1 #ifndef ALIHADCORRTASK_H
2 #define ALIHADCORRTASK_H
10 class AliEmcalParticle;
13 #include "AliAnalysisTaskEmcal.h"
15 class AliHadCorrTask : public AliAnalysisTaskEmcal {
19 AliHadCorrTask(const char *name);
20 AliHadCorrTask(const char *name, Bool_t histo);
21 virtual ~AliHadCorrTask();
23 void UserCreateOutputObjects();
24 void Terminate(Option_t *);
26 void SetEexcl(Double_t Emin) { fEexclCell = Emin ; }
27 void SetEtaMatch(Double_t eta) { fEtaMatch = eta ; }
28 void SetHadCorr(Double_t c) { fHadCorr = c ; }
29 void SetOutClusName(const char *n) { fOutCaloName = n ; }
30 void SetPhiMatch(Double_t phi) { fPhiMatch = phi ; }
31 void SetTrackClus(Int_t c) { fDoTrackClus = c ; }
34 Double_t ApplyHadCorrOneTrack(AliEmcalParticle *emccluster, Double_t hadCorr);
35 Double_t ApplyHadCorrAllTracks(AliEmcalParticle *emccluster, Double_t hadCorr);
36 void DoMatchedTracksLoop(AliEmcalParticle *emccluster, Double_t &totalTrkP, Int_t &Nmatches);
37 void DoTrackClusLoop() ;
38 Double_t GetEtaSigma(Int_t pbin) const;
39 Int_t GetMomBin(Double_t pt) const;
40 Double_t GetPhiMean(Int_t pbin, Int_t centbin) const;
41 Double_t GetPhiSigma(Int_t pbin, Int_t centbin) const;
42 virtual Bool_t Run() ;
44 TString fOutCaloName; // name of output clusters
45 Double_t fPhiMatch; // phi match value (pp=0.050)
46 Double_t fEtaMatch; // eta match value (pp=0.025)
47 Int_t fDoTrackClus; // loop over tracks first
48 Double_t fHadCorr; // hadronic correction (fraction)
49 Double_t fEexclCell; // energy/cell that we cannot subtract from the clusters
50 TClonesArray *fOutClusters; //!output cluster collection
51 TH2F *fHistMatchEtaPhi[8][9][2]; //!output histograms
52 TH2F *fHistMatchEvsP[4]; //!output histograms
53 TH2F *fHistNMatchEnergy[4]; //!output histograms
54 TH2F *fHistNCellsEnergy[8][4]; //!output histograms
55 TH2F *fHistMatchdRvsEP[4]; //!output histograms
56 TH1F *fHistNclusvsCent; //!output histograms
57 TH1F *fHistNclusMatchvsCent; //!output histograms
58 TH1F *fHistEbefore; //!output histograms
59 TH1F *fHistEafter; //!output histograms
60 TH2F *fHistEoPCent; //!output histograms
61 TH2F *fHistNMatchCent; //!output histograms
62 TH2F *fHistNMatchCent_trk; //!output histograms
63 TH1F *fHistEsubPch[8]; //!output histograms
64 TH2F *fHistEsubPchRat[8]; //!output histograms
65 TH1F *fHistCentrality; //!output histograms
66 TH2F *fHistNoMatchEtaPhi; //!output histograms
69 AliHadCorrTask(const AliHadCorrTask&); // not implemented
70 AliHadCorrTask &operator=(const AliHadCorrTask&); // not implemented
72 ClassDef(AliHadCorrTask, 9) // Hadronic correction task