-#ifndef AliEMCALHistoUtilities_H
-#define AliEMCALHistoUtilities_H
-/* Copyright(c) 1998-2002, ALICE Experiment at CERN, All rights reserved. *
+#ifndef ALIEMCALHISTOUTILITIES_H
+#define ALIEMCALHISTOUTILITIES_H
+/* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
/* $Id$ */
//*-- Authors: J.L. Klay (LLNL) & Aleksei Pavlinov (WSU)
#include <TNamed.h>
+
class TList;
+class TString;
class TH1;
+class TGraph;
+class TGraphErrors;
+class TF1;
+class TLatex;
+class TChain;
+//class TLorentzVector;
+class TArrayF;
+
+//class AliESDCaloCluster;
+//class AliEMCALRecPoint;
+//class AliRunLoader;
class AliEMCALHistoUtilities: public TNamed {
- public:
- AliEMCALHistoUtilities(const char *name="emcalUtilitiesRoutines",
+ public: AliEMCALHistoUtilities(const char *name="emcalUtilitiesRoutines",
const char *tit="EMCAL utility routines");
AliEMCALHistoUtilities(const AliEMCALHistoUtilities &) : TNamed("", ""){
Fatal("cpy ctor", "not implemented") ; }
virtual ~AliEMCALHistoUtilities();
// service routine
- static TList *MoveHistsToList(const char* name="ListOfHists", Bool_t putToBrowser=kTRUE);
- static void FillH1(TList *l=0, Int_t ind=0, Double_t x=-99999., Double_t w=1.);
+ static TList *MoveHistsToList(const char* name="ListOfHists", Bool_t putToBrowser=kTRUE, Bool_t setOwner=kFALSE);
+ static void FillH1(TList *l=0, Int_t ind=0, Double_t x=-99999., Double_t w=1., Double_t error=0.);
static void FillH2(TList *l=0, Int_t ind=0, Double_t x=-99999., Double_t y=-99999., Double_t w=1.);
+ static void FillHProf(TList *l=0, Int_t ind=0, Double_t x=-99999., Double_t y=-99999., Double_t w=1.);
+ static void FillHnSparse(TList *l=0, Int_t ind=0, Double_t* x=0, Double_t w=1.);
static int SaveListOfHists(TList *mylist=0, const char* name="test", Bool_t kSingleKey=kFALSE,
const char* opt="RECREATE");
- static void AddToNameAndTitle(TH1 *h=0, const char *name=0, const char *title=0);
+ static void AddToNameAndTitle(TNamed *h=0, const char *name=0, const char *title=0);
static void AddToNameAndTitleToList(TList *l=0, const char *name=0, const char *title=0);
- //
+ static void ResetListOfHists(TList *l);
+ static void Titles(TH1 *hid=0, const char *titx="",const char *tity="");
+ static TList* CreateProjectionsX(TList *l=0, const Int_t ind=2, const Char_t* name="prox");
+
+ static TLatex *Lat(const char *text="", Float_t x=0.0,Float_t y=0.0, Int_t align=22, Float_t tsize=0.05, short tcolor = 1);
+ static TGraph *DrawGraph(Int_t n=4, Double_t *x=0, Double_t *y=0, Int_t markerColor=4,
+ Int_t markerStyle=4, const char* opt="", const char* tit="", const char* xTit=" jet E_{t} ",
+ const char* yTit="", Int_t ifun=0, const char *optFit="W+", const char *fun="");
+ static TGraphErrors *DrawGraphErrors(const Int_t n=4,Double_t *x=0,Double_t *y=0,Double_t *ex=0,
+ Double_t *ey=0, Int_t markerColor=4,Int_t markerStyle=4, const char* opt="",
+ const char* tit="", const char* xTit=" jet E_{t} ",
+ const char* yTit="", Int_t ifun=0, const char *optFit="W+", const char *fun="");
+ // TChain
+ static void InitChain(TChain *chain=0, const char* nameListOfFiles=0, Int_t nFileMax=0);
+ //static AliRunLoader* InitKinematics(const Int_t nev=0, const char* galiceName="galice.root");
+ //static AliRunLoader* GetRunLoader(const Int_t nev, const Char_t* galiceName,
+ // const Char_t* eventFolderName, AliRunLoader* rlOld);
+ //
+ static Double_t GetMomentum(const char* nameListOfFiles);
static int ParseString(const TString &topt, TObjArray &Opt);
+ // Analysis utilites
+ //static Bool_t GetLorentzVectorFromESDCluster(TLorentzVector &v, const AliESDCaloCluster *cl);
+ //static Bool_t GetLorentzVectorFromRecPoint(TLorentzVector &v, const AliEMCALRecPoint *rp);
+ // Drawing
+ static void DrawHist(TH1* hid=0,int lineWidth=1,int lineColor=1,const char* opt="",int lineStyle=1);
+ // Fitting:
+ static TF1* Gausi(const char *addName, double xmi, double xma, double N, double mean, double sig,
+ double width);
+ static TF1* Gausi(const char *addName, double xmi, double xma, TH1 *h);
+
+ static TF1* GausiPol2(const char *addName, double xmi, double xma, TF1 *g, TF1* bg);
+ //
+ static Double_t Gi(Double_t *x, Double_t *par);
+ static Double_t GiPol2(Double_t *x, Double_t *par);
+ // Calibration stuff
+ static Double_t GetCorrectionCoefficientForGamma1(const Double_t eRec);
+ static Double_t GetCorrectedEnergyForGamma1(const Double_t eRec);
+ static TF1* GetResolutionFunction(const char *opt, TString &latexName);
+ //
+ // Analysis
+ //
+ // Trigger
+ static TList* GetTriggersListOfHists(const Int_t scale=0, const Int_t nTrig=5, const Bool_t toBrowser=kFALSE);
+ static void FillTriggersListOfHists(TList *l=0, TArrayF *triggerPosition=0, TArrayF *triggerAmplitudes=0);
+ // Jet(s) kinematics
+ static TList* GetJetsListOfHists(Int_t njet=2, Bool_t toBrowser=kFALSE);
+ //static void FillJetKineListOfHists(TList *l, AliRunLoader* rl, TLorentzVector &goodJet);
AliEMCALHistoUtilities & operator = (const AliEMCALHistoUtilities &) {
Fatal("operator =", "not implemented") ; return *this ; }
ClassDef(AliEMCALHistoUtilities,1) // EMCAL utility routines
};
-#endif // AliEMCALHistoUtilities_H
+#endif // ALIEMCALHISTOUTILITIES_H