2 #include "Rivet/Analysis.hh" 3 #include "Rivet/Projections/FinalState.hh" 4 #include "Rivet/Projections/FastJets.hh" 5 #include "Rivet/Projections/WFinder.hh" 6 #include "Rivet/Projections/ZFinder.hh" 7 #include "fastjet/tools/Filter.hh" 8 #include "fastjet/tools/Pruner.hh" 21 : Analysis(
"CMS_2013_I1224539_DIJET"),
37 FinalState fs(-2.4, 2.4, 0*
GeV);
38 addProjection(fs,
"FS");
41 addProjection(FastJets(fs, FastJets::ANTIKT, 0.7),
"JetsAK7");
42 addProjection(FastJets(fs, FastJets::CAM, 0.8),
"JetsCA8");
43 addProjection(FastJets(fs, FastJets::CAM, 1.2),
"JetsCA12");
58 const double ptBins_dj[
N_PT_BINS_dj+1] = { 220.0, 300.0, 450.0, 500.0, 600.0, 800.0, 1000.0, 1500.0};
60 if (inRange(ptJ, ptBins_dj[ibin], ptBins_dj[ibin+1]))
return ibin;
68 const double weight =
event.weight();
71 const PseudoJets& psjetsAK7 = applyProjection<FastJets>(
event,
"JetsAK7").pseudoJetsByPt( 50.0*
GeV );
72 if (psjetsAK7.size() < 2) vetoEvent;
75 const fastjet::PseudoJet& j0 = psjetsAK7[0];
76 const fastjet::PseudoJet& j1 = psjetsAK7[1];
77 double ptAvg = 0.5 * (j0.pt() + j1.pt());
84 fastjet::PseudoJet filtered0 =
_filter(j0);
85 fastjet::PseudoJet filtered1 =
_filter(j1);
86 fastjet::PseudoJet trimmed0 =
_trimmer(j0);
87 fastjet::PseudoJet trimmed1 =
_trimmer(j1);
88 fastjet::PseudoJet pruned0 =
_pruner(j0);
89 fastjet::PseudoJet pruned1 =
_pruner(j1);
101 const double normalizationVal = 1000;
const fastjet::Filter _trimmer
Histo1DPtr _h_filteredAvgJetMass_dj[N_PT_BINS_dj]
Histo1DPtr _h_ungroomedJet1pt
size_t findPtBin(double ptJ)
const fastjet::Filter _filter
Histo1DPtr _h_ungroomedAvgJetMass_dj[N_PT_BINS_dj]
Histo1DPtr _h_prunedAvgJetMass_dj[N_PT_BINS_dj]
const fastjet::Pruner _pruner
DECLARE_RIVET_PLUGIN(CMS_2013_I1224539_DIJET)
void analyze(const Event &event) override
Perform the per-event analysis.
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
enum Rivet::CMS_2013_I1224539_DIJET::@758 BINS_dj
void init() override
Book histograms and initialise projections before the run.
Histo1DPtr _h_ungroomedJet0pt
void finalize() override
Normalise histograms etc., after the run.
CMS_2013_I1224539_DIJET()
Constructor.
Histo1DPtr _h_trimmedAvgJetMass_dj[N_PT_BINS_dj]