@@ -230,33 +230,33 @@ struct OnTheFlyTofPid {
230230
231231 std::string particleNames[kParticles ] = {" #it{e}" , " #it{#mu}" , " #it{#pi}" , " #it{K}" , " #it{p}" };
232232 std::string particleNames2[kParticles ] = {" Elec" , " Muon" , " Pion" , " Kaon" , " Prot" };
233- for (int i_true = 0 ; i_true < kParticles ; i_true ++) {
233+ for (int iTrue = 0 ; iTrue < kParticles ; iTrue ++) {
234234 auto addHistogram = [&](const std::string& name, const AxisSpec& axis) {
235235 return histos.add <TH2>(name, " " , kTH2F , {axisMomentum, axis});
236236 };
237237
238- const AxisSpec axisTrackTimeRes{plotsConfig.nBinsTimeRes , 0 .0f , +200 .0f , " Track time resolution - " + particleNames[i_true ] + " (ps)" };
239- h2dInnerTimeResTrack[i_true ] = addHistogram (" iTOF/res/h2dInnerTimeResTrack" + particleNames2[i_true ] + " VsP" , axisTrackTimeRes);
240- h2dOuterTimeResTrack[i_true ] = addHistogram (" oTOF/res/h2dOuterTimeResTrack" + particleNames2[i_true ] + " VsP" , axisTrackTimeRes);
241- const AxisSpec axisTotalTimeRes{plotsConfig.nBinsTimeRes , 0 .0f , +200 .0f , " Total time resolution - " + particleNames[i_true ] + " (ps)" };
242- h2dInnerTimeResTotal[i_true ] = addHistogram (" iTOF/res/h2dInnerTimeResTotal" + particleNames2[i_true ] + " VsP" , axisTotalTimeRes);
243- h2dOuterTimeResTotal[i_true ] = addHistogram (" oTOF/res/h2dOuterTimeResTotal" + particleNames2[i_true ] + " VsP" , axisTotalTimeRes);
244- for (int i_hyp = 0 ; i_hyp < kParticles ; i_hyp ++) {
245- std::string nameTitleInner = " h2dInnerNsigmaTrue" + particleNames2[i_true ] + " Vs" + particleNames2[i_hyp ] + " Hypothesis" ;
246- std::string nameTitleOuter = " h2dOuterNsigmaTrue" + particleNames2[i_true ] + " Vs" + particleNames2[i_hyp ] + " Hypothesis" ;
247- std::string nameTitleInnerDelta = " h2dInnerDeltaTrue" + particleNames2[i_true ] + " Vs" + particleNames2[i_hyp ] + " Hypothesis" ;
248- std::string nameTitleOuterDelta = " h2dOuterDeltaTrue" + particleNames2[i_true ] + " Vs" + particleNames2[i_hyp ] + " Hypothesis" ;
238+ const AxisSpec axisTrackTimeRes{plotsConfig.nBinsTimeRes , 0 .0f , +200 .0f , " Track time resolution - " + particleNames[iTrue ] + " (ps)" };
239+ h2dInnerTimeResTrack[iTrue ] = addHistogram (" iTOF/res/h2dInnerTimeResTrack" + particleNames2[iTrue ] + " VsP" , axisTrackTimeRes);
240+ h2dOuterTimeResTrack[iTrue ] = addHistogram (" oTOF/res/h2dOuterTimeResTrack" + particleNames2[iTrue ] + " VsP" , axisTrackTimeRes);
241+ const AxisSpec axisTotalTimeRes{plotsConfig.nBinsTimeRes , 0 .0f , +200 .0f , " Total time resolution - " + particleNames[iTrue ] + " (ps)" };
242+ h2dInnerTimeResTotal[iTrue ] = addHistogram (" iTOF/res/h2dInnerTimeResTotal" + particleNames2[iTrue ] + " VsP" , axisTotalTimeRes);
243+ h2dOuterTimeResTotal[iTrue ] = addHistogram (" oTOF/res/h2dOuterTimeResTotal" + particleNames2[iTrue ] + " VsP" , axisTotalTimeRes);
244+ for (int iHyp = 0 ; iHyp < kParticles ; iHyp ++) {
245+ std::string nameTitleInner = " h2dInnerNsigmaTrue" + particleNames2[iTrue ] + " Vs" + particleNames2[iHyp ] + " Hypothesis" ;
246+ std::string nameTitleOuter = " h2dOuterNsigmaTrue" + particleNames2[iTrue ] + " Vs" + particleNames2[iHyp ] + " Hypothesis" ;
247+ std::string nameTitleInnerDelta = " h2dInnerDeltaTrue" + particleNames2[iTrue ] + " Vs" + particleNames2[iHyp ] + " Hypothesis" ;
248+ std::string nameTitleOuterDelta = " h2dOuterDeltaTrue" + particleNames2[iTrue ] + " Vs" + particleNames2[iHyp ] + " Hypothesis" ;
249249 const AxisSpec axisX{plotsConfig.doSeparationVsPt .value ? axisPt : axisMomentum};
250- const AxisSpec axisNsigmaCorrect{plotsConfig.nBinsNsigmaCorrectSpecies , plotsConfig.minNsigmaRange , plotsConfig.maxNsigmaRange , " N#sigma - True " + particleNames[i_true ] + " vs " + particleNames[i_hyp ] + " hypothesis" };
251- const AxisSpec axisDeltaCorrect{plotsConfig.nBinsDeltaCorrectSpecies , plotsConfig.minDeltaRange , plotsConfig.maxDeltaRange , " #Delta - True " + particleNames[i_true ] + " vs " + particleNames[i_hyp ] + " hypothesis" };
252- const AxisSpec axisNsigmaWrong{plotsConfig.nBinsNsigmaWrongSpecies , plotsConfig.minNsigmaRange , plotsConfig.maxNsigmaRange , " N#sigma - True " + particleNames[i_true ] + " vs " + particleNames[i_hyp ] + " hypothesis" };
253- const AxisSpec axisDeltaWrong{plotsConfig.nBinsDeltaWrongSpecies , plotsConfig.minDeltaRange , plotsConfig.maxDeltaRange , " #Delta - True " + particleNames[i_true ] + " vs " + particleNames[i_hyp ] + " hypothesis" };
254- const AxisSpec axisNSigma{i_true == i_hyp ? axisNsigmaCorrect : axisNsigmaWrong};
255- const AxisSpec axisDelta{i_true == i_hyp ? axisDeltaCorrect : axisDeltaWrong};
256- h2dInnerNsigmaTrue[i_true][i_hyp ] = histos.add <TH2>(" iTOF/nsigma/h2dInnerNsigmaTrue" + particleNames2[i_true ] + " Vs" + particleNames2[i_hyp ] + " Hypothesis" , " " , kTH2F , {axisX, axisNSigma});
257- h2dOuterNsigmaTrue[i_true][i_hyp ] = histos.add <TH2>(" oTOF/nsigma/h2dOuterNsigmaTrue" + particleNames2[i_true ] + " Vs" + particleNames2[i_hyp ] + " Hypothesis" , " " , kTH2F , {axisX, axisNSigma});
258- h2dInnerDeltaTrue[i_true][i_hyp ] = histos.add <TH2>(" iTOF/delta/h2dInnerDeltaTrue" + particleNames2[i_true ] + " Vs" + particleNames2[i_hyp ] + " Hypothesis" , " " , kTH2F , {axisX, axisDelta});
259- h2dOuterDeltaTrue[i_true][i_hyp ] = histos.add <TH2>(" oTOF/delta/h2dOuterDeltaTrue" + particleNames2[i_true ] + " Vs" + particleNames2[i_hyp ] + " Hypothesis" , " " , kTH2F , {axisX, axisDelta});
250+ const AxisSpec axisNsigmaCorrect{plotsConfig.nBinsNsigmaCorrectSpecies , plotsConfig.minNsigmaRange , plotsConfig.maxNsigmaRange , " N#sigma - True " + particleNames[iTrue ] + " vs " + particleNames[iHyp ] + " hypothesis" };
251+ const AxisSpec axisDeltaCorrect{plotsConfig.nBinsDeltaCorrectSpecies , plotsConfig.minDeltaRange , plotsConfig.maxDeltaRange , " #Delta - True " + particleNames[iTrue ] + " vs " + particleNames[iHyp ] + " hypothesis" };
252+ const AxisSpec axisNsigmaWrong{plotsConfig.nBinsNsigmaWrongSpecies , plotsConfig.minNsigmaRange , plotsConfig.maxNsigmaRange , " N#sigma - True " + particleNames[iTrue ] + " vs " + particleNames[iHyp ] + " hypothesis" };
253+ const AxisSpec axisDeltaWrong{plotsConfig.nBinsDeltaWrongSpecies , plotsConfig.minDeltaRange , plotsConfig.maxDeltaRange , " #Delta - True " + particleNames[iTrue ] + " vs " + particleNames[iHyp ] + " hypothesis" };
254+ const AxisSpec axisNSigma{iTrue == iHyp ? axisNsigmaCorrect : axisNsigmaWrong};
255+ const AxisSpec axisDelta{iTrue == iHyp ? axisDeltaCorrect : axisDeltaWrong};
256+ h2dInnerNsigmaTrue[iTrue][iHyp ] = histos.add <TH2>(" iTOF/nsigma/h2dInnerNsigmaTrue" + particleNames2[iTrue ] + " Vs" + particleNames2[iHyp ] + " Hypothesis" , " " , kTH2F , {axisX, axisNSigma});
257+ h2dOuterNsigmaTrue[iTrue][iHyp ] = histos.add <TH2>(" oTOF/nsigma/h2dOuterNsigmaTrue" + particleNames2[iTrue ] + " Vs" + particleNames2[iHyp ] + " Hypothesis" , " " , kTH2F , {axisX, axisNSigma});
258+ h2dInnerDeltaTrue[iTrue][iHyp ] = histos.add <TH2>(" iTOF/delta/h2dInnerDeltaTrue" + particleNames2[iTrue ] + " Vs" + particleNames2[iHyp ] + " Hypothesis" , " " , kTH2F , {axisX, axisDelta});
259+ h2dOuterDeltaTrue[iTrue][iHyp ] = histos.add <TH2>(" oTOF/delta/h2dOuterDeltaTrue" + particleNames2[iTrue ] + " Vs" + particleNames2[iHyp ] + " Hypothesis" , " " , kTH2F , {axisX, axisDelta});
260260 }
261261 }
262262 }
0 commit comments