1 #ifndef ALIJETPARTICLESREADERKINE_H
2 #define ALIJETPARTICLESREADERKINE_H
6 //_______________________________________________________________________
7 /////////////////////////////////////////////////////////////////////////
9 // class AliJetParticlesReaderKine
11 // Reader for Kinematics
13 // loizides@ikf.uni-frankfurt.de
15 /////////////////////////////////////////////////////////////////////////
17 #include <Riostream.h>
19 #include <TParticle.h>
20 #include "AliJetParticlesReader.h"
23 class AliJetParticlesReaderKine: public AliJetParticlesReader
26 AliJetParticlesReaderKine();
27 AliJetParticlesReaderKine(TString&);
28 AliJetParticlesReaderKine(TObjArray*, const Char_t *filename="galice.root");
29 virtual ~AliJetParticlesReaderKine();
32 void SetNeutral(Bool_t b){fNeutral=b;}
33 void SetCharged(Bool_t b){fCharged=b;}
34 void SetEM(Bool_t b){fEM=b;}
35 void SetUseTracks(Bool_t b){fUseTracks=b;}
39 Int_t OpenFile(Int_t n);
40 Bool_t IsAcceptedParticle(TParticle *p) const;
42 TString fFileName; //file name of galice
43 AliRunLoader* fRunLoader; //!pointer to loader
45 Bool_t fNeutral; //neutral cut
46 Bool_t fCharged; //charged cut
47 Bool_t fEM; //em (e+,e-,gamma) cut
48 Bool_t fUseTracks; // use ntracks instead of nprimaries
49 ClassDef(AliJetParticlesReaderKine,1)
52 inline Bool_t AliJetParticlesReaderKine::IsAcceptedParticle(TParticle *p) const
54 Int_t pcode=TMath::Abs(p->GetPdgCode());
57 if ((pcode==11)||(pcode==22)) {
58 if(!fEM) return kFALSE;
60 /*this is slighly misleading, needs fixing*/
61 if(pcode!=211 && pcode!=321 && pcode!=2212 &&
62 pcode!=111 && pcode!=311 && pcode!=2112 &&
63 pcode!=3122 && pcode!=213 && pcode!=113 &&
64 pcode!=130 &&pcode!=310) {
68 TParticlePDG *pdg=p->GetPDG();
69 Int_t ch=(Int_t)pdg->Charge();
70 if((!fCharged)&&(ch)) return kFALSE;
71 if((!fNeutral)&&(!ch)) return kFALSE;
74 if(pcode==11 && pcode!=22 && pcode <= 110) return kFALSE;
75 if(pcode==990 || pcode > 5224) return kFALSE;
76 //if(pcode!=11 && pcode!=22 && pcode!=211 && pcode!=321 &&
77 // pcode!=2212 && pcode!=111 && pcode!=211 && pcode!=2112)
78 // return kFALSE; /*keep only e-, gammas, pi, ks, p, n*/
80 if(fEM && (pcode==11 || pcode==22 || pcode==111)){
81 // em particles are accepted
83 TParticlePDG *pdg=p->GetPDG();
84 Int_t ch=(Int_t)pdg->Charge();
85 if((!fCharged)&&(ch)) return kFALSE;
86 if((!fNeutral)&&(!ch)) return kFALSE;
95 if(pabs-TMath::Abs(pz)>1.e-8) eta=0.5*TMath::Log((pabs+pz)/(pabs-pz));
98 if((eta<fEtaMin)||(eta>fEtaMax)) return kFALSE;
100 Float_t phi=p->Phi();
101 if((phi<fPhiMin)||(phi>fPhiMax)) return kFALSE;
104 if((pt<fPtMin)||(pt>fPtMax)) return kFALSE;