Skip to content

Commit c3c09d5

Browse files
giovannimalfattoregiovannimalfattore
authored
[PWGLF] NucleiTask - Fix bugged hosts (#13474)
1 parent bf13110 commit c3c09d5

File tree

1 file changed

+106
-25
lines changed

1 file changed

+106
-25
lines changed

‎PWGLF/Tasks/Nuspex/LFNucleiBATask.cxx‎

Lines changed: 106 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -330,17 +330,29 @@ struct LFNucleiBATask {
330330
LOG(fatal) << "Can't enable processData and processMCReco in the same time, pick one!";
331331
}
332332
if (doprocessEvSgLossMC) {
333-
evLossHistos.add<TH1>("evLoss/hEvent", "Event loss histograms; ; counts", HistType::kTH1F, {{3, 0., 3.}});
333+
evLossHistos.add<TH1>("evLoss/hEvent", "Event loss histograms; ; counts", HistType::kTH1F, {{4, 0., 4.}});
334334
evLossHistos.get<TH1>(HIST("evLoss/hEvent"))->GetXaxis()->SetBinLabel(1, "All Gen.");
335335
evLossHistos.get<TH1>(HIST("evLoss/hEvent"))->GetXaxis()->SetBinLabel(2, "TVX (reco.)");
336-
evLossHistos.get<TH1>(HIST("evLoss/hEvent"))->GetXaxis()->SetBinLabel(3, "Sel8 (reco.)");
336+
evLossHistos.get<TH1>(HIST("evLoss/hEvent"))->GetXaxis()->SetBinLabel(3, "MC Sel8 (TVX + NoTFB) (reco.)");
337+
evLossHistos.get<TH1>(HIST("evLoss/hEvent"))->GetXaxis()->SetBinLabel(4, "Sel8 (reco.)");
337338

338339
evLossHistos.add<TH1>("evLoss/pt/hDeuteronTriggeredTVX", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
339340
evLossHistos.add<TH1>("evLoss/pt/hDeuteronTriggeredSel8", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
341+
evLossHistos.add<TH1>("evLoss/pt/hDeuteronTriggeredMCSel8", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
340342
evLossHistos.add<TH1>("evLoss/pt/hDeuteronGen", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
341343
evLossHistos.add<TH1>("evLoss/pt/hAntiDeuteronTriggeredTVX", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
344+
evLossHistos.add<TH1>("evLoss/pt/hAntiDeuteronTriggeredMCSel8", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
342345
evLossHistos.add<TH1>("evLoss/pt/hAntiDeuteronTriggeredSel8", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
343346
evLossHistos.add<TH1>("evLoss/pt/hAntiDeuteronGen", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
347+
348+
evLossHistos.add<TH1>("evLoss/pt/hHeliumTriggeredTVX", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
349+
evLossHistos.add<TH1>("evLoss/pt/hHeliumTriggeredSel8", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
350+
evLossHistos.add<TH1>("evLoss/pt/hHeliumTriggeredMCSel8", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
351+
evLossHistos.add<TH1>("evLoss/pt/hHeliumGen", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
352+
evLossHistos.add<TH1>("evLoss/pt/hAntiHeliumTriggeredTVX", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
353+
evLossHistos.add<TH1>("evLoss/pt/hAntiHeliumTriggeredSel8", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
354+
evLossHistos.add<TH1>("evLoss/pt/hAntiHeliumTriggeredMCSel8", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
355+
evLossHistos.add<TH1>("evLoss/pt/hAntiHeliumGen", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{100, 0., 5.}});
344356
}
345357
if (doprocessMCRecoLfPidEv) {
346358
spectraGen.add<TH1>("LfEv/pT_nocut", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
@@ -356,15 +368,27 @@ struct LFNucleiBATask {
356368
spectraGen.add<TH1>("LfEv/helium/pT_ITSROFBorder_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
357369
spectraGen.add<TH1>("LfEv/helium/pT_sel8_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
358370
spectraGen.add<TH1>("LfEv/helium/pT_MCsel8_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
359-
spectraGen.add<TH1>("LfEv/helium/pT_MCsel8_HePrim", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
371+
372+
spectraGen.add<TH1>("LfEv/helium/prim/pT_nocut_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
373+
spectraGen.add<TH1>("LfEv/helium/prim/pT_TVXtrigger_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
374+
spectraGen.add<TH1>("LfEv/helium/prim/pT_TFrameBorder_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
375+
spectraGen.add<TH1>("LfEv/helium/prim/pT_ITSROFBorder_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
376+
spectraGen.add<TH1>("LfEv/helium/prim/pT_sel8_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
377+
spectraGen.add<TH1>("LfEv/helium/prim/pT_MCsel8_He", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
360378

361379
spectraGen.add<TH1>("LfEv/helium/pT_nocut_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
362380
spectraGen.add<TH1>("LfEv/helium/pT_TVXtrigger_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
363381
spectraGen.add<TH1>("LfEv/helium/pT_TFrameBorder_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
364382
spectraGen.add<TH1>("LfEv/helium/pT_ITSROFBorder_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
365383
spectraGen.add<TH1>("LfEv/helium/pT_sel8_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
366384
spectraGen.add<TH1>("LfEv/helium/pT_MCsel8_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
367-
spectraGen.add<TH1>("LfEv/helium/pT_MCsel8_antiHePrim", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
385+
386+
spectraGen.add<TH1>("LfEv/helium/prim/pT_nocut_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
387+
spectraGen.add<TH1>("LfEv/helium/prim/pT_TVXtrigger_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
388+
spectraGen.add<TH1>("LfEv/helium/prim/pT_TFrameBorder_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
389+
spectraGen.add<TH1>("LfEv/helium/prim/pT_ITSROFBorder_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
390+
spectraGen.add<TH1>("LfEv/helium/prim/pT_sel8_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
391+
spectraGen.add<TH1>("LfEv/helium/prim/pT_MCsel8_antiHe", "Track #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); counts", HistType::kTH1F, {{ptHeAxis}});
368392
}
369393

370394
if (enableDebug) {
@@ -419,7 +443,7 @@ struct LFNucleiBATask {
419443
}
420444

421445
if (enableCentrality)
422-
histos.add<TH2>("event/h1VtxZ", "V_{z};V_{z} (in cm); counts", HistType::kTH1F, {{1500, -15, 15}, {binsPercentile, "Centrality FT0M"}});
446+
histos.add<TH2>("event/h1VtxZ", "V_{z};V_{z} (in cm); counts", HistType::kTH2F, {{1500, -15, 15}, {binsPercentile, "Centrality FT0M"}});
423447
else
424448
histos.add<TH1>("event/h1VtxZ", "V_{z};V_{z} (in cm); counts", HistType::kTH1F, {{1500, -15, 15}});
425449

@@ -2162,14 +2186,14 @@ struct LFNucleiBATask {
21622186
spectraGen.add("helium/histGenPtHe", "generated particles", HistType::kTH1F, {ptHeAxis});
21632187
spectraGen.add("helium/histGenPtHePrim", "generated particles", HistType::kTH1F, {ptHeAxis});
21642188
if (enableCentrality)
2165-
spectraGen.add("helium/histGenPtHePrim", "generated particles", HistType::kTH2F, {{ptHeAxis}, {binsPercentile}});
2189+
spectraGen.add("helium/histGenPtHePrimVsMult", "generated particles", HistType::kTH2F, {{ptHeAxis}, {binsPercentile}});
21662190
spectraGen.add("helium/histGenPtHeSec", "generated particles", HistType::kTH1F, {ptHeAxis});
21672191
spectraGen.add("helium/histSecTransportPtHe", "generated particles", HistType::kTH1F, {ptHeAxis});
21682192

21692193
spectraGen.add("helium/histGenPtantiHe", "generated particles", HistType::kTH1F, {ptHeAxis});
21702194
spectraGen.add("helium/histGenPtantiHePrim", "generated particles", HistType::kTH1F, {ptHeAxis});
21712195
if (enableCentrality)
2172-
spectraGen.add("helium/histGenPtantiHePrim", "generated particles", HistType::kTH2F, {{ptHeAxis}, {binsPercentile}});
2196+
spectraGen.add("helium/histGenPtantiHePrimVsMult", "generated particles", HistType::kTH2F, {{ptHeAxis}, {binsPercentile}});
21732197
spectraGen.add("helium/histGenPtantiHeSec", "generated particles", HistType::kTH1F, {ptHeAxis});
21742198
spectraGen.add("helium/histSecTransportPtantiHe", "generated particles", HistType::kTH1F, {ptHeAxis});
21752199
}
@@ -5752,7 +5776,7 @@ struct LFNucleiBATask {
57525776
if (isWeakDecay) {
57535777
histos.fill(HIST("tracks/helium/h1HeliumSpectraTrueSec_Z2"), 2 * hePt);
57545778
if (enableCentrality)
5755-
histos.fill(HIST("tracks/helium/h1HeliumSpectraTrueSecVsMult_Z2"), 2 * hePt);
5779+
histos.fill(HIST("tracks/helium/h2HeliumSpectraTrueSecVsMult_Z2"), 2 * hePt, event.centFT0M());
57565780
if (outFlagOptions.makeDCAAfterCutPlots) {
57575781
histos.fill(HIST("tracks/helium/dca/after/hDCAxyVsPtHeliumTrueSec"), hePt, track.dcaXY());
57585782
histos.fill(HIST("tracks/helium/dca/after/hDCAzVsPtHeliumTrueSec"), hePt, track.dcaZ());
@@ -6042,55 +6066,83 @@ struct LFNucleiBATask {
60426066

60436067
// No cut
60446068
spectraGen.fill(HIST("LfEv/pT_nocut"), pt);
6045-
if (pdg == PDGHelium)
6069+
if (pdg == PDGHelium) {
60466070
spectraGen.fill(HIST("LfEv/helium/pT_nocut_He"), pt);
6047-
if (pdg == -PDGHelium)
6071+
if (isPhysPrim)
6072+
spectraGen.fill(HIST("LfEv/helium/prim/pT_nocut_He"), pt);
6073+
}
6074+
if (pdg == -PDGHelium) {
60486075
spectraGen.fill(HIST("LfEv/helium/pT_nocut_antiHe"), pt);
6076+
if (isPhysPrim)
6077+
spectraGen.fill(HIST("LfEv/helium/prim/pT_nocut_antiHe"), pt);
6078+
}
60496079
// Trigger TVX
60506080
if (hasTVX) {
60516081
spectraGen.fill(HIST("LfEv/pT_TVXtrigger"), pt);
6052-
if (pdg == PDGHelium)
6082+
if (pdg == PDGHelium) {
60536083
spectraGen.fill(HIST("LfEv/helium/pT_TVXtrigger_He"), pt);
6054-
if (pdg == -PDGHelium)
6084+
if (isPhysPrim)
6085+
spectraGen.fill(HIST("LfEv/helium/prim/pT_TVXtrigger_He"), pt);
6086+
}
6087+
if (pdg == -PDGHelium) {
60556088
spectraGen.fill(HIST("LfEv/helium/pT_TVXtrigger_antiHe"), pt);
6089+
if (isPhysPrim)
6090+
spectraGen.fill(HIST("LfEv/helium/prim/pT_TVXtrigger_antiHe"), pt);
6091+
}
60566092
}
60576093
// No Time Frame Border
60586094
if (hasNoTFB) {
60596095
spectraGen.fill(HIST("LfEv/pT_TFrameBorder"), pt);
6060-
if (pdg == PDGHelium)
6096+
if (pdg == PDGHelium) {
60616097
spectraGen.fill(HIST("LfEv/helium/pT_TFrameBorder_He"), pt);
6062-
if (pdg == -PDGHelium)
6098+
if (isPhysPrim)
6099+
spectraGen.fill(HIST("LfEv/helium/prim/pT_TFrameBorder_He"), pt);
6100+
}
6101+
if (pdg == -PDGHelium) {
60636102
spectraGen.fill(HIST("LfEv/helium/pT_TFrameBorder_antiHe"), pt);
6103+
if (isPhysPrim)
6104+
spectraGen.fill(HIST("LfEv/helium/prim/pT_TFrameBorder_antiHe"), pt);
6105+
}
60646106
}
60656107
// No ITS ROF Frame Border
60666108
if (hasNoItsRofFB) {
60676109
spectraGen.fill(HIST("LfEv/pT_ITSROFBorder"), pt);
6068-
if (pdg == PDGHelium)
6110+
if (pdg == PDGHelium) {
60696111
spectraGen.fill(HIST("LfEv/helium/pT_ITSROFBorder_He"), pt);
6070-
if (pdg == -PDGHelium)
6112+
if (isPhysPrim)
6113+
spectraGen.fill(HIST("LfEv/helium/prim/pT_ITSROFBorder_He"), pt);
6114+
}
6115+
if (pdg == -PDGHelium) {
60716116
spectraGen.fill(HIST("LfEv/helium/pT_ITSROFBorder_antiHe"), pt);
6117+
if (isPhysPrim)
6118+
spectraGen.fill(HIST("LfEv/helium/prim/pT_ITSROFBorder_antiHe"), pt);
6119+
}
60726120
}
60736121
// Sel8 MC
60746122
if (hasTVX && hasNoTFB) {
60756123
spectraGen.fill(HIST("LfEv/pT_MCsel8"), pt);
60766124
if (pdg == PDGHelium) {
60776125
spectraGen.fill(HIST("LfEv/helium/pT_MCsel8_He"), pt);
60786126
if (isPhysPrim)
6079-
spectraGen.fill(HIST("LfEv/helium/pT_MCsel8_HePrim"), pt);
6127+
spectraGen.fill(HIST("LfEv/helium/prim/pT_MCsel8_He"), pt);
60806128
}
60816129
if (pdg == -PDGHelium) {
60826130
spectraGen.fill(HIST("LfEv/helium/pT_MCsel8_antiHe"), pt);
60836131
if (isPhysPrim)
6084-
spectraGen.fill(HIST("LfEv/helium/pT_MCsel8_antiHePrim"), pt);
6132+
spectraGen.fill(HIST("LfEv/helium/prim/pT_MCsel8_antiHe"), pt);
60856133
}
60866134
}
60876135
// Sel8 tag
60886136
if (hasSel8) {
60896137
spectraGen.fill(HIST("LfEv/pT_sel8"), pt);
60906138
if (pdg == PDGHelium)
60916139
spectraGen.fill(HIST("LfEv/helium/pT_sel8_He"), pt);
6140+
if (isPhysPrim)
6141+
spectraGen.fill(HIST("LfEv/helium/prim/pT_sel8_He"), pt);
60926142
if (pdg == -PDGHelium)
60936143
spectraGen.fill(HIST("LfEv/helium/pT_sel8_antiHe"), pt);
6144+
if (isPhysPrim)
6145+
spectraGen.fill(HIST("LfEv/helium/prim/pT_sel8_antiHe"), pt);
60946146
}
60956147
}
60966148
}
@@ -6304,10 +6356,9 @@ struct LFNucleiBATask {
63046356
if (pdgCode == PDGHelium) {
63056357
spectraGen.fill(HIST("helium/histGenPtHe"), ptMC);
63066358
if (isPhysPrim) {
6359+
spectraGen.fill(HIST("helium/histGenPtHePrim"), ptMC);
63076360
if (enableCentrality)
6308-
spectraGen.fill(HIST("helium/histGenPtHePrim"), ptMC, mcCollision.centFT0M());
6309-
else
6310-
spectraGen.fill(HIST("helium/histGenPtHePrim"), ptMC);
6361+
spectraGen.fill(HIST("helium/histGenPtHePrimVsMult"), ptMC, mcCollision.centFT0M());
63116362
}
63126363
if (!isPhysPrim && isProdByGen) {
63136364
{
@@ -6324,10 +6375,9 @@ struct LFNucleiBATask {
63246375
if (pdgCode == -PDGHelium) {
63256376
spectraGen.fill(HIST("helium/histGenPtantiHe"), ptMC);
63266377
if (isPhysPrim) {
6378+
spectraGen.fill(HIST("helium/histGenPtantiHePrim"), ptMC);
63276379
if (enableCentrality)
6328-
spectraGen.fill(HIST("helium/histGenPtantiHePrim"), ptMC, mcCollision.centFT0M());
6329-
else
6330-
spectraGen.fill(HIST("helium/histGenPtantiHePrim"), ptMC);
6380+
spectraGen.fill(HIST("helium/histGenPtantiHePrimVsMult"), ptMC, mcCollision.centFT0M());
63316381
}
63326382
if (!isPhysPrim && isProdByGen) {
63336383
{
@@ -6385,6 +6435,8 @@ struct LFNucleiBATask {
63856435

63866436
bool isSel8Event = false;
63876437
bool isTvxEvent = false;
6438+
bool isNoTFBEvent = false;
6439+
bool isMCSel8Event = false;
63886440

63896441
// Check if there is an event reconstructed for a generated event
63906442
for (const auto& recoColl : recoColls) {
@@ -6394,14 +6446,21 @@ struct LFNucleiBATask {
63946446
isTvxEvent = true;
63956447
if (recoColl.sel8())
63966448
isSel8Event = true;
6397-
if (isTvxEvent && isSel8Event)
6449+
if (recoColl.selection_bit(aod::evsel::kNoTimeFrameBorder))
6450+
isNoTFBEvent = true;
6451+
if (isTvxEvent && isSel8Event && isNoTFBEvent)
63986452
break; // Optimize loop
63996453
}
64006454

6455+
if (isTvxEvent && isNoTFBEvent)
6456+
isMCSel8Event = true;
6457+
64016458
if (isTvxEvent)
64026459
evLossHistos.fill(HIST("evLoss/hEvent"), 1.5);
64036460
if (isSel8Event)
64046461
evLossHistos.fill(HIST("evLoss/hEvent"), 2.5);
6462+
if (isMCSel8Event)
6463+
evLossHistos.fill(HIST("evLoss/hEvent"), 3.5);
64056464

64066465
// Loop over all the Generated level particles
64076466
for (const auto& mcPart : mcParticles) {
@@ -6419,13 +6478,35 @@ struct LFNucleiBATask {
64196478
evLossHistos.fill(HIST("evLoss/pt/hDeuteronTriggeredTVX"), pt);
64206479
if (isSel8Event)
64216480
evLossHistos.fill(HIST("evLoss/pt/hDeuteronTriggeredSel8"), pt);
6481+
if (isMCSel8Event)
6482+
evLossHistos.fill(HIST("evLoss/pt/hDeuteronTriggeredMCSel8"), pt);
64226483
}
64236484
if (pdg == -PDGDeuteron) {
64246485
evLossHistos.fill(HIST("evLoss/pt/hAntiDeuteronGen"), pt);
64256486
if (isTvxEvent)
64266487
evLossHistos.fill(HIST("evLoss/pt/hAntiDeuteronTriggeredTVX"), pt);
64276488
if (isSel8Event)
64286489
evLossHistos.fill(HIST("evLoss/pt/hAntiDeuteronTriggeredSel8"), pt);
6490+
if (isMCSel8Event)
6491+
evLossHistos.fill(HIST("evLoss/pt/hAntiDeuteronTriggeredMCSel8"), pt);
6492+
}
6493+
if (pdg == PDGHelium) {
6494+
evLossHistos.fill(HIST("evLoss/pt/hHeliumGen"), pt);
6495+
if (isTvxEvent)
6496+
evLossHistos.fill(HIST("evLoss/pt/hHeliumTriggeredTVX"), pt);
6497+
if (isSel8Event)
6498+
evLossHistos.fill(HIST("evLoss/pt/hHeliumTriggeredSel8"), pt);
6499+
if (isMCSel8Event)
6500+
evLossHistos.fill(HIST("evLoss/pt/hHeliumTriggeredMCSel8"), pt);
6501+
}
6502+
if (pdg == -PDGHelium) {
6503+
evLossHistos.fill(HIST("evLoss/pt/hAntiHeliumGen"), pt);
6504+
if (isTvxEvent)
6505+
evLossHistos.fill(HIST("evLoss/pt/hAntiHeliumTriggeredTVX"), pt);
6506+
if (isSel8Event)
6507+
evLossHistos.fill(HIST("evLoss/pt/hAntiHeliumTriggeredSel8"), pt);
6508+
if (isMCSel8Event)
6509+
evLossHistos.fill(HIST("evLoss/pt/hAntiHeliumTriggeredMCSel8"), pt);
64296510
}
64306511
}
64316512
}

0 commit comments

Comments
 (0)