@@ -820,25 +820,32 @@ struct FlowGfwLightIons {
820820 fFCptgenFull ->fillPtProfiles (centmult, rndm);
821821 fFCptgenFull ->fillCMProfiles (centmult, rndm);
822822 fFCptgenFull ->fillSubeventPtProfiles (centmult, rndm);
823+ fFCptgenFull ->fillCMSubeventProfiles (centmult, rndm);
823824 } else {
824825 fFCpt ->fillPtProfiles (centmult, rndm);
825826 fFCpt ->fillCMProfiles (centmult, rndm);
826827 fFCptFull ->fillPtProfiles (centmult, rndm);
827828 fFCptFull ->fillSubeventPtProfiles (centmult, rndm);
828829 fFCptFull ->fillCMProfiles (centmult, rndm);
830+ fFCptFull ->fillCMSubeventProfiles (centmult, rndm);
829831 }
830832 for (uint l_ind = 0 ; l_ind < corrconfigs.size (); ++l_ind) {
831833 if (!corrconfigs.at (l_ind).pTDif ) {
834+ uint8_t vnptmask = o2::analysis::gfw::configs.GetpTCorrMasks ()[l_ind];
832835 auto dnx = fGFW ->Calculate (corrconfigs.at (l_ind), 0 , kTRUE ).real ();
833- if (dnx == 0 )
836+ if (dnx == 0 ) {
837+ (dt == kGen ) ? fFCptgen ->skipVnPtProfiles (vnptmask) : fFCpt ->skipVnPtProfiles (vnptmask);
834838 continue ;
839+ }
835840 auto val = fGFW ->Calculate (corrconfigs.at (l_ind), 0 , kFALSE ).real () / dnx;
836841 if (std::abs (val) < 1 ) {
837842 (dt == kGen ) ? fFCgen ->FillProfile (corrconfigs.at (l_ind).Head .c_str (), centmult, val, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 , rndm) : fFC ->FillProfile (corrconfigs.at (l_ind).Head .c_str (), centmult, val, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 , rndm);
838- (dt == kGen ) ? fFCptgen ->fillVnPtProfiles (centmult, val, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 , rndm, o2::analysis::gfw::configs. GetpTCorrMasks ()[l_ind]) : fFCpt ->fillVnPtProfiles (centmult, val, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 , rndm, o2::analysis::gfw::configs. GetpTCorrMasks ()[l_ind] );
843+ (dt == kGen ) ? fFCptgen ->fillVnPtProfiles (centmult, val, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 , rndm, vnptmask) : fFCpt ->fillVnPtProfiles (centmult, val, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 , rndm, vnptmask );
839844 if (cfgRunByRun && cfgFillFlowRunByRun && dt != kGen && l_ind == 0 ) {
840845 tpfsList[run][pfCorr22]->Fill (centmult, val, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 );
841846 }
847+ } else {
848+ (dt == kGen ) ? fFCptgen ->skipVnPtProfiles (vnptmask) : fFCpt ->skipVnPtProfiles (vnptmask);
842849 }
843850 continue ;
844851 }
0 commit comments