Line data    Source code

       1             : #ifndef ALIEMCALPID_H
       2             : #define ALIEMCALPID_H
       3             : 
       4             : //_________________________________________________________________________
       5             : /// \class AliEMCALPID
       6             : /// \brief Compute cluster PID weights 
       7             : ///
       8             : ///   Compute PID weights for all the clusters that are in AliESDs.root file
       9             : ///   the AliESDs.root have to be in the same directory as the class.
      10             : ///   Wrapper class to be used during reconstruction, the main algorithm sits in 
      11             : ///   EMCALUtils/AliEMCALPIDUtils.
      12             : ///
      13             : ///   Implementation based on simulations before 2009, executed in reconstruction but
      14             : ///   not really used after 2009. To be revisited.
      15             : ///
      16             : ///   Do:    
      17             : ///   AliEMCALPID *pid = new AliEMCALPID(kFALSE); // this calls the constructor which avoids the call to recparam 
      18             : ///   pid->SetReconstructor(kFALSE);
      19             : ///   pid->SetPrintInfo(kTRUE);
      20             : ///   pid->SetHighFluxParam(); //   pid->SetLowFluxParam(); 
      21             : ///   
      22             : ///   then in cluster loop do
      23             : ///   pid->ComputePID(energy, lambda0);
      24             : ///     
      25             : ///   Get back the probabilities with 
      26             : ///   pid->GetPIDFinal(pidFinal) 
      27             : ///
      28             : ///   where Double_t pidFinal[AliPID::kSPECIESCN] is the standard PID for :
      29             : ///
      30             : ///       kElectron :  fPIDFinal[0]
      31             : ///       kMuon     :  fPIDFinal[1]
      32             : ///       kPion     :  fPIDFinal[2]
      33             : ///       kKaon     :  fPIDFinal[3]
      34             : ///       kProton   :  fPIDFinal[4]
      35             : ///       kPhoton   :  fPIDFinal[5]
      36             : ///       kPi0      :  fPIDFinal[6]
      37             : ///       kNeutron  :  fPIDFinal[7]
      38             : ///       kKaon0    :  fPIDFinal[8]
      39             : ///       kEleCon   :  fPIDFinal[9]
      40             : ///       kUnknown  :  fPIDFinal[10]
      41             : ///
      42             : /// \author Genole Bourdaud, SUBATECH
      43             : /// First implementation, 2007
      44             : /// \author Marie Germain, <Marie.Germain@subatech.in2p3.fr>, SUBATECH 
      45             : /// New parametrization for low and high flux environment, 07/2009
      46             : /// \author Gustavo Conesa, <Gustavo.Conesa.Balbastre@cern.ch>, LNF, 08/2009 
      47             : /// Divide class in AliEMCALPID and AliEMCALPIDUtils, PIDUtils belong to library EMCALUtils 
      48             : ///
      49             : 
      50             : // Root includes
      51             : class TArrayD ;
      52             : 
      53             : // AliRoot includes
      54             : class AliESDEvent ;
      55             : #include "AliEMCALPIDUtils.h" 
      56             : 
      57          48 : class AliEMCALPID : public AliEMCALPIDUtils {
      58             : 
      59             : public:
      60             :   
      61             :   AliEMCALPID();
      62             :   
      63             :   AliEMCALPID(Bool_t reconstructor);
      64             :   
      65             :   //virtual ~AliEMCALPID() { }
      66             :   
      67             :   void    RunPID(AliESDEvent *esd);
      68             :   
      69             :   void    InitParameters();
      70             :   
      71           8 :   void    SetReconstructor(Bool_t yesno) { fReconstructor = yesno ; }
      72             :         
      73             :  private:
      74             :   
      75             :   Bool_t   fReconstructor ;  ///< Fill AliESDCaloCluster object when called from AliEMCALReconstructor
      76             :   
      77             :   /// \cond CLASSIMP
      78          42 :   ClassDef(AliEMCALPID, 5) ;
      79             :   /// \endcond
      80             : 
      81             : };
      82             : 
      83             : #endif // ALIEMCALPID_H
      84             : 
      85             :