PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.h

   1 #ifndef ALIANAPARTICLEISOLATION_H
   2 #define ALIANAPARTICLEISOLATION_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice     */
   5
   6 //_________________________________________________________________________
   7
   8 // Class for the analysis of particle isolation
   9 // Input is selected particles put in AOD branch (AliAODPWG4ParticleCorrelation)
  10 //
  11 //  Class created from old AliPHOSGammaJet
  12 //  (see AliRoot versions previous Release 4-09)
  13
  14 //-- Author: Gustavo Conesa (INFN-LNF)
  15
  16 // --- ROOT system ---
  17 class TH2F;
  18 class TList ;
  19 class TObjString;
  20
  21 // --- ANALYSIS system ---
  22 #include "AliAnaCaloTrackCorrBaseClass.h"
  23 class AliAODPWG4Particle;
  24 class AliAODPWG4ParticleCorrelation ;
  25
  26
  27 class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
  28
  29  public:
  30   AliAnaParticleIsolation() ; // default ctor
  31   virtual ~AliAnaParticleIsolation() { ; } //virtual dtor
  32
  33   // Main general methods
  34
  35   void         CalculateCaloUEBand    (AliAODPWG4ParticleCorrelation * pCandidate,
  36                                        Float_t & etaBand, Float_t & phiBand) ;
  37   void         CalculateCaloCellUEBand(AliAODPWG4ParticleCorrelation * pCandidate,
  38                                        Float_t & etaBand, Float_t & phiBand) ;
  39   void         CalculateTrackUEBand   (AliAODPWG4ParticleCorrelation * pCandidate,
  40                                        Float_t & etaBand, Float_t & phiBand) ;
  41
  42   void         CalculateCaloSignalInCone    (AliAODPWG4ParticleCorrelation * aodParticle, Float_t & coneptsumCluster) ;
  43   void         CalculateCaloCellSignalInCone(AliAODPWG4ParticleCorrelation * aodParticle, Float_t & coneptsumCell) ;
  44   void         CalculateTrackSignalInCone   (AliAODPWG4ParticleCorrelation * aodParticle, Float_t & coneptsumTrack  ) ;
  45
  46
  47   void         CalculateNormalizeUEBandPerUnitArea(AliAODPWG4ParticleCorrelation * pCandidate,
  48                                                    Float_t coneptsumCluster,       Float_t coneptsumCell,     Float_t coneptsumTrack,
  49                                                    Float_t &etaBandptsumTrackNorm, Float_t &etaBandptsumClusterNorm ) ;
  50
  51   TObjString * GetAnalysisCuts() ;
  52
  53   TList      * GetCreateOutputObjects() ;
  54
  55   void         Init() ;
  56
  57   void         InitParameters() ;
  58
  59   void         MakeAnalysisFillAOD()  ;
  60
  61   void         MakeAnalysisFillHistograms() ;
  62
  63
  64   void         Print( const Option_t * opt ) const ;
  65
  66   //Analysis specific methods
  67
  68   void         FillPileUpHistograms(Int_t clusterID) ;
  69
  70   void         FillAcceptanceHistograms();
  71
  72   void         FillTrackMatchingShowerShapeControlHistograms(AliAODPWG4ParticleCorrelation  * pCandidate,
  73                                                              AliCaloPID * pid) ;
  74
  75   void         MakeSeveralICAnalysis( AliAODPWG4ParticleCorrelation * ph ) ;
  76
  77   // Analysis Setters and Getters
  78
  79   TString      GetCalorimeter()                const { return fCalorimeter       ; }
  80   Int_t        GetNCones()                     const { return fNCones            ; }
  81   Int_t        GetNPtThresFrac()               const { return fNPtThresFrac      ; }
  82   Float_t      GetConeSizes(Int_t i)           const { return fConeSizes[i]      ; }
  83   Float_t      GetPtThresholds(Int_t i)        const { return fPtThresholds[i]   ; }
  84   Float_t      GetSumPtThresholds(Int_t i)     const { return fSumPtThresholds[i]; }
  85   Float_t      GetPtFractions(Int_t i)         const { return fPtFractions[i]    ; }
  86
  87   void         SetCalorimeter(TString & det)         { fCalorimeter     = det    ; }
  88   void         SetNCones(Int_t ncs)                  { fNCones          = ncs    ; }
  89   void         SetNPtThresFrac(Int_t npt)            { fNPtThresFrac    = npt    ; }
  90   void         SetConeSizes(Int_t i, Float_t r)      { fConeSizes[i]    = r      ; }
  91   void         SetPtThresholds(Int_t i, Float_t pt)  { fPtThresholds[i] = pt     ; }
  92   void         SetPtFractions(Int_t i, Float_t pt)   { fPtFractions[i]  = pt     ; }
  93   void         SetSumPtThresholds(Int_t i, Float_t pt){ fSumPtThresholds[i] = pt ; }
  94
  95
  96   Bool_t       IsReIsolationOn()               const { return fReMakeIC          ; }
  97   void         SwitchOnReIsolation()                 { fReMakeIC      = kTRUE    ; }
  98   void         SwitchOffReIsolation()                { fReMakeIC      = kFALSE   ; }
  99
 100   Bool_t       IsSeveralIsolationOn()          const { return fMakeSeveralIC     ; }
 101   void         SwitchOnSeveralIsolation()            { fMakeSeveralIC = kTRUE    ; }
 102   void         SwitchOffSeveralIsolation()           { fMakeSeveralIC = kFALSE   ; }
 103
 104   void         SwitchOnFillPileUpHistograms()        { fFillPileUpHistograms = kTRUE  ; }
 105   void         SwitchOffFillPileUpHistograms()       { fFillPileUpHistograms = kFALSE ; }
 106
 107   void         SwitchOnTMHistoFill()                 { fFillTMHisto   = kTRUE    ; }
 108   void         SwitchOffTMHistoFill()                { fFillTMHisto   = kFALSE   ; }
 109
 110   void         SwitchOnSSHistoFill()                 { fFillSSHisto   = kTRUE    ; }
 111   void         SwitchOffSSHistoFill()                { fFillSSHisto   = kFALSE   ; }
 112
 113   //Histogrammes setters and getters
 114
 115   void         SetHistoPtSumRangeAndNBins(Float_t min, Float_t max, Int_t n)       {
 116     fHistoNPtSumBins = n ;    fHistoPtSumMax = max ;    fHistoPtSumMin = min ;     }
 117
 118   Int_t        GetHistoNPtSumBins()            const { return fHistoNPtSumBins   ; }
 119   Float_t      GetHistoPtSumMin()              const { return fHistoPtSumMin     ; }
 120   Float_t      GetHistoPtSumMax()              const { return fHistoPtSumMax     ; }
 121
 122   void         SetHistoPtInConeRangeAndNBins(Float_t min, Float_t max, Int_t n)    {
 123     fHistoNPtInConeBins = n ; fHistoPtInConeMax = max ; fHistoPtInConeMin = min  ; }
 124
 125   Int_t        GetHistoNPtInConeBins()         const { return fHistoNPtInConeBins; }
 126   Float_t      GetHistoPtInConeMin()           const { return fHistoPtInConeMin  ; }
 127   Float_t      GetHistoPtInConeMax()           const { return fHistoPtInConeMax  ; }
 128
 129  // For histograms in arrays, index in the array, corresponding to a particle
 130   enum mcPTypes   { kmcPPhoton = 0,       kmcPPi0Decay = 1,       kmcPOtherDecay = 2,  kmcPOther = 3,
 131                     kmcPPrompt = 4,       kmcPFragmentation = 5,  kmcPISR = 6           };
 132
 133  private:
 134
 135   TString  fCalorimeter ;                         // Calorimeter where neutral particles in cone for isolation are;
 136   Bool_t   fReMakeIC ;                            // Do isolation analysis
 137   Bool_t   fMakeSeveralIC ;                       // Do analysis for different IC
 138   Bool_t   fFillPileUpHistograms;                 // Fill pile-up related histograms
 139   Bool_t   fFillTMHisto;                          // Fill track matching plots
 140   Bool_t   fFillSSHisto;                          // Fill Shower shape plots
 141
 142   // Analysis data members for multiple cones and pt thresholds
 143   Int_t    fNCones ;                              //! Number of cone sizes to test
 144   Int_t    fNPtThresFrac ;                        //! Number of ptThres and ptFrac to test
 145
 146   Float_t  fConeSizes[5] ;                        //! Array with cones to test
 147   Float_t  fPtThresholds[5] ;                     //! Array with pt thresholds to test
 148   Float_t  fPtFractions[5] ;                      //! Array with pt thresholds to test frac
 149   Float_t  fSumPtThresholds[5] ;                  //! Array with pt thresholds to test frac
 150
 151   //Histograms
 152
 153   TH1F *   fhEIso ;                               //! Number of isolated particles vs energy
 154   TH1F *   fhPtIso ;                              //! Number of isolated particles vs pT
 155   TH2F *   fhPtCentralityIso ;                    //! centrality vs pT
 156   TH2F *   fhPtEventPlaneIso ;                    //! event plane angle vs pT
 157   TH2F *   fhPtNLocMaxIso ;                       //! Number of isolated particles vs NLM in cluster
 158   TH2F *   fhPhiIso ;                             //! Phi of isolated particles
 159   TH2F *   fhEtaIso ;                             //! eta of isolated particles
 160   TH2F *   fhEtaPhiIso ;                          //! eta vs phi of isolated particles
 161   TH2F *   fhEtaPhiNoIso ;                        //! eta vs phi of not isolated leading particles
 162   TH1F *   fhENoIso ;                             //! Number of not isolated leading particles vs Energy
 163   TH1F *   fhPtNoIso ;                            //! Number of not isolated leading particles vs pT
 164   TH2F *   fhPtNLocMaxNoIso ;                     //! Number of not isolated particles vs NLM in cluster
 165   TH1F *   fhPtDecayIso ;                         //! Number of isolated Pi0 decay particles (invariant mass tag)
 166   TH1F *   fhPtDecayNoIso ;                       //! Number of not isolated Pi0 decay leading particles (invariant mass tag)
 167   TH2F *   fhEtaPhiDecayIso ;                     //! eta vs phi of isolated Pi0 decay particles
 168   TH2F *   fhEtaPhiDecayNoIso ;                   //! eta vs phi of not isolated leading Pi0 decay particles
 169
 170   TH2F *   fhPtInCone ;                           //! Cluster/track Pt in the cone
 171   TH2F *   fhPtClusterInCone ;                    //! Cluster Pt in the cone
 172   TH2F *   fhPtCellInCone ;                       //! Cell amplitude in the cone
 173   TH2F *   fhPtTrackInCone ;                      //! Track Pt in the cone
 174   TH2F *   fhPtTrackInConeOtherBC ;               //! Track Pt in the cone, tracks out of main BC Time window
 175   TH2F *   fhPtTrackInConeOtherBCPileUpSPD ;      //! Track Pt in the cone, tracks out of main BC Time window
 176   TH2F *   fhPtTrackInConeBC0 ;                   //! Track Pt in the cone, tracks in BC=0
 177   TH2F *   fhPtTrackInConeVtxBC0 ;                //! Track Pt in the cone, tracks in BC=0
 178   TH2F *   fhPtTrackInConeBC0PileUpSPD ;          //! Track Pt in the cone, tracks in BC=0
 179   TH2F *   fhPtInConePileUp[7] ;                  //! Particle Pt in the cone, if event is from pile-up (SPD method)
 180   TH2F *   fhPtInConeCent ;                       //! Particle Pt in the cone versus centrality
 181   TH2F *   fhPerpConeSumPt ;                      //! Sum Pt in cone at the perpendicular phi region to trigger axis  (phi +90)
 182   TH2F *   fhPtInPerpCone ;                       //! Particle Pt  in cone at the perpendicular phi region to trigger axis  (phi +90)
 183
 184   TH2F *   fhEtaPhiInConeCluster ;                //! Eta vs. phi of clusters in cone
 185   TH2F *   fhEtaPhiCluster ;                      //! Eta vs. phi of all clusters
 186   TH2F *   fhEtaPhiInConeTrack ;                  //! Eta vs. phi of tracks in cone
 187   TH2F *   fhEtaPhiTrack ;                        //! Eta vs. phi of all tracks
 188
 189   TH2F *   fhEtaBandCluster ;                     //! Accumulated pT in Eta band to estimate UE in cone, only clusters
 190   TH2F *   fhPhiBandCluster ;                     //! Accumulated pT in Phi band to estimate UE in cone, only clusters
 191   TH2F *   fhEtaBandTrack   ;                     //! Accumulated pT in Eta band to estimate UE in cone, only tracks
 192   TH2F *   fhPhiBandTrack   ;                     //! Accumulated pT in Phi band to estimate UE in cone, only tracks
 193   TH2F *   fhEtaBandCell ;                        //! Accumulated pT in Eta band to estimate UE in cone, only cells
 194   TH2F *   fhPhiBandCell ;                        //! Accumulated pT in Phi band to estimate UE in cone, only cells
 195
 196   TH2F *   fhConeSumPt ;                          //! Cluster and tracks Sum Pt Sum Pt in the cone
 197   TH2F *   fhConeSumPtCellTrack ;                 //! Cells and tracks Sum Pt Sum Pt in the cone
 198   TH2F *   fhConeSumPtCell ;                      //! Cells Sum Pt Sum Pt in the cone
 199   TH2F *   fhConeSumPtCluster ;                   //! Clusters Sum Pt Sum Pt in the cone
 200   TH2F *   fhConeSumPtTrack ;                     //! Tracks Sum Pt Sum Pt in the cone
 201   TH2F *   fhConeSumPtEtaBandUECluster;           //! Cluster Sum Pt in the eta band for clusters, before normalization
 202   TH2F *   fhConeSumPtPhiBandUECluster;           //! Cluster Sum Pt in the phi band for clusters, before normalization
 203   TH2F *   fhConeSumPtEtaBandUETrack;             //! Track Sum Pt in the eta band for tracks  , before normalization
 204   TH2F *   fhConeSumPtPhiBandUETrack;             //! Track Sum Pt in the phi badn for tracks  , before normalization
 205   TH2F *   fhConeSumPtEtaBandUECell;              //! Cell Sum amplitude in the eta band for cells, before normalization
 206   TH2F *   fhConeSumPtPhiBandUECell;              //! Cell Sum amplitude in the phi band for cells, before normalization
 207
 208   TH2F *   fhConeSumPtTrigEtaPhi ;                //! Cluster and tracks Sum Pt Sum Pt in the cone, per eta-phi bin of trigger,
 209   TH2F *   fhConeSumPtCellTrackTrigEtaPhi ;       //! Cell and tracks Sum Pt Sum Pt in the cone, per eta-phi bin of trigger,
 210   TH2F *   fhConeSumPtEtaBandUEClusterTrigEtaPhi; //! Cluster Sum Pt in the eta band for clusters, per eta-phi bin of trigger,before normalization
 211   TH2F *   fhConeSumPtPhiBandUEClusterTrigEtaPhi; //! Cluster Sum Pt in the phi band for clusters, per eta-phi bin of trigger, before normalization
 212   TH2F *   fhConeSumPtEtaBandUETrackTrigEtaPhi;   //! Track Sum Pt in the eta band for tracks  , per eta-phi bin of trigger, before normalization
 213   TH2F *   fhConeSumPtPhiBandUETrackTrigEtaPhi;   //! Track Sum Pt in the phi badn for tracks  , per eta-phi bin of trigger, before normalization
 214   TH2F *   fhConeSumPtEtaBandUECellTrigEtaPhi;    //! Cluster Sum amplitude in the eta band for cells, per eta-phi bin of trigger, before normalization
 215   TH2F *   fhConeSumPtPhiBandUECellTrigEtaPhi;    //! Cluster Sum amplitude in the phi band for cells, per eta-phi bin of trigger, before normalization
 216
 217   TH2F *   fhConeSumPtEtaUESub;                   //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs pT trigger
 218   TH2F *   fhConeSumPtPhiUESub;                   //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs pT trigger
 219   TH2F *   fhConeSumPtEtaUESubTrigEtaPhi;         //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs eta-phi trigger
 220   TH2F *   fhConeSumPtPhiUESubTrigEtaPhi;         //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs eta-phi trigger
 221
 222   TH2F *   fhConeSumPtEtaUESubTrackCell;          //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs pT trigger
 223   TH2F *   fhConeSumPtPhiUESubTrackCell;          //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs pT trigger
 224   TH2F *   fhConeSumPtEtaUESubTrackCellTrigEtaPhi;//! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs eta-phi trigger
 225   TH2F *   fhConeSumPtPhiUESubTrackCellTrigEtaPhi;//! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs eta-phi trigger
 226
 227   TH2F *   fhConeSumPtEtaUESubCluster;            //! Cluster Sum Pt in the cone after bkg subtraction, vs pT trigger
 228   TH2F *   fhConeSumPtPhiUESubCluster;            //! Cluster Sum Pt in the cone after bkg subtraction, vs pT trigger
 229   TH2F *   fhConeSumPtEtaUESubClusterTrigEtaPhi;  //! Cluster Sum Pt in the cone after bkg subtraction, vs eta-phi trigger
 230   TH2F *   fhConeSumPtPhiUESubClusterTrigEtaPhi;  //! Cluster Sum Pt in the cone after bkg subtraction, vs eta-phi trigger
 231
 232   TH2F *   fhConeSumPtEtaUESubCell;               //! Cell Sum amplitude in the cone after bkg subtraction, vs pT trigger
 233   TH2F *   fhConeSumPtPhiUESubCell;               //! Cell Sum amplitude in the cone after bkg subtraction, vs pT trigger
 234   TH2F *   fhConeSumPtEtaUESubCellTrigEtaPhi;     //! Cell Sum amplitude in the cone after bkg subtraction, vs eta-phi trigger
 235   TH2F *   fhConeSumPtPhiUESubCellTrigEtaPhi;     //! Cell Sum amplitude in the cone after bkg subtraction, vs eta-phi trigger
 236
 237   TH2F *   fhConeSumPtEtaUESubTrack;              //! Track Sum Pt in the cone after bkg subtraction, vs pT trigger
 238   TH2F *   fhConeSumPtPhiUESubTrack;              //! Track Sum Pt in the cone after bkg subtraction, vs pT trigger
 239   TH2F *   fhConeSumPtEtaUESubTrackTrigEtaPhi;    //! Track Sum Pt in the cone after bkg subtraction, vs eta-phi trigger
 240   TH2F *   fhConeSumPtPhiUESubTrackTrigEtaPhi;    //! Track Sum Pt in the cone after bkg subtraction, vs eta-phi trigger
 241
 242   TH2F *   fhFractionTrackOutConeEta;             //! Fraction of cone out of tracks acceptance in eta
 243   TH2F *   fhFractionTrackOutConeEtaTrigEtaPhi;   //! Fraction of cone out of tracks acceptance in eta, vs trigger eta-phi
 244   TH2F *   fhFractionClusterOutConeEta;           //! Fraction of cone out of clusters acceptance in eta
 245   TH2F *   fhFractionClusterOutConeEtaTrigEtaPhi; //! Fraction of cone out of clusters acceptance in eta, vs trigger eta-phi
 246   TH2F *   fhFractionClusterOutConePhi;           //! Fraction of cone out of clusters acceptance in phi
 247   TH2F *   fhFractionClusterOutConePhiTrigEtaPhi; //! Fraction of cone out of clusters acceptance in phi, vs trigger eta-phi
 248
 249   TH2F *   fhFractionCellOutConeEta;              //! Fraction of cone out of cells acceptance in eta
 250   TH2F *   fhFractionCellOutConeEtaTrigEtaPhi;    //! Fraction of cone out of cells acceptance in eta, vs trigger eta-phi
 251   TH2F *   fhFractionCellOutConePhi;              //! Fraction of cone out of cells acceptance in phi
 252   TH2F *   fhFractionCellOutConePhiTrigEtaPhi;    //! Fraction of cone out of cells acceptance in phi, vs trigger eta-phi
 253
 254   TH2F *   fhConeSumPtClustervsTrack ;            //! Cluster vs tracks Sum Pt Sum Pt in the cone
 255   TH2F *   fhConeSumPtEtaUESubClustervsTrack ;    //! Cluster vs tracks Sum Pt Sum Pt in the cone, after subtraction in eta band
 256   TH2F *   fhConeSumPtPhiUESubClustervsTrack ;    //! Cluster vs tracks Sum Pt Sum Pt in the cone, after subtraction in phi band
 257   TH2F *   fhConeSumPtCellvsTrack;                //! Cell vs tracks Sum Pt Sum Pt in the cone
 258   TH2F *   fhConeSumPtEtaUESubCellvsTrack ;       //! Cell vs tracks Sum Pt Sum Pt in the cone, after subtraction in eta band
 259   TH2F *   fhConeSumPtPhiUESubCellvsTrack ;       //! Cell vs tracks Sum Pt Sum Pt in the cone, after subtraction in phi band
 260
 261   TH2F *   fhEtaBandClustervsTrack ;              //! Accumulated pT in Eta band to estimate UE in cone, clusters vs tracks
 262   TH2F *   fhPhiBandClustervsTrack ;              //! Accumulated pT in Phi band to estimate UE in cone, clusters vs tracks
 263   TH2F *   fhEtaBandNormClustervsTrack ;          //! Accumulated pT in Eta band to estimate UE in cone, normalized to cone size, clusters vs tracks
 264   TH2F *   fhPhiBandNormClustervsTrack ;          //! Accumulated pT in Phi band to estimate UE in cone, normalized to cone size, clusters vs tracks
 265   TH2F *   fhEtaBandCellvsTrack ;                 //! Accumulated pT in Eta band to estimate UE in cone, cells vs tracks
 266   TH2F *   fhPhiBandCellvsTrack ;                 //! Accumulated pT in Phi band to estimate UE in cone, cells vs tracks
 267   TH2F *   fhEtaBandNormCellvsTrack ;             //! Accumulated pT cell in Eta band to estimate UE in cone, normalized to cone size, clusters vs tracks
 268   TH2F *   fhPhiBandNormCellvsTrack ;             //! Accumulated pT cell in Phi band to estimate UE in cone, normalized to cone
 269
 270   TH2F *   fhConeSumPtSubvsConeSumPtTotPhiTrack;       //! Tracks, phi band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub
 271   TH2F *   fhConeSumPtSubNormvsConeSumPtTotPhiTrack;   //! Tracks, phi band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub
 272   TH2F *   fhConeSumPtSubvsConeSumPtTotEtaTrack;       //! Tracks, eta band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub
 273   TH2F *   fhConeSumPtSubNormvsConeSumPtTotEtaTrack;   //! Tracks, eta band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub
 274   TH2F *   fhConeSumPtSubvsConeSumPtTotPhiCluster;     //! Clusters, phi band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub
 275   TH2F *   fhConeSumPtSubNormvsConeSumPtTotPhiCluster; //! Clusters, phi band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub
 276   TH2F *   fhConeSumPtSubvsConeSumPtTotEtaCluster;     //! Clusters, eta band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub
 277   TH2F *   fhConeSumPtSubNormvsConeSumPtTotEtaCluster; //! Clusters, eta band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub
 278   TH2F *   fhConeSumPtSubvsConeSumPtTotPhiCell;        //! Cells, phi band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub
 279   TH2F *   fhConeSumPtSubNormvsConeSumPtTotPhiCell;    //! Cells, phi band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub
 280   TH2F *   fhConeSumPtSubvsConeSumPtTotEtaCell;        //! Cells, eta band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub
 281   TH2F *   fhConeSumPtSubNormvsConeSumPtTotEtaCell;    //! Cells, eta band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub
 282   TH2F *   fhConeSumPtVSUETracksEtaBand;               //! fhConeSumPtVSUETracksEtaBand
 283   TH2F *   fhConeSumPtVSUETracksPhiBand;               //! fhConeSumPtVSUETracksPhiBand
 284   TH2F *   fhConeSumPtVSUEClusterEtaBand;              //! fhConeSumPtVSUEClusterEtaBand
 285   TH2F *   fhConeSumPtVSUEClusterPhiBand;              //! fhConeSumPtVSUEClusterPhiBand
 286
 287   //MC
 288
 289   TH2F *   fhEtaPrimMC  [7];                      //! Pt vs Eta of generated photon
 290   TH1F *   fhEPrimMC    [7];                      //! Number of generated photon vs E
 291   TH1F *   fhPtPrimMCiso[7];                      //! Number of generated isolated photon vs pT
 292   TH1F *   fhPtIsoPrompt;                         //! Number of isolated prompt gamma
 293   TH2F *   fhPhiIsoPrompt;                        //! Phi of isolated prompt gamma
 294   TH2F *   fhEtaIsoPrompt;                        //! eta of isolated prompt gamma
 295   TH1F *   fhPtThresIsolatedPrompt[5][5];         //! Isolated prompt gamma with pt threshold
 296   TH1F *   fhPtFracIsolatedPrompt[5][5];          //! Isolated prompt gamma with pt frac
 297   TH2F *   fhPtSumIsolatedPrompt[5];              //! Isolated prompt gamma with threshold on cone pt sume
 298   TH1F *   fhPtIsoFragmentation;                  //! Number of isolated fragmentation gamma
 299   TH2F *   fhPhiIsoFragmentation;                 //! Phi of isolated fragmentation gamma
 300   TH2F *   fhEtaIsoFragmentation;                 //! eta of isolated fragmentation gamma
 301   TH1F *   fhPtThresIsolatedFragmentation[5][5];  //! Isolated fragmentation gamma with pt threshold
 302   TH1F *   fhPtFracIsolatedFragmentation[5][5];   //! Isolated fragmentation gamma with pt frac
 303   TH2F *   fhPtSumIsolatedFragmentation[5];       //! Isolated fragmentation gamma with threshold on cone pt sume
 304   TH1F *   fhPtIsoPi0;                            //! Number of isolated pi0 (2 gamma)
 305   TH2F *   fhPhiIsoPi0;                           //! Phi of isolated pi0 (2 gamma)
 306   TH2F *   fhEtaIsoPi0;                           //! eta of isolated pi0 (2 gamma)
 307   TH1F *   fhPtThresIsolatedPi0[5][5];            //! Isolated pi0 (2 gamma) with pt threshold
 308   TH1F *   fhPtFracIsolatedPi0[5][5];             //! Isolated pi0 (2 gamma) with pt frac
 309   TH2F *   fhPtSumIsolatedPi0[5];                 //! Isolated pi0 (2 gamma) with threshold on cone pt sum
 310   TH1F *   fhPtIsoPi0Decay;                       //! Number of isolated pi0 decay gamma
 311   TH2F *   fhPhiIsoPi0Decay;                      //! Phi of isolated pi0 decay gamma
 312   TH2F *   fhEtaIsoPi0Decay;                      //! eta of isolated pi0 decay gamma
 313   TH1F *   fhPtThresIsolatedPi0Decay[5][5];       //! Isolated pi0 decay gamma with pt threshold
 314   TH1F *   fhPtFracIsolatedPi0Decay[5][5];        //! Isolated pi0 decay gamma with pt frac
 315   TH2F *   fhPtSumIsolatedPi0Decay[5];            //! Isolated pi0 decay gamma with threshold on cone pt sume
 316   TH1F *   fhPtIsoEtaDecay;                       //! Number of isolated eta decay gamma
 317   TH2F *   fhPhiIsoEtaDecay;                      //! Phi of isolated eta decay gamma
 318   TH2F *   fhEtaIsoEtaDecay;                      //! eta of isolated eta decay gamma
 319   TH1F *   fhPtThresIsolatedEtaDecay[5][5];       //! Isolated eta decay gamma with pt threshold
 320   TH1F *   fhPtFracIsolatedEtaDecay[5][5];        //! Isolated eta decay gamma with pt frac
 321   TH2F *   fhPtSumIsolatedEtaDecay[5];            //! Isolated eta fecay gamma with threshold on cone pt sume
 322   TH1F *   fhPtIsoOtherDecay;                     //! Number of isolated other decay gamma
 323   TH2F *   fhPhiIsoOtherDecay;                    //! Phi of isolated other decay gamma
 324   TH2F *   fhEtaIsoOtherDecay;                    //! eta of isolated other decay gamma
 325   TH1F *   fhPtThresIsolatedOtherDecay[5][5];     //! Isolated OtherDecay gamma with pt threshold
 326   TH1F *   fhPtFracIsolatedOtherDecay[5][5];      //! Isolated OtherDecay gamma with pt frac
 327   TH2F *   fhPtSumIsolatedOtherDecay[5];          //! Isolated OtherDecay gamma with threshold on cone pt sume
 328 //  TH1F *   fhPtIsoConversion;                     //! Number of isolated Conversion gamma
 329 //  TH2F *   fhPhiIsoConversion;                    //! Phi of isolated Conversion gamma
 330 //  TH2F *   fhEtaIsoConversion;                    //! eta of isolated Conversion gamma
 331 //  TH1F *   fhPtThresIsolatedConversion[5][5];     //! Isolated Conversion gamma with pt threshold
 332 //  TH1F *   fhPtFracIsolatedConversion[5][5];      //! Isolated Conversion gamma with pt frac
 333 //  TH2F *   fhPtSumIsolatedConversion[5];          //! Isolated Conversion gamma with threshold on cone pt sume
 334   TH1F *   fhPtIsoHadron;                         //! Number of isolated Hadron
 335   TH2F *   fhPhiIsoHadron;                        //! Phi of isolated Hadron
 336   TH2F *   fhEtaIsoHadron;                        //! eta of isolated Hadron
 337   TH1F *   fhPtThresIsolatedHadron[5][5];         //! Isolated Hadron gamma with pt threshold
 338   TH1F *   fhPtFracIsolatedHadron[5][5];          //! Isolated Hadron gamma with pt frac
 339   TH2F *   fhPtSumIsolatedHadron[5];              //! Isolated Hadron gamma with threshold on cone pt sume
 340
 341   // Multi Cut analysis Several IC
 342   TH1F *   fhPtNoIsoPi0;                          //! Number of not isolated leading pi0 (2 gamma)
 343   TH1F *   fhPtNoIsoPi0Decay;                     //! Number of not isolated leading pi0 decay gamma
 344   TH1F *   fhPtNoIsoEtaDecay;                     //! Number of not isolated leading eta decay gamma
 345   TH1F *   fhPtNoIsoOtherDecay;                   //! Number of not isolated leading other decay gamma
 346   TH1F *   fhPtNoIsoPrompt;                       //! Number of not isolated leading prompt gamma
 347   TH1F *   fhPtIsoMCPhoton;                       //! Number of isolated leading gamma
 348   TH1F *   fhPtNoIsoMCPhoton;                     //! Number of not isolated leading gamma
 349 //  TH1F *   fhPtNoIsoConversion;                   //! Number of not isolated leading conversion gamma
 350   TH1F *   fhPtNoIsoFragmentation;                //! Number of not isolated leading fragmentation gamma
 351   TH1F *   fhPtNoIsoHadron;                       //! Number of not isolated leading hadrons
 352
 353   TH2F *   fhSumPtLeadingPt[5] ;                  //! Sum Pt in the cone
 354   TH2F *   fhPtLeadingPt[5] ;                     //! Particle Pt in the cone
 355   TH2F *   fhPerpSumPtLeadingPt[5] ;              //! Sum Pt in the cone at the perpendicular phi region to trigger axis  (phi +90)
 356   TH2F *   fhPerpPtLeadingPt[5];                  //! Sum Pt in the cone at the perpendicular phi region to trigger axis  (phi +90)
 357
 358   TH1F *   fhPtThresIsolated[5][5] ;              //! Isolated particle with pt threshold
 359   TH1F *   fhPtFracIsolated[5][5] ;               //! Isolated particle with pt threshold frac
 360   TH1F *   fhPtSumIsolated[5][5] ;                //! Isolated particle with threshold on cone pt sum
 361
 362   TH2F *   fhEtaPhiPtThresIso[5][5] ;             //! eta vs phi of isolated particles with pt threshold
 363   TH2F *   fhEtaPhiPtThresDecayIso[5][5] ;        //! eta vs phi of isolated particles with pt threshold
 364   TH1F *   fhPtPtThresDecayIso[5][5] ;            //! Number of isolated Pi0 decay particles (invariant mass tag) with pt threshold
 365
 366   TH2F *   fhEtaPhiPtFracIso[5][5] ;              //! eta vs phi of isolated particles with pt frac
 367   TH2F *   fhEtaPhiPtFracDecayIso[5][5] ;         //! eta vs phi of isolated particles with pt frac
 368   TH1F *   fhPtPtFracDecayIso[5][5] ;             //! Number of isolated Pi0 decay particles (invariant mass tag) with pt fra
 369
 370   TH2F *   fhEtaPhiPtSumIso[5][5] ;               //! eta vs phi of isolated particles with pt sum
 371   TH2F *   fhEtaPhiPtSumDecayIso[5][5] ;          //! eta vs phi of isolated particles with pt sum
 372   TH1F *   fhPtPtSumDecayIso[5][5] ;              //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum
 373
 374   TH2F *   fhEtaPhiSumDensityIso[5][5];           //! Isolated particle with threshold on cone sum density
 375   TH2F *   fhEtaPhiSumDensityDecayIso[5][5];      //! Isolated particle with threshold on cone sum density
 376   TH1F *   fhPtSumDensityIso[5][5];               //! Isolated particle with threshold on cone sum density
 377   TH1F *   fhPtSumDensityDecayIso[5][5];          //! Isolated decay particle with threshold on cone sum density
 378
 379   TH1F *   fhPtFracPtSumIso[5][5] ;               //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum
 380   TH1F *   fhPtFracPtSumDecayIso[5][5] ;          //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum
 381   TH2F *   fhEtaPhiFracPtSumIso[5][5];            //! Isolated particle with threshold on cone sum density
 382   TH2F *   fhEtaPhiFracPtSumDecayIso[5][5];       //! Isolated particle with threshold on cone sum density
 383
 384   // Track matching studies
 385   TH2F *   fhTrackMatchedDEta[2]     ;            //! Eta distance between track and cluster vs cluster E
 386   TH2F *   fhTrackMatchedDPhi[2]     ;            //! Phi distance between track and cluster vs cluster E
 387   TH2F *   fhTrackMatchedDEtaDPhi[2] ;            //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV
 388   TH2F *   fhdEdx[2]  ;                           //! matched track dEdx vs cluster E
 389   TH2F *   fhEOverP[2];                           //! matched track E cluster over P track vs cluster E, after dEdx cut
 390   TH2F *   fhTrackMatchedMCParticle[2];           //! Trace origin of matched particle
 391
 392   // Shower Shape histograms
 393   TH2F *   fhELambda0[2];                         //! Shower shape of (non) isolated photons (do not apply SS cut previously)
 394   TH2F *   fhPtLambda0[2];                        //! Shower shape of (non) isolated photons (do not apply SS cut previously)
 395   TH2F *   fhELambda1[2];                         //! Shower shape of (non) isolated photons (do not apply SS cut previously)
 396   TH2F *   fhELambda0SSBkg;                       //! Shower shape of non isolated photons close to isolation threshold (do not apply SS cut previously)
 397   TH2F *   fhELambda0TRD[2];                      //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
 398   TH2F *   fhPtLambda0TRD[2];                     //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
 399   TH2F *   fhELambda1TRD[2];                      //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
 400   TH2F *   fhELambda0MCPhoton[2];                 //! Shower shape of (non) isolated photon candidates originated by particle photon not decay (do not apply SS cut previously)
 401   TH2F *   fhPtLambda0MCPhotonPrompt[2];          //! Shower shape of (non) isolated photon candidates originated by particle photon not decay (do not apply SS cut previously)
 402   TH2F *   fhPtLambda0MCPhotonFrag[2];            //! Shower shape of (non) isolated photon candidates originated by particle photon not decay (do not apply SS cut previously)
 403   TH2F *   fhELambda0MCPi0[2];                    //! Shower shape of (non) isolated photon candidates originated by particle 2 merged pi0 photons (do not apply SS cut previously)
 404   TH2F *   fhELambda0MCPi0Decay[2];               //! Shower shape of (non) isolated photon candidates originated by particle pi0 decay photon (do not apply SS cut previously)
 405   TH2F *   fhELambda0MCEtaDecay[2];               //! Shower shape of (non) isolated photon candidates originated by particle eta decay photon (do not apply SS cut previously)
 406   TH2F *   fhELambda0MCOtherDecay[2];             //! Shower shape of (non) isolated photon candidates originated by particle other decay photon (do not apply SS cut previously)
 407   TH2F *   fhELambda0MCHadron[2];                 //! Shower shape of (non) isolated photon candidates originated by particle other hadrons (do not apply SS cut previously)
 408
 409   // Local maxima
 410   TH2F *   fhNLocMax[2];                          //! number of maxima in selected clusters
 411   TH2F *   fhELambda0LocMax1[2] ;                 //! E vs lambda0 of selected cluster, 1 local maxima in cluster
 412   TH2F *   fhELambda1LocMax1[2] ;                 //! E vs lambda1 of selected cluster, 1 local maxima in cluster
 413   TH2F *   fhELambda0LocMax2[2] ;                 //! E vs lambda0 of selected cluster, 2 local maxima in cluster
 414   TH2F *   fhELambda1LocMax2[2] ;                 //! E vs lambda1 of selected cluster, 2 local maxima in cluster
 415   TH2F *   fhELambda0LocMaxN[2] ;                 //! E vs lambda0 of selected cluster, N>2 local maxima in cluster
 416   TH2F *   fhELambda1LocMaxN[2] ;                 //! E vs lambda1 of selected cluster, N>2 local maxima in cluster
 417
 418   // Pile-up
 419   TH1F *   fhEIsoPileUp[7] ;                      //! Number of isolated particles
 420   TH1F *   fhPtIsoPileUp[7] ;                     //! Number of isolated particles
 421   TH1F *   fhENoIsoPileUp[7] ;                    //! Number of not isolated particles
 422   TH1F *   fhPtNoIsoPileUp[7] ;                   //! Number of not isolated particles
 423   TH2F *   fhTimeENoCut;                          //! time of cluster vs E, no cut
 424   TH2F *   fhTimeESPD;                            //! time of cluster vs E, IsSPDPileUp
 425   TH2F *   fhTimeESPDMulti;                       //! time of cluster vs E, IsSPDPileUpMulti
 426   TH2F *   fhTimeNPileUpVertSPD;                  //! time of cluster vs n pile-up vertices from SPD
 427   TH2F *   fhTimeNPileUpVertTrack;                //! time of cluster vs n pile-up vertices from Tracks
 428   TH2F *   fhTimeNPileUpVertContributors;         //! time of cluster vs n pile-up vertex from SPD contributors
 429   TH2F *   fhTimePileUpMainVertexZDistance;       //! time of cluster vs difference of z main vertex and pile-up vertex
 430   TH2F *   fhTimePileUpMainVertexZDiamond;        //! time of cluster vs difference of z diamond and pile-up vertex
 431
 432   //Histograms settings
 433   Int_t    fHistoNPtSumBins;                      // Number of bins in PtSum histograms
 434   Float_t  fHistoPtSumMax;                        // PtSum maximum in histogram
 435   Float_t  fHistoPtSumMin;                            // PtSum minimum in histogram
 436   Int_t    fHistoNPtInConeBins;                   // Number of bins in PtInCone histogram
 437   Float_t  fHistoPtInConeMax;                     // PtInCone maximum in histogram
 438   Float_t  fHistoPtInConeMin;                     // PtInCone maximum in histogram
 439
 440   AliAnaParticleIsolation(              const AliAnaParticleIsolation & iso) ; // cpy ctor
 441   AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & iso) ; // cpy assignment
 442
 443   ClassDef(AliAnaParticleIsolation,25)
 444 } ;
 445
 446
 447 #endif //ALIANAPARTICLEISOLATION_H
 448
 449
 450