AliAnalysisTaskEMCALTriggerQA.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEMCALTRIGGERQA_H
#define ALIANALYSISTASKEMCALTRIGGERQA_H

//------------------------------------------------------------------------
//------------------------------------------------------------------------//

//--- Root ---
class TList;
class TH1F;
class TH2I;
class TH2F;
class AliEMCALGeometry;
class TProfile2D;

//--- AliRoot ---
class AliEMCALRecoUtils;
#include "AliEMCALGeoParams.h"
#include "AliAnalysisTaskSE.h"

class AliAnalysisTaskEMCALTriggerQA : public AliAnalysisTaskSE 
{
public:

  AliAnalysisTaskEMCALTriggerQA();                   // default constructor
  
  AliAnalysisTaskEMCALTriggerQA(const char *name);   // named constructor
  
  virtual ~AliAnalysisTaskEMCALTriggerQA() { ; }     // destructor
  
  void   ClusterAnalysis();
  
  void   FillCellMaps();
  
  void   FillTriggerPatchMaps(TString triggerclasses);
  
  void   FillClusterHistograms(Int_t triggerNumber, Bool_t maxCluster,
                               Float_t e,Float_t eta,Float_t phi,
                               Float_t ietamax,Float_t iphimax, Int_t sm,
                               Float_t centrality, Float_t v0AC);
  
  void   FillCorrelationHistograms();
  
  void   FillEventCounterHistogram();
  
  void   FillL1GammaPatchHistograms();
  
  void   FillL1JetPatchHistograms();
  
  void   FillMapHistograms();
  
  void   FillV0Histograms();
  
  void   Init() ;

  void   InitHistogramArrays() ;

  void   InitCellPatchMaps();
  
  void   LocalInit()                     { Init()                       ; }

  void   UserCreateOutputObjects();    
  
  void   UserExec(Option_t *option);   
  
  AliEMCALRecoUtils* GetRecoUtils()      { if(!fRecoUtils) fRecoUtils = new AliEMCALRecoUtils ;
                                           return fRecoUtils            ; }
  
  void   SetEtaPhiEnMin(Float_t en)      { fEtaPhiEnMin       = en      ; }

  void   SetTriggerEventBit(TString list) ;
  
  // OADB and geometry settings
  
  void   InitGeometry();
  
  void   SetGeometryName(TString name)   { fGeoName           = name    ; }
    
  void   SetEventTriggerL1Bit(Int_t ega, Int_t eje)
                                         { fBitEGA   = ega ; fBitEJE = eje; }

  void   AccessOADB() ;
  
  void   SwitchOnEMCALOADB()             { fAccessOADB        = kTRUE   ; }
  void   SwitchOffEMCALOADB()            { fAccessOADB        = kFALSE  ; }

  void   SwitchOnMCData()                { fMCData            = kTRUE   ; }
  void   SwitchOffMCData()               { fMCData            = kFALSE  ; }

  void   SwitchOnV0SignalHistograms()    { fFillV0SigHisto    = kTRUE   ; }
  void   SwitchOffV0SignalHistograms()   { fFillV0SigHisto    = kFALSE  ; }

  void   SwitchOnCentralityHistograms()  { fFillCenHisto      = kTRUE   ; }
  void   SwitchOffCentralityHistograms() { fFillCenHisto      = kFALSE  ; }
  
  void   SwitchOnAliCentrality ()        { fUseAliCentrality  = kTRUE  ; }
  void   SwitchOffAliCentrality()        { fUseAliCentrality  = kFALSE ; }
  
  void   SwitchOnClusterAcceptanceHistograms()  { fFillClusAcceptHisto = kTRUE   ; }
  void   SwitchOffClusterAcceptanceHistograms() { fFillClusAcceptHisto = kFALSE  ; }
  
  void   SetSuperModulesRange(Int_t min, Int_t max) { fFirstSM = min ; if(max < 20) fLastSM = max ; else fLastSM = 19 ; }
  
  void   SetCentralityEstimator(TString est) { fCentEstimator  = est    ; }
  
  void   SetOADBFilePath(TString path)       { fOADBFilePath   = path   ; }
  
  // Histogram setters

  void   SetTRUTotalSignalHistogramsRange(Int_t nbins,  Float_t max) { fNBinsTRUSignal   = nbins; fMaxTRUSignal   = max ; }
  void   SetSTUTotalSignalHistogramsRange(Int_t nbins,  Float_t max) { fNBinsSTUSignal   = nbins; fMaxSTUSignal   = max ; }
  void   SetV0TotalSignalHistogramsRange (Int_t nbins,  Float_t max) { fNBinsV0Signal    = nbins; fMaxV0Signal    = max ; }
  void   SetSTUFEERatioHistogramsRange   (Int_t nbins,  Float_t max) { fNBinsSTUFEERatio = nbins; fMaxSTUFEERatio = max ; }
  void   SetSTUTRURatioHistogramsRange   (Int_t nbins,  Float_t max) { fNBinsSTUTRURatio = nbins; fMaxSTUFEERatio = max ; }
  void   SetClusterEHistogramsRange      (Int_t nbins,  Float_t max) { fNBinsClusterE    = nbins; fMaxClusterE    = max ; }
  void   SetClusterEtaHistogramsRange    (Int_t nbins,  Float_t max) { fNBinsClusterEta  = nbins; fMaxClusterEta  = max ; }
  void   SetClusterPhiHistogramsRange    (Int_t nbins,  Float_t max, Float_t min)
  { fNBinsClusterPhi  = nbins; fMaxClusterPhi  = max ;  fMinClusterPhi = min ; }
  
private:
    
  TList            *fOutputList;      
  
  AliEMCALRecoUtils *fRecoUtils;      

  Bool_t            fGeoSet  ;        
  AliEMCALGeometry *fGeometry;        
  TString           fGeoName;         
  
  Bool_t            fOADBSet ;        
  Bool_t            fAccessOADB ;     
  TString           fOADBFilePath ;   
  
  Int_t             fBitEGA;          
  Int_t             fBitEJE;          
  
  Float_t           fEtaPhiEnMin;     
  
  Int_t             fSTUTotal;        
  Float_t           fTRUTotal;        
  Float_t           fV0Trigger;       
  Float_t           fV0A;             
  Float_t           fV0C;             
  
  Bool_t            fFillV0SigHisto;  
  Bool_t            fFillCenHisto;    
  Bool_t            fFillClusAcceptHisto; 
  Bool_t            fMCData;          
  Int_t             fFirstSM;         
  Int_t             fLastSM;          
  
  TString           fCentEstimator;   
  Bool_t            fUseAliCentrality;

  // Event by event trigger recognition bit
  Bool_t            fEventMB   ;      
  Bool_t            fEventL0   ;      
  Bool_t            fEventL0D  ;      
  Bool_t            fEventL1G  ;      
  Bool_t            fEventL1GD  ;     
  Bool_t            fEventL1G2 ;      
  Bool_t            fEventL1G2D ;     
  Bool_t            fEventL1J  ;      
  Bool_t            fEventL1JD  ;     
  Bool_t            fEventL1J2 ;      
  Bool_t            fEventL1J2D ;     

  Bool_t            fEventCen  ;      
  Bool_t            fEventSem  ;      
  
  TLorentzVector    fMomentum ;       
    
  // Histograms
  
  TH1F             *fhNEvents;                
  TH2F             *fhFORAmp;                 
  TH2F             *fhFORAmpL1G;              
  TH2F             *fhFORAmpL1G2;             
  TH2F             *fhFORAmpL1J;              
  TH2F             *fhFORAmpL1J2;             
  TH2F             *fhL0Amp;                  
  TH2F             *fhL0AmpL1G;               
  TH2F             *fhL0AmpL1J;               
  TH2F             *fhL1Amp;                  
  TH2F             *fhL1GAmp;                 
  TH2F             *fhL1G2Amp;                
  TH2F             *fhL1JAmp;                 
  TH2F             *fhL1J2Amp;                
  TH2F             *fhL1FOREnergy;            
  
  TH2F             *fhL0Patch;                
  TH2F             *fhL1GPatch;               
  TH2F             *fhL1G2Patch;              
  TH2F             *fhL1GPatchNotFake;        
  TH2F             *fhL1GPatchFake;           
  TH2F             *fhL1GPatchNotAllFake;     
  TH2F             *fhL1GPatchAllFake;        
  TH2F             *fhL1GPatchNotAllFakeMax;  
  TH2F             *fhL1GPatchAllFakeMax;     
  TH1F             *fhL1GPatchNotAllFakeMaxE; 
  TH1F             *fhL1GPatchAllFakeMaxE;    
  TH1F             *fhL1GPatchNotAllFakeE;    
  TH1F             *fhL1GPatchAllFakeE;       
  TH1F             *fhL1GPatchFakeE;          
  TH1F             *fhL1GPatchNotFakeE;       
  TH2F             *fhNPatchFake;             
  TH2F             *fhNPatchNotFake;          
  
  TH2F             *fhL1JPatch;               
  TH2F             *fhL1J2Patch;              
  TH2F             *fhFEESTU;                 
  TH2F             *fhTRUSTU;                 
  TH2I             *fhV0STU;                  
  
  TH2F             *fhGPMaxVV0TT;             
  TH2F             *fhJPMaxVV0TT;             
  TProfile2D       *fhFORMeanAmp;             
  TProfile2D       *fhL0MeanAmp;              
  TProfile2D       *fhL1MeanAmp;              
  TH2F             *fhL1GPatchMax;            
  TH2F             *fhL1G2PatchMax;           
  TH2F             *fhL1JPatchMax;            
  TH2F             *fhL1J2PatchMax;           
  
  enum triggerType{ kMBTrig                = 0,  
                    kL0Trig                = 1,  kL0TrigD                = 2,
                    kL1GammaTrig           = 3,  kL1GammaTrigD           = 4, 
                    kL1GammaTrig2          = 5,  kL1GammaTrig2D          = 6,
                    kL1JetTrig             = 7,  kL1JetTrigD             = 8,
                    kL1JetTrig2            = 9,  kL1JetTrig2D            = 10,
                    kL1GammaOnlyTrig       = 11, kL1GammaOnlyTrigD       = 12,  
                    kL1JetOnlyTrig         = 13, kL1JetOnlyTrigD         = 14,
                    kL1Gamma2OnlyGammaTrig = 15, kL1Gamma2OnlyGammaTrigD = 16,  
                    kL1Jet2OnlyJetTrig     = 17, kL1Jet2OnlyJetTrigD     = 18,
                    kL0TrigPureEMC         = 19, kL0TrigPureDMC          = 20,
                    kL1GTrigPureEMC        = 21, kL1GTrigPureDMC         = 22,
                    kL1JTrigPureEMC        = 23, kL1JTrigPureDMC         = 24,
                    kCentralTrig           = 25, kSemiCentralTrig        = 26 };
  
  TH1F             *fhClusMBPure[3];                                 
  TH1F             *fhClusMaxMBPure[3];                              
  
  static const int  fgkTriggerCombi = 27;                            
  
  TH1F             *fhClus   [fgkTriggerCombi];                      
  TH1F             *fhClusMax[fgkTriggerCombi];                      

  TH2F             *fhClusCen   [fgkTriggerCombi];                   
  TH2F             *fhClusCenMax[fgkTriggerCombi];                   
  
  TH2F             *fhClusV0   [fgkTriggerCombi];                    
  TH2F             *fhClusV0Max[fgkTriggerCombi];                    

  TH1F             *fhClusSM   [fgkTriggerCombi][20];                
  TH2F             *fhClusCenSM[fgkTriggerCombi][20];                
  TH2F             *fhClusV0SM [fgkTriggerCombi][20];                
  
  TH2F             *fhClusEta   [fgkTriggerCombi];                   
  TH2F             *fhClusEtaMax[fgkTriggerCombi];                   

  TH2F             *fhClusPhi   [fgkTriggerCombi];                   
  TH2F             *fhClusPhiMax[fgkTriggerCombi];                   

  TH2F             *fhClusEtaPhiHigh      [fgkTriggerCombi];         
  TH2F             *fhClusEtaPhiHighCluMax[fgkTriggerCombi];         
 
  TH2F             *fhClusEtaPhiHighCellMax      [fgkTriggerCombi];  
  TH2F             *fhClusEtaPhiHighCellMaxCluMax[fgkTriggerCombi];  
 
  TH2F             *fhClusEtaPhiLow      [fgkTriggerCombi];          
  TH2F             *fhClusEtaPhiLowCluMax[fgkTriggerCombi];          
 
  TH2F             *fhClusEtaPhiLowCellMax      [fgkTriggerCombi];   
  TH2F             *fhClusEtaPhiLowCellMaxCluMax[fgkTriggerCombi];   

  TH1F             *fhV0[fgkTriggerCombi];                           
  TH1F             *fhCentrality[fgkTriggerCombi];                   

  // Histograms bins
  
  Int_t             fNBinsSTUSignal   ;     // Number of bins for STU total signal histograms
  Float_t           fMaxSTUSignal     ;     // Maximum value for TRU total signal histograms
  Int_t             fNBinsTRUSignal   ;     // Number of bins for TRU total signal histograms
  Float_t           fMaxTRUSignal     ;     // Maximum value for TRU total signal histograms
  Int_t             fNBinsV0Signal    ;     // Number of bins for V0 total signal histograms
  Float_t           fMaxV0Signal      ;     // Maximum value for V0 total signal histograms
  Int_t             fNBinsSTUFEERatio ;     // Number of bins for STU/FEE ratios histograms
  Float_t           fMaxSTUFEERatio   ;     // Maximum value for STU/FEE ratios histograms
  Int_t             fNBinsSTUTRURatio ;     // Number of bins for STU/TRU ratios histograms
  Float_t           fMaxSTUTRURatio   ;     // Maximum value for STU/TRU ratios histograms
  Int_t             fNBinsClusterE    ;     // Number of bins for E cluster histograms
  Float_t           fMaxClusterE      ;     // Maximum value for E cluster histograms
  Int_t             fNBinsClusterPhi  ;     // Number of bins for Phi cluster histograms
  Float_t           fMaxClusterPhi    ;     // Maximum value for Phi cluster histograms
  Float_t           fMinClusterPhi    ;     // Maximum value for Phi cluster histograms
  Int_t             fNBinsClusterEta  ;     // Number of bins for Eta cluster histograms
  Float_t           fMaxClusterEta    ;     // Maximum value for Eta cluster histograms
    
  //static const int  fgkFALTRORows = AliEMCALGeoParams::fgkEMCALRows*(AliEMCALGeoParams::fgkEMCALModules-7)/2;   // total number
    
  static const int  fgkFALTRORows = 104; // 60 //AliEMCALGeoParams::fgkEMCALSTURows-4;
  
  static const int  fgkFALTROCols = AliEMCALGeoParams::fgkEMCALSTUCols;
  
  // cell, patch maps
  Double_t fMapCell     [fgkFALTRORows][fgkFALTROCols]; // Cell map
  Double_t fMapCellL1G  [fgkFALTRORows][fgkFALTROCols]; // Cell map for L1G
  Double_t fMapCellL1G2 [fgkFALTRORows][fgkFALTROCols]; // Cell map for L1G2
  Double_t fMapCellL1J  [fgkFALTRORows][fgkFALTROCols]; // Cell map for L1J
  Double_t fMapCellL1J2 [fgkFALTRORows][fgkFALTROCols]; // Cell map for L1J2
  Double_t fMapTrigL0   [fgkFALTRORows][fgkFALTROCols]; // Patch map for L0
  Double_t fMapTrigL1   [fgkFALTRORows][fgkFALTROCols]; // Patch map for L1
  Double_t fMapTrigL0L1G[fgkFALTRORows][fgkFALTROCols]; // Patch map for L0L1G
  Double_t fMapTrigL0L1J[fgkFALTRORows][fgkFALTROCols]; // Patch map for L0L1J
  Double_t fMapTrigL1G  [fgkFALTRORows][fgkFALTROCols]; // Patch map for L1G
  Double_t fMapTrigL1G2 [fgkFALTRORows][fgkFALTROCols]; // Patch map for L1G2
  Double_t fMapTrigL1J  [fgkFALTRORows][fgkFALTROCols]; // Patch map for L1J
  Double_t fMapTrigL1J2 [fgkFALTRORows][fgkFALTROCols]; // Patch map for L1J2

  AliAnalysisTaskEMCALTriggerQA           (const AliAnalysisTaskEMCALTriggerQA&) ;
  
  AliAnalysisTaskEMCALTriggerQA& operator=(const AliAnalysisTaskEMCALTriggerQA&) ; 
  
  ClassDef(AliAnalysisTaskEMCALTriggerQA, 18) ;

};

#endif