Page MenuHomeHEPForge
Files F7879145

No OneTemporary
Actions

View Options

diff --git a/src/Analyses/CDF_2015_1388868.cc b/src/Analyses/CDF_2015_1388868.cc
--- a/src/Analyses/CDF_2015_1388868.cc
+++ b/src/Analyses/CDF_2015_1388868.cc
@@ -1,132 +1,132 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Tools/Logging.hh"
namespace Rivet {
/// @brief CDF leading track underlying event at 300, 900 and 1960 GeV
/// @author Orestes Tumbarell Aranda (Havana), Hannes Jung (DESY)
class CDF_2015_I1388868 : public Analysis {
public:
/// Constructor
DEFAULT_RIVET_ANALYSIS_CTOR(CDF_2015_I1388868);
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
// Energy selection
double isqrts = 0;
if (fuzzyEquals(sqrtS()/GeV, 300, 1E-3)) {
isqrts = 3;
} else if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) {
isqrts = 2;
} else if (fuzzyEquals(sqrtS()/GeV, 1960, 1E-3)) {
isqrts = 1;
} else {
throw UserError("Unexpected sqrtS ! Only 300, 900, 1960 GeV is supported by CDF_2015_I1388868");
}
MSG_DEBUG("CDF Tevatron UE: running with " << sqrtS()/GeV);
// Book projection
const ChargedFinalState cfs(Cuts::abseta < 0.8 && Cuts::pT > 0.5*GeV);
- addProjection(cfs, "Tracks");
+ declare(cfs, "Tracks");
// Book profile histos
_NchgPDFden1 = bookProfile1D(isqrts,1,1);
_NchgPMNden1 = bookProfile1D(isqrts,1,2);
_NchgPMXden1 = bookProfile1D(isqrts,1,3);
_NchgPden1 = bookProfile1D(isqrts,1,4);
_PTsumPDFden1 = bookProfile1D(isqrts,1,6);
_PTsumPMNden1 = bookProfile1D(isqrts,1,7);
_PTsumPMXden1 = bookProfile1D(isqrts,1,8);
_PTsumPden1 = bookProfile1D(isqrts,1,9);
}
/// Perform the per-event analysis
void analyze(const Event& event) {
// Require at least one track in the event with pT >= 0.5 GeV
- const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(event, "Tracks");
+ const ChargedFinalState& cfs = apply<ChargedFinalState>(event, "Tracks");
if (cfs.empty()) vetoEvent;
const Particles trks = cfs.particlesByPt();
// Get lead track
const Particle p_lead = trks[0];
const double philead = p_lead.phi();
const double ptlead = p_lead.pT();
// Loop over tracks and compute variables
double NchgP1 = 0, NchgP2 = 0, PTsumP1 = 0, PTsumP2 = 0;
for (const Particle& p : trks) {
// Region definition -- if not in transverse region, ignore
const double dphi = mapAngle0To2Pi(p.phi() - philead);
if (!inRange(dphi, PI/3, 2*PI/3) && !inRange(dphi, 4*PI/3, 5*PI/3)) continue;
// Transverse region 1
if (inRange(dphi, PI/3, 2*PI/3)) {
NchgP1 += 1;
PTsumP1 += p.pT();
}
// Transverse region 2
else if (inRange(dphi, 4*PI/3, 5*PI/3)) {
NchgP2 += 1;
PTsumP2 += p.pT();
}
}
// Calculate total variables
const double NchgPtot = (NchgP1 + NchgP2)/2;
const double NchgPmax = max(NchgP1,NchgP2);
const double NchgPmin = min(NchgP1,NchgP2);
const double PTsumPtot = (PTsumP1 + PTsumP2)/2;
const double PTsumPmax = max(PTsumP1,PTsumP2);
const double PTsumPmin = min(PTsumP1,PTsumP2);
//
const double PTsumPMXden = PTsumPmax/AREA;
const double PTsumPMNden = PTsumPmin/AREA;
const double NchgPMXden = NchgPmax/AREA;
const double NchgPMNden = NchgPmin/AREA;
//
const double NchgPDFden = NchgPMXden - NchgPMNden;
const double PTsumPDFden = PTsumPMXden - PTsumPMNden;
// Fill histograms
const double weight = event.weight();
_NchgPden1 ->fill(ptlead/GeV, NchgPtot/AREA, weight );
_NchgPMXden1->fill(ptlead/GeV, NchgPmax/AREA, weight );
_NchgPMNden1->fill(ptlead/GeV, NchgPmin/AREA, weight );
_NchgPDFden1->fill(ptlead/GeV, NchgPDFden , weight );
_PTsumPden1 ->fill(ptlead/GeV, PTsumPtot/AREA, weight );
_PTsumPMXden1->fill(ptlead/GeV, PTsumPmax/AREA, weight );
_PTsumPMNden1->fill(ptlead/GeV, PTsumPmin/AREA, weight );
_PTsumPDFden1->fill(ptlead/GeV, PTsumPDFden , weight );
}
//@}
/// eta-phi area of the transverse region
constexpr static double AREA = 2*0.8 * M_PI/3;
/// @name Histograms
Profile1DPtr _NchgPden1, _NchgPMXden1,_NchgPMNden1,_NchgPDFden1,_PTsumPden1,_PTsumPMXden1,_PTsumPMNden1,_PTsumPDFden1;
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(CDF_2015_I1388868);
}

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 7:39 PM (1 d, 9 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3801056
Default Alt Text
(4 KB)

Event Timeline