Rivet 3.1.9
MC_Cent_pPb.hh
1// -*- C++ -*-
2#ifndef RIVET_MC_Cent_pPb_HH
3#define RIVET_MC_Cent_pPb_HH
4
5#include "Rivet/Projections/FinalState.hh"
6#include "Rivet/Projections/ChargedFinalState.hh"
7#include "Rivet/Projections/SingleValueProjection.hh"
8#include "Rivet/Projections/TriggerProjection.hh"
9
10namespace Rivet {
11
15public:
16
19 declare(FinalState(Cuts::eta < -3.2 && Cuts::eta > -4.9 && Cuts::pT > 0.1*GeV),
20 "FSSumETFwdCentrality");
21 }
22
25
27 using Projection::operator =;
28
29protected:
30
32 void project(const Event& e) {
33 clear();
34 const FinalState & fsfwd =
35 apply<FinalState>(e, "FSSumETFwdCentrality");
36 double estimate = 0.0;
37 for ( const Particle & p : fsfwd.particles() ) {
38 estimate += p.Et();
39 }
40 set(estimate);
41 }
42
44 CmpState compare(const Projection& p) const {
45 return mkNamedPCmp(p, "FSSumETFwdCentrality");
46 }
47
48};
49
50
54public:
55
58 declare(FinalState(Cuts::eta < -3.2 && Cuts::eta > -4.9 && Cuts::pT > 0.1*GeV),
59 "FSSumETFwdCentrality");
60 declare(ChargedFinalState(Cuts::eta > 2.09 && Cuts::eta < 3.84 &&
61 Cuts::pT > 0.1*GeV), "MBB");
62 declare(ChargedFinalState(Cuts::eta < -2.09 && Cuts::eta > -3.84 &&
63 Cuts::pT > 0.1*GeV), "MBF");
64 }
65
68
70 using Projection::operator =;
71
72protected:
73
75 void project(const Event& event) {
76 pass();
77 if ( applyProjection<FinalState>(event,"MBF").particles().empty() ||
78 applyProjection<FinalState>(event,"MBB").particles().empty() )
79 fail();
80 }
81
83 CmpState compare(const Projection& p) const {
84 return mkNamedPCmp(p, "MBF") || mkNamedPCmp(p, "MBB");
85 }
86
87};
88
89
90}
91
92#endif
Project only charged final state particles.
Definition ChargedFinalState.hh:11
Representation of a HepMC event, and enabler of Projection caching.
Definition Event.hh:22
Project out all final-state particles in an event. Probably the most important projection in Rivet!
Definition FinalState.hh:12
Definition MC_Cent_pPb.hh:14
void project(const Event &e)
Perform the projection on the Event.
Definition MC_Cent_pPb.hh:32
CmpState compare(const Projection &p) const
Compare projections.
Definition MC_Cent_pPb.hh:44
MC_SumETFwdPbCentrality()
Constructor.
Definition MC_Cent_pPb.hh:18
DEFAULT_RIVET_PROJ_CLONE(MC_SumETFwdPbCentrality)
Clone on the heap.
Definition MC_Cent_pPb.hh:53
void project(const Event &event)
Perform the projection on the Event.
Definition MC_Cent_pPb.hh:75
DEFAULT_RIVET_PROJ_CLONE(MC_pPbMinBiasTrigger)
Clone on the heap.
CmpState compare(const Projection &p) const
Compare projections.
Definition MC_Cent_pPb.hh:83
MC_pPbMinBiasTrigger()
Constructor.
Definition MC_Cent_pPb.hh:57
virtual const Particles & particles() const
Get the particles in no particular order, with no cuts.
Definition ParticleFinder.hh:52
Particle representation, either from a HepMC::GenEvent or reconstructed.
Definition Particle.hh:53
const PROJ & declare(const PROJ &proj, const std::string &name)
Register a contained projection (user-facing version)
Definition ProjectionApplier.hh:170
Base class for all Rivet projections.
Definition Projection.hh:29
Cmp< Projection > mkNamedPCmp(const Projection &otherparent, const std::string &pname) const
Base class for projections returning a single floating point value.
Definition SingleValueProjection.hh:17
void set(double v)
Definition SingleValueProjection.hh:55
void clear()
Unset the value.
Definition SingleValueProjection.hh:58
Base class for projections returning a bool corresponding to a trigger.
Definition TriggerProjection.hh:17
void fail()
Indicate that the event has failed the trigger.
Definition TriggerProjection.hh:56
void pass()
Indicate that the event has passed the trigger.
Definition TriggerProjection.hh:46
double p(const ParticleBase &p)
Unbound function access to p.
Definition ParticleBaseUtils.hh:684
Definition MC_Cent_pPb.hh:10