AliAnaChargedParticles.h

Go to the documentation of this file.

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

//_________________________________________________________________________
//_________________________________________________________________________

// Root system
class TH2F; 

// Analysis system
#include "AliAnaCaloTrackCorrBaseClass.h"
 
class AliAnaChargedParticles : public AliAnaCaloTrackCorrBaseClass {
  
 public:
    
  AliAnaChargedParticles() ;
    
  virtual ~AliAnaChargedParticles() { ; }
  
  void    FillPrimaryHistograms();
  
  TList * GetCreateOutputObjects();
    
  Bool_t  GetTrackSide  (Float_t trackEta) const { if ( trackEta >=0 ) return kTRUE ; else return kFALSE ; }
  Int_t   GetTrackSector(Float_t trackPhi) const { return Int_t (trackPhi / 0.349065850398865896) ; } // 20*pi/180
  
  void    Init();
  
  void    InitParameters();
  
  void    Print(const Option_t * opt) const;
  
  void    MakeAnalysisFillAOD()  ;
  
  void    MakeAnalysisFillHistograms() ;
  
  void    SwitchOnFillTrackBCHistograms()        { fFillTrackBCHistograms   = kTRUE  ; }
  void    SwitchOffFillTrackBCHistograms()       { fFillTrackBCHistograms   = kFALSE ; }

  void    SwitchOnFillVertexBC0Histograms()      { fFillVertexBC0Histograms = kTRUE  ; }
  void    SwitchOffFillVertexBC0Histograms()     { fFillVertexBC0Histograms = kFALSE ; }

  void    SwitchOnFillEtaPhiRegionsHistograms()  { fFillEtaPhiRegionHistograms = kTRUE  ; }
  void    SwitchOffFillEtaPhiRegionsHistograms() { fFillEtaPhiRegionHistograms = kFALSE ; }

  void    SwitchOnFillTrackMultiplicityHistograms()  { fFillTrackMultHistograms = kTRUE  ; }
  void    SwitchOffFillTrackMultiplicityHistograms() { fFillTrackMultHistograms = kFALSE ; }
  
  void    SwitchOnFillTrackDCAHistograms()       { fFillTrackDCAHistograms = kTRUE  ; }
  void    SwitchOffFillTrackDCAHistograms()      { fFillTrackDCAHistograms = kFALSE ; }
  
 private:
  
  Bool_t  fFillTrackBCHistograms;           
  Bool_t  fFillVertexBC0Histograms;         
  Bool_t  fFillEtaPhiRegionHistograms;      
  Bool_t  fFillTrackMultHistograms;         
  Bool_t  fFillTrackDCAHistograms;          
  
  TLorentzVector fMomentum;                 
  
  //Histograms
  TH2F * fhNTracks;                         
  TH2F * fhSumPtTracks;                     
  
  TH2F * fhPtTrackNTracks    [10];          
  TH2F * fhPtTrackSumPtTracks[10];          
  
  TH1F * fhPt;                              
  TH1F * fhPtNoCut;                         
  TH1F * fhPtCutDCA;                        
  TH1F * fhPtCutDCABCOK;                    

  TH1F * fhPtPerRegion   [18][2];           
  TH1F * fhSumPtPerRegion[18][2];           
  
  TH1F * fhPtNotPrimary;                    
  TH1F * fhPtNotSharedClusterCut;           
  
  TH1F * fhPtPileUp[7];                     
  TH2F * fhPhiNeg;                          
  TH2F * fhEtaNeg;                          
  TH2F * fhPhiPos;                          
  TH2F * fhEtaPos;                          
  TH2F * fhEtaPhiPos;                       
  TH2F * fhEtaPhiNeg;                       
  
  TH2F * fhTrackResolution;                 
  
  TH1F * fhPtVtxOutBC0;                     
  TH2F * fhEtaPhiVtxOutBC0;                 
  TH1F * fhPtVtxInBC0;                      
  TH2F * fhEtaPhiVtxInBC0;                  

  TH1F * fhPtSPDRefit;                      
  TH1F * fhPtNoSPDRefit;                    
  TH1F * fhPtNoSPDNoRefit;                  
  
  TH2F * fhEtaPhiSPDRefit;                  
  TH2F * fhEtaPhiNoSPDRefit;                
  TH2F * fhEtaPhiNoSPDNoRefit;              

  TH2F * fhEtaPhiSPDRefitPt02;              
  TH2F * fhEtaPhiNoSPDRefitPt02;            
  TH2F * fhEtaPhiNoSPDNoRefitPt02;          

  TH2F * fhEtaPhiSPDRefitPt3;               
  TH2F * fhEtaPhiNoSPDRefitPt3;             
  TH2F * fhEtaPhiNoSPDNoRefitPt3;           

  // MC
  enum mvType { kmcPion = 0, kmcProton = 1, kmcKaon = 2, kmcMuon = 3, kmcElectron = 4, kmcUnknown = 5 };

  TH1F * fhPtMCPart [6];                    
  TH2F * fhPhiMCPart[6];                    
  TH2F * fhEtaMCPart[6];                    
  
  TH1F * fhPtMCPrimPart [6];                
  TH2F * fhPhiMCPrimPart[6];                
  TH2F * fhEtaMCPrimPart[6];                

  // TOF and BC
  TH1F * fhTOFSignal;                       
  TH1F * fhTOFSignalPtCut;                  
  TH1F * fhTOFSignalBCOK;                   
  TH2F * fhPtTOFSignal;                     
  TH2F * fhPtTOFSignalDCACut;               
  TH2F * fhPtTOFSignalPileUp[7];            
  TH2F * fhPtTOFSignalVtxOutBC0;            
  TH2F * fhPtTOFSignalVtxOutBC0PileUp[7];   
  TH2F * fhPtTOFSignalVtxInBC0;             
  TH2F * fhPtTOFSignalVtxInBC0PileUp[7];    
  TH1F * fhPtTOFStatus0;                    
  TH2F * fhEtaPhiTOFStatus0;                
  TH2F * fhEtaPhiTOFBC0;                    
  TH2F * fhEtaPhiTOFBCPlus;                 
  TH2F * fhEtaPhiTOFBCMinus;                
  TH2F * fhEtaPhiTOFBC0PileUpSPD;           
  TH2F * fhEtaPhiTOFBCPlusPileUpSPD;        
  TH2F * fhEtaPhiTOFBCMinusPileUpSPD;       
  
  TH1F * fhProductionVertexBC;              
  TH1F * fhProductionVertexBCPileUp[7];     

  // DCA
  TH2F * fhPtDCA[3];                        
  
  TH2F * fhPtDCASPDRefit[3];                
  TH2F * fhPtDCANoSPDRefit[3];              
  TH2F * fhPtDCANoSPDNoRefit[3];            

  TH2F * fhPtDCAVtxOutBC0[3];               
  TH2F * fhPtDCAVtxInBC0[3];                
  TH2F * fhPtDCAPileUp[3];                  
  TH2F * fhPtDCAVtxOutBC0PileUp[3];         
  TH2F * fhPtDCAVtxInBC0PileUp[3];          

  TH2F * fhPtDCATOFBC0[3];                  
  TH2F * fhPtDCAPileUpTOFBC0[3];            
  TH2F * fhPtDCATOFBCOut[3];                

  TH2F * fhPtDCANoTOFHit[3];                
  TH2F * fhPtDCAVtxOutBC0NoTOFHit[3];       
  TH2F * fhPtDCAVtxInBC0NoTOFHit[3];        
  TH2F * fhPtDCAPileUpNoTOFHit[3];          
  TH2F * fhPtDCAVtxOutBC0PileUpNoTOFHit[3]; 
  TH2F * fhPtDCAVtxInBC0PileUpNoTOFHit[3];  
  
  // pile-up
  TH2F * fhPtNPileUpSPDVtx;                 
  TH2F * fhPtNPileUpTrkVtx;                 
  TH2F * fhPtNPileUpSPDVtxBC0;              
  TH2F * fhPtNPileUpTrkVtxBC0;              

  AliAnaChargedParticles(              const AliAnaChargedParticles & ch) ;

  AliAnaChargedParticles & operator = (const AliAnaChargedParticles & ch) ;
  
  ClassDef(AliAnaChargedParticles,12) ;

} ;


#endif //ALIANACHARGEDPARTICLES_H