Skip to content

Commit 1bba6d8

Browse files
dsekihatdsekihat
authored
[PWGEM/PhotonMeson] update pcm data table and tasks (#10425)
1 parent 229e8d1 commit 1bba6d8

File tree

8 files changed

+20
-7
lines changed

8 files changed

+20
-7
lines changed

‎PWGEM/PhotonMeson/DataModel/gammaTables.h‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -175,7 +175,7 @@ DECLARE_SOA_TABLE(V0Legs, "AOD", "V0LEG", //!
175175
track::TPCChi2NCl, track::TPCInnerParam,
176176
track::TPCSignal, pidtpc::TPCNSigmaEl, pidtpc::TPCNSigmaPi,
177177
track::ITSClusterSizes, track::ITSChi2NCl, track::DetectorMap,
178-
track::X, track::Y, track::Z, track::Tgl,
178+
track::X, track::Y, track::Z, track::Tgl, track::C1Pt21Pt2,
179179

180180
// dynamic column
181181
v0leg::P<v0leg::Px, v0leg::Py, v0leg::Pz>,
@@ -284,7 +284,7 @@ using V0PhotonKFPrefilterBitDerived = V0PhotonsKFPrefilterBitDerived::iterator;
284284
DECLARE_SOA_TABLE(EMPrimaryElectronsFromDalitz, "AOD", "EMPRIMARYELDA", //!
285285
o2::soa::Index<>, emprimaryelectron::CollisionId,
286286
emprimaryelectron::TrackId, emprimaryelectron::Sign,
287-
track::Pt, track::Eta, track::Phi, track::DcaXY, track::DcaZ, track::CYY, track::CZY, track::CZZ,
287+
track::Pt, track::Eta, track::Phi, track::DcaXY, track::DcaZ, track::CYY, track::CZY, track::CZZ, track::C1Pt21Pt2,
288288
track::TPCNClsFindable, track::TPCNClsFindableMinusFound, track::TPCNClsFindableMinusCrossedRows, track::TPCNClsShared,
289289
track::TPCChi2NCl, track::TPCInnerParam,
290290
track::TPCSignal, pidtpc::TPCNSigmaEl, pidtpc::TPCNSigmaPi,

‎PWGEM/PhotonMeson/TableProducer/photonconversionbuilder.cxx‎

Lines changed: 6 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -166,6 +166,7 @@ struct PhotonConversionBuilder {
166166
{"V0/hPCA_diffX", "PCA vs. trackiu X - R_{xy};distance btween 2 legs (cm);min trackiu X - R_{xy} (cm)", {HistType::kTH2F, {{500, 0.0f, 5.f}, {100, -50.0, 50.0f}}}},
167167
{"V0Leg/hPt", "pT of leg at SV;p_{T,e} (GeV/c)", {HistType::kTH1F, {{1000, 0.0f, 10.0f}}}},
168168
{"V0Leg/hEtaPhi", "#eta vs. #varphi of leg at SV;#varphi (rad.);#eta", {HistType::kTH2F, {{72, 0.0f, 2 * M_PI}, {200, -1, +1}}}},
169+
{"V0Leg/hRelDeltaPt", "pT resolution;p_{T} (GeV/c);#Deltap_{T}/p_{T}", {HistType::kTH2F, {{1000, 0.f, 10.f}, {100, 0, 1}}}},
169170
{"V0Leg/hDCAxyz", "DCA xy vs. z to PV;DCA_{xy} (cm);DCA_{z} (cm)", {HistType::kTH2F, {{200, -50.f, 50.f}, {200, -50.f, +50.f}}}},
170171
{"V0Leg/hdEdx_Pin", "TPC dE/dx vs. p_{in};p_{in} (GeV/c);TPC dE/dx", {HistType::kTH2F, {{1000, 0.f, 10.f}, {200, 0.f, 200.f}}}},
171172
{"V0Leg/hTPCNsigmaEl", "TPC dE/dx vs. p_{in};p_{in} (GeV/c);n #sigma_{e}^{TPC}", {HistType::kTH2F, {{1000, 0.f, 10.f}, {100, -5.f, +5.f}}}},
@@ -368,7 +369,7 @@ struct PhotonConversionBuilder {
368369
track.tpcChi2NCl(), track.tpcInnerParam(), track.tpcSignal(),
369370
track.tpcNSigmaEl(), track.tpcNSigmaPi(),
370371
track.itsClusterSizes(), track.itsChi2NCl(), track.detectorMap(),
371-
shiftedtrack.getX(), shiftedtrack.getY(), shiftedtrack.getZ(), shiftedtrack.getTgl());
372+
shiftedtrack.getX(), shiftedtrack.getY(), shiftedtrack.getZ(), shiftedtrack.getTgl(), shiftedtrack.getSigma1Pt2());
372373
}
373374

374375
template <bool isMC, class TBCs, class TCollisions, class TTracks, typename TV0>
@@ -647,8 +648,10 @@ struct PhotonConversionBuilder {
647648
registry.fill(HIST("V0Leg/hdEdx_Pin"), leg.tpcInnerParam(), leg.tpcSignal());
648649
registry.fill(HIST("V0Leg/hTPCNsigmaEl"), leg.tpcInnerParam(), leg.tpcNSigmaEl());
649650
} // end of leg loop
650-
registry.fill(HIST("V0Leg/hXZ"), pTrack.getZ(), pTrack.getX());
651-
registry.fill(HIST("V0Leg/hXZ"), nTrack.getZ(), nTrack.getX());
651+
for (auto& leg : {pTrack, nTrack}) {
652+
registry.fill(HIST("V0Leg/hXZ"), leg.getZ(), leg.getX());
653+
registry.fill(HIST("V0Leg/hRelDeltaPt"), leg.getPt(), leg.getPt() * std::sqrt(leg.getSigma1Pt2()));
654+
} // end of leg loop
652655
registry.fill(HIST("V0Leg/hDCAxyz"), posdcaXY, posdcaZ);
653656
registry.fill(HIST("V0Leg/hDCAxyz"), eledcaXY, eledcaZ);
654657

‎PWGEM/PhotonMeson/TableProducer/skimmerGammaConversion.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -190,7 +190,7 @@ struct skimmerGammaConversion {
190190
theTrack.tpcChi2NCl(), theTrack.tpcInnerParam(), theTrack.tpcSignal(),
191191
theTrack.tpcNSigmaEl(), theTrack.tpcNSigmaPi(),
192192
theTrack.itsClusterSizes(), theTrack.itsChi2NCl(), theTrack.detectorMap(),
193-
theTrack.x(), theTrack.y(), theTrack.z(), theTrack.tgl());
193+
theTrack.x(), theTrack.y(), theTrack.z(), theTrack.tgl(), theTrack.c1Pt21Pt2());
194194
}
195195

196196
template <typename TTRACK>

‎PWGEM/PhotonMeson/TableProducer/skimmerPrimaryElectronFromDalitzEE.cxx‎

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -103,6 +103,7 @@ struct skimmerPrimaryElectronFromDalitzEE {
103103
fRegistry.add("Track/hPt", "pT;p_{T} (GeV/c)", kTH1F, {{1000, 0.0f, 10}}, false);
104104
fRegistry.add("Track/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {400, -2.0f, 2.0f}}, false);
105105
fRegistry.add("Track/hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", kTH1F, {{400, -20, 20}}, false);
106+
fRegistry.add("Track/hRelDeltaPt", "pT resolution;p_{T} (GeV/c);#Deltap_{T}/p_{T}", kTH2F, {{1000, 0, 10}, {100, 0, 1}}, false);
106107
fRegistry.add("Track/hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", kTH2F, {{200, -1.0f, 1.0f}, {200, -1.0f, 1.0f}}, false);
107108
fRegistry.add("Track/hDCAxy_Pt", "DCA_{xy} vs. pT;p_{T} (GeV/c);DCA_{xy} (cm)", kTH2F, {{200, 0, 10}, {200, -1, 1}}, false);
108109
fRegistry.add("Track/hDCAz_Pt", "DCA_{z} vs. pT;p_{T} (GeV/c);DCA_{z} (cm)", kTH2F, {{200, 0, 10}, {200, -1, 1}}, false);
@@ -267,7 +268,7 @@ struct skimmerPrimaryElectronFromDalitzEE {
267268
{
268269
if (std::find(stored_trackIds.begin(), stored_trackIds.end(), std::make_pair(collision.globalIndex(), track.globalIndex())) == stored_trackIds.end()) {
269270
emprimaryelectrons(collision.globalIndex(), track.globalIndex(), track.sign(),
270-
track.pt(), track.eta(), track.phi(), track.dcaXY(), track.dcaZ(), track.cYY(), track.cZY(), track.cZZ(),
271+
track.pt(), track.eta(), track.phi(), track.dcaXY(), track.dcaZ(), track.cYY(), track.cZY(), track.cZZ(), track.c1Pt21Pt2(),
271272
track.tpcNClsFindable(), track.tpcNClsFindableMinusFound(), track.tpcNClsFindableMinusCrossedRows(), track.tpcNClsShared(),
272273
track.tpcChi2NCl(), track.tpcInnerParam(),
273274
track.tpcSignal(), track.tpcNSigmaEl(), track.tpcNSigmaPi(),
@@ -277,6 +278,7 @@ struct skimmerPrimaryElectronFromDalitzEE {
277278
fRegistry.fill(HIST("Track/hPt"), track.pt());
278279
fRegistry.fill(HIST("Track/hEtaPhi"), track.phi(), track.eta());
279280
fRegistry.fill(HIST("Track/hQoverPt"), track.sign() / track.pt());
281+
fRegistry.fill(HIST("Track/hRelDeltaPt"), track.pt(), track.sigma1Pt() * track.pt());
280282
fRegistry.fill(HIST("Track/hDCAxyz"), track.dcaXY(), track.dcaZ());
281283
fRegistry.fill(HIST("Track/hDCAxy_Pt"), track.pt(), track.dcaXY());
282284
fRegistry.fill(HIST("Track/hDCAz_Pt"), track.pt(), track.dcaZ());

‎PWGEM/PhotonMeson/Tasks/dalitzEEQC.cxx‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -185,6 +185,7 @@ struct DalitzEEQC {
185185

186186
fRegistry.add("Track/hPt", "pT;p_{T} (GeV/c)", kTH1F, {{1000, 0.0f, 10}}, false);
187187
fRegistry.add("Track/hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", kTH1F, {{400, -20, 20}}, false);
188+
fRegistry.add("Track/hRelDeltaPt", "pT resolution;p_{T} (GeV/c);#Deltap_{T}/p_{T}", kTH2F, {{1000, 0, 10}, {200, 0, 0.2}}, false);
188189
fRegistry.add("Track/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {40, -2.0f, 2.0f}}, false);
189190
fRegistry.add("Track/hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", kTH2F, {{200, -1.0f, 1.0f}, {200, -1.0f, 1.0f}}, false);
190191
fRegistry.add("Track/hDCAxyzSigma", "DCA xy vs. z;DCA_{xy} (#sigma);DCA_{z} (#sigma)", kTH2F, {{200, -10.0f, 10.0f}, {200, -10.0f, 10.0f}}, false);
@@ -314,6 +315,7 @@ struct DalitzEEQC {
314315
{
315316
fRegistry.fill(HIST("Track/hPt"), track.pt());
316317
fRegistry.fill(HIST("Track/hQoverPt"), track.sign() / track.pt());
318+
fRegistry.fill(HIST("Track/hRelDeltaPt"), track.pt(), std::sqrt(track.c1Pt21Pt2()) * track.pt());
317319
fRegistry.fill(HIST("Track/hEtaPhi"), track.phi(), track.eta());
318320
fRegistry.fill(HIST("Track/hDCAxyz"), track.dcaXY(), track.dcaZ());
319321
fRegistry.fill(HIST("Track/hDCAxyzSigma"), track.dcaXY() / sqrt(track.cYY()), track.dcaZ() / sqrt(track.cZZ()));

‎PWGEM/PhotonMeson/Tasks/dalitzEEQCMC.cxx‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -158,6 +158,7 @@ struct DalitzEEQCMC {
158158
// track info
159159
fRegistry.add("Track/hPt", "pT;p_{T} (GeV/c)", kTH1F, {{1000, 0.0f, 10}}, false);
160160
fRegistry.add("Track/hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", kTH1F, {{400, -20, 20}}, false);
161+
fRegistry.add("Track/hRelDeltaPt", "pT resolution;p_{T} (GeV/c);#Deltap_{T}/p_{T}", kTH2F, {{1000, 0, 10}, {200, 0, 0.2}}, false);
161162
fRegistry.add("Track/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {40, -2.0f, 2.0f}}, false);
162163
fRegistry.add("Track/hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", kTH2F, {{200, -1.0f, 1.0f}, {200, -1.0f, 1.0f}}, false);
163164
fRegistry.add("Track/hDCAxyzSigma", "DCA xy vs. z;DCA_{xy} (#sigma);DCA_{z} (#sigma)", kTH2F, {{200, -10.0f, 10.0f}, {200, -10.0f, 10.0f}}, false);
@@ -418,6 +419,7 @@ struct DalitzEEQCMC {
418419
{
419420
fRegistry.fill(HIST("Track/hPt"), track.pt());
420421
fRegistry.fill(HIST("Track/hQoverPt"), track.sign() / track.pt());
422+
fRegistry.fill(HIST("Track/hRelDeltaPt"), track.pt(), std::sqrt(track.c1Pt21Pt2()) * track.pt());
421423
fRegistry.fill(HIST("Track/hEtaPhi"), track.phi(), track.eta());
422424
fRegistry.fill(HIST("Track/hDCAxyz"), track.dcaXY(), track.dcaZ());
423425
fRegistry.fill(HIST("Track/hDCAxyzSigma"), track.dcaXY() / sqrt(track.cYY()), track.dcaZ() / sqrt(track.cZZ()));

‎PWGEM/PhotonMeson/Tasks/pcmQC.cxx‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -162,6 +162,7 @@ struct PCMQC {
162162
// v0leg info
163163
fRegistry.add("V0Leg/hPt", "pT;p_{T,e} (GeV/c)", kTH1F, {{1000, 0.0f, 10}}, false);
164164
fRegistry.add("V0Leg/hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", kTH1F, {{1000, -50, 50}}, false);
165+
fRegistry.add("V0Leg/hRelDeltaPt", "pT resolution;p_{T} (GeV/c);#Deltap_{T}/p_{T}", kTH2F, {{1000, 0, 10}, {200, 0, 0.2}}, false);
165166
fRegistry.add("V0Leg/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{90, 0, 2 * M_PI}, {200, -1.0f, 1.0f}}, false);
166167
fRegistry.add("V0Leg/hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", kTH2F, {{200, -50.0f, 50.0f}, {200, -50.0f, 50.0f}}, false);
167168
fRegistry.add("V0Leg/hNclsTPC", "number of TPC clusters", kTH1F, {{161, -0.5, 160.5}}, false);
@@ -311,6 +312,7 @@ struct PCMQC {
311312
{
312313
fRegistry.fill(HIST("V0Leg/hPt"), leg.pt());
313314
fRegistry.fill(HIST("V0Leg/hQoverPt"), leg.sign() / leg.pt());
315+
fRegistry.fill(HIST("V0Leg/hRelDeltaPt"), leg.pt(), std::sqrt(leg.c1Pt21Pt2()) * leg.pt());
314316
fRegistry.fill(HIST("V0Leg/hEtaPhi"), leg.phi(), leg.eta());
315317
fRegistry.fill(HIST("V0Leg/hDCAxyz"), leg.dcaXY(), leg.dcaZ());
316318
fRegistry.fill(HIST("V0Leg/hNclsITS"), leg.itsNCls());

‎PWGEM/PhotonMeson/Tasks/pcmQCMC.cxx‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -216,6 +216,7 @@ struct PCMQCMC {
216216
// v0leg info
217217
fRegistry.add("V0Leg/primary/hPt", "pT;p_{T,e} (GeV/c)", kTH1F, {{1000, 0.0f, 10}}, false);
218218
fRegistry.add("V0Leg/primary/hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", kTH1F, {{1000, -50, 50}}, false);
219+
fRegistry.add("V0Leg/primary/hRelDeltaPt", "pT resolution;p_{T} (GeV/c);#Deltap_{T}/p_{T}", kTH2F, {{1000, 0, 10}, {200, 0, 0.2}}, false);
219220
fRegistry.add("V0Leg/primary/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{90, 0, 2 * M_PI}, {200, -1.0f, 1.0f}}, false);
220221
fRegistry.add("V0Leg/primary/hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", kTH2F, {{200, -50.0f, 50.0f}, {200, -50.0f, 50.0f}}, false);
221222
fRegistry.add("V0Leg/primary/hNclsTPC", "number of TPC clusters", kTH1F, {{161, -0.5, 160.5}}, false);
@@ -380,6 +381,7 @@ struct PCMQCMC {
380381
fRegistry.fill(HIST("V0Leg/") + HIST(mcphoton_types[mctype]) + HIST("hPt"), leg.pt());
381382
fRegistry.fill(HIST("V0Leg/") + HIST(mcphoton_types[mctype]) + HIST("hQoverPt"), leg.sign() / leg.pt());
382383
fRegistry.fill(HIST("V0Leg/") + HIST(mcphoton_types[mctype]) + HIST("hEtaPhi"), leg.phi(), leg.eta());
384+
fRegistry.fill(HIST("V0Leg/") + HIST(mcphoton_types[mctype]) + HIST("hRelDeltaPt"), leg.pt(), std::sqrt(leg.c1Pt21Pt2()) * leg.pt());
383385
fRegistry.fill(HIST("V0Leg/") + HIST(mcphoton_types[mctype]) + HIST("hDCAxyz"), leg.dcaXY(), leg.dcaZ());
384386
fRegistry.fill(HIST("V0Leg/") + HIST(mcphoton_types[mctype]) + HIST("hNclsITS"), leg.itsNCls());
385387
fRegistry.fill(HIST("V0Leg/") + HIST(mcphoton_types[mctype]) + HIST("hNclsTPC"), leg.tpcNClsFound());

0 commit comments

Comments
 (0)