@@ -841,14 +841,15 @@ struct FlowGfwV02 {
841841 template <DataType dt>
842842 void fillOutputContainers (const float & centmult, const int & multiplicity, const double & rndm, const int & /* run*/ = 0 )
843843 {
844+ double threshold = 1.01 ;
844845 for (uint l_ind = 0 ; l_ind < corrconfigs.size (); ++l_ind) {
845846 if (!corrconfigs.at (l_ind).pTDif ) {
846847 auto dnx = fGFW ->Calculate (corrconfigs.at (l_ind), 0 , kTRUE ).real ();
847848 if (dnx == 0 )
848849 continue ;
849850 auto val = fGFW ->Calculate (corrconfigs.at (l_ind), 0 , kFALSE ).real () / dnx;
850851
851- if (std::abs (val) < 1 ) {
852+ if (std::abs (val) < threshold ) {
852853 fFC ->FillProfile (corrconfigs.at (l_ind).Head .c_str (), centmult, val, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 , rndm);
853854 }
854855 continue ;
@@ -871,7 +872,6 @@ struct FlowGfwV02 {
871872 ebyeWeight *= pidStates.hPtMid [PidCharged]->Integral ();
872873 }
873874 double ptFraction = 0 ;
874- double threshold = 1.01 ;
875875 int normIndex = (cfgNormalizeByCharged) ? PidCharged : pidInd; // Configured to normalize by charged particles or the selected particle
876876 if (pidStates.hPtMid [normIndex]->Integral () > 0 ) {
877877 ptFraction = pidStates.hPtMid [pidInd]->GetBinContent (i) / pidStates.hPtMid [normIndex]->Integral ();
@@ -1002,7 +1002,6 @@ struct FlowGfwV02 {
10021002 double ptFraction = 0 ;
10031003 if (pidStates.hPtMid [PidCharged]->Integral () > 0 ) {
10041004 ptFraction = pidStates.hPtMid [PidCharged]->GetBinContent (i) / pidStates.hPtMid [PidCharged]->Integral ();
1005- double threshold = 1.01 ;
10061005 if (std::abs (val) < threshold)
10071006 registry.fill (HIST (" v02pt" ), fSecondAxis ->GetBinCenter (i), centmult, multiplicity, val * ptFraction, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0 );
10081007 registry.fill (HIST (" nchMid" ), fSecondAxis ->GetBinCenter (i), centmult, multiplicity, ptFraction);
0 commit comments