]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGDQ/dielectron/AliDielectronVarManager.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
* remove unused var
[u/mrichter/AliRoot.git] / PWGDQ / dielectron / AliDielectronVarManager.h
CommitLineData
b2a297fa 1#ifndef ALIDIELECTRONVARMANAGER_H
2#define ALIDIELECTRONVARMANAGER_H
3/* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
5
6//#############################################################
7//# #
8//# Class AliDielectronVarManager #
9//# Class for management of available variables #
10//# #
11//# Authors: #
12//# Anton Andronic, GSI / A.Andronic@gsi.de #
13//# Ionut C. Arsene, GSI / I.C.Arsene@gsi.de #
14//# Julian Book, Uni Ffm / Julian.Book@cern.ch #
572b0139 15//# Markus Köhler, GSI / M.Koehler@gsi.de #
6551594b 16//# Frederick Kramer, Uni Ffm / Frederick.Kramer@cern.ch #
b2a297fa 17//# Magnus Mager, CERN / Magnus.Mager@cern.ch #
18//# WooJin J. Park, GSI / W.J.Park@gsi.de #
19//# Jens Wiechula, Uni HD / Jens.Wiechula@cern.ch #
20//# #
21//#############################################################
22
b2a297fa 23#include <TNamed.h>
5720c765 24#include <TProfile.h>
a823f01b 25#include <TProfile2D.h>
5720c765 26#include <TH3D.h>
27#include <TFile.h>
28#include <TDatabasePDG.h>
29#include <TKey.h>
29091f0d 30#include <TBits.h>
b2a297fa 31
32#include <AliVEvent.h>
33#include <AliESDEvent.h>
6551594b 34#include <AliAODEvent.h>
35#include <AliMCEvent.h>
36#include <AliVVertex.h>
37#include <AliESDVertex.h>
99345a64 38#include <AliAODVertex.h>
5720c765 39#include <AliEventplane.h>
40
41#include <AliESDVZERO.h>
42#include <AliAODVZERO.h>
b2a297fa 43
44#include <AliVParticle.h>
b2a297fa 45#include <AliESDtrack.h>
46#include <AliAODTrack.h>
572b0139 47#include <AliAODPid.h>
b2a297fa 48#include <AliKFParticle.h>
572b0139 49#include <AliKFVertex.h>
b2a297fa 50#include <AliMCParticle.h>
51#include <AliAODMCParticle.h>
52#include <AliVTrack.h> // ?
53
6551594b 54#include <AliExternalTrackParam.h>
55#include <AliESDpid.h>
4649c75b 56#include <AliCentrality.h>
2a14a7b1 57#include <AliAODpidUtil.h>
8df8e382 58#include <AliPID.h>
ba15fdfb 59#include <AliPIDResponse.h>
6551594b 60
b2a297fa 61#include "AliDielectronPair.h"
8df8e382 62#include "AliDielectronMC.h"
48609e3d 63#include "AliDielectronPID.h"
ba15fdfb 64#include "AliDielectronHelper.h"
b2a297fa 65
8d42b5b9 66#include "AliAnalysisManager.h"
67#include "AliVZEROEPSelectionTask.h"
68
b2a297fa 69class AliVEvent;
70
71//________________________________________________________________
72class AliDielectronVarManager : public TNamed {
73
74public:
8df8e382 75
b2a297fa 76 // Particle specific variables
77 enum ValueTypes {
78 kPx = 0, // px
79 kPy, // py
80 kPz, // pz
81 kPt, // transverse momentum
82 kP, // momentum
83 kXv, // vertex position in x
84 kYv, // vertex position in y
85 kZv, // vertex position in z
86 kOneOverPt, // 1/pt
87 kPhi, // phi angle
88 kTheta, // theta angle
89 kEta, // pseudo-rapidity
90 kY, // rapidity
91 kE, // energy
92 kM, // mass
93 kCharge, // charge
94 kNclsITS, // number of clusters assigned in the ITS
5720c765 95 kITSchi2Cl, // chi2/cl in the ITS
b2a297fa 96 kNclsTPC, // number of clusters assigned in the TPC
5720c765 97 kNclsSTPC, // number of shared clusters assigned in the TPC
98 kNclsSFracTPC, // fraction of shared clusters assigned in the TPC
554e40f8 99 kNclsTPCiter1, // number of clusters assigned in the TPC after first iteration
b2a297fa 100 kNFclsTPC, // number of findable clusters in the TPC
09835c8f 101 kNFclsTPCr, // number of findable clusters(crossed rows) in the TPC with more robust definition
fb7d2d99 102 kNFclsTPCrFrac, // number of found/findable clusters in the TPC with more robust definition
fdda39b6 103 kNFclsTPCfCross, // fraction crossed rows/findable clusters in the TPC, as done in AliESDtrackCuts
b2a297fa 104 kTPCsignalN, // number of points used for dEdx
ba15fdfb 105 kTPCsignalNfrac, // fraction of points used for dEdx / cluster used for tracking
61d106d3 106 kTPCchi2Cl, // chi2/cl in TPC
5720c765 107 kTPCclsDiff, // TPC cluster difference
29091f0d 108 kTPCclsSegments, // TPC cluster segments
61d106d3 109 kTrackStatus, // track status bits
07480df6 110 kFilterBit, // AOD filter bits
111
b2a297fa 112 kNclsTRD, // number of clusters assigned in the TRD
113 kTRDntracklets, // number of TRD tracklets used for tracking/PID TODO: correct getter
114 kTRDpidQuality, // number of TRD tracklets used for PID
67fd1119 115 kTRDchi2, // chi2 in TRD
8df8e382 116 kTRDprobEle, // TRD electron pid probability
117 kTRDprobPio, // TRD electron pid probability
07480df6 118 kTRDprob2DEle, // TRD electron pid probability 2D LQ
119 kTRDprob2DPio, // TRD electron pid probability 2D LQ
5720c765 120 kTRDphi, // Phi angle of the track at the entrance of the TRD
121 kTRDpidEffLeg, // TRD pid efficiency from conversion electrons
35363fb9 122 kTRDsignal, // TRD signal
8df8e382 123
b2a297fa 124 kImpactParXY, // Impact parameter in XY plane
125 kImpactParZ, // Impact parameter in Z
126 kTrackLength, // Track length
8df8e382 127
b2ad74d0 128
5720c765 129 kPdgCode, // PDG code
fdda39b6 130 kPdgCodeMother,
5720c765 131 kPdgCodeGrandMother, // PDG code of the grandmother
1bb1fef1 132 kHasCocktailMother, // true if particle is added via MC generator cocktail (AliDielectronSignal::kDirect)
133 kHasCocktailGrandMother, // true if particle is added via MC generator cocktail (AliDielectronSignal::kDirect)
8df8e382 134 kNumberOfDaughters, // number of daughters
135 kHaveSameMother, // check that particles have the same mother (MC)
fb7d2d99 136 kIsJpsiPrimary, // check if the particle is primary (MC)
5720c765 137 kNumberOfJPsis, // number of generated inclusive jpsis per event (MC)
138 kNumberOfJPsisPrompt, // number of generated prompt jpsis per event (MC)
139 kNumberOfJPsisNPrompt, // number of generated non-prompt jpsis per event (MC)
140
141 kITSsignal, // ITS dE/dx signal
142 kITSsignalSSD1, // SSD1 dE/dx signal
143 kITSsignalSSD2, // SSD2 dE/dx signal
144 kITSsignalSDD1, // SDD1 dE/dx signal
145 kITSsignalSDD2, // SDD2 dE/dx signal
146 kITSclusterMap, // ITS cluster map
147 kITSLayerFirstCls, // No of innermost ITS layer with a cluster of a track
9143d69f 148 kITSnSigmaEle, // number of sigmas to the dE/dx electron line in the ITS
149 kITSnSigmaPio, // number of sigmas to the dE/dx pion line in the ITS
150 kITSnSigmaMuo, // number of sigmas to the dE/dx muon line in the ITS
151 kITSnSigmaKao, // number of sigmas to the dE/dx kaon line in the ITS
152 kITSnSigmaPro, // number of sigmas to the dE/dx proton line in the ITS
8df8e382 153
b2a297fa 154 kPIn, // momentum at inner wall of TPC (if available), used for PID
5720c765 155 kPOut, // momentum at outer wall of TPC, used for TRD studies
156 kYsignedIn, // signed local y at inner wall of TPC
b2a297fa 157 kTPCsignal, // TPC dE/dx signal
ba15fdfb 158
5720c765 159 kTOFsignal, // TOF signal
160 kTOFbeta, // TOF beta
b2ad74d0 161 kTOFPIDBit, // TOF PID bit (1:set, 0:TOF not available)a
35363fb9 162 kTOFmismProb, // and mismatchPorbability as explain in TOF-twiki
ba15fdfb 163
6551594b 164 kTPCnSigmaEle, // number of sigmas to the dE/dx electron line in the TPC
572b0139 165 kTPCnSigmaPio, // number of sigmas to the dE/dx pion line in the TPC
166 kTPCnSigmaMuo, // number of sigmas to the dE/dx muon line in the TPC
167 kTPCnSigmaKao, // number of sigmas to the dE/dx kaon line in the TPC
168 kTPCnSigmaPro, // number of sigmas to the dE/dx proton line in the TPC
169
5720c765 170 kTOFnSigmaEle, // number of sigmas to the electron line in the TOF
572b0139 171 kTOFnSigmaPio, // number of sigmas to the pion line in the TOF
172 kTOFnSigmaMuo, // number of sigmas to the muon line in the TOF
173 kTOFnSigmaKao, // number of sigmas to the kaon line in the TOF
174 kTOFnSigmaPro, // number of sigmas to the proton line in the TOF
8df8e382 175
5720c765 176 kEMCALnSigmaEle, // number of sigmas to the proton line in the TOF
a823f01b 177 kEMCALEoverP, // E over P from EMCAL
88204efa 178 kEMCALE, // E from EMCAL
a823f01b 179 kEMCALNCells, // NCells from EMCAL
180 kEMCALM02, // M02 showershape parameter
181 kEMCALM20, // M20 showershape parameter
182 kEMCALDispersion, // Dispersion paramter
5720c765 183
88e9bae6 184 kV0Index0, // v0 index 0
ffbede40 185 kKinkIndex0, // kink index 0
186
b2a297fa 187 kParticleMax, //
188 // TODO: kRNClusters ??
5720c765 189 // AliDielectronPair specific variables
b2a297fa 190 kChi2NDF = kParticleMax, // Chi^2/NDF
191 kDecayLength, // decay length
192 kR, // distance to the origin
193 kOpeningAngle, // opening angle
2e02dba4 194 kCosPointingAngle, // cosine of the pointing angle
99345a64 195 kArmAlpha, // Armenteros-Podolanski alpha
196 kArmPt, // Armenteros-Podolanski pt
8df8e382 197 // helicity picture: Z-axis is considered the direction of the mother's 3-momentum vector
198 kThetaHE, // theta in mother's rest frame in the helicity picture
199 kPhiHE, // phi in mother's rest frame in the helicity picture
1762f182 200 kThetaSqHE, // squared value of kThetaHE
5720c765 201 kCos2PhiHE, // Cosine of 2*phi in mother's rest frame in the helicity picture
8d7a07b9 202 kCosTilPhiHE, // Shifted phi depending on kThetaHE
8df8e382 203 // Collins-Soper picture: Z-axis is considered the direction of the vectorial difference between
204 // the 3-mom vectors of target and projectile beams
205 kThetaCS, // theta in mother's rest frame in Collins-Soper picture
206 kPhiCS, // phi in mother's rest frame in Collins-Soper picture
1762f182 207 kThetaSqCS, // squared value of kThetaCS
236e1bda 208 kPsiPair, // phi in mother's rest frame in Collins-Soper picture
b2ad74d0 209 kPhivPair, // angle between ee plane and the magnetic field (can be useful for conversion rejection)
210 kPairPlanev0rpH2Angle, // angle between ee plane and VZERO-C reaction plane
211 kPairPlaneMagAngle, // angle between ee plane and strong magnetic field
151ddcbc 212 kPairPlaneAngle, // angle between ee plane and strong magnetic field
213 kRotPairx, //ee plane vector
214 kRotPairy, //ee plane vector
215 kRotPairz, //ee plane vector
216 kCos2PhiCS, // Cosine of 2*phi in mother's rest frame in the Collins-Soper picture
8d7a07b9 217 kCosTilPhiCS, // Shifted phi depending on kThetaCS
5720c765 218 kDeltaPhiV0ArpH2, // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-A
219 kDeltaPhiV0CrpH2, // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-C
220 kDeltaPhiV0ACrpH2, // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-A + V0-C
221 kV0ArpH2FlowV2, // v2 coefficient with respect to the 2nd order reaction plane from V0-A
222 kV0CrpH2FlowV2, // v2 coefficient with respect to the 2nd order reaction plane from V0-C
223 kV0ACrpH2FlowV2, // v2 coefficient with respect to the 2nd order reaction plane from V0-A + V0-C
c315310a 224 kDeltaPhiv0ArpH2, // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-A (EPtask)
225 kDeltaPhiv0CrpH2, // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-C
226 kDeltaPhiv0ACrpH2, // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-AC
227 kv0ArpH2FlowV2, // v2 coefficient with respect to the 2nd order reaction plane from V0-A (EPtask)
228 kv0CrpH2FlowV2, // v2 coefficient with respect to the 2nd order reaction plane from V0-C
229 kv0ACrpH2FlowV2, // v2 coefficient with respect to the 2nd order reaction plane from V0-A + V0-C
230
572b0139 231 kLegDist, // distance of the legs
232 kLegDistXY, // distance of the legs in XY
61d106d3 233 kDeltaEta, // Absolute value of Delta Eta for the legs
234 kDeltaPhi, // Absolute value of Delta Phi for the legs
b2a297fa 235 kMerr, // error of mass calculation
236 kDCA, // distance of closest approach TODO: not implemented yet
237 kPairType, // type of the pair, like like sign ++ unlikesign ...
ba15fdfb 238 kPseudoProperTime, // pseudo proper time
e4339752 239 kPseudoProperTimeErr, // pseudo proper time error
5720c765 240 kPseudoProperTimeResolution, // resolution for pseudo proper decay time (reconstructed - MC truth)
241 kPseudoProperTimePull, // normalizd resolution for pseudo proper time = (reco - MC truth)/dReco
242 kTRDpidEffPair, // TRD pid efficieny from conversion electrons
36d19280 243 kMomAsymDau1, // momentum fraction of daughter1
244 kMomAsymDau2, // momentum fraction of daughter2
b2a297fa 245 kPairMax, //
246 // Event specific variables
6551594b 247 kXvPrim=kPairMax, // prim vertex
248 kYvPrim, // prim vertex
249 kZvPrim, // prim vertex
250 kXRes, // primary vertex x-resolution
251 kYRes, // primary vertex y-resolution
252 kZRes, // primary vertex z-resolution
20b6b7d4 253 kPhiMaxPt, // phi angle of the track with maximum pt
254 kMaxPt, // track with maximum pt
5720c765 255
8d42b5b9 256 //// v0 reaction plane quantities from AliEPSelectionTaks, angles interval [-pi,+pi]
fd6ebd85 257 kv0ArpH2, // VZERO-A reaction plane of the Q vector for 2nd harmonic
258 kv0CrpH2, // reaction plane
259 kv0ACrpH2, // VZERO-AC reaction plane of the Q vector for 2nd harmonic
260 kv0AxH2, // VZERO-A x-component of the Q vector for 2nd harmonic
261 kv0AyH2, // VZERO-A y-component of the Q vector for 2nd harmonic
262 kv0CxH2, // VZERO-C x-component of the Q vector for 2nd harmonic
263 kv0CyH2, // VZERO-C y-component of the Q vector for 2nd harmonic
264 kv0ACxH2, // VZERO-AC x-component of the Q vector for 2nd harmonic
265 kv0ACyH2, // VZERO-AC y-component of the Q vector for 2nd harmonic
20b6b7d4 266 kv0AmagH2, // VZERO-A the Q vectors magnitude for 2nd harmonic
267 kv0CmagH2, // VZERO-A the Q vectors magnitude for 2nd harmonic
268 kv0ACmagH2, // VZERO-A the Q vectors magnitude for 2nd harmonic
fd6ebd85 269 kv0A0rpH2, // VZERO-A 1st ring reaction plane of the Q vector for 2nd harmonic
270 kv0A3rpH2, // VZERO-A last ring reaction plane of the Q vector for 2nd harmonic
271 kv0C0rpH2, // VZERO-C 1st ring reaction plane of the Q vector for 2nd harmonic
272 kv0C3rpH2, // VZERO-C last ring reaction plane of the Q vector for 2nd harmonic
5720c765 273 kv0ATPCDiffH2, // V0A-TPC reaction plane difference for 2nd harmonic
274 kv0CTPCDiffH2, // V0C-TPC reaction plane difference for 2nd harmonic
275 kv0Av0CDiffH2, // V0A-V0C reaction plane difference for 2nd harmonic
1d2d723c 276 kv0Av0C0DiffH2, // V0A-ring 0 ofV0C reaction plane difference for 2nd harmonic
277 kv0Av0C3DiffH2, // V0A-ring 3 ofV0C reaction plane difference for 2nd harmonic
278 kv0Cv0A0DiffH2, // V0C-ring 0 ofV0A reaction plane difference for 2nd harmonic
279 kv0Cv0A3DiffH2, // V0C-ring 3 ofV0A reaction plane difference for 2nd harmonic
280 kv0A0v0A3DiffH2, // V0C-ring 0 ofV0A reaction plane difference for 2nd harmonic
281 kv0C0v0C3DiffH2, // V0C-ring 0 ofV0A reaction plane difference for 2nd harmonic
5720c765 282
283 kMultV0A, // VZERO multiplicity and ADC amplitudes
284 kMultV0C,
285 kMultV0,
286 kAdcV0A,
287 kAdcV0C,
288 kAdcV0,
289 kVZEROchMult,
290 // VZERO reaction plane quantities
291 kV0AxH2=kVZEROchMult+64, // VZERO-A x-component of the Q vector for 2nd harmonic
292 kV0AyH2, // VZERO-A y-component of the Q vector for 2nd harmonic
293 kV0ArpH2, // VZERO-A reaction plane of the Q vector for 2nd harmonic
294 kV0CxH2, // VZERO-C x-component of the Q vector for 2nd harmonic
295 kV0CyH2, // y-component
296 kV0CrpH2, // reaction plane
297 kV0ACxH2, // VZERO-AC x-component of the Q vector for 2nd harmonic
298 kV0ACyH2, // VZERO-AC y-component of the Q vector for 2nd harmonic
299 kV0ACrpH2, // VZERO-AC reaction plane of the Q vector for 2nd harmonic
300 kV0ArpResH2, // 2nd harmonic reaction plane resolution for V0A
301 kV0CrpResH2, // V0C
302 kV0ACrpResH2, // V0A+V0C
303 kV0XaXcH2, // Correlation quantities to check V0 reaction plane quality
304 kV0XaYaH2,
305 kV0XaYcH2,
306 kV0YaXcH2,
307 kV0YaYcH2,
308 kV0XcYcH2,
309 kV0ATPCDiffH2, // V0A-TPC reaction plane difference for 2nd harmonic
310 kV0CTPCDiffH2, // V0C-TPC reaction plane difference for 2nd harmonic
311 kV0AV0CDiffH2, // V0A-V0C reaction plane difference for 2nd harmonic
8d42b5b9 312 // TPC reaction plane quantities, angle interval [0,+pi]
5720c765 313 kTPCxH2, // TPC x-component of the Q vector for 2nd harmonic (corrected)
314 kTPCyH2, // TPC y-component of the Q vector for 2nd harmonic (corrected)
20b6b7d4 315 kTPCmagH2, // TPC reaction plane the Q vectors magnitude for 2nd harmonic (corrected)
316 kTPCrpH2, // TPC reaction plane angle of the Q vector for 2nd harmonic (corrected)
5720c765 317 kTPCsub1xH2, // TPC x-component of the Q vector for 2nd harmonic (corrected, sub event 1)
318 kTPCsub1yH2, // TPC y-component of the Q vector for 2nd harmonic (corrected, sub event 1)
319 kTPCsub1rpH2, // TPC reaction plane of the Q vector for 2nd harmonic (corrected, sub event 1)
320 kTPCsub2xH2, // TPC x-component of the Q vector for 2nd harmonic (corrected, sub event 2)
321 kTPCsub2yH2, // TPC y-component of the Q vector for 2nd harmonic (corrected, sub event 2)
322 kTPCsub2rpH2, // TPC reaction plane of the Q vector for 2nd harmonic (corrected, sub event 2)
323 kTPCsub12DiffH2, // TPC reaction plane difference of sub event 1,2 for 2nd harmonic
324 kTPCsub12DiffH2Sin, // TPC reaction plane difference of sub event 1,2 for 2nd harmonic, sinus term
325
326 kTPCxH2uc, // TPC x-component of the Q vector for 2nd harmonic (uncorrected)
327 kTPCyH2uc, // TPC y-component of the Q vector for 2nd harmonic (uncorrected)
20b6b7d4 328 kTPCmagH2uc, // TPC reaction plane the Q vectors magnitude for 2nd harmonic (uncorrected)
329 kTPCrpH2uc, // TPC reaction plane angle of the Q vector for 2nd harmonic (uncorrected)
5720c765 330 kTPCsub1xH2uc, // TPC x-component of the Q vector for 2nd harmonic (uncorrected, sub event 1)
331 kTPCsub1yH2uc, // TPC y-component of the Q vector for 2nd harmonic (uncorrected, sub event 1)
332 kTPCsub1rpH2uc, // TPC reaction plane of the Q vector for 2nd harmonic (uncorrected, sub event 1)
333 kTPCsub2xH2uc, // TPC x-component of the Q vector for 2nd harmonic (uncorrected, sub event 2)
334 kTPCsub2yH2uc, // TPC y-component of the Q vector for 2nd harmonic (uncorrected, sub event 2)
335 kTPCsub2rpH2uc, // TPC reaction plane of the Q vector for 2nd harmonic (uncorrected, sub event 2)
336 kTPCsub12DiffH2uc, // TPC reaction plane difference of sub event 1,2 for 2nd harmonic (uncorrected)
337
ba15fdfb 338 kNTrk, // number of tracks (or tracklets) TODO: ambiguous
339 kTracks, // ESD tracks TODO: ambiguous
a823f01b 340 kNVtxContrib, // number of primary vertex contibutors
6bd73568 341 kNVtxContribTPC, // number of TPC vertex contibutors
ba15fdfb 342 kNacc, // Number of accepted tracks
41a64bfb 343 kMatchEffITSTPC, // ruff estimate on the ITS-TPC matching efficiceny
5720c765 344 kNaccTrcklts, // number of accepted SPD tracklets in |eta|<1.6
345 kNaccTrcklts0916, // number of accepted SPD tracklets in 0.9<|eta|<1.6
346
347 kNaccTrckltsEsd05, // number of accepted SPD tracklets in |eta|<0.5 (AliESDEvent::EstimateMultiplicity())
348 kNaccTrckltsEsd10, // number of accepted SPD tracklets in |eta|<1.0 (AliESDEvent::EstimateMultiplicity())
349 kNaccTrckltsEsd16, // number of accepted SPD tracklets in |eta|<1.6 (AliESDEvent::EstimateMultiplicity())
350 kNaccTrckltsEsd05Corr, //
351 kNaccTrckltsEsd10Corr, //
352 kNaccTrckltsEsd16Corr, //
353 kNaccItsTpcEsd05, // ITS-TPC tracks + ITS SA complementary tracks + tracklets from unassigned tracklets in |eta|<0.5 (AliESDEvent::EstimateMultiplicity())
354 kNaccItsTpcEsd10, // ITS-TPC tracks + ITS SA complementary tracks + tracklets from unassigned tracklets in |eta|<1.0 (AliESDEvent::EstimateMultiplicity())
355 kNaccItsTpcEsd16, // ITS-TPC tracks + ITS SA complementary tracks + tracklets from unassigned tracklets in |eta|<1.6 (AliESDEvent::EstimateMultiplicity())
356 kNaccItsTpcEsd05Corr, //
357 kNaccItsTpcEsd10Corr, //
358 kNaccItsTpcEsd16Corr, //
359
360 kNaccItsPureEsd05, // ITS SA tracks + tracklets from unassigned tracklets in |eta|<0.5 (AliESDEvent::EstimateMultiplicity())
361 kNaccItsPureEsd10, // ITS SA tracks + tracklets from unassigned tracklets in |eta|<1.0 (AliESDEvent::EstimateMultiplicity())
362 kNaccItsPureEsd16, // ITS SA tracks + tracklets from unassigned tracklets in |eta|<1.6 (AliESDEvent::EstimateMultiplicity())
363 kNaccItsPureEsd05Corr, //
364 kNaccItsPureEsd10Corr, //
365 kNaccItsPureEsd16Corr, //
7a1ddc4b 366 kRefMult, // reference multiplicity (only in AODs) should be Ntrk w/o double counts
367 kRefMultTPConly, // TPC only Reference Multiplicty (AliESDtrackCuts::GetReferenceMultiplicity(&esd, kTRUE))
5720c765 368
369 kNch, // MC true number of charged particles in |eta|<1.6
370 kNch05, // MC true number of charged particles in |eta|<0.5
371 kNch10, // MC true number of charged particles in |eta|<1.0
372
ffbede40 373 kCentrality, // event centrality fraction
a823f01b 374 kCentralitySPD, // centrality using SPD
6551594b 375 kNevents, // event counter
5720c765 376 kRunNumber, // run number
d327d9cd 377 kMixingBin,
b2a297fa 378 kNMaxValues //
379 // TODO: (for A+A) ZDCEnergy, impact parameter, Iflag??
380 };
381
382
383 AliDielectronVarManager();
384 AliDielectronVarManager(const char* name, const char* title);
385 virtual ~AliDielectronVarManager();
386 static void Fill(const TObject* particle, Double_t * const values);
ba15fdfb 387 static void FillVarMCParticle2(const AliVParticle *p1, const AliVParticle *p2, Double_t * const values);
61749cc7 388 static void FillVarVParticle(const AliVParticle *particle, Double_t * const values);
b2a297fa 389
6551594b 390 static void InitESDpid(Int_t type=0);
2a14a7b1 391 static void InitAODpidUtil(Int_t type=0);
5720c765 392 static void InitEstimatorAvg(const Char_t* filename);
393 static void InitTRDpidEffHistograms(const Char_t* filename);
a823f01b 394 static void SetVZEROCalibrationFile(const Char_t* filename) {fgVZEROCalibrationFile = filename;}
395
396 static void SetVZERORecenteringFile(const Char_t* filename) {fgVZERORecenteringFile = filename;}
ba15fdfb 397 static void SetPIDResponse(AliPIDResponse *pidResponse) {fgPIDResponse=pidResponse;}
5720c765 398 static AliPIDResponse* GetPIDResponse() { return fgPIDResponse; }
572b0139 399 static void SetEvent(AliVEvent * const ev);
5720c765 400 static void SetEventData(const Double_t data[AliDielectronVarManager::kNMaxValues]);
2a14a7b1 401 static Bool_t GetDCA(const AliAODTrack *track, Double_t d0z0[2]);
5720c765 402 static void SetTPCEventPlane(AliEventplane *const evplane);
a823f01b 403 static void GetVzeroRP(const AliVEvent* event, Double_t* qvec, Int_t sideOption); // 0- V0A; 1- V0C; 2- V0A+V0C
99345a64 404 static AliAODVertex* GetVertex(const AliAODEvent *event, AliAODVertex::AODVtx_t vtype);
5720c765 405
406 static TProfile* GetEstimatorHistogram(Int_t period, Int_t type) {return fgMultEstimatorAvg[period][type];}
407 static Double_t GetTRDpidEfficiency(Int_t runNo, Double_t centrality, Double_t eta, Double_t trdPhi, Double_t pout, Double_t& effErr);
ba15fdfb 408
409 static const AliKFVertex* GetKFVertex() {return fgKFVertex;}
410
90c1ae2d 411 static const char* GetValueName(Int_t i) { return (i>=0&&i<kNMaxValues)?fgkParticleNames[i][0]:""; }
412 static const char* GetValueLabel(Int_t i) { return (i>=0&&i<kNMaxValues)?fgkParticleNames[i][1]:""; }
413 static const char* GetValueUnit(Int_t i) { return (i>=0&&i<kNMaxValues)?fgkParticleNames[i][2]:""; }
07480df6 414 static UInt_t GetValueType(const char* valname);
ba15fdfb 415 static const Double_t* GetData() {return fgData;}
99345a64 416 static AliVEvent* GetCurrentEvent() {return fgEvent;}
ba15fdfb 417
d327d9cd 418 static Double_t GetValue(ValueTypes var) {return fgData[var];}
419 static void SetValue(ValueTypes var, Double_t val) { fgData[var]=val; }
420
b2a297fa 421private:
422
90c1ae2d 423 static const char* fgkParticleNames[kNMaxValues][3]; //variable names
b2a297fa 424
61749cc7 425
b2a297fa 426 static void FillVarESDtrack(const AliESDtrack *particle, Double_t * const values);
427 static void FillVarAODTrack(const AliAODTrack *particle, Double_t * const values);
8df8e382 428 static void FillVarMCParticle(const AliMCParticle *particle, Double_t * const values);
b2a297fa 429 static void FillVarAODMCParticle(const AliAODMCParticle *particle, Double_t * const values);
430 static void FillVarDielectronPair(const AliDielectronPair *pair, Double_t * const values);
1201a1a9 431 static void FillVarKFParticle(const AliKFParticle *pair, Double_t * const values);
432
6551594b 433 static void FillVarVEvent(const AliVEvent *event, Double_t * const values);
434 static void FillVarESDEvent(const AliESDEvent *event, Double_t * const values);
435 static void FillVarAODEvent(const AliAODEvent *event, Double_t * const values);
436 static void FillVarMCEvent(const AliMCEvent *event, Double_t * const values);
5720c765 437 static void FillVarTPCEventPlane(const AliEventplane *evplane, Double_t * const values);
438
a823f01b 439 static void InitVZEROCalibrationHistograms(Int_t runNo);
440 static void InitVZERORecenteringHistograms(Int_t runNo);
441
ba15fdfb 442 static AliPIDResponse *fgPIDResponse; // PID response object
443 static AliVEvent *fgEvent; // current event pointer
8d42b5b9 444 static AliEventplane *fgTPCEventPlane; // current event tpc plane pointer
ba15fdfb 445 static AliKFVertex *fgKFVertex; // kf vertex
5720c765 446 static TProfile *fgMultEstimatorAvg[4][9]; // multiplicity estimator averages (4 periods x 9 estimators)
447 static Double_t fgTRDpidEffCentRanges[10][4]; // centrality ranges for the TRD pid efficiency histograms
448 static TH3D *fgTRDpidEff[10][4]; // TRD pid efficiencies from conversion electrons
a823f01b 449 static TString fgVZEROCalibrationFile; // file with VZERO channel-by-channel calibrations
450 static TString fgVZERORecenteringFile; // file with VZERO Q-vector averages needed for event plane recentering
451 static TProfile2D *fgVZEROCalib[64]; // 1 histogram per VZERO channel
452 static TProfile2D *fgVZERORecentering[2][2]; // 2 VZERO sides x 2 Q-vector components
2e02dba4 453 static Int_t fgCurrentRun; // current run number
a823f01b 454
ba15fdfb 455 static Double_t fgData[kNMaxValues]; //! data
572b0139 456
b2a297fa 457 AliDielectronVarManager(const AliDielectronVarManager &c);
458 AliDielectronVarManager &operator=(const AliDielectronVarManager &c);
459
460 ClassDef(AliDielectronVarManager,1);
461};
462
463
464//Inline functions
6551594b 465inline void AliDielectronVarManager::Fill(const TObject* object, Double_t * const values)
b2a297fa 466{
467 //
468 // Main function to fill all available variables according to the type of particle
469 //
40875e45 470 if (!object) return;
6551594b 471 if (object->IsA() == AliESDtrack::Class()) FillVarESDtrack(static_cast<const AliESDtrack*>(object), values);
472 else if (object->IsA() == AliAODTrack::Class()) FillVarAODTrack(static_cast<const AliAODTrack*>(object), values);
473 else if (object->IsA() == AliMCParticle::Class()) FillVarMCParticle(static_cast<const AliMCParticle*>(object), values);
474 else if (object->IsA() == AliAODMCParticle::Class()) FillVarAODMCParticle(static_cast<const AliAODMCParticle*>(object), values);
475 else if (object->IsA() == AliDielectronPair::Class()) FillVarDielectronPair(static_cast<const AliDielectronPair*>(object), values);
1201a1a9 476 else if (object->IsA() == AliKFParticle::Class()) FillVarKFParticle(static_cast<const AliKFParticle*>(object),values);
6551594b 477 // Main function to fill all available variables according to the type of event
478
479 else if (object->IsA() == AliVEvent::Class()) FillVarVEvent(static_cast<const AliVEvent*>(object), values);
480 else if (object->IsA() == AliESDEvent::Class()) FillVarESDEvent(static_cast<const AliESDEvent*>(object), values);
481 else if (object->IsA() == AliAODEvent::Class()) FillVarAODEvent(static_cast<const AliAODEvent*>(object), values);
482 else if (object->IsA() == AliMCEvent::Class()) FillVarMCEvent(static_cast<const AliMCEvent*>(object), values);
5720c765 483 else if (object->IsA() == AliEventplane::Class()) FillVarTPCEventPlane(static_cast<const AliEventplane*>(object), values);
1201a1a9 484// else printf(Form("AliDielectronVarManager::Fill: Type %s is not supported by AliDielectronVarManager!", object->ClassName())); //TODO: implement without object needed
b2a297fa 485}
486
487inline void AliDielectronVarManager::FillVarVParticle(const AliVParticle *particle, Double_t * const values)
488{
489 //
490 // Fill track information available in AliVParticle into an array
491 //
492 values[AliDielectronVarManager::kPx] = particle->Px();
493 values[AliDielectronVarManager::kPy] = particle->Py();
494 values[AliDielectronVarManager::kPz] = particle->Pz();
495 values[AliDielectronVarManager::kPt] = particle->Pt();
496 values[AliDielectronVarManager::kP] = particle->P();
497
498 values[AliDielectronVarManager::kXv] = particle->Xv();
499 values[AliDielectronVarManager::kYv] = particle->Yv();
500 values[AliDielectronVarManager::kZv] = particle->Zv();
501
ba15fdfb 502 values[AliDielectronVarManager::kOneOverPt] = (particle->Pt()>1.0e-3 ? particle->OneOverPt() : 0.0);
b2a297fa 503 values[AliDielectronVarManager::kPhi] = particle->Phi();
504 values[AliDielectronVarManager::kTheta] = particle->Theta();
505 values[AliDielectronVarManager::kEta] = particle->Eta();
506 values[AliDielectronVarManager::kY] = particle->Y();
8df8e382 507
b2a297fa 508 values[AliDielectronVarManager::kE] = particle->E();
509 values[AliDielectronVarManager::kM] = particle->M();
510 values[AliDielectronVarManager::kCharge] = particle->Charge();
2a14a7b1 511
512 values[AliDielectronVarManager::kPdgCode] = particle->PdgCode();
513
ba15fdfb 514// if ( fgEvent ) AliDielectronVarManager::Fill(fgEvent, values);
515 for (Int_t i=AliDielectronVarManager::kPairMax; i<AliDielectronVarManager::kNMaxValues; ++i)
516 values[i]=fgData[i];
b2a297fa 517}
518
519inline void AliDielectronVarManager::FillVarESDtrack(const AliESDtrack *particle, Double_t * const values)
520{
521 //
522 // Fill track information available for histogramming into an array
523 //
524
525 // Fill common AliVParticle interface information
526 FillVarVParticle(particle, values);
527
5720c765 528 AliESDtrack *esdTrack=0x0;
529 Double_t origdEdx=particle->GetTPCsignal();
530
531 // apply ETa correction, remove once this is in the tender
d327d9cd 532 esdTrack=const_cast<AliESDtrack*>(particle);
9395992d 533 if (!esdTrack) return;
d327d9cd 534 esdTrack->SetTPCsignal(origdEdx/AliDielectronPID::GetEtaCorr(esdTrack)/AliDielectronPID::GetCorrValdEdx(),esdTrack->GetTPCsignalSigma(),esdTrack->GetTPCsignalN());
5720c765 535
8df8e382 536 Double_t pidProbs[AliPID::kSPECIES];
b2a297fa 537 // Fill AliESDtrack interface specific information
164bfb53 538 Double_t tpcNcls=particle->GetTPCNcls();
5720c765 539 Double_t tpcNclsS = particle->GetTPCnclsS();
540 Double_t itsNcls=particle->GetNcls(0);
ba15fdfb 541 Double_t tpcSignalN=particle->GetTPCsignalN();
9afa0f7e 542 Double_t tpcClusFindable=particle->GetTPCNclsF();
5720c765 543 values[AliDielectronVarManager::kNclsITS] = itsNcls; // TODO: get rid of the plain numbers
61d106d3 544 values[AliDielectronVarManager::kNclsTPC] = tpcNcls; // TODO: get rid of the plain numbers
5720c765 545 values[AliDielectronVarManager::kNclsSTPC] = tpcNclsS;
546 values[AliDielectronVarManager::kNclsSFracTPC] = tpcNcls>0?tpcNclsS/tpcNcls:0;
ffbede40 547 values[AliDielectronVarManager::kNclsTPCiter1] = particle->GetTPCNclsIter1(); // TODO: get rid of the plain numbers
9afa0f7e 548 values[AliDielectronVarManager::kNFclsTPC] = tpcClusFindable;
549 values[AliDielectronVarManager::kNFclsTPCr] = particle->GetTPCClusterInfo(2,1);
550 values[AliDielectronVarManager::kNFclsTPCrFrac] = particle->GetTPCClusterInfo(2);
551 values[AliDielectronVarManager::kNFclsTPCfCross]= (tpcClusFindable>0)?(particle->GetTPCClusterInfo(2,1)/tpcClusFindable):0;
ba15fdfb 552 values[AliDielectronVarManager::kTPCsignalN] = tpcSignalN;
553 values[AliDielectronVarManager::kTPCsignalNfrac]= tpcNcls>0?tpcSignalN/tpcNcls:0;
b2a297fa 554 values[AliDielectronVarManager::kNclsTRD] = particle->GetNcls(2); // TODO: get rid of the plain numbers
555 values[AliDielectronVarManager::kTRDntracklets] = particle->GetTRDntracklets(); // TODO: GetTRDtracklets/GetTRDntracklets?
e4339752 556 values[AliDielectronVarManager::kTRDpidQuality] = particle->GetTRDntrackletsPID();
67fd1119 557 values[AliDielectronVarManager::kTRDchi2] = particle->GetTRDchi2();
35363fb9 558 values[AliDielectronVarManager::kTRDsignal] = particle->GetTRDsignal();
5720c765 559 values[AliDielectronVarManager::kTPCclsDiff] = tpcSignalN-tpcNcls;
29091f0d 560 values[AliDielectronVarManager::kTPCclsSegments] = 0.0;
561 const UChar_t threshold = 5;
562 TBits tpcClusterMap = particle->GetTPCClusterMap();
563 UChar_t n=0; UChar_t j=0;
564 for(UChar_t i=0; i<8; ++i) {
565 n=0;
566 for(j=i*20; j<(i+1)*20 && j<159; ++j) n+=tpcClusterMap.TestBitNumber(j);
567 if(n>=threshold) values[AliDielectronVarManager::kTPCclsSegments] += 1.0;
568 }
61d106d3 569 values[AliDielectronVarManager::kTrackStatus] = (Double_t)particle->GetStatus();
07480df6 570 values[AliDielectronVarManager::kFilterBit] = 0;
61d106d3 571
572 values[AliDielectronVarManager::kTPCchi2Cl] = -1;
573 if (tpcNcls>0) values[AliDielectronVarManager::kTPCchi2Cl] = particle->GetTPCchi2() / tpcNcls;
5720c765 574 values[AliDielectronVarManager::kITSchi2Cl] = -1;
575 if (itsNcls>0) values[AliDielectronVarManager::kITSchi2Cl] = particle->GetITSchi2() / itsNcls;
8df8e382 576 //TRD pidProbs
577 particle->GetTRDpid(pidProbs);
578 values[AliDielectronVarManager::kTRDprobEle] = pidProbs[AliPID::kElectron];
579 values[AliDielectronVarManager::kTRDprobPio] = pidProbs[AliPID::kPion];
ffbede40 580
88e9bae6 581 values[AliDielectronVarManager::kV0Index0] = particle->GetV0Index(0);
ffbede40 582 values[AliDielectronVarManager::kKinkIndex0] = particle->GetKinkIndex(0);
8df8e382 583
b2a297fa 584 Float_t impactParXY, impactParZ;
585 particle->GetImpactParameters(impactParXY, impactParZ);
586 values[AliDielectronVarManager::kImpactParXY] = impactParXY;
587 values[AliDielectronVarManager::kImpactParZ] = impactParZ;
588
8df8e382 589
ba15fdfb 590 values[AliDielectronVarManager::kPdgCode]=-1;
591 values[AliDielectronVarManager::kPdgCodeMother]=-1;
592 values[AliDielectronVarManager::kPdgCodeGrandMother]=-1;
1bb1fef1 593 values[AliDielectronVarManager::kHasCocktailMother]=0;
594 values[AliDielectronVarManager::kHasCocktailGrandMother]=0;
ba15fdfb 595
8df8e382 596 values[AliDielectronVarManager::kNumberOfDaughters]=-999;
597
598 AliDielectronMC *mc=AliDielectronMC::Instance();
599
600 if (mc->HasMC()){
1bb1fef1 601 if (mc->GetMCTrack(particle)) {
ba15fdfb 602 values[AliDielectronVarManager::kPdgCode]=mc->GetMCTrack(particle)->PdgCode();
1bb1fef1 603 Int_t trkLbl = mc->GetMCTrack(particle)->GetLabel();
604 values[AliDielectronVarManager::kHasCocktailMother]=mc->CheckParticleSource(trkLbl, AliDielectronSignalMC::kDirect);
605 }
ba15fdfb 606 AliMCParticle *motherMC=mc->GetMCTrackMother(particle); //mother
607 if (motherMC){
608 values[AliDielectronVarManager::kPdgCodeMother]=motherMC->PdgCode();
1bb1fef1 609 Int_t motherLbl = motherMC->GetLabel();
610 values[AliDielectronVarManager::kHasCocktailGrandMother]=mc->CheckParticleSource(motherLbl, AliDielectronSignalMC::kDirect);
ba15fdfb 611
612 motherMC=mc->GetMCTrackMother(motherMC); //grand motherMC
613 if (motherMC) values[AliDielectronVarManager::kPdgCodeGrandMother]=motherMC->PdgCode();;
614 }
615
8df8e382 616 values[AliDielectronVarManager::kNumberOfDaughters]=mc->NumberOfDaughters(particle);
617 } //if(mc->HasMC())
618
619
620
621 values[AliDielectronVarManager::kITSsignal] = particle->GetITSsignal();
622
623 Double_t itsdEdx[4];
624 particle->GetITSdEdxSamples(itsdEdx);
625
626 values[AliDielectronVarManager::kITSsignalSSD1] = itsdEdx[0];
627 values[AliDielectronVarManager::kITSsignalSSD2] = itsdEdx[1];
628 values[AliDielectronVarManager::kITSsignalSDD1] = itsdEdx[2];
629 values[AliDielectronVarManager::kITSsignalSDD2] = itsdEdx[3];
3505bfad 630 values[AliDielectronVarManager::kITSclusterMap] = particle->GetITSClusterMap();
5720c765 631 values[AliDielectronVarManager::kITSLayerFirstCls] = -1.;
632
633 for (Int_t iC=0; iC<6; iC++) {
634 if (((particle->GetITSClusterMap()) & (1<<(iC))) > 0) {
635 values[AliDielectronVarManager::kITSLayerFirstCls] = iC;
636 break;
637 }
638 }
639
3505bfad 640
b2a297fa 641 values[AliDielectronVarManager::kTrackLength] = particle->GetIntegratedLength();
b2a297fa 642 //dEdx information
643 Double_t mom = particle->GetP();
644 const AliExternalTrackParam *in=particle->GetInnerParam();
5720c765 645 Double_t ysignedIn=-100;
646 if (in) {
647 mom = in->GetP();
648 ysignedIn=particle->Charge()*in->GetY();
649 }
b2a297fa 650 values[AliDielectronVarManager::kPIn]=mom;
5720c765 651 values[AliDielectronVarManager::kYsignedIn]=ysignedIn;
652 const AliExternalTrackParam *out=particle->GetOuterParam();
653 if(out) values[AliDielectronVarManager::kPOut] = out->GetP();
654 else values[AliDielectronVarManager::kPOut] = mom;
655 if(out && fgEvent) {
656 Double_t localCoord[3]={0.0};
657 Bool_t localCoordGood = out->GetXYZAt(298.0, ((AliESDEvent*)fgEvent)->GetMagneticField(), localCoord);
658 values[AliDielectronVarManager::kTRDphi] = (localCoordGood && TMath::Abs(localCoord[0])>1.0e-6 && TMath::Abs(localCoord[1])>1.0e-6 ? TMath::ATan2(localCoord[1], localCoord[0]) : -999.);
659 }
660 if(mc->HasMC() && fgTRDpidEff[0][0]) {
661 Int_t runNo = (fgEvent ? fgEvent->GetRunNumber() : -1);
662 Float_t centrality=-1.0;
663 AliCentrality *esdCentrality = (fgEvent ? fgEvent->GetCentrality() : 0x0);
664 if(esdCentrality) centrality = esdCentrality->GetCentralityPercentile("V0M");
665 Double_t effErr=0.0;
666 values[kTRDpidEffLeg] = GetTRDpidEfficiency(runNo, centrality, values[AliDielectronVarManager::kEta],
667 values[AliDielectronVarManager::kTRDphi],
668 values[AliDielectronVarManager::kPOut], effErr);
669 }
b2a297fa 670 values[AliDielectronVarManager::kTPCsignal]=particle->GetTPCsignal();
ba15fdfb 671
672 values[AliDielectronVarManager::kTOFsignal]=particle->GetTOFsignal();
673
674 Double_t l = particle->GetIntegratedLength(); // cm
675 Double_t t = particle->GetTOFsignal();
676 Double_t t0 = fgPIDResponse->GetTOFResponse().GetTimeZero(); // ps
677
678 if( (l < 360. || l > 800.) || (t <= 0.) || (t0 >999990.0) ) {
679 values[AliDielectronVarManager::kTOFbeta]=0.0;
680 }
681 else {
682 t -= t0; // subtract the T0
683 l *= 0.01; // cm ->m
684 t *= 1e-12; //ps -> s
685
686 Double_t v = l / t;
687 Float_t beta = v / TMath::C();
688 values[AliDielectronVarManager::kTOFbeta]=beta;
689 }
5720c765 690 values[AliDielectronVarManager::kTOFPIDBit]=(particle->GetStatus()&AliESDtrack::kTOFpid? 1: 0);
b2ad74d0 691
692 values[AliDielectronVarManager::kTOFmismProb] = fgPIDResponse->GetTOFMismatchProbability(particle);
ba15fdfb 693
6551594b 694 // nsigma to Electron band
695 // TODO: for the moment we set the bethe bloch parameters manually
696 // this should be changed in future!
61749cc7 697 values[AliDielectronVarManager::kTPCnSigmaEle]=(fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kElectron)-AliDielectronPID::GetCorrVal()-AliDielectronPID::GetCntrdCorr(particle)) / AliDielectronPID::GetWdthCorr(particle);
ba15fdfb 698 values[AliDielectronVarManager::kTPCnSigmaPio]=fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kPion);
699 values[AliDielectronVarManager::kTPCnSigmaMuo]=fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kMuon);
700 values[AliDielectronVarManager::kTPCnSigmaKao]=fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kKaon);
701 values[AliDielectronVarManager::kTPCnSigmaPro]=fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kProton);
702
703 values[AliDielectronVarManager::kITSnSigmaEle]=fgPIDResponse->NumberOfSigmasITS(particle,AliPID::kElectron);
704 values[AliDielectronVarManager::kITSnSigmaPio]=fgPIDResponse->NumberOfSigmasITS(particle,AliPID::kPion);
705 values[AliDielectronVarManager::kITSnSigmaMuo]=fgPIDResponse->NumberOfSigmasITS(particle,AliPID::kMuon);
706 values[AliDielectronVarManager::kITSnSigmaKao]=fgPIDResponse->NumberOfSigmasITS(particle,AliPID::kKaon);
707 values[AliDielectronVarManager::kITSnSigmaPro]=fgPIDResponse->NumberOfSigmasITS(particle,AliPID::kProton);
708
709 values[AliDielectronVarManager::kTOFnSigmaEle]=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kElectron);
710 values[AliDielectronVarManager::kTOFnSigmaPio]=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kPion);
711 values[AliDielectronVarManager::kTOFnSigmaMuo]=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kMuon);
712 values[AliDielectronVarManager::kTOFnSigmaKao]=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kKaon);
713 values[AliDielectronVarManager::kTOFnSigmaPro]=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kProton);
5720c765 714
a823f01b 715 //EMCAL PID information
716 Double_t eop=0;
717 Double_t showershape[4]={0.,0.,0.,0.};
718// values[AliDielectronVarManager::kEMCALnSigmaEle] = fgPIDResponse->NumberOfSigmasEMCAL(particle,AliPID::kElectron);
719 values[AliDielectronVarManager::kEMCALnSigmaEle] = fgPIDResponse->NumberOfSigmasEMCAL(particle,AliPID::kElectron,eop,showershape);
720 values[AliDielectronVarManager::kEMCALEoverP] = eop;
88204efa 721 values[AliDielectronVarManager::kEMCALE] = eop*values[AliDielectronVarManager::kP];
a823f01b 722 values[AliDielectronVarManager::kEMCALNCells] = showershape[0];
723 values[AliDielectronVarManager::kEMCALM02] = showershape[1];
724 values[AliDielectronVarManager::kEMCALM20] = showershape[2];
725 values[AliDielectronVarManager::kEMCALDispersion] = showershape[3];
726
727
5720c765 728 //restore TPC signal if it was changed
729 if (esdTrack) esdTrack->SetTPCsignal(origdEdx,esdTrack->GetTPCsignalSigma(),esdTrack->GetTPCsignalN());
b2a297fa 730}
731
732inline void AliDielectronVarManager::FillVarAODTrack(const AliAODTrack *particle, Double_t * const values)
733{
734 //
735 // Fill track information available for histogramming into an array
736 //
737
738 // Fill common AliVParticle interface information
739 FillVarVParticle(particle, values);
ba15fdfb 740 Double_t tpcNcls=particle->GetTPCNcls();
5720c765 741
742 //GetNclsS not present in AODtrack
743 //Replace with method as soon as available
744 TBits tpcSharedMap = particle->GetTPCSharedMap();
745 Double_t tpcNclsS= tpcSharedMap.CountBits(0)-tpcSharedMap.CountBits(159);
9afa0f7e 746 Double_t tpcClusFindable=particle->GetTPCNclsF();
747
572b0139 748 // Reset AliESDtrack interface specific information
1762f182 749 values[AliDielectronVarManager::kNclsITS] = particle->GetITSNcls();
5720c765 750 values[AliDielectronVarManager::kITSchi2Cl] = -1;
ba15fdfb 751 values[AliDielectronVarManager::kNclsTPC] = tpcNcls;
e4339752 752 values[AliDielectronVarManager::kNclsSTPC] = tpcNclsS;
753 values[AliDielectronVarManager::kNclsSFracTPC] = tpcNcls>0?tpcNclsS/tpcNcls:0;
ba15fdfb 754 values[AliDielectronVarManager::kNclsTPCiter1] = tpcNcls; // not really available in AOD
9afa0f7e 755 values[AliDielectronVarManager::kNFclsTPC] = tpcClusFindable;
756 values[AliDielectronVarManager::kNFclsTPCr] = particle->GetTPCClusterInfo(2,1);
757 values[AliDielectronVarManager::kNFclsTPCrFrac] = particle->GetTPCClusterInfo(2);
758 values[AliDielectronVarManager::kNFclsTPCfCross]= (tpcClusFindable>0)?(particle->GetTPCClusterInfo(2,1)/tpcClusFindable):0;
572b0139 759 values[AliDielectronVarManager::kNclsTRD] = 0;
760 values[AliDielectronVarManager::kTRDntracklets] = 0;
b863bc83 761 values[AliDielectronVarManager::kTRDpidQuality] = particle->GetTRDntrackletsPID();
762 values[AliDielectronVarManager::kTRDchi2] = (particle->GetTRDntrackletsPID()!=0.?particle->GetTRDchi2():-1);
35363fb9 763 values[AliDielectronVarManager::kTRDsignal] = particle->GetTRDsignal();
1ca9f651 764 values[AliDielectronVarManager::kTPCclsSegments] = 0.0;
765 const UChar_t threshold = 5;
766 TBits tpcClusterMap = particle->GetTPCClusterMap();
767 UChar_t n=0; UChar_t j=0;
768 for(UChar_t i=0; i<8; ++i) {
769 n=0;
770 for(j=i*20; j<(i+1)*20 && j<159; ++j) n+=tpcClusterMap.TestBitNumber(j);
771 if(n>=threshold) values[AliDielectronVarManager::kTPCclsSegments] += 1.0;
772 }
773
07480df6 774 values[AliDielectronVarManager::kTPCchi2Cl] = (tpcNcls>0)?particle->Chi2perNDF()*(tpcNcls-5)/tpcNcls:-1.; // it is stored as normalized to tpcNcls-5 (see AliAnalysisTaskESDfilter)
61d106d3 775 values[AliDielectronVarManager::kTrackStatus] = (Double_t)particle->GetStatus();
07480df6 776 values[AliDielectronVarManager::kFilterBit] = (Double_t)particle->GetFilterMap();
777
8df8e382 778 //TRD pidProbs
779 //TODO: set correctly
780 values[AliDielectronVarManager::kTRDprobEle] = 0;
781 values[AliDielectronVarManager::kTRDprobPio] = 0;
67fd1119 782
2a14a7b1 783 values[AliDielectronVarManager::kTPCsignalN] = 0;
ba15fdfb 784 values[AliDielectronVarManager::kTPCsignalNfrac]= 0;
b863bc83 785
2a14a7b1 786 // Fill AliAODTrack interface information
787 //
0c09cae4 788 Int_t v0Index=-1;
789 Int_t kinkIndex=-1;
0c09cae4 790 if (particle->GetProdVertex()) {
0c09cae4 791 v0Index = particle->GetProdVertex()->GetType()==AliAODVertex::kV0 ? 1 : 0;
792 kinkIndex = particle->GetProdVertex()->GetType()==AliAODVertex::kKink ? 1 : 0;
793 }
794 values[AliDielectronVarManager::kV0Index0] = v0Index;
795 values[AliDielectronVarManager::kKinkIndex0] = kinkIndex;
a823f01b 796
2a14a7b1 797 Double_t d0z0[2];
798 GetDCA(particle, d0z0);
799 values[AliDielectronVarManager::kImpactParXY] = d0z0[0];
800 values[AliDielectronVarManager::kImpactParZ] = d0z0[1];
572b0139 801
0c09cae4 802 values[AliDielectronVarManager::kPIn] = 0.;
803 values[AliDielectronVarManager::kTPCsignal] = 0.;
804 values[AliDielectronVarManager::kTPCsignalN] = -1.;
805 values[AliDielectronVarManager::kTPCsignalNfrac] = -1.;
806 values[AliDielectronVarManager::kTPCclsDiff] = -999.;
ba15fdfb 807
808 values[AliDielectronVarManager::kTOFsignal]=0;
29091f0d 809 //values[AliDielectronVarManager::kTOFbeta]=0;
572b0139 810
811 values[AliDielectronVarManager::kTPCnSigmaEle]=0;
812 values[AliDielectronVarManager::kTPCnSigmaPio]=0;
813 values[AliDielectronVarManager::kTPCnSigmaMuo]=0;
814 values[AliDielectronVarManager::kTPCnSigmaKao]=0;
815 values[AliDielectronVarManager::kTPCnSigmaPro]=0;
b2a297fa 816
0c09cae4 817 values[AliDielectronVarManager::kTOFnSigmaEle]=0;
818 values[AliDielectronVarManager::kTOFnSigmaPio]=0;
819 values[AliDielectronVarManager::kTOFnSigmaMuo]=0;
820 values[AliDielectronVarManager::kTOFnSigmaKao]=0;
821 values[AliDielectronVarManager::kTOFnSigmaPro]=0;
822
3505bfad 823 values[AliDielectronVarManager::kITSclusterMap] = particle->GetITSClusterMap();
e4339752 824 values[AliDielectronVarManager::kITSLayerFirstCls] = -1.;
825 for (Int_t iC=0; iC<6; iC++) {
826 if (((particle->GetITSClusterMap()) & (1<<(iC))) > 0) {
827 values[AliDielectronVarManager::kITSLayerFirstCls] = iC;
828 break;
829 }
830 }
831
0c09cae4 832 AliAODPid *pid=const_cast<AliAODPid*>(particle->GetDetPid());
b863bc83 833 if (pid) {
0c09cae4 834 Double_t origdEdx=pid->GetTPCsignal();
835 //overwrite signal
836 pid->SetTPCsignal(origdEdx/AliDielectronPID::GetEtaCorr(particle)/AliDielectronPID::GetCorrValdEdx());
837
572b0139 838 Double_t mom =pid->GetTPCmomentum();
ba15fdfb 839 Double_t tpcSignalN=pid->GetTPCsignalN();
0c09cae4 840 values[AliDielectronVarManager::kTPCsignalN] = tpcSignalN;
841 values[AliDielectronVarManager::kTPCsignalNfrac] = tpcNcls>0?tpcSignalN/tpcNcls:0;
842 values[AliDielectronVarManager::kTPCclsDiff] = tpcSignalN-tpcNcls;
61749cc7 843 Double_t tpcNsigmaEle=(fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kElectron)-AliDielectronPID::GetCntrdCorr(particle))/AliDielectronPID::GetWdthCorr(particle);
ba15fdfb 844 Double_t tpcNsigmaPio=fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kPion);
845 Double_t tpcNsigmaMuo=fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kMuon);
846 Double_t tpcNsigmaKao=fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kKaon);
847 Double_t tpcNsigmaPro=fgPIDResponse->NumberOfSigmasTPC(particle,AliPID::kProton);
572b0139 848
164bfb53 849 values[AliDielectronVarManager::kPIn]=mom;
572b0139 850 values[AliDielectronVarManager::kTPCsignal]=pid->GetTPCsignal();
851
852 values[AliDielectronVarManager::kTPCnSigmaEle]=tpcNsigmaEle;
853 values[AliDielectronVarManager::kTPCnSigmaPio]=tpcNsigmaPio;
854 values[AliDielectronVarManager::kTPCnSigmaMuo]=tpcNsigmaMuo;
855 values[AliDielectronVarManager::kTPCnSigmaKao]=tpcNsigmaKao;
856 values[AliDielectronVarManager::kTPCnSigmaPro]=tpcNsigmaPro;
857
e4339752 858 Double_t prob[AliPID::kSPECIES];
859 fgPIDResponse->ComputeTRDProbability(particle,AliPID::kSPECIES,prob);
35363fb9 860 values[AliDielectronVarManager::kTRDprobEle] = prob[AliPID::kElectron];
861 values[AliDielectronVarManager::kTRDprobPio] = prob[AliPID::kPion];
862 // fgPIDResponse->ComputeTRDProbability(particle,AliPID::kSPECIES,prob, AliTRDPIDResponse::kLQ2D);
863 values[AliDielectronVarManager::kTRDprob2DEle] = prob[AliPID::kElectron];
864 values[AliDielectronVarManager::kTRDprob2DPio] = prob[AliPID::kPion];
e4339752 865
866 values[AliDielectronVarManager::kTOFsignal]=pid->GetTOFsignal();
29091f0d 867 Double32_t expt[5];
868 particle->GetIntegratedTimes(expt); // ps
869 Double_t l = TMath::C()* expt[0]*1e-12; // m
870 Double_t t = pid->GetTOFsignal(); // ps start time subtracted (until v5-02-Rev09)
871 AliTOFHeader* tofH=0x0; // from v5-02-Rev10 on subtract the start time
872 if(fgEvent) tofH = (AliTOFHeader*)fgEvent->GetTOFHeader();
873 if(tofH) t -= fgPIDResponse->GetTOFResponse().GetStartTime(particle->P()); // ps
874
875 if( (l < 360.e-2 || l > 800.e-2) || (t <= 0.) ) {
876 values[AliDielectronVarManager::kTOFbeta] =0;
877 }
878 else {
879 t *= 1e-12; //ps -> s
880
881 Double_t v = l / t;
882 Float_t beta = v / TMath::C();
883 values[AliDielectronVarManager::kTOFbeta]=beta;
884 }
5720c765 885
886 Double_t tofNsigmaEle=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kElectron);
887 Double_t tofNsigmaPio=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kPion);
888 Double_t tofNsigmaMuo=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kMuon);
889 Double_t tofNsigmaKao=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kKaon);
890 Double_t tofNsigmaPro=fgPIDResponse->NumberOfSigmasTOF(particle,AliPID::kProton);
891
892 values[AliDielectronVarManager::kTOFnSigmaEle]=tofNsigmaEle;
893 values[AliDielectronVarManager::kTOFnSigmaPio]=tofNsigmaPio;
894 values[AliDielectronVarManager::kTOFnSigmaMuo]=tofNsigmaMuo;
895 values[AliDielectronVarManager::kTOFnSigmaKao]=tofNsigmaKao;
896 values[AliDielectronVarManager::kTOFnSigmaPro]=tofNsigmaPro;
b2ad74d0 897
61749cc7 898 values[AliDielectronVarManager::kTOFmismProb] = fgPIDResponse->GetTOFMismatchProbability(particle);
b2ad74d0 899
61749cc7 900 pid->SetTPCsignal(origdEdx);
572b0139 901 }
ba15fdfb 902
a823f01b 903 //EMCAL PID information
904 Double_t eop=0;
905 Double_t showershape[4]={0.,0.,0.,0.};
906// values[AliDielectronVarManager::kEMCALnSigmaEle] = fgPIDResponse->NumberOfSigmasEMCAL(particle,AliPID::kElectron);
907 values[AliDielectronVarManager::kEMCALnSigmaEle] = fgPIDResponse->NumberOfSigmasEMCAL(particle,AliPID::kElectron,eop,showershape);
908 values[AliDielectronVarManager::kEMCALEoverP] = eop;
88204efa 909 values[AliDielectronVarManager::kEMCALE] = eop*values[AliDielectronVarManager::kP];
a823f01b 910 values[AliDielectronVarManager::kEMCALNCells] = showershape[0];
911 values[AliDielectronVarManager::kEMCALM02] = showershape[1];
912 values[AliDielectronVarManager::kEMCALM20] = showershape[2];
913 values[AliDielectronVarManager::kEMCALDispersion] = showershape[3];
914
ba15fdfb 915 values[AliDielectronVarManager::kPdgCode]=-1;
916 values[AliDielectronVarManager::kPdgCodeMother]=-1;
917 values[AliDielectronVarManager::kPdgCodeGrandMother]=-1;
1bb1fef1 918 values[AliDielectronVarManager::kHasCocktailMother]=0;
919 values[AliDielectronVarManager::kHasCocktailGrandMother]=0;
920
07480df6 921 values[AliDielectronVarManager::kNumberOfDaughters]=-1;
ba15fdfb 922
923 AliDielectronMC *mc=AliDielectronMC::Instance();
924
925 if (mc->HasMC()){
1bb1fef1 926 if (mc->GetMCTrack(particle)) {
ba15fdfb 927 values[AliDielectronVarManager::kPdgCode]=mc->GetMCTrack(particle)->PdgCode();
1bb1fef1 928 Int_t trkLbl = mc->GetMCTrack(particle)->GetLabel();
929 // printf("trklbl %d for %p->%p \n",trkLbl,particle,mc->GetMCTrack(particle));
930 values[AliDielectronVarManager::kHasCocktailMother]=mc->CheckParticleSource(trkLbl, AliDielectronSignalMC::kDirect);
931 }
ba15fdfb 932 AliAODMCParticle *motherMC=mc->GetMCTrackMother(particle); //mother
933 if (motherMC){
934 values[AliDielectronVarManager::kPdgCodeMother]=motherMC->PdgCode();
1bb1fef1 935 Int_t motherLbl = motherMC->GetLabel();
936 values[AliDielectronVarManager::kHasCocktailGrandMother]=mc->CheckParticleSource(motherLbl, AliDielectronSignalMC::kDirect);
937
ba15fdfb 938 motherMC=mc->GetMCTrackMother(motherMC); //grand motherMC
939 if (motherMC) values[AliDielectronVarManager::kPdgCodeGrandMother]=motherMC->PdgCode();;
940 }
941
942 values[AliDielectronVarManager::kNumberOfDaughters]=mc->NumberOfDaughters(particle);
943 } //if(mc->HasMC())
944
5720c765 945 values[AliDielectronVarManager::kTOFPIDBit]=(particle->GetStatus()&AliESDtrack::kTOFpid? 1: 0);
b2a297fa 946}
947
948inline void AliDielectronVarManager::FillVarMCParticle(const AliMCParticle *particle, Double_t * const values)
949{
950 //
951 // Fill track information available for histogramming into an array
952 //
953
554e40f8 954 values[AliDielectronVarManager::kNclsITS] = 0;
5720c765 955 values[AliDielectronVarManager::kITSchi2Cl] = 0;
554e40f8 956 values[AliDielectronVarManager::kNclsTPC] = 0;
07480df6 957 values[AliDielectronVarManager::kNclsSTPC] = 0;
958 values[AliDielectronVarManager::kNclsSFracTPC] = 0;
554e40f8 959 values[AliDielectronVarManager::kNclsTPCiter1] = 0;
960 values[AliDielectronVarManager::kNFclsTPC] = 0;
fb7d2d99 961 values[AliDielectronVarManager::kNFclsTPCr] = 0;
962 values[AliDielectronVarManager::kNFclsTPCrFrac] = 0;
554e40f8 963 values[AliDielectronVarManager::kNclsTRD] = 0;
964 values[AliDielectronVarManager::kTRDntracklets] = 0;
965 values[AliDielectronVarManager::kTRDpidQuality] = 0;
966 values[AliDielectronVarManager::kTPCchi2Cl] = 0;
967 values[AliDielectronVarManager::kTrackStatus] = 0;
07480df6 968 values[AliDielectronVarManager::kFilterBit] = 0;
554e40f8 969 values[AliDielectronVarManager::kTRDprobEle] = 0;
970 values[AliDielectronVarManager::kTRDprobPio] = 0;
971 values[AliDielectronVarManager::kTPCsignalN] = 0;
5720c765 972 values[AliDielectronVarManager::kTPCclsDiff] = 0;
ba15fdfb 973 values[AliDielectronVarManager::kTPCsignalNfrac] = 0;
554e40f8 974 values[AliDielectronVarManager::kImpactParXY] = 0;
975 values[AliDielectronVarManager::kImpactParZ] = 0;
976 values[AliDielectronVarManager::kPIn] = 0;
5720c765 977 values[AliDielectronVarManager::kYsignedIn] = 0;
554e40f8 978 values[AliDielectronVarManager::kTPCsignal] = 0;
ba15fdfb 979 values[AliDielectronVarManager::kTOFsignal] = 0;
980 values[AliDielectronVarManager::kTOFbeta] = 0;
554e40f8 981 values[AliDielectronVarManager::kTPCnSigmaEle] = 0;
982 values[AliDielectronVarManager::kTPCnSigmaPio] = 0;
983 values[AliDielectronVarManager::kTPCnSigmaMuo] = 0;
984 values[AliDielectronVarManager::kTPCnSigmaKao] = 0;
985 values[AliDielectronVarManager::kTPCnSigmaPro] = 0;
986 values[AliDielectronVarManager::kITSclusterMap] = 0;
987
ba15fdfb 988 values[AliDielectronVarManager::kPdgCode] = -1;
989 values[AliDielectronVarManager::kPdgCodeMother] = -1;
990 values[AliDielectronVarManager::kPdgCodeGrandMother] = -1;
1bb1fef1 991 values[AliDielectronVarManager::kHasCocktailMother]=0;
992 values[AliDielectronVarManager::kHasCocktailGrandMother]=0;
993
b2a297fa 994 // Fill common AliVParticle interface information
995 FillVarVParticle(particle, values);
8df8e382 996
997 AliDielectronMC *mc=AliDielectronMC::Instance();
998
b2a297fa 999 // Fill AliMCParticle interface specific information
1000 values[AliDielectronVarManager::kPdgCode] = particle->PdgCode();
1bb1fef1 1001 Int_t trkLbl = particle->GetLabel();
1002 values[AliDielectronVarManager::kHasCocktailMother]=mc->CheckParticleSource(trkLbl, AliDielectronSignalMC::kDirect);
8df8e382 1003
ba15fdfb 1004 AliMCParticle *motherMC = mc->GetMCTrackMother(particle);
1005 if (motherMC){
1006 values[AliDielectronVarManager::kPdgCodeMother]=motherMC->PdgCode();
1bb1fef1 1007 Int_t motherLbl = motherMC->GetLabel();
1008 values[AliDielectronVarManager::kHasCocktailGrandMother]=mc->CheckParticleSource(motherLbl, AliDielectronSignalMC::kDirect);
1009
ba15fdfb 1010 motherMC=mc->GetMCTrackMother(motherMC); //grand mother
1011 if (motherMC) values[AliDielectronVarManager::kPdgCodeGrandMother]=motherMC->PdgCode();;
1012 }
1013
fb7d2d99 1014 values[AliDielectronVarManager::kIsJpsiPrimary] = mc->IsJpsiPrimary(particle);
8df8e382 1015 values[AliDielectronVarManager::kNumberOfDaughters]=mc->NumberOfDaughters(particle);
b2a297fa 1016}
1017
ba15fdfb 1018
1019inline void AliDielectronVarManager::FillVarMCParticle2(const AliVParticle *p1, const AliVParticle *p2, Double_t * const values) {
1020 //
1021 // fill 2 track information starting from MC legs
1022 //
1023
1024 values[AliDielectronVarManager::kNclsITS] = 0;
5720c765 1025 values[AliDielectronVarManager::kITSchi2Cl] = -1;
ba15fdfb 1026 values[AliDielectronVarManager::kNclsTPC] = 0;
5720c765 1027 values[AliDielectronVarManager::kNclsSTPC] = 0;
1028 values[AliDielectronVarManager::kNclsSFracTPC] = 0;
ba15fdfb 1029 values[AliDielectronVarManager::kNclsTPCiter1] = 0;
1030 values[AliDielectronVarManager::kNFclsTPC] = 0;
1031 values[AliDielectronVarManager::kNFclsTPCr] = 0;
1032 values[AliDielectronVarManager::kNFclsTPCrFrac] = 0;
1033 values[AliDielectronVarManager::kNclsTRD] = 0;
1034 values[AliDielectronVarManager::kTRDntracklets] = 0;
1035 values[AliDielectronVarManager::kTRDpidQuality] = 0;
1036 values[AliDielectronVarManager::kTPCchi2Cl] = 0;
1037 values[AliDielectronVarManager::kTrackStatus] = 0;
07480df6 1038 values[AliDielectronVarManager::kFilterBit] = 0;
ba15fdfb 1039 values[AliDielectronVarManager::kTRDprobEle] = 0;
1040 values[AliDielectronVarManager::kTRDprobPio] = 0;
1041 values[AliDielectronVarManager::kTPCsignalN] = 0;
5720c765 1042 values[AliDielectronVarManager::kTPCclsDiff] = 0;
ba15fdfb 1043 values[AliDielectronVarManager::kTPCsignalNfrac] = 0;
1044 values[AliDielectronVarManager::kImpactParXY] = 0;
1045 values[AliDielectronVarManager::kImpactParZ] = 0;
1046 values[AliDielectronVarManager::kPIn] = 0;
5720c765 1047 values[AliDielectronVarManager::kYsignedIn] = 0;
ba15fdfb 1048 values[AliDielectronVarManager::kTPCsignal] = 0;
1049 values[AliDielectronVarManager::kTPCnSigmaEle] = 0;
1050 values[AliDielectronVarManager::kTPCnSigmaPio] = 0;
1051 values[AliDielectronVarManager::kTPCnSigmaMuo] = 0;
1052 values[AliDielectronVarManager::kTPCnSigmaKao] = 0;
1053 values[AliDielectronVarManager::kTPCnSigmaPro] = 0;
1054 values[AliDielectronVarManager::kITSclusterMap] = 0;
1055
1056 values[AliDielectronVarManager::kPdgCode] = 0;
1057 values[AliDielectronVarManager::kPdgCodeMother] = 0;
1bb1fef1 1058 values[AliDielectronVarManager::kHasCocktailMother]=0;
ba15fdfb 1059
1060 AliDielectronMC *mc=AliDielectronMC::Instance();
1061 AliVParticle* mother=0x0;
1062 Int_t mLabel1 = mc->GetMothersLabel(p1->GetLabel());
1063 Int_t mLabel2 = mc->GetMothersLabel(p2->GetLabel());
1064 if(mLabel1==mLabel2)
1065 mother = mc->GetMCTrackFromMCEvent(mLabel1);
5720c765 1066
1067 values[AliDielectronVarManager::kPseudoProperTime] = -2e10;
1068 if(mother) { // same mother
ba15fdfb 1069 FillVarVParticle(mother, values);
4533e78e 1070 Double_t vtxX, vtxY, vtxZ;
1071 mc->GetPrimaryVertex(vtxX,vtxY,vtxZ);
1072 Double_t lxy = ((mother->Xv()- vtxX) * mother->Px() +
1073 (mother->Yv()- vtxY) * mother->Py() )/mother->Pt();
5720c765 1074 values[AliDielectronVarManager::kPseudoProperTime] = lxy*(TDatabasePDG::Instance()->GetParticle(443)->Mass())/mother->Pt();
1075 }
213ceb97 1076 // AliVParticle part
1077 values[AliDielectronVarManager::kPx] = p1->Px()+p2->Px();
1078 values[AliDielectronVarManager::kPy] = p1->Py()+p2->Py();
1079 values[AliDielectronVarManager::kPz] = p1->Pz()+p2->Pz();
1080 values[AliDielectronVarManager::kPt] = TMath::Sqrt(values[AliDielectronVarManager::kPx]*
1081 values[AliDielectronVarManager::kPx]+
1082 values[AliDielectronVarManager::kPy]*
1083 values[AliDielectronVarManager::kPy]);
1084 values[AliDielectronVarManager::kP] = TMath::Sqrt(values[AliDielectronVarManager::kPt]*
1085 values[AliDielectronVarManager::kPt]+
1086 values[AliDielectronVarManager::kPz]*
1087 values[AliDielectronVarManager::kPz]);
ba15fdfb 1088
213ceb97 1089 values[AliDielectronVarManager::kXv] = 0;
1090 values[AliDielectronVarManager::kYv] = 0;
1091 values[AliDielectronVarManager::kZv] = 0;
ba15fdfb 1092
213ceb97 1093 values[AliDielectronVarManager::kOneOverPt] = (values[AliDielectronVarManager::kPt]>1.0e-6 ? 1.0/values[AliDielectronVarManager::kPt] : 0.0);
1094 values[AliDielectronVarManager::kPhi] = TMath::ATan2(values[AliDielectronVarManager::kPy],values[AliDielectronVarManager::kPx]);
1095 values[AliDielectronVarManager::kTheta] = TMath::ATan2(values[AliDielectronVarManager::kPt],values[AliDielectronVarManager::kPz]);
1096 values[AliDielectronVarManager::kEta] = ((values[AliDielectronVarManager::kP]-values[AliDielectronVarManager::kPz])>1.0e-6 && (values[AliDielectronVarManager::kP]+values[AliDielectronVarManager::kPz])>1.0e-6 ? 0.5*TMath::Log((values[AliDielectronVarManager::kP]+values[AliDielectronVarManager::kPz])/(values[AliDielectronVarManager::kP]-values[AliDielectronVarManager::kPz])) : -9999.);
1097 values[AliDielectronVarManager::kE] = p1->E()+p2->E();
1098 values[AliDielectronVarManager::kY] = ((values[AliDielectronVarManager::kE]-values[AliDielectronVarManager::kPz])>1.0e-6 && (values[AliDielectronVarManager::kE]+values[AliDielectronVarManager::kPz])>1.0e-6 ? 0.5*TMath::Log((values[AliDielectronVarManager::kE]+values[AliDielectronVarManager::kPz])/(values[AliDielectronVarManager::kE]-values[AliDielectronVarManager::kPz])) : -9999.);
1099 values[AliDielectronVarManager::kCharge] = p1->Charge()+p2->Charge();
1100
88e9bae6 1101 values[AliDielectronVarManager::kM] = p1->M()*p1->M()+p2->M()*p2->M()+
1102 2.0*(p1->E()*p2->E()-p1->Px()*p2->Px()-p1->Py()*p2->Py()-p1->Pz()*p2->Pz());
1103 values[AliDielectronVarManager::kM] = (values[AliDielectronVarManager::kM]>1.0e-8 ? TMath::Sqrt(values[AliDielectronVarManager::kM]) : -1.0);
ba15fdfb 1104
213ceb97 1105 if ( fgEvent ) AliDielectronVarManager::Fill(fgEvent, values);
a823f01b 1106
5720c765 1107 values[AliDielectronVarManager::kThetaHE] = AliDielectronPair::ThetaPhiCM(p1,p2,kTRUE, kTRUE);
1108 values[AliDielectronVarManager::kPhiHE] = AliDielectronPair::ThetaPhiCM(p1,p2,kTRUE, kFALSE);
1762f182 1109 values[AliDielectronVarManager::kThetaSqHE] = values[AliDielectronVarManager::kThetaHE] * values[AliDielectronVarManager::kThetaHE];
5720c765 1110 values[AliDielectronVarManager::kCos2PhiHE] = TMath::Cos(2*values[AliDielectronVarManager::kPhiHE]);
1111 values[AliDielectronVarManager::kThetaCS] = AliDielectronPair::ThetaPhiCM(p1,p2,kFALSE, kTRUE);
1112 values[AliDielectronVarManager::kPhiCS] = AliDielectronPair::ThetaPhiCM(p1,p2,kFALSE, kFALSE);
1762f182 1113 values[AliDielectronVarManager::kThetaSqCS] = values[AliDielectronVarManager::kThetaCS] * values[AliDielectronVarManager::kThetaCS];
5720c765 1114 values[AliDielectronVarManager::kCos2PhiCS] = TMath::Cos(2*values[AliDielectronVarManager::kPhiCS]);
8d7a07b9 1115 values[AliDielectronVarManager::kCosTilPhiHE] = (values[AliDielectronVarManager::kThetaHE]>0)?(TMath::Cos(values[AliDielectronVarManager::kPhiHE]-TMath::Pi()/4.)):(TMath::Cos(values[AliDielectronVarManager::kPhiHE]-3*TMath::Pi()/4.));
1116 values[AliDielectronVarManager::kCosTilPhiCS] = (values[AliDielectronVarManager::kThetaCS]>0)?(TMath::Cos(values[AliDielectronVarManager::kPhiCS]-TMath::Pi()/4.)):(TMath::Cos(values[AliDielectronVarManager::kPhiCS]-3*TMath::Pi()/4.));
ba15fdfb 1117}
1118
1119
b2a297fa 1120inline void AliDielectronVarManager::FillVarAODMCParticle(const AliAODMCParticle *particle, Double_t * const values)
1121{
1122 //
1123 // Fill track information available for histogramming into an array
1124 //
1125
554e40f8 1126 values[AliDielectronVarManager::kNclsITS] = 0;
5720c765 1127 values[AliDielectronVarManager::kITSchi2Cl] = -1;
554e40f8 1128 values[AliDielectronVarManager::kNclsTPC] = 0;
5720c765 1129 values[AliDielectronVarManager::kNclsSTPC] = 0;
1130 values[AliDielectronVarManager::kNclsSFracTPC] = 0;
554e40f8 1131 values[AliDielectronVarManager::kNclsTPCiter1] = 0;
1132 values[AliDielectronVarManager::kNFclsTPC] = 0;
1133 values[AliDielectronVarManager::kNclsTRD] = 0;
1134 values[AliDielectronVarManager::kTRDntracklets] = 0;
1135 values[AliDielectronVarManager::kTRDpidQuality] = 0;
1136 values[AliDielectronVarManager::kTPCchi2Cl] = 0;
1137 values[AliDielectronVarManager::kTrackStatus] = 0;
07480df6 1138 values[AliDielectronVarManager::kFilterBit] = 0;
554e40f8 1139 values[AliDielectronVarManager::kTRDprobEle] = 0;
1140 values[AliDielectronVarManager::kTRDprobPio] = 0;
1141 values[AliDielectronVarManager::kTPCsignalN] = 0;
5720c765 1142 values[AliDielectronVarManager::kTPCclsDiff] = 0;
ba15fdfb 1143 values[AliDielectronVarManager::kTPCsignalNfrac]= 0;
554e40f8 1144 values[AliDielectronVarManager::kImpactParXY] = 0;
1145 values[AliDielectronVarManager::kImpactParZ] = 0;
1146 values[AliDielectronVarManager::kPIn] = 0;
5720c765 1147 values[AliDielectronVarManager::kYsignedIn] = 0;
554e40f8 1148 values[AliDielectronVarManager::kTPCsignal] = 0;
1149 values[AliDielectronVarManager::kTPCnSigmaEle] = 0;
1150 values[AliDielectronVarManager::kTPCnSigmaPio] = 0;
1151 values[AliDielectronVarManager::kTPCnSigmaMuo] = 0;
1152 values[AliDielectronVarManager::kTPCnSigmaKao] = 0;
1153 values[AliDielectronVarManager::kTPCnSigmaPro] = 0;
1154 values[AliDielectronVarManager::kITSclusterMap] = 0;
1155
ba15fdfb 1156 values[AliDielectronVarManager::kPdgCode] = -1;
1157 values[AliDielectronVarManager::kPdgCodeMother] = -1;
1158 values[AliDielectronVarManager::kPdgCodeGrandMother] = -1;
1bb1fef1 1159 values[AliDielectronVarManager::kHasCocktailMother]=0;
1160 values[AliDielectronVarManager::kHasCocktailGrandMother]=0;
1161
b2a297fa 1162 // Fill common AliVParticle interface information
1163 FillVarVParticle(particle, values);
554e40f8 1164
8df8e382 1165 AliDielectronMC *mc=AliDielectronMC::Instance();
1166
1167
b2a297fa 1168 // Fill AliAODMCParticle interface specific information
8df8e382 1169 values[AliDielectronVarManager::kPdgCode] = particle->PdgCode();
1bb1fef1 1170 Int_t trkLbl = particle->GetLabel();
1171 values[AliDielectronVarManager::kHasCocktailMother]=mc->CheckParticleSource(trkLbl, AliDielectronSignalMC::kDirect);
8df8e382 1172
ba15fdfb 1173 AliAODMCParticle *motherMC = mc->GetMCTrackMother(particle);
1174 if (motherMC){
1175 values[AliDielectronVarManager::kPdgCodeMother]=motherMC->PdgCode();
1bb1fef1 1176 Int_t motherLbl = motherMC->GetLabel();
1177 values[AliDielectronVarManager::kHasCocktailGrandMother]=mc->CheckParticleSource(motherLbl, AliDielectronSignalMC::kDirect);
1178
ba15fdfb 1179 motherMC=mc->GetMCTrackMother(motherMC); //grand mother
1180 if (motherMC) values[AliDielectronVarManager::kPdgCodeGrandMother]=motherMC->PdgCode();;
1181 }
1182
1183
fb7d2d99 1184 values[AliDielectronVarManager::kIsJpsiPrimary] = mc->IsJpsiPrimary(particle);
8df8e382 1185
1186 values[AliDielectronVarManager::kNumberOfDaughters]=mc->NumberOfDaughters(particle);
b2a297fa 1187}
1188
1189inline void AliDielectronVarManager::FillVarDielectronPair(const AliDielectronPair *pair, Double_t * const values)
1190{
1191 //
1192 // Fill pair information available for histogramming into an array
1193 //
2a14a7b1 1194
ba15fdfb 1195 values[AliDielectronVarManager::kPdgCode]=-1;
1196 values[AliDielectronVarManager::kPdgCodeMother]=-1;
1197 values[AliDielectronVarManager::kPdgCodeGrandMother]=-1;
1bb1fef1 1198 values[AliDielectronVarManager::kHasCocktailMother]=0;
1199 values[AliDielectronVarManager::kHasCocktailGrandMother]=0;
08b801a6 1200
1201 Double_t errPseudoProperTime2 = -1;
b2a297fa 1202 // Fill common AliVParticle interface information
1203 FillVarVParticle(pair, values);
1204
1205 // Fill AliDielectronPair specific information
1206 const AliKFParticle &kfPair = pair->GetKFParticle();
61d106d3 1207
1208 Double_t thetaHE=0;
1209 Double_t phiHE=0;
1210 Double_t thetaCS=0;
1211 Double_t phiCS=0;
1212
1213 pair->GetThetaPhiCM(thetaHE,phiHE,thetaCS,phiCS);
5720c765 1214
b2a297fa 1215 values[AliDielectronVarManager::kChi2NDF] = kfPair.GetChi2()/kfPair.GetNDF();
1216 values[AliDielectronVarManager::kDecayLength] = kfPair.GetDecayLength();
1217 values[AliDielectronVarManager::kR] = kfPair.GetR();
1218 values[AliDielectronVarManager::kOpeningAngle] = pair->OpeningAngle();
2e02dba4 1219 values[AliDielectronVarManager::kCosPointingAngle] = fgEvent ? pair->GetCosPointingAngle(fgEvent->GetPrimaryVertex()) : -1;
61d106d3 1220 values[AliDielectronVarManager::kThetaHE] = thetaHE;
1221 values[AliDielectronVarManager::kPhiHE] = phiHE;
1762f182 1222 values[AliDielectronVarManager::kThetaSqHE] = thetaHE * thetaHE;
5720c765 1223 values[AliDielectronVarManager::kCos2PhiHE] = TMath::Cos(2.0*phiHE);
8d7a07b9 1224 values[AliDielectronVarManager::kCosTilPhiHE] = (thetaHE>0)?(TMath::Cos(phiHE-TMath::Pi()/4.)):(TMath::Cos(phiHE-3*TMath::Pi()/4.));
61d106d3 1225 values[AliDielectronVarManager::kThetaCS] = thetaCS;
1226 values[AliDielectronVarManager::kPhiCS] = phiCS;
1762f182 1227 values[AliDielectronVarManager::kThetaSqCS] = thetaCS * thetaCS;
5720c765 1228 values[AliDielectronVarManager::kCos2PhiCS] = TMath::Cos(2.0*phiCS);
8d7a07b9 1229 values[AliDielectronVarManager::kCosTilPhiCS] = (thetaCS>0)?(TMath::Cos(phiCS-TMath::Pi()/4.)):(TMath::Cos(phiCS-3*TMath::Pi()/4.));
572b0139 1230 values[AliDielectronVarManager::kLegDist] = pair->DistanceDaughters();
1231 values[AliDielectronVarManager::kLegDistXY] = pair->DistanceDaughtersXY();
61d106d3 1232 values[AliDielectronVarManager::kDeltaEta] = pair->DeltaEta();
1233 values[AliDielectronVarManager::kDeltaPhi] = pair->DeltaPhi();
8df8e382 1234 values[AliDielectronVarManager::kMerr] = kfPair.GetErrMass()>1e-30&&kfPair.GetMass()>1e-30?kfPair.GetErrMass()/kfPair.GetMass():1000000;
b2a297fa 1235 values[AliDielectronVarManager::kPairType] = pair->GetType();
99345a64 1236 // Armenteros-Podolanski quantities
1237 values[AliDielectronVarManager::kArmAlpha] = pair->GetArmAlpha();
1238 values[AliDielectronVarManager::kArmPt] = pair->GetArmPt();
8df8e382 1239
91487037 1240 values[AliDielectronVarManager::kPsiPair] = fgEvent ? pair->PsiPair(fgEvent->GetMagneticField()) : -5;
805bd069 1241 values[AliDielectronVarManager::kPhivPair] = fgEvent ? pair->PhivPair(fgEvent->GetMagneticField()) : -5;
e4339752 1242 values[AliDielectronVarManager::kPseudoProperTime] = fgEvent ? kfPair.GetPseudoProperDecayTime(*(fgEvent->GetPrimaryVertex()), TDatabasePDG::Instance()->GetParticle(443)->Mass(), &errPseudoProperTime2 ) : -1e10;
5720c765 1243 // values[AliDielectronVarManager::kPseudoProperTime] = fgEvent ? pair->GetPseudoProperTime(fgEvent->GetPrimaryVertex()): -1e10;
e4339752 1244 values[AliDielectronVarManager::kPseudoProperTimeErr] = (errPseudoProperTime2 > 0) ? TMath::Sqrt(errPseudoProperTime2) : -1e10;
5720c765 1245
08b801a6 1246
1247 if (!(pair->GetKFUsage())) {
1248 //if KF Pairing is not enabled, overwrite values that can be easily derived from legs
1249 //use the INDIVIDUAL KF particles as source, which should be a copy of the corresponding properties
1250 //the ESDtrack, the reference to the ESDtrack is not (always) accessible in Mixing, while KF
1251 //particles are copied in the Pair-Object
1252 static const Double_t mElectron = AliPID::ParticleMass(AliPID::kElectron); // MeV
1253
1254 const AliKFParticle& fD1 = pair->GetKFFirstDaughter();
1255 const AliKFParticle& fD2 = pair->GetKFSecondDaughter();
1256
1257 //Define local buffer variables for leg properties
1258 Double_t px1=-9999.,py1=-9999.,pz1=-9999.;
1259 Double_t px2=-9999.,py2=-9999.,pz2=-9999.;
1260 Double_t e1 =-9999.,e2 =-9999.;
20b6b7d4 1261 Double_t feta1=-9999.;//,fphi1=-9999.;
1262 Double_t feta2=-9999.;//,fphi2=-9999.;
08b801a6 1263
1264 px1 = fD1.GetPx();
1265 py1 = fD1.GetPy();
1266 pz1 = fD1.GetPz();
1267 feta1 = fD1.GetEta();
20b6b7d4 1268 // fphi1 = fD1.GetPhi();
08b801a6 1269
1270 px2 = fD2.GetPx();
1271 py2 = fD2.GetPy();
1272 pz2 = fD2.GetPz();
1273 feta2 = fD2.GetEta();
20b6b7d4 1274 // fphi2 = fD2.GetPhi();
08b801a6 1275
1276 //Calculate Energy per particle by hand
1277 e1 = TMath::Sqrt(mElectron*mElectron+px1*px1+py1*py1+pz1*pz1);
1278 e2 = TMath::Sqrt(mElectron*mElectron+px2*px2+py2*py2+pz2*pz2);
1279
1280 //Now Create TLorentzVector:
1281 TLorentzVector lv1,lv2;
1282 lv1.SetPxPyPzE(px1,py1,pz1,e1);
1283 lv2.SetPxPyPzE(px2,py2,pz2,e2);
1284
1285 values[AliDielectronVarManager::kPx] = (lv1+lv2).Px();
1286 values[AliDielectronVarManager::kPy] = (lv1+lv2).Py();
1287 values[AliDielectronVarManager::kPz] = (lv1+lv2).Pz();
1288
1289 values[AliDielectronVarManager::kPt] = (lv1+lv2).Pt();
1290
1291 values[AliDielectronVarManager::kP] = (lv1+lv2).P();
1292
1293 //Not overwritten, could take event vertex in next iteration
1294 values[AliDielectronVarManager::kXv] = (lv1+lv2).X();
1295 values[AliDielectronVarManager::kYv] = (lv1+lv2).Y();
1296 values[AliDielectronVarManager::kZv] = (lv1+lv2).Z();
1297
1298 values[AliDielectronVarManager::kE] = (lv1+lv2).E();
1299
1300
1301 values[AliDielectronVarManager::kM] = (lv1+lv2).M();
1302
1303 values[AliDielectronVarManager::kOpeningAngle] = lv1.Angle(lv2.Vect());
1304
1305 values[AliDielectronVarManager::kOneOverPt] = (values[AliDielectronVarManager::kPt]>0. ? 1./values[AliDielectronVarManager::kPt] : -9999.);
1306 values[AliDielectronVarManager::kPhi] = (lv1+lv2).Phi();
1307 values[AliDielectronVarManager::kEta] = (lv1+lv2).Eta();
1308
1309 values[AliDielectronVarManager::kY] = (lv1+lv2).Rapidity();
1310
1311 for (Int_t i=AliDielectronVarManager::kPairMax; i<AliDielectronVarManager::kNMaxValues; ++i)
1312 values[i]=fgData[i];
1313
1314 // Fill AliDielectronPair specific information
1315 values[AliDielectronVarManager::kDeltaEta] = TMath::Abs(feta1 -feta2 );
1316 values[AliDielectronVarManager::kDeltaPhi] = lv1.DeltaPhi(lv2);
1317 values[AliDielectronVarManager::kPairType] = pair->GetType();
1318
1319 /*
1320 //Also not overwritten, still coming from KF particle
1321 //where needed to be replaced by independent determination
1322 values[AliDielectronVarManager::kCharge] = 0.;
1323 values[AliDielectronVarManager::kPdgCode] = 0.;
1324 values[AliDielectronVarManager::kChi2NDF] = 0.;
1325 values[AliDielectronVarManager::kDecayLength] = 0.;
1326 values[AliDielectronVarManager::kR] = 0.;
1327 values[AliDielectronVarManager::kCosPointingAngle] = 0.;
1328 values[AliDielectronVarManager::kThetaHE] = 0.;
1329 values[AliDielectronVarManager::kPhiHE] = 0.;
1330 values[AliDielectronVarManager::kThetaSqHE] = 0.;
1331 values[AliDielectronVarManager::kCos2PhiHE] = 0.;
1332 values[AliDielectronVarManager::kCosTilPhiHE] = 0.;
1333 values[AliDielectronVarManager::kThetaCS] = 0.;
1334 values[AliDielectronVarManager::kPhiCS] = 0.;
1335 values[AliDielectronVarManager::kThetaSqCS] = 0.;
1336 values[AliDielectronVarManager::kCos2PhiCS] = 0.;
1337 values[AliDielectronVarManager::kCosTilPhiCS] = 0.;
1338 values[AliDielectronVarManager::kLegDist] = 0.;
1339 values[AliDielectronVarManager::kLegDistXY] = 0.;
1340 values[AliDielectronVarManager::kMerr] = 0.;
1341 values[AliDielectronVarManager::kPseudoProperTime] = 0.;
1342 values[AliDielectronVarManager::kPseudoProperTimeErr] = 0.;
1343 //Fill in Taku's PhiV?
1344 values[AliDielectronVarManager::kPsiPair] = 0.;
1345
1346 */
08b801a6 1347 }
1348 //common, regardless of calculation method
1349 // Flow quantities
5720c765 1350 Double_t delta=0.0;
1351 // v2 with respect to VZERO-A event plane
1352 delta = values[AliDielectronVarManager::kPhi] - fgData[AliDielectronVarManager::kV0ArpH2];
1353 if(delta>TMath::Pi()) delta -= 2.0*TMath::Pi(); // keep the [-pi,+pi] interval
1354 if(delta<-1.0*TMath::Pi()) delta += 2.0*TMath::Pi();
1355 values[AliDielectronVarManager::kV0ArpH2FlowV2] = TMath::Cos(2.0*delta); // 2nd harmonic flow coefficient
1356 values[AliDielectronVarManager::kDeltaPhiV0ArpH2] = delta;
1357 // v2 with respect to VZERO-C event plane
1358 delta = values[AliDielectronVarManager::kPhi] - fgData[AliDielectronVarManager::kV0CrpH2];
1359 if(delta>TMath::Pi()) delta -= 2.0*TMath::Pi(); // keep the [-pi,+pi] interval
1360 if(delta<-1.0*TMath::Pi()) delta += 2.0*TMath::Pi();
1361 values[AliDielectronVarManager::kV0CrpH2FlowV2] = TMath::Cos(2.0*delta); // 2nd harmonic flow coefficient
1362 values[AliDielectronVarManager::kDeltaPhiV0CrpH2] = delta;
1363 // v2 with respect to the combined VZERO-A and VZERO-C event plane
1364 delta = values[AliDielectronVarManager::kPhi] - fgData[AliDielectronVarManager::kV0ACrpH2];
1365 if(delta>TMath::Pi()) delta -= 2.0*TMath::Pi(); // keep the [-pi,+pi] interval
1366 if(delta<-1.0*TMath::Pi()) delta += 2.0*TMath::Pi();
1367 values[AliDielectronVarManager::kV0ACrpH2FlowV2] = TMath::Cos(2.0*delta); // 2nd harmonic flow coefficient
1368 values[AliDielectronVarManager::kDeltaPhiV0ACrpH2] = delta;
1369
a3cde8ca 1370
1371 // quantities using the values of AliEPSelectionTask
c315310a 1372 values[AliDielectronVarManager::kDeltaPhiv0ArpH2] = values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0ArpH2];
1373 values[AliDielectronVarManager::kDeltaPhiv0CrpH2] = values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0CrpH2];
a3cde8ca 1374 values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] = values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0ACrpH2];
c315310a 1375 values[AliDielectronVarManager::kv0ACrpH2FlowV2] = TMath::Cos( 2*(values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0ACrpH2]) );
1376 values[AliDielectronVarManager::kv0ArpH2FlowV2] = TMath::Cos( 2*(values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0ArpH2]) );
1377 values[AliDielectronVarManager::kv0CrpH2FlowV2] = TMath::Cos( 2*(values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0CrpH2]) );
1378
1379 // keep the interval [-pi,+pi]
a3cde8ca 1380 if ( values[AliDielectronVarManager::kDeltaPhiv0ArpH2] > TMath::Pi() )
1381 values[AliDielectronVarManager::kDeltaPhiv0ArpH2] -= TMath::TwoPi();
1382 if ( values[AliDielectronVarManager::kDeltaPhiv0CrpH2] > TMath::Pi() )
1383 values[AliDielectronVarManager::kDeltaPhiv0CrpH2] -= TMath::TwoPi();
1384 if ( values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] > TMath::Pi() )
1385 values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] -= TMath::TwoPi();
1386
1387 if ( values[AliDielectronVarManager::kDeltaPhiv0ArpH2] < -1.*TMath::Pi() )
1388 values[AliDielectronVarManager::kDeltaPhiv0ArpH2] += TMath::TwoPi();
1389 if ( values[AliDielectronVarManager::kDeltaPhiv0CrpH2] < -1.*TMath::Pi() )
1390 values[AliDielectronVarManager::kDeltaPhiv0CrpH2] += TMath::TwoPi();
1391 if ( values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] < -1.*TMath::Pi() )
1392 values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] += TMath::TwoPi();
1393
b2ad74d0 1394 //angle between ee plane and Mag/Reaction plane
1395 values[AliDielectronVarManager::kPairPlanev0rpH2Angle] = pair->PairPlanev0rpH2Angle(values[AliDielectronVarManager::kv0CrpH2]);
1396 values[AliDielectronVarManager::kPairPlaneMagAngle] = pair->PairPlaneMagAngle(values[AliDielectronVarManager::kv0CrpH2]);
151ddcbc 1397 values[AliDielectronVarManager::kPairPlaneAngle] = pair->PairPlaneAngle(values[AliDielectronVarManager::kv0CrpH2]);
1398 //ee plane vector
1399 Double_t RotPairx = 0;
1400 Double_t RotPairy = 0;
1401 Double_t RotPairz = 0;
1402 pair->GetRotPair(RotPairx,RotPairy,RotPairz);
1403 values[AliDielectronVarManager::kRotPairx]=RotPairx;
1404 values[AliDielectronVarManager::kRotPairy]=RotPairy;
1405 values[AliDielectronVarManager::kRotPairz]=RotPairz;
b2ad74d0 1406
c315310a 1407
8df8e382 1408 AliDielectronMC *mc=AliDielectronMC::Instance();
1409
1410 if (mc->HasMC()){
5720c765 1411 values[AliDielectronVarManager::kPseudoProperTimeResolution] = -10.0e+10;
ba15fdfb 1412 Bool_t samemother = mc->HaveSameMother(pair);
1413 values[AliDielectronVarManager::kIsJpsiPrimary] = mc->IsJpsiPrimary(pair);
1414 values[AliDielectronVarManager::kHaveSameMother] = samemother ;
5720c765 1415
1bb1fef1 1416 // fill kPseudoProperTimeResolution
5720c765 1417 values[AliDielectronVarManager::kPseudoProperTimeResolution] = -1e10;
1418 // values[AliDielectronVarManager::kPseudoProperTimePull] = -1e10;
1419 if(samemother && fgEvent) {
1420 if(pair->GetFirstDaughter()->GetLabel() > 0) {
4533e78e 1421 const AliVParticle *motherMC = 0x0;
1422 if(fgEvent->IsA() == AliESDEvent::Class()) motherMC = (AliMCParticle*)mc->GetMCTrackMother((AliESDtrack*)pair->GetFirstDaughter());
1423 else if(fgEvent->IsA() == AliAODEvent::Class()) motherMC = (AliAODMCParticle*)mc->GetMCTrackMother((AliAODTrack*)pair->GetFirstDaughter());
1424 Double_t vtxX, vtxY, vtxZ;
1425 if(motherMC && mc->GetPrimaryVertex(vtxX,vtxY,vtxZ)) {
1bb1fef1 1426 Int_t motherLbl = motherMC->GetLabel();
1427 values[AliDielectronVarManager::kHasCocktailMother]=mc->CheckParticleSource(motherLbl, AliDielectronSignalMC::kDirect);
4533e78e 1428 const Double_t lxyMC = ( (motherMC->Xv() - vtxX) * motherMC->Px() +
1429 (motherMC->Yv() - vtxY) * motherMC->Py() ) / motherMC->Pt();
5720c765 1430 const Double_t pseudoMC = lxyMC * (TDatabasePDG::Instance()->GetParticle(443)->Mass())/motherMC->Pt();
1431 values[AliDielectronVarManager::kPseudoProperTimeResolution] = values[AliDielectronVarManager::kPseudoProperTime] - pseudoMC;
1432 if (errPseudoProperTime2 > 0)
1433 values[AliDielectronVarManager::kPseudoProperTimePull] = values[AliDielectronVarManager::kPseudoProperTimeResolution]/sqrt(errPseudoProperTime2);
1bb1fef1 1434 }
5720c765 1435 }
1436 }
2e02dba4 1437
696b423a 1438 values[AliDielectronVarManager::kTRDpidEffPair] = 0.;
1439 if (fgTRDpidEff[0][0]){
1440 Double_t valuesLeg1[AliDielectronVarManager::kNMaxValues];
1441 Double_t valuesLeg2[AliDielectronVarManager::kNMaxValues];
1442 AliVParticle* leg1 = pair->GetFirstDaughter();
1443 AliVParticle* leg2 = pair->GetSecondDaughter();
1444 if (leg1 && leg2){
1445 Fill(leg1, valuesLeg1);
1446 Fill(leg2, valuesLeg2);
1447 values[AliDielectronVarManager::kTRDpidEffPair] = valuesLeg1[AliDielectronVarManager::kTRDpidEffLeg]*valuesLeg2[AliDielectronVarManager::kTRDpidEffLeg];
1448 }
1449 }
b2ad74d0 1450
1451
8df8e382 1452 }//if (mc->HasMC())
1453
696b423a 1454 AliVParticle* leg1 = pair->GetFirstDaughter();
1455 AliVParticle* leg2 = pair->GetSecondDaughter();
1456 if (leg1)
1457 values[AliDielectronVarManager::kMomAsymDau1] = (values[AliDielectronVarManager::kP] != 0)? leg1->P() / values[AliDielectronVarManager::kP]: 0;
1458 else
1459 values[AliDielectronVarManager::kMomAsymDau1] = -9999.;
1460 if (leg2)
1461 values[AliDielectronVarManager::kMomAsymDau2] = (values[AliDielectronVarManager::kP] != 0)? leg2->P() / values[AliDielectronVarManager::kP]: 0;
1462 else
1463 values[AliDielectronVarManager::kMomAsymDau2] = -9999.;
b2a297fa 1464}
1465
1201a1a9 1466inline void AliDielectronVarManager::FillVarKFParticle(const AliKFParticle *particle, Double_t * const values)
1467{
1468 //
1469 // Fill track information available in AliVParticle into an array
1470 //
1471 values[AliDielectronVarManager::kPx] = particle->GetPx();
1472 values[AliDielectronVarManager::kPy] = particle->GetPy();
1473 values[AliDielectronVarManager::kPz] = particle->GetPz();
1474 values[AliDielectronVarManager::kPt] = particle->GetPt();
1475 values[AliDielectronVarManager::kP] = particle->GetP();
1476
1477 values[AliDielectronVarManager::kXv] = particle->GetX();
1478 values[AliDielectronVarManager::kYv] = particle->GetY();
1479 values[AliDielectronVarManager::kZv] = particle->GetZ();
1480
1481 values[AliDielectronVarManager::kOneOverPt] = 0;
1482 values[AliDielectronVarManager::kPhi] = particle->GetPhi();
1483 values[AliDielectronVarManager::kTheta] = 0.;
1484 values[AliDielectronVarManager::kEta] = particle->GetEta();
1485 values[AliDielectronVarManager::kY] = ((particle->GetE()*particle->GetE()-particle->GetPx()*particle->GetPx()-particle->GetPy()*particle->GetPy()-particle->GetPz()*particle->GetPz())>0.) ? TLorentzVector(particle->GetPx(),particle->GetPy(),particle->GetPz(),particle->GetE()).Rapidity() : -1111.;
1486
1487 values[AliDielectronVarManager::kE] = particle->GetE();
1488 values[AliDielectronVarManager::kM] = particle->GetMass();
1489 values[AliDielectronVarManager::kCharge] = particle->GetQ();
1490
1491 values[AliDielectronVarManager::kNclsITS] = 0;
5720c765 1492 values[AliDielectronVarManager::kITSchi2Cl] = -1;
1201a1a9 1493 values[AliDielectronVarManager::kNclsTPC] = 0;
5720c765 1494 values[AliDielectronVarManager::kNclsSTPC] = 0;
1495 values[AliDielectronVarManager::kNclsSFracTPC] = 0;
1201a1a9 1496 values[AliDielectronVarManager::kNclsTPCiter1] = 0;
1497 values[AliDielectronVarManager::kNFclsTPC] = 0;
1498 values[AliDielectronVarManager::kNclsTRD] = 0;
1499 values[AliDielectronVarManager::kTRDntracklets] = 0;
1500 values[AliDielectronVarManager::kTRDpidQuality] = 0;
1501 values[AliDielectronVarManager::kTPCchi2Cl] = 0;
1502 values[AliDielectronVarManager::kTrackStatus] = 0;
07480df6 1503 values[AliDielectronVarManager::kFilterBit] = 0;
1201a1a9 1504 values[AliDielectronVarManager::kTRDprobEle] = 0;
1505 values[AliDielectronVarManager::kTRDprobPio] = 0;
1506 values[AliDielectronVarManager::kTPCsignalN] = 0;
5720c765 1507 values[AliDielectronVarManager::kTPCclsDiff] = 0;
ba15fdfb 1508 values[AliDielectronVarManager::kTPCsignalNfrac]= 0;
1201a1a9 1509 values[AliDielectronVarManager::kImpactParXY] = 0;
1510 values[AliDielectronVarManager::kImpactParZ] = 0;
1511 values[AliDielectronVarManager::kPIn] = 0;
5720c765 1512 values[AliDielectronVarManager::kYsignedIn] = 0;
1201a1a9 1513 values[AliDielectronVarManager::kTPCsignal] = 0;
ba15fdfb 1514 values[AliDielectronVarManager::kTOFsignal] = 0;
1515 values[AliDielectronVarManager::kTOFbeta] = 0;
1201a1a9 1516 values[AliDielectronVarManager::kTPCnSigmaEle] = 0;
1517 values[AliDielectronVarManager::kTPCnSigmaPio] = 0;
1518 values[AliDielectronVarManager::kTPCnSigmaMuo] = 0;
1519 values[AliDielectronVarManager::kTPCnSigmaKao] = 0;
1520 values[AliDielectronVarManager::kTPCnSigmaPro] = 0;
1521 values[AliDielectronVarManager::kITSclusterMap] = 0;
1522
ba15fdfb 1523 values[AliDielectronVarManager::kPdgCode] = -1;
1524 values[AliDielectronVarManager::kPdgCodeMother] = -1;
1525 values[AliDielectronVarManager::kPdgCodeGrandMother] = -1;
1bb1fef1 1526 values[AliDielectronVarManager::kHasCocktailMother]=0;
1527 values[AliDielectronVarManager::kHasCocktailGrandMother]=0;
1201a1a9 1528
5720c765 1529// if ( fgEvent ) AliDielectronVarManager::Fill(fgEvent, values);
1530 for (Int_t i=AliDielectronVarManager::kPairMax; i<AliDielectronVarManager::kNMaxValues; ++i)
1531 values[i]=fgData[i];
1532
1201a1a9 1533}
6551594b 1534
1535inline void AliDielectronVarManager::FillVarVEvent(const AliVEvent *event, Double_t * const values)
1536{
1537 //
1538 // Fill event information available for histogramming into an array
1539 //
5720c765 1540 values[AliDielectronVarManager::kRunNumber] = event->GetRunNumber();
a823f01b 1541 if(fgCurrentRun!=event->GetRunNumber()) {
1542 if(fgVZEROCalibrationFile.Contains(".root")) InitVZEROCalibrationHistograms(event->GetRunNumber());
1543 if(fgVZERORecenteringFile.Contains(".root")) InitVZERORecenteringHistograms(event->GetRunNumber());
1544 fgCurrentRun=event->GetRunNumber();
1545 }
0c09cae4 1546 values[AliDielectronVarManager::kMixingBin]=0;
1547
6551594b 1548 const AliVVertex *primVtx = event->GetPrimaryVertex();
8df8e382 1549
1550 values[AliDielectronVarManager::kXvPrim] = 0;
1551 values[AliDielectronVarManager::kYvPrim] = 0;
1552 values[AliDielectronVarManager::kZvPrim] = 0;
0c09cae4 1553 values[AliDielectronVarManager::kNVtxContrib] = 0;
ba15fdfb 1554// values[AliDielectronVarManager::kChi2NDF] = 0; //This is the pair value!!!
5720c765 1555
0c09cae4 1556 values[AliDielectronVarManager::kNTrk] = 0;
1557 values[AliDielectronVarManager::kNVtxContrib] = 0;
1558 values[AliDielectronVarManager::kNacc] = 0;
1559 values[AliDielectronVarManager::kNaccTrcklts] = 0;
5720c765 1560 values[AliDielectronVarManager::kNaccTrcklts0916] = 0;
0c09cae4 1561 values[AliDielectronVarManager::kNevents] = 0; //always fill bin 0;
7a1ddc4b 1562 values[AliDielectronVarManager::kRefMult] = 0;
1563 values[AliDielectronVarManager::kRefMultTPConly] = 0;
8df8e382 1564
0c09cae4 1565 if (primVtx){
1566 values[AliDielectronVarManager::kXvPrim] = primVtx->GetX();
1567 values[AliDielectronVarManager::kYvPrim] = primVtx->GetY();
1568 values[AliDielectronVarManager::kZvPrim] = primVtx->GetZ();
1569 values[AliDielectronVarManager::kNVtxContrib] = primVtx->GetNContributors();
1570 }
1571 // values[AliDielectronVarManager::kChi2NDF] = primVtx->GetChi2perNDF(); //this is the pair value
7a1ddc4b 1572
5720c765 1573 values[AliDielectronVarManager::kNTrk] = event->GetNumberOfTracks();
1574 values[AliDielectronVarManager::kNacc] = AliDielectronHelper::GetNacc(event);
41a64bfb 1575 values[AliDielectronVarManager::kMatchEffITSTPC] = AliDielectronHelper::GetITSTPCMatchEff(event);
5720c765 1576 values[AliDielectronVarManager::kNaccTrcklts] = AliDielectronHelper::GetNaccTrcklts(event); // etaRange = 1.6 (default)
1577 values[AliDielectronVarManager::kNaccTrcklts0916] = AliDielectronHelper::GetNaccTrcklts(event,1.6)-AliDielectronHelper::GetNaccTrcklts(event,.9);
1578 // values[AliDielectronVarManager::kNaccTrcklts05] = AliDielectronHelper::GetNaccTrcklts(event, 0.5);
1579 // values[AliDielectronVarManager::kNaccTrcklts10] = AliDielectronHelper::GetNaccTrcklts(event, 1.0);
1580 // values[AliDielectronVarManager::kNaccTrckltsCorr] = AliDielectronHelper::GetNaccTrckltsCorrected(event, values[AliDielectronVarManager::kNaccTrcklts], values[AliDielectronVarManager::kZvPrim]);
1581
20b6b7d4 1582 values[AliDielectronVarManager::kPhiMaxPt] = 0;
1583 values[AliDielectronVarManager::kMaxPt] = 0;
1584 Double_t ptMaxEv = -1., phiptMaxEv= -1.;
1585 for(Int_t itrk=0; itrk<event->GetNumberOfTracks(); itrk++) {
1586 AliVParticle *part= event->GetTrack(itrk);
1587 if(part->Pt() > ptMaxEv) {
1588 ptMaxEv = part->Pt();
1589 phiptMaxEv = part->Phi();
1590 }
1591 }
1592 values[AliDielectronVarManager::kPhiMaxPt] = phiptMaxEv;
1593 values[AliDielectronVarManager::kMaxPt] = ptMaxEv;
1594
5720c765 1595
07480df6 1596 // event plane quantities from the AliEPSelectionTask
882407d7 1597 for(Int_t ivar=AliDielectronVarManager::kv0ArpH2; ivar<=kv0C0v0C3DiffH2; ivar++) values[ivar] = 0.0; // v0 variables
07480df6 1598 for(Int_t ivar=AliDielectronVarManager::kTPCxH2; ivar<=kTPCsub12DiffH2uc; ivar++) values[ivar] = 0.0; // tpc variables
8d42b5b9 1599
1600 // ep angle interval [todo, fill]
99345a64 1601 AliEventplane *ep = const_cast<AliVEvent*>(event)->GetEventplane();
07480df6 1602 if(ep) {
8d42b5b9 1603
1604 // TPC event plane quantities (uncorrected)
07480df6 1605 TVector2 *qstd = ep->GetQVector(); // This is the "standard" Q-Vector for TPC
1606 TVector2 *qsub1 = ep->GetQsub1(); // random subevent plane
1607 TVector2 *qsub2 = ep->GetQsub2();
1608 if(qstd && qsub1 && qsub2) {
8d42b5b9 1609 values[AliDielectronVarManager::kTPCxH2uc] = qstd->X();
1610 values[AliDielectronVarManager::kTPCyH2uc] = qstd->Y();
1611 values[AliDielectronVarManager::kTPCmagH2uc] = qstd->Mod();
1612 values[AliDielectronVarManager::kTPCrpH2uc] = ((TMath::Abs(qstd->X())>1.0e-10) ? TMath::ATan2(qstd->Y(),qstd->X())/2.0 : 0.0);
07480df6 1613 values[AliDielectronVarManager::kTPCsub1xH2uc] = qsub1->X();
1614 values[AliDielectronVarManager::kTPCsub1yH2uc] = qsub1->Y();
5c2b3da2 1615 values[AliDielectronVarManager::kTPCsub1rpH2uc] = ((TMath::Abs(qsub1->X())>1.0e-10) ? TMath::ATan2(qsub1->Y(),qsub1->X())/2.0 : 0.0);
07480df6 1616 values[AliDielectronVarManager::kTPCsub2xH2uc] = qsub2->X();
1617 values[AliDielectronVarManager::kTPCsub2yH2uc] = qsub2->Y();
5c2b3da2 1618 values[AliDielectronVarManager::kTPCsub2rpH2uc] = ((TMath::Abs(qsub2->X())>1.0e-10) ? TMath::ATan2(qsub2->Y(),qsub2->X())/2.0 : 0.0);
07480df6 1619
1620 values[AliDielectronVarManager::kTPCsub12DiffH2uc] = TMath::Cos( 2.*(values[AliDielectronVarManager::kTPCsub1rpH2uc] -
1621 values[AliDielectronVarManager::kTPCsub2rpH2uc]) );
1622 }
20b6b7d4 1623
8d42b5b9 1624 // VZERO event plane
1625 TVector2 qvec;
1626 Double_t qx = 0, qy = 0;
1627 ep->CalculateVZEROEventPlane(event,10, 2, qx, qy); qvec.Set(qx,qy);
1628 values[AliDielectronVarManager::kv0ACrpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1629 values[AliDielectronVarManager::kv0ACxH2] = qvec.X();
1630 values[AliDielectronVarManager::kv0ACyH2] = qvec.Y();
1631 values[AliDielectronVarManager::kv0ACmagH2] = qvec.Mod();
1632 ep->CalculateVZEROEventPlane(event, 8, 2, qx, qy); qvec.Set(qx,qy);
1633 values[AliDielectronVarManager::kv0ArpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1634 values[AliDielectronVarManager::kv0AxH2] = qvec.X();
1635 values[AliDielectronVarManager::kv0AyH2] = qvec.Y();
1636 values[AliDielectronVarManager::kv0AmagH2] = qvec.Mod();
1637 ep->CalculateVZEROEventPlane(event, 9, 2, qx, qy); qvec.Set(qx,qy);
1638 values[AliDielectronVarManager::kv0CrpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1639 values[AliDielectronVarManager::kv0CxH2] = qvec.X();
1640 values[AliDielectronVarManager::kv0CyH2] = qvec.Y();
1641 values[AliDielectronVarManager::kv0CmagH2] = qvec.Mod();
1642 ep->CalculateVZEROEventPlane(event, 0, 0, 2, qx, qy); qvec.Set(qx,qy);
1643 values[AliDielectronVarManager::kv0C0rpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1644 ep->CalculateVZEROEventPlane(event, 3, 3, 2, qx, qy); qvec.Set(qx,qy);
1645 values[AliDielectronVarManager::kv0C3rpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1646 ep->CalculateVZEROEventPlane(event, 4, 4, 2, qx, qy); qvec.Set(qx,qy);
1647 values[AliDielectronVarManager::kv0A0rpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1648 ep->CalculateVZEROEventPlane(event, 7, 7, 2, qx, qy); qvec.Set(qx,qy);
1649 values[AliDielectronVarManager::kv0A3rpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
07480df6 1650 } //if: eventplane
1651
5720c765 1652 // ESD VZERO information
1653 AliVVZERO* vzeroData = event->GetVZEROData();
1654 values[AliDielectronVarManager::kMultV0A] = 0.0;
1655 values[AliDielectronVarManager::kMultV0C] = 0.0;
1656 values[AliDielectronVarManager::kAdcV0A] = 0.0;
1657 values[AliDielectronVarManager::kAdcV0C] = 0.0;
1658 for(Int_t i=0; i<32; ++i) {
a823f01b 1659 values[AliDielectronVarManager::kVZEROchMult+i] = vzeroData->GetMultiplicity(i);
1660 values[AliDielectronVarManager::kVZEROchMult+32+i] = vzeroData->GetMultiplicity(i+32);
1661 //values[AliDielectronVarManager::kVZEROchMult+i] = event->GetVZEROEqMultiplicity(i);
1662 //values[AliDielectronVarManager::kVZEROchMult+32+i] = event->GetVZEROEqMultiplicity(i+32);
5720c765 1663 values[AliDielectronVarManager::kMultV0A] += vzeroData->GetMultiplicityV0A(i);
1664 values[AliDielectronVarManager::kMultV0C] += vzeroData->GetMultiplicityV0C(i);
1665 //values[AliDielectronVarManager::kAdcV0A] += vzeroData->GetAdcV0A(i);
1666 //values[AliDielectronVarManager::kAdcV0C] += vzeroData->GetAdcV0C(i);
1667 }
1668 values[AliDielectronVarManager::kMultV0] = values[AliDielectronVarManager::kMultV0A] + values[AliDielectronVarManager::kMultV0C];
1669 values[AliDielectronVarManager::kAdcV0] = values[AliDielectronVarManager::kAdcV0A] + values[AliDielectronVarManager::kAdcV0C];
1670 // VZERO event plane quantities
1671 Double_t qvec[3]={0.0};
1672 GetVzeroRP(event, qvec,0); // V0-A
1673 values[AliDielectronVarManager::kV0AxH2] = qvec[0]; values[AliDielectronVarManager::kV0AyH2] = qvec[1];
1674 values[AliDielectronVarManager::kV0ArpH2] = qvec[2];
1675 qvec[0]=0.0; qvec[1]=0.0; qvec[2]=0.0;
1676 GetVzeroRP(event, qvec,1); // V0-C
1677 values[AliDielectronVarManager::kV0CxH2] = qvec[0]; values[AliDielectronVarManager::kV0CyH2] = qvec[1];
1678 values[AliDielectronVarManager::kV0CrpH2] = qvec[2];
1679 qvec[0]=0.0; qvec[1]=0.0; qvec[2]=0.0;
1680 GetVzeroRP(event, qvec,2); // V0-A and V0-C combined
1681 values[AliDielectronVarManager::kV0ACxH2] = qvec[0]; values[AliDielectronVarManager::kV0ACyH2] = qvec[1];
1682 values[AliDielectronVarManager::kV0ACrpH2] = qvec[2];
1683 // VZERO event plane resolution
1684 values[AliDielectronVarManager::kV0ArpResH2] = 1.0;
1685 values[AliDielectronVarManager::kV0CrpResH2] = 1.0;
1686 values[AliDielectronVarManager::kV0ACrpResH2] = 1.0;
1687 // Q vector components correlations
1688 values[AliDielectronVarManager::kV0XaXcH2] = values[AliDielectronVarManager::kV0AxH2]*values[AliDielectronVarManager::kV0CxH2];
1689 values[AliDielectronVarManager::kV0XaYaH2] = values[AliDielectronVarManager::kV0AxH2]*values[AliDielectronVarManager::kV0AyH2];
1690 values[AliDielectronVarManager::kV0XaYcH2] = values[AliDielectronVarManager::kV0AxH2]*values[AliDielectronVarManager::kV0CyH2];
1691 values[AliDielectronVarManager::kV0YaXcH2] = values[AliDielectronVarManager::kV0AyH2]*values[AliDielectronVarManager::kV0CxH2];
1692 values[AliDielectronVarManager::kV0YaYcH2] = values[AliDielectronVarManager::kV0AyH2]*values[AliDielectronVarManager::kV0CyH2];
1693 values[AliDielectronVarManager::kV0XcYcH2] = values[AliDielectronVarManager::kV0CxH2]*values[AliDielectronVarManager::kV0CyH2];
1694
07480df6 1695
1696 // event plane differences used for EP resolution calculation
5720c765 1697 values[AliDielectronVarManager::kV0ATPCDiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kV0ArpH2] -
1698 values[AliDielectronVarManager::kTPCrpH2]) );
1699
1700 values[AliDielectronVarManager::kV0CTPCDiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kV0CrpH2] -
1701 values[AliDielectronVarManager::kTPCrpH2]) );
1702
1703 values[AliDielectronVarManager::kV0AV0CDiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kV0ArpH2] -
1704 values[AliDielectronVarManager::kV0CrpH2]) );
1705
1706 values[AliDielectronVarManager::kv0ATPCDiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0ArpH2] -
1707 values[AliDielectronVarManager::kTPCrpH2]) );
1708
1709 values[AliDielectronVarManager::kv0CTPCDiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0CrpH2] -
1710 values[AliDielectronVarManager::kTPCrpH2]) );
1711
1712 values[AliDielectronVarManager::kv0Av0CDiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0ArpH2] -
1713 values[AliDielectronVarManager::kv0CrpH2]) );
d327d9cd 1714
1d2d723c 1715 values[AliDielectronVarManager::kv0Av0C0DiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0ArpH2] -
1716 values[AliDielectronVarManager::kv0C0rpH2]) );
1717
1718 values[AliDielectronVarManager::kv0Av0C3DiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0ArpH2] -
1719 values[AliDielectronVarManager::kv0C3rpH2]) );
1720
1721 values[AliDielectronVarManager::kv0Cv0A0DiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0CrpH2] -
1722 values[AliDielectronVarManager::kv0A0rpH2]) );
1723
1724 values[AliDielectronVarManager::kv0Cv0A3DiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0CrpH2] -
1725 values[AliDielectronVarManager::kv0A3rpH2]) );
1726
1727 values[AliDielectronVarManager::kv0A0v0A3DiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0A0rpH2] -
1728 values[AliDielectronVarManager::kv0A3rpH2]) );
1729
1730 values[AliDielectronVarManager::kv0C0v0C3DiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kv0C0rpH2] -
1731 values[AliDielectronVarManager::kv0C3rpH2]) );
1732
1733
6551594b 1734}
1735
1736inline void AliDielectronVarManager::FillVarESDEvent(const AliESDEvent *event, Double_t * const values)
1737{
1738 //
1739 // Fill event information available for histogramming into an array
1740 //
1741
1742 // Fill common AliVEvent interface information
1743 FillVarVEvent(event, values);
ffbede40 1744
a823f01b 1745 Double_t centralityF=-1; Double_t centralitySPD=-1;
4649c75b 1746 AliCentrality *esdCentrality = const_cast<AliESDEvent*>(event)->GetCentrality();
ffbede40 1747 if (esdCentrality) centralityF = esdCentrality->GetCentralityPercentile("V0M");
a823f01b 1748 if (esdCentrality) centralitySPD = esdCentrality->GetCentralityPercentile("CL1");
ffbede40 1749
6551594b 1750 // Fill AliESDEvent interface specific information
1751 const AliESDVertex *primVtx = event->GetPrimaryVertex();
1752 values[AliDielectronVarManager::kXRes] = primVtx->GetXRes();
1753 values[AliDielectronVarManager::kYRes] = primVtx->GetYRes();
1754 values[AliDielectronVarManager::kZRes] = primVtx->GetZRes();
ffbede40 1755 values[AliDielectronVarManager::kCentrality] = centralityF;
a823f01b 1756 values[AliDielectronVarManager::kCentralitySPD] = centralitySPD;
6bd73568 1757
1758 const AliESDVertex *vtxTPC = event->GetPrimaryVertexTPC();
1759 values[AliDielectronVarManager::kNVtxContribTPC] = (vtxTPC ? vtxTPC->GetNContributors() : 0);
1760
5720c765 1761 // Event multiplicity estimators
1762 Int_t nTrSPD05=0; Int_t nTrITSTPC05=0; Int_t nTrITSSA05=0;
1763 event->EstimateMultiplicity(nTrSPD05, nTrITSTPC05, nTrITSSA05, 0.5);
1764 values[AliDielectronVarManager::kNaccTrckltsEsd05] = nTrSPD05;
1765 values[AliDielectronVarManager::kNaccItsTpcEsd05] = nTrITSTPC05;
1766 values[AliDielectronVarManager::kNaccItsPureEsd05] = nTrITSSA05;
1767 values[AliDielectronVarManager::kNaccTrckltsEsd05Corr] =
1768 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrSPD05),values[AliDielectronVarManager::kZvPrim],0);
1769 values[AliDielectronVarManager::kNaccItsTpcEsd05Corr] =
1770 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrITSTPC05),values[AliDielectronVarManager::kZvPrim],3);
1771 values[AliDielectronVarManager::kNaccItsPureEsd05Corr] =
1772 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrITSSA05),values[AliDielectronVarManager::kZvPrim],6);
1773
1774 Int_t nTrSPD10=0; Int_t nTrITSTPC10=0; Int_t nTrITSSA10=0;
1775 event->EstimateMultiplicity(nTrSPD10, nTrITSTPC10, nTrITSSA10, 1.0);
1776 values[AliDielectronVarManager::kNaccTrckltsEsd10] = nTrSPD10;
1777 values[AliDielectronVarManager::kNaccItsTpcEsd10] = nTrITSTPC10;
1778 values[AliDielectronVarManager::kNaccItsPureEsd10] = nTrITSSA10;
1779 values[AliDielectronVarManager::kNaccTrckltsEsd10Corr] =
1780 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrSPD10),values[AliDielectronVarManager::kZvPrim],1);
1781 values[AliDielectronVarManager::kNaccItsTpcEsd10Corr] =
1782 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrITSTPC10),values[AliDielectronVarManager::kZvPrim],4);
1783 values[AliDielectronVarManager::kNaccItsPureEsd10Corr] =
1784 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrITSSA10),values[AliDielectronVarManager::kZvPrim],7);
1785
1786 Int_t nTrSPD16=0; Int_t nTrITSTPC16=0; Int_t nTrITSSA16=0;
1787 event->EstimateMultiplicity(nTrSPD16, nTrITSTPC16, nTrITSSA16, 1.6);
1788 values[AliDielectronVarManager::kNaccTrckltsEsd16] = nTrSPD16;
1789 values[AliDielectronVarManager::kNaccItsTpcEsd16] = nTrITSTPC16;
1790 values[AliDielectronVarManager::kNaccItsPureEsd16] = nTrITSSA16;
1791 values[AliDielectronVarManager::kNaccTrckltsEsd16Corr] =
1792 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrSPD16),values[AliDielectronVarManager::kZvPrim],2);
1793 values[AliDielectronVarManager::kNaccItsTpcEsd16Corr] =
1794 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrITSTPC16),values[AliDielectronVarManager::kZvPrim],5);
1795 values[AliDielectronVarManager::kNaccItsPureEsd16Corr] =
1796 AliDielectronHelper::GetNaccTrckltsCorrected(event,Double_t(nTrITSSA16),values[AliDielectronVarManager::kZvPrim],8);
1797
6551594b 1798}
8d42b5b9 1799
6551594b 1800inline void AliDielectronVarManager::FillVarAODEvent(const AliAODEvent *event, Double_t * const values)
1801{
1802 //
1803 // Fill event information available for histogramming into an array
20b6b7d4 1804 //
6551594b 1805
1806 // Fill common AliVEvent interface information
1807 FillVarVEvent(event, values);
1808
1809 // Fill AliAODEvent interface specific information
d327d9cd 1810 AliAODHeader *header = event->GetHeader();
20b6b7d4 1811
a823f01b 1812 Double_t centralityF=-1; Double_t centralitySPD=-1;
d327d9cd 1813 AliCentrality *aodCentrality = header->GetCentralityP();
1814 if (aodCentrality) centralityF = aodCentrality->GetCentralityPercentile("V0M");
a823f01b 1815 if (aodCentrality) centralitySPD = aodCentrality->GetCentralityPercentile("CL1");
d327d9cd 1816 values[AliDielectronVarManager::kCentrality] = centralityF;
a823f01b 1817 values[AliDielectronVarManager::kCentralitySPD] = centralitySPD;
20b6b7d4 1818
7a1ddc4b 1819 values[AliDielectronVarManager::kRefMult] = header->GetRefMultiplicity(); // similar to Ntrk
1820 values[AliDielectronVarManager::kRefMultTPConly] = header->GetTPConlyRefMultiplicity(); // similar to Nacc
1821
20b6b7d4 1822 // nanoAODs (w/o AliCentrality branch) should have the VOM centrality stored in the header
1823 if(!header->GetCentralityP())
1824 values[AliDielectronVarManager::kCentrality] = header->GetCentrality();
1825 // nanoAODs (w/o AliEventPlane branch) should have the tpc event plane angle stored in the header
7a1ddc4b 1826 if(!header->GetEventplaneP()) {
61749cc7 1827
8d42b5b9 1828 // values[AliDielectronVarManager::kNTrk] = header->GetRefMultiplicity(); // overwritten datamembers in "our" nanoAODs
1829 // values[AliDielectronVarManager::kNacc] = header->GetRefMultiplicityPos(); // overwritten datamembers in "our" nanoAODs
61749cc7 1830
8d42b5b9 1831 TVector2 qvec;
123e7cdf 1832 // TPC
8d42b5b9 1833 qvec.Set(header->GetEventplaneQx(), header->GetEventplaneQy());
1834 values[AliDielectronVarManager::kTPCxH2uc] = qvec.X();
1835 values[AliDielectronVarManager::kTPCyH2uc] = qvec.Y();
1836 values[AliDielectronVarManager::kTPCmagH2uc] = qvec.Mod();
1837 values[AliDielectronVarManager::kTPCrpH2uc] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1838
1839 // VZERO
1840 AliEventplane ep2;
1841 // get event plane corrections from the VZERO EP selection task
1842 AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
1843 AliVZEROEPSelectionTask *eptask = dynamic_cast<AliVZEROEPSelectionTask *>(man->GetTask("AliVZEROEPSelectionTask"));
1844 if(eptask) eptask->SetEventplaneParams(&ep2,centralityF);
1845 else printf("no VZERO event plane selection task added! \n");
1846
1847 Double_t qx = 0, qy = 0;
1848 ep2.CalculateVZEROEventPlane(event,10, 2, qx, qy); qvec.Set(qx,qy);
1849 values[AliDielectronVarManager::kv0ACrpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1850 values[AliDielectronVarManager::kv0ACxH2] = qvec.X();
1851 values[AliDielectronVarManager::kv0ACyH2] = qvec.Y();
1852 values[AliDielectronVarManager::kv0ACmagH2] = qvec.Mod();
1853 ep2.CalculateVZEROEventPlane(event, 8, 2, qx, qy); qvec.Set(qx,qy);
1854 values[AliDielectronVarManager::kv0ArpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1855 values[AliDielectronVarManager::kv0AxH2] = qvec.X();
1856 values[AliDielectronVarManager::kv0AyH2] = qvec.Y();
1857 values[AliDielectronVarManager::kv0AmagH2] = qvec.Mod();
1858 ep2.CalculateVZEROEventPlane(event, 9, 2, qx, qy); qvec.Set(qx,qy);
1859 values[AliDielectronVarManager::kv0CrpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1860 values[AliDielectronVarManager::kv0CxH2] = qvec.X();
1861 values[AliDielectronVarManager::kv0CyH2] = qvec.Y();
1862 values[AliDielectronVarManager::kv0CmagH2] = qvec.Mod();
1863 ep2.CalculateVZEROEventPlane(event, 0, 0, 2, qx, qy); qvec.Set(qx,qy);
1864 values[AliDielectronVarManager::kv0C0rpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1865 ep2.CalculateVZEROEventPlane(event, 3, 3, 2, qx, qy); qvec.Set(qx,qy);
1866 values[AliDielectronVarManager::kv0C3rpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1867 ep2.CalculateVZEROEventPlane(event, 4, 4, 2, qx, qy); qvec.Set(qx,qy);
1868 values[AliDielectronVarManager::kv0A0rpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
1869 ep2.CalculateVZEROEventPlane(event, 7, 7, 2, qx, qy); qvec.Set(qx,qy);
1870 values[AliDielectronVarManager::kv0A3rpH2] = ((TMath::Abs(qvec.X())>1.0e-10) ? TMath::ATan2(qvec.Y(),qvec.X())/2.0 : 0.0);
7a1ddc4b 1871 }
20b6b7d4 1872
99345a64 1873 const AliAODVertex *vtxtpc = GetVertex(event, AliAODVertex::kMainTPC);
1874 values[AliDielectronVarManager::kNVtxContribTPC] = (vtxtpc ? vtxtpc->GetNContributors() : 0);
35363fb9 1875
6551594b 1876}
1877
1878inline void AliDielectronVarManager::FillVarMCEvent(const AliMCEvent *event, Double_t * const values)
35363fb9 1879{
6551594b 1880 //
1881 // Fill event information available for histogramming into an array
35363fb9 1882 //
1883
6551594b 1884 // Fill common AliVEvent interface information
5720c765 1885 // FillVarVEvent(event, values);
1886 const AliVVertex* vtx = event->GetPrimaryVertex();
1887 values[AliDielectronVarManager::kXvPrim] = (vtx ? vtx->GetX() : 0.0);
1888 values[AliDielectronVarManager::kYvPrim] = (vtx ? vtx->GetY() : 0.0);
1889 values[AliDielectronVarManager::kZvPrim] = (vtx ? vtx->GetZ() : 0.0);
6551594b 1890 // Fill AliMCEvent interface specific information
5720c765 1891 values[AliDielectronVarManager::kNch] = AliDielectronHelper::GetNch(event, 1.6);
1892 values[AliDielectronVarManager::kNch05] = AliDielectronHelper::GetNch(event, 0.5);
1893 values[AliDielectronVarManager::kNch10] = AliDielectronHelper::GetNch(event, 1.0);
1894
1895 values[AliDielectronVarManager::kNumberOfJPsis] = AliDielectronHelper::GetNMothers(event, 0.9, 443, 11);
1896 values[AliDielectronVarManager::kNumberOfJPsisPrompt] = AliDielectronHelper::GetNMothers(event, 0.9, 443, 11, 1);
1897 values[AliDielectronVarManager::kNumberOfJPsisNPrompt] = AliDielectronHelper::GetNMothers(event, 0.9, 443, 11, 0);
1898}
1899
1900inline void AliDielectronVarManager::FillVarTPCEventPlane(const AliEventplane *evplane, Double_t * const values)
1901{
1902 //
1903 // Fill TPC event plane information after correction
1904 //
07480df6 1905 if(evplane) {
1906 TVector2 *qcorr = const_cast<AliEventplane *>(evplane)->GetQVector(); // This is the "corrected" Q-Vector
1907 TVector2 *qcsub1 = const_cast<AliEventplane *>(evplane)->GetQsub1();
1908 TVector2 *qcsub2 = const_cast<AliEventplane *>(evplane)->GetQsub2();
4533e78e 1909 if(qcorr) {
07480df6 1910 values[AliDielectronVarManager::kTPCxH2] = qcorr->X();
1911 values[AliDielectronVarManager::kTPCyH2] = qcorr->Y();
6bd73568 1912 values[AliDielectronVarManager::kTPCmagH2] = qcorr->Mod();
5c2b3da2 1913 values[AliDielectronVarManager::kTPCrpH2] = ((TMath::Abs(qcorr->X())>1.0e-10) ? TMath::ATan2(qcorr->Y(),qcorr->X())/2.0 : 0.0);
4533e78e 1914 }
1915 if(qcsub1 && qcsub2) {
07480df6 1916 values[AliDielectronVarManager::kTPCsub1xH2] = qcsub1->X();
1917 values[AliDielectronVarManager::kTPCsub1yH2] = qcsub1->Y();
5c2b3da2 1918 values[AliDielectronVarManager::kTPCsub1rpH2] = ((TMath::Abs(qcsub1->X())>1.0e-10) ? TMath::ATan2(qcsub1->Y(),qcsub1->X())/2.0 : 0.0);
07480df6 1919
1920 values[AliDielectronVarManager::kTPCsub2xH2] = qcsub2->X();
1921 values[AliDielectronVarManager::kTPCsub2yH2] = qcsub2->Y();
5c2b3da2 1922 values[AliDielectronVarManager::kTPCsub2rpH2] = ((TMath::Abs(qcsub2->X())>1.0e-10) ? TMath::ATan2(qcsub2->Y(),qcsub2->X())/2.0 : 0.0);
07480df6 1923
1924 values[AliDielectronVarManager::kTPCsub12DiffH2] = TMath::Cos( 2.*(values[AliDielectronVarManager::kTPCsub1rpH2] -
1925 values[AliDielectronVarManager::kTPCsub2rpH2]) );
1926 values[AliDielectronVarManager::kTPCsub12DiffH2Sin] = TMath::Sin( 2.*(values[AliDielectronVarManager::kTPCsub1rpH2] -
1927 values[AliDielectronVarManager::kTPCsub2rpH2]) );
1928 }
1929 }
5720c765 1930}
a823f01b 1931
6551594b 1932inline void AliDielectronVarManager::InitESDpid(Int_t type)
1933{
1934 //
1935 // initialize PID parameters
1936 // type=0 is simulation
1937 // type=1 is data
1938
ba15fdfb 1939 if (!fgPIDResponse) fgPIDResponse=new AliESDpid((Bool_t)(type==0));
6551594b 1940 Double_t alephParameters[5];
1941 // simulation
1942 alephParameters[0] = 2.15898e+00/50.;
1943 alephParameters[1] = 1.75295e+01;
1944 alephParameters[2] = 3.40030e-09;
1945 alephParameters[3] = 1.96178e+00;
1946 alephParameters[4] = 3.91720e+00;
ba15fdfb 1947 fgPIDResponse->GetTOFResponse().SetTimeResolution(80.);
6551594b 1948
1949 // data
83e37742 1950 if (type==1){
1951 alephParameters[0] = 0.0283086/0.97;
6551594b 1952 alephParameters[1] = 2.63394e+01;
1953 alephParameters[2] = 5.04114e-11;
1954 alephParameters[3] = 2.12543e+00;
1955 alephParameters[4] = 4.88663e+00;
ba15fdfb 1956 fgPIDResponse->GetTOFResponse().SetTimeResolution(130.);
1957 fgPIDResponse->GetTPCResponse().SetMip(50.);
6551594b 1958 }
1959
ba15fdfb 1960 fgPIDResponse->GetTPCResponse().SetBetheBlochParameters(
6551594b 1961 alephParameters[0],alephParameters[1],alephParameters[2],
1962 alephParameters[3],alephParameters[4]);
1963
ba15fdfb 1964 fgPIDResponse->GetTPCResponse().SetSigma(3.79301e-03, 2.21280e+04);
6551594b 1965}
1966
2a14a7b1 1967inline void AliDielectronVarManager::InitAODpidUtil(Int_t type)
1968{
ba15fdfb 1969 if (!fgPIDResponse) fgPIDResponse=new AliAODpidUtil;
2a14a7b1 1970 Double_t alephParameters[5];
1971 // simulation
1972 alephParameters[0] = 2.15898e+00/50.;
1973 alephParameters[1] = 1.75295e+01;
1974 alephParameters[2] = 3.40030e-09;
1975 alephParameters[3] = 1.96178e+00;
1976 alephParameters[4] = 3.91720e+00;
ba15fdfb 1977 fgPIDResponse->GetTOFResponse().SetTimeResolution(80.);
2a14a7b1 1978
1979 // data
1980 if (type==1){
1981 alephParameters[0] = 0.0283086/0.97;
1982 alephParameters[1] = 2.63394e+01;
1983 alephParameters[2] = 5.04114e-11;
1984 alephParameters[3] = 2.12543e+00;
1985 alephParameters[4] = 4.88663e+00;
ba15fdfb 1986 fgPIDResponse->GetTOFResponse().SetTimeResolution(130.);
1987 fgPIDResponse->GetTPCResponse().SetMip(50.);
2a14a7b1 1988 }
1989
ba15fdfb 1990 fgPIDResponse->GetTPCResponse().SetBetheBlochParameters(
2a14a7b1 1991 alephParameters[0],alephParameters[1],alephParameters[2],
1992 alephParameters[3],alephParameters[4]);
1993
ba15fdfb 1994 fgPIDResponse->GetTPCResponse().SetSigma(3.79301e-03, 2.21280e+04);
2a14a7b1 1995}
1996
572b0139 1997
5720c765 1998inline void AliDielectronVarManager::InitEstimatorAvg(const Char_t* filename)
1999{
2000 //
2001 // initialize the profile histograms neccessary for the correction of the multiplicity estimators in pp collisions
2002 //
2003
2004 const Char_t* estimatorNames[9] = {"SPDmult05","SPDmult10","SPDmult16",
2005 "ITSTPC05", "ITSTPC10", "ITSTPC16",
2006 "ITSSA05", "ITSSA10", "ITSSA16"};
2007 const Char_t* periodNames[4] = {"LHC10b", "LHC10c", "LHC10d", "LHC10e"};
2008 TFile* file=TFile::Open(filename);
2009 if(!file) return;
2010
2011 for(Int_t ip=0; ip<4; ++ip) {
2012 for(Int_t ie=0; ie<9; ++ie) {
2013 fgMultEstimatorAvg[ip][ie] = (TProfile*)(file->Get(Form("%s_%s",estimatorNames[ie],periodNames[ip]))->Clone(Form("%s_%s_clone",estimatorNames[ie],periodNames[ip])));
2014 }
2015 }
2016}
2017
2018
2019inline void AliDielectronVarManager::InitTRDpidEffHistograms(const Char_t* filename)
2020{
2021 //
2022 // initialize the 3D histograms with the TRD pid efficiency histograms
2023 //
2024
2025 // reset the centrality ranges and the efficiency histograms
2026 for(Int_t i=0; i<10; ++i) { // centrality ranges
2027 for(Int_t j=0; j<4; ++j) fgTRDpidEffCentRanges[i][j] = -1.;
2028 if(fgTRDpidEff[i][0]) {
2029 delete fgTRDpidEff[i][0];
2030 fgTRDpidEff[i][0] = 0x0;
2031 }
2032 if(fgTRDpidEff[i][1]) {
2033 delete fgTRDpidEff[i][1];
2034 fgTRDpidEff[i][1] = 0x0;
2035 }
2036 }
2037
2038 TFile* file=TFile::Open(filename);
2039 TList* keys=file->GetListOfKeys();
2040 Int_t idxp=0; Int_t idxn=0;
2041 for(Int_t i=0; i<keys->GetEntries(); ++i) {
2042 if(idxp>=10) continue;
2043 if(idxn>=10) continue;
2044 TString name=((TKey*)keys->At(i))->ReadObj()->GetName();
2045 // Name of histograms should be in the format:
2046 // TRDeff<field>_cent_<centLow>_<centHigh>
2047 // <field> is either "BPLUS" or "BMINUS"
2048 if(!(name.Contains("BPLUS") || name.Contains("BMINUS"))) continue;
2049 TObjArray* arr = name.Tokenize("_");
2050 Bool_t isBplus = kTRUE;
2051 if(name.Contains("BMINUS")) isBplus = kFALSE;
99345a64 2052 TString centMinStr = arr->At(2)->GetName();
09d5920f 2053 TString centMaxStr = arr->At(3)->GetName();
2054 delete arr;
5720c765 2055 if(isBplus) {
2056 fgTRDpidEffCentRanges[idxp][2] = centMinStr.Atof();
2057 fgTRDpidEffCentRanges[idxp][3] = centMaxStr.Atof();
2058 fgTRDpidEff[idxp][1] = (TH3D*)(file->Get(name.Data())->Clone(Form("%s_clone",name.Data())));
2059 ++idxp;
2060 }
2061 else {
2062 fgTRDpidEffCentRanges[idxn][0] = centMinStr.Atof();
2063 fgTRDpidEffCentRanges[idxn][1] = centMaxStr.Atof();
2064 fgTRDpidEff[idxn][0] = (TH3D*)(file->Get(name.Data())->Clone(Form("%s_clone",name.Data())));
2065 ++idxn;
2066 }
2067 }
2068}
2069
2070
a823f01b 2071inline void AliDielectronVarManager::InitVZEROCalibrationHistograms(Int_t runNo) {
2072 //
2073 // Initialize the VZERO channel-by-channel calibration histograms
2074 //
2075
2076 //initialize only once
2077 if(fgVZEROCalib[0]) return;
2078
2079 for(Int_t i=0; i<64; ++i)
2080 if(fgVZEROCalib[i]) {
2081 delete fgVZEROCalib[i];
2082 fgVZEROCalib[i] = 0x0;
2083 }
2084
2085 TFile file(fgVZEROCalibrationFile.Data());
2086
2087 for(Int_t i=0; i<64; ++i){
2088 fgVZEROCalib[i] = (TProfile2D*)(file.Get(Form("RUN%d_ch%d_VtxCent", runNo, i)));
2089 if (fgVZEROCalib[i]) fgVZEROCalib[i]->SetDirectory(0x0);
2090 }
2091}
2092
2093
2094inline void AliDielectronVarManager::InitVZERORecenteringHistograms(Int_t runNo) {
2095 //
2096 // Initialize the VZERO event plane recentering histograms
2097 //
2098
2099 //initialize only once
2100 if(fgVZERORecentering[0][0]) return;
2101
2102 for(Int_t i=0; i<2; ++i)
2103 for(Int_t j=0; j<2; ++j)
2104 if(fgVZERORecentering[i][j]) {
2105 delete fgVZERORecentering[i][j];
2106 fgVZERORecentering[i][j] = 0x0;
2107 }
2108
2109 TFile file(fgVZERORecenteringFile.Data());
2110 if (!file.IsOpen()) return;
2111
2112 fgVZERORecentering[0][0] = (TProfile2D*)(file.Get(Form("RUN%d_QxA_CentVtx", runNo)));
2113 fgVZERORecentering[0][1] = (TProfile2D*)(file.Get(Form("RUN%d_QyA_CentVtx", runNo)));
2114 fgVZERORecentering[1][0] = (TProfile2D*)(file.Get(Form("RUN%d_QxC_CentVtx", runNo)));
2115 fgVZERORecentering[1][1] = (TProfile2D*)(file.Get(Form("RUN%d_QyC_CentVtx", runNo)));
2116
2117 if (fgVZERORecentering[0][0]) fgVZERORecentering[0][0]->SetDirectory(0x0);
2118 if (fgVZERORecentering[0][1]) fgVZERORecentering[0][1]->SetDirectory(0x0);
2119 if (fgVZERORecentering[1][0]) fgVZERORecentering[1][0]->SetDirectory(0x0);
2120 if (fgVZERORecentering[1][1]) fgVZERORecentering[1][1]->SetDirectory(0x0);
2121
2122}
2123
2124
5720c765 2125inline Double_t AliDielectronVarManager::GetTRDpidEfficiency(Int_t runNo, Double_t centrality,
2126 Double_t eta, Double_t trdPhi, Double_t pout,
2127 Double_t& effErr) {
2128 //
2129 // return the efficiency in the given phase space cell
2130 //
2131 // LHC10h data----------------------------------------------
2132 Bool_t isBplus = kTRUE;
2133 if(runNo<=138275) isBplus = kFALSE;
2134 // TODO: check magnetic polarity for runs in 2011 data
2135 // ---------------------------------------------------------
2136 Int_t centIdx = -1;
2137 for(Int_t icent=0; icent<10; ++icent) {
2138 if(isBplus) {
2139 if(centrality>=fgTRDpidEffCentRanges[icent][2] && centrality<fgTRDpidEffCentRanges[icent][3]) {
2140 centIdx = icent;
2141 break;
2142 }
2143 }
2144 else {
2145 if(centrality>=fgTRDpidEffCentRanges[icent][0] && centrality<fgTRDpidEffCentRanges[icent][1]) {
2146 centIdx = icent;
2147 break;
2148 }
2149 }
2150 }
ac390e40 2151 //TODO: chek logick
2152 if (centIdx<0) return 1;
2153
5720c765 2154 TH3D* effH = fgTRDpidEff[centIdx][(isBplus ? 1 : 0)];
2155 if(!effH) {effErr=0x0; return 1.0;}
2156 Int_t etaBin = effH->GetXaxis()->FindBin(eta);
2157 if(eta<effH->GetXaxis()->GetXmin()) etaBin=1;
2158 if(eta>effH->GetXaxis()->GetXmax()) etaBin=effH->GetXaxis()->GetNbins();
2159 Int_t phiBin = effH->GetYaxis()->FindBin(trdPhi);
2160 if(trdPhi<effH->GetYaxis()->GetXmin()) phiBin=1;
2161 if(trdPhi>effH->GetYaxis()->GetXmax()) phiBin=effH->GetYaxis()->GetNbins();
2162 Int_t poutBin = effH->GetZaxis()->FindBin(pout);
2163 if(pout<effH->GetZaxis()->GetXmin()) poutBin=1;
2164 if(pout>effH->GetZaxis()->GetXmax()) poutBin=effH->GetZaxis()->GetNbins();
2165 Double_t eff = effH->GetBinContent(etaBin, phiBin, poutBin);
2166 effErr = effH->GetBinError(etaBin, phiBin, poutBin);
2167 if(eff<-0.0001) {
2168 effErr = 0.0;
2169 eff = 1.0;
2170 }
2171 return eff;
2172}
2173
2174
572b0139 2175inline void AliDielectronVarManager::SetEvent(AliVEvent * const ev)
2176{
2177
2178 fgEvent = ev;
2179 if (fgKFVertex) delete fgKFVertex;
2180 fgKFVertex=0x0;
2e02dba4 2181 if (!ev) return;
2182 if (ev->GetPrimaryVertex()) fgKFVertex=new AliKFVertex(*ev->GetPrimaryVertex());
ba15fdfb 2183
2184 for (Int_t i=0; i<AliDielectronVarManager::kNMaxValues;++i) fgData[i]=0.;
2185 AliDielectronVarManager::Fill(fgEvent, fgData);
572b0139 2186}
2a14a7b1 2187
5720c765 2188inline void AliDielectronVarManager::SetEventData(const Double_t data[AliDielectronVarManager::kNMaxValues])
2189{
2190 for (Int_t i=0; i<kNMaxValues;++i) fgData[i]=0.;
2191 for (Int_t i=kPairMax; i<kNMaxValues;++i) fgData[i]=data[i];
2192}
2193
2a14a7b1 2194
2195inline Bool_t AliDielectronVarManager::GetDCA(const AliAODTrack *track, Double_t d0z0[2])
2196{
2197 if(track->TestBit(AliAODTrack::kIsDCA)){
2198 d0z0[0]=track->DCA();
2199 d0z0[1]=track->ZAtDCA();
2200 return kTRUE;
2201 }
2202
67fd1119 2203 Bool_t ok=kFALSE;
2204 if(fgEvent) {
2205 Double_t covd0z0[3];
31b0e8c5 2206 //AliAODTrack copy(*track);
2207 AliExternalTrackParam etp; etp.CopyFromVTrack(track);
2208
2209 Float_t xstart = etp.GetX();
2210 if(xstart>3.) {
2211 d0z0[0]=-999.;
2212 d0z0[1]=-999.;
2213 //printf("This method can be used only for propagation inside the beam pipe \n");
2214 return kFALSE;
2215 }
2216
2217
67fd1119 2218 AliAODVertex *vtx =(AliAODVertex*)(fgEvent->GetPrimaryVertex());
2219 Double_t fBzkG = fgEvent->GetMagneticField(); // z componenent of field in kG
31b0e8c5 2220 ok = etp.PropagateToDCA(vtx,fBzkG,kVeryBig,d0z0,covd0z0);
2221 //ok = copy.PropagateToDCA(vtx,fBzkG,kVeryBig,d0z0,covd0z0);
67fd1119 2222 }
2a14a7b1 2223 if(!ok){
2224 d0z0[0]=-999.;
2225 d0z0[1]=-999.;
2226 }
2227 return ok;
2228}
2229
5720c765 2230inline void AliDielectronVarManager::SetTPCEventPlane(AliEventplane *const evplane)
2231{
2232
2233 fgTPCEventPlane = evplane;
c315310a 2234 FillVarTPCEventPlane(evplane,fgData);
5720c765 2235 // for (Int_t i=0; i<AliDielectronVarManager::kNMaxValues;++i) fgData[i]=0.;
2236 // AliDielectronVarManager::Fill(fgEvent, fgData);
2237}
2238
2239
2240//_________________________________________________________________
2241inline void AliDielectronVarManager::GetVzeroRP(const AliVEvent* event, Double_t* qvec, Int_t sideOption) {
2242 //
2243 // Get the reaction plane from the VZERO detector for a given harmonic
2244 //
2245 // sideOption = 0- V0A, 1- V0C, 2-both
2246 // Q{x,y} = SUM_i mult(i) * {cos(n*phi_i), sin(n*phi_i)}
2247 // phi_i - phi angle of the VZERO sector i
2248 // Each sector covers 45 degrees(8 sectors per ring). Middle of sector 0 is at 45/2
2249 // channel 0: 22.5
2250 // 1: 22.5+45
2251 // 2: 22.5+45*2
2252 // ...
2253 // at the next ring continues the same
2254 // channel 8: 22.5
2255 // channel 9: 22.5 + 45
2256 // ...
2257 const Double_t kX[8] = {0.92388, 0.38268, -0.38268, -0.92388, -0.92388, -0.38268, 0.38268, 0.92388}; // cosines of the angles of the VZERO sectors (8 per ring)
2258 const Double_t kY[8] = {0.38268, 0.92388, 0.92388, 0.38268, -0.38268, -0.92388, -0.92388, -0.38268}; // sines -- " --
2259 Int_t phi;
2260 Float_t mult;
5720c765 2261
a823f01b 2262 // get centrality and vertex for this event
2263 Double_t centralitySPD = -1; Double_t vtxZ = -999.;
2264 if(event->IsA() == AliESDEvent::Class()) {
2265 const AliESDEvent* esdEv = static_cast<const AliESDEvent*>(event);
2266 AliCentrality *esdCentrality = const_cast<AliESDEvent*>(esdEv)->GetCentrality();
2267 if(esdCentrality) centralitySPD = esdCentrality->GetCentralityPercentile("CL1");
2268 }
2269 if(event->IsA() == AliAODEvent::Class()) {
2270 const AliAODEvent* aodEv = static_cast<const AliAODEvent*>(event);
2271 AliAODHeader *header = aodEv->GetHeader();
2272 AliCentrality *aodCentrality = header->GetCentralityP();
2273 if(aodCentrality) centralitySPD = aodCentrality->GetCentralityPercentile("CL1");
2274 }
2275 const AliVVertex *primVtx = event->GetPrimaryVertex();
2276 if(!primVtx) return;
2277 vtxZ = primVtx->GetZ();
2278 if(TMath::Abs(vtxZ)>10.) return;
2279 if(centralitySPD<0. || centralitySPD>80.) return;
2280
2281 Int_t binCent = -1; Int_t binVtx = -1;
2282 if(fgVZEROCalib[0]) {
2283 binVtx = fgVZEROCalib[0]->GetXaxis()->FindBin(vtxZ);
2284 binCent = fgVZEROCalib[0]->GetYaxis()->FindBin(centralitySPD);
2285 }
2286
2287 AliVVZERO* vzero = event->GetVZEROData();
2288 Double_t average = 0.0;
5720c765 2289 for(Int_t iChannel=0; iChannel<64; ++iChannel) {
2290 if(iChannel<32 && sideOption==0) continue;
2291 if(iChannel>=32 && sideOption==1) continue;
2292 phi=iChannel%8;
a823f01b 2293 mult = vzero->GetMultiplicity(iChannel);
2294 if(fgVZEROCalib[iChannel])
2295 average = fgVZEROCalib[iChannel]->GetBinContent(binVtx, binCent);
2296 if(average>1.0e-10 && mult>0.5)
2297 mult /= average;
2298 else
2299 mult = 0.0;
5720c765 2300 // 2nd harmonic
2301 qvec[0] += mult*(2.0*TMath::Power(kX[phi],2.0)-1);
2302 qvec[1] += mult*(2.0*kX[phi]*kY[phi]);
2303 } // end loop over channels
a823f01b 2304
2305 // do recentering
2306 if(fgVZERORecentering[0][0]) {
2307// printf("vzero: %p\n",fgVZERORecentering[0][0]);
2308 Int_t binCentRecenter = -1; Int_t binVtxRecenter = -1;
2309 binCentRecenter = fgVZERORecentering[0][0]->GetXaxis()->FindBin(centralitySPD);
2310 binVtxRecenter = fgVZERORecentering[0][0]->GetYaxis()->FindBin(vtxZ);
2311 if(sideOption==0) { // side A
2312 qvec[0] -= fgVZERORecentering[0][0]->GetBinContent(binCentRecenter, binVtxRecenter);
2313 qvec[1] -= fgVZERORecentering[0][1]->GetBinContent(binCentRecenter, binVtxRecenter);
2314 }
2315 if(sideOption==1) { // side C
2316 qvec[0] -= fgVZERORecentering[1][0]->GetBinContent(binCentRecenter, binVtxRecenter);
2317 qvec[1] -= fgVZERORecentering[1][1]->GetBinContent(binCentRecenter, binVtxRecenter);
2318 }
2319 if(sideOption==2) { // side A and C together
2320 qvec[0] -= fgVZERORecentering[0][0]->GetBinContent(binCentRecenter, binVtxRecenter);
2321 qvec[0] -= fgVZERORecentering[1][0]->GetBinContent(binCentRecenter, binVtxRecenter);
2322 qvec[1] -= fgVZERORecentering[0][1]->GetBinContent(binCentRecenter, binVtxRecenter);
2323 qvec[1] -= fgVZERORecentering[1][1]->GetBinContent(binCentRecenter, binVtxRecenter);
2324 }
2325 }
2326
2327 // calculate the reaction plane
5720c765 2328 if(TMath::Abs(qvec[0])>1.0e-10)
2329 qvec[2] = TMath::ATan2(qvec[1],qvec[0])/2.0;
2330}
2331
99345a64 2332//______________________________________________________________________________
2333inline AliAODVertex* AliDielectronVarManager::GetVertex(const AliAODEvent* event, AliAODVertex::AODVtx_t vtype) {
2334 // Get vertex
2335 Int_t nVertices=event->GetNumberOfVertices();
2336 for(Int_t iVert=0; iVert<nVertices; iVert++){
2337 AliAODVertex *v=event->GetVertex(iVert);
2338 // printf(" vtx %d contrib %d daughters %d \n ",v->GetType(),v->GetNContributors(), v->GetNDaughters());
2339 if(v->GetType()==vtype) return v;
2340 }
2341 return 0;
2342}
5720c765 2343
b2a297fa 2344/*
2345inline void AliDielectronVarManager::FillValues(const TParticle *particle, Double_t *values)
2346{
2347 //
2348 // Fill track information available for histogramming into an array
2349 //
2350
2351 // Fill TParticle interface information
2352 values[AliDielectronVarManager::kPx] = particle->Px();
2353 values[AliDielectronVarManager::kPy] = particle->Py();
2354 values[AliDielectronVarManager::kPz] = particle->Pz();
2355 values[AliDielectronVarManager::kPt] = particle->Pt();
2356 values[AliDielectronVarManager::kP] = particle->P();
2357
2358 values[AliDielectronVarManager::kXv] = particle->Vx();
2359 values[AliDielectronVarManager::kYv] = particle->Vy();
2360 values[AliDielectronVarManager::kZv] = particle->Vz();
2361
2362 values[AliDielectronVarManager::kOneOverPt] = 1./particle->Pt();
2363 values[AliDielectronVarManager::kPhi] = particle->Phi();
2364 values[AliDielectronVarManager::kTheta] =
2365 values[AliDielectronVarManager::kEta] = particle->Eta();
2366 values[AliDielectronVarManager::kY] =
2367
2368 values[AliDielectronVarManager::kE] = particle->Energy();
2369 values[AliDielectronVarManager::kM] = particle->GetMass();
2370
2371 values[AliDielectronVarManager::kCharge] = particle->GetPDG()->Charge()/3; // uggly
2372
2373}*/
2374
2375#endif
fb7d2d99 2376