@@ -1359,12 +1359,23 @@ struct Photonhbt {
13591359 float qside_lcms = q3_lcms.Dot (uv_side);
13601360 float qlong_lcms = q3_lcms.Dot (uv_long);
13611361 if (hbtanalysis.cfgDo3D ) {
1362- fRegistryPairMC .fill (HIST (" Pair/same/MC/NoLabel/CF_3D" ), std::fabs (qout_lcms), std::fabs (qside_lcms), std::fabs (qlong_lcms), kt);
1363- if (hbtanalysis.cfgDo2D ) {
1364- fRegistryPairMC .fill (HIST (" Pair/same/MC/NoLabel/CF_2D" ), std::fabs (qout_lcms), std::fabs (qinv), kt);
1362+ if constexpr (ev_id == 0 ) {
1363+ fRegistryCF .fill (HIST (" Pair/same/CF_3D" ), std::fabs (qout_lcms), std::fabs (qside_lcms), std::fabs (qlong_lcms), kt, weight);
1364+ if (hbtanalysis.cfgDo2D ) {
1365+ fRegistryCF .fill (HIST (" Pair/same/CF_2D" ), std::fabs (qout_lcms), std::fabs (qinv), kt, weight);
1366+ }
1367+ } else {
1368+ fRegistryCF .fill (HIST (" Pair/mix/CF_3D" ), std::fabs (qout_lcms), std::fabs (qside_lcms), std::fabs (qlong_lcms), kt, weight);
1369+ if (hbtanalysis.cfgDo2D ) {
1370+ fRegistryCF .fill (HIST (" Pair/mix/CF_2D" ), std::fabs (qout_lcms), std::fabs (qinv), kt, weight);
1371+ }
13651372 }
13661373 } else {
1367- fRegistryPairMC .fill (HIST (" Pair/same/MC/NoLabel/CF_1D" ), hbtanalysis.cfgUseLCMS ? qabs_lcms : qinv, kt);
1374+ if constexpr (ev_id == 0 ) {
1375+ fRegistryCF .fill (HIST (" Pair/same/CF_1D" ), hbtanalysis.cfgUseLCMS ? qabs_lcms : qinv, kt, weight);
1376+ } else {
1377+ fRegistryCF .fill (HIST (" Pair/mix/CF_1D" ), hbtanalysis.cfgUseLCMS ? qabs_lcms : qinv, kt, weight);
1378+ }
13681379 }
13691380 float deta_pair = v1.Eta () - v2.Eta ();
13701381 float dphi_pair = v1.Phi () - v2.Phi ();
@@ -2706,8 +2717,8 @@ struct Photonhbt {
27062717 }
27072718 }
27082719 } // end cfgDoTruthMix
2709- } // end collision loop
2710- } // end runTruthEfficiency
2720+ } // end collision loop
2721+ } // end runTruthEfficiency
27112722
27122723 /* ************************************************/
27132724 // PROCESS FUNCTIONS
0 commit comments