Skip to content

Commit a03a13b

Browse files
dsekihatdsekihat
authored
[PWGEM/PhotonMeson] add histograms for misidentified ee as V0 photons in MC (#10466)
1 parent cd054ed commit a03a13b

File tree

8 files changed

+68
-48
lines changed

8 files changed

+68
-48
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, track::C1Pt21Pt2,
178+
track::X, track::Y, track::Z, track::Tgl,
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, track::C1Pt21Pt2,
287+
track::Pt, track::Eta, track::Phi, track::DcaXY, track::DcaZ, track::CYY, track::CZY, track::CZZ,
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: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -369,7 +369,7 @@ struct PhotonConversionBuilder {
369369
track.tpcChi2NCl(), track.tpcInnerParam(), track.tpcSignal(),
370370
track.tpcNSigmaEl(), track.tpcNSigmaPi(),
371371
track.itsClusterSizes(), track.itsChi2NCl(), track.detectorMap(),
372-
shiftedtrack.getX(), shiftedtrack.getY(), shiftedtrack.getZ(), shiftedtrack.getTgl(), shiftedtrack.getSigma1Pt2());
372+
shiftedtrack.getX(), shiftedtrack.getY(), shiftedtrack.getZ(), shiftedtrack.getTgl());
373373
}
374374

375375
template <bool isMC, class TBCs, class TCollisions, class TTracks, typename TV0>

‎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(), theTrack.c1Pt21Pt2());
193+
theTrack.x(), theTrack.y(), theTrack.z(), theTrack.tgl());
194194
}
195195

196196
template <typename TTRACK>

‎PWGEM/PhotonMeson/TableProducer/skimmerPrimaryElectronFromDalitzEE.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -268,7 +268,7 @@ struct skimmerPrimaryElectronFromDalitzEE {
268268
{
269269
if (std::find(stored_trackIds.begin(), stored_trackIds.end(), std::make_pair(collision.globalIndex(), track.globalIndex())) == stored_trackIds.end()) {
270270
emprimaryelectrons(collision.globalIndex(), track.globalIndex(), track.sign(),
271-
track.pt(), track.eta(), track.phi(), track.dcaXY(), track.dcaZ(), track.cYY(), track.cZY(), track.cZZ(), track.c1Pt21Pt2(),
271+
track.pt(), track.eta(), track.phi(), track.dcaXY(), track.dcaZ(), track.cYY(), track.cZY(), track.cZZ(),
272272
track.tpcNClsFindable(), track.tpcNClsFindableMinusFound(), track.tpcNClsFindableMinusCrossedRows(), track.tpcNClsShared(),
273273
track.tpcChi2NCl(), track.tpcInnerParam(),
274274
track.tpcSignal(), track.tpcNSigmaEl(), track.tpcNSigmaPi(),

‎PWGEM/PhotonMeson/Tasks/dalitzEEQC.cxx‎

Lines changed: 3 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -142,7 +142,7 @@ struct DalitzEEQC {
142142
if (d_bz_input > -990) {
143143
d_bz = d_bz_input;
144144
o2::parameters::GRPMagField grpmag;
145-
if (fabs(d_bz) > 1e-5) {
145+
if (std::fabs(d_bz) > 1e-5) {
146146
grpmag.setL3Current(30000.f / (d_bz / 5.0f));
147147
}
148148
mRunNumber = collision.runNumber();
@@ -186,9 +186,8 @@ struct DalitzEEQC {
186186

187187
fRegistry.add("Track/hPt", "pT;p_{T} (GeV/c)", kTH1F, {{1000, 0.0f, 10}}, false);
188188
fRegistry.add("Track/hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", kTH1F, {{400, -20, 20}}, false);
189-
fRegistry.add("Track/hRelDeltaPt", "pT resolution;p_{T} (GeV/c);#Deltap_{T}/p_{T}", kTH2F, {{1000, 0, 10}, {200, 0, 0.2}}, false);
190189
fRegistry.add("Track/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {40, -2.0f, 2.0f}}, false);
191-
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);
190+
fRegistry.add("Track/hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", kTH2F, {{200, -0.1f, 0.1f}, {200, -0.1f, 0.1f}}, false);
192191
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);
193192
fRegistry.add("Track/hDCAxyRes_Pt", "DCA_{xy} resolution vs. pT;p_{T} (GeV/c);DCA_{xy} resolution (#mum)", kTH2F, {{200, 0, 10}, {500, 0., 500}}, false);
194193
fRegistry.add("Track/hDCAzRes_Pt", "DCA_{z} resolution vs. pT;p_{T} (GeV/c);DCA_{z} resolution (#mum)", kTH2F, {{200, 0, 10}, {500, 0., 500}}, false);
@@ -277,7 +276,7 @@ struct DalitzEEQC {
277276
ROOT::Math::PtEtaPhiMVector v1(t1.pt(), t1.eta(), t1.phi(), o2::constants::physics::MassElectron);
278277
ROOT::Math::PtEtaPhiMVector v2(t2.pt(), t2.eta(), t2.phi(), o2::constants::physics::MassElectron);
279278
ROOT::Math::PtEtaPhiMVector v12 = v1 + v2;
280-
if (abs(v12.Rapidity()) > maxY) {
279+
if (std::fabs(v12.Rapidity()) > maxY) {
281280
return false;
282281
}
283282
float phiv = o2::aod::pwgem::dilepton::utils::pairutil::getPhivPair(t1.px(), t1.py(), t1.pz(), t2.px(), t2.py(), t2.pz(), t1.sign(), t2.sign(), d_bz);
@@ -317,7 +316,6 @@ struct DalitzEEQC {
317316
{
318317
fRegistry.fill(HIST("Track/hPt"), track.pt());
319318
fRegistry.fill(HIST("Track/hQoverPt"), track.sign() / track.pt());
320-
fRegistry.fill(HIST("Track/hRelDeltaPt"), track.pt(), std::sqrt(track.c1Pt21Pt2()) * track.pt());
321319
fRegistry.fill(HIST("Track/hEtaPhi"), track.phi(), track.eta());
322320
fRegistry.fill(HIST("Track/hDCAxyz"), track.dcaXY(), track.dcaZ());
323321
fRegistry.fill(HIST("Track/hDCAxyzSigma"), track.dcaXY() / sqrt(track.cYY()), track.dcaZ() / sqrt(track.cZZ()));

‎PWGEM/PhotonMeson/Tasks/dalitzEEQCMC.cxx‎

Lines changed: 8 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -159,9 +159,8 @@ struct DalitzEEQCMC {
159159
// track info
160160
fRegistry.add("Track/hPt", "pT;p_{T} (GeV/c)", kTH1F, {{1000, 0.0f, 10}}, false);
161161
fRegistry.add("Track/hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", kTH1F, {{400, -20, 20}}, false);
162-
fRegistry.add("Track/hRelDeltaPt", "pT resolution;p_{T} (GeV/c);#Deltap_{T}/p_{T}", kTH2F, {{1000, 0, 10}, {200, 0, 0.2}}, false);
163162
fRegistry.add("Track/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {40, -2.0f, 2.0f}}, false);
164-
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);
163+
fRegistry.add("Track/hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", kTH2F, {{200, -0.1f, 0.1f}, {200, -0.1f, 0.1f}}, false);
165164
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);
166165
fRegistry.add("Track/hDCAxyRes_Pt", "DCA_{xy} resolution vs. pT;p_{T} (GeV/c);DCA_{xy} resolution (#mum)", kTH2F, {{200, 0, 10}, {500, 0., 500}}, false);
167166
fRegistry.add("Track/hDCAzRes_Pt", "DCA_{z} resolution vs. pT;p_{T} (GeV/c);DCA_{z} resolution (#mum)", kTH2F, {{200, 0, 10}, {500, 0., 500}}, false);
@@ -213,7 +212,7 @@ struct DalitzEEQCMC {
213212
if (d_bz_input > -990) {
214213
d_bz = d_bz_input;
215214
o2::parameters::GRPMagField grpmag;
216-
if (fabs(d_bz) > 1e-5) {
215+
if (std::fabs(d_bz) > 1e-5) {
217216
grpmag.setL3Current(30000.f / (d_bz / 5.0f));
218217
}
219218
mRunNumber = collision.runNumber();
@@ -307,7 +306,7 @@ struct DalitzEEQCMC {
307306
template <typename T>
308307
bool isInAcceptance(T const& t1)
309308
{
310-
if ((mctrackcuts.min_mcPt < t1.pt() && t1.pt() < mctrackcuts.max_mcPt) && abs(t1.eta()) < mctrackcuts.max_mcEta) {
309+
if ((mctrackcuts.min_mcPt < t1.pt() && t1.pt() < mctrackcuts.max_mcPt) && std::fabs(t1.eta()) < mctrackcuts.max_mcEta) {
311310
return true;
312311
} else {
313312
return false;
@@ -335,7 +334,7 @@ struct DalitzEEQCMC {
335334
ROOT::Math::PtEtaPhiMVector v1(t1.pt(), t1.eta(), t1.phi(), o2::constants::physics::MassElectron);
336335
ROOT::Math::PtEtaPhiMVector v2(t2.pt(), t2.eta(), t2.phi(), o2::constants::physics::MassElectron);
337336
ROOT::Math::PtEtaPhiMVector v12 = v1 + v2;
338-
if (abs(v12.Rapidity()) > maxY) {
337+
if (std::fabs(v12.Rapidity()) > maxY) {
339338
return false;
340339
}
341340
float phiv = o2::aod::pwgem::dilepton::utils::pairutil::getPhivPair(t1.px(), t1.py(), t1.pz(), t2.px(), t2.py(), t2.pz(), t1.sign(), t2.sign(), d_bz);
@@ -348,7 +347,7 @@ struct DalitzEEQCMC {
348347

349348
if (mcmother.isPhysicalPrimary() || mcmother.producedByGenerator()) {
350349
if ((t1mc.isPhysicalPrimary() || t1mc.producedByGenerator()) && (t2mc.isPhysicalPrimary() || t2mc.producedByGenerator())) {
351-
switch (abs(mcmother.pdgCode())) {
350+
switch (std::abs(mcmother.pdgCode())) {
352351
case 111:
353352
fRegistry.fill(HIST("Pair/sm/Pi0/hMvsPt"), v12.M(), v12.Pt());
354353
fRegistry.fill(HIST("Pair/sm/Pi0/hMvsPhiV"), phiv, v12.M());
@@ -377,7 +376,7 @@ struct DalitzEEQCMC {
377376
break;
378377
}
379378
} else if (!(t1mc.isPhysicalPrimary() || t1mc.producedByGenerator()) && !(t2mc.isPhysicalPrimary() || t2mc.producedByGenerator())) {
380-
switch (abs(mcmother.pdgCode())) {
379+
switch (std::abs(mcmother.pdgCode())) {
381380
case 22:
382381
fRegistry.fill(HIST("Pair/sm/Photon/hMvsPt"), v12.M(), v12.Pt());
383382
fRegistry.fill(HIST("Pair/sm/Photon/hMvsPhiV"), phiv, v12.M());
@@ -421,7 +420,6 @@ struct DalitzEEQCMC {
421420
{
422421
fRegistry.fill(HIST("Track/hPt"), track.pt());
423422
fRegistry.fill(HIST("Track/hQoverPt"), track.sign() / track.pt());
424-
fRegistry.fill(HIST("Track/hRelDeltaPt"), track.pt(), std::sqrt(track.c1Pt21Pt2()) * track.pt());
425423
fRegistry.fill(HIST("Track/hEtaPhi"), track.phi(), track.eta());
426424
fRegistry.fill(HIST("Track/hDCAxyz"), track.dcaXY(), track.dcaZ());
427425
fRegistry.fill(HIST("Track/hDCAxyzSigma"), track.dcaXY() / sqrt(track.cYY()), track.dcaZ() / sqrt(track.cZZ()));
@@ -542,15 +540,15 @@ struct DalitzEEQCMC {
542540
ROOT::Math::PtEtaPhiMVector v2(t2.pt(), t2.eta(), t2.phi(), o2::constants::physics::MassElectron);
543541
ROOT::Math::PtEtaPhiMVector v12 = v1 + v2;
544542

545-
if (abs(v12.Rapidity()) > maxY) {
543+
if (std::fabs(v12.Rapidity()) > maxY) {
546544
continue;
547545
}
548546

549547
if (mother_id > -1) {
550548
auto mcmother = mcparticles.iteratorAt(mother_id);
551549
if (mcmother.isPhysicalPrimary() || mcmother.producedByGenerator()) {
552550

553-
switch (abs(mcmother.pdgCode())) {
551+
switch (std::abs(mcmother.pdgCode())) {
554552
case 111:
555553
fRegistry.fill(HIST("Generated/sm/Pi0/hMvsPt"), v12.M(), v12.Pt());
556554
break;

‎PWGEM/PhotonMeson/Tasks/pcmQC.cxx‎

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -162,7 +162,6 @@ 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);
166165
fRegistry.add("V0Leg/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{90, 0, 2 * M_PI}, {200, -1.0f, 1.0f}}, false);
167166
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);
168167
fRegistry.add("V0Leg/hNclsTPC", "number of TPC clusters", kTH1F, {{161, -0.5, 160.5}}, false);
@@ -312,7 +311,6 @@ struct PCMQC {
312311
{
313312
fRegistry.fill(HIST("V0Leg/hPt"), leg.pt());
314313
fRegistry.fill(HIST("V0Leg/hQoverPt"), leg.sign() / leg.pt());
315-
fRegistry.fill(HIST("V0Leg/hRelDeltaPt"), leg.pt(), std::sqrt(leg.c1Pt21Pt2()) * leg.pt());
316314
fRegistry.fill(HIST("V0Leg/hEtaPhi"), leg.phi(), leg.eta());
317315
fRegistry.fill(HIST("V0Leg/hDCAxyz"), leg.dcaXY(), leg.dcaZ());
318316
fRegistry.fill(HIST("V0Leg/hNclsITS"), leg.itsNCls());

0 commit comments

Comments
 (0)