]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.h
change location of geometry.root file
[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / AliAnaParticleIsolation.h
index 1a1ed8ff8a54e8eda855083737e98ac7296164ba..913b6b66a8fd08267a49385306f0c7778f01c959 100755 (executable)
@@ -49,9 +49,17 @@ class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
  
   //Analysis specific methods 
   
-  void         FillTrackMatchingShowerShapeControlHistograms(const Int_t clusterID, 
-                                                             const Int_t nLocMax,
-                                                             const Int_t mcTag ) ;
+  void         FillTrackMatchingShowerShapeControlHistograms(const Bool_t isolated,
+                                                             const Int_t  clusterID, 
+                                                             const Int_t  nLocMax,
+                                                             const Int_t  mcTag, 
+                                                            const TObjArray * plCTS,
+                                                            const TObjArray * plNe,
+                                                            AliAODPWG4ParticleCorrelation  * pCandidate,
+                                                            const AliCaloTrackReader * reader, 
+                                                            const AliCaloPID * pid) ;
+  
+
   
   void         MakeSeveralICAnalysis( AliAODPWG4ParticleCorrelation * ph ) ; 
   
@@ -136,8 +144,15 @@ class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
   TH2F *   fhEtaPhiDecayNoIso ;                   //! eta vs phi of not isolated leading Pi0 decay particles
   TH2F *   fhConeSumPt ;                          //! Sum Pt in the cone
   TH2F *   fhPtInCone ;                           //! Particle Pt in the cone
+  TH2F *   fhPtInConeCent ;                       //! Particle Pt in the cone versus centrality
   TH2F *   fhFRConeSumPt ;                        //! Sum Pt in the forward region cone (phi +90)
   TH2F *   fhPtInFRCone ;                         //! Particle Pt in the forward region cone (phi +90 ) 
+  TH2F *   fhPhiUEConeSumPt ;                     //! UE Pt sum in phi band around cone
+  TH2F *   fhEtaUEConeSumPt ;                     //! UE Pt sum in eta band around cone
+  TH2F *   fhEtaBand ;                            //! Eta band to estimate UE in cone
+  TH2F *   fhPhiBand ;                            //! Phi band to estimate UE in cone
+  TH2F *   fhConeSumPtEtaUESub;                   //! Sum Pt in the cone after bkg subtraction
+  TH2F *   fhConeSumPtPhiUESub;                   //! Sum Pt in the cone after bkg subtraction
   
   //MC
   TH1F *   fhPtIsoPrompt;                         //! Number of isolated prompt gamma 
@@ -152,6 +167,12 @@ class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
   TH1F *   fhPtThresIsolatedFragmentation[5][5];  //! Isolated fragmentation gamma with pt threshold 
   TH1F *   fhPtFracIsolatedFragmentation[5][5];   //! Isolated fragmentation gamma with pt frac
   TH2F *   fhPtSumIsolatedFragmentation[5];       //! Isolated fragmentation gamma with threshold on cone pt sume
+  TH1F *   fhPtIsoPi0;                            //! Number of isolated pi0 (2 gamma)
+  TH2F *   fhPhiIsoPi0;                           //! Phi of isolated pi0 (2 gamma)
+  TH2F *   fhEtaIsoPi0;                           //! eta of isolated pi0 (2 gamma)
+  TH1F *   fhPtThresIsolatedPi0[5][5];            //! Isolated pi0 (2 gamma) with pt threshold 
+  TH1F *   fhPtFracIsolatedPi0[5][5];             //! Isolated pi0 (2 gamma) with pt frac
+  TH2F *   fhPtSumIsolatedPi0[5];                 //! Isolated pi0 (2 gamma) with threshold on cone pt sum
   TH1F *   fhPtIsoPi0Decay;                       //! Number of isolated pi0 decay gamma 
   TH2F *   fhPhiIsoPi0Decay;                      //! Phi of isolated pi0 decay gamma
   TH2F *   fhEtaIsoPi0Decay;                      //! eta of isolated pi0 decay gamma
@@ -170,29 +191,36 @@ class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
   TH1F *   fhPtThresIsolatedOtherDecay[5][5];     //! Isolated OtherDecay gamma with pt threshold 
   TH1F *   fhPtFracIsolatedOtherDecay[5][5];      //! Isolated OtherDecay gamma with pt frac
   TH2F *   fhPtSumIsolatedOtherDecay[5];          //! Isolated OtherDecay gamma with threshold on cone pt sume 
-  TH1F *   fhPtIsoConversion;                     //! Number of isolated Conversion gamma 
-  TH2F *   fhPhiIsoConversion;                    //! Phi of isolated Conversion gamma
-  TH2F *   fhEtaIsoConversion;                    //! eta of isolated Conversion gamma
-  TH1F *   fhPtThresIsolatedConversion[5][5];     //! Isolated Conversion gamma with pt threshold 
-  TH1F *   fhPtFracIsolatedConversion[5][5];      //! Isolated Conversion gamma with pt frac
-  TH2F *   fhPtSumIsolatedConversion[5];          //! Isolated Conversion gamma with threshold on cone pt sume
-  TH1F *   fhPtIsoUnknown;                        //! Number of isolated Unknow
-  TH2F *   fhPhiIsoUnknown;                       //! Phi of isolated Unknown
-  TH2F *   fhEtaIsoUnknown;                       //! eta of isolated Unknown
-  TH1F *   fhPtThresIsolatedUnknown[5][5];        //! Isolated Unknown gamma with pt threshold 
-  TH1F *   fhPtFracIsolatedUnknown[5][5];         //! Isolated Unknown gamma with pt frac
-  TH2F *   fhPtSumIsolatedUnknown[5];             //! Isolated Unknown gamma with threshold on cone pt sume
+//  TH1F *   fhPtIsoConversion;                     //! Number of isolated Conversion gamma 
+//  TH2F *   fhPhiIsoConversion;                    //! Phi of isolated Conversion gamma
+//  TH2F *   fhEtaIsoConversion;                    //! eta of isolated Conversion gamma
+//  TH1F *   fhPtThresIsolatedConversion[5][5];     //! Isolated Conversion gamma with pt threshold 
+//  TH1F *   fhPtFracIsolatedConversion[5][5];      //! Isolated Conversion gamma with pt frac
+//  TH2F *   fhPtSumIsolatedConversion[5];          //! Isolated Conversion gamma with threshold on cone pt sume
+  TH1F *   fhPtIsoHadron;                        //! Number of isolated Hadro
+  TH2F *   fhPhiIsoHadron;                       //! Phi of isolated Hadron
+  TH2F *   fhEtaIsoHadron;                       //! eta of isolated Hadron
+  TH1F *   fhPtThresIsolatedHadron[5][5];        //! Isolated Hadron gamma with pt threshold 
+  TH1F *   fhPtFracIsolatedHadron[5][5];         //! Isolated Hadron gamma with pt frac
+  TH2F *   fhPtSumIsolatedHadron[5];             //! Isolated Hadron gamma with threshold on cone pt sume
 
+  // Multi Cut analysis Several IC
+  TH1F *   fhPtNoIsoPi0;                          //! Number of not isolated leading pi0 (2 gamma)
   TH1F *   fhPtNoIsoPi0Decay;                     //! Number of not isolated leading pi0 decay gamma 
   TH1F *   fhPtNoIsoEtaDecay;                     //! Number of not isolated leading eta decay gamma 
   TH1F *   fhPtNoIsoOtherDecay;                   //! Number of not isolated leading other decay gamma 
   TH1F *   fhPtNoIsoPrompt;                       //! Number of not isolated leading prompt gamma 
   TH1F *   fhPtIsoMCPhoton;                       //! Number of isolated leading gamma 
   TH1F *   fhPtNoIsoMCPhoton;                     //! Number of not isolated leading gamma 
-  TH1F *   fhPtNoIsoConversion;                   //! Number of not isolated leading conversion gamma 
+//  TH1F *   fhPtNoIsoConversion;                   //! Number of not isolated leading conversion gamma 
   TH1F *   fhPtNoIsoFragmentation;                //! Number of not isolated leading fragmentation gamma 
-  TH1F *   fhPtNoIsoUnknown;                      //! Number of not isolated leading hadrons 
+  TH1F *   fhPtNoIsoHadron;                      //! Number of not isolated leading hadrons 
   
+  TH2F *   fhSumPtLeadingPt[5] ;                  //! Sum Pt in the cone
+  TH2F *   fhPtLeadingPt[5] ;                     //! Particle Pt in the cone
+  TH2F *   fhFRSumPtLeadingPt[5] ;                //! Sum Pt in the forward region cone (phi +90)
+  TH2F *   fhFRPtLeadingPt[5];                    //! Particle Pt in the forward region cone (phi +90 ) 
+
   TH1F *   fhPtThresIsolated[5][5] ;              //! Isolated particle with pt threshold 
   TH1F *   fhPtFracIsolated[5][5] ;               //! Isolated particle with pt threshold frac
   TH1F *   fhPtSumIsolated[5][5] ;                //! Isolated particle with threshold on cone pt sum
@@ -220,32 +248,39 @@ class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
   TH2F *   fhEtaPhiFracPtSumDecayIso[5][5];       //! Isolated particle with threshold on cone sum density
  
   // Track matching studies
-  TH2F *   fhTrackMatchedDEta     ;               //! Eta distance between track and cluster vs cluster E
-  TH2F *   fhTrackMatchedDPhi     ;               //! Phi distance between track and cluster vs cluster E
-  TH2F *   fhTrackMatchedDEtaDPhi ;               //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV
-  TH2F *   fhdEdx  ;                              //! matched track dEdx vs cluster E 
-  TH2F *   fhEOverP;                              //! matched track E cluster over P track vs cluster E, after dEdx cut 
-  TH2F *   fhTrackMatchedMCParticle;              //! Trace origin of matched particle
+  TH2F *   fhTrackMatchedDEta[2]     ;            //! Eta distance between track and cluster vs cluster E
+  TH2F *   fhTrackMatchedDPhi[2]     ;            //! Phi distance between track and cluster vs cluster E
+  TH2F *   fhTrackMatchedDEtaDPhi[2] ;            //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV
+  TH2F *   fhdEdx[2]  ;                           //! matched track dEdx vs cluster E 
+  TH2F *   fhEOverP[2];                           //! matched track E cluster over P track vs cluster E, after dEdx cut 
+  TH2F *   fhTrackMatchedMCParticle[2];           //! Trace origin of matched particle
 
   // Shower Shape histograms
-  TH2F *   fhELambda0;                            //! Shower shape of isolated photons
-  TH2F *   fhELambda1;                            //! Shower shape of isolated photons
-  TH2F *   fhELambda0TRD;                         //! Shower shape of isolated photons, SM behind TRD
-  TH2F *   fhELambda1TRD;                         //! Shower shape of isolated photons, SM behind TRD
-  
+  TH2F *   fhELambda0[2];                         //! Shower shape of (non) isolated photons (do not apply SS cut previously)
+  TH2F *   fhELambda1[2];                         //! Shower shape of (non) isolated photons (do not apply SS cut previously)
+  TH2F *   fhELambda0SSBkg;                       //! Shower shape of non isolated photons close to isolation threshold (do not apply SS cut previously)
+  TH2F *   fhELambda0TRD[2];                      //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
+  TH2F *   fhELambda1TRD[2];                      //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
+  TH2F *   fhELambda0MCPhoton[2];                 //! Shower shape of (non) isolated photon candidates originated by particle photon not decay (do not apply SS cut previously)
+  TH2F *   fhELambda0MCPi0[2];                    //! Shower shape of (non) isolated photon candidates originated by particle 2 merged pi0 photons (do not apply SS cut previously)
+  TH2F *   fhELambda0MCPi0Decay[2];               //! Shower shape of (non) isolated photon candidates originated by particle pi0 decay photon (do not apply SS cut previously)
+  TH2F *   fhELambda0MCEtaDecay[2];               //! Shower shape of (non) isolated photon candidates originated by particle eta decay photon (do not apply SS cut previously)
+  TH2F *   fhELambda0MCOtherDecay[2];             //! Shower shape of (non) isolated photon candidates originated by particle other decay photon (do not apply SS cut previously)
+  TH2F *   fhELambda0MCHadron[2];                 //! Shower shape of (non) isolated photon candidates originated by particle other hadrons (do not apply SS cut previously)
+
   // Local maxima
-  TH2F *   fhNLocMax;                             //! number of maxima in selected clusters
-  TH2F *   fhELambda0LocMax1 ;                    //! E vs lambda0 of selected cluster, 1 local maxima in cluster 
-  TH2F *   fhELambda1LocMax1 ;                    //! E vs lambda1 of selected cluster, 1 local maxima in cluster 
-  TH2F *   fhELambda0LocMax2 ;                    //! E vs lambda0 of selected cluster, 2 local maxima in cluster 
-  TH2F *   fhELambda1LocMax2 ;                    //! E vs lambda1 of selected cluster, 2 local maxima in cluster
-  TH2F *   fhELambda0LocMaxN ;                    //! E vs lambda0 of selected cluster, N>2 local maxima in cluster 
-  TH2F *   fhELambda1LocMaxN ;                    //! E vs lambda1 of selected cluster, N>2 local maxima in cluster 
+  TH2F *   fhNLocMax[2];                          //! number of maxima in selected clusters
+  TH2F *   fhELambda0LocMax1[2] ;                 //! E vs lambda0 of selected cluster, 1 local maxima in cluster 
+  TH2F *   fhELambda1LocMax1[2] ;                 //! E vs lambda1 of selected cluster, 1 local maxima in cluster 
+  TH2F *   fhELambda0LocMax2[2] ;                 //! E vs lambda0 of selected cluster, 2 local maxima in cluster 
+  TH2F *   fhELambda1LocMax2[2] ;                 //! E vs lambda1 of selected cluster, 2 local maxima in cluster
+  TH2F *   fhELambda0LocMaxN[2] ;                 //! E vs lambda0 of selected cluster, N>2 local maxima in cluster 
+  TH2F *   fhELambda1LocMaxN[2] ;                 //! E vs lambda1 of selected cluster, N>2 local maxima in cluster 
   
   //Histograms settings
   Int_t    fHistoNPtSumBins;                      // Number of bins in PtSum histograms
   Float_t  fHistoPtSumMax;                        // PtSum maximum in histogram
-  Float_t  fHistoPtSumMin;                           // PtSum minimum in histogram
+  Float_t  fHistoPtSumMin;                       // PtSum minimum in histogram
   Int_t    fHistoNPtInConeBins;                   // Number of bins in PtInCone histogram
   Float_t  fHistoPtInConeMax;                     // PtInCone maximum in histogram
   Float_t  fHistoPtInConeMin;                     // PtInCone maximum in histogram 
@@ -253,7 +288,7 @@ class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
   AliAnaParticleIsolation(              const AliAnaParticleIsolation & iso) ; // cpy ctor
   AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & iso) ; // cpy assignment
   
-  ClassDef(AliAnaParticleIsolation,13)
+  ClassDef(AliAnaParticleIsolation,16)
 } ;