22 PI = 3.141592653589793238462643383279502884197169399375105820974945;
34 return mcisocut.IsIsolatedPP(pp);
42 if ( pp.
pdgId() != 22)
54 cout <<
" found a prompt photon" << endl;
65 {
PI = 3.141592653589793238462643383279502884197169399375105820974945;
73 PI = 3.141592653589793238462643383279502884197169399375105820974945;
86 double PtMaxHadron = 0;
89 for(
int i=0;
i < maxindex ; ++
i) {
99 if(apid>11 && apid<20)
108 bool isHadron =
false;
109 if ( fabs(pp.
pdgId())>100 && fabs(pp.
pdgId()) != 310)
112 if(apid>100 && apid!=310 && pt>PtMaxHadron)
114 if ( (isHadron ==
true) && (fabs(pp.
pt()-
pt)<0.001) && (pp.
pdgId()==p.
pdgId()) )
124 if(EtCone-EtPhoton > 5+EtPhoton/20)
126 if(PtMaxHadron > 4.5+EtPhoton/40)
137 using namespace reco;
146 for(
int i=0;
i < maxindex ; ++
i) {
171 if(EtCone-EtPhoton > 5)
183 using namespace reco;
192 for(
int i=0;
i < maxindex ; ++
i) {
213 if(EtCone-EtPhoton > 5)
224 while(dPhi >=
PI) dPhi -= (2.0*
PI);
225 while(dPhi < (-1.0*
PI)) dPhi += (2.0*
PI);
227 return sqrt(dEta*dEta+dPhi*dPhi);
234 while(dPhi >=
PI) dPhi -= (2.0*
PI);
235 while(dPhi < (-1.0*
PI)) dPhi += (2.0*
PI);
std::vector< GenParticle > GenParticleCollection
collection of GenParticles
double getDeltaR(const reco::Candidate &track1, const reco::Candidate &track2)
int pdgId() const final
PDG identifier.
double pt() const final
transverse momentum
bool IsIsolated(const reco::GenParticle &pp)
double et() const final
transverse energy
bool IsIsolatedJP(const reco::GenParticle &pp)
virtual int pdgId() const =0
PDG identifier.
Abs< T >::type abs(const T &t)
virtual double eta() const =0
momentum pseudorapidity
T const * product() const
bool IsIsolatedPP(const reco::GenParticle &pp)
HiPhotonType(edm::Handle< reco::GenParticleCollection > inputHandle)
bool IsPrompt(const reco::GenParticle &pp)
et
define resolution functions of each parameter
double getDeltaPhi(const reco::Candidate &track1, const reco::Candidate &track2)
int status() const final
status word
virtual double phi() const =0
momentum azimuthal angle
bool IsIsolated(const reco::GenParticle &pp)
const reco::GenParticleCollection * fSigParticles
const Candidate * mother(size_type=0) const override
return mother at a given position, i = 0, ... numberOfMothers() - 1 (read only mode) ...