Skip to content

Commit e2ca4e6

Browse files
dsekihatdsekihat
authored
[PWGEM/Dilepton] add impact parameters in taggingHFE (#16197)
1 parent 14faea5 commit e2ca4e6

3 files changed

Lines changed: 112 additions & 88 deletions

File tree

‎PWGEM/Dilepton/DataModel/lmeeMLTables.h‎

Lines changed: 25 additions & 22 deletions
Original file line numberDiff line numberDiff line change
@@ -219,11 +219,11 @@ DECLARE_SOA_INDEX_COLUMN(EMMLLepton, emmllepton); //! most propable emeventId
219219
DECLARE_SOA_COLUMN(RapidityV0, rapidityV0, float); //! rapidity of associated V0
220220
DECLARE_SOA_COLUMN(RapidityC, rapidityC, float); //! rapidity of associated Cascade
221221

222-
DECLARE_SOA_COLUMN(ImpParXY, impParXY, float); //! impact parameter for V0/Cascade in XY plane
223-
DECLARE_SOA_COLUMN(ImpParZ, impParZ, float); //! impact parameter for V0/Cascade in XY plane
224-
DECLARE_SOA_COLUMN(ImpParCYY, impParCYY, float); //! sigma of impact parameter for V0/Cascade in XY plane
225-
DECLARE_SOA_COLUMN(ImpParCZY, impParCZY, float); //! sigma of impact parameter for V0/Cascade in XY plane
226-
DECLARE_SOA_COLUMN(ImpParCZZ, impParCZZ, float); //! sigma of impact parameter for V0/Cascade in XY plane
222+
DECLARE_SOA_COLUMN(ImpParHXY, impParHXY, float); //! impact parameter for V0/Cascade in XY plane
223+
DECLARE_SOA_COLUMN(ImpParHZ, impParHZ, float); //! impact parameter for V0/Cascade in Z plane
224+
DECLARE_SOA_COLUMN(ImpParHCYY, impParHCYY, float); //! sigma of impact parameter for V0/Cascade in XY
225+
DECLARE_SOA_COLUMN(ImpParHCZY, impParHCZY, float); //! sigma of impact parameter for V0/Cascade, correlaion term
226+
DECLARE_SOA_COLUMN(ImpParHCZZ, impParHCZZ, float); //! sigma of impact parameter for V0/Cascade in Z
227227

228228
DECLARE_SOA_COLUMN(V0CosPA, v0copa, float); //! cosPA of V0
229229
DECLARE_SOA_COLUMN(V0CosPAXY, v0copaXY, float); //! cosPA of V0 in XY plane
@@ -236,12 +236,18 @@ DECLARE_SOA_COLUMN(DcaLH, dcalh, float); //! DCA between lepton and hadron
236236
DECLARE_SOA_COLUMN(CosPA, cospa, float); //! flight direction of LH pair
237237
DECLARE_SOA_COLUMN(CosPAXY, cospaXY, float); //! flight direction of LH pair in XY
238238

239-
DECLARE_SOA_COLUMN(Lxyz, lxyz, float); //! decay length of LH pair
240-
DECLARE_SOA_COLUMN(LxyzSigma, lxyzSigma, float); //! decay length resolution of LH pair
241-
DECLARE_SOA_COLUMN(Lxy, lxy, float); //! decay length of LH pair
242-
DECLARE_SOA_COLUMN(LxySigma, lxySigma, float); //! decay length resolution of LH pair
243-
DECLARE_SOA_COLUMN(Lz, lz, float); //! decay length of LH pair
244-
DECLARE_SOA_COLUMN(LzSigma, lzSigma, float); //! decay length resolution of LH pair
239+
DECLARE_SOA_COLUMN(Lxy, lxy, float); //! decay length of LH pair
240+
DECLARE_SOA_COLUMN(Lz, lz, float); //! decay length of LH pair
241+
DECLARE_SOA_COLUMN(Lxyz, lxyz, float); //! decay length of LH pair
242+
DECLARE_SOA_COLUMN(LxyErr, lxyErr, float); //! decay length resolution of LH pair
243+
DECLARE_SOA_COLUMN(LzErr, lzErr, float); //! decay length resolution of LH pair
244+
DECLARE_SOA_COLUMN(LxyzErr, lxyzErr, float); //! decay length resolution of LH pair
245+
246+
DECLARE_SOA_COLUMN(ImpParXY, impParXY, float); //! impact parameter for LH in XY plane
247+
DECLARE_SOA_COLUMN(ImpParZ, impParZ, float); //! impact parameter for LH in Z plane
248+
DECLARE_SOA_COLUMN(ImpParCYY, impParCYY, float); //! sigma of impact parameter for LH in XY
249+
DECLARE_SOA_COLUMN(ImpParCZY, impParCZY, float); //! sigma of impact parameter for LH, correlation term
250+
DECLARE_SOA_COLUMN(ImpParCZZ, impParCZZ, float); //! sigma of impact parameter for LH in Z
245251

246252
DECLARE_SOA_COLUMN(PdgCodeH, pdgCodeH, int); //! pdg code of associated hadron
247253
DECLARE_SOA_COLUMN(PdgCodeIM, pdgCodeIM, int); //! pdg code of intermediate hadron from HF hadrons. e.g K*, D*
@@ -256,9 +262,8 @@ DECLARE_SOA_TABLE(EMMLLHPairs, "AOD", "EMMLLHPAIR", //!
256262
pidtpc::TPCNSigmaKa, pidtof::TOFNSigmaKa,
257263
pidtpc::TPCNSigmaPr, pidtof::TOFNSigmaPr,
258264
emmllhpair::Mass, emmllhpair::DcaLH, emmllhpair::CosPA, emmllhpair::CosPAXY,
259-
emmllhpair::Lxyz, emmllhpair::LxyzSigma,
260-
emmllhpair::Lxy, emmllhpair::LxySigma,
261-
emmllhpair::Lz, emmllhpair::LzSigma,
265+
emmllhpair::Lxy, emmllhpair::Lz, emmllhpair::Lxyz, emmllhpair::LxyErr, emmllhpair::LzErr, emmllhpair::LxyzErr,
266+
emmllhpair::ImpParXY, emmllhpair::ImpParZ, emmllhpair::ImpParCYY, emmllhpair::ImpParCZY, emmllhpair::ImpParCZZ,
262267
emmllhpair::PdgCodeH, emmllhpair::PdgCodeIM, emmllhpair::FoundCommonMother);
263268
// iterators
264269
using EMMLLHPair = EMMLLHPairs::iterator;
@@ -267,11 +272,10 @@ DECLARE_SOA_TABLE(EMMLLV0Pairs, "AOD", "EMMLLV0PAIR", //!
267272
emmllhpair::EMMLLeptonId,
268273
track::Pt, emmllhpair::RapidityV0,
269274
emmllhpair::V0CosPA, emmllhpair::V0CosPAXY,
270-
emmllhpair::ImpParXY, emmllhpair::ImpParZ, emmllhpair::ImpParCYY, emmllhpair::ImpParCZY, emmllhpair::ImpParCZZ,
275+
emmllhpair::ImpParHXY, emmllhpair::ImpParHZ, emmllhpair::ImpParHCYY, emmllhpair::ImpParHCZY, emmllhpair::ImpParHCZZ,
271276
emmllhpair::Mass, emmllhpair::DcaLH, emmllhpair::CosPA, emmllhpair::CosPAXY,
272-
emmllhpair::Lxyz, emmllhpair::LxyzSigma,
273-
emmllhpair::Lxy, emmllhpair::LxySigma,
274-
emmllhpair::Lz, emmllhpair::LzSigma,
277+
emmllhpair::Lxy, emmllhpair::Lz, emmllhpair::Lxyz, emmllhpair::LxyErr, emmllhpair::LzErr, emmllhpair::LxyzErr,
278+
emmllhpair::ImpParXY, emmllhpair::ImpParZ, emmllhpair::ImpParCYY, emmllhpair::ImpParCZY, emmllhpair::ImpParCZZ,
275279
emmllhpair::PdgCodeH, emmllhpair::PdgCodeIM, emmllhpair::FoundCommonMother);
276280
// iterators
277281
using EMMLLV0Pair = EMMLLV0Pairs::iterator;
@@ -280,11 +284,10 @@ DECLARE_SOA_TABLE(EMMLLCascPairs, "AOD", "EMMLLCAPAIR", //!
280284
emmllhpair::EMMLLeptonId,
281285
track::Pt, emmllhpair::RapidityC,
282286
emmllhpair::CascCosPA, emmllhpair::CascCosPAXY,
283-
emmllhpair::ImpParXY, emmllhpair::ImpParZ, emmllhpair::ImpParCYY, emmllhpair::ImpParCZY, emmllhpair::ImpParCZZ,
287+
emmllhpair::ImpParHXY, emmllhpair::ImpParHZ, emmllhpair::ImpParHCYY, emmllhpair::ImpParHCZY, emmllhpair::ImpParHCZZ,
284288
emmllhpair::Mass, emmllhpair::DcaLH, emmllhpair::CosPA, emmllhpair::CosPAXY,
285-
emmllhpair::Lxyz, emmllhpair::LxyzSigma,
286-
emmllhpair::Lxy, emmllhpair::LxySigma,
287-
emmllhpair::Lz, emmllhpair::LzSigma,
289+
emmllhpair::Lxy, emmllhpair::Lz, emmllhpair::Lxyz, emmllhpair::LxyErr, emmllhpair::LzErr, emmllhpair::LxyzErr,
290+
emmllhpair::ImpParXY, emmllhpair::ImpParZ, emmllhpair::ImpParCYY, emmllhpair::ImpParCZY, emmllhpair::ImpParCZZ,
288291
emmllhpair::PdgCodeH, emmllhpair::PdgCodeIM, emmllhpair::FoundCommonMother);
289292
// iterators
290293
using EMMLLCascPair = EMMLLCascPairs::iterator;

‎PWGEM/Dilepton/Tasks/taggingHFE.cxx‎

Lines changed: 24 additions & 36 deletions
Original file line numberDiff line numberDiff line change
@@ -1710,9 +1710,8 @@ struct taggingHFE {
17101710
kaon.tpcNSigmaKa(), tofNSigmaKa,
17111711
kaon.tpcNSigmaPr(), tofNSigmaPr,
17121712
eKpair.mass, eKpair.dca2legs, eKpair.cospa, eKpair.cospaXY,
1713-
eKpair.lxyz, eKpair.lxyzErr,
1714-
eKpair.lxy, eKpair.lxyErr,
1715-
eKpair.lz, eKpair.lzErr,
1713+
eKpair.lxy, eKpair.lz, eKpair.lxyz, eKpair.lxyErr, eKpair.lzErr, eKpair.lxyzErr,
1714+
eKpair.impParXY, eKpair.impParZ, eKpair.impParCYY, eKpair.impParCZY, eKpair.impParCZZ,
17161715
mckaon.pdgCode(), pdgCodeIM, foundCommonMother);
17171716

17181717
} // end of kaon loop
@@ -1774,9 +1773,8 @@ struct taggingHFE {
17741773
kaon.tpcNSigmaKa(), tofNSigmaKa,
17751774
kaon.tpcNSigmaPr(), tofNSigmaPr,
17761775
eKpair.mass, eKpair.dca2legs, eKpair.cospa, eKpair.cospaXY,
1777-
eKpair.lxyz, eKpair.lxyzErr,
1778-
eKpair.lxy, eKpair.lxyErr,
1779-
eKpair.lz, eKpair.lzErr,
1776+
eKpair.lxy, eKpair.lz, eKpair.lxyz, eKpair.lxyErr, eKpair.lzErr, eKpair.lxyzErr,
1777+
eKpair.impParXY, eKpair.impParZ, eKpair.impParCYY, eKpair.impParCZY, eKpair.impParCZZ,
17801778
mckaon.pdgCode(), pdgCodeIM, foundCommonMother);
17811779

17821780
} // end of kaon loop
@@ -1851,9 +1849,8 @@ struct taggingHFE {
18511849
RecoDecay::cpaXY(std::array<float, 3>{collision.posX(), collision.posY(), collision.posZ()}, std::array<float, 3>{v0.x(), v0.y(), v0.z()}, std::array<float, 3>{v0.px(), v0.py(), v0.pz()}),
18521850
impactParameterV0.getY(), impactParameterV0.getZ(), impactParameterV0.getSigmaY2(), impactParameterV0.getSigmaYZ(), impactParameterV0.getSigmaZ2(),
18531851
eV0pair.mass, eV0pair.dca2legs, eV0pair.cospa, eV0pair.cospaXY,
1854-
eV0pair.lxyz, eV0pair.lxyzErr,
1855-
eV0pair.lxy, eV0pair.lxyErr,
1856-
eV0pair.lz, eV0pair.lzErr,
1852+
eV0pair.lxy, eV0pair.lz, eV0pair.lxyz, eV0pair.lxyErr, eV0pair.lzErr, eV0pair.lxyzErr,
1853+
eV0pair.impParXY, eV0pair.impParZ, eV0pair.impParCYY, eV0pair.impParCZY, eV0pair.impParCZZ,
18571854
pdgCodeV0, pdgCodeIM, foundCommonMother);
18581855

18591856
} // end of K0S loop
@@ -1912,9 +1909,8 @@ struct taggingHFE {
19121909
RecoDecay::cpaXY(std::array<float, 3>{collision.posX(), collision.posY(), collision.posZ()}, std::array<float, 3>{v0.x(), v0.y(), v0.z()}, std::array<float, 3>{v0.px(), v0.py(), v0.pz()}),
19131910
impactParameterV0.getY(), impactParameterV0.getZ(), impactParameterV0.getSigmaY2(), impactParameterV0.getSigmaYZ(), impactParameterV0.getSigmaZ2(),
19141911
eV0pair.mass, eV0pair.dca2legs, eV0pair.cospa, eV0pair.cospaXY,
1915-
eV0pair.lxyz, eV0pair.lxyzErr,
1916-
eV0pair.lxy, eV0pair.lxyErr,
1917-
eV0pair.lz, eV0pair.lzErr,
1912+
eV0pair.lxy, eV0pair.lz, eV0pair.lxyz, eV0pair.lxyErr, eV0pair.lzErr, eV0pair.lxyzErr,
1913+
eV0pair.impParXY, eV0pair.impParZ, eV0pair.impParCYY, eV0pair.impParCZY, eV0pair.impParCZZ,
19181914
pdgCodeV0, pdgCodeIM, foundCommonMother);
19191915

19201916
} // end of Lambda loop
@@ -1978,9 +1974,8 @@ struct taggingHFE {
19781974
RecoDecay::cpaXY(std::array<float, 3>{collision.posX(), collision.posY(), collision.posZ()}, std::array<float, 3>{cascade.x(), cascade.y(), cascade.z()}, std::array<float, 3>{cascade.px(), cascade.py(), cascade.pz()}),
19791975
impactParameterCasc.getY(), impactParameterCasc.getZ(), impactParameterCasc.getSigmaY2(), impactParameterCasc.getSigmaYZ(), impactParameterCasc.getSigmaZ2(),
19801976
eCpair.mass, eCpair.dca2legs, eCpair.cospa, eCpair.cospaXY,
1981-
eCpair.lxyz, eCpair.lxyzErr,
1982-
eCpair.lxy, eCpair.lxyErr,
1983-
eCpair.lz, eCpair.lzErr,
1977+
eCpair.lxy, eCpair.lz, eCpair.lxyz, eCpair.lxyErr, eCpair.lzErr, eCpair.lxyzErr,
1978+
eCpair.impParXY, eCpair.impParZ, eCpair.impParCYY, eCpair.impParCZY, eCpair.impParCZZ,
19841979
pdgCodeCascade, pdgCodeIM, foundCommonMother);
19851980

19861981
} // end of Xi- loop
@@ -2044,9 +2039,8 @@ struct taggingHFE {
20442039
RecoDecay::cpaXY(std::array<float, 3>{collision.posX(), collision.posY(), collision.posZ()}, std::array<float, 3>{cascade.x(), cascade.y(), cascade.z()}, std::array<float, 3>{cascade.px(), cascade.py(), cascade.pz()}),
20452040
impactParameterCasc.getY(), impactParameterCasc.getZ(), impactParameterCasc.getSigmaY2(), impactParameterCasc.getSigmaYZ(), impactParameterCasc.getSigmaZ2(),
20462041
eCpair.mass, eCpair.dca2legs, eCpair.cospa, eCpair.cospaXY,
2047-
eCpair.lxyz, eCpair.lxyzErr,
2048-
eCpair.lxy, eCpair.lxyErr,
2049-
eCpair.lz, eCpair.lzErr,
2042+
eCpair.lxy, eCpair.lz, eCpair.lxyz, eCpair.lxyErr, eCpair.lzErr, eCpair.lxyzErr,
2043+
eCpair.impParXY, eCpair.impParZ, eCpair.impParCYY, eCpair.impParCZY, eCpair.impParCZZ,
20502044
pdgCodeCascade, pdgCodeIM, foundCommonMother);
20512045

20522046
} // end of Omega- loop
@@ -2137,9 +2131,8 @@ struct taggingHFE {
21372131
kaon.tpcNSigmaKa(), tofNSigmaKa,
21382132
kaon.tpcNSigmaPr(), tofNSigmaPr,
21392133
eKpair.mass, eKpair.dca2legs, eKpair.cospa, eKpair.cospaXY,
2140-
eKpair.lxyz, eKpair.lxyzErr,
2141-
eKpair.lxy, eKpair.lxyErr,
2142-
eKpair.lz, eKpair.lzErr,
2134+
eKpair.lxy, eKpair.lz, eKpair.lxyz, eKpair.lxyErr, eKpair.lzErr, eKpair.lxyzErr,
2135+
eKpair.impParXY, eKpair.impParZ, eKpair.impParCYY, eKpair.impParCZY, eKpair.impParCZZ,
21432136
mckaon.pdgCode(), pdgCodeIM, foundCommonMother);
21442137

21452138
} // end of kaon loop
@@ -2202,9 +2195,8 @@ struct taggingHFE {
22022195
kaon.tpcNSigmaKa(), tofNSigmaKa,
22032196
kaon.tpcNSigmaPr(), tofNSigmaPr,
22042197
eKpair.mass, eKpair.dca2legs, eKpair.cospa, eKpair.cospaXY,
2205-
eKpair.lxyz, eKpair.lxyzErr,
2206-
eKpair.lxy, eKpair.lxyErr,
2207-
eKpair.lz, eKpair.lzErr,
2198+
eKpair.lxy, eKpair.lz, eKpair.lxyz, eKpair.lxyErr, eKpair.lzErr, eKpair.lxyzErr,
2199+
eKpair.impParXY, eKpair.impParZ, eKpair.impParCYY, eKpair.impParCZY, eKpair.impParCZZ,
22082200
mckaon.pdgCode(), pdgCodeIM, foundCommonMother);
22092201

22102202
} // end of kaon loop
@@ -2278,9 +2270,8 @@ struct taggingHFE {
22782270
RecoDecay::cpaXY(std::array<float, 3>{collision.posX(), collision.posY(), collision.posZ()}, std::array<float, 3>{v0.x(), v0.y(), v0.z()}, std::array<float, 3>{v0.px(), v0.py(), v0.pz()}),
22792271
impactParameterV0.getY(), impactParameterV0.getZ(), impactParameterV0.getSigmaY2(), impactParameterV0.getSigmaYZ(), impactParameterV0.getSigmaZ2(),
22802272
eV0pair.mass, eV0pair.dca2legs, eV0pair.cospa, eV0pair.cospaXY,
2281-
eV0pair.lxyz, eV0pair.lxyzErr,
2282-
eV0pair.lxy, eV0pair.lxyErr,
2283-
eV0pair.lz, eV0pair.lzErr,
2273+
eV0pair.lxy, eV0pair.lz, eV0pair.lxyz, eV0pair.lxyErr, eV0pair.lzErr, eV0pair.lxyzErr,
2274+
eV0pair.impParXY, eV0pair.impParZ, eV0pair.impParCYY, eV0pair.impParCZY, eV0pair.impParCZZ,
22842275
pdgCodeV0, pdgCodeIM, foundCommonMother);
22852276

22862277
} // end of K0S loop
@@ -2339,9 +2330,8 @@ struct taggingHFE {
23392330
RecoDecay::cpaXY(std::array<float, 3>{collision.posX(), collision.posY(), collision.posZ()}, std::array<float, 3>{v0.x(), v0.y(), v0.z()}, std::array<float, 3>{v0.px(), v0.py(), v0.pz()}),
23402331
impactParameterV0.getY(), impactParameterV0.getZ(), impactParameterV0.getSigmaY2(), impactParameterV0.getSigmaYZ(), impactParameterV0.getSigmaZ2(),
23412332
eV0pair.mass, eV0pair.dca2legs, eV0pair.cospa, eV0pair.cospaXY,
2342-
eV0pair.lxyz, eV0pair.lxyzErr,
2343-
eV0pair.lxy, eV0pair.lxyErr,
2344-
eV0pair.lz, eV0pair.lzErr,
2333+
eV0pair.lxy, eV0pair.lz, eV0pair.lxyz, eV0pair.lxyErr, eV0pair.lzErr, eV0pair.lxyzErr,
2334+
eV0pair.impParXY, eV0pair.impParZ, eV0pair.impParCYY, eV0pair.impParCZY, eV0pair.impParCZZ,
23452335
pdgCodeV0, pdgCodeIM, foundCommonMother);
23462336

23472337
} // end of Anti-Lambda loop
@@ -2405,9 +2395,8 @@ struct taggingHFE {
24052395
RecoDecay::cpaXY(std::array<float, 3>{collision.posX(), collision.posY(), collision.posZ()}, std::array<float, 3>{cascade.x(), cascade.y(), cascade.z()}, std::array<float, 3>{cascade.px(), cascade.py(), cascade.pz()}),
24062396
impactParameterCasc.getY(), impactParameterCasc.getZ(), impactParameterCasc.getSigmaY2(), impactParameterCasc.getSigmaYZ(), impactParameterCasc.getSigmaZ2(),
24072397
eCpair.mass, eCpair.dca2legs, eCpair.cospa, eCpair.cospaXY,
2408-
eCpair.lxyz, eCpair.lxyzErr,
2409-
eCpair.lxy, eCpair.lxyErr,
2410-
eCpair.lz, eCpair.lzErr,
2398+
eCpair.lxy, eCpair.lz, eCpair.lxyz, eCpair.lxyErr, eCpair.lzErr, eCpair.lxyzErr,
2399+
eCpair.impParXY, eCpair.impParZ, eCpair.impParCYY, eCpair.impParCZY, eCpair.impParCZZ,
24112400
pdgCodeCascade, pdgCodeIM, foundCommonMother);
24122401

24132402
} // end of Xi+ loop
@@ -2471,9 +2460,8 @@ struct taggingHFE {
24712460
RecoDecay::cpaXY(std::array<float, 3>{collision.posX(), collision.posY(), collision.posZ()}, std::array<float, 3>{cascade.x(), cascade.y(), cascade.z()}, std::array<float, 3>{cascade.px(), cascade.py(), cascade.pz()}),
24722461
impactParameterCasc.getY(), impactParameterCasc.getZ(), impactParameterCasc.getSigmaY2(), impactParameterCasc.getSigmaYZ(), impactParameterCasc.getSigmaZ2(),
24732462
eCpair.mass, eCpair.dca2legs, eCpair.cospa, eCpair.cospaXY,
2474-
eCpair.lxyz, eCpair.lxyzErr,
2475-
eCpair.lxy, eCpair.lxyErr,
2476-
eCpair.lz, eCpair.lzErr,
2463+
eCpair.lxy, eCpair.lz, eCpair.lxyz, eCpair.lxyErr, eCpair.lzErr, eCpair.lxyzErr,
2464+
eCpair.impParXY, eCpair.impParZ, eCpair.impParCYY, eCpair.impParCZY, eCpair.impParCZZ,
24772465
pdgCodeCascade, pdgCodeIM, foundCommonMother);
24782466

24792467
} // end of Omega+ loop

0 commit comments

Comments
 (0)