Removed old jetfinder
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALJetMicroDst.h
index 55bf53f97ecec21f078b8d1bbe95f83381ea9a35..49642b73755b1385274fe6390ea2a6108940f1a7 100644 (file)
@@ -9,10 +9,12 @@
 //                  
 //*-- Author: Aleksei Pavlinov (WSU)
 #include <TNamed.h>
-#include <TFile.h>
-#include <TTree.h>
-#include <TH1.h>
-#include <TH2.h>
+class TFile;
+class TTree;
+class TH1;
+class TH1F;
+class TH2F;
+class TH2;
 
 class AliGenHijingEventHeader;
 class AliRun;
@@ -22,49 +24,12 @@ class TBrowser;
 
 class AliEMCALJetMicroDst: public TNamed {
 
-  private:
-  Int_t   fDebug;
-
-  public:
-  TFile*  fFile;
-  TTree*  fTree;
-  TString fName;
-  TList*  fListHist;    //!
-  TString fFileName;    // for convenience
-
-  Float_t decone;   //! for EMCAL
-  Float_t ptcone;   //! for ch.particles 
-  // For partons after hard scattering
-  Int_t   npart;
-  Float_t xpt[4];  //[npart]
-  Float_t xeta[4]; //[npart]
-  Float_t xphi[4]; //[npart]
-  // Jet 
-  Int_t   njet;
-  Float_t jet[10];   //[njet]
-  Float_t jetal[10]; //[njet]
-  Float_t jphil[10]; //[njet]
-  Float_t jetaw[10]; //[njet]
-  Float_t jphiw[10]; //[njet]
-  // Charge particle in jet ??
-  // eT in EMCAL itself - 24-jan-2003
-  Int_t   ncell;         // 96*144 =13824 
-  Int_t   idcell[13824]; //[ncell]
-  Float_t etcell[13824]; //[ncell] : de = det*sf
-  // eT in EMCAL grid for jet finder 
-  Int_t   ngrid;         // 96*144 =13824 
-  Int_t   idgrid[13824]; //[ngrid]
-  Float_t etgrid[13824]; //[ngrid]
-  // charge particle which hit to EMCAL - 28-jan-2003
-  Int_t   nchp;
-  Int_t   pid[20000];  //[nchp]
-  Float_t ppt[20000];  //[nchp]
-  Float_t peta[20000]; //[nchp]
-  Float_t pphi[20000]; //[nchp]
 
   public:
   AliEMCALJetMicroDst(char *name="jetMicroDst",
   char *tit="jet Micro Dst for preparation of proposal");
+  AliEMCALJetMicroDst(const  AliEMCALJetMicroDst &) : TNamed("", ""){
+    Fatal("cpy ctor", "not implemented") ; }
   virtual ~AliEMCALJetMicroDst();
   Bool_t  Create(TFile *file);
   Bool_t  Create(const char  *fname);
@@ -86,38 +51,94 @@ class AliEMCALJetMicroDst: public TNamed {
   void    Print(Option_t* option="") const;                          // *MENU* 
   Int_t   GetEntry(Int_t entry);
   void    Test();
-  Int_t   GetNpart() {return npart;}
-  Bool_t  GetParton(Int_t i, Float_t& pt, Float_t& eta, Float_t& phi);
-  Bool_t  GetParton(Int_t i, TVector3& vec);
-  Int_t   GetNjet() {return njet;} 
-  Bool_t  GetJet(Int_t i,Int_t mode, Float_t& pt,Float_t& eta,Float_t& phi);
-  Bool_t  GetJet(Int_t i,Int_t mode, TVector3& vec);
+  Int_t   GetNpart() const {return fnpart;}
+  Bool_t  GetParton(Int_t i, Float_t& pt, Float_t& eta, Float_t& phi) const ;
+  Bool_t  GetParton(Int_t i, TVector3& vec) const ;
+  Int_t   GetNjet() const {return fnjet;} 
+  Bool_t  GetJet(Int_t i,Int_t mode, Float_t& pt,Float_t& eta,Float_t& phi) const ;
+  Bool_t  GetJet(Int_t i,Int_t mode, TVector3& vec) const ;
   static  void FillVector(Float_t pt, Float_t eta, Float_t phi, TVector3& vec);
-  void    GetEtaPhi(Int_t id, Double_t &eta, Double_t &phi);
-  TVector3& GetCellVector(Int_t i);
-  TVector3& GetGridVector(Int_t i);
+  void    GetEtaPhi(Int_t id, Double_t &eta, Double_t &phi) const ;
+  TVector3& GetCellVector(Int_t i) const ;
+  TVector3& GetGridVector(Int_t i) const ;
   // 13-apr-2003
-  Double_t GetSumInCone(TVector3 &jet, Int_t nc, Float_t *et,Float_t *eta,Float_t *phi, Double_t cellEtCut, Double_t rJet);
-  Double_t GetEmcalEtInCone(TVector3 &jet, Double_t cellEtCut=0.0, Double_t rJet=0.5);
-  Double_t GetTpcPtInCone(TVector3 &jet, Double_t cellEtCut=0.0, Double_t rJet=0.5);
-  Double_t GetSum(Int_t n, Float_t *ar, Double_t cut=0.0);
-  Double_t GetSumEmcal(Double_t cut=0.0) {return GetSum(ncell, etcell, cut);}
-  Double_t GetSumTpc(Double_t cut=0.0) {return GetSum(nchp, ppt, cut);}
+  Double_t GetSumInCone(TVector3 &jet, Int_t nc, Float_t *et,Float_t *eta,Float_t *phi, Double_t cellEtCut, Double_t rJet) const ;
+  Double_t GetEmcalEtInCone(TVector3 &jet, Double_t cellEtCut=0.0, Double_t rJet=0.5) ;
+  Double_t GetTpcPtInCone(TVector3 &jet, Double_t cellEtCut=0.0, Double_t rJet=0.5) ;
+  Double_t GetSum(Int_t n, Float_t *ar, Double_t cut=0.0) const ;
+  Double_t GetSumEmcal(Double_t cut=0.0) {return GetSum(fncell, fetcell, cut);}
+  Double_t GetSumTpc(Double_t cut=0.0) {return GetSum(fnchp, fppt, cut);}
 
   void    SetDebug(Int_t flag) {fDebug = flag;}
   Float_t GetDebug() const  {return fDebug;}
 
-  TTree* GetTree() {return fTree;}
-  TFile* GetFile() {return fFile;}
+  TTree* GetTree() const {return fTree;}
+  TFile* GetFile() const {return fFile;}
   void   Close();
 
-  Bool_t  IsPythiaDst();
+  Bool_t  IsPythiaDst() const ;
   virtual Bool_t  IsFolder() const;
-  virtual void Browse(TBrowser* b);
+  virtual void Browse(TBrowser* b) const ;
 
   static TList *MoveHistsToList(char* name="List of Hist", Bool_t putToBrowser=kTRUE);
 
-  ClassDef(AliEMCALJetMicroDst,1) // Micro Dst for jet analysis
+  AliEMCALJetMicroDst & operator = (const AliEMCALJetMicroDst &) {
+    Fatal("operator =", "not implemented") ; return *this ; }
+  
+  private:
+  
+  Float_t fpphi[20000]; //[nchp]
+  Int_t   fDebug;      // debug flag
+  TFile*  fFile;       // filename
+  TTree*  fTree;       // Tree pointer
+  TList*  fListHist;    //!
+  TString fFileName;    // for convenience
+
+  Float_t fdecone;   //! for EMCAL
+  Float_t fptcone;   //! for ch.particles 
+  // For partons after hard scattering
+  Int_t   fnpart;      //npartons
+  Float_t fxpt[4];  //[npart]
+  Float_t fxeta[4]; //[npart]
+  Float_t fxphi[4]; //[npart]
+  // Jet 
+  Int_t   fnjet;       // number of jets
+  Float_t fjet[10];   //[njet]
+  Float_t fjetal[10]; //[njet]
+  Float_t fjphil[10]; //[njet]
+  Float_t fjetaw[10]; //[njet]
+  Float_t fjphiw[10]; //[njet]
+  // Charge particle in jet ??
+  // eT in EMCAL itself - 24-jan-2003
+  Int_t   fncell;         // 96*144 =13824 
+  Int_t   fidcell[13824]; //[ncell]
+  Float_t fetcell[13824]; //[ncell] : de = det*sf
+  // eT in EMCAL grid for jet finder 
+  Int_t   fngrid;         // 96*144 =13824 
+  Int_t   fidgrid[13824]; //[ngrid]
+  Float_t fetgrid[13824]; //[ngrid]
+  // charge particle which hit to EMCAL - 28-jan-2003
+  Int_t   fnchp;       //number of charged particles
+  Int_t   fpid[20000];  //[nchp]
+  Float_t fppt[20000];  //[nchp]
+  Float_t fpeta[20000]; //[nchp]
+  TH1F*  fhPtPart; //hist
+  TH1F*         fhNJet;        //hist
+  TH1F*  fhPtJet;      //hist
+  TH2F*  fhEtaPhiPart; //hist
+  TH2F*  fhEtaPhiJet;  //hist
+  TH1F*  fhNcell;      //hist
+  TH1F*         fhCellId;      //hist
+  TH1F*         fhCellEt;      //hist
+       
+  TH1F*         fhSumEt;       //hist
+  TH1F*  fhNgrid;      //hist
+  TH1F*         fhGridId;      //hist
+  TH1F*  fhGridEt;     //hist
+  TH1F*  fhSumEtGrForJF;       //hist
+
+
+  ClassDef(AliEMCALJetMicroDst,3) // Micro Dst for jet analysis
 };
 
 #endif // AliEMCALJETMICRODST_H