Skip to content

Commit ac188d9

Browse files
KotliarovArKotliarovAr
authored
[PWGJE] Fix variable shadowing (#13771)
1 parent 1f4c127 commit ac188d9

File tree

1 file changed

+57
-46
lines changed

1 file changed

+57
-46
lines changed

‎PWGJE/Tasks/recoilJets.cxx‎

Lines changed: 57 additions & 46 deletions
Original file line numberDiff line numberDiff line change
@@ -100,31 +100,21 @@ struct RecoilJets {
100100
Configurable<float> meanFT0CPartLevel{"meanFT0CPartLevel", -1., "Mean number of charged part. within FT0C acceptance"};
101101

102102
// Parameters for recoil jet selection
103-
Configurable<float> ptTTrefMin{"ptTTrefMin", 5., "Minimum pT of reference TT"};
104-
Configurable<float> ptTTrefMax{"ptTTrefMax", 7., "Maximum pT of reference TT"};
105-
Configurable<float> ptTTsigMin{"ptTTsigMin", 20., "Minimum pT of signal TT"};
106-
Configurable<float> ptTTsigMax{"ptTTsigMax", 50., "Maximum pT of signal TT"};
103+
Configurable<std::vector<float>> ptTTref{"ptTTref", {5., 7}, "Transverse momentum (min, max) range for reference TT"};
104+
Configurable<std::vector<float>> ptTTsig{"ptTTsig", {20., 50}, "Transverse momentum (min, max) range for signal TT"};
107105
Configurable<float> recoilRegion{"recoilRegion", 0.6, "Width of recoil acceptance"};
108-
Configurable<float> minPhiForTTSelection{"minPhiForTTSelection", 0.0, "Min rectriction of phi angle for TT if phi is non-uniform"};
109-
Configurable<float> maxPhiForTTSelection{"maxPhiForTTSelection", 6.3, "Max rectriction of phi angle for TT if phi is non-uniform"};
106+
Configurable<std::vector<float>> phiRestrTTSelection{"phiRestrTTSelection", {0., 6.3}, "Restriction on phi angle (min, max) to search for TT"};
110107

111108
// Leading track and associated track
112-
Configurable<std::vector<float>> pTLeadTrack{"pTLeadTrack", {1., 3.}, "Transverse momenturm range for leading tracks"};
113-
Configurable<std::vector<float>> pTAssociatTrack{"pTAssociatTrack", {1., 3.}, "Transverse momenturm range for associated tracks"};
109+
Configurable<std::vector<float>> pTLeadTrack{"pTLeadTrack", {1., 3.}, "Transverse momenturm range (min, max) for leading tracks"};
110+
Configurable<std::vector<float>> pTAssociatTrack{"pTAssociatTrack", {1., 3.}, "Transverse momenturm range (min, max) for associated tracks"};
114111

115112
// List of configurable parameters for histograms
116-
Configurable<uint16_t> histJetPt{"histJetPt", 100, "Maximum value of jet pT shown in histograms"};
113+
Configurable<uint16_t> histJetPt{"histJetPt", 100, "Maximum value of jet pT stored in histograms"};
117114
Configurable<uint16_t> histMultBins{"histMultBins", 1000, "Number of bins for scaled FT0M multiplicity"};
118115
Configurable<uint16_t> histZDCTimeBins{"histZDCTimeBins", 240, "Number of bins for ZDC timing histograms"};
119116

120117
// Axes specification
121-
AxisSpec phiAngle{40, 0.0, constants::math::TwoPI, "#it{#varphi} (rad)"};
122-
AxisSpec deltaPhiAngle{52, 0.0, constants::math::PI, "#Delta#it{#varphi} (rad)"};
123-
AxisSpec pseudorap{40, -1., 1., "#it{#eta}"};
124-
AxisSpec pseudorapJets{20, -0.5, 0.5, "#it{#eta}_{jet}"};
125-
AxisSpec jetArea{50, 0.0, 5., "Area_{jet}"};
126-
AxisSpec rhoArea{60, 0.0, 60., "#it{#rho} #times Area_{jet}"};
127-
AxisSpec rho{50, 0.0, 50., "#it{#rho}"};
128118
ConfigurableAxis multFT0CThresh{"multFT0CThresh", {VARIABLE_WIDTH, 0.0, 0.133, 0.233, 0.367, 0.567, 0.767, 1.067, 1.4, 1.867, 2.5, 3.9, 5.4, 6.9, 20.}, "Percentiles of scaled FT0C: 100%, 90%, 80%, 70%, 60%, 50%, 40%, 30%, 20%, 10%, 1%, 0.1%, 0.01%"}; // default values for raw data
129119
ConfigurableAxis multFT0MThresh{"multFT0MThresh", {VARIABLE_WIDTH, 0.0, 0.167, 0.267, 0.4, 0.567, 0.8, 1.067, 1.4, 1.833, 2.433, 3.667, 5.1, 6.433, 20.}, "Percentiles of scaled FT0M: 100%, 90%, 80%, 70%, 60%, 50%, 40%, 30%, 20%, 10%, 1%, 0.1%, 0.01%"}; // default values for raw data
130120

@@ -165,6 +155,13 @@ struct RecoilJets {
165155
AxisSpec scaledFT0M{histMultBins, 0.0, 20., "FT0M^{*}"};
166156
AxisSpec zdcTiming{histZDCTimeBins, -30., 30., ""};
167157

158+
AxisSpec phiAngle{40, 0.0, constants::math::TwoPI, "#it{#varphi} (rad)"};
159+
AxisSpec deltaPhiAngle{52, 0.0, constants::math::PI, "#Delta#it{#varphi} (rad)"};
160+
AxisSpec pseudorap{40, -1., 1., "#it{#eta}"};
161+
AxisSpec pseudorapJets{20, -0.5, 0.5, "#it{#eta}_{jet}"};
162+
AxisSpec jetArea{50, 0.0, 5., "Area_{jet}"};
163+
AxisSpec rho{50, 0.0, 50., "#it{#rho}"};
164+
168165
std::string nameFT0Caxis = "FT0C / #LT FT0C #GT";
169166
std::string nameFT0Maxis = "FT0M^{*}";
170167

@@ -256,11 +253,11 @@ struct RecoilJets {
256253
spectra.add("hScaledFT0M_TTSig", "Events w. TT_{Sig}: scaled FT0M", kTH1F, {scaledFT0M});
257254

258255
// Rectricted phi range for TT selection
259-
spectra.add("hScaledFT0C_DPhi_JetPt_Corr_TTRef_RectrictedPhi", Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled FT0C & #Delta#varphi & #it{p}_{T, jet}^{ch}", static_cast<float>(minPhiForTTSelection), static_cast<float>(maxPhiForTTSelection)), kTH3F, {{multFT0CThresh, nameFT0Caxis}, deltaPhiAngle, jetPTcorr});
260-
spectra.add("hScaledFT0M_DPhi_JetPt_Corr_TTRef_RectrictedPhi", Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled FT0M & #Delta#varphi & #it{p}_{T, jet}^{ch}", static_cast<float>(minPhiForTTSelection), static_cast<float>(maxPhiForTTSelection)), kTH3F, {{multFT0MThresh, nameFT0Maxis}, deltaPhiAngle, jetPTcorr});
256+
spectra.add("hScaledFT0C_DPhi_JetPt_Corr_TTRef_RectrictedPhi", Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled FT0C & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiRestrTTSelection->at(0), phiRestrTTSelection->at(1)), kTH3F, {{multFT0CThresh, nameFT0Caxis}, deltaPhiAngle, jetPTcorr});
257+
spectra.add("hScaledFT0M_DPhi_JetPt_Corr_TTRef_RectrictedPhi", Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled FT0M & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiRestrTTSelection->at(0), phiRestrTTSelection->at(1)), kTH3F, {{multFT0MThresh, nameFT0Maxis}, deltaPhiAngle, jetPTcorr});
261258

262-
spectra.add("hScaledFT0C_DPhi_JetPt_Corr_TTSig_RectrictedPhi", Form("Events w. TT_{Sig} #in #varphi (%.2f, %.2f): scaled FT0C & #Delta#varphi & #it{p}_{T, jet}^{ch}", static_cast<float>(minPhiForTTSelection), static_cast<float>(maxPhiForTTSelection)), kTH3F, {{multFT0CThresh, nameFT0Caxis}, deltaPhiAngle, jetPTcorr});
263-
spectra.add("hScaledFT0M_DPhi_JetPt_Corr_TTSig_RectrictedPhi", Form("Events w. TT_{Sig} #in #varphi (%.2f, %.2f): scaled FT0M & #Delta#varphi & #it{p}_{T, jet}^{ch}", static_cast<float>(minPhiForTTSelection), static_cast<float>(maxPhiForTTSelection)), kTH3F, {{multFT0MThresh, nameFT0Maxis}, deltaPhiAngle, jetPTcorr});
259+
spectra.add("hScaledFT0C_DPhi_JetPt_Corr_TTSig_RectrictedPhi", Form("Events w. TT_{Sig} #in #varphi (%.2f, %.2f): scaled FT0C & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiRestrTTSelection->at(0), phiRestrTTSelection->at(1)), kTH3F, {{multFT0CThresh, nameFT0Caxis}, deltaPhiAngle, jetPTcorr});
260+
spectra.add("hScaledFT0M_DPhi_JetPt_Corr_TTSig_RectrictedPhi", Form("Events w. TT_{Sig} #in #varphi (%.2f, %.2f): scaled FT0M & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiRestrTTSelection->at(0), phiRestrTTSelection->at(1)), kTH3F, {{multFT0MThresh, nameFT0Maxis}, deltaPhiAngle, jetPTcorr});
264261
}
265262

266263
// List of MC particle level distributions
@@ -337,11 +334,11 @@ struct RecoilJets {
337334
spectra.add("hScaledFT0M_TTSig_Part", "Events w. TT_{Sig}: scaled FT0M", kTH1F, {scaledFT0M});
338335

339336
// Rectricted phi range for TT selection
340-
spectra.add("hScaledFT0C_DPhi_JetPt_Corr_TTRef_RectrictedPhi_Part", Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled FT0C & #Delta#varphi & #it{p}_{T, jet}^{ch}", static_cast<float>(minPhiForTTSelection), static_cast<float>(maxPhiForTTSelection)), kTH3F, {{multFT0CThreshPartLevel, nameFT0Caxis}, deltaPhiAngle, jetPTcorr});
341-
spectra.add("hScaledFT0M_DPhi_JetPt_Corr_TTRef_RectrictedPhi_Part", Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled FT0M & #Delta#varphi & #it{p}_{T, jet}^{ch}", static_cast<float>(minPhiForTTSelection), static_cast<float>(maxPhiForTTSelection)), kTH3F, {{multFT0MThreshPartLevel, nameFT0Maxis}, deltaPhiAngle, jetPTcorr});
337+
spectra.add("hScaledFT0C_DPhi_JetPt_Corr_TTRef_RectrictedPhi_Part", Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled FT0C & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiRestrTTSelection->at(0), phiRestrTTSelection->at(1)), kTH3F, {{multFT0CThreshPartLevel, nameFT0Caxis}, deltaPhiAngle, jetPTcorr});
338+
spectra.add("hScaledFT0M_DPhi_JetPt_Corr_TTRef_RectrictedPhi_Part", Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled FT0M & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiRestrTTSelection->at(0), phiRestrTTSelection->at(1)), kTH3F, {{multFT0MThreshPartLevel, nameFT0Maxis}, deltaPhiAngle, jetPTcorr});
342339

343-
spectra.add("hScaledFT0C_DPhi_JetPt_Corr_TTSig_RectrictedPhi_Part", Form("Events w. TT_{Sig} #in #varphi (%.2f, %.2f): scaled FT0C & #Delta#varphi & #it{p}_{T, jet}^{ch}", static_cast<float>(minPhiForTTSelection), static_cast<float>(maxPhiForTTSelection)), kTH3F, {{multFT0CThreshPartLevel, nameFT0Caxis}, deltaPhiAngle, jetPTcorr});
344-
spectra.add("hScaledFT0M_DPhi_JetPt_Corr_TTSig_RectrictedPhi_Part", Form("Events w. TT_{Sig} #in #varphi (%.2f, %.2f): scaled FT0M & #Delta#varphi & #it{p}_{T, jet}^{ch}", static_cast<float>(minPhiForTTSelection), static_cast<float>(maxPhiForTTSelection)), kTH3F, {{multFT0MThreshPartLevel, nameFT0Maxis}, deltaPhiAngle, jetPTcorr});
340+
spectra.add("hScaledFT0C_DPhi_JetPt_Corr_TTSig_RectrictedPhi_Part", Form("Events w. TT_{Sig} #in #varphi (%.2f, %.2f): scaled FT0C & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiRestrTTSelection->at(0), phiRestrTTSelection->at(1)), kTH3F, {{multFT0CThreshPartLevel, nameFT0Caxis}, deltaPhiAngle, jetPTcorr});
341+
spectra.add("hScaledFT0M_DPhi_JetPt_Corr_TTSig_RectrictedPhi_Part", Form("Events w. TT_{Sig} #in #varphi (%.2f, %.2f): scaled FT0M & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiRestrTTSelection->at(0), phiRestrTTSelection->at(1)), kTH3F, {{multFT0MThreshPartLevel, nameFT0Maxis}, deltaPhiAngle, jetPTcorr});
345342
}
346343

347344
// Jet matching: part. vs. det.
@@ -470,6 +467,11 @@ struct RecoilJets {
470467
spectra.get<TH2>(HIST("hIsFT0SignalComeFromCollPerRun"))->GetXaxis()->SetBinLabel(3, "Coll. w. BC");
471468
spectra.get<TH2>(HIST("hIsFT0SignalComeFromCollPerRun"))->GetXaxis()->SetBinLabel(4, "Coll. w/o BC");
472469

470+
// FT0 signal for the case when there is no associated BC
471+
spectra.add("hFT0AsignalWithoutBC", "", kTH2F, {{2000, 0.0, 40000., "FT0A"}, {nRunsOO, 0., nRunsOO * 1.}});
472+
spectra.add("hFT0CsignalWithoutBC", "", kTH2F, {{2000, 0.0, 40000., "FT0C"}, {nRunsOO, 0., nRunsOO * 1.}});
473+
spectra.add("hFT0MsignalWithoutBC", "", kTH2F, {{2000, 0.0, 40000., "FT0M"}, {nRunsOO, 0., nRunsOO * 1.}});
474+
473475
// Rename Y axis with Run numbers
474476
for (int iRun = 0; iRun < nRunsOO; ++iRun) {
475477
spectra.get<TH3>(HIST("hScaledFT0APerRunPerSetOfFlags"))->GetYaxis()->SetBinLabel(iRun + 1, runNumbersOO[iRun]);
@@ -481,6 +483,10 @@ struct RecoilJets {
481483
spectra.get<TH3>(HIST("hFT0MPerRunPerSetOfFlags"))->GetYaxis()->SetBinLabel(iRun + 1, runNumbersOO[iRun]);
482484

483485
spectra.get<TH2>(HIST("hIsFT0SignalComeFromCollPerRun"))->GetYaxis()->SetBinLabel(iRun + 1, runNumbersOO[iRun]);
486+
487+
spectra.get<TH2>(HIST("hFT0AsignalWithoutBC"))->GetYaxis()->SetBinLabel(iRun + 1, runNumbersOO[iRun]);
488+
spectra.get<TH2>(HIST("hFT0CsignalWithoutBC"))->GetYaxis()->SetBinLabel(iRun + 1, runNumbersOO[iRun]);
489+
spectra.get<TH2>(HIST("hFT0MsignalWithoutBC"))->GetYaxis()->SetBinLabel(iRun + 1, runNumbersOO[iRun]);
484490
}
485491

486492
// Rename Z axis with event selection flags
@@ -502,9 +508,9 @@ struct RecoilJets {
502508
}
503509

504510
// Fill histograms with raw or MC det. level data
505-
template <typename Collision, typename Jets, typename Tracks>
506-
void fillHistograms(Collision const& collision, Jets const& jets,
507-
Tracks const& tracks, float weight = 1.)
511+
template <typename JCollision, typename Jets, typename JTracks>
512+
void fillHistograms(JCollision const& collision, Jets const& jets,
513+
JTracks const& tracks, float weight = 1.)
508514
{
509515
bool bSigEv = false;
510516
std::vector<double> vPhiOfTT;
@@ -538,13 +544,13 @@ struct RecoilJets {
538544
spectra.fill(HIST("hScaledFT0MTrackPtEtaPhi"), scaledFT0M, trackPt, track.eta(), trackPhi, weight);
539545

540546
// Search for TT candidate
541-
if (bSigEv && (trackPt > ptTTsigMin && trackPt < ptTTsigMax)) {
547+
if (bSigEv && (trackPt > ptTTsig->at(0) && trackPt < ptTTsig->at(1))) {
542548
vPhiOfTT.push_back(trackPhi);
543549
spectra.fill(HIST("hTTSig_pT"), trackPt, weight);
544550
++nTT;
545551
}
546552

547-
if (!bSigEv && (trackPt > ptTTrefMin && trackPt < ptTTrefMax)) {
553+
if (!bSigEv && (trackPt > ptTTref->at(0) && trackPt < ptTTref->at(1))) {
548554
vPhiOfTT.push_back(trackPhi);
549555
++nTT;
550556
}
@@ -601,7 +607,7 @@ struct RecoilJets {
601607
spectra.fill(HIST("hScaledFT0M_DPhi_JetPt_TTSig"), scaledFT0M, dphi, jetPt, weight);
602608
spectra.fill(HIST("hJetArea_JetPt_Rho_TTSig"), jetArea, jetPt, rho, weight);
603609

604-
if (phiTT > minPhiForTTSelection && phiTT < maxPhiForTTSelection) {
610+
if (phiTT > phiRestrTTSelection->at(0) && phiTT < phiRestrTTSelection->at(1)) {
605611
spectra.fill(HIST("hScaledFT0C_DPhi_JetPt_Corr_TTSig_RectrictedPhi"), scaledFT0C, dphi, jetPtCorr, weight);
606612
spectra.fill(HIST("hScaledFT0M_DPhi_JetPt_Corr_TTSig_RectrictedPhi"), scaledFT0M, dphi, jetPtCorr, weight);
607613
}
@@ -621,7 +627,7 @@ struct RecoilJets {
621627
spectra.fill(HIST("hScaledFT0M_DPhi_JetPt_TTRef"), scaledFT0M, dphi, jetPt, weight);
622628
spectra.fill(HIST("hJetArea_JetPt_Rho_TTRef"), jetArea, jetPt, rho, weight);
623629

624-
if (phiTT > minPhiForTTSelection && phiTT < maxPhiForTTSelection) {
630+
if (phiTT > phiRestrTTSelection->at(0) && phiTT < phiRestrTTSelection->at(1)) {
625631
spectra.fill(HIST("hScaledFT0C_DPhi_JetPt_Corr_TTRef_RectrictedPhi"), scaledFT0C, dphi, jetPtCorr, weight);
626632
spectra.fill(HIST("hScaledFT0M_DPhi_JetPt_Corr_TTRef_RectrictedPhi"), scaledFT0M, dphi, jetPtCorr, weight);
627633
}
@@ -638,9 +644,9 @@ struct RecoilJets {
638644
}
639645
}
640646

641-
template <typename Collision, typename Jets, typename Particles>
642-
void fillMCPHistograms(Collision const& collision, Jets const& jets,
643-
Particles const& particles, float weight = 1.)
647+
template <typename JCollision, typename Jets, typename JParticles>
648+
void fillMCPHistograms(JCollision const& collision, Jets const& jets,
649+
JParticles const& particles, float weight = 1.)
644650
{
645651
bool bSigEv = false;
646652
std::vector<double> vPhiOfTT;
@@ -676,12 +682,12 @@ struct RecoilJets {
676682
spectra.fill(HIST("hScaledFT0CPartPtEtaPhi"), scaledFT0C, particlePt, particle.eta(), particlePhi, weight);
677683
spectra.fill(HIST("hScaledFT0MPartPtEtaPhi"), scaledFT0M, particlePt, particle.eta(), particlePhi, weight);
678684

679-
if (bSigEv && (particlePt > ptTTsigMin && particlePt < ptTTsigMax)) {
685+
if (bSigEv && (particlePt > ptTTsig->at(0) && particlePt < ptTTsig->at(1))) {
680686
vPhiOfTT.push_back(particlePhi);
681687
++nTT;
682688
}
683689

684-
if (!bSigEv && (particlePt > ptTTrefMin && particlePt < ptTTrefMax)) {
690+
if (!bSigEv && (particlePt > ptTTref->at(0) && particlePt < ptTTref->at(1))) {
685691
vPhiOfTT.push_back(particlePhi);
686692
++nTT;
687693
}
@@ -736,7 +742,7 @@ struct RecoilJets {
736742

737743
spectra.fill(HIST("hJetArea_JetPt_Rho_TTSig_Part"), jetArea, jetPt, rho, weight);
738744

739-
if (phiTT > minPhiForTTSelection && phiTT < maxPhiForTTSelection) {
745+
if (phiTT > phiRestrTTSelection->at(0) && phiTT < phiRestrTTSelection->at(1)) {
740746
spectra.fill(HIST("hScaledFT0C_DPhi_JetPt_Corr_TTSig_RectrictedPhi_Part"), scaledFT0C, dphi, jetPtCorr, weight);
741747
spectra.fill(HIST("hScaledFT0M_DPhi_JetPt_Corr_TTSig_RectrictedPhi_Part"), scaledFT0M, dphi, jetPtCorr, weight);
742748
}
@@ -759,7 +765,7 @@ struct RecoilJets {
759765

760766
spectra.fill(HIST("hJetArea_JetPt_Rho_TTRef_Part"), jetArea, jetPt, rho, weight);
761767

762-
if (phiTT > minPhiForTTSelection && phiTT < maxPhiForTTSelection) {
768+
if (phiTT > phiRestrTTSelection->at(0) && phiTT < phiRestrTTSelection->at(1)) {
763769
spectra.fill(HIST("hScaledFT0C_DPhi_JetPt_Corr_TTRef_RectrictedPhi_Part"), scaledFT0C, dphi, jetPtCorr, weight);
764770
spectra.fill(HIST("hScaledFT0M_DPhi_JetPt_Corr_TTRef_RectrictedPhi_Part"), scaledFT0M, dphi, jetPtCorr, weight);
765771
}
@@ -776,8 +782,8 @@ struct RecoilJets {
776782
}
777783
}
778784

779-
template <typename TracksTable, typename JetsBase, typename JetsTag>
780-
void fillMatchedHistograms(TracksTable const& tracks,
785+
template <typename JTracksTable, typename JetsBase, typename JetsTag>
786+
void fillMatchedHistograms(JTracksTable const& tracks,
781787
JetsBase const& jetsBase, JetsTag const& jetsTag,
782788
float weight = 1.)
783789
{
@@ -788,7 +794,7 @@ struct RecoilJets {
788794
if (skipTrack(track))
789795
continue;
790796

791-
if (track.pt() > ptTTsigMin && track.pt() < ptTTsigMax) {
797+
if (track.pt() > ptTTsig->at(0) && track.pt() < ptTTsig->at(1)) {
792798
vPhiOfTT.push_back(track.phi());
793799
}
794800
}
@@ -805,8 +811,8 @@ struct RecoilJets {
805811
}
806812
}
807813

808-
template <typename Collision>
809-
void fillMultiplicityHistogramsOO(Collision const& collision,
814+
template <typename JCollision>
815+
void fillMultiplicityHistogramsOO(JCollision const& collision,
810816
float weight = 1.)
811817
{
812818
float multFT0A = collision.multFT0A();
@@ -861,8 +867,8 @@ struct RecoilJets {
861867
spectra.fill(HIST("hScaleMultFT0M_vs_ZPA_vs_ZPC"), scaledFT0M, multZPA, multZPC, weight);
862868
}
863869

864-
template <typename CollisionMC>
865-
void fillMultiplicityHistogramsPartLevelMC(CollisionMC const& collision,
870+
template <typename JCollisionMC>
871+
void fillMultiplicityHistogramsPartLevelMC(JCollisionMC const& collision,
866872
float weight = 1.)
867873
{
868874
spectra.fill(HIST("hMultFT0APartLevel"), collision.multFT0A(), weight);
@@ -976,9 +982,14 @@ struct RecoilJets {
976982
spectra.fill(HIST("hScaledFT0M_TracksPV"), scaledFT0M, multNContribs, weight);
977983

978984
if (foundBC.foundFT0Id() > 0) // -1 if does not
985+
{
979986
spectra.fill(HIST("hIsFT0SignalComeFromCollPerRun"), 0.5, fillNumber, weight);
980-
else
987+
} else {
981988
spectra.fill(HIST("hIsFT0SignalComeFromCollPerRun"), 1.5, fillNumber, weight);
989+
spectra.fill(HIST("hFT0AsignalWithoutBC"), multFT0A, fillNumber, weight);
990+
spectra.fill(HIST("hFT0CsignalWithoutBC"), multFT0C, fillNumber, weight);
991+
spectra.fill(HIST("hFT0MsignalWithoutBC"), multFT0M, fillNumber, weight);
992+
}
982993

983994
if (collision.foundBCId() > 0)
984995
spectra.fill(HIST("hIsFT0SignalComeFromCollPerRun"), 2.5, fillNumber, weight);

0 commit comments

Comments
 (0)