Skip to content

Commit a694bb3

Browse files
dsekihatdsekihat
committed
PWGEM/Dilepton: update for ITSsa
1 parent 61234b3 commit a694bb3

2 files changed

Lines changed: 252 additions & 106 deletions

File tree

‎PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx‎

Lines changed: 152 additions & 68 deletions
Original file line numberDiff line numberDiff line change
@@ -50,7 +50,7 @@ using namespace o2::constants::physics;
5050
using MyCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::EMEvSels>;
5151
using MyCollisionsWithSWT = soa::Join<MyCollisions, aod::EMSWTriggerInfosTMP>;
5252

53-
using MyTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TracksCov,
53+
using MyTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksCov,
5454
aod::pidTPCFullEl, aod::pidTPCFullMu, aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr,
5555
aod::pidTOFFullEl, aod::pidTOFFullMu, aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr, aod::pidTOFbeta>;
5656
using MyTrack = MyTracks::iterator;
@@ -81,7 +81,8 @@ struct skimmerPrimaryElectron {
8181
Configurable<int> min_ncluster_itsib{"min_ncluster_itsib", 1, "min ncluster itsib"};
8282
Configurable<float> maxchi2tpc{"maxchi2tpc", 5.0, "max. chi2/NclsTPC"};
8383
Configurable<float> maxchi2its{"maxchi2its", 6.0, "max. chi2/NclsITS"};
84-
Configurable<float> minpt{"minpt", 0.15, "min pt for track"};
84+
Configurable<float> minpt_itstpc{"minpt_itstpc", 0.1, "min pt for ITS-TPC track"};
85+
Configurable<float> minpt_itssa{"minpt_itssa", 0.05, "min pt for ITSsa track"};
8586
Configurable<float> maxeta{"maxeta", 0.9, "eta acceptance"};
8687
Configurable<float> dca_xy_max{"dca_xy_max", 1.0, "max DCAxy in cm"};
8788
Configurable<float> dca_z_max{"dca_z_max", 1.0, "max DCAz in cm"};
@@ -99,6 +100,9 @@ struct skimmerPrimaryElectron {
99100
Configurable<float> min_pin_for_pion_rejection{"min_pin_for_pion_rejection", 0.0, "pion rejection is applied above this pin"}; // this is used only in TOFreq
100101
Configurable<float> max_pin_for_pion_rejection{"max_pin_for_pion_rejection", 0.5, "pion rejection is applied below this pin"};
101102
Configurable<float> max_frac_shared_clusters_tpc{"max_frac_shared_clusters_tpc", 999.f, "max fraction of shared clusters in TPC"};
103+
Configurable<bool> includeITSsa{"includeITSsa", false, "Flag to include ITSsa tracks"};
104+
Configurable<float> maxpt_itssa{"maxpt_itssa", 0.15, "max pt for ITSsa track"};
105+
Configurable<float> maxMeanITSClusterSize{"maxMeanITSClusterSize", 16, "max <ITS cluster size> x cos(lambda)"};
102106

103107
// configuration for PID ML
104108
Configurable<bool> usePIDML{"usePIDML", false, "Flag to use PID ML"};
@@ -166,7 +170,9 @@ struct skimmerPrimaryElectron {
166170
fRegistry.add("Track/hTOFNsigmaPi", "TOF n sigma pi;p_{in} (GeV/c);n #sigma_{#pi}^{TOF}", kTH2F, {{1000, 0, 10}, {100, -5, +5}}, false);
167171
fRegistry.add("Track/hTOFNsigmaKa", "TOF n sigma ka;p_{in} (GeV/c);n #sigma_{K}^{TOF}", kTH2F, {{1000, 0, 10}, {100, -5, +5}}, false);
168172
fRegistry.add("Track/hTOFNsigmaPr", "TOF n sigma pr;p_{in} (GeV/c);n #sigma_{p}^{TOF}", kTH2F, {{1000, 0, 10}, {100, -5, +5}}, false);
169-
fRegistry.add("Track/hMeanClusterSizeITS", "mean cluster size ITS;p_{pv} (GeV/c);<cluster size> on ITS #times cos(#lambda)", kTH2F, {{1000, 0, 10}, {150, 0, 15}}, false);
173+
fRegistry.add("Track/hMeanClusterSizeITS", "mean cluster size ITS;p_{pv} (GeV/c);<ITS cluster size> #times cos(#lambda)", kTH2F, {{1000, 0, 10}, {150, 0, 15}}, false);
174+
fRegistry.add("Track/hMeanClusterSizeITSib", "mean cluster size ITSib;p_{pv} (GeV/c);<ITSib cluster size> #times cos(#lambda)", kTH2F, {{1000, 0, 10}, {150, 0, 15}}, false);
175+
fRegistry.add("Track/hMeanClusterSizeITSob", "mean cluster size ITSob;p_{pv} (GeV/c);<ITSob cluster size> #times cos(#lambda)", kTH2F, {{1000, 0, 10}, {150, 0, 15}}, false);
170176
fRegistry.add("Track/hITSNsigmaEl", "ITS n sigma el;p_{pv} (GeV/c);n #sigma_{e}^{ITS}", kTH2F, {{1000, 0, 10}, {100, -5, +5}}, false);
171177
fRegistry.add("Track/hITSNsigmaMu", "ITS n sigma mu;p_{pv} (GeV/c);n #sigma_{#mu}^{ITS}", kTH2F, {{1000, 0, 10}, {100, -5, +5}}, false);
172178
fRegistry.add("Track/hITSNsigmaPi", "ITS n sigma pi;p_{pv} (GeV/c);n #sigma_{#pi}^{ITS}", kTH2F, {{1000, 0, 10}, {100, -5, +5}}, false);
@@ -270,42 +276,44 @@ struct skimmerPrimaryElectron {
270276
return false;
271277
}
272278

273-
if (track.tpcChi2NCl() > maxchi2tpc) {
279+
if (!track.hasITS()) {
274280
return false;
275281
}
276282

277283
if (track.itsChi2NCl() > maxchi2its) {
278284
return false;
279285
}
280-
281-
if (!track.hasITS() || !track.hasTPC()) {
282-
return false;
283-
}
284286
if (track.itsNCls() < min_ncluster_its) {
285287
return false;
286288
}
287289
if (track.itsNClsInnerBarrel() < min_ncluster_itsib) {
288290
return false;
289291
}
290292

291-
if (track.tpcNClsFound() < min_ncluster_tpc) {
293+
if (!includeITSsa && (!track.hasITS() || !track.hasTPC())) {
292294
return false;
293295
}
294296

295-
if (track.tpcNClsCrossedRows() < mincrossedrows) {
296-
return false;
297-
}
297+
if (track.hasTPC()) {
298+
if (track.tpcChi2NCl() > maxchi2tpc) {
299+
return false;
300+
}
298301

299-
if (track.tpcCrossedRowsOverFindableCls() < min_tpc_cr_findable_ratio) {
300-
return false;
301-
}
302+
if (track.tpcNClsFound() < min_ncluster_tpc) {
303+
return false;
304+
}
302305

303-
if (track.tpcFractionSharedCls() > max_frac_shared_clusters_tpc) {
304-
return false;
305-
}
306+
if (track.tpcNClsCrossedRows() < mincrossedrows) {
307+
return false;
308+
}
306309

307-
if (track.hasTOF() && (maxTOFNsigmaEl < std::fabs(track.tofNSigmaEl()))) {
308-
return false;
310+
if (track.tpcCrossedRowsOverFindableCls() < min_tpc_cr_findable_ratio) {
311+
return false;
312+
}
313+
314+
if (track.tpcFractionSharedCls() > max_frac_shared_clusters_tpc) {
315+
return false;
316+
}
309317
}
310318

311319
o2::dataformats::DCA mDcaInfoCov;
@@ -324,10 +332,6 @@ struct skimmerPrimaryElectron {
324332
return false;
325333
}
326334

327-
if (track_par_cov_recalc.getPt() < minpt || std::fabs(track_par_cov_recalc.getEta()) > maxeta) {
328-
return false;
329-
}
330-
331335
float dca_3d = 999.f;
332336
float det = track_par_cov_recalc.getSigmaY2() * track_par_cov_recalc.getSigmaZ2() - track_par_cov_recalc.getSigmaZY() * track_par_cov_recalc.getSigmaZY();
333337
if (det < 0) {
@@ -340,12 +344,43 @@ struct skimmerPrimaryElectron {
340344
return false;
341345
}
342346

347+
if (std::fabs(track_par_cov_recalc.getEta()) > maxeta) {
348+
return false;
349+
}
350+
351+
if ((track.hasITS() && track.hasTPC()) && track_par_cov_recalc.getPt() < minpt_itstpc) {
352+
return false;
353+
}
354+
355+
if ((track.hasITS() && !track.hasTPC() && !track.hasTOF() && !track.hasTRD()) && (track_par_cov_recalc.getPt() < minpt_itssa || maxpt_itssa < track_par_cov_recalc.getPt())) {
356+
return false;
357+
}
358+
359+
if (track.hasITS() && !track.hasTPC() && !track.hasTRD() && !track.hasTOF()) {
360+
int total_cluster_size = 0, nl = 0;
361+
for (unsigned int layer = 0; layer < 7; layer++) {
362+
int cluster_size_per_layer = track.itsClsSizeInLayer(layer);
363+
if (cluster_size_per_layer > 0) {
364+
nl++;
365+
}
366+
total_cluster_size += cluster_size_per_layer;
367+
}
368+
369+
if (maxMeanITSClusterSize < static_cast<float>(total_cluster_size) / static_cast<float>(nl) * std::cos(std::atan(track_par_cov_recalc.getTgl()))) {
370+
return false;
371+
}
372+
}
373+
343374
return true;
344375
}
345376

346377
template <typename TCollision, typename TTrack>
347378
bool isElectron(TCollision const& collision, TTrack const& track)
348379
{
380+
if (includeITSsa && (track.hasITS() && !track.hasTPC() && !track.hasTRD() && !track.hasTOF())) {
381+
return true;
382+
}
383+
349384
if (usePIDML) {
350385
if (track.tpcNSigmaEl() < minTPCNsigmaEl || maxTPCNsigmaEl < track.tpcNSigmaEl()) {
351386
return false;
@@ -451,16 +486,34 @@ struct skimmerPrimaryElectron {
451486
stored_trackIds.emplace_back(std::pair<int, int>{collision.globalIndex(), track.globalIndex()});
452487

453488
if (fillQAHistogram) {
454-
uint32_t itsClusterSizes = track.itsClusterSizes();
489+
// uint32_t itsClusterSizes = track.itsClusterSizes();
455490
int total_cluster_size = 0, nl = 0;
456491
for (unsigned int layer = 0; layer < 7; layer++) {
457-
int cluster_size_per_layer = (itsClusterSizes >> (layer * 4)) & 0xf;
492+
int cluster_size_per_layer = track.itsClsSizeInLayer(layer);
458493
if (cluster_size_per_layer > 0) {
459494
nl++;
460495
}
461496
total_cluster_size += cluster_size_per_layer;
462497
}
463498

499+
int total_cluster_size_ib = 0, nl_ib = 0;
500+
for (unsigned int layer = 0; layer < 3; layer++) {
501+
int cluster_size_per_layer = track.itsClsSizeInLayer(layer);
502+
if (cluster_size_per_layer > 0) {
503+
nl_ib++;
504+
}
505+
total_cluster_size_ib += cluster_size_per_layer;
506+
}
507+
508+
int total_cluster_size_ob = 0, nl_ob = 0;
509+
for (unsigned int layer = 3; layer < 7; layer++) {
510+
int cluster_size_per_layer = track.itsClsSizeInLayer(layer);
511+
if (cluster_size_per_layer > 0) {
512+
nl_ob++;
513+
}
514+
total_cluster_size_ob += cluster_size_per_layer;
515+
}
516+
464517
fRegistry.fill(HIST("Track/hPt"), pt_recalc);
465518
fRegistry.fill(HIST("Track/hQoverPt"), track.sign() / pt_recalc);
466519
fRegistry.fill(HIST("Track/hEtaPhi"), phi_recalc, eta_recalc);
@@ -484,26 +537,27 @@ struct skimmerPrimaryElectron {
484537
fRegistry.fill(HIST("Track/hTPCNsigmaPi"), track.tpcInnerParam(), track.tpcNSigmaPi());
485538
fRegistry.fill(HIST("Track/hTPCNsigmaKa"), track.tpcInnerParam(), track.tpcNSigmaKa());
486539
fRegistry.fill(HIST("Track/hTPCNsigmaPr"), track.tpcInnerParam(), track.tpcNSigmaPr());
487-
fRegistry.fill(HIST("Track/hTOFbeta"), track.p(), track.beta());
540+
fRegistry.fill(HIST("Track/hTOFbeta"), track_par_cov_recalc.getP(), track.beta());
488541
fRegistry.fill(HIST("Track/hTOFNsigmaEl"), track.tpcInnerParam(), track.tofNSigmaEl());
489542
fRegistry.fill(HIST("Track/hTOFNsigmaMu"), track.tpcInnerParam(), track.tofNSigmaMu());
490543
fRegistry.fill(HIST("Track/hTOFNsigmaPi"), track.tpcInnerParam(), track.tofNSigmaPi());
491544
fRegistry.fill(HIST("Track/hTOFNsigmaKa"), track.tpcInnerParam(), track.tofNSigmaKa());
492545
fRegistry.fill(HIST("Track/hTOFNsigmaPr"), track.tpcInnerParam(), track.tofNSigmaPr());
493-
fRegistry.fill(HIST("Track/hMeanClusterSizeITS"), track.p(), static_cast<float>(total_cluster_size) / static_cast<float>(nl) * std::cos(std::atan(track.tgl())));
494-
fRegistry.fill(HIST("Track/hITSNsigmaEl"), track.p(), track.itsNSigmaEl());
495-
fRegistry.fill(HIST("Track/hITSNsigmaMu"), track.p(), track.itsNSigmaMu());
496-
fRegistry.fill(HIST("Track/hITSNsigmaPi"), track.p(), track.itsNSigmaPi());
497-
fRegistry.fill(HIST("Track/hITSNsigmaKa"), track.p(), track.itsNSigmaKa());
498-
fRegistry.fill(HIST("Track/hITSNsigmaPr"), track.p(), track.itsNSigmaPr());
546+
fRegistry.fill(HIST("Track/hMeanClusterSizeITS"), track_par_cov_recalc.getP(), static_cast<float>(total_cluster_size) / static_cast<float>(nl) * std::cos(std::atan(track_par_cov_recalc.getTgl())));
547+
fRegistry.fill(HIST("Track/hMeanClusterSizeITSib"), track_par_cov_recalc.getP(), static_cast<float>(total_cluster_size_ib) / static_cast<float>(nl_ib) * std::cos(std::atan(track_par_cov_recalc.getTgl())));
548+
fRegistry.fill(HIST("Track/hMeanClusterSizeITSob"), track_par_cov_recalc.getP(), static_cast<float>(total_cluster_size_ob) / static_cast<float>(nl_ob) * std::cos(std::atan(track_par_cov_recalc.getTgl())));
549+
fRegistry.fill(HIST("Track/hITSNsigmaEl"), track_par_cov_recalc.getP(), track.itsNSigmaEl());
550+
fRegistry.fill(HIST("Track/hITSNsigmaMu"), track_par_cov_recalc.getP(), track.itsNSigmaMu());
551+
fRegistry.fill(HIST("Track/hITSNsigmaPi"), track_par_cov_recalc.getP(), track.itsNSigmaPi());
552+
fRegistry.fill(HIST("Track/hITSNsigmaKa"), track_par_cov_recalc.getP(), track.itsNSigmaKa());
553+
fRegistry.fill(HIST("Track/hITSNsigmaPr"), track_par_cov_recalc.getP(), track.itsNSigmaPr());
499554
}
500555
}
501556
}
502557

503558
Preslice<aod::TrackAssoc> trackIndicesPerCollision = aod::track_association::collisionId;
504559
std::vector<std::pair<int, int>> stored_trackIds;
505-
Filter trackFilter = o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& o2::aod::track::tpcChi2NCl < maxchi2tpc&& o2::aod::track::itsChi2NCl < maxchi2its&& ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::ITS) == true && ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::TPC) == true;
506-
Filter pidFilter = minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl;
560+
Filter trackFilter = o2::aod::track::pt > minpt_itssa&& nabs(o2::aod::track::eta) < maxeta&& o2::aod::track::itsChi2NCl < maxchi2its&& ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::ITS) == true;
507561
using MyFilteredTracks = soa::Filtered<MyTracks>;
508562

509563
Partition<MyFilteredTracks> posTracks = o2::aod::track::signed1Pt > 0.f;
@@ -724,7 +778,8 @@ struct prefilterPrimaryElectron {
724778
Configurable<bool> fillQAHistogram{"fillQAHistogram", false, "flag to fill QA histograms"};
725779
Configurable<float> max_dcaxy{"max_dcaxy", 0.3, "DCAxy To PV for loose track sample"};
726780
Configurable<float> max_dcaz{"max_dcaz", 0.3, "DCAz To PV for loose track sample"};
727-
Configurable<float> minpt{"minpt", 0.1, "min pt for track for loose track sample"};
781+
Configurable<float> minpt_itstpc{"minpt_itstpc", 0.1, "min pt for ITS-TPC track"};
782+
Configurable<float> minpt_itssa{"minpt_itssa", 0.05, "min pt for ITSsa track"};
728783
Configurable<float> maxeta{"maxeta", 1.2, "eta acceptance for loose track sample"};
729784
Configurable<int> min_ncluster_tpc{"min_ncluster_tpc", 0, "min ncluster tpc"};
730785
Configurable<int> mincrossedrows{"mincrossedrows", 70, "min crossed rows"};
@@ -738,6 +793,9 @@ struct prefilterPrimaryElectron {
738793
Configurable<float> maxTPCNsigmaEl{"maxTPCNsigmaEl", 3.0, "max. TPC n sigma for electron inclusion"};
739794
Configurable<float> slope{"slope", 0.0185, "slope for m vs. phiv"};
740795
Configurable<float> intercept{"intercept", -0.0280, "intercept for m vs. phiv"};
796+
Configurable<bool> includeITSsa{"includeITSsa", false, "Flag to include ITSsa tracks"};
797+
Configurable<float> maxpt_itssa{"maxpt_itssa", 0.15, "mix pt for ITSsa track"};
798+
Configurable<float> maxMeanITSClusterSize{"maxMeanITSClusterSize", 16, "max <ITS cluster size> x cos(lambda)"};
741799

742800
Configurable<std::vector<float>> max_mee_vec{"max_mee_vec", std::vector<float>{0.06, 0.08, 0.10}, "vector fo max mee for prefilter in ULS. Please sort this by increasing order."}; // currently, 3 thoresholds are allowed.
743801

@@ -852,26 +910,29 @@ struct prefilterPrimaryElectron {
852910
return false;
853911
}
854912

855-
if (!track.hasTPC()) {
856-
return false;
857-
}
858-
if (track.tpcNSigmaEl() < minTPCNsigmaEl || maxTPCNsigmaEl < track.tpcNSigmaEl()) {
859-
return false;
860-
}
861-
if (track.tpcNClsFound() < min_ncluster_tpc) {
862-
return false;
863-
}
864-
if (track.tpcNClsCrossedRows() < mincrossedrows) {
865-
return false;
866-
}
867-
if (track.tpcCrossedRowsOverFindableCls() < min_tpc_cr_findable_ratio) {
868-
return false;
869-
}
870-
if (track.tpcFractionSharedCls() > max_frac_shared_clusters_tpc) {
913+
if (!includeITSsa && (!track.hasITS() || !track.hasTPC())) {
871914
return false;
872915
}
873-
if (track.tpcChi2NCl() > maxchi2its) {
874-
return false;
916+
917+
if (track.hasTPC()) {
918+
if (track.tpcNSigmaEl() < minTPCNsigmaEl || maxTPCNsigmaEl < track.tpcNSigmaEl()) {
919+
return false;
920+
}
921+
if (track.tpcNClsFound() < min_ncluster_tpc) {
922+
return false;
923+
}
924+
if (track.tpcNClsCrossedRows() < mincrossedrows) {
925+
return false;
926+
}
927+
if (track.tpcCrossedRowsOverFindableCls() < min_tpc_cr_findable_ratio) {
928+
return false;
929+
}
930+
if (track.tpcFractionSharedCls() > max_frac_shared_clusters_tpc) {
931+
return false;
932+
}
933+
if (track.tpcChi2NCl() > maxchi2tpc) {
934+
return false;
935+
}
875936
}
876937

877938
o2::dataformats::DCA mDcaInfoCov;
@@ -888,10 +949,33 @@ struct prefilterPrimaryElectron {
888949
return false;
889950
}
890951

891-
if (track_par_cov_recalc.getPt() < minpt || std::fabs(track_par_cov_recalc.getEta()) > maxeta) {
952+
if (std::fabs(track_par_cov_recalc.getEta()) > maxeta) {
892953
return false;
893954
}
894955

956+
if ((track.hasITS() && track.hasTPC()) && track_par_cov_recalc.getPt() < minpt_itstpc) {
957+
return false;
958+
}
959+
960+
if ((track.hasITS() && !track.hasTPC() && !track.hasTOF() && !track.hasTRD()) && (track_par_cov_recalc.getPt() < minpt_itssa || maxpt_itssa < track_par_cov_recalc.getPt())) {
961+
return false;
962+
}
963+
964+
if (track.hasITS() && !track.hasTPC() && !track.hasTOF() && !track.hasTRD()) {
965+
int total_cluster_size = 0, nl = 0;
966+
for (unsigned int layer = 0; layer < 7; layer++) {
967+
int cluster_size_per_layer = track.itsClsSizeInLayer(layer);
968+
if (cluster_size_per_layer > 0) {
969+
nl++;
970+
}
971+
total_cluster_size += cluster_size_per_layer;
972+
}
973+
974+
if (maxMeanITSClusterSize < static_cast<float>(total_cluster_size) / static_cast<float>(nl) * std::cos(std::atan(track_par_cov_recalc.getTgl()))) {
975+
return false;
976+
}
977+
}
978+
895979
return true;
896980
}
897981

@@ -936,7 +1020,7 @@ struct prefilterPrimaryElectron {
9361020

9371021
Preslice<aod::TrackAssoc> trackIndicesPerCollision = aod::track_association::collisionId;
9381022

939-
Filter trackFilter = o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::ITS) == true && ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::TPC) == true;
1023+
Filter trackFilter = o2::aod::track::pt > minpt_itssa&& nabs(o2::aod::track::eta) < maxeta&& o2::aod::track::itsChi2NCl < maxchi2its&& ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::ITS) == true;
9401024
using MyFilteredTracks = soa::Filtered<MyTracks>;
9411025
Partition<MyFilteredTracks> posTracks = o2::aod::track::signed1Pt > 0.f;
9421026
Partition<MyFilteredTracks> negTracks = o2::aod::track::signed1Pt < 0.f;
@@ -981,9 +1065,9 @@ struct prefilterPrimaryElectron {
9811065
}
9821066

9831067
for (const auto& ele : negTracks_per_coll) {
984-
if (!checkTrack(collision, ele)) {
985-
continue;
986-
}
1068+
// if (!checkTrack(collision, ele)) {
1069+
// continue;
1070+
// }
9871071
o2::dataformats::DCA mDcaInfoCov;
9881072
mDcaInfoCov.set(999, 999, 999, 999, 999);
9891073
std::array<float, 3> pVec_recalc = {0, 0, 0}; // px, py, pz
@@ -1026,9 +1110,9 @@ struct prefilterPrimaryElectron {
10261110
} // end of loose electron loop
10271111

10281112
for (const auto& pos : posTracks_per_coll) {
1029-
if (!checkTrack(collision, pos)) { // track cut is applied to loose sample
1030-
continue;
1031-
}
1113+
// if (!checkTrack(collision, pos)) { // track cut is applied to loose sample
1114+
// continue;
1115+
// }
10321116
o2::dataformats::DCA mDcaInfoCov;
10331117
mDcaInfoCov.set(999, 999, 999, 999, 999);
10341118
std::array<float, 3> pVec_recalc = {0, 0, 0}; // px, py, pz
@@ -1069,9 +1153,9 @@ struct prefilterPrimaryElectron {
10691153
} // end of loose positon loop
10701154

10711155
for (const auto& pos : posTracks_per_coll) {
1072-
if (!checkTrack(collision, pos)) { // track cut is applied to loose sample
1073-
continue;
1074-
}
1156+
// if (!checkTrack(collision, pos)) { // track cut is applied to loose sample
1157+
// continue;
1158+
// }
10751159
o2::dataformats::DCA mDcaInfoCov;
10761160
mDcaInfoCov.set(999, 999, 999, 999, 999);
10771161
std::array<float, 3> pVec_recalc = {0, 0, 0}; // px, py, pz
@@ -1098,9 +1182,9 @@ struct prefilterPrimaryElectron {
10981182
} // end of loose positon loop
10991183

11001184
for (const auto& ele : negTracks_per_coll) {
1101-
if (!checkTrack(collision, ele)) {
1102-
continue;
1103-
}
1185+
// if (!checkTrack(collision, ele)) {
1186+
// continue;
1187+
// }
11041188
o2::dataformats::DCA mDcaInfoCov;
11051189
mDcaInfoCov.set(999, 999, 999, 999, 999);
11061190
std::array<float, 3> pVec_recalc = {0, 0, 0}; // px, py, pz

0 commit comments

Comments
 (0)