]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG4/PartCorrDep/AliAnaPi0EbE.h
correct binning of hNCellsMod histogram, set remaining TH3 histogram to be filled...
[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / AliAnaPi0EbE.h
index 5a34a946efcf4ee7c7bd04a9e9864967f0408164..ecac2fba274ae17edca8d1c87a8f448d9ee0ea38 100755 (executable)
-#ifndef ALIANAPI0EBE_H\r
-#define ALIANAPI0EBE_H\r
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
- * See cxx source for full Copyright notice     */\r
-/* $Id: AliAnaPi0EbE.h 27413 2008-07-18 13:28:12Z gconesab $ */\r
-\r
-//_________________________________________________________________________\r
-//\r
-// Class for the analysis of high pT pi0 event by event\r
-// Pi0 identified by one of the following:\r
-//  -Invariant mass of 2 cluster in calorimeter\r
-//  -Shower shape analysis in calorimeter\r
-//  -Invariant mass of one cluster in calorimeter and one photon reconstructed in TPC (in near future)\r
-//\r
-//-- Author: Gustavo Conesa (INFN-LNF)  &  Raphaelle Ichou (SUBATECH)\r
-//_________________________________________________________________________\r
-\r
-\r
-// --- ROOT system ---\r
-#include <TH2F.h>\r
-#include <TString.h>\r
-\r
-// --- ANALYSIS system ---\r
-#include "AliAnaPartCorrBaseClass.h"\r
-class AliAODPWG4ParticleCorrelation ;\r
-\r
-class TList ;\r
-\r
-class AliAnaPi0EbE : public AliAnaPartCorrBaseClass {\r
-\r
-public: \r
-\r
-       AliAnaPi0EbE() ; // default ctor\r
-       AliAnaPi0EbE(const AliAnaPi0EbE & g) ; // cpy ctor\r
-       AliAnaPi0EbE & operator = (const AliAnaPi0EbE & g) ;//cpy assignment\r
-       virtual ~AliAnaPi0EbE() ; //virtual dtor\r
-       \r
-       enum anaTypes {kIMCalo, kSSCalo, kIMCaloTracks};\r
-       \r
-       TList *  GetCreateOutputObjects();\r
-       \r
-       void Init();
-       void InitParameters();\r
-       \r
-       void MakeAnalysisFillAOD()  ;\r
-       void MakeAnalysisFillHistograms() ; \r
-       \r
-       void MakeInvMassInCalorimeter() ;\r
-       void MakeInvMassInCalorimeterAndCTS() ;\r
-       void MakeShowerShapeIdentification() ;\r
-       \r
-       void Print(const Option_t * opt)const;\r
-       \r
-       anaTypes GetAnalysisType()   const {return fAnaType ; }\r
-       void SetAnalysisType(anaTypes ana)    {fAnaType = ana ; }\r
-       \r
-       TString GetInputAODGammaConvName()   const {return fInputAODGammaConvName ; }\r
-       void SetInputAODGammaConvName(TString name)    {fInputAODGammaConvName = name ; }       \r
-       \r
-       //Only for pi0 SS identification case\r
-       void SetCalorimeter(TString det)    {fCalorimeter = det ; }\r
-\r
-       void SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) {\r
-               fMinDist = m1;\r
-               fMinDist2 = m2;\r
-               fMinDist3 = m3;\r
-       }\r
-       \r
-private:\r
\r
-       anaTypes fAnaType; //Select analysis type\r
-       \r
-       //Only for pi0 SS identification case, kSSCalo\r
-       TString fCalorimeter ; // Calorimeter where the gamma is searched;\r
-       Float_t fMinDist ;     // Minimal distance to bad channel to accept cluster\r
-       Float_t fMinDist2;     // Cuts on Minimal distance to study acceptance evaluation\r
-       Float_t fMinDist3;     // One more cut on distance used for acceptance-efficiency study\r
-       \r
-       //Only for combination of calorimeter and conversion photons, kIMCaloTracks\r
-       TClonesArray * fInputAODGammaConv; //AOD array with conversion photons reconstructed in CTS\r
-       TString fInputAODGammaConvName;    //Name of AOD branch with conversion photons\r
-       \r
-       //Histograms\r
-       TH1F * fhPtPi0   ;         //! Number of identified  pi0\r
-       TH2F * fhPhiPi0  ;         //! Phi of identified  pi0\r
-       TH2F * fhEtaPi0  ;         //! eta of identified  pi0  \r
-       \r
-       //MC\r
-       TH1F * fhPtMCNoPi0;   //! Number of identified pi0, not coming from pi0\r
-       TH2F * fhPhiMCNoPi0;  //! Phi of identified pi0, not coming from pi0\r
-       TH2F * fhEtaMCNoPi0;  //! eta of identified  pi0, not coming from pi0\r
-       TH1F * fhPtMCPi0;     //! Number of identified pi0, coming from pi0\r
-       TH2F * fhPhiMCPi0;    //! Phi of identified pi0, coming from pi0\r
-       TH2F * fhEtaMCPi0;    //! eta of identified pi0, coming from pi0\r
-       \r
-       ClassDef(AliAnaPi0EbE,1)\r
-} ;\r
\r
-\r
-#endif //ALIANAPI0EBE_H\r
-\r
-\r
-\r
+#ifndef ALIANAPI0EBE_H
+#define ALIANAPI0EBE_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice     */
+/* $Id: AliAnaPi0EbE.h 27413 2008-07-18 13:28:12Z gconesab $ */
+
+//_________________________________________________________________________
+//
+// Class for the analysis of high pT pi0 event by event
+// Pi0 identified by one of the following:
+//  -Invariant mass of 2 cluster in calorimeter
+//  -Shower shape analysis in calorimeter
+//  -Invariant mass of one cluster in calorimeter and one photon reconstructed in TPC (in near future)
+//
+//-- Author: Gustavo Conesa (INFN-LNF)  &  Raphaelle Ichou (SUBATECH)
+//_________________________________________________________________________
+
+
+// --- ROOT system ---
+class TH3F ; 
+class TList ;
+class TObjString;
+
+// --- ANALYSIS system ---
+#include "AliAnaPartCorrBaseClass.h"
+
+class AliAnaPi0EbE : public AliAnaPartCorrBaseClass {
+
+ public: 
+  AliAnaPi0EbE() ; // default ctor
+  virtual ~AliAnaPi0EbE() ; //virtual dtor
+ private:
+  AliAnaPi0EbE(const AliAnaPi0EbE & g) ; // cpy ctor
+  AliAnaPi0EbE & operator = (const AliAnaPi0EbE & g) ;//cpy assignment
+
+ public:  
+       
+  //General
+  
+  TObjString *   GetAnalysisCuts();
+  
+  TList      *   GetCreateOutputObjects();
+  
+  void           Init();
+  
+  void           InitParameters();
+
+  void           MakeAnalysisFillAOD()  ;
+   
+  void           MakeAnalysisFillHistograms() ; 
+  
+  void           Print(const Option_t * opt) const;
+  
+  // Main
+  
+  void           MakeInvMassInCalorimeter() ;
+  
+  void           MakeInvMassInCalorimeterAndCTS() ;
+  
+  void           MakeShowerShapeIdentification() ;
+  
+  //Setters Getters
+  
+  //Analysis types
+  enum anaTypes  {kIMCalo, kSSCalo, kIMCaloTracks};  
+  anaTypes       GetAnalysisType()                     const { return fAnaType               ; }
+  void           SetAnalysisType(anaTypes ana)               { fAnaType = ana                ; }
+  
+  TString        GetInputAODGammaConvName()            const { return fInputAODGammaConvName ; }
+  void           SetInputAODGammaConvName(TString name)      { fInputAODGammaConvName = name ; }       
+  
+  //Only for pi0 SS identification case
+  void           SetCalorimeter(TString & det)               { fCalorimeter = det            ; }
+  
+  void           SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) {
+                  fMinDist = m1; fMinDist2 = m2; fMinDist3 = m3                              ; }
+
+  //For histograms
+  enum mcTypes   { mcPhoton = 0, mcConversion = 1, mcPi0    = 2,  
+                   mcEta    = 3, mcElectron   = 4, mcHadron = 5 };
+
+ private:
+  
+  anaTypes       fAnaType; //Select analysis type
+  
+  //Only for pi0 SS identification case, kSSCalo
+  TString        fCalorimeter ;            // Calorimeter where the gamma is searched;
+  Float_t        fMinDist ;                // Minimal distance to bad channel to accept cluster
+  Float_t        fMinDist2;                // Cuts on Minimal distance to study acceptance evaluation
+  Float_t        fMinDist3;                // One more cut on distance used for acceptance-efficiency study
+  
+  //Only for combination of calorimeter and conversion photons, kIMCaloTracks
+  TClonesArray * fInputAODGammaConv;       //! AOD array with conversion photons reconstructed in CTS
+  TString        fInputAODGammaConvName;   //  Name of AOD branch with conversion photons
+  
+  //Histograms
+  
+  TH1F         * fhPtPi0  ;                //! Number of identified  pi0 vs pT
+  TH1F         * fhEPi0   ;                //! Number of identified  pi0 vs E
+  TH3F         * fhEEtaPhiPi0  ;           //! E vs eta phi of identified  pi0 
+  
+  TH2F         * fhEDispersion ;           //! E vs disp of pi0 pairs
+  TH2F         * fhELambda0 ;              //! E vs lambda0 of pi0 pairs 
+  TH2F         * fhELambda1 ;              //! E vs lambda1 of pi0 pairs 
+  TH2F         * fhELambda0NoTRD ;         //! E vs lambda0 of pi0 pairs, not behind TRD 
+  TH2F         * fhELambda0FracMaxCellCut ;//! E vs lambda0 of pi0 pairs, fraction of cluster energy in max cell cut 
+  TH2F         * fhEFracMaxCell ;          //! E vs frac max cell of cluster  
+  TH2F         * fhEFracMaxCellNoTRD ;     //! E vs frac max cell of cluster, not behind TRD  
+
+  TH2F         * fhClusterPairDiffTimeE;   //! Pair of clusters time difference vs E
+  TH2F         * fhClusterPairDiffTimeAsy; //! Pair of clusters time difference vs Asymmetry
+  
+  //MC histograms
+  
+  TH2F         * fhEMCLambda0[6] ;         //! E vs lambda0 of pi0 pairs but really from MC particle
+  TH2F         * fhEMCLambda1[6] ;         //! E vs lambda1 of pi0 pairs but really from MC particle
+  TH2F         * fhEMCDispersion[6] ;      //! E vs dispersion of pi0 pairs but really from MC particle
+  TH2F         * fhEMCLambda0NoTRD[6] ;         //! E vs lambda0 of pi0 pairs but really from MC particle, not behind TRD
+  TH2F         * fhEMCLambda0FracMaxCellCut[6] ;//! E vs lambda0 of pi0 pairs but really from MC particle, fraction of cluster energy in max cell cut
+  TH2F         * fhEMCFracMaxCell[6] ;     //! E vs fraction of max cell 
+  
+  TH1F         * fhPtMCNoPi0;              //! Number of identified pi0, not coming from pi0
+  TH2F         * fhPhiMCNoPi0;             //! Phi of identified pi0, not coming from pi0
+  TH2F         * fhEtaMCNoPi0;             //! eta of identified  pi0, not coming from pi0
+  TH1F         * fhPtMCPi0;                //! Number of identified pi0, coming from pi0
+  TH2F         * fhPhiMCPi0;               //! Phi of identified pi0, coming from pi0
+  TH2F         * fhEtaMCPi0;               //! eta of identified pi0, coming from pi0
+  
+  ClassDef(AliAnaPi0EbE,7)
+} ;
+
+
+#endif //ALIANAPI0EBE_H
+
+
+