AliPhysics  aaf9c62 (aaf9c62)
AliAnalysisTaskJetPP.h
Go to the documentation of this file.
1 
10 #ifndef ALIANALYSISTASKJETPP_H
11 #define ALIANALYSISTASKJETPP_H
12 
13 
14 class TH1I;
15 class TH1F;
16 class TH2F;
17 class TH2D;
18 class TH1D;
19 class TArrayD;
20 class THnSparse;
21 class TProfile;
22 class TList;
23 class AliEmcalList;
24 class TClonesArray;
25 class TString;
26 class AliEmcalJet;
27 class AliRhoParameter;
28 class AliVParticle;
29 class AliLog;
30 class AliAnalysisUtils;
31 class TRandom3;
32 class AliJetContainer;
35 
36 #include <vector>
38 
39 using std::vector;
40 
41 
50  public:
51 
52  enum MyContainer {
57  };
58 
59 
60 
61  // ######### CONTRUCTORS/DESTRUCTORS AND STD FUNCTIONS
63  AliAnalysisTaskJetPP(const char *name);
64  virtual ~AliAnalysisTaskJetPP();
66  void Terminate(Option_t *);
67 
68 
69  // ######### SETTERS/GETTERS
70 
72  void SetUsePileUpCut (Bool_t val) {fUsePileUpCut = val;}
73  void SetSignalJetMinArea(Double_t minArea) {fMinJetArea = minArea;}
74  void SetMC(Bool_t a) {fIsMC=a;}
75 
76  void SetAcceptanceWindows(Double_t trackEta, Double_t signalJetRadius){
77  fTrackEtaWindow = trackEta;
78  fSignalJetRadius = signalJetRadius;
80  }
81 
82  void SetCentralityType(const char* type){
83  fCentralityType = type;
84  }
85 
86  void SetVertexCut(Double_t vz){ fZVertexCut = vz; }
87  void SetMinTrackPt(Double_t mpt){ fMinTrackPt = mpt;}
88 
89  void SetASPDCvsTCut(Float_t a) { fASPDCvsTCut = a; }//PILEUP
91 
93  Bool_t Run();
95 
96  private:
97 
98  Double_t EstimateBgKT(AliJetContainer *jetCont); // median p/A of kt jets
99  Double_t EstimateLocalBg(AliJetContainer *jetCont, AliParticleContainer *trkCont); // local bg in a perpendicular cone
100  Double_t EstimateBgKTCMS(AliJetContainer *jetCont); //CMS-style background form kT jets
101 
102  Double_t Convert(Double_t input);
103  // ######### CHECK FUNCTIONS
104  Bool_t IsTrackInAcceptance(AliVParticle* track, Bool_t isprimary=0);
105  Bool_t IsEventInAcceptance(AliVEvent* event);
106  Bool_t IsSPDClusterVsTrackletBG(AliVEvent *event);
107  Bool_t IsSignalJetInAcceptance(AliEmcalJet* jet, Bool_t suppressGhost=1);
108 
109 
110  // ######### STANDARD FUNCTIONS
111  void ExecOnceLocal();
112 
113  // ########## USAGE TRIGGERS
114 
118 
119  // ########## JET/DIJET/RC PROPERTIES
121  // ########## CUTS
127 
131 
132  // ########## GENERAL ////VARS
133  AliAnalysisUtils* fHelperClass;
135 
136 
142 
149 
155 
161 
166 
171 
176 
181 
185 
186 
189 
191  ClassDef(AliAnalysisTaskJetPP, 4); // Charged jet analysis for pA,,, increase the last parametre by 1 after a change
193 };
194 #endif
void ExecOnceLocal()
Perform steps needed to initialize the analysis.
AliAnalysisTaskJetPP & operator=(const AliAnalysisTaskJetPP &)
Double_t EstimateBgKT(AliJetContainer *jetCont)
Bool_t IsSignalJetInAcceptance(AliEmcalJet *jet, Bool_t suppressGhost=1)
void SetASPDCvsTCut(Float_t a)
double Double_t
Definition: External.C:58
Definition: External.C:236
TH1D * fhXVertex
! X vertex
TH2D * fhRemx
! Response matrix
Bool_t IsTrackInAcceptance(AliVParticle *track, Bool_t isprimary=0)
TH2D * fhTrackEtaPhi
! Track eta-phi distribution
TH2D * fhGenTrackEtaPt
! Generated track eta-pt distribution
Bool_t IsSPDClusterVsTrackletBG(AliVEvent *event)
void SetCentralityType(const char *type)
TH1D * fhZVertex
! Z vertex
TH1D * fhJetPtConeRho
! pt spectrum of AKT jets without local bgk
TH2D * fhInvPtQVsEta[2]
! Inverse pt vs eta
Bool_t Run()
Run the nalysis.
TH2D * fhTrackPhiCG
! Global tracks
TH2D * fhTrackPhiTPCG
! complementary tracks
Double_t fMinJetArea
Min jet area to be accepted.
Double_t fSignalJetRadius
Radius for the signal jets.
TH1D * fhRecTrkPt
! Pt spectrum of correctly reconstructed tracks
void SetBSPDCvsTCut(Float_t b)
TH1D * fhPrimGenTrkPt
! Pt spectrum of MC particle
TH2D * fhTrackEtaPt
! Track eta-pt distribution
TH1D * fhMult
! Multiplicity
TH2D * fhInvPtQVsPhiCSide[2]
! Inverse pt vs pi in close-side rapidity
TH1D * fhTrackPt
! Track pt
Double_t EstimateLocalBg(AliJetContainer *jetCont, AliParticleContainer *trkCont)
void SetMinTrackPt(Double_t mpt)
TH1D * fhJetConstituentPt
! Jet constituent pt
Double_t fSignalJetEtaWindow
+- window in eta for signal jets
Container for particles within the EMCAL framework.
TH1D * fhCMSRho
! CMS KT bgk
Pure virtual base class for AliAnalysisTaskJetPP.
TH1D * fhJetPt
! pt spectrum of AKT jets
void SetVertexCut(Double_t vz)
TH2D * fhAktJetEtaPhi
! Jet et-phi distribution
Double_t Convert(Double_t input)
TH1D * fhAtimesRhoCone
! cone jet background times area
Float_t fBSPDCvsTCut
SPD cluster-vs-tracklet cut.
Definition: External.C:204
Double_t fTrackEtaWindow
gc +- window in eta for tracks
TH2D * fhTrackPhiPt
! Track phi-pt distribution
TH2D * fhJetAreaPt
! Jet Area-pt distribution
float Float_t
Definition: External.C:68
void SetUseDefaultVertexCut(Bool_t val)
Double_t fZVertexCut
zvertex cut
TH2D * fhKtJetEtaPhi
! Jet et-phi distribution
TH1D * fhJetPtRho
! pt spectrum of AKT jets without kt bgk
Definition: External.C:228
void UserCreateOutputObjects()
Overloads base class method. Creates output objects.
Definition: External.C:212
TH2D * fhInvPtQVsPhiASide[2]
! Inverse pt vs phi in away-side rapidity
TH1D * fhConeRho
! local KT bgk
TH1D * fhAtimesRhoCMS
! CMS jet background times area
TH2D * fhJetPhiPt
! Jet phi-pt distribution
TH1D * fhYVertex
! Y vertex
Bool_t fInitializedLocal
! gc trigger if tracks/jets are loaded initiates calling ExecOnce
TH1D * fhKTJetPt
! pt spectrum of KT jets
TH2D * fhSigmaPtOverPtVsPt[2]
! Pt vs sigmaPt/Pt
Enhanced TList-derived class that implements correct merging for pt_hard binned production.
Definition: AliEmcalList.h:25
TH1D * fhCuts
! histogram for pilup/vertex cuts
TH2D * fhInvPtQVsPhi[2]
! Inverse pt vs phi
AliAnalysisUtils * fHelperClass
! gc Vertex selection helper
Bool_t RetrieveEventObjects()
Retrieve common objects from event.
Double_t EstimateBgKTCMS(AliJetContainer *jetCont)
void SetUsePileUpCut(Bool_t val)
Double_t fMinTrackPt
gc Min track pt to be accepted
TH1D * fhJetPtCMSRho
! pt spectrum of AKT jets without CMS bgk
Float_t fASPDCvsTCut
SPD cluster-vs-tracklet cut.
TH1D * fhAtimesRhoMedian
! kt jet background times area
Base task in the EMCAL jet framework.
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:51
TH1D * fhZVertexBC
! Z vertex before cut
const char Option_t
Definition: External.C:48
Bool_t IsEventInAcceptance(AliVEvent *event)
bool Bool_t
Definition: External.C:53
void SetSignalJetMinArea(Double_t minArea)
TH1D * fhGenJetPt
! Pt spectrum of MC jets
Bool_t fUseDefaultVertexCut
trigger if automatic vertex cut from helper class should be done
Container structure for EMCAL clusters.
TString fCentralityType
gc Used centrality estimate (V0A, V0C, V0M, ...)
Bool_t fUsePileUpCut
trigger if pileup cut should be done
Container for jet within the EMCAL jet framework.
void SetAcceptanceWindows(Double_t trackEta, Double_t signalJetRadius)
AliAnalysisTaskJetPP()
Default constructor.
TH1D * fhFakeTrkPt
! Pt spectrum of fake tracks
TH2D * fhJetEtaPt
! Jet eta-pt distribution