@@ -155,8 +155,8 @@ struct AnalysisEnergyCorrelator {
155155 MixingHandler* fMixHandler = nullptr ;
156156
157157 AnalysisCompositeCut* fEventCut = nullptr ;
158- std::vector<AnalysisCompositeCut *> fTrackCuts ; // Electron cuts
159- std::vector<AnalysisCompositeCut *> fHadronCuts ; // Hadron cuts
158+ std::vector<AnalysisCut *> fTrackCuts ; // Electron cuts
159+ std::vector<AnalysisCut *> fHadronCuts ; // Hadron cuts
160160
161161 std::vector<TString> fTrackCutNames ;
162162 std::vector<TString> fHadronCutNames ;
@@ -182,12 +182,12 @@ struct AnalysisEnergyCorrelator {
182182 using MixingBinning = ColumnBinningPolicy<aod::collision::PosZ, aod::mult::MultNTracksPV>;
183183 std::unique_ptr<MixingBinning> fMixingBinning ;
184184
185- TH2F * hAcceptance_rec;
186- TH2F * hAcceptance_gen;
187- TH2F * hEfficiency_dilepton;
188- TH2F * hEfficiency_hadron;
189- TH1F * hMasswindow;
190- TH2F * hReweighthadron;
185+ TH2F * hAcceptance_rec = nullptr ;
186+ TH2F * hAcceptance_gen = nullptr ;
187+ TH2F * hEfficiency_dilepton = nullptr ;
188+ TH2F * hEfficiency_hadron = nullptr ;
189+ TH1F * hMasswindow = nullptr ;
190+ TH2F * hReweighthadron = nullptr ;
191191
192192 void init (o2::framework::InitContext& context)
193193 {
@@ -240,7 +240,7 @@ struct AnalysisEnergyCorrelator {
240240 if (trackCutsJSON != " " ) {
241241 std::vector<AnalysisCut*> addTrackCuts = dqcuts::GetCutsFromJSON (trackCutsJSON.Data ());
242242 for (auto & t : addTrackCuts) {
243- fTrackCuts .push_back (dynamic_cast <AnalysisCompositeCut*>(t) );
243+ fTrackCuts .push_back (t );
244244 fTrackCutNames .push_back (t->GetName ());
245245 trackCutStr += Form (" ,%s" , t->GetName ());
246246 }
@@ -313,7 +313,7 @@ struct AnalysisEnergyCorrelator {
313313 if (hadronCutsJSON != " " ) {
314314 std::vector<AnalysisCut*> addHadronCuts = dqcuts::GetCutsFromJSON (hadronCutsJSON.Data ());
315315 for (auto & t : addHadronCuts) {
316- fHadronCuts .push_back (dynamic_cast <AnalysisCompositeCut*>(t) );
316+ fHadronCuts .push_back (t );
317317 fHadronCutNames .push_back (t->GetName ());
318318 hadronCutStr += Form (" ,%s" , t->GetName ());
319319 }
@@ -651,26 +651,26 @@ struct AnalysisEnergyCorrelator {
651651 VarManager::FillPair<TPairType, TTrackFillMap>(t1, t2);
652652 VarManager::FillPairVertexing<TPairType, gkEventFillMapWithMults, TTrackFillMap>(event, t1, t2, VarManager::fgValues);
653653
654- if (sign1 * sign2 < 0 ) { // +- pairs
655- fHistMan ->FillHistClass (histNames[iEleCut][0 ].Data (), VarManager::fgValues); // reconstructed, unmatched
656- for (size_t isig = 0 ; isig < fRecMCSignals .size (); isig ++) { // loop over MC signals
657- if (mcDecision & (static_cast <uint32_t >(1 ) << isig )) {
658- fHistMan ->FillHistClass (histNamesMC[iEleCut * fRecMCSignals .size () + isig ][0 ].Data (), VarManager::fgValues); // matched signal
654+ if (sign1 * sign2 < 0 ) { // +- pairs
655+ fHistMan ->FillHistClass (histNames[iEleCut][0 ].Data (), VarManager::fgValues); // reconstructed, unmatched
656+ for (size_t mcRecSignal = 0 ; mcRecSignal < fRecMCSignals .size (); mcRecSignal ++) { // loop over MC signals
657+ if (mcDecision & (static_cast <uint32_t >(1 ) << mcRecSignal )) {
658+ fHistMan ->FillHistClass (histNamesMC[iEleCut * fRecMCSignals .size () + mcRecSignal ][0 ].Data (), VarManager::fgValues); // matched signal
659659 }
660660 }
661661 } else {
662662 if (sign1 > 0 ) { // ++ pairs
663663 fHistMan ->FillHistClass (histNames[iEleCut][1 ].Data (), VarManager::fgValues);
664- for (size_t isig = 0 ; isig < fRecMCSignals .size (); isig ++) { // loop over MC signals
665- if (mcDecision & (static_cast <uint32_t >(1 ) << isig )) {
666- fHistMan ->FillHistClass (histNamesMC[iEleCut * fRecMCSignals .size () + isig ][1 ].Data (), VarManager::fgValues);
664+ for (size_t mcRecSignal = 0 ; mcRecSignal < fRecMCSignals .size (); mcRecSignal ++) { // loop over MC signals
665+ if (mcDecision & (static_cast <uint32_t >(1 ) << mcRecSignal )) {
666+ fHistMan ->FillHistClass (histNamesMC[iEleCut * fRecMCSignals .size () + mcRecSignal ][1 ].Data (), VarManager::fgValues);
667667 }
668668 }
669669 } else { // -- pairs
670670 fHistMan ->FillHistClass (histNames[iEleCut][2 ].Data (), VarManager::fgValues);
671- for (unsigned int isig = 0 ; isig < fRecMCSignals .size (); isig ++) { // loop over MC signals
672- if (mcDecision & (static_cast <uint32_t >(1 ) << isig )) {
673- fHistMan ->FillHistClass (histNamesMC[iEleCut * fRecMCSignals .size () + isig ][2 ].Data (), VarManager::fgValues);
671+ for (unsigned int mcRecSignal = 0 ; mcRecSignal < fRecMCSignals .size (); mcRecSignal ++) { // loop over MC signals
672+ if (mcDecision & (static_cast <uint32_t >(1 ) << mcRecSignal )) {
673+ fHistMan ->FillHistClass (histNamesMC[iEleCut * fRecMCSignals .size () + mcRecSignal ][2 ].Data (), VarManager::fgValues);
674674 }
675675 }
676676 }
0 commit comments