JETAN/jetan2004/AliJFJetTrigger.cxx

   1 // $Id$
   2
   3 #include <Riostream.h>
   4
   5 #include <TClonesArray.h>
   6 #include <TIterator.h>
   7 #include <TParticle.h>
   8
   9 #include "AliJFJetTrigger.h"
  10
  11 ClassImp(AliJFJetTrigger)
  12
  13 AliJFJetTrigger::AliJFJetTrigger(Int_t n) : AliJFJetFinder(n)
  14 {
  15   fParticles=new TClonesArray("TParticle",100000);
  16 }
  17
  18 AliJFJetTrigger::~AliJFJetTrigger()
  19 {
  20   delete fParticles;
  21 }
  22
  23 Bool_t AliJFJetTrigger::IsAcceptedParticle(TParticle *p)
  24 {
  25   if(p->GetStatusCode()%100!=1) return kFALSE;
  26
  27   Int_t pcode=p->GetPdgCode();
  28
  29   if((!fEM) && ((pcode==11)||(pcode==-11)||(pcode==22))) return kFALSE;
  30
  31   TParticlePDG *pdg=p->GetPDG();
  32   Float_t ch=pdg->Charge();
  33   if((!fCharged)&&(ch)) return kFALSE;
  34   if((!fNeutral)&&(!ch)) return kFALSE;
  35
  36   Float_t eta=p->Eta();
  37   if((eta<fEtaMin)||(eta>fEtaMax)) return kFALSE;
  38
  39   Float_t phi=p->Phi();
  40   if((phi<fPhiMin)||(phi>fPhiMax)) return kFALSE;
  41
  42   Float_t pt=p->Pt();
  43   if((pt<fPtMin)||(pt>fPtMax)) return kFALSE;
  44
  45   return kTRUE;
  46 }
  47