added pi0 calibration, linearity, shower profile

[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecPointsQaESDSelector.h

diff --git a/EMCAL/AliEMCALRecPointsQaESDSelector.h b/EMCAL/AliEMCALRecPointsQaESDSelector.h
index 9837810324a38932e483139d0892915e563488b2..a447ee5d0e519d3abca87072196d397b8101366c 100644 (file)
--- a/EMCAL/AliEMCALRecPointsQaESDSelector.h
+++ b/EMCAL/AliEMCALRecPointsQaESDSelector.h
@@ -9,60 +9,132 @@
 
 #include "AliSelector.h"
 
+#include <TObjArray.h>
+
 class AliEMCALFolder;
+class AliRunLoader;
+class AliEMCALRecPoint;
+class AliEMCALCellInfo;
+class cellInfo;
 
 class TList;
+class TCanvas;
+class TH1F;
+class TH2F;
 class TBrowser;
 class TChain;
-class TObjectSet;
+class TFolder;
+class TArrayI;
 
 class AliEMCALRecPointsQaESDSelector :  public AliSelector {
   public:
     AliEMCALRecPointsQaESDSelector();
     virtual ~AliEMCALRecPointsQaESDSelector();
 
-    virtual void    Begin(TTree* tree);
-    virtual void    Init(TTree* tree);
-    virtual void    SlaveBegin(TTree *tree);
+    virtual void    Begin(TTree*);
+    virtual void    SlaveBegin(TTree* tree);
+    virtual void    Init(TTree *tree);
+    virtual Bool_t  Notify();
     virtual Bool_t  Process(Long64_t entry);
     virtual void    SlaveTerminate();
     virtual void    Terminate();
   //
   void InitStructure(Int_t it);
-  static TList *DefineHistsOfRP(const char *name="RP");
-  // static TList *DefineHistsOfTowers(const char *name="towers");
+  static TList *DefineHistsOfRP(const char *name="RP",  Double_t p=110.0, Int_t keyOpt=0);
+  static TList *DefineHistsOfKineVsRP(const char *name="KineVsRP",  Double_t p=110.0, Int_t keyOpt=0);
+  static TList *DefineHistsForShowerProfile(const char *name="ProfY", Double_t p=1.);
+  static void   FillHistsOfKineVsRP(TList *l, AliRunLoader* RL, TClonesArray &lvM);
+  static void   FillHistsForShowerProfile(TList *l, AliEMCALRecPoint *rp, AliEMCALCellInfo* t);
+  static void   EvalLocalPhiPosition(const Double_t wlog, const AliEMCALRecPoint *rp, const AliEMCALCellInfo* t, Double_t &xcog, Int_t &phiSize, cellInfo &rMax); 
+  //
+  TList *GetListKineVsRP()      {return fLKineVsRP;}
+  TList *GetListShowerProfile() {return fLShowerProfile;}
+  //  static TList *DefineHistsOfTowers(const char *name="towers");
   // 
-  void FitEffMassHist(); //*MENU*  
+  void FitEffMassHist(); // *MENU*  
   void PrintInfo();      // *MENU*
   //
   void    SetChain(TChain *chain)  {fChain = chain;}
   TChain* GetChain()               {return fChain;}
+  void    SetMomentum(Double_t p);
+  Double_t GetMomentum() {return fPmom;}
 
   AliEMCALFolder* CreateEmcalFolder(const Int_t it);
   void SetEmcalFolder(AliEMCALFolder* folder); 
   void SetEmcalOldFolder(AliEMCALFolder* folder); 
-  AliEMCALFolder* GetEmcalFolder() {return fEMCAL;}
-  AliEMCALFolder* GetEmcalOldFolder() {return fEMCALOld;}
   AliEMCALFolder* GetEmcalOldFolder(const Int_t nsm);
-  // 
+  AliEMCALCellInfo* GetCellsInfo() {return fCellsInfo;}
+  //
+  void      SetStringOfRunOpts(const char *st) {fRunOpts = st;}
+  TObjArray GetOptsArray() const {return fArrOpts;}
+  Int_t     GetKeyOptsValue(Int_t key);  // *MENU*
+  void      CheckRunOpts();
+  //
   virtual void Browse(TBrowser* b);
   virtual Bool_t  IsFolder() const;
+  //
+  void    Save(Int_t ver=0, const char *optIO="NEW");   // *MENU*
+  static  AliEMCALRecPointsQaESDSelector* ReadSelector(const char* nf = "/home/pavlinov/ALICE/SHISHKEBAB/RF/CALIB/PROF/PROFILE_0.root");
+  
 
   //
   //// Pictures staf - Jun 26, 2007
   //
   void ReadAllEmcalFolders();
   void PictVsIterNumber(const Int_t ind=0, const Int_t nsm=0);
-
+  // Gamma
+  TH1F* FitHistOfRecPointEnergy(const char *opt="CLONE");
+  static TCanvas *Linearity(TList *l, Int_t ifun=3);
+  static TCanvas *DrawKineVsRP(TList *l);
+  // Profile
+  static TCanvas *DrawMeanVsLog(TH2F *h2);
+  // Geometry staff
+  TCanvas *DrawPhiEtaAnglesDistribution(const char *gn="SHISH_TRD1_CURRENT_2X2"); // *MENU*
+  // Geometry constants 
+  TCanvas *DrawDeffVsEnergy();  // *MENU*
+  TCanvas *DrawDeffVsEnergy2(const char *opt="fit1"); // *MENU*
+  void     ReadParsDeffAndW0(const char *dirName="/data/r22b/ALICE/CALIB/FIT/",
+          double *deff=0, double *edeff=0, double *w0=0, double *ew0=0, const Int_t pri=0);
+  TCanvas *DrawSpaceResolution();
  protected:
+  Double_t fPmom; // positive if defined
   //
   TChain* fChain; //! chain if ESD files
-  TList* fLofHistsPC; //! list of histograms of pseudo clusters 
-  TList* fLofHistsRP; //! list of histograms of rec.points 
+  TList* fLofHistsPC; // list of histograms of pseudo clusters 
+  TList* fLofHistsRP; // list of histograms of rec.points 
+  TList* fLKineVsRP;  // list of histograms kinematics vs rec.points 
+  TList* fLShowerProfile;  // list of histograms for shower profile business
+  //
+  AliEMCALCellInfo *fCellsInfo; //
+  TFolder*         fEmcalPool;  //
+  //
+  // Options - Jul 10, 2007
   //
-  TObjectSet*     fEmcalPool;
-  AliEMCALFolder* fEMCAL;    //! current  EMCAL object
-  AliEMCALFolder* fEMCALOld; //! previous EMCAL object
+  TString   fRunOpts;        // String of running options
+  TObjArray fArrOpts;        // Array of options 
+  // Options keys
+  TArrayI  *fKeyOpts;
+  // Static parameters
+ public:
+  static AliEMCALFolder* GetEmcalFolder() {return fEMCAL;}
+  static AliEMCALFolder* GetEmcalOldFolder() {return fEMCALOld;}
+  static AliEMCALFolder*  fEMCAL;      // current  EMCAL object
+  static AliEMCALFolder*  fEMCALOld;   // previous EMCAL object
+  //
+  static Double_t fDistEff;  // effective depth of electromagnetic shower
+  static Double_t fW0;       // 
+  static Double_t fSlopePhiShift; // phi shift of cluster = fSlopePhiShift * phi
+  static void SetFitParameters(Double_t deff, Double_t w0, Double_t slope) 
+  {
+    fDistEff = deff; fW0 = w0; fSlopePhiShift = slope;
+  }
+  static void GetFitParameters(Double_t &deff, Double_t &w0, Double_t &slope)
+  {
+    deff = fDistEff; w0 = fW0; slope = fSlopePhiShift;
+  } 
+  void ResetAllListOfHists();
+  void ReloadChain(Long64_t entry=0);
+  void GetInitialParsForFit(const Int_t var, Double_t &deff, Double_t &w0, Double_t &phislope, const int phiCase=0);
 
  private:
   AliEMCALRecPointsQaESDSelector(const AliEMCALRecPointsQaESDSelector&);
@@ -70,5 +142,4 @@ class AliEMCALRecPointsQaESDSelector :  public AliSelector {
 
   ClassDef(AliEMCALRecPointsQaESDSelector, 1);
 };
-
 #endif