/* $Id$ */
+/* History of cvs commits:
+ *
+ * $Log$
+ * Revision 1.36 2005/05/28 14:19:04 schutz
+ * Compilation warnings fixed by T.P.
+ *
+ */
+
//_________________________________________________________________________
// A Particle modified by PHOS response and produced by AliPHOSvFast
+// This is also a base class for AliPHOSRecParticle produced by AliPHOSPIDv1
+// Defines the particle type
// To become a general class of AliRoot ?
-//
+//--
//*-- Author: Yves Schutz (SUBATECH)
// --- ROOT system ---
+class TClonesArray;
#include "TParticle.h"
-#include "TVector3.h"
-
+#include "AliPID.h"
// --- Standard library ---
// --- AliRoot header files ---
AliPHOSFastRecParticle(const AliPHOSFastRecParticle & rp) ; // ctor
AliPHOSFastRecParticle(const TParticle & p) ; // ctor
- virtual ~AliPHOSFastRecParticle(){
- // dtor
- }
+ virtual ~AliPHOSFastRecParticle(){ } //dtor
+ AliPHOSFastRecParticle & operator = (const AliPHOSFastRecParticle & /*rp*/);
+
virtual Int_t DistancetoPrimitive(Int_t px, Int_t py) ;
virtual void Draw(Option_t *option) ;
virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py) ;
// returns the index of this in the list
return fIndexInList ;
}
- virtual const Int_t GetNPrimaries() const {return 0 ;}
- virtual const TParticle * GetPrimary(Int_t index=0) const {return 0 ;}
- const Int_t GetType() const {
+ virtual Int_t GetNPrimaries() const {return 0 ;}
+ virtual const TParticle * GetPrimary(Int_t) const {return 0 ;}
+ Int_t GetType() const {
// returns the type of the particle
return fType ;
}
- void SetPIDBit(UInt_t fSet)
- {
- fType |= (1<<fSet) ;
- }
-
- Bool_t TestPIDBit(UInt_t fTest){
- if (fType & (1<<fTest) )
- return kTRUE ;
- else
- return kFALSE ;
- }
-
- Bool_t IsPhotonHiPuLoEf() {
- Bool_t pid=kFALSE ;
- if(TestPIDBit(8)&& //PCA
- TestPIDBit(5)&& //TOF
- TestPIDBit(2)) //RCPV
- pid = kTRUE;
- return pid ;
- }
-
- Bool_t IsPhotonMedPuEf(){
- Bool_t pid=kFALSE ;
- if(TestPIDBit(7)&& //PCA
- TestPIDBit(4)&& //TOF
- TestPIDBit(1)) //RCPV
- pid = kTRUE ;
- return pid ;
- }
+ void SetPIDBit(UInt_t fSet) {
+ // Set PID bit number fSet
+ fType |= (1<<fSet) ;
+ }
- Bool_t IsPhotonHiEfLoPu() {
- Bool_t pid=kFALSE ;
- if(TestPIDBit(6)&& //PCA
- TestPIDBit(3)&& //TOF
- TestPIDBit(0)) //RCPV
- pid = kTRUE ;
- return pid ;
- }
- Bool_t IsPhoton() {
- Bool_t pid=kFALSE ;
- if(IsPhotonHiEfLoPu()) pid = kTRUE ;
- return pid ;
+ Bool_t TestPIDBit(UInt_t fTest) const {
+ // Check PID bit number fTest
+ if (fType & (1<<fTest) ) return kTRUE ;
+ else return kFALSE ;
}
- Bool_t IsFastChargedHadron() {
- Bool_t pid=kFALSE ;
- if(TestPIDBit(5)&&TestPIDBit(4)&&TestPIDBit(3)) //TOF
- pid = kTRUE ;
- return pid ;
- }
- Bool_t IsSlowChargedHadron() {
- Bool_t pid=kFALSE ;
- if(TestPIDBit(1)||TestPIDBit(0)) //CPV
- pid = kTRUE ;
- return pid ;
- }
- Bool_t IsFastNeutralHadron() {
- Bool_t pid=kFALSE ;
- if(TestPIDBit(5)&&TestPIDBit(4)&&TestPIDBit(3)&& //TOF
- TestPIDBit(2)&&TestPIDBit(1)&&TestPIDBit(0))//RCPV
- pid = kTRUE ;
- return pid ;
- }
- Bool_t IsSlowNeutralHadron() {
- Bool_t pid=kFALSE ;
- if(TestPIDBit(2)&&TestPIDBit(1)&&TestPIDBit(0))//RCPV
- pid = kTRUE ;
- return pid ;
- }
+ Bool_t IsPhoton (TString purity = "low") const;
+ Bool_t IsPi0 (TString purity = "low") const;
+ Bool_t IsElectron (TString purity = "low") const;
+ Bool_t IsHardPhoton () const;
+ Bool_t IsHardPi0 () const;
+ Bool_t IsHadron () const;
+ Bool_t IsChargedHadron () const;
+ Bool_t IsNeutralHadron () const;
+ Bool_t IsFastChargedHadron() const;
+ Bool_t IsSlowChargedHadron() const;
+ Bool_t IsFastNeutralHadron() const;
+ Bool_t IsSlowNeutralHadron() const;
+ Bool_t IsEleCon(TString purity = "low") const;
- Bool_t IsFastChargedEM() {
- Bool_t pid=kFALSE ;
- if( (TestPIDBit(8)||TestPIDBit(7)||TestPIDBit(6))&& // PCA
- (TestPIDBit(5)||TestPIDBit(4)||TestPIDBit(3))&& // TOF
- !(TestPIDBit(2)||TestPIDBit(1)||TestPIDBit(0))) // !CPV
- pid = kTRUE ;
- return pid ;
- }
-
- Bool_t IsSlowChargedEM() {
- Bool_t pid=kFALSE ;
- if(TestPIDBit(8)||TestPIDBit(7)||TestPIDBit(6))
- pid = kTRUE ;
- return pid ;
- }
-
TString Name() const ;
virtual void Paint(Option_t * option="");
- virtual void Print(Option_t * option = "") const ;
+ virtual void Print(const Option_t * = "") const ;
+ void SetTof(Float_t tof) { fTof = tof ; }
+ Float_t ToF() const { return fTof ; }
void SetType(Int_t type) ;
void SetIndexInList(Int_t val) {
fIndexInList = val ;
}
//This has to disappear
- enum EParticleType { kUNDEFINED=-1,
+ enum EParticleType { kTYPE = 8,
+ kUNDEFINED=-1,
kNEUTRALEMFAST, kNEUTRALHAFAST, kNEUTRALEMSLOW, kNEUTRALHASLOW,
kCHARGEDEMFAST, kCHARGEDHAFAST, kCHARGEDEMSLOW, kCHARGEDHASLOW } ;
-
typedef TClonesArray FastRecParticlesList ;
protected:
Int_t fIndexInList ; // the index of this RecParticle in the list stored in TreeR (to be set by analysis)
+ Float_t fTof ; // time of fliht
Int_t fType ; // particle type obtained by "virtual" reconstruction
- private:
+ Float_t fPID[AliPID::kSPECIESCN] ; // PID probability densities
+ private:
ClassDef(AliPHOSFastRecParticle,2) // Reconstructed Particle produced by the fast simulation