AliPhysics  d497547 (d497547)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnalysisTaskV0sInJets.h
Go to the documentation of this file.
1 #ifndef AliAnalysisTaskV0sInJets_cxx
2 #define AliAnalysisTaskV0sInJets_cxx
3 
4 // task for analysis of V0s (K0S, (anti-)Lambda) in charged jets
5 // Author: Vit Kucera (vit.kucera@cern.ch)
6 
7 class TH1D;
8 class TH2D;
9 class THnSparse;
10 class TRandom;
11 class TClonesArray;
12 
13 class AliAODv0;
14 class AliAODVertex;
15 class AliAODJet;
16 
17 #include "AliAnalysisTaskSE.h"
18 #include "THnSparse.h"
19 //#include "AuxFunctions.h"
20 
22 {
23 public:
24  AliAnalysisTaskV0sInJets(); // Default constructor
25  AliAnalysisTaskV0sInJets(const char* name); // Constructor
26  virtual ~AliAnalysisTaskV0sInJets(); // Destructor
27  virtual void UserCreateOutputObjects();
28  virtual void UserExec(Option_t* option);
29  virtual void Terminate(Option_t*) {}
30 
31  void SetTypeAOD(Int_t type = 1) {fiAODAnalysis = type;}
32  void SetIsPbPb(Bool_t val = 1) {fbIsPbPb = val;}
33  void SetJetBranchName(char* line) {fsJetBranchName = line;}
34  void SetJetBgBranchName(char* line) {fsJetBgBranchName = line;}
35  void SetCuts(Double_t z = 10, Double_t r = 1, Double_t cL = 0, Double_t cH = 80) {fdCutVertexZ = z; fdCutVertexR2 = r * r; fdCutCentLow = cL; fdCutCentHigh = cH;}
38  void SetJetRadius(Double_t r = 0.4) {fdRadiusJet = r;}
39  void SetJetRadiusBg(Double_t r = 0.4) {fdRadiusJetBg = r;}
40  void SetJetSelection(Bool_t select = kTRUE) {fbJetSelection = select;}
41  void SetMCAnalysis(Bool_t select = kTRUE) {fbMCAnalysis = select;}
42 // void SetTreeOutput(Bool_t select = kTRUE){fbTreeOutput = select;}
43  void FillQAHistogramV0(AliAODVertex* vtx, const AliAODv0* vZero, Int_t iIndexHisto, Bool_t IsCandK0s, Bool_t IsCandLambda, Bool_t IsInPeakK0s, Bool_t IsInPeakLambda);
44 // virtual Double_t MassPeakSigma(Double_t pt, Int_t particle);
45 // virtual Double_t MassPeakSigma(Int_t iCent, Double_t pt, Int_t particle);
46  void FillCandidates(Double_t mK, Double_t mL, Double_t mAL, Bool_t isK, Bool_t isL, Bool_t isAL, Int_t iCut, Int_t iCent);
47  Bool_t IsParticleInCone(const AliVParticle* part1, const AliVParticle* part2, Double_t dRMax) const; // decides whether a particle is inside a jet cone
48  Bool_t OverlapWithJets(const TClonesArray* array, const AliVParticle* cone, Double_t dDistance) const; // decides whether a cone overlaps with other jets
49  AliAODJet* GetRandomCone(const TClonesArray* array, Double_t dEtaConeMax, Double_t dDistance) const; // generate a random cone which does not overlap with selected jets
50  AliAODJet* GetMedianCluster(const TClonesArray* array, Double_t dEtaConeMax) const; // get median kt cluster
51  Double_t AreaCircSegment(Double_t dRadius, Double_t dDistance) const; // area of circular segment
52 
56  void SetCutCPAMin(Double_t val = 0.998) {fdCutCPAMin = val;}
57  void SetCutNTauMax(Double_t val = 5.) {fdCutNTauMax = val;}
58 
59  Bool_t IsSelectedForJets(AliAODEvent* fAOD, Double_t dVtxZCut, Double_t dVtxR2Cut, Double_t dCentCutLo, Double_t dCentCutUp, Bool_t bCutDeltaZ = kFALSE, Double_t dDeltaZMax = 100.);
64  static bool CompareClusters(const std::vector<Double_t> cluster1, const std::vector<Double_t> cluster2); // compare clusters by their pt/area
65 
66  // upper edges of centrality bins
67  static const Int_t fgkiNBinsCent = 1; // number of centrality bins
68  static const Int_t fgkiCentBinRanges[fgkiNBinsCent]; // upper edges of centrality bins
69  // axis: pT of V0
70  static const Double_t fgkdBinsPtV0[2]; // [GeV/c] minimum and maximum or desired binning of the axis (intended for the rebinned axis)
71  static const Int_t fgkiNBinsPtV0; // number of bins (intended for the rebinned axis)
72  static const Int_t fgkiNBinsPtV0Init; // initial number of bins (uniform binning)
73  // axis: pT of jets
74  static const Double_t fgkdBinsPtJet[2]; // [GeV/c] minimum and maximum or desired binning of the axis (intended for the rebinned axis)
75  static const Int_t fgkiNBinsPtJet; // number of bins (intended for the rebinned axis)
76  static const Int_t fgkiNBinsPtJetInit; // initial number of bins (uniform binning)
77  // axis: K0S invariant mass
78  static const Int_t fgkiNBinsMassK0s; // number of bins (uniform binning)
79  static const Double_t fgkdMassK0sMin; // minimum
80  static const Double_t fgkdMassK0sMax; // maximum
81  // axis: Lambda invariant mass
82  static const Int_t fgkiNBinsMassLambda; // number of bins (uniform binning)
83  static const Double_t fgkdMassLambdaMin; // minimum
84  static const Double_t fgkdMassLambdaMax; // maximum
85 
86 private:
93 // TTree* ftreeOut; //! output tree
94 
95  Int_t fiAODAnalysis; // switch for input AOD/ESD
96  Bool_t fbIsPbPb; // switch Pb-Pb / p-p collisions
97 
98  // V0 selection
99  Double_t fdCutDCAToPrimVtxMin; // [cm] min DCA of daughters to the prim vtx
100  Double_t fdCutDCADaughtersMax; // [sigma of TPC tracking] max DCA between daughters
101  Double_t fdCutNSigmadEdxMax; // [sigma dE/dx] max difference between measured and expected signal of dE/dx in the TPC
102  Double_t fdCutCPAMin; // min cosine of the pointing angle
103  Double_t fdCutNTauMax; // [tau] max proper lifetime in multiples of the mean lifetime
104  // jet selection
105  TString fsJetBranchName; // name of the branch with jets
106  TString fsJetBgBranchName; // name of the branch with kt clusters used for the rho calculation
107  Double_t fdCutPtJetMin; // [GeV/c] minimum jet pt
108  Double_t fdCutPtTrackMin; // [GeV/c] minimum pt of leading jet-track
109  Double_t fdRadiusJet; // R of jet finder used for finding V0s in the jet cone
110  Double_t fdRadiusJetBg; // R of kt jet finder used for reconstruction of bg clusters
111  Bool_t fbJetSelection; // switch for the analysis of V0s in jets
112 
113  Bool_t fbMCAnalysis; // switch for the analysis of simulated data
114 // Bool_t fbTreeOutput; // switch for the output tree
115  TRandom* fRandom;
116 
117  // event cuts
118  Double_t fdCutVertexZ; // [cm] maximum |z| of primary vertex
119  Double_t fdCutVertexR2; // [cm^2] maximum r^2 of primary vertex
120  Double_t fdCutCentLow; // [%] minimum centrality
121  Double_t fdCutCentHigh; // [%] maximum centrality
122  /*
123  // output branches
124  TClonesArray* fBranchV0Rec; //! output branch for reconstructed V0s
125  TClonesArray* fBranchV0Gen; //! output branch for generated V0s
126  TClonesArray* fBranchJet; //! output branch for selected jets
127  AliEventInfoObject* fEventInfo; //! class to store info about events
128  */
130 
131  // event histograms
142 
143  // jet histograms
154 
155  static const Int_t fgkiNCategV0 = 17; // number of V0 selection steps
156 
157  // QA histograms
158  static const Int_t fgkiNQAIndeces = 2; // 0 - before cuts, 1 - after cuts
170 
171  // K0s
182  // K0s Inclusive
184  // K0s Cones
191 
195 // TH1D* fh1DeltaZK0s[fgkiNBinsCent]; //! z-distance between V0 vertex and primary vertex, in centrality bins
196  // MC histograms
197  // inclusive
201  // inclusive eta-pT efficiency
204  // MC daughter eta inclusive
205 // THnSparse* fhnV0K0sInclDaughterEtaPtPtMCGen[fgkiNBinsCent]; //! eta_daughter-pt_daughter-pt_V0 generated
207  // in jets
210  // in jets eta-pT efficiency
213  // MC daughter eta in JC
214 // THnSparse* fhnV0K0sInJetsDaughterEtaPtPtMCGen[fgkiNBinsCent]; //! eta_daughter-pt_daughter-pt_V0 generated
216 
217  // resolution
220 
221  // Lambda
232  // Lambda Inclusive
234  // Lambda Cones
241 
245 // TH1D* fh1DeltaZLambda[fgkiNBinsCent]; //!
246  // MC histograms
247  // inclusive
251  // inclusive eta-pT efficiency
254  // MC daughter eta inclusive
255 // THnSparse* fhnV0LambdaInclDaughterEtaPtPtMCGen[fgkiNBinsCent]; //! eta_daughter-pt_daughter-pt_V0 generated
257  // in jets
260  // in jets eta-pT efficiency
263  // MC daughter eta in JC
264 // THnSparse* fhnV0LambdaInJetsDaughterEtaPtPtMCGen[fgkiNBinsCent]; //! eta_daughter-pt_daughter-pt_V0 generated
266 
267  // resolution
270  // feed-down
275 
276  // ALambda
288  // ALambda Inclusive
290  // ALambda Cones
297 
301 // TH1D* fh1DeltaZALambda[fgkiNBinsCent]; //!
302  // MC histograms
303  // inclusive
308  // inclusive eta-pT efficiency
311  // MC daughter eta inclusive
312 // THnSparse* fhnV0ALambdaInclDaughterEtaPtPtMCGen[fgkiNBinsCent]; //! eta_daughter-pt_daughter-pt_V0 generated
314  // in jets
318  // in jets eta-pT efficiency
321  // MC daughter eta in JC
322 // THnSparse* fhnV0ALambdaInJetsDaughterEtaPtPtMCGen[fgkiNBinsCent]; //! eta_daughter-pt_daughter-pt_V0 generated
324 
325  // resolution
328  // feed-down
333 
345  THnSparse* fh3CCMassCorrelBoth;
346  THnSparse* fh3CCMassCorrelKNotL;
347  THnSparse* fh3CCMassCorrelLNotK;
348 
349  // Cut tuning
350  // crossed/findable, daughter pt, dca, cpa, r, pseudorapidity, y, decay length, PID sigma
351  /*
352  TH2D* fh2CutTPCRowsK0s[fgkiNQAIndeces]; //! inv mass vs TPC rows
353  TH2D* fh2CutTPCRowsLambda[fgkiNQAIndeces]; //!
354  TH2D* fh2CutPtPosK0s[fgkiNQAIndeces]; //! inv mass vs pt of positive daughter
355  TH2D* fh2CutPtNegK0s[fgkiNQAIndeces]; //! inv mass vs pt of negative daughter
356  TH2D* fh2CutPtPosLambda[fgkiNQAIndeces]; //!
357  TH2D* fh2CutPtNegLambda[fgkiNQAIndeces]; //!
358  TH2D* fh2CutDCAVtx[fgkiNQAIndeces]; //! inv mass vs DCA of daughters to prim vtx
359  TH2D* fh2CutDCAV0[fgkiNQAIndeces]; //! inv mass vs DCA between daughters
360  TH2D* fh2CutCos[fgkiNQAIndeces]; //! inv mass vs CPA
361  TH2D* fh2CutR[fgkiNQAIndeces]; //! inv mass vs R
362  TH2D* fh2CutEtaK0s[fgkiNQAIndeces]; //! inv mass vs pseudorapidity
363  TH2D* fh2CutEtaLambda[fgkiNQAIndeces]; //!
364  TH2D* fh2CutRapK0s[fgkiNQAIndeces]; //! inv mass vs rapidity
365  TH2D* fh2CutRapLambda[fgkiNQAIndeces]; //!
366  TH2D* fh2CutCTauK0s[fgkiNQAIndeces]; //! inv mass vs lifetime
367  TH2D* fh2CutCTauLambda[fgkiNQAIndeces]; //!
368  TH2D* fh2CutPIDPosK0s[fgkiNQAIndeces]; //! inv mass vs number of dE/dx sigmas for positive daughter
369  TH2D* fh2CutPIDNegK0s[fgkiNQAIndeces]; //! inv mass vs number of dE/dx sigmas for negative daughter
370  TH2D* fh2CutPIDPosLambda[fgkiNQAIndeces]; //!
371  TH2D* fh2CutPIDNegLambda[fgkiNQAIndeces]; //!
372 
373  TH2D* fh2Tau3DVs2D[fgkiNQAIndeces]; //! pt vs ratio 3D lifetime / 2D lifetime
374  */
375 
376  AliAnalysisTaskV0sInJets(const AliAnalysisTaskV0sInJets&); // not implemented
377  AliAnalysisTaskV0sInJets& operator=(const AliAnalysisTaskV0sInJets&); // not implemented
378 
379  ClassDef(AliAnalysisTaskV0sInJets, 3) // example of analysis
380 };
381 
382 #endif
TH1D * fh1V0ALambdaPtMCRecFalse[fgkiNBinsCent]
Bool_t IsParticleInCone(const AliVParticle *part1, const AliVParticle *part2, Double_t dRMax) const
THnSparseD * fhnV0ALambdaInJetsMCFD[fgkiNBinsCent]
TH1D * fh1V0CandPerEventCentALambda[fgkiNBinsCent]
TList * fOutputListMC
Output list for checking cuts.
THnSparse * fh3V0ALambdaEtaPtMassMCRec[fgkiNBinsCent]
TH1D * fh1V0AXiPtMCGen[fgkiNBinsCent]
TH1D * fh1QAV0R[fgkiNQAIndeces]
cosine of pointing angle (CPA)
Int_t GetCentralityBinIndex(Double_t centrality)
THnSparse * fh3V0K0sInJetPtMassMCRec[fgkiNBinsCent]
pt spectrum of generated K0s in jet
THnSparse * fhnV0OutJetK0s[fgkiNBinsCent]
V0 invariant mass vs V0 pt vs jet pt, in centrality bins.
THnSparse * fhnV0ALambdaInclDaughterEtaPtPtMCRec[fgkiNBinsCent]
TH1D * fh1NMedConeCent
random cone eta-pT
TH1D * fh1DCAOutALambda[fgkiNBinsCent]
double Double_t
Definition: External.C:58
THnSparse * fhnV0InPerpK0s[fgkiNBinsCent]
V0 invariant mass vs V0 pt vs jet pt, in centrality bins.
TH2D * fh2V0LambdaPtMassMCRec[fgkiNBinsCent]
TH2D * fh2QAV0EtaPtALambdaPeak[fgkiNQAIndeces]
THnSparse * fhnV0InJetK0s[fgkiNBinsCent]
V0 inv mass vs pt before and after cuts, in centrality bins.
static const Int_t fgkiNBinsPtV0Init
THnSparse * fhnV0K0sInJetsDaughterEtaPtPtMCRec[fgkiNBinsCent]
mass-eta-pt spectrum of successfully reconstructed K0s in jet
THnSparse * fhnV0InRndALambda[fgkiNBinsCent]
THnSparse * fh3CCMassCorrelKNotL
mass correlation of candidates
TH1D * fh1V0K0sPtMCGen[fgkiNBinsCent]
DCA between daughters of V0 outside jets, in centrality bins.
TH2D * fh2QAV0PtRows[fgkiNQAIndeces]
pseudorapidity vs TPC rows
AliAODEvent * fAODOut
Input AOD event.
TH1D * fh1V0CounterCentALambda[fgkiNBinsCent]
virtual void Terminate(Option_t *)
TH1D * fh1DCAInLambda[fgkiNBinsCent]
TH1D * fh1V0LambdaPtMCRecFalse[fgkiNBinsCent]
THnSparse * fhnV0InMedLambda[fgkiNBinsCent]
TH1D * fh1PhiJet[fgkiNBinsCent]
jet eta-pT
void SetCutNTauMax(Double_t val=5.)
void SetCutNSigmadEdxMax(Double_t val=3.)
void SetJetRadius(Double_t r=0.4)
TList * fOutputListStd
Output AOD event.
static const Double_t fgkdBinsPtJet[2]
virtual void UserExec(Option_t *option)
TH2D * fh2ArmPodALambda[fgkiNQAIndeces]
TH1D * fh1V0CandPerEventCentLambda[fgkiNBinsCent]
THnSparse * fh3V0LambdaEtaPtMassMCRec[fgkiNBinsCent]
centrality
TList * fOutputListQA
Output list for standard analysis results.
static bool CompareClusters(const std::vector< Double_t > cluster1, const std::vector< Double_t > cluster2)
TH1D * fh1NJetPerEvent[fgkiNBinsCent]
jet phi
THnSparse * fh3V0K0sInJetEtaPtMCGen[fgkiNBinsCent]
mass-pt spectrum of successfully reconstructed K0s in jet
static const Double_t fgkdMassK0sMax
Bool_t OverlapWithJets(const TClonesArray *array, const AliVParticle *cone, Double_t dDistance) const
TH1D * fh1QAV0Pt[fgkiNQAIndeces]
TH1D * fh1V0ALambdaPtMCRec[fgkiNBinsCent]
THnSparseD * fhnV0ALambdaBulkMCFD[fgkiNBinsCent]
THnSparse * fh4V0ALambdaInJetEtaPtMassMCRec[fgkiNBinsCent]
static const Int_t fgkiNBinsMassLambda
static const Double_t fgkdMassLambdaMax
TH1D * fh1V0InvMassLambdaCent[fgkiNBinsCent]
TH1D * fh1QAV0TPCFindable[fgkiNQAIndeces]
crossed TPC pad rows
THnSparseD * fhnV0LambdaInclMCFD[fgkiNBinsCent]
TH2D * fh2CCLambda
K0s candidates in Lambda peak.
TH1D * fh1V0InvMassK0sAll[fgkiNCategV0]
number of K0s candidates after various cuts
TH2D * fh2V0PtJetAngleLambda[fgkiNBinsCent]
AliAODJet * GetRandomCone(const TClonesArray *array, Double_t dEtaConeMax, Double_t dDistance) const
Double_t ptMin
TH1D * fh1EventCent
number of events for different selection steps and different centralities
void FillCandidates(Double_t mK, Double_t mL, Double_t mAL, Bool_t isK, Bool_t isL, Bool_t isAL, Int_t iCut, Int_t iCent)
TH2D * fh2V0ALambdaInJetPtMCRec[fgkiNBinsCent]
TH1D * fh1QAV0DCAVtx[fgkiNQAIndeces]
charge
TH2D * fh2ArmPodLambda[fgkiNQAIndeces]
Double_t fdCutVertexZ
random-number generator
THnSparse * fhnV0InRndK0s[fgkiNBinsCent]
V0 invariant mass vs V0 pt vs jet pt, in centrality bins.
THnSparseD * fhnV0LambdaInJetsMCFD[fgkiNBinsCent]
THnSparse * fh3V0LambdaInJetEtaPtMCGen[fgkiNBinsCent]
TH1D * fh1QAV0Status[fgkiNQAIndeces]
TH1D * fh1DCAOutLambda[fgkiNBinsCent]
THnSparse * fh3CCMassCorrelLNotK
mass correlation of candidates
TH1D * fh1V0ALambdaPtMCGen[fgkiNBinsCent]
THnSparse * fhnV0InJetLambda[fgkiNBinsCent]
THnSparse * fhnV0NoJetALambda[fgkiNBinsCent]
void SetCutCPAMin(Double_t val=0.998)
TH1D * fh1VtxZ[fgkiNBinsCent]
number of tracks vs centrality
TH1D * fh1V0InvMassLambdaAll[fgkiNCategV0]
number of Lambda candidates after various cuts
TH2D * fh2V0ALambdaInJetPtMCGen[fgkiNBinsCent]
eta_daughter-pt_daughter-pt_V0 reconstructed
TH2D * fh2QAV0PtPtK0sPeak[fgkiNQAIndeces]
V0 rapidity.
Double_t AreaCircSegment(Double_t dRadius, Double_t dDistance) const
TH1D * fh1PtJet[fgkiNBinsCent]
number of V0 cand per event
TH1D * fh1V0CandPerEventCentK0s[fgkiNBinsCent]
Armenteros-Podolanski.
TH2D * fh2QAV0NClRows[fgkiNQAIndeces]
azimuth vs TPC rows
TH2D * fh2V0K0sMCPtGenPtRec[fgkiNBinsCent]
K0s mass resolution vs pt.
TH2D * fh2V0K0sPtMassMCRec[fgkiNBinsCent]
pt spectrum of all generated K0s in event
THnSparse * fhnV0OutJetALambda[fgkiNBinsCent]
THnSparse * fh4V0LambdaInJetEtaPtMassMCRec[fgkiNBinsCent]
TH1D * fh1DCAInK0s[fgkiNBinsCent]
pt jet vs angle V0-jet, in centrality bins
static const Int_t fgkiNCategV0
area of excluded cones for outside-cones V0s
TH2D * fh2V0K0sInJetPtMCGen[fgkiNBinsCent]
eta_daughter-pt_daughter-pt_V0 reconstructed
THnSparse * fhnV0InclusiveK0s[fgkiNBinsCent]
V0 invariant mass, in centrality bins.
void SetCuts(Double_t z=10, Double_t r=1, Double_t cL=0, Double_t cH=80)
THnSparse * fhnV0InclusiveLambda[fgkiNBinsCent]
Int_t cH
Definition: Combine.C:26
TH2D * fh2EtaPhiRndCone[fgkiNBinsCent]
number of generated random cones in centrality bins
TH2D * fh2V0ALambdaPtMassMCRec[fgkiNBinsCent]
TH1D * fh1V0ALambdaPt[fgkiNBinsCent]
int Int_t
Definition: External.C:63
TH1D * fh1V0XiPtMCGen[fgkiNBinsCent]
THnSparse * fhnV0InMedALambda[fgkiNBinsCent]
TH1D * fh1QAV0RapALambda[fgkiNQAIndeces]
THnSparse * fh3CCMassCorrelBoth
Lambda candidates in K0s peak.
TH1D * fh1QACTau3D[fgkiNQAIndeces]
lifetime calculated in xy
TH1D * fh1NRndConeCent
number of jets per event
TH2D * fh2QAV0PtPtALambdaPeak[fgkiNQAIndeces]
THnSparse * fhnV0InPerpALambda[fgkiNBinsCent]
TH2D * fh2ArmPod[fgkiNQAIndeces]
lifetime calculated in xyz
TH2D * fh2V0K0sMCResolMPt[fgkiNBinsCent]
eta_daughter-pt_daughter-pt_V0 reconstructed
void SetJetSelection(Bool_t select=kTRUE)
Definition: External.C:228
Definition: External.C:212
THnSparse * fhnV0LambdaInclDaughterEtaPtPtMCRec[fgkiNBinsCent]
TH2D * fh2QAV0EtaEtaK0s[fgkiNQAIndeces]
daughters pseudorapidity vs V0 pt, in mass peak
THnSparse * fhnV0InclusiveALambda[fgkiNBinsCent]
THnSparse * fhnV0InPerpLambda[fgkiNBinsCent]
THnSparse * fhnV0OutJetLambda[fgkiNBinsCent]
void SetPtTrackMin(Double_t ptMin=0)
TH1D * fh1QACTau2D[fgkiNQAIndeces]
radial distance between prim vtx and decay vertex
TH2D * fh2QAV0PtPtLambdaPeak[fgkiNQAIndeces]
TH2D * fh2EtaPtJet[fgkiNBinsCent]
jet eta
void SetCutDCADaughtersMax(Double_t val=1.)
void SetMCAnalysis(Bool_t select=kTRUE)
TH1D * fh1AreaExcluded
median-cluster cone eta-phi
TH2D * fh2QAV0PhiRows[fgkiNQAIndeces]
pt vs TPC rows
Int_t fiAODAnalysis
Output list for MC related results.
THnSparse * fhnV0NoJetK0s[fgkiNBinsCent]
V0 invariant mass vs V0 pt, in centrality bins.
void SetPtJetMin(Double_t ptMin=0)
THnSparse * fh4V0K0sInJetEtaPtMassMCRec[fgkiNBinsCent]
eta-pt spectrum of generated K0s in jet
TH2D * fh2V0LambdaEtaPtMCGen[fgkiNBinsCent]
THnSparse * fhnV0LambdaInJetsDaughterEtaPtPtMCRec[fgkiNBinsCent]
static const Int_t fgkiCentBinRanges[fgkiNBinsCent]
TH2D * fh2QAV0EtaEtaALambda[fgkiNQAIndeces]
void SetJetRadiusBg(Double_t r=0.4)
THnSparse * fhnV0NoJetLambda[fgkiNBinsCent]
TH1D * fh1QAV0DCAV0[fgkiNQAIndeces]
DCA of daughters to prim vtx.
THnSparseD * fhnV0LambdaBulkMCFD[fgkiNBinsCent]
TH2D * fh2V0LambdaInJetPtMCGen[fgkiNBinsCent]
eta_daughter-pt_daughter-pt_V0 reconstructed
TH1D * fh1QAV0RapLambda[fgkiNQAIndeces]
TH2D * fh2QAV0PhiPhiLambda[fgkiNQAIndeces]
TH1D * fh1QAV0Cos[fgkiNQAIndeces]
DCA between daughters.
TH2D * fh2V0PtJetAngleK0s[fgkiNBinsCent]
V0 invariant mass vs V0 pt, in centrality bins.
TH2D * fh2EtaPhiMedCone[fgkiNBinsCent]
number of found median-cluster cones in centrality bins
TH2D * fh2QAV0EtaEtaLambda[fgkiNQAIndeces]
THnSparse * fhnV0InJetALambda[fgkiNBinsCent]
TH2D * fh2V0PtJetAngleALambda[fgkiNBinsCent]
TH1D * fh1V0CandPerEvent
xy coordinates of the primary vertex
TH1D * fh1DCAOutK0s[fgkiNBinsCent]
DCA between daughters of V0 inside jets, in centrality bins.
TH2D * fh2QAV0PhiPhiALambda[fgkiNQAIndeces]
static const Double_t fgkdBinsPtV0[2]
TH1D * fh1V0LambdaPtMCGen[fgkiNBinsCent]
TH1D * fh1EventCent2Jets
number of events for different centralities
Double_t MassPeakSigmaOld(Double_t pt, Int_t particle)
THnSparse * fhnV0InMedK0s[fgkiNBinsCent]
V0 invariant mass vs V0 pt vs jet pt, in centrality bins.
TH1D * fh1EtaJet[fgkiNBinsCent]
pt spectra of jets for normalisation of in-jet V0 spectra
TH1D * fh1QAV0TPCRowsFind[fgkiNQAIndeces]
findable clusters
TH1D * fh1EventCent2NoJets
number of events for different centralities
TH1D * fh1V0InvMassALambdaAll[fgkiNCategV0]
number of ALambda candidates after various cuts
TH1D * fh1QAV0RapK0s[fgkiNQAIndeces]
daughters azimuth vs azimuth
TH1D * fh1V0InvMassK0sCent[fgkiNBinsCent]
number of K0s candidates per event, in centrality bins
THnSparse * fh3V0LambdaInJetPtMassMCRec[fgkiNBinsCent]
TH2D * fh2CCK0s
Armenteros-Podolanski.
static const Double_t fgkdMassK0sMin
TH1D * fh1V0K0sPtMCRecFalse[fgkiNBinsCent]
pt-mass spectrum of successfully reconstructed K0s in event
TList * fOutputListCuts
Output list for quality assurance.
TH2D * fh2V0ALambdaMCResolMPt[fgkiNBinsCent]
eta_daughter-pt_daughter-pt_V0 reconstructed
TH2D * fh2ArmPodK0s[fgkiNQAIndeces]
daughters pt vs pt, in mass peak
void FillQAHistogramV0(AliAODVertex *vtx, const AliAODv0 *vZero, Int_t iIndexHisto, Bool_t IsCandK0s, Bool_t IsCandLambda, Bool_t IsInPeakK0s, Bool_t IsInPeakLambda)
const char Option_t
Definition: External.C:48
TH2D * fh2QAV0EtaPtK0sPeak[fgkiNQAIndeces]
V0 invariant mass, selection steps.
THnSparse * fhnV0K0sInclDaughterEtaPtPtMCRec[fgkiNBinsCent]
eta-pt-mass spectrum of successfully reconstructed K0s in event
TH1D * fh1V0CounterCentLambda[fgkiNBinsCent]
K0s generated pt vs reconstructed pt.
THnSparseD * fhnV0ALambdaInclMCFD[fgkiNBinsCent]
TH2D * fh2VtxXY[fgkiNBinsCent]
z coordinate of the primary vertex
THnSparse * fh3V0ALambdaInJetEtaPtMCGen[fgkiNBinsCent]
bool Bool_t
Definition: External.C:53
TH2D * fh2V0ALambdaEtaPtMCGen[fgkiNBinsCent]
THnSparse * fhnV0ALambdaInJetsDaughterEtaPtPtMCRec[fgkiNBinsCent]
TH1D * fh1QAV0Eta[fgkiNQAIndeces]
ratio rows/clusters
TH1D * fh1V0CounterCentK0s[fgkiNBinsCent]
pseudorapidity vs clusters
AliAODJet * GetMedianCluster(const TClonesArray *array, Double_t dEtaConeMax) const
TH2D * fh2V0LambdaMCPtGenPtRec[fgkiNBinsCent]
TH1D * fh1V0InvMassALambdaCent[fgkiNBinsCent]
THnSparse * fhnV0InRndLambda[fgkiNBinsCent]
TH2D * fh2V0ALambdaMCPtGenPtRec[fgkiNBinsCent]
TH1D * fh1QAV0Charge[fgkiNQAIndeces]
pt
TH2D * fh2QAV0EtaNCl[fgkiNQAIndeces]
clusters vs TPC rows
TH1D * fh1QAV0TPCRefit[fgkiNQAIndeces]
online vs offline reconstructed V0 candidates
TH1D * fh1EventCounterCutCent[fgkiNBinsCent]
number of events for different selection steps
TH1D * fh1QAV0TPCRows[fgkiNQAIndeces]
TPC refit on vs off.
Bool_t IsSelectedForJets(AliAODEvent *fAOD, Double_t dVtxZCut, Double_t dVtxR2Cut, Double_t dCentCutLo, Double_t dCentCutUp, Bool_t bCutDeltaZ=kFALSE, Double_t dDeltaZMax=100.)
static const Double_t fgkdMassLambdaMin
THnSparse * fh3V0K0sEtaPtMassMCRec[fgkiNBinsCent]
eta-pt spectrum of all generated K0s in event
TH2D * fh2QAV0EtaRows[fgkiNQAIndeces]
pseudorapidity
TH1D * fh1EventCent2
number of events for different centralities
TH2D * fh2QAV0PhiPhiK0s[fgkiNQAIndeces]
daughters pseudorapidity vs pseudorapidity
static const Int_t fgkiNBinsPtJetInit
TH1D * fh1DCAInALambda[fgkiNBinsCent]
TH2D * fh2V0LambdaMCResolMPt[fgkiNBinsCent]
eta_daughter-pt_daughter-pt_V0 reconstructed
TH2D * fh2QAV0EtaPtLambdaPeak[fgkiNQAIndeces]
TH2D * fh2V0K0sEtaPtMCGen[fgkiNBinsCent]
pt spectrum of false reconstructed K0s in event
AliAnalysisTaskV0sInJets & operator=(const AliAnalysisTaskV0sInJets &)
TH2D * fh2EventCentTracks
number of events for different centralities
void SetCutDCAToPrimVtxMin(Double_t val=0.1)
THnSparse * fh3V0ALambdaInJetPtMassMCRec[fgkiNBinsCent]