@@ -973,19 +973,23 @@ struct MultiparticleCorrelationsMei // this name is used in lower-case format to
973973 } // end of for (int64_t nTrack = 0; nTrack < tracks.size(); ++nTrack) {
974974
975975 std::vector<TComplex> resultMultCorr (2 , TComplex (0 ., 0 .));
976- if (cfMasterCutSwitch && passEventCutsReal && passTechnicalCut) {
977- resultMultCorr = two (qVectorsTableReal, n2);
978- if (noneZeroDenom (resultMultCorr)) {
979- obs.fProfTwo [eRec]->Fill (thisCollCent, (resultMultCorr[0 ] / resultMultCorr[1 ].Re ()).Re ());
980- }
981-
982- if (rs == eRecAndSim && passEventCutsMC) {
983- resultMultCorr = two (qVectorsTableMC, n2);
976+ if (cfMasterCutSwitch) {
977+ if (passEventCutsReal && passTechnicalCut) {
978+ resultMultCorr = two (qVectorsTableReal, n2);
984979 if (noneZeroDenom (resultMultCorr)) {
985- obs.fProfTwo [eSim ]->Fill (thisCollCent, (resultMultCorr[0 ] / resultMultCorr[1 ].Re ()).Re ());
980+ obs.fProfTwo [eRec ]->Fill (thisCollCent, (resultMultCorr[0 ] / resultMultCorr[1 ].Re ()).Re ());
986981 }
987982 }
988- }
983+
984+ if constexpr (rs == eRecAndSim) {
985+ if (passEventCutsMC) {
986+ resultMultCorr = two (qVectorsTableMC, n2);
987+ if (noneZeroDenom (resultMultCorr)) {
988+ obs.fProfTwo [eSim]->Fill (thisCollCent, (resultMultCorr[0 ] / resultMultCorr[1 ].Re ()).Re ());
989+ }
990+ } // end of if (passEventCutsMC) {
991+ } // end of if constexpr (rs == eRecAndSim) {
992+ } // end of if (cfMasterCutSwitch) {
989993
990994 isFirstCollision = false ; // Now the first collision ends
991995 } // end of template <ERecSim rs, typename T1, typename T2> void steer(T1 const& collision, T2 const& tracks) {
0 commit comments