]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PHOS/AliPHOSFastRecParticle.h
make sure EMCal trigger-cluster matching is done when event is EMCal triggered, in...
[u/mrichter/AliRoot.git] / PHOS / AliPHOSFastRecParticle.h
index 9c270cef3e8a4b40f54d6d20ca479944c6240c4c..d9a55e165f449f4573e0d36af61cdb0b743e5514 100644 (file)
@@ -5,17 +5,27 @@
 
 /* $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 ---
@@ -28,9 +38,9 @@ class AliPHOSFastRecParticle : public TParticle {
   
   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) ;
@@ -38,125 +48,66 @@ class AliPHOSFastRecParticle : public TParticle {
     // 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 IsPhotonHiPu_LoEf()  {
-    Bool_t pid=kFALSE ;
-    if(TestPIDBit(8)&& //PCA
-       TestPIDBit(5)&& //TOF
-       TestPIDBit(2))  //RCPV
-      pid = kTRUE;
-    return pid ;
-  }
-  
-  Bool_t IsPhotonMed_Pu_Ef(){
-    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 IsPhotonHiEf_LoPu()  {
-    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(IsPhotonHiEf_LoPu()) 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))&&
-       TestPIDBit(5)&&TestPIDBit(4)&&TestPIDBit(3))//TOF
-      pid = kTRUE ;
-    return pid ;
-  }
-  
-  Bool_t IsSlowChargedEM()  {
-    Bool_t pid=kFALSE ;
-    if(TestPIDBit(8)||TestPIDBit(7)||TestPIDBit(6))
-      pid = kTRUE ;
-    return pid ;
-  }
-  
-  TString Name() ; 
+  TString Name() const ; 
   virtual void Paint(Option_t * option="");
-  virtual void Print(const char * opt) ; 
+  virtual void Print(const Option_t * = "") const ; 
   
-  void SetType(Int_t type) { 
-    // sets the particle type 
-    fType = type ; 
-  }        
+  void SetTof(Float_t tof) { fTof = tof ; } 
+  Float_t ToF() const { return fTof ; } 
+  void SetType(Int_t type) ;
   
   void SetIndexInList(Int_t val) { 
     // sets the value of the index in the list 
     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