Adding statistical function (Marian)
[u/mrichter/AliRoot.git] / PHOS / AliPHOSGammaJet.h
index 9fb3a5c..32ce4ca 100644 (file)
@@ -4,6 +4,20 @@
  * See cxx source for full Copyright notice     */
 /* $Id$ */
 
+/* History of cvs commits:
+ *
+ * $Log$
+ * Revision 1.11  2006/01/12 17:01:00  schutz
+ * Correct few coding violations, missing comments
+ *
+ * Revision 1.9  2005/12/20 07:31:51  schutz
+ * added data members
+ *
+ * Revision 1.7  2005/05/28 14:19:04  schutz
+ * Compilation warnings fixed by T.P.
+ *
+ */
+
 //_________________________________________________________________________
 //  Class for the analysis of gamma-jet correlations.     
 //  Basically it seaches for a prompt photon in the PHOS acceptance, 
 #include "TTask.h"
 #include "TRandom.h"
 #include "TArrayD.h"
-class AliPHOSGeometry ;
+#include "AliESD.h"
+#include "TChain.h"
+
 class AliPHOSFastGlobalReconstruction ;
+
 //#include "../PYTHIA6/AliGenPythia.h"
+
 // --- AliRoot header files ---
 
 class AliPHOSGammaJet : public TTask {
@@ -29,8 +47,7 @@ public:
   AliPHOSGammaJet() ; // default ctor
   AliPHOSGammaJet(const TString inputfilename) ; //ctor 
   AliPHOSGammaJet(const AliPHOSGammaJet & gj) ; // cpy ctor
-  AliPHOSGammaJet & operator = (const AliPHOSGammaJet & /*rvalue*/) 
-    { return *this ;} //assignement operator requested by coding convention but not needed
+
   virtual ~AliPHOSGammaJet() ; //virtual dtor
   virtual void   Exec(Option_t *option); 
   void List() const; 
@@ -62,6 +79,11 @@ public:
    TString  GetConeName(Int_t i) const {return fNameCones[i] ; }
    TString  GetPtThresName(Int_t i) const {return fNamePtThres[i] ; }
    Bool_t   GetTPCCutsLikeEMCAL() const {return fTPCCutsLikeEMCAL ; }
+   Bool_t   IsESDdata() const {return fESDdata ; }
+
+   TString  GetDirName() const {return fDirName ; }
+   TString  GetESDTreeName() const {return fESDTree ; }
+   char *   GetDirPattern()  const {return fPattern ; }
 
    Bool_t   IsAnyConeOrPt() const {return fAnyConeOrPt ; }
    Bool_t   IsFastReconstruction() const {return fOptFast ; }
@@ -69,7 +91,7 @@ public:
    Bool_t   IsOnlyCharged() const {return fOnlyCharged ; }
 
   void Plot(TString what="all", Option_t *option="") const;
-  void Print(char * opt);
+  void Print(const Option_t * opt)const;
 
   void SetAngleMaxParam(Int_t i, Double_t par)
   {fAngleMaxParam.AddAt(par,i) ; }
@@ -112,16 +134,24 @@ public:
   void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
   {fRatioMaxCut = ratiomax;  fRatioMinCut = ratiomin;}
   void SetTPCCutsLikeEMCAL(Bool_t b){ fTPCCutsLikeEMCAL= b ; }
+  void SetESDdata(Bool_t b){ fESDdata= b ; }
  
+  void  SetDirName(TString dn) { fDirName  = dn ; }
+  void  SetESDTreeName(TString en) { fESDTree = en ; }
+  void  SetDirPattern(char * dp)  { fPattern = dp ; }
 
  private:
+  AliPHOSGammaJet & operator = (const AliPHOSGammaJet & /*rvalue*/) 
+  { return *this ;}
+
+
 //   void AddHIJINGToList(TList & particleList, TList & particleListCh, 
 //                    TList & particleListNe, const Int_t iEvent, 
 //                    const TLorentzVector gamma, Double_t & rot ); 
  
   void AddHIJINGToList(Int_t iEvent, TClonesArray * particleList, 
                       TClonesArray * plCh, TClonesArray * plNe, 
-                      TClonesArray * plNePHOS, const AliPHOSGeometry * geom); 
+                      TClonesArray * plNePHOS); 
 
 
   Double_t CalculateJetRatioLimit(const Double_t ptg, const Double_t *param, 
@@ -129,8 +159,11 @@ public:
 
   void CreateParticleList(Int_t iEvent, TClonesArray * particleList, 
                          TClonesArray * plCh, TClonesArray * plNe, 
-                         TClonesArray * plNePHOS, 
-                         const AliPHOSGeometry * geom );
+                         TClonesArray * plNePHOS); 
+  void CreateParticleListFromESD(TClonesArray * particleList, 
+                                TClonesArray * plCh, TClonesArray * plNe, 
+                                TClonesArray * plNePHOS,  
+                                const AliESD * esd );//, Int_t iEvent); 
   
   void FillJetHistos(TClonesArray * pl, Double_t ptg, TString conf, TString type);
 
@@ -145,7 +178,7 @@ public:
               Double_t ptl, Double_t phil, Double_t etal, 
               TString  type, TLorentzVector & jet); 
   void MakeJetAnyConeOrPt(TClonesArray * particleList, Double_t ptg, 
-                         Double_t phil, Double_t ptl, Double_t phil, 
+                         Double_t phig, Double_t ptl, Double_t phil, 
                          Double_t etal, TString  type); 
   void GetGammaJet(TClonesArray * pl,  Double_t &pt, 
                   Double_t &phi, Double_t &eta, Bool_t &Is)  const;
@@ -165,9 +198,19 @@ public:
  
   void Pi0Decay(Double_t mPi0, TLorentzVector &p0, 
                TLorentzVector &p1, TLorentzVector &p2, Double_t &angle) ;
+
+  TChain * ReadESDfromdisk(const UInt_t eventsToRead,
+                             const TString dirName, 
+                             const TString esdTreeName = "esdTree",
+                             const char *  pattern     = "Evt") ;
+  TChain * ReadESD(const UInt_t eventsToRead = 1,
+                     TString fileName  = "",
+                     const TString esdTreeName = "esdTree",
+                     const char *  pattern     = "Evt" ) ;
+  
   Double_t SigmaE(Double_t energy) ;
   Double_t SigmaP(Double_t energy) ;
-
+  
   void SetJet(TParticle * part, Bool_t & b, Float_t cone, Double_t eta, 
              Double_t phi);
 
@@ -180,6 +223,7 @@ public:
   TString    fInputFileName;   //!
   TString    fHIJINGFileName;  //!
   Bool_t     fHIJING;          // Add HIJING event to PYTHIA event?
+  Bool_t     fESDdata ;        // Read ESD?      
   Double_t   fEtaCut ;         // Eta cut
   Bool_t     fOnlyCharged ;    // Only jets of charged particles
   Double_t   fPhiEMCALCut[2] ; // Phi cut maximum
@@ -195,6 +239,11 @@ public:
   Double_t   fRatioMinCut ;    // Leading particle/gamma Ratio cut minimum
   Bool_t     fTPCCutsLikeEMCAL ; //Same jet energy ratio limits for both conf.
 
+  //Read ESD Paramenters
+  TString fDirName ; //Name of the directory where data is
+  TString fESDTree ; //Name of the ESD Tree
+  char *  fPattern ; //Pattern followed by  directory data
+
   //Jet selection parameters
   //Fixed cuts (old)
   Double_t   fJetTPCRatioMaxCut ; // Leading particle/gamma Ratio cut maximum
@@ -237,7 +286,7 @@ public:
   TArrayD    fAngleMaxParam ; //Max opening angle selection parameters
   Bool_t fSelect  ;  //Select jet within limits
 
-  ClassDef(AliPHOSGammaJet,2)
+  ClassDef(AliPHOSGammaJet,3)
 } ;