[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisHadEtCorrections.h

//Created by Christine Nattrass, Rebecca Scott, Irakli Martashvili
//University of Tennessee at Knoxville
// This is a container class for the correction factors for the hadronic 
// component of transverse energy
// It is filled by the output of AliAnalysisTaskHadEt from spinning over Monte 
// Carlo data (using AliAnalysisHadEtMonteCarlo)
// It is used by AliAnalysisTaskHadEt while spinning over reconstructed data 
// (using AliAnalysisHadEtReconstructed)
//Please see https://twiki.cern.ch/twiki/bin/view/ALICE/ETCaloAnalysis
#ifndef ALIANALYSISHADETCORRECTIONS_H
#define ALIANALYSISHADETCORRECTIONS_H

#include "TString.h"
#include "TNamed.h"
class Rtypes;
class TH1D;
class TNamed;
class TObjArray;

class AliAnalysisHadEtCorrections : public TNamed
{
public:
   
  AliAnalysisHadEtCorrections();
    virtual ~AliAnalysisHadEtCorrections();
    virtual void Init();


    Float_t GetEtaCut() const {return fEtaCut;}
    Float_t GetAcceptanceCorrectionFull() const {return fAcceptanceCorrectionFull;}
    Float_t GetAcceptanceCorrectionEMCAL() const {return fAcceptanceCorrectionEMCAL;}
    Float_t GetAcceptanceCorrectionPHOS() const {return fAcceptanceCorrectionPHOS;}
    Float_t GetNeutralCorrection() const {return fNeutralCorrection;}
    Float_t GetNotHadronicCorrection() const {return fNotHadronicCorrection;}
    Float_t GetpTCutCorrectionTPC() const {return fpTcutCorrectionTPC;}
    Float_t GetpTCutCorrectionITS() const {return fpTcutCorrectionITS;}
    Float_t GetNotIDConstCorrectionTPC() const {return fNotIDConstTPC;}
    Float_t GetNotIDConstCorrectionITS() const {return fNotIDConstITS;}
    Float_t GetNotIDConstCorrectionTPCNoID() const {return fNotIDConstTPCNoID;}
    Float_t GetNotIDConstCorrectionITSNoID() const {return fNotIDConstITSNoID;}
    Float_t GetNeutralCorrectionLowBound() const {return fNeutralCorrectionLow;}
    Float_t GetNotHadronicCorrectionLowBound() const {return fNotHadronicCorrectionLow;}
    Float_t GetpTCutCorrectionTPCLowBound() const {return ffpTcutCorrectionTPCLow;}
    Float_t GetpTCutCorrectionITSLowBound() const {return ffpTcutCorrectionITSLow;}
    Float_t GetNeutralCorrectionHighBound() const {return fNeutralCorrectionHigh;}
    Float_t GetNotHadronicCorrectionHighBound() const {return fNotHadronicCorrectionHigh;}
    Float_t GetpTCutCorrectionTPCHighBound() const {return ffpTcutCorrectionTPCHigh;}
    Float_t GetpTCutCorrectionITSHighBound() const {return ffpTcutCorrectionITSHigh;}
    Float_t GetNotIDConstCorrectionTPCLowBound() const {return fNotIDConstTPCLow;}
    Float_t GetNotIDConstCorrectionITSLowBound() const {return fNotIDConstITSLow;}
    Float_t GetNotIDConstCorrectionTPCNoIDLowBound() const {return fNotIDConstTPCNoIDLow;}
    Float_t GetNotIDConstCorrectionITSNoIDLowBound() const {return fNotIDConstITSNoIDLow;}
    Float_t GetNotIDConstCorrectionTPCHighBound() const {return fNotIDConstTPCHigh;}
    Float_t GetNotIDConstCorrectionITSHighBound() const {return fNotIDConstITSHigh;}
    Float_t GetNotIDConstCorrectionTPCNoIDHighBound() const {return fNotIDConstTPCNoIDHigh;}
    Float_t GetNotIDConstCorrectionITSNoIDHighBound() const {return fNotIDConstITSNoIDHigh;}
    TH1D *GetNotIDCorrectionTPC() const {return fnotIDTPC;}
    TH1D *GetNotIDCorrectionITS() const {return fnotIDITS;}
    TH1D *GetNotIDCorrectionNoPID() const {return fnotIDNoID;}
    TH1D *GetEfficiencyPionTPC(const int cb = -1);// const {return fEfficiencyPionTPC;}
    TH1D *GetEfficiencyKaonTPC(const int cb = -1);// const {return fEfficiencyKaonTPC;}
    TH1D *GetEfficiencyProtonTPC(const int cb = -1);// const {return fEfficiencyProtonTPC;}
    TH1D *GetEfficiencyHadronTPC(const int cb = -1);// const {return fEfficiencyHadronTPC;}
    TH1D *GetEfficiencyPionITS(const int cb = -1);// const {return fEfficiencyPionITS;}
    TH1D *GetEfficiencyKaonITS(const int cb = -1);// const {return fEfficiencyKaonITS;}
    TH1D *GetEfficiencyProtonITS(const int cb = -1);// const {return fEfficiencyProtonITS;}
    TH1D *GetEfficiencyHadronITS(const int cb = -1);// const {return fEfficiencyHadronITS;}
    TH1D *GetBackgroundCorrectionTPC() const {return fBackgroundTPC;}
    TH1D *GetBackgroundCorrectionITS() const {return fBackgroundITS;}

    //This is stored as the inverse of the correction
    Float_t GetNotIDCorrectionTPC(const float pT);//{return 1.0/(fnotIDTPC->GetBinContent(fnotIDTPC->FindBin(pT)));}
    Float_t GetNotIDCorrectionITS(const float pT);//{return 1.0/(fnotIDITS->GetBinContent(fnotIDITS->FindBin(pT)));}
    Float_t GetNotIDCorrectionNoPID(const float pT);//{return 1.0/(fnotIDNoID->GetBinContent(fnotIDNoID->FindBin(pT)));}
    //As is this...
    Float_t GetTPCEfficiencyCorrectionPion(const float pT, const int cb = -1);
    Float_t GetTPCEfficiencyCorrectionKaon(const float pT, const int cb = -1);
    Float_t GetTPCEfficiencyCorrectionProton(const float pT, const int cb = -1);
    Float_t GetTPCEfficiencyCorrectionHadron(const float pT, const int cb = -1);
    Float_t GetITSEfficiencyCorrectionPion(const float pT, const int cb = -1);
    Float_t GetITSEfficiencyCorrectionKaon(const float pT, const int cb = -1);
    Float_t GetITSEfficiencyCorrectionProton(const float pT, const int cb = -1);
    Float_t GetITSEfficiencyCorrectionHadron(const float pT, const int cb = -1);
    //...and these guys are too
    Float_t GetBackgroundCorrectionTPC(const float pT);//{return (1.0-fBackgroundTPC->GetBinContent(fBackgroundTPC->FindBin(pT)));}
    Float_t GetBackgroundCorrectionITS(const float pT);//{return (1.0-fBackgroundITS->GetBinContent(fBackgroundITS->FindBin(pT)));}
    Float_t GetEfficiencyErrorLowBound() const {return fEfficiencyErrorLow;}
    Float_t GetEfficiencyErrorHighBound() const {return fEfficiencyErrorHigh;}
    Float_t GetBackgroundErrorLowBound() const {return fBackgroundErrorLow;}
    Float_t GetBackgroundErrorHighBound() const {return fBackgroundErrorHigh;}


    void SetEtaCut(const Float_t val){fEtaCut=val;}
    void SetAcceptanceCorrectionFull(const Float_t val){fAcceptanceCorrectionFull=val;}
    void SetAcceptanceCorrectionEMCAL(const Float_t val){fAcceptanceCorrectionEMCAL=val;}
    void SetAcceptanceCorrectionPHOS(const Float_t val){fAcceptanceCorrectionPHOS=val;}
    void SetNeutralCorrection(const Float_t val){fNeutralCorrection=val;}
    void SetNotHadronicCorrection(const Float_t val){fNotHadronicCorrection=val;}
    void SetpTCutCorrectionTPC(const Float_t val){fpTcutCorrectionTPC=val;}
    void SetpTCutCorrectionITS(const Float_t val){fpTcutCorrectionITS=val;}
    void SetNotIDConstCorrectionTPC(const Float_t val){fNotIDConstTPC=val;}
    void SetNotIDConstCorrectionITS(const Float_t val){fNotIDConstITS=val;}
    void SetNotIDConstCorrectionTPCNoID(const Float_t val){fNotIDConstTPCNoID=val;}
    void SetNotIDConstCorrectionITSNoID(const Float_t val){fNotIDConstITSNoID=val;}
    void SetNeutralCorrectionLowBound(const Float_t val){fNeutralCorrectionLow=val;}
    void SetNotHadronicCorrectionLowBound(const Float_t val){fNotHadronicCorrectionLow=val;}
    void SetpTCutCorrectionTPCLowBound(const Float_t val){ffpTcutCorrectionTPCLow=val;}
    void SetpTCutCorrectionITSLowBound(const Float_t val){ffpTcutCorrectionITSLow=val;}
    void SetNeutralCorrectionHighBound(const Float_t val){fNeutralCorrectionHigh=val;}
    void SetNotHadronicCorrectionHighBound(const Float_t val){fNotHadronicCorrectionHigh=val;}
    void SetpTCutCorrectionTPCHighBound(const Float_t val){ffpTcutCorrectionTPCHigh=val;}
    void SetpTCutCorrectionITSHighBound(const Float_t val){ffpTcutCorrectionITSHigh=val;}
    void SetNotIDConstCorrectionTPCLowBound(const Float_t val){fNotIDConstTPCLow=val;}
    void SetNotIDConstCorrectionITSLowBound(const Float_t val){fNotIDConstITSLow=val;}
    void SetNotIDConstCorrectionTPCNoIDLowBound(const Float_t val){fNotIDConstTPCNoIDLow=val;}
    void SetNotIDConstCorrectionITSNoIDLowBound(const Float_t val){fNotIDConstITSNoIDLow=val;}
    void SetNotIDConstCorrectionTPCHighBound(const Float_t val){fNotIDConstTPCHigh=val;}
    void SetNotIDConstCorrectionITSHighBound(const Float_t val){fNotIDConstITSHigh=val;}
    void SetNotIDConstCorrectionTPCNoIDHighBound(const Float_t val){fNotIDConstTPCNoIDHigh=val;}
    void SetNotIDConstCorrectionITSNoIDHighBound(const Float_t val){fNotIDConstITSNoIDHigh=val;}
    void SetNotIDCorrectionTPC(const TH1D *histo){fnotIDTPC=(TH1D*) histo;}
    void SetNotIDCorrectionITS(const TH1D *histo){fnotIDITS=(TH1D*) histo;}
    void SetNotIDCorrectionNoPID(const TH1D *histo){fnotIDNoID=(TH1D*) histo;}
    void SetEfficiencyPionTPC(const TH1D *histo){fEfficiencyPionTPC=(TH1D*) histo;}
    void SetEfficiencyKaonTPC(const TH1D *histo){fEfficiencyKaonTPC=(TH1D*) histo;}
    void SetEfficiencyProtonTPC(const TH1D *histo){fEfficiencyProtonTPC=(TH1D*) histo;}
    void SetEfficiencyHadronTPC(const TH1D *histo){fEfficiencyHadronTPC=(TH1D*) histo;}
    void SetEfficiencyPionITS(const TH1D *histo){fEfficiencyPionITS=(TH1D*) histo;}
    void SetEfficiencyKaonITS(const TH1D *histo){fEfficiencyKaonITS=(TH1D*) histo;}
    void SetEfficiencyProtonITS(const TH1D *histo){fEfficiencyProtonITS=(TH1D*) histo;}
    void SetEfficiencyHadronITS(const TH1D *histo){fEfficiencyHadronITS=(TH1D*) histo;}
    void SetEfficiencyPionTPC(TH1D *histo, const int cb);
    void SetEfficiencyKaonTPC(TH1D *histo, const int cb);
    void SetEfficiencyProtonTPC(TH1D *histo, const int cb);
    void SetEfficiencyHadronTPC(TH1D *histo, const int cb);
    void SetEfficiencyPionITS(TH1D *histo, const int cb);
    void SetEfficiencyKaonITS(TH1D *histo, const int cb);
    void SetEfficiencyProtonITS(TH1D *histo, const int cb);
    void SetEfficiencyHadronITS(TH1D *histo, const int cb);
    void SetBackgroundCorrectionTPC(const TH1D *histo){fBackgroundTPC=(TH1D*) histo;}
    void SetBackgroundCorrectionITS(const TH1D *histo){fBackgroundITS=(TH1D*) histo;}
    void SetEfficiencyErrorLowBound(const Float_t val){fEfficiencyErrorLow=val;}
    void SetEfficiencyErrorHighBound(const Float_t val){fEfficiencyErrorHigh=val;}
    void SetBackgroundErrorLowBound(const Float_t val){fBackgroundErrorLow=val;}
    void SetBackgroundErrorHighBound(const Float_t val){fBackgroundErrorHigh=val;}
    void IsEMCal(Bool_t val){fIsEMCal=val;}
    void IsData(Bool_t val){fIsData=val;}
    void SetDataSet(Int_t val){fDataSet=val;}
    void SetProduction(char *prod){fProduction = prod;}
    void SetProductionDescription(char *prod){fProductionDescription = prod;}
    void Report();//Gives a report on the status of all of the corrections

    //Returns the factor one needs to multiply by to get the corrected et for all constant (not pt dependent) factors
    Float_t GetConstantCorrections(Bool_t totEt, Float_t ptcut, TString type) const;
    Float_t GetSystematicErrorBound(Float_t et,Bool_t isLowBound, Bool_t isHadronic, Bool_t isTPC) const;


    AliAnalysisHadEtCorrections(const AliAnalysisHadEtCorrections *g) ; // cpy ctor
    //AliAnalysisHadEtCorrections & operator = (const AliAnalysisHadEtCorrections & g) ;//cpy assignment

protected:

    Float_t fEtaCut;//the eta cut used for this analysis
    Float_t fAcceptanceCorrectionFull;//the acceptance correction for full azimuthal acceptance
    Float_t fAcceptanceCorrectionEMCAL;//the acceptance correction for the EMCal azimuthal acceptance
    Float_t fAcceptanceCorrectionPHOS;//the acceptance correction for the PHOS azimuthal acceptance
    //Systematic errors - low value, mean value, high value
    Float_t fNeutralCorrection;//the correction for the fraction of energy from neutral particles (for using both the calorimeters and the tracking detectors)
    Float_t fNotHadronicCorrection;//the correction for the fraction of energy which is not measured by the tracking detectors 
    Float_t fpTcutCorrectionTPC;//the correction for the momentum cut for the tpc (150 MeV/c)
    Float_t fpTcutCorrectionITS;//the correction for the momentum cut for the ITS (100 MeV/c)
    Float_t fNotIDConstTPC;//the correction for the constant correction for unidentified particles with the TPC momentum cut
    Float_t fNotIDConstITS;//the correction for the constant correction for unidentified particles with the ITS momentum cut
    Float_t fNotIDConstTPCNoID;//the correction for the constant correction for unidentified particles with the TPC momentum cut if no PID was done
    Float_t fNotIDConstITSNoID;//the correction for the constant correction for unidentified particles with the ITS momentum cut if no PID was done
    Float_t fNeutralCorrectionLow;//the low bound on the neutral energy fraction correction
    Float_t fNotHadronicCorrectionLow;//the low bound on the hadronic energy fraction correction
    Float_t ffpTcutCorrectionTPCLow;//the low bound on the TPC momentum cut correction
    Float_t ffpTcutCorrectionITSLow;//the low bound on the ITS momentum cut correction
    Float_t fNeutralCorrectionHigh;//the high bound on the neutral energy correcton
    Float_t fNotHadronicCorrectionHigh;//the high bound on the hadronic energy correction
    Float_t ffpTcutCorrectionTPCHigh;//the high bound on the TPC momentum cut correction
    Float_t ffpTcutCorrectionITSHigh;//the high bound on the ITS momentum cut correction
    Float_t fNotIDConstTPCLow;//the correction for the constant correction for unidentified particles with the TPC momentum cut
    Float_t fNotIDConstITSLow;//the correction for the constant correction for unidentified particles with the ITS momentum cut
    Float_t fNotIDConstTPCNoIDLow;//the correction for the constant correction for unidentified particles with the TPC momentum cut if no PID was done
    Float_t fNotIDConstITSNoIDLow;//the correction for the constant correction for unidentified particles with the ITS momentum cut if no PID was done
    Float_t fNotIDConstTPCHigh;//the correction for the constant correction for unidentified particles with the TPC momentum cut
    Float_t fNotIDConstITSHigh;//the correction for the constant correction for unidentified particles with the ITS momentum cut
    Float_t fNotIDConstTPCNoIDHigh;//the correction for the constant correction for unidentified particles with the TPC momentum cut if no PID was done
    Float_t fNotIDConstITSNoIDHigh;//the correction for the constant correction for unidentified particles with the ITS momentum cut if no PID was done

    //Histograms with the pT dependent fCorrections
    TH1D *fnotIDTPC;//correction for unidentified tracks in the TPC
    TH1D *fnotIDITS;//correction for unidentified tracks in the ITS
    TH1D *fnotIDNoID;//correction for unidentified tracks assuming no ID
    TH1D *fEfficiencyPionTPC;//efficiency correction for pions in the TPC
    TH1D *fEfficiencyKaonTPC;//efficiency correction for kaons in the TPC
    TH1D *fEfficiencyProtonTPC;//efficiency correction for protons in the TPC
    TH1D *fEfficiencyHadronTPC;//efficiency correction for unidentified hadrons in the TPC
    TH1D *fEfficiencyPionITS;//efficiency correction for pions in the ITS
    TH1D *fEfficiencyKaonITS;//efficiency correction for kaons in the ITS
    TH1D *fEfficiencyProtonITS;//efficiency correction for protons in the ITS
    TH1D *fEfficiencyHadronITS;//efficiency correction for unidentified hadrons in the ITS
    TObjArray *fEfficiencyTPC;//TList containing efficiencies for ITS standalone tracks for different centrality bins
    TObjArray *fEfficiencyITS;//TList containing efficiencies for ITS standalone tracks for different centrality bins
    Float_t fEfficiencyErrorLow;//Relative error on efficiency, lower bound
    Float_t fEfficiencyErrorHigh;//Relative error on efficiency, upper bound
    Float_t fBackgroundErrorLow;//Relative error on efficiency, lower bound
    Float_t fBackgroundErrorHigh;//Relative error on efficiency, upper bound
    TH1D *fBackgroundTPC;//background correction for the TPC
    TH1D *fBackgroundITS;//background correction for the ITS
    Bool_t fIsEMCal;//boolean to keep track of whether this is for EMCal or PHOS acceptance
    Bool_t fIsData;//boolean to keep track of whether this is for data or simulation acceptance
    Int_t fDataSet;//integer to keep track of data set
    TString fProduction;//Short version name of production
    TString fProductionDescription;//Long description of production


private:
  //Declare it private to avoid compilation warning
    AliAnalysisHadEtCorrections & operator = (const AliAnalysisHadEtCorrections & g) ;//cpy assignment
    AliAnalysisHadEtCorrections(const AliAnalysisHadEtCorrections & g) ; // cpy ctor

    ClassDef(AliAnalysisHadEtCorrections, 1);
};

#endif // ALIANALYSISHADETCORRECTIONS_H
Commit	Line	Data
cf6522d1	1	//Created by Christine Nattrass, Rebecca Scott, Irakli Martashvili
020abff0	2	//University of Tennessee at Knoxville
cf6522d1	3	// This is a container class for the correction factors for the hadronic
	4	// component of transverse energy
	5	// It is filled by the output of AliAnalysisTaskHadEt from spinning over Monte
	6	// Carlo data (using AliAnalysisHadEtMonteCarlo)
	7	// It is used by AliAnalysisTaskHadEt while spinning over reconstructed data
	8	// (using AliAnalysisHadEtReconstructed)
020abff0	9	//Please see https://twiki.cern.ch/twiki/bin/view/ALICE/ETCaloAnalysis
	10	#ifndef ALIANALYSISHADETCORRECTIONS_H
	11	#define ALIANALYSISHADETCORRECTIONS_H
	12
020abff0	13	#include "TString.h"
45e5a6ec	14	#include "TNamed.h"
cf6522d1	15	class Rtypes;
45e5a6ec	16	class TH1D;
cf6522d1	17	class TNamed;
226c54bf	18	class TObjArray;
020abff0	19
	20	class AliAnalysisHadEtCorrections : public TNamed
	21	{
	22	public:
	23
	24	AliAnalysisHadEtCorrections();
	25	virtual ~AliAnalysisHadEtCorrections();
45e5a6ec	26	virtual void Init();
45e5a6ec	27
020abff0	28
	29	Float_t GetEtaCut() const {return fEtaCut;}
	30	Float_t GetAcceptanceCorrectionFull() const {return fAcceptanceCorrectionFull;}
	31	Float_t GetAcceptanceCorrectionEMCAL() const {return fAcceptanceCorrectionEMCAL;}
	32	Float_t GetAcceptanceCorrectionPHOS() const {return fAcceptanceCorrectionPHOS;}
	33	Float_t GetNeutralCorrection() const {return fNeutralCorrection;}
	34	Float_t GetNotHadronicCorrection() const {return fNotHadronicCorrection;}
	35	Float_t GetpTCutCorrectionTPC() const {return fpTcutCorrectionTPC;}
	36	Float_t GetpTCutCorrectionITS() const {return fpTcutCorrectionITS;}
0e866ddc	37	Float_t GetNotIDConstCorrectionTPC() const {return fNotIDConstTPC;}
	38	Float_t GetNotIDConstCorrectionITS() const {return fNotIDConstITS;}
	39	Float_t GetNotIDConstCorrectionTPCNoID() const {return fNotIDConstTPCNoID;}
	40	Float_t GetNotIDConstCorrectionITSNoID() const {return fNotIDConstITSNoID;}
020abff0	41	Float_t GetNeutralCorrectionLowBound() const {return fNeutralCorrectionLow;}
	42	Float_t GetNotHadronicCorrectionLowBound() const {return fNotHadronicCorrectionLow;}
	43	Float_t GetpTCutCorrectionTPCLowBound() const {return ffpTcutCorrectionTPCLow;}
	44	Float_t GetpTCutCorrectionITSLowBound() const {return ffpTcutCorrectionITSLow;}
	45	Float_t GetNeutralCorrectionHighBound() const {return fNeutralCorrectionHigh;}
	46	Float_t GetNotHadronicCorrectionHighBound() const {return fNotHadronicCorrectionHigh;}
	47	Float_t GetpTCutCorrectionTPCHighBound() const {return ffpTcutCorrectionTPCHigh;}
	48	Float_t GetpTCutCorrectionITSHighBound() const {return ffpTcutCorrectionITSHigh;}
0e866ddc	49	Float_t GetNotIDConstCorrectionTPCLowBound() const {return fNotIDConstTPCLow;}
	50	Float_t GetNotIDConstCorrectionITSLowBound() const {return fNotIDConstITSLow;}
	51	Float_t GetNotIDConstCorrectionTPCNoIDLowBound() const {return fNotIDConstTPCNoIDLow;}
	52	Float_t GetNotIDConstCorrectionITSNoIDLowBound() const {return fNotIDConstITSNoIDLow;}
825704b5	53	Float_t GetNotIDConstCorrectionTPCHighBound() const {return fNotIDConstTPCHigh;}
	54	Float_t GetNotIDConstCorrectionITSHighBound() const {return fNotIDConstITSHigh;}
	55	Float_t GetNotIDConstCorrectionTPCNoIDHighBound() const {return fNotIDConstTPCNoIDHigh;}
	56	Float_t GetNotIDConstCorrectionITSNoIDHighBound() const {return fNotIDConstITSNoIDHigh;}
020abff0	57	TH1D *GetNotIDCorrectionTPC() const {return fnotIDTPC;}
	58	TH1D *GetNotIDCorrectionITS() const {return fnotIDITS;}
	59	TH1D *GetNotIDCorrectionNoPID() const {return fnotIDNoID;}
45f9be5c	60	TH1D *GetEfficiencyPionTPC(const int cb = -1);// const {return fEfficiencyPionTPC;}
	61	TH1D *GetEfficiencyKaonTPC(const int cb = -1);// const {return fEfficiencyKaonTPC;}
	62	TH1D *GetEfficiencyProtonTPC(const int cb = -1);// const {return fEfficiencyProtonTPC;}
	63	TH1D *GetEfficiencyHadronTPC(const int cb = -1);// const {return fEfficiencyHadronTPC;}
	64	TH1D *GetEfficiencyPionITS(const int cb = -1);// const {return fEfficiencyPionITS;}
	65	TH1D *GetEfficiencyKaonITS(const int cb = -1);// const {return fEfficiencyKaonITS;}
	66	TH1D *GetEfficiencyProtonITS(const int cb = -1);// const {return fEfficiencyProtonITS;}
	67	TH1D *GetEfficiencyHadronITS(const int cb = -1);// const {return fEfficiencyHadronITS;}
020abff0	68	TH1D *GetBackgroundCorrectionTPC() const {return fBackgroundTPC;}
	69	TH1D *GetBackgroundCorrectionITS() const {return fBackgroundITS;}
	70
	71	//This is stored as the inverse of the correction
0e866ddc	72	Float_t GetNotIDCorrectionTPC(const float pT);//{return 1.0/(fnotIDTPC->GetBinContent(fnotIDTPC->FindBin(pT)));}
	73	Float_t GetNotIDCorrectionITS(const float pT);//{return 1.0/(fnotIDITS->GetBinContent(fnotIDITS->FindBin(pT)));}
	74	Float_t GetNotIDCorrectionNoPID(const float pT);//{return 1.0/(fnotIDNoID->GetBinContent(fnotIDNoID->FindBin(pT)));}
020abff0	75	//As is this...
226c54bf	76	Float_t GetTPCEfficiencyCorrectionPion(const float pT, const int cb = -1);
	77	Float_t GetTPCEfficiencyCorrectionKaon(const float pT, const int cb = -1);
	78	Float_t GetTPCEfficiencyCorrectionProton(const float pT, const int cb = -1);
	79	Float_t GetTPCEfficiencyCorrectionHadron(const float pT, const int cb = -1);
	80	Float_t GetITSEfficiencyCorrectionPion(const float pT, const int cb = -1);
	81	Float_t GetITSEfficiencyCorrectionKaon(const float pT, const int cb = -1);
	82	Float_t GetITSEfficiencyCorrectionProton(const float pT, const int cb = -1);
	83	Float_t GetITSEfficiencyCorrectionHadron(const float pT, const int cb = -1);
020abff0	84	//...and these guys are too
45e5a6ec	85	Float_t GetBackgroundCorrectionTPC(const float pT);//{return (1.0-fBackgroundTPC->GetBinContent(fBackgroundTPC->FindBin(pT)));}
45e5a6ec	86	Float_t GetBackgroundCorrectionITS(const float pT);//{return (1.0-fBackgroundITS->GetBinContent(fBackgroundITS->FindBin(pT)));}
57ea58bf	87	Float_t GetEfficiencyErrorLowBound() const {return fEfficiencyErrorLow;}
	88	Float_t GetEfficiencyErrorHighBound() const {return fEfficiencyErrorHigh;}
	89	Float_t GetBackgroundErrorLowBound() const {return fBackgroundErrorLow;}
	90	Float_t GetBackgroundErrorHighBound() const {return fBackgroundErrorHigh;}
020abff0	91
020abff0	92
cf6522d1	93	void SetEtaCut(const Float_t val){fEtaCut=val;}
	94	void SetAcceptanceCorrectionFull(const Float_t val){fAcceptanceCorrectionFull=val;}
	95	void SetAcceptanceCorrectionEMCAL(const Float_t val){fAcceptanceCorrectionEMCAL=val;}
	96	void SetAcceptanceCorrectionPHOS(const Float_t val){fAcceptanceCorrectionPHOS=val;}
	97	void SetNeutralCorrection(const Float_t val){fNeutralCorrection=val;}
	98	void SetNotHadronicCorrection(const Float_t val){fNotHadronicCorrection=val;}
	99	void SetpTCutCorrectionTPC(const Float_t val){fpTcutCorrectionTPC=val;}
	100	void SetpTCutCorrectionITS(const Float_t val){fpTcutCorrectionITS=val;}
0e866ddc	101	void SetNotIDConstCorrectionTPC(const Float_t val){fNotIDConstTPC=val;}
	102	void SetNotIDConstCorrectionITS(const Float_t val){fNotIDConstITS=val;}
	103	void SetNotIDConstCorrectionTPCNoID(const Float_t val){fNotIDConstTPCNoID=val;}
	104	void SetNotIDConstCorrectionITSNoID(const Float_t val){fNotIDConstITSNoID=val;}
cf6522d1	105	void SetNeutralCorrectionLowBound(const Float_t val){fNeutralCorrectionLow=val;}
	106	void SetNotHadronicCorrectionLowBound(const Float_t val){fNotHadronicCorrectionLow=val;}
	107	void SetpTCutCorrectionTPCLowBound(const Float_t val){ffpTcutCorrectionTPCLow=val;}
	108	void SetpTCutCorrectionITSLowBound(const Float_t val){ffpTcutCorrectionITSLow=val;}
	109	void SetNeutralCorrectionHighBound(const Float_t val){fNeutralCorrectionHigh=val;}
	110	void SetNotHadronicCorrectionHighBound(const Float_t val){fNotHadronicCorrectionHigh=val;}
	111	void SetpTCutCorrectionTPCHighBound(const Float_t val){ffpTcutCorrectionTPCHigh=val;}
	112	void SetpTCutCorrectionITSHighBound(const Float_t val){ffpTcutCorrectionITSHigh=val;}
0e866ddc	113	void SetNotIDConstCorrectionTPCLowBound(const Float_t val){fNotIDConstTPCLow=val;}
	114	void SetNotIDConstCorrectionITSLowBound(const Float_t val){fNotIDConstITSLow=val;}
	115	void SetNotIDConstCorrectionTPCNoIDLowBound(const Float_t val){fNotIDConstTPCNoIDLow=val;}
	116	void SetNotIDConstCorrectionITSNoIDLowBound(const Float_t val){fNotIDConstITSNoIDLow=val;}
	117	void SetNotIDConstCorrectionTPCHighBound(const Float_t val){fNotIDConstTPCHigh=val;}
	118	void SetNotIDConstCorrectionITSHighBound(const Float_t val){fNotIDConstITSHigh=val;}
	119	void SetNotIDConstCorrectionTPCNoIDHighBound(const Float_t val){fNotIDConstTPCNoIDHigh=val;}
	120	void SetNotIDConstCorrectionITSNoIDHighBound(const Float_t val){fNotIDConstITSNoIDHigh=val;}
cf6522d1	121	void SetNotIDCorrectionTPC(const TH1D histo){fnotIDTPC=(TH1D) histo;}
	122	void SetNotIDCorrectionITS(const TH1D histo){fnotIDITS=(TH1D) histo;}
	123	void SetNotIDCorrectionNoPID(const TH1D histo){fnotIDNoID=(TH1D) histo;}
	124	void SetEfficiencyPionTPC(const TH1D histo){fEfficiencyPionTPC=(TH1D) histo;}
	125	void SetEfficiencyKaonTPC(const TH1D histo){fEfficiencyKaonTPC=(TH1D) histo;}
	126	void SetEfficiencyProtonTPC(const TH1D histo){fEfficiencyProtonTPC=(TH1D) histo;}
	127	void SetEfficiencyHadronTPC(const TH1D histo){fEfficiencyHadronTPC=(TH1D) histo;}
	128	void SetEfficiencyPionITS(const TH1D histo){fEfficiencyPionITS=(TH1D) histo;}
	129	void SetEfficiencyKaonITS(const TH1D histo){fEfficiencyKaonITS=(TH1D) histo;}
	130	void SetEfficiencyProtonITS(const TH1D histo){fEfficiencyProtonITS=(TH1D) histo;}
	131	void SetEfficiencyHadronITS(const TH1D histo){fEfficiencyHadronITS=(TH1D) histo;}
226c54bf	132	void SetEfficiencyPionTPC(TH1D *histo, const int cb);
	133	void SetEfficiencyKaonTPC(TH1D *histo, const int cb);
	134	void SetEfficiencyProtonTPC(TH1D *histo, const int cb);
	135	void SetEfficiencyHadronTPC(TH1D *histo, const int cb);
	136	void SetEfficiencyPionITS(TH1D *histo, const int cb);
	137	void SetEfficiencyKaonITS(TH1D *histo, const int cb);
	138	void SetEfficiencyProtonITS(TH1D *histo, const int cb);
	139	void SetEfficiencyHadronITS(TH1D *histo, const int cb);
cf6522d1	140	void SetBackgroundCorrectionTPC(const TH1D histo){fBackgroundTPC=(TH1D) histo;}
cf6522d1	141	void SetBackgroundCorrectionITS(const TH1D histo){fBackgroundITS=(TH1D) histo;}
45f9be5c	142	void SetEfficiencyErrorLowBound(const Float_t val){fEfficiencyErrorLow=val;}
	143	void SetEfficiencyErrorHighBound(const Float_t val){fEfficiencyErrorHigh=val;}
	144	void SetBackgroundErrorLowBound(const Float_t val){fBackgroundErrorLow=val;}
	145	void SetBackgroundErrorHighBound(const Float_t val){fBackgroundErrorHigh=val;}
226c54bf	146	void IsEMCal(Bool_t val){fIsEMCal=val;}
	147	void IsData(Bool_t val){fIsData=val;}
	148	void SetDataSet(Int_t val){fDataSet=val;}
	149	void SetProduction(char *prod){fProduction = prod;}
	150	void SetProductionDescription(char *prod){fProductionDescription = prod;}
	151	void Report();//Gives a report on the status of all of the corrections
020abff0	152
3ce6b879	153	//Returns the factor one needs to multiply by to get the corrected et for all constant (not pt dependent) factors
464aa50c	154	Float_t GetConstantCorrections(Bool_t totEt, Float_t ptcut, TString type) const;
d263b18f	155	Float_t GetSystematicErrorBound(Float_t et,Bool_t isLowBound, Bool_t isHadronic, Bool_t isTPC) const;
3ce6b879	156
020abff0	157
	158	AliAnalysisHadEtCorrections(const AliAnalysisHadEtCorrections *g) ; // cpy ctor
	159	//AliAnalysisHadEtCorrections & operator = (const AliAnalysisHadEtCorrections & g) ;//cpy assignment
	160
	161	protected:
	162
	163	Float_t fEtaCut;//the eta cut used for this analysis
	164	Float_t fAcceptanceCorrectionFull;//the acceptance correction for full azimuthal acceptance
	165	Float_t fAcceptanceCorrectionEMCAL;//the acceptance correction for the EMCal azimuthal acceptance
	166	Float_t fAcceptanceCorrectionPHOS;//the acceptance correction for the PHOS azimuthal acceptance
	167	//Systematic errors - low value, mean value, high value
	168	Float_t fNeutralCorrection;//the correction for the fraction of energy from neutral particles (for using both the calorimeters and the tracking detectors)
	169	Float_t fNotHadronicCorrection;//the correction for the fraction of energy which is not measured by the tracking detectors
	170	Float_t fpTcutCorrectionTPC;//the correction for the momentum cut for the tpc (150 MeV/c)
	171	Float_t fpTcutCorrectionITS;//the correction for the momentum cut for the ITS (100 MeV/c)
0e866ddc	172	Float_t fNotIDConstTPC;//the correction for the constant correction for unidentified particles with the TPC momentum cut
	173	Float_t fNotIDConstITS;//the correction for the constant correction for unidentified particles with the ITS momentum cut
	174	Float_t fNotIDConstTPCNoID;//the correction for the constant correction for unidentified particles with the TPC momentum cut if no PID was done
	175	Float_t fNotIDConstITSNoID;//the correction for the constant correction for unidentified particles with the ITS momentum cut if no PID was done
020abff0	176	Float_t fNeutralCorrectionLow;//the low bound on the neutral energy fraction correction
	177	Float_t fNotHadronicCorrectionLow;//the low bound on the hadronic energy fraction correction
	178	Float_t ffpTcutCorrectionTPCLow;//the low bound on the TPC momentum cut correction
	179	Float_t ffpTcutCorrectionITSLow;//the low bound on the ITS momentum cut correction
	180	Float_t fNeutralCorrectionHigh;//the high bound on the neutral energy correcton
	181	Float_t fNotHadronicCorrectionHigh;//the high bound on the hadronic energy correction
	182	Float_t ffpTcutCorrectionTPCHigh;//the high bound on the TPC momentum cut correction
	183	Float_t ffpTcutCorrectionITSHigh;//the high bound on the ITS momentum cut correction
0e866ddc	184	Float_t fNotIDConstTPCLow;//the correction for the constant correction for unidentified particles with the TPC momentum cut
	185	Float_t fNotIDConstITSLow;//the correction for the constant correction for unidentified particles with the ITS momentum cut
	186	Float_t fNotIDConstTPCNoIDLow;//the correction for the constant correction for unidentified particles with the TPC momentum cut if no PID was done
	187	Float_t fNotIDConstITSNoIDLow;//the correction for the constant correction for unidentified particles with the ITS momentum cut if no PID was done
	188	Float_t fNotIDConstTPCHigh;//the correction for the constant correction for unidentified particles with the TPC momentum cut
	189	Float_t fNotIDConstITSHigh;//the correction for the constant correction for unidentified particles with the ITS momentum cut
	190	Float_t fNotIDConstTPCNoIDHigh;//the correction for the constant correction for unidentified particles with the TPC momentum cut if no PID was done
	191	Float_t fNotIDConstITSNoIDHigh;//the correction for the constant correction for unidentified particles with the ITS momentum cut if no PID was done
020abff0	192
464aa50c	193	//Histograms with the pT dependent fCorrections
020abff0	194	TH1D *fnotIDTPC;//correction for unidentified tracks in the TPC
	195	TH1D *fnotIDITS;//correction for unidentified tracks in the ITS
	196	TH1D *fnotIDNoID;//correction for unidentified tracks assuming no ID
	197	TH1D *fEfficiencyPionTPC;//efficiency correction for pions in the TPC
	198	TH1D *fEfficiencyKaonTPC;//efficiency correction for kaons in the TPC
	199	TH1D *fEfficiencyProtonTPC;//efficiency correction for protons in the TPC
	200	TH1D *fEfficiencyHadronTPC;//efficiency correction for unidentified hadrons in the TPC
	201	TH1D *fEfficiencyPionITS;//efficiency correction for pions in the ITS
	202	TH1D *fEfficiencyKaonITS;//efficiency correction for kaons in the ITS
	203	TH1D *fEfficiencyProtonITS;//efficiency correction for protons in the ITS
	204	TH1D *fEfficiencyHadronITS;//efficiency correction for unidentified hadrons in the ITS
226c54bf	205	TObjArray *fEfficiencyTPC;//TList containing efficiencies for ITS standalone tracks for different centrality bins
226c54bf	206	TObjArray *fEfficiencyITS;//TList containing efficiencies for ITS standalone tracks for different centrality bins
45f9be5c	207	Float_t fEfficiencyErrorLow;//Relative error on efficiency, lower bound
	208	Float_t fEfficiencyErrorHigh;//Relative error on efficiency, upper bound
	209	Float_t fBackgroundErrorLow;//Relative error on efficiency, lower bound
	210	Float_t fBackgroundErrorHigh;//Relative error on efficiency, upper bound
020abff0	211	TH1D *fBackgroundTPC;//background correction for the TPC
020abff0	212	TH1D *fBackgroundITS;//background correction for the ITS
226c54bf	213	Bool_t fIsEMCal;//boolean to keep track of whether this is for EMCal or PHOS acceptance
	214	Bool_t fIsData;//boolean to keep track of whether this is for data or simulation acceptance
	215	Int_t fDataSet;//integer to keep track of data set
	216	TString fProduction;//Short version name of production
	217	TString fProductionDescription;//Long description of production
020abff0	218
	219
	220
	221	private:
	222	//Declare it private to avoid compilation warning
	223	AliAnalysisHadEtCorrections & operator = (const AliAnalysisHadEtCorrections & g) ;//cpy assignment
	224	AliAnalysisHadEtCorrections(const AliAnalysisHadEtCorrections & g) ; // cpy ctor
	225
	226	ClassDef(AliAnalysisHadEtCorrections, 1);
	227	};
	228
	229	#endif // ALIANALYSISHADETCORRECTIONS_H