Skip to content

Commit 99ca939

Browse files
dsekihatdsekihat
authored
[PWGEM/Dilepton] update matchingMFT.cxx (#11004)
1 parent bccdb6a commit 99ca939

File tree

1 file changed

+40
-52
lines changed

1 file changed

+40
-52
lines changed

‎PWGEM/Dilepton/Tasks/matchingMFT.cxx‎

Lines changed: 40 additions & 52 deletions
Original file line numberDiff line numberDiff line change
@@ -76,7 +76,8 @@ struct matchingMFT {
7676
Configurable<bool> requireTrueAssociation{"requireTrueAssociation", false, "flag to require true mc collision association"};
7777
Configurable<float> maxDEta{"maxDEta", 1e+10f, "max. deta between MFT-MCH-MID and MCH-MID"};
7878
Configurable<float> maxDPhi{"maxDPhi", 1e+10f, "max. dphi between MFT-MCH-MID and MCH-MID"};
79-
Configurable<std::vector<int>> requireMFTDisks{"requireMFTDisks", std::vector<int>{0}, "hit map on MFT disks [0,1,2,3,4]. logical-OR of each double-sided disk"};
79+
Configurable<bool> requireMFTHitMap{"requireMFTHitMap", false, "flag to require MFT hit map"};
80+
Configurable<std::vector<int>> requiredMFTDisks{"requiredMFTDisks", std::vector<int>{0}, "hit map on MFT disks [0,1,2,3,4]. logical-OR of each double-sided disk"};
8081

8182
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2"};
8283
Configurable<float> cfgCentMin{"cfgCentMin", -1.f, "min. centrality"};
@@ -148,10 +149,12 @@ struct matchingMFT {
148149
hMuonType->GetXaxis()->SetBinLabel(4, "MCH-MID");
149150
hMuonType->GetXaxis()->SetBinLabel(5, "MCH standalone");
150151

152+
const AxisSpec axis_pt{{0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2, 3, 4, 5, 6, 7, 8, 9, 10}, "p_{T}^{gl} (GeV/c)"};
153+
151154
fRegistry.add("MFTMCHMID/primary/correct/hPt", "pT;p_{T} (GeV/c)", kTH1F, {{100, 0.0f, 10}}, false);
152-
fRegistry.add("MFTMCHMID/primary/correct/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {100, -6.f, -1.f}}, false);
153-
fRegistry.add("MFTMCHMID/primary/correct/hEtaPhi_MatchedMCHMID", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {100, -6.f, -1.f}}, false);
154-
fRegistry.add("MFTMCHMID/primary/correct/hDeltaEtaDeltaPhi", "#Delta#eta vs. #Delta#varphi;#Delta#varphi (rad.);#Delta#eta", kTH2F, {{180, -M_PI / 2, M_PI / 2}, {200, -1.f, +1.f}}, false);
155+
fRegistry.add("MFTMCHMID/primary/correct/hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {80, -5.f, -1.f}}, false);
156+
fRegistry.add("MFTMCHMID/primary/correct/hEtaPhi_MatchedMCHMID", "#eta vs. #varphi;#varphi (rad.);#eta", kTH2F, {{180, 0, 2 * M_PI}, {80, -5.f, -1.f}}, false);
157+
fRegistry.add("MFTMCHMID/primary/correct/hsDelta", "diff. between GL and associated SA;p_{T}^{gl} (GeV/c);#Delta#eta;#Delta#varphi (rad.);", kTHnSparseF, {axis_pt, {100, -0.5, +0.5}, {90, -M_PI / 4, M_PI / 4}}, false);
155158
fRegistry.add("MFTMCHMID/primary/correct/hDiffCollId", "difference in collision index;collisionId_{TTCA} - collisionId_{MP}", kTH1F, {{41, -20.5, +20.5}}, false);
156159
fRegistry.add("MFTMCHMID/primary/correct/hSign", "sign;sign", kTH1F, {{3, -1.5, +1.5}}, false);
157160
fRegistry.add("MFTMCHMID/primary/correct/hNclusters", "Nclusters;Nclusters", kTH1F, {{21, -0.5f, 20.5}}, false);
@@ -171,9 +174,6 @@ struct matchingMFT {
171174
fRegistry.add("MFTMCHMID/primary/correct/hDCAxResolutionvsPt", "DCA_{x} resolution vs. p_{T};p_{T} (GeV/c);DCA_{x} resolution (#mum);", kTH2F, {{100, 0, 10.f}, {500, 0, 500}}, false);
172175
fRegistry.add("MFTMCHMID/primary/correct/hDCAyResolutionvsPt", "DCA_{y} resolution vs. p_{T};p_{T} (GeV/c);DCA_{y} resolution (#mum);", kTH2F, {{100, 0, 10.f}, {500, 0, 500}}, false);
173176
fRegistry.add("MFTMCHMID/primary/correct/hDCAxyResolutionvsPt", "DCA_{xy} resolution vs. p_{T};p_{T} (GeV/c);DCA_{xy} resolution (#mum);", kTH2F, {{100, 0, 10.f}, {500, 0, 500}}, false);
174-
fRegistry.add("MFTMCHMID/primary/correct/hDCAxResolutionvsP", "DCA_{x} resolution vs. p;p (GeV/c);DCA_{x} resolution (#mum);", kTH2F, {{100, 0, 100.f}, {500, 0, 500}}, false);
175-
fRegistry.add("MFTMCHMID/primary/correct/hDCAyResolutionvsP", "DCA_{y} resolution vs. p;p (GeV/c);DCA_{y} resolution (#mum);", kTH2F, {{100, 0, 100.f}, {500, 0, 500}}, false);
176-
fRegistry.add("MFTMCHMID/primary/correct/hDCAxyResolutionvsP", "DCA_{xy} resolution vs. p;p (GeV/c);DCA_{xy} resolution (#mum);", kTH2F, {{100, 0, 100.f}, {500, 0, 500}}, false);
177177

178178
fRegistry.add("MFTMCHMID/primary/correct/hProdVtxZ", "prod. vtx Z of muon;V_{z} (cm)", kTH1F, {{200, -100, 100}}, false);
179179
fRegistry.add("MFTMCHMID/primary/correct/hRelDeltaPt", "pT resolution;p_{T}^{gen} (GeV/c);(p_{T}^{rec} - p_{T}^{gen})/p_{T}^{gen}", kTH2F, {{100, 0, 10}, {400, -1, +1}}, false);
@@ -183,8 +183,7 @@ struct matchingMFT {
183183
fRegistry.add("MFTMCHMID/primary/correct/hDeltaPhi_Neg", "#varphi resolution;p_{T}^{gen} (GeV/c);#varphi^{rec} - #varphi^{gen} (rad.)", kTH2F, {{100, 0, 10}, {400, -0.2, +0.2}}, false);
184184
fRegistry.addClone("MFTMCHMID/primary/correct/", "MFTMCHMID/primary/wrong/");
185185
fRegistry.addClone("MFTMCHMID/primary/", "MFTMCHMID/secondary/");
186-
187-
fRegistry.add("Generated/primary/hs", "gen. info;p_{T} (GeV/c);#eta;#varphi (rad.)", kTHnSparseF, {{100, 0.0f, 10}, {60, -5, -2}, {90, 0, 2.f * M_PI}}, false);
186+
// fRegistry.add("Generated/primary/hs", "gen. info;p_{T} (GeV/c);#eta;#varphi (rad.)", kTHnSparseF, {{100, 0.0f, 10}, {60, -5, -2}, {90, 0, 2.f * M_PI}}, false);
188187
}
189188

190189
bool isSelected(const float pt, const float eta, const float rAtAbsorberEnd, const float pDCA, const float chi2, const uint8_t trackType, const float dcaXY)
@@ -296,7 +295,6 @@ struct matchingMFT {
296295
o2::dataformats::GlobalFwdTrack propmuonAtPV = propagateMuon(fwdtrack, collision, propagationPoint::kToVertex);
297296
o2::dataformats::GlobalFwdTrack propmuonAtDCA = propagateMuon(fwdtrack, collision, propagationPoint::kToDCA);
298297

299-
float p = propmuonAtPV.getP();
300298
float pt = propmuonAtPV.getPt();
301299
float eta = propmuonAtPV.getEta();
302300
float phi = propmuonAtPV.getPhi();
@@ -355,10 +353,12 @@ struct matchingMFT {
355353
return;
356354
}
357355

358-
std::vector<bool> hasMFTs{hasMFT<0, 1>(mfttrack), hasMFT<2, 3>(mfttrack), hasMFT<4, 5>(mfttrack), hasMFT<6, 7>(mfttrack), hasMFT<8, 9>(mfttrack)};
359-
for (int i = 0; i < static_cast<int>(requireMFTDisks->size()); i++) {
360-
if (!hasMFTs[requireMFTDisks->at(i)]) {
361-
return;
356+
if (requireMFTHitMap) {
357+
std::vector<bool> hasMFTs{hasMFT<0, 1>(mfttrack), hasMFT<2, 3>(mfttrack), hasMFT<4, 5>(mfttrack), hasMFT<6, 7>(mfttrack), hasMFT<8, 9>(mfttrack)};
358+
for (int i = 0; i < static_cast<int>(requiredMFTDisks->size()); i++) {
359+
if (!hasMFTs[requiredMFTDisks->at(i)]) {
360+
return;
361+
}
362362
}
363363
}
364364

@@ -368,7 +368,7 @@ struct matchingMFT {
368368
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hPt"), pt);
369369
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hEtaPhi"), phi, eta);
370370
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hEtaPhi_MatchedMCHMID"), phiMatchedMCHMID, etaMatchedMCHMID);
371-
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hDeltaEtaDeltaPhi"), dphi, deta);
371+
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hsDelta"), pt, deta, dphi);
372372
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hDiffCollId"), collision.globalIndex() - fwdtrack.collisionId());
373373
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hSign"), fwdtrack.sign());
374374
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hNclusters"), fwdtrack.nClusters());
@@ -388,9 +388,6 @@ struct matchingMFT {
388388
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hDCAxResolutionvsPt"), pt, std::sqrt(cXXatDCA) * 1e+4); // convert cm to um
389389
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hDCAyResolutionvsPt"), pt, std::sqrt(cYYatDCA) * 1e+4); // convert cm to um
390390
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hDCAxyResolutionvsPt"), pt, sigma_dcaXY * 1e+4); // convert cm to um
391-
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hDCAxResolutionvsP"), p, std::sqrt(cXXatDCA) * 1e+4); // convert cm to um
392-
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hDCAyResolutionvsP"), p, std::sqrt(cYYatDCA) * 1e+4); // convert cm to um
393-
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hDCAxyResolutionvsP"), p, sigma_dcaXY * 1e+4); // convert cm to um
394391
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hProdVtxZ"), mcParticle_MFTMCHMID.vz());
395392

396393
fRegistry.fill(HIST("MFTMCHMID/primary/correct/hRelDeltaPt"), mcParticle_MFTMCHMID.pt(), (pt - mcParticle_MFTMCHMID.pt()) / mcParticle_MFTMCHMID.pt());
@@ -405,7 +402,7 @@ struct matchingMFT {
405402
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hPt"), pt);
406403
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hEtaPhi"), phi, eta);
407404
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hEtaPhi_MatchedMCHMID"), phiMatchedMCHMID, etaMatchedMCHMID);
408-
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hDeltaEtaDeltaPhi"), dphi, deta);
405+
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hsDelta"), pt, deta, dphi);
409406
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hDiffCollId"), collision.globalIndex() - fwdtrack.collisionId());
410407
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hSign"), fwdtrack.sign());
411408
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hNclusters"), fwdtrack.nClusters());
@@ -425,9 +422,6 @@ struct matchingMFT {
425422
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hDCAxResolutionvsPt"), pt, std::sqrt(cXXatDCA) * 1e+4); // convert cm to um
426423
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hDCAyResolutionvsPt"), pt, std::sqrt(cYYatDCA) * 1e+4); // convert cm to um
427424
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hDCAxyResolutionvsPt"), pt, sigma_dcaXY * 1e+4); // convert cm to um
428-
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hDCAxResolutionvsP"), p, std::sqrt(cXXatDCA) * 1e+4); // convert cm to um
429-
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hDCAyResolutionvsP"), p, std::sqrt(cYYatDCA) * 1e+4); // convert cm to um
430-
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hDCAxyResolutionvsP"), p, sigma_dcaXY * 1e+4); // convert cm to um
431425
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hProdVtxZ"), mcParticle_MFTMCHMID.vz());
432426
fRegistry.fill(HIST("MFTMCHMID/primary/wrong/hRelDeltaPt"), mcParticle_MFTMCHMID.pt(), (pt - mcParticle_MFTMCHMID.pt()) / mcParticle_MFTMCHMID.pt());
433427
if (mcParticle_MFTMCHMID.pdgCode() > 0) {
@@ -443,7 +437,7 @@ struct matchingMFT {
443437
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hPt"), pt);
444438
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hEtaPhi"), phi, eta);
445439
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hEtaPhi_MatchedMCHMID"), phiMatchedMCHMID, etaMatchedMCHMID);
446-
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hDeltaEtaDeltaPhi"), dphi, deta);
440+
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hsDelta"), pt, deta, dphi);
447441
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hDiffCollId"), collision.globalIndex() - fwdtrack.collisionId());
448442
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hSign"), fwdtrack.sign());
449443
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hNclusters"), fwdtrack.nClusters());
@@ -463,9 +457,6 @@ struct matchingMFT {
463457
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hDCAxResolutionvsPt"), pt, std::sqrt(cXXatDCA) * 1e+4); // convert cm to um
464458
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hDCAyResolutionvsPt"), pt, std::sqrt(cYYatDCA) * 1e+4); // convert cm to um
465459
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hDCAxyResolutionvsPt"), pt, sigma_dcaXY * 1e+4); // convert cm to um
466-
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hDCAxResolutionvsP"), p, std::sqrt(cXXatDCA) * 1e+4); // convert cm to um
467-
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hDCAyResolutionvsP"), p, std::sqrt(cYYatDCA) * 1e+4); // convert cm to um
468-
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hDCAxyResolutionvsP"), p, sigma_dcaXY * 1e+4); // convert cm to um
469460
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hProdVtxZ"), mcParticle_MFTMCHMID.vz());
470461
fRegistry.fill(HIST("MFTMCHMID/secondary/correct/hRelDeltaPt"), mcParticle_MFTMCHMID.pt(), (pt - mcParticle_MFTMCHMID.pt()) / mcParticle_MFTMCHMID.pt());
471462
if (mcParticle_MFTMCHMID.pdgCode() > 0) {
@@ -479,7 +470,7 @@ struct matchingMFT {
479470
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hPt"), pt);
480471
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hEtaPhi"), phi, eta);
481472
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hEtaPhi_MatchedMCHMID"), phiMatchedMCHMID, etaMatchedMCHMID);
482-
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hDeltaEtaDeltaPhi"), dphi, deta);
473+
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hsDelta"), pt, deta, dphi);
483474
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hDiffCollId"), collision.globalIndex() - fwdtrack.collisionId());
484475
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hSign"), fwdtrack.sign());
485476
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hNclusters"), fwdtrack.nClusters());
@@ -499,9 +490,6 @@ struct matchingMFT {
499490
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hDCAxResolutionvsPt"), pt, std::sqrt(cXXatDCA) * 1e+4); // convert cm to um
500491
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hDCAyResolutionvsPt"), pt, std::sqrt(cYYatDCA) * 1e+4); // convert cm to um
501492
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hDCAxyResolutionvsPt"), pt, sigma_dcaXY * 1e+4); // convert cm to um
502-
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hDCAxResolutionvsP"), p, std::sqrt(cXXatDCA) * 1e+4); // convert cm to um
503-
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hDCAyResolutionvsP"), p, std::sqrt(cYYatDCA) * 1e+4); // convert cm to um
504-
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hDCAxyResolutionvsP"), p, sigma_dcaXY * 1e+4); // convert cm to um
505493
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hProdVtxZ"), mcParticle_MFTMCHMID.vz());
506494
fRegistry.fill(HIST("MFTMCHMID/secondary/wrong/hRelDeltaPt"), mcParticle_MFTMCHMID.pt(), (pt - mcParticle_MFTMCHMID.pt()) / mcParticle_MFTMCHMID.pt());
507495
if (mcParticle_MFTMCHMID.pdgCode() > 0) {
@@ -529,24 +517,24 @@ struct matchingMFT {
529517
fRegistry.fill(HIST("Event/hMultFT0CvsMultNTracksPV"), collision.multFT0C(), collision.multNTracksPV());
530518
}
531519

532-
template <typename TMCParticles>
533-
void runGen(TMCParticles const& mcParticles)
534-
{
535-
for (const auto& mcParticle : mcParticles) {
536-
if (std::abs(mcParticle.pdgCode()) != 13) { // select true muon
537-
continue;
538-
}
539-
if (!(mcParticle.isPhysicalPrimary() || mcParticle.producedByGenerator())) {
540-
continue;
541-
}
542-
if (mcParticle.eta() < minEtaGL || maxEtaGL < mcParticle.eta()) {
543-
continue;
544-
}
545-
546-
fRegistry.fill(HIST("Generated/primary/hs"), mcParticle.pt(), mcParticle.eta(), mcParticle.phi());
547-
548-
} // end of mc particles
549-
}
520+
// template <typename TMCParticles>
521+
// void runGen(TMCParticles const& mcParticles)
522+
// {
523+
// for (const auto& mcParticle : mcParticles) {
524+
// if (std::abs(mcParticle.pdgCode()) != 13) { // select true muon
525+
// continue;
526+
// }
527+
// if (!(mcParticle.isPhysicalPrimary() || mcParticle.producedByGenerator())) {
528+
// continue;
529+
// }
530+
// if (mcParticle.eta() < minEtaGL || maxEtaGL < mcParticle.eta()) {
531+
// continue;
532+
// }
533+
534+
// fRegistry.fill(HIST("Generated/primary/hs"), mcParticle.pt(), mcParticle.eta(), mcParticle.phi());
535+
536+
// } // end of mc particles
537+
// }
550538

551539
SliceCache cache;
552540
PresliceUnsorted<aod::FwdTracks> perMFTTrack = o2::aod::fwdtrack::matchMFTTrackId;
@@ -558,7 +546,7 @@ struct matchingMFT {
558546
Filter collisionFilter_centrality = (cfgCentMin < o2::aod::cent::centFT0M && o2::aod::cent::centFT0M < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0A && o2::aod::cent::centFT0A < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0C && o2::aod::cent::centFT0C < cfgCentMax);
559547
using FilteredMyCollisions = soa::Filtered<MyCollisions>;
560548

561-
void processWithoutFTTCA(FilteredMyCollisions const& collisions, MyFwdTracks const& fwdtracks, MyMFTTracks const& mfttracks, aod::BCsWithTimestamps const&, aod::McParticles const& mcParticles)
549+
void processWithoutFTTCA(FilteredMyCollisions const& collisions, MyFwdTracks const& fwdtracks, MyMFTTracks const& mfttracks, aod::BCsWithTimestamps const&, aod::McParticles const&)
562550
{
563551
for (const auto& collision : collisions) {
564552
const auto& bc = collision.template bc_as<aod::BCsWithTimestamps>();
@@ -582,11 +570,11 @@ struct matchingMFT {
582570
fillHistograms(collision, fwdtrack, fwdtracks, mfttracks);
583571
} // end of fwdtrack loop
584572
} // end of collision loop
585-
runGen(mcParticles);
573+
// runGen(mcParticles);
586574
}
587575
PROCESS_SWITCH(matchingMFT, processWithoutFTTCA, "process without FTTCA", false);
588576

589-
void processWithFTTCA(FilteredMyCollisions const& collisions, MyFwdTracks const& fwdtracks, MyMFTTracks const& mfttracks, aod::BCsWithTimestamps const&, aod::FwdTrackAssoc const& fwdtrackIndices, aod::McParticles const& mcParticles)
577+
void processWithFTTCA(FilteredMyCollisions const& collisions, MyFwdTracks const& fwdtracks, MyMFTTracks const& mfttracks, aod::BCsWithTimestamps const&, aod::FwdTrackAssoc const& fwdtrackIndices, aod::McParticles const&)
590578
{
591579
for (const auto& collision : collisions) {
592580
const auto& bc = collision.template bc_as<aod::BCsWithTimestamps>();
@@ -611,7 +599,7 @@ struct matchingMFT {
611599
fillHistograms(collision, fwdtrack, fwdtracks, mfttracks);
612600
} // end of fwdtrack loop
613601
} // end of collision loop
614-
runGen(mcParticles);
602+
// runGen(mcParticles);
615603
}
616604
PROCESS_SWITCH(matchingMFT, processWithFTTCA, "process with FTTCA", true);
617605
};

0 commit comments

Comments
 (0)