X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PWGHF%2Fhfe%2FAliAnalysisTaskHFECal.h;h=86bf16e4a7f85365815755902f1329fbb4e896e3;hb=307e1ff7ab30795aa778e42133ddafbc306485f4;hp=6f51813774ffef037c75cde35764635b42876bf9;hpb=d19af75d6c56e34087bb84b708064dcae13e5364;p=u%2Fmrichter%2FAliRoot.git diff --git a/PWGHF/hfe/AliAnalysisTaskHFECal.h b/PWGHF/hfe/AliAnalysisTaskHFECal.h index 6f51813774f..86bf16e4a7f 100644 --- a/PWGHF/hfe/AliAnalysisTaskHFECal.h +++ b/PWGHF/hfe/AliAnalysisTaskHFECal.h @@ -51,11 +51,19 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE { void SetHFECuts(AliHFEcuts * const cuts) { fCuts = cuts; }; void SetOpeningAngleCut (Double_t openingAngle) {fOpeningAngleCut = openingAngle;}; + void SetMimpTassCut (Double_t MimpTassCut) {fMimpTassCut = MimpTassCut;}; + void SetMimNsigassCut (Double_t MimNsigassCut) {fMimNsigassCut = MimNsigassCut;}; void SetInvariantMassCut (Double_t invmass) {fInvmassCut = invmass;}; + void SetMassConstraint (Bool_t MassConstraint) { fSetMassConstraint = MassConstraint; }; + void SetMassWidthCut (Bool_t MassWidthCut) { fSetMassWidthCut = MassWidthCut; }; + void SetMassNonlinear (Bool_t MassNonlinear) { fSetMassNonlinear = MassNonlinear; }; + void SetMassCalMethod (Bool_t KFpart) { fSetKFpart = KFpart; }; void SetQAHist (int qahist) {fqahist = qahist;}; AliHFEpid *GetPID() const { return fPID; } void SetRejectKinkMother(Bool_t rejectKinkMother = kFALSE) { fRejectKinkMother = rejectKinkMother; }; void SelectPhotonicElectron(Int_t itrack, Double_t cent, AliESDtrack *track, Bool_t &fFlagPhotonicElec, Bool_t &fFlagConvinatElec, Double_t nSig, Double_t shower, Double_t ep, Double_t mce, Double_t w, Int_t ibgevent, Bool_t tagpi0, Bool_t tageta); + void SelectPhotonicElectron2(Int_t itrack, Double_t cent, AliESDtrack *track, Bool_t &fFlagPhotonicElec, Bool_t &fFlagConvinatElec, Double_t nSig, Double_t shower, Double_t ep, Double_t mce, Double_t w, Int_t ibgevent, Bool_t tagpi0, Bool_t tageta); + double SumpT(Int_t itrack, AliESDtrack *track); void FindMother(TParticle* part, int &label, int &pid); double GetMCweight(double mcPi0pT); double GetMCweightEta(double mcEtapT); @@ -85,7 +93,13 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE { AliHFEpid *fPID; //! PID AliHFEpidQAmanager *fPIDqa; //! PID QA manager Double_t fOpeningAngleCut; //openingAngle cut value + Double_t fMimpTassCut; //openingAngle cut value + Double_t fMimNsigassCut; //openingAngle cut value Double_t fInvmassCut; //invariant mass cut value + Bool_t fSetMassConstraint; // Set mass constraint + Bool_t fSetMassWidthCut; // Set mass constraint + Bool_t fSetMassNonlinear; // Set mass constraint + Bool_t fSetKFpart; // Set mass constraint int ftriggers[48][60];//! int ftriggersCut[48][60];//! @@ -106,6 +120,8 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE { THnSparseD *fInvmassULS; //! Inv mass of ULS (e,e) THnSparseD *fInvmassLSmc; //! Inv mass of LS (e,e) THnSparseD *fInvmassULSmc; //! Inv mass of ULS (e,e) + TH2D *fInvmassLSreco; //! Inv mass of LS (e,e) + TH2D *fInvmassULSreco; //! Inv mass of ULS (e,e) TH2D *fInvmassLSmc0; //! Inv mass of ULS (e,e) TH2D *fInvmassLSmc1; //! Inv mass of ULS (e,e) TH2D *fInvmassLSmc2; //! Inv mass of ULS (e,e) @@ -128,6 +144,7 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE { TH1F *fCent; //! centrality THnSparseD *fEleInfo; //! EMC acc + THnSparseD *fElenSigma; //! EMC acc /* //<---- trigger info TH1F *fClsEBftTrigCut; //Cluster E before trigger selection @@ -166,6 +183,9 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE { THnSparseD *fIncpTMCpho_pi0e_TPC; //! MC HFE pid electron vs centrality THnSparseD *fPhoElecPtMC_pi0e_TPC; //! Pho inclusive ele pt THnSparseD *fSameElecPtMC_pi0e_TPC; //! Same inclusive ele pt + THnSparseD *fIncpTMCpho_eta_TPC; //! MC HFE pid electron vs centrality + THnSparseD *fPhoElecPtMC_eta_TPC; //! Pho inclusive ele pt + THnSparseD *fSameElecPtMC_eta_TPC; //! Same inclusive ele pt TH1D *CheckNclust; TH1D *CheckNits; THnSparseD *Hpi0pTcheck; @@ -180,7 +200,14 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE { TH2D *fFakeTrk0; TH2D *fFakeTrk1; TH2D *ftimingEle; - + TH2D *fIncMaxE; + TH2D *fIncReco; + TH2D *fPhoReco; + TH2D *fSamReco; + TH2D *fIncRecoMaxE; + TH2D *fPhoRecoMaxE; + TH2D *fSamRecoMaxE; + //<----- correction TGraphErrors *fnSigEtaCorr[7];