AliAnaElectron.h

Go to the documentation of this file.

#ifndef ALIANAELECTRON_H
#define ALIANAELECTRON_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice     */

//_________________________________________________________________________
//_________________________________________________________________________

// --- ROOT system ---
class TH2F ;
class TH1F;
class TH3D;
class TObjString;

// --- ANALYSIS system ---
#include "AliAnaCaloTrackCorrBaseClass.h"

class TList ;

class AliAnaElectron : public AliAnaCaloTrackCorrBaseClass {

 public: 
  
  AliAnaElectron() ;
  
  virtual ~AliAnaElectron() { ; }
  
  //---------------------------------------
  // General analysis frame methods
  //---------------------------------------
  
  TObjString * GetAnalysisCuts();
  
  TList      * GetCreateOutputObjects();
  
  void         Init();

  void         InitParameters();

  void         MakeAnalysisFillAOD()  ;

  void         MakeAnalysisFillHistograms() ; 
  
  void         Print(const Option_t * opt)const;
    
  
  // Analysis methods
  
  Bool_t       ClusterSelected(AliVCluster* cl, Int_t nMaxima) ;
  
  void         FillShowerShapeHistograms( AliVCluster* cluster, Int_t mcTag , Int_t pidTag) ;
  
  void         SwitchOnFillShowerShapeHistograms()    { fFillSSHistograms = kTRUE  ; }
  void         SwitchOffFillShowerShapeHistograms()   { fFillSSHistograms = kFALSE ; }  
    
  void         WeightHistograms(AliVCluster *clus);
  
  void         SwitchOnFillWeightHistograms()         { fFillWeightHistograms = kTRUE  ; }
  void         SwitchOffFillWeightHistograms()        { fFillWeightHistograms = kFALSE ; }  
  
  //---------------------------------------
  // Analysis parameters setters getters
  //---------------------------------------
  
  // ** Cluster selection methods **
  
  void         SetdEdxCut(Double_t min, Double_t max) { fdEdxMin    = min ; 
                                                        fdEdxMax    = max          ; }
  
  void         SetEOverP(Double_t min, Double_t max)  { fEOverPMin  = min ; 
                                                        fEOverPMax  = max          ; }

  
  void         SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) {
                fMinDist = m1; fMinDist2 = m2; fMinDist3 = m3; }

  void         SetTimeCut(Double_t min, Double_t max) { fTimeCutMin = min; 
                                                        fTimeCutMax = max          ; }
  Double_t     GetTimeCutMin()                  const { return fTimeCutMin         ; }
  Double_t     GetTimeCutMax()                  const { return fTimeCutMax         ; }  
  
  void         SetNCellCut(Int_t n)                   { fNCellsCut = n             ; }
  Double_t     GetNCellCut()                    const { return fNCellsCut          ; }
  
  void         SetNLMCut(Int_t min, Int_t max)        { fNLMCutMin = min;
                                                        fNLMCutMax = max                ; }
  Int_t        GetNLMCutMin()                   const { return fNLMCutMin               ; }
  Int_t        GetNLMCutMax()                   const { return fNLMCutMax               ; }
  
  void         FillNOriginHistograms(Int_t n)         { fNOriginHistograms = n ; 
                                                        if(n > 10) fNOriginHistograms = 10; }

  
  void         FillAODWithElectrons()                 { fAODParticle = AliCaloPID::kElectron      ; }
  void         FillAODWithHadrons()                   { fAODParticle = AliCaloPID::kChargedHadron ; }
  void         FillAODWithAny()                       { fAODParticle = 0 ; }

  void         SwitchOnOnlySimpleSSHistoFill()        { fFillOnlySimpleSSHisto = kTRUE  ; }
  void         SwitchOffOnlySimpleHistoFill()         { fFillOnlySimpleSSHisto = kFALSE ; }
  
  enum mcTypes    { kmcPhoton = 0,        kmcPi0Decay = 1,       kmcOtherDecay = 2,  
                    kmcPi0 = 3,           kmcEta = 4,            kmcElectron = 5,       
                    kmcConversion = 6,    kmcOther = 7,          kmcAntiNeutron = 8,    
                    kmcAntiProton = 9                                                 };    
  
  enum mcssTypes  { kmcssPhoton = 0,      kmcssOther = 1,       kmcssPi0 = 2,         
                    kmcssEta = 3,         kmcssConversion = 4,  kmcssElectron = 5       };  
  
  private:
 
  Float_t  fMinDist ;                           
  Float_t  fMinDist2;                           
  Float_t  fMinDist3;                           
  Double_t fTimeCutMin  ;                       
  Double_t fTimeCutMax  ;                       
  Int_t    fNCellsCut ;                         
  Int_t    fNLMCutMin  ;                        
  Int_t    fNLMCutMax  ;                        
  Bool_t   fFillSSHistograms ;                  
  Bool_t   fFillOnlySimpleSSHisto;              
  Bool_t   fFillWeightHistograms ;              
  Int_t    fNOriginHistograms;                  

  Float_t  fdEdxMin;                            
  Float_t  fdEdxMax;                            
  Float_t  fEOverPMin;                          
  Float_t  fEOverPMax;                          

  Int_t    fAODParticle;                        
  
  TLorentzVector fMomentum;                     
  TLorentzVector fMomentumMC;                   
  TVector3       fProdVertex;                   

  //Histograms
  TH2F * fhdEdxvsE;                             
  TH2F * fhdEdxvsP;                             
  TH2F * fhEOverPvsE;                           
  TH2F * fhEOverPvsP;                           

  TH2F * fhdEdxvsECutM02;                       
  TH2F * fhdEdxvsPCutM02;                       
  TH2F * fhEOverPvsECutM02;                     
  TH2F * fhEOverPvsPCutM02;                     
  
  TH2F * fhdEdxvsECutEOverP;                    
  TH2F * fhdEdxvsPCutEOverP;                    
  TH2F * fhEOverPvsECutM02CutdEdx;              
  TH2F * fhEOverPvsPCutM02CutdEdx;              

  TH2F * fhMCdEdxvsE[10];                       
  TH2F * fhMCdEdxvsP[10];                       
  TH2F * fhMCEOverPvsE[10];                     
  TH2F * fhMCEOverPvsP[10];                     
  
  TH2F * fhNCellsE[2];                          
  TH2F * fhNLME[2];                             
  TH2F * fhMaxCellDiffClusterE[2];              
  TH2F * fhTimeE[2];                            

  TH1F * fhE[2]    ;                            
  TH1F * fhPt[2]   ;                            
  TH2F * fhPhi[2]  ;                            
  TH2F * fhEta[2]  ;                            
  TH2F * fhEtaPhi[2]  ;                         
  TH2F * fhEtaPhi05[2]  ;                       
  
  //Shower shape
  
  TH2F * fhDispE[2];                            
  TH2F * fhLam0E[2];                            
  TH2F * fhLam1E[2];                            

  TH2F * fhDispETRD[2];                         
  TH2F * fhLam0ETRD[2];                         
  TH2F * fhLam1ETRD[2];                         
  
  TH2F * fhNCellsLam0LowE[2];                   
  TH2F * fhNCellsLam0HighE[2];                  

  TH2F * fhEtaLam0LowE[2];                      
  TH2F * fhPhiLam0LowE[2];                      
  TH2F * fhEtaLam0HighE[2];                     
  TH2F * fhPhiLam0HighE[2];                     
    
  TH2F * fhDispEtaE[2] ;                        
  TH2F * fhDispPhiE[2] ;                        
  TH2F * fhSumEtaE[2] ;                         
  TH2F * fhSumPhiE[2] ;                         
  TH2F * fhSumEtaPhiE[2] ;                      
  TH2F * fhDispEtaPhiDiffE[2] ;                 
  TH2F * fhSphericityE[2] ;                     
  TH2F * fhDispEtaDispPhiEBin[2][5] ;           

 //  Weight studies
  
  TH2F * fhECellClusterRatio;                   
  TH2F * fhECellClusterLogRatio;                
  TH2F * fhEMaxCellClusterRatio;                
  TH2F * fhEMaxCellClusterLogRatio;             
  TH2F * fhLambda0ForW0[14];                    
//TH2F * fhLambda1ForW0[14];                    //!<! L1 for 7 defined w0= 3, 3.5 ... 6 for selected electrons
  
  // Fill MC dependent histograms, Origin of this cluster is ...

  TH2F * fhMCDeltaE[2][10]  ;                   
  TH2F * fhMC2E[2][10]  ;                       
  
  TH1F * fhMCE[2][10];                          
  TH1F * fhMCPt[2][10];                         
  TH2F * fhMCPhi[2][10];                        
  TH2F * fhMCEta[2][10];                        
  
  // Shower Shape MC

  TH2F * fhMCELambda0[2][6] ;                   
  
  TH2F * fhMCEDispEta[2][6] ;                   
  TH2F * fhMCEDispPhi[2][6] ;                   
  TH2F * fhMCESumEtaPhi[2][6] ;                 
  TH2F * fhMCEDispEtaPhiDiff[2][6] ;            
  TH2F * fhMCESphericity[2][6] ;                

  TH2F * fhMCElectronELambda0NoOverlap ;        
  TH2F * fhMCElectronELambda0TwoOverlap ;       
  TH2F * fhMCElectronELambda0NOverlap ;         
  
  //Embedding
  TH2F * fhEmbeddedSignalFractionEnergy ;       
  
  TH2F * fhEmbedElectronELambda0FullSignal ;    
  TH2F * fhEmbedElectronELambda0MostlySignal ;  
  TH2F * fhEmbedElectronELambda0MostlyBkg ;     
  TH2F * fhEmbedElectronELambda0FullBkg ;       
  
  AliAnaElectron(              const AliAnaElectron & el) ;
    
  AliAnaElectron & operator = (const AliAnaElectron & el) ;
  
  ClassDef(AliAnaElectron,5) ;

} ;
 

#endif//ALIANAELECTRON_H