@@ -127,8 +127,8 @@ struct OnTheFlyTofPid {
127127 // for handling basic QA histograms if requested
128128 HistogramRegistry histos{" Histos" , {}, OutputObjHandlingPolicy::AnalysisObject};
129129 static constexpr int kParticles = 5 ;
130- std::string particle_names1 [kParticles ] = {" #it{e}" , " #it{#mu}" , " #it{#pi}" , " #it{K}" , " #it{p}" };
131- std::string particle_names2 [kParticles ] = {" Elec" , " Muon" , " Pion" , " Kaon" , " Prot" };
130+ std::string particleNames [kParticles ] = {" #it{e}" , " #it{#mu}" , " #it{#pi}" , " #it{K}" , " #it{p}" };
131+ std::string particleNames2 [kParticles ] = {" Elec" , " Muon" , " Pion" , " Kaon" , " Prot" };
132132 static constexpr int kIdPion = 2 ;
133133
134134 void init (o2::framework::InitContext&)
@@ -197,38 +197,38 @@ struct OnTheFlyTofPid {
197197 histos.add (" h2dRelativeEtaResolution" , " h2dRelativeEtaResolution" , kTH2F , {axisEta, axisRelativeEta});
198198
199199 for (int i_true = 0 ; i_true < kParticles ; i_true++) {
200- std::string name_title_inner_track_res = " h2dInnerTimeResTrack" + particle_names2 [i_true] + " VsP" ;
201- std::string name_title_inner_total_res = " h2dInnerTimeResTotal" + particle_names2 [i_true] + " VsP" ;
202- std::string name_title_outer_track_res = " h2dOuterTimeResTrack" + particle_names2 [i_true] + " VsP" ;
203- std::string name_title_outer_total_res = " h2dOuterTimeResTotal" + particle_names2 [i_true] + " VsP" ;
204- const AxisSpec axisTrackTimeRes{static_cast <int >(plotsConfig.nBinsTimeRes ), 0 .0f , +200 .0f , " Track time resolution - " + particle_names1 [i_true] + " (ps)" };
205- const AxisSpec axisTotalTimeRes{static_cast <int >(plotsConfig.nBinsTimeRes ), 0 .0f , +200 .0f , " Total time resolution - " + particle_names1 [i_true] + " (ps)" };
206- h2dInnerTimeResTrack[i_true] = histos.add <TH2>(name_title_inner_track_res .c_str (), name_title_inner_track_res .c_str (), kTH2F , {axisMomentum, axisTrackTimeRes});
207- h2dInnerTimeResTotal[i_true] = histos.add <TH2>(name_title_inner_total_res .c_str (), name_title_inner_total_res .c_str (), kTH2F , {axisMomentum, axisTotalTimeRes});
208- h2dOuterTimeResTrack[i_true] = histos.add <TH2>(name_title_outer_track_res .c_str (), name_title_outer_track_res .c_str (), kTH2F , {axisMomentum, axisTrackTimeRes});
209- h2dOuterTimeResTotal[i_true] = histos.add <TH2>(name_title_outer_total_res .c_str (), name_title_outer_total_res .c_str (), kTH2F , {axisMomentum, axisTotalTimeRes});
200+ std::string nameTitleInnerTrackRes = " h2dInnerTimeResTrack" + particleNames2 [i_true] + " VsP" ;
201+ std::string nameTitleTotalRes = " h2dInnerTimeResTotal" + particleNames2 [i_true] + " VsP" ;
202+ std::string nameTitleOuterTrackRes = " h2dOuterTimeResTrack" + particleNames2 [i_true] + " VsP" ;
203+ std::string nameTitleOuterTotalRes = " h2dOuterTimeResTotal" + particleNames2 [i_true] + " VsP" ;
204+ const AxisSpec axisTrackTimeRes{static_cast <int >(plotsConfig.nBinsTimeRes ), 0 .0f , +200 .0f , " Track time resolution - " + particleNames [i_true] + " (ps)" };
205+ const AxisSpec axisTotalTimeRes{static_cast <int >(plotsConfig.nBinsTimeRes ), 0 .0f , +200 .0f , " Total time resolution - " + particleNames [i_true] + " (ps)" };
206+ h2dInnerTimeResTrack[i_true] = histos.add <TH2>(nameTitleInnerTrackRes .c_str (), nameTitleInnerTrackRes .c_str (), kTH2F , {axisMomentum, axisTrackTimeRes});
207+ h2dInnerTimeResTotal[i_true] = histos.add <TH2>(nameTitleTotalRes .c_str (), nameTitleTotalRes .c_str (), kTH2F , {axisMomentum, axisTotalTimeRes});
208+ h2dOuterTimeResTrack[i_true] = histos.add <TH2>(nameTitleOuterTrackRes .c_str (), nameTitleOuterTrackRes .c_str (), kTH2F , {axisMomentum, axisTrackTimeRes});
209+ h2dOuterTimeResTotal[i_true] = histos.add <TH2>(nameTitleOuterTotalRes .c_str (), nameTitleOuterTotalRes .c_str (), kTH2F , {axisMomentum, axisTotalTimeRes});
210210 for (int i_hyp = 0 ; i_hyp < kParticles ; i_hyp++) {
211- std::string name_title_inner = " h2dInnerNsigmaTrue" + particle_names2 [i_true] + " Vs" + particle_names2 [i_hyp] + " Hypothesis" ;
212- std::string name_title_outer = " h2dOuterNsigmaTrue" + particle_names2 [i_true] + " Vs" + particle_names2 [i_hyp] + " Hypothesis" ;
213- std::string name_title_inner_delta = " h2dInnerDeltaTrue" + particle_names2 [i_true] + " Vs" + particle_names2 [i_hyp] + " Hypothesis" ;
214- std::string name_title_outer_delta = " h2dOuterDeltaTrue" + particle_names2 [i_true] + " Vs" + particle_names2 [i_hyp] + " Hypothesis" ;
211+ std::string nameTitle = " h2dInnerNsigmaTrue" + particleNames2 [i_true] + " Vs" + particleNames2 [i_hyp] + " Hypothesis" ;
212+ std::string nameTitleOuter = " h2dOuterNsigmaTrue" + particleNames2 [i_true] + " Vs" + particleNames2 [i_hyp] + " Hypothesis" ;
213+ std::string nameTitle_delta = " h2dInnerDeltaTrue" + particleNames2 [i_true] + " Vs" + particleNames2 [i_hyp] + " Hypothesis" ;
214+ std::string nameTitleOuter_delta = " h2dOuterDeltaTrue" + particleNames2 [i_true] + " Vs" + particleNames2 [i_hyp] + " Hypothesis" ;
215215 const AxisSpec axisX{plotsConfig.doSeparationVsPt .value ? axisPt : axisMomentum};
216216 if (i_true == i_hyp) {
217- const AxisSpec axisNsigmaCorrect{static_cast <int >(plotsConfig.nBinsNsigmaCorrectSpecies ), plotsConfig.minNsigmaRange , plotsConfig.maxNsigmaRange , " N#sigma - True " + particle_names1 [i_true] + " vs " + particle_names1 [i_hyp] + " hypothesis" };
218- h2dInnerNsigmaTrue[i_true][i_hyp] = histos.add <TH2>(name_title_inner .c_str (), name_title_inner .c_str (), kTH2F , {axisX, axisNsigmaCorrect});
219- h2dOuterNsigmaTrue[i_true][i_hyp] = histos.add <TH2>(name_title_outer .c_str (), name_title_outer .c_str (), kTH2F , {axisX, axisNsigmaCorrect});
217+ const AxisSpec axisNsigmaCorrect{static_cast <int >(plotsConfig.nBinsNsigmaCorrectSpecies ), plotsConfig.minNsigmaRange , plotsConfig.maxNsigmaRange , " N#sigma - True " + particleNames [i_true] + " vs " + particleNames [i_hyp] + " hypothesis" };
218+ h2dInnerNsigmaTrue[i_true][i_hyp] = histos.add <TH2>(nameTitle .c_str (), nameTitle .c_str (), kTH2F , {axisX, axisNsigmaCorrect});
219+ h2dOuterNsigmaTrue[i_true][i_hyp] = histos.add <TH2>(nameTitleOuter .c_str (), nameTitleOuter .c_str (), kTH2F , {axisX, axisNsigmaCorrect});
220220
221- const AxisSpec axisDeltaCorrect{static_cast <int >(plotsConfig.nBinsDeltaCorrectSpecies ), plotsConfig.minDeltaRange , plotsConfig.maxDeltaRange , " #Delta - True " + particle_names1 [i_true] + " vs " + particle_names1 [i_hyp] + " hypothesis" };
222- h2dInnerDeltaTrue[i_true][i_hyp] = histos.add <TH2>(name_title_inner_delta .c_str (), name_title_inner_delta .c_str (), kTH2F , {axisX, axisDeltaCorrect});
223- h2dOuterDeltaTrue[i_true][i_hyp] = histos.add <TH2>(name_title_outer_delta .c_str (), name_title_outer_delta .c_str (), kTH2F , {axisX, axisDeltaCorrect});
221+ const AxisSpec axisDeltaCorrect{static_cast <int >(plotsConfig.nBinsDeltaCorrectSpecies ), plotsConfig.minDeltaRange , plotsConfig.maxDeltaRange , " #Delta - True " + particleNames [i_true] + " vs " + particleNames [i_hyp] + " hypothesis" };
222+ h2dInnerDeltaTrue[i_true][i_hyp] = histos.add <TH2>(nameTitle_delta .c_str (), nameTitle_delta .c_str (), kTH2F , {axisX, axisDeltaCorrect});
223+ h2dOuterDeltaTrue[i_true][i_hyp] = histos.add <TH2>(nameTitleOuter_delta .c_str (), nameTitleOuter_delta .c_str (), kTH2F , {axisX, axisDeltaCorrect});
224224 } else {
225- const AxisSpec axisNsigmaWrong{static_cast <int >(plotsConfig.nBinsNsigmaWrongSpecies ), plotsConfig.minNsigmaRange , plotsConfig.maxNsigmaRange , " N#sigma - True " + particle_names1 [i_true] + " vs " + particle_names1 [i_hyp] + " hypothesis" };
226- h2dInnerNsigmaTrue[i_true][i_hyp] = histos.add <TH2>(name_title_inner .c_str (), name_title_inner .c_str (), kTH2F , {axisX, axisNsigmaWrong});
227- h2dOuterNsigmaTrue[i_true][i_hyp] = histos.add <TH2>(name_title_outer .c_str (), name_title_outer .c_str (), kTH2F , {axisX, axisNsigmaWrong});
225+ const AxisSpec axisNsigmaWrong{static_cast <int >(plotsConfig.nBinsNsigmaWrongSpecies ), plotsConfig.minNsigmaRange , plotsConfig.maxNsigmaRange , " N#sigma - True " + particleNames [i_true] + " vs " + particleNames [i_hyp] + " hypothesis" };
226+ h2dInnerNsigmaTrue[i_true][i_hyp] = histos.add <TH2>(nameTitle .c_str (), nameTitle .c_str (), kTH2F , {axisX, axisNsigmaWrong});
227+ h2dOuterNsigmaTrue[i_true][i_hyp] = histos.add <TH2>(nameTitleOuter .c_str (), nameTitleOuter .c_str (), kTH2F , {axisX, axisNsigmaWrong});
228228
229- const AxisSpec axisDeltaWrong{static_cast <int >(plotsConfig.nBinsDeltaWrongSpecies ), plotsConfig.minDeltaRange , plotsConfig.maxDeltaRange , " #Delta - True " + particle_names1 [i_true] + " vs " + particle_names1 [i_hyp] + " hypothesis" };
230- h2dInnerDeltaTrue[i_true][i_hyp] = histos.add <TH2>(name_title_inner_delta .c_str (), name_title_inner_delta .c_str (), kTH2F , {axisX, axisDeltaWrong});
231- h2dOuterDeltaTrue[i_true][i_hyp] = histos.add <TH2>(name_title_outer_delta .c_str (), name_title_outer_delta .c_str (), kTH2F , {axisX, axisDeltaWrong});
229+ const AxisSpec axisDeltaWrong{static_cast <int >(plotsConfig.nBinsDeltaWrongSpecies ), plotsConfig.minDeltaRange , plotsConfig.maxDeltaRange , " #Delta - True " + particleNames [i_true] + " vs " + particleNames [i_hyp] + " hypothesis" };
230+ h2dInnerDeltaTrue[i_true][i_hyp] = histos.add <TH2>(nameTitle_delta .c_str (), nameTitle_delta .c_str (), kTH2F , {axisX, axisDeltaWrong});
231+ h2dOuterDeltaTrue[i_true][i_hyp] = histos.add <TH2>(nameTitleOuter_delta .c_str (), nameTitleOuter_delta .c_str (), kTH2F , {axisX, axisDeltaWrong});
232232 }
233233 }
234234 }
0 commit comments