AliAnaCalorimeterQA.h

Go to the documentation of this file.

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

//_________________________________________________________________________

// --- Root system ---
class TH3F;
class TH2F;
class TH1F;
class TObjString;
class TObjArray;

// --- Analysis system --- 
class AliVCaloCells;
class AliVCaloCluster;
class AliVTrack;

#include "AliAnaCaloTrackCorrBaseClass.h"
 
class AliAnaCalorimeterQA : public AliAnaCaloTrackCorrBaseClass {
  
public:
    
  AliAnaCalorimeterQA() ;
    
  virtual ~AliAnaCalorimeterQA() { ; }
    
  // General methods
  
  TObjString * GetAnalysisCuts();

  TList      * GetCreateOutputObjects();
  
  void         Init();
  
  void         InitParameters();
    
  void         MakeAnalysisFillHistograms() ;
  
  void         Print(const Option_t * opt) const;
    
  // Main methods
  
  void         BadClusterHistograms(AliVCluster* clus, const TObjArray *caloClusters,  AliVCaloCells * cells,
                                    Int_t absIdMax,    Double_t maxCellFraction, Float_t eCrossFrac, Double_t tmax);
    
  void         CalculateAverageTime(AliVCluster *clus, AliVCaloCells *cells, Double_t timeAverages[2]);
  
  void         CellHistograms(AliVCaloCells * cells);

  void         CellInClusterPositionHistograms(AliVCluster* cluster);
    
  void         ClusterAsymmetryHistograms(AliVCluster* clus, Int_t absIdMax, Bool_t goodCluster );
  
  void         ClusterHistograms(AliVCluster* cluster, const TObjArray *caloClusters,  AliVCaloCells * cells, 
                                 Int_t absIdMax, Double_t maxCellFraction, Float_t eCrossFrac, Double_t tmax);
  
  void         ClusterLoopHistograms(const TObjArray * clusters, AliVCaloCells * cells);
  
  Bool_t       ClusterMCHistograms(Bool_t matched, const Int_t * labels, Int_t nLabels, Int_t & pdg );

  void         ClusterMatchedWithTrackHistograms(AliVCluster* clus, Bool_t mcOK, Int_t pdg);

  void         Correlate();
  
  void         ExoticHistograms(Int_t absIdMax, Float_t ampMax,
                                AliVCluster *clus, AliVCaloCells* cells);
  
  Float_t      GetECross(Int_t absId, AliVCaloCells* cells,Float_t dtcut = 10000);
  
  void         InvariantMassHistograms(Int_t iclus, Int_t nModule, const TObjArray* caloClusters, AliVCaloCells * cells);

  Bool_t       IsGoodCluster(Int_t absIdMax, Float_t m02, Int_t nCellsPerCluster, AliVCaloCells *cells);

  void         MCHistograms();  
  
  void         WeightHistograms(AliVCluster *clus, AliVCaloCells* cells);

  
  // Setters and getters
  
  Float_t      GetEMCALCellAmpMin()      const  { return fEMCALCellAmpMin    ; }
  void         SetEMCALCellAmpMin(Float_t amp)  { fEMCALCellAmpMin = amp     ; }
  
  Float_t      GetPHOSCellAmpMin()       const  { return fPHOSCellAmpMin     ; }
  void         SetPHOSCellAmpMin (Float_t amp)  { fPHOSCellAmpMin  = amp     ; }

  Int_t        GetEMCALNCellsPerClusterMin() const  { return fEMCALClusterNCellMin ; }
  void         SetEMCALNCellsPerClusterMin(Int_t n) { fEMCALClusterNCellMin = n    ; }
  
  Int_t        GetPHOSNCellsPerClusterMin()  const  { return fPHOSClusterNCellMin  ; }
  void         SetPHOSNCellsPerClusterMin (Int_t n) { fPHOSClusterNCellMin  = n    ; }
  
  Float_t      GetInvMassMinECut()       const  { return fInvMassMinECut     ; }
  void         SetInvMassMinECut(Float_t cut)   { fInvMassMinECut = cut      ; }

  Float_t      GetInvMassMaxECut()       const  { return fInvMassMaxECut     ; }
  void         SetInvMassMaxECut(Float_t cut)   { fInvMassMaxECut = cut      ; }

  Float_t      GetInvMassMinM02Cut()     const  { return fInvMassMinM02Cut   ; }
  void         SetInvMassMinM02Cut(Float_t cut) { fInvMassMinM02Cut = cut    ; }
  
  Float_t      GetInvMassMaxM02Cut()     const  { return fInvMassMaxM02Cut   ; }
  void         SetInvMassMaxM02Cut(Float_t cut) { fInvMassMaxM02Cut = cut    ; }
  
  Float_t      GetInvMassMaxOpenAngle()  const  { return fInvMassMaxOpenAngle; }
  void         SetInvMassMaxOpenAngle(Float_t c){ fInvMassMaxOpenAngle = c   ; }
  
  Float_t      GetInvMassMaxTimeDifference()  const   { return fInvMassMaxTimeDifference; }
  void         SetInvMassMaxTimeDifference(Float_t c) { fInvMassMaxTimeDifference = c   ; }  
  
  Double_t     GetTimeCutMin()           const  { return fTimeCutMin         ; }
  Double_t     GetTimeCutMax()           const  { return fTimeCutMax         ; }
  void         SetTimeCut(Double_t min, Double_t max) {
                          fTimeCutMin = min ; fTimeCutMax = max              ; }
    
  // Histogram switchs
  
  void SwitchOnFillAllCellTimeHisto()           { fFillAllCellTimeHisto = kTRUE  ; }
  void SwitchOffFillAllCellTimeHisto()          { fFillAllCellTimeHisto = kFALSE ; }
  
  void SwitchOnFillAllPositionHistogram()       { fFillAllPosHisto  = kTRUE  ; }
  void SwitchOffFillAllPositionHistogram()      { fFillAllPosHisto  = kFALSE ; }
  
  void SwitchOnFillAllPositionHistogram2()      { fFillAllPosHisto2 = kTRUE  ; }
  void SwitchOffFillAllPositionHistogram2()     { fFillAllPosHisto2 = kFALSE ; }
  
  void SwitchOnFillAllTH3Histogram()            { fFillAllTH3       = kTRUE  ; }
  void SwitchOffFillAllTH3Histogram()           { fFillAllTH3       = kFALSE ; }
  
  void SwitchOnFillAllTrackMatchingHistogram()  { fFillAllTMHisto   = kTRUE  ; }
  void SwitchOffFillAllTrackMatchingHistogram() { fFillAllTMHisto   = kFALSE ; }
  
  void SwitchOnFillAllPi0Histogram()            { fFillAllPi0Histo  = kTRUE  ; }
  void SwitchOffFillAllPi0Histogram()           { fFillAllPi0Histo  = kFALSE ; }

  void SwitchOnFillInvMassOpAngleHistogram()    { fFillInvMassOpenAngle = kTRUE  ; }
  void SwitchOffFillInvMassOpAngleHistogram()   { fFillInvMassOpenAngle = kFALSE ; }  

  void SwitchOnFillPi0PairDiffTimeHistogram()   { fFillPi0PairDiffTime = kTRUE  ; }
  void SwitchOffFillPi0PairDiffTimeHistogram()  { fFillPi0PairDiffTime = kFALSE ; }

  void SwitchOnFillInvMassReducedEMCALHistogram()  { fFillInvMassInEMCALWithPHOSDCalAcc = kTRUE  ; }
  void SwitchOffFillInvMassReducedEMCALHistogram() { fFillInvMassInEMCALWithPHOSDCalAcc = kFALSE ; }

  void SwitchOnCorrelation()                    { fCorrelate        = kTRUE  ; }
  void SwitchOffCorrelation()                   { fCorrelate        = kFALSE ; }

  void SwitchOnStudyBadClusters()               { fStudyBadClusters = kTRUE  ; }
  void SwitchOffStudyBadClusters()              { fStudyBadClusters = kFALSE ; }
  
  // Analysis not to be used in QA
  
  void SwitchOnStudyClustersAsymmetry()         { fStudyClustersAsymmetry = kTRUE  ; }
  void SwitchOffStudyClustersAsymmetry()        { fStudyClustersAsymmetry = kFALSE ; }

  void SwitchOnStudyWeight()                    { fStudyWeight      = kTRUE  ; }
  void SwitchOffStudyWeight()                   { fStudyWeight      = kFALSE ; }
  
  void SwitchOnStudyExotic()                    { fStudyExotic      = kTRUE  ; }
  void SwitchOffStudyExotic()                   { fStudyExotic      = kFALSE ; }
  
  void SetNECrossCuts(Int_t n)                  { fExoNECrossCuts   = n      ; }
  void SetNDTimeCuts (Int_t n)                  { fExoNDTimeCuts    = n      ; }
  
  void SetExoECrossCuts (Int_t i, Float_t c)    { if (i<fExoNECrossCuts) fExoECrossCuts[i] = c ; }
  void SetExoDTimeCuts  (Int_t i, Float_t c)    { if (i<fExoNDTimeCuts ) fExoDTimeCuts [i] = c ; }
  
  void SetConstantTimeShift(Float_t shift)      { fConstantTimeShift     = shift  ; }

  
 private:
  
  // Switches
    
  Bool_t   fFillAllCellTimeHisto;               
  Bool_t   fFillAllPosHisto;                    
  Bool_t   fFillAllPosHisto2;                   
  Bool_t   fFillAllTH3 ;                        
  Bool_t   fFillAllTMHisto ;                    
  
  Bool_t   fFillAllPi0Histo ;                   
  Bool_t   fFillInvMassOpenAngle;               
  Bool_t   fFillPi0PairDiffTime;                
  Bool_t   fFillInvMassInEMCALWithPHOSDCalAcc;  
  
  Bool_t   fCorrelate   ;                       
  Bool_t   fStudyBadClusters;                   
  
  // Analysis not to be used in QA
  
  Bool_t   fStudyClustersAsymmetry;             
  Bool_t   fStudyExotic;                        
  Bool_t   fStudyWeight;                        
  
  // Parameters
    
  Int_t    fNModules    ;                       
  Int_t    fNRCU        ;                       
  Int_t    fNMaxCols    ;                       
  Int_t    fNMaxRows    ;                       
  
  // Cuts
  
  Double_t fTimeCutMin  ;                       
  Double_t fTimeCutMax  ;                       
  Float_t  fCellAmpMin;                         
  Float_t  fEMCALCellAmpMin;                    
  Float_t  fPHOSCellAmpMin ;                    
  Int_t    fEMCALClusterNCellMin;               
  Int_t    fPHOSClusterNCellMin ;               

  // Invariant mass analysis
  
  Float_t  fInvMassMinECut;                     
  Float_t  fInvMassMaxECut;                     
  Float_t  fInvMassMinM02Cut;                   
  Float_t  fInvMassMaxM02Cut;                   
  Float_t  fInvMassMaxOpenAngle;                
  Float_t  fInvMassMaxTimeDifference;           
  // Exotic studies
    
  Float_t  fExoNECrossCuts   ;                  
  Float_t  fExoECrossCuts[10];                  
  Float_t  fExoNDTimeCuts    ;                  
  Float_t  fExoDTimeCuts[5]  ;                  
  
  TLorentzVector fClusterMomentum;              
  TLorentzVector fClusterMomentum2;             
  TLorentzVector fPrimaryMomentum;              
  
  Float_t  fConstantTimeShift;                  

  
  // Calorimeter Clusters
    
  TH1F *   fhE  ;                               
  TH1F *   fhPt ;                               
  TH1F *   fhPhi;                               
  TH1F *   fhEta;                               
  TH2F *   fhEtaPhi;                            
  TH3F *   fhEtaPhiE  ;                         
  TH1F *   fhECharged  ;                        
  TH1F *   fhPtCharged ;                        
  TH1F *   fhPhiCharged;                        
  TH1F *   fhEtaCharged;                        
  TH2F *   fhEtaPhiCharged;                     
  TH3F *   fhEtaPhiECharged;                    
    
  TH2F *   fhIM;                                
  TH2F *   fhIMSame;                            
  TH2F *   fhIMDiff;                            
  TH2F *   fhIMDCAL;                            
  TH2F *   fhIMDCALSame;                        
  TH2F *   fhIMDCALDiff;                        
  TH2F *   fhIMDCALPHOS;                        
  TH2F *   fhIMDCALPHOSSame;                    
  TH2F *   fhIMEMCALPHOS;                       
  TH2F *   fhIMEMCALPHOSSame;                   
  
  TH2F *   fhAsym;                              
  
  TH2F*    fhOpAngle;                           
  TH2F*    fhIMvsOpAngle;                       
  
  TH2F *   fhClusterPairDiffTimeEPi0Mass;       
  TH2F *   fhClusterPairDiffTimeEPi0MassSame;   
  
  TH2F *   fhClusterPairDiffTimeEPi0MassDCal;     
  TH2F *   fhClusterPairDiffTimeEPi0MassDCalSame; 
  
  TH2F *   fhNCellsPerCluster;                  
  TH2F *   fhNCellsPerClusterNoCut;             

  TH1F *   fhNClusters;                         

  TH2F *   fhClusterTimeEnergy;                 
  TH2F *   fhCellTimeSpreadRespectToCellMax;    
  TH1F *   fhCellIdCellLargeTimeSpread;         
  TH2F *   fhClusterPairDiffTimeE;              
  TH2F *   fhClusterPairDiffTimeESameMod;       
  
  TH2F *   fhClusterMaxCellCloseCellRatio;      
  TH2F *   fhClusterMaxCellCloseCellDiff;       
  TH2F *   fhClusterMaxCellDiff;                
  TH2F *   fhClusterMaxCellDiffNoCut;           
  
//  TH2F *   fhClusterMaxCellDiffAverageTime;     //!<! Difference between cluster average time and time of cell with more energy
//  TH2F *   fhClusterMaxCellDiffWeightedTime;    //!<! Difference between cluster weighted time and time of cell with more energy
  TH2F *   fhClusterMaxCellECross;              
  
  TH2F *   fhLambda0;                           
  TH2F *   fhLambda1;                           
//  TH2F *   fhDispersion;                        //!<! Cluster Dispersion vs Energy
  
  // Bad clusters histograms
  
  TH1F *   fhBadClusterEnergy;                  
  TH2F *   fhBadClusterTimeEnergy;              
  TH2F *   fhBadClusterEtaPhi;                  
  TH2F *   fhBadClusterPairDiffTimeE;           
  TH2F *   fhBadCellTimeSpreadRespectToCellMax; 
  
  TH2F *   fhBadClusterMaxCellCloseCellRatio;   
  TH2F *   fhBadClusterMaxCellCloseCellDiff ;   
  TH2F *   fhBadClusterMaxCellDiff;             
  
  TH2F *   fhBadClusterMaxCellDiffAverageTime;  
  TH2F *   fhBadClusterMaxCellDiffWeightedTime; 
  TH2F *   fhBadClusterMaxCellECross;           

  TH2F *   fhBadClusterLambda0;                 
  TH2F *   fhBadClusterLambda1;                 

  // Cluster cell size
    
  TH2F *   fhDeltaIEtaDeltaIPhiE0[2];           
  TH2F *   fhDeltaIEtaDeltaIPhiE2[2];           
  TH2F *   fhDeltaIEtaDeltaIPhiE6[2];           
  TH2F *   fhDeltaIA[2];                        
  TH2F *   fhDeltaIAL0[2];                      
  TH2F *   fhDeltaIAL1[2];                      
  TH2F *   fhDeltaIANCells[2] ;                 
  TH2F *   fhDeltaIAMC[4];                      
  TH2F *   fhBadClusterDeltaIEtaDeltaIPhiE0;    
  TH2F *   fhBadClusterDeltaIEtaDeltaIPhiE2;    
  TH2F *   fhBadClusterDeltaIEtaDeltaIPhiE6;    
  TH2F *   fhBadClusterDeltaIA;                 
  
  // Cluster/cell Position
    
  TH2F *   fhRNCells ;                          
  TH2F *   fhXNCells ;                          
  TH2F *   fhYNCells ;                          
  TH2F *   fhZNCells ;                          
  
  TH2F *   fhRE ;                               
  TH2F *   fhXE ;                               
  TH2F *   fhYE ;                               
  TH2F *   fhZE ;                               
  TH3F *   fhXYZ;                               
  
  TH2F *   fhRCellE ;                           
  TH2F *   fhXCellE ;                           
  TH2F *   fhYCellE ;                           
  TH2F *   fhZCellE ;                           
  TH3F *   fhXYZCell;                           
  
  TH2F *   fhDeltaCellClusterRNCells ;          
  TH2F *   fhDeltaCellClusterXNCells ;          
  TH2F *   fhDeltaCellClusterYNCells ;          
  TH2F *   fhDeltaCellClusterZNCells ;          
  
  TH2F *   fhDeltaCellClusterRE ;               
  TH2F *   fhDeltaCellClusterXE ;               
  TH2F *   fhDeltaCellClusterYE ;               
  TH2F *   fhDeltaCellClusterZE ;               
  
  // Calorimeter cells
    
  TH1F *   fhNCells;                            
  TH1F *   fhNCellsCutAmpMin;                   
  TH1F *   fhAmplitude;                         
  TH2F *   fhAmpId;                             
  TH3F *   fhEtaPhiAmpCell;                     
  TH2F *   fhEtaPhiCell;                        
   
  TH1F *   fhTime;                              
 //TH2F *   fhTimeVz;                            //!<! Time measured in towers/crystals vs vertex z component, for E > 0.5
  TH2F *   fhTimeId;                            
  TH2F *   fhTimeAmp;                           
  
  TH2F *   fhAmpIdLowGain;                      
  TH2F *   fhTimeIdLowGain;                     
  TH2F *   fhTimeAmpLowGain;                    

  TH2F *   fhCellECross;                        
  
  // Calorimeters Correlation
  
  TH2F *   fhEMCALPHOSCorrNClusters;            
  TH2F *   fhEMCALPHOSCorrEClusters;            
  TH2F *   fhEMCALPHOSCorrNCells;               
  TH2F *   fhEMCALPHOSCorrECells;               

  TH2F *   fhEMCALDCALCorrNClusters;            
  TH2F *   fhEMCALDCALCorrEClusters;            
  TH2F *   fhEMCALDCALCorrNCells;               
  TH2F *   fhEMCALDCALCorrECells;               

  TH2F *   fhDCALPHOSCorrNClusters;             
  TH2F *   fhDCALPHOSCorrEClusters;             
  TH2F *   fhDCALPHOSCorrNCells;                
  TH2F *   fhDCALPHOSCorrECells;                
  
  // V0 Correlation
  
  TH2F *   fhCaloV0SCorrNClusters;              
  TH2F *   fhCaloV0SCorrEClusters;              
  TH2F *   fhCaloV0SCorrNCells;                 
  TH2F *   fhCaloV0SCorrECells;                 
  TH2F *   fhCaloV0MCorrNClusters;              
  TH2F *   fhCaloV0MCorrEClusters;              
  TH2F *   fhCaloV0MCorrNCells;                 
  TH2F *   fhCaloV0MCorrECells;                 
  
  // Track Correlation
  
  TH2F *   fhCaloTrackMCorrNClusters;           
  TH2F *   fhCaloTrackMCorrEClusters;           
  TH2F *   fhCaloTrackMCorrNCells;              
  TH2F *   fhCaloTrackMCorrECells;              
  
  // Centrality
  
  TH2F *   fhCaloCenNClusters;                  
  TH2F *   fhCaloCenEClusters;                  
  TH2F *   fhCaloCenNCells;                     
  TH2F *   fhCaloCenECells;                     

  // Event plane
  
  TH2F *   fhCaloEvPNClusters;                  
  TH2F *   fhCaloEvPEClusters;                  
  TH2F *   fhCaloEvPNCells;                     
  TH2F *   fhCaloEvPECells;                     
  
  // Module histograms
  
  TH2F *   fhEMod  ;                            
  TH2F *   fhAmpMod ;                           
  TH2F *   fhEWeirdMod  ;                       
  TH2F *   fhAmpWeirdMod ;                      
  TH2F *   fhTimeMod ;                          
  TH2F *   fhNClustersMod ;                     
  TH2F *   fhNCellsMod ;                        
  TH2F *   fhSumClustersEnergyMod ;             
  TH2F *   fhSumCellsAmpMod ;                   
  TH2F **  fhNCellsSumAmpPerMod ;               
  TH2F **  fhNClustersSumEnergyPerMod ;         
  TH2F **  fhNCellsPerClusterMod ;              
  TH2F **  fhNCellsPerClusterModNoCut ;         
  TH2F *   fhNCellsPerClusterWeirdMod ;         
  TH2F *   fhNCellsPerClusterWeirdModNoCut ;    
  TH2F *   fhGridCells ;                        
  TH2F *   fhGridCellsE ;                       
  TH2F *   fhGridCellsTime ;                    
  TH2F *   fhGridCellsLowGain ;                 
  TH2F *   fhGridCellsELowGain ;                
  TH2F *   fhGridCellsTimeLowGain ;             
  TH2F **  fhTimeAmpPerRCU;                     
  TH2F **  fhIMMod;                             
  
  // Weight studies
  
  TH2F *   fhECellClusterRatio;                 
  TH2F *   fhECellClusterLogRatio;              
  TH2F *   fhEMaxCellClusterRatio;              
  TH2F *   fhEMaxCellClusterLogRatio;           
  
  TH2F *   fhLambda0ForW0AndCellCuts    [12][4][3]; 
  TH2F *   fhLambda0ForW0AndCellCutsEta0[12][4][3]; 
//TH2F *   fhLambda1ForW0AndCellCuts    [12][4][3]; //!<! L1 for different w0 and cell cuts

  TH2F *   fhLambda0ForW0MC[12][5];             
//TH2F *   fhLambda1ForW0MC[12][5];             //!<! L1 for different w0, depending on the particle of origin
  
  TH2F *   fhECellTotalRatio;                   
  TH2F *   fhECellTotalLogRatio;                
  TH2F **  fhECellTotalRatioMod;                
  TH2F **  fhECellTotalLogRatioMod;             

  // Exotic studies
  
  TH2F *   fhExoNCell  [10][5] ;                
  TH2F *   fhExoL0     [10][5] ;                
  TH2F *   fhExoL1     [10][5] ;                
  TH2F *   fhExoECross [10][5] ;                
  TH2F *   fhExoTime   [10][5] ;                
  TH2F *   fhExoDTime  [10]    ;                
  TH2F *   fhExoL0NCell[10][5] ;                
  TH2F *   fhExoL0ECross      ;                 
  TH2F *   fhExoL1NCell[10][5] ;                
  TH2F *   fhExoL1ECross      ;                 
   
  // Pure MC histograms

  enum mcTypes { kmcPhoton   = 0, kmcPi0        = 1, kmcEta = 2,
                 kmcElectron = 3, kmcPhotonConv = 4,
                 kmcNeHadron = 5, kmcChHadron   = 6             };
  
  TH2F *   fhRecoMCE[7][2]  ;                   
  TH2F *   fhRecoMCPhi[7][2] ;                  
  TH2F *   fhRecoMCEta[7][2] ;                  
  TH2F *   fhRecoMCDeltaE[7][2]  ;              
  TH2F *   fhRecoMCRatioE[7][2]  ;              
  TH2F *   fhRecoMCDeltaPhi[7][2];              
  TH2F *   fhRecoMCDeltaEta[7][2];              
  
  TH1F *   fhGenMCE [4]     ;                   
  TH1F *   fhGenMCPt[4]     ;                   
  TH2F *   fhGenMCEtaPhi[4] ;                   
  TH1F *   fhGenMCAccE [4]     ;                
  TH1F *   fhGenMCAccPt[4]     ;                
  TH2F *   fhGenMCAccEtaPhi[4] ;                
  
  TH2F *   fhEMVxyz    ;                        
  TH2F *   fhEMR       ;                        
  TH2F *   fhHaVxyz    ;                        
  TH2F *   fhHaR       ;                        
  
  // Histograms for MC track-matching
  
  TH2F *   fh1EOverP;                           
  TH2F *   fh2dR;                               
  TH2F *   fh2EledEdx;                          
  TH2F *   fh2MatchdEdx;                        
  TH2F *   fh1EOverPR02;                        
  TH2F *   fh1EleEOverP;                        

  TH2F *   fh1EOverPMod;                        
  TH2F *   fh2dRMod;                            
  TH2F *   fh2EledEdxMod;                       
  TH2F *   fh2MatchdEdxMod;                     
  TH2F *   fh1EOverPR02Mod;                     
  TH2F *   fh1EleEOverPMod;                     

  TH2F *   fhMCEle1EOverP;                      
  TH1F *   fhMCEle1dR;                          
  TH2F *   fhMCEle2MatchdEdx;                   
  
  TH2F *   fhMCChHad1EOverP;                    
  TH1F *   fhMCChHad1dR;                        
  TH2F *   fhMCChHad2MatchdEdx;                 
  
  TH2F *   fhMCNeutral1EOverP;                  
  TH1F *   fhMCNeutral1dR;                      
  TH2F *   fhMCNeutral2MatchdEdx;               
  
  TH2F *   fhMCEle1EOverPR02;                   
  TH2F *   fhMCChHad1EOverPR02;                 
  TH2F *   fhMCNeutral1EOverPR02;               

  TH2F *   fhMCEle1EleEOverP;                   
  TH2F *   fhMCChHad1EleEOverP;                 
  TH2F *   fhMCNeutral1EleEOverP;               

  TH2F *   fhTrackMatchedDEtaNeg;               
  TH2F *   fhTrackMatchedDPhiNeg;               
  TH2F *   fhTrackMatchedDEtaDPhiNeg;           
  
  TH2F *   fhTrackMatchedDEtaPos;               
  TH2F *   fhTrackMatchedDPhiPos;               
  TH2F *   fhTrackMatchedDEtaDPhiPos;           

  TH2F *   fhTrackMatchedDEtaNegMod;            
  TH2F *   fhTrackMatchedDPhiNegMod;            
  TH2F *   fhTrackMatchedDEtaPosMod;            
  TH2F *   fhTrackMatchedDPhiPosMod;            
  
  AliAnaCalorimeterQA & operator = (const AliAnaCalorimeterQA & qa) ;
    
  AliAnaCalorimeterQA(              const AliAnaCalorimeterQA & qa) ;
  
  ClassDef(AliAnaCalorimeterQA,33) ;

} ;

#endif //ALIANACALORIMETERQA_H