@@ -106,6 +106,7 @@ struct TrackEfficiency {
106106
107107 std::vector<int > eventSelectionBits;
108108 int trackSelection = -1 ;
109+ float pTHatSettingSentinelValue = 999.0 ;
109110
110111 enum AcceptSplitCollisionsOptions {
111112 NonSplitOnly = 0 ,
@@ -518,7 +519,7 @@ struct TrackEfficiency {
518519 }
519520 registry.fill (HIST (" hMcCollCutsCounts" ), 5.5 ); // at least one of the reconstructed collisions associated with this mcCollision is selected with regard to centrality
520521
521- float pTHat = mcCollision.ptHard () < 999 . 0f ? mcCollision.ptHard () : simPtRef / (std::pow (mcCollision.weight (), 1.0 / pTHatExponent));
522+ float pTHat = mcCollision.ptHard () < pTHatSettingSentinelValue ? mcCollision.ptHard () : simPtRef / (std::pow (mcCollision.weight (), 1.0 / pTHatExponent));
522523 if (pTHat < ptHatMin || pTHat > ptHatMax) { // only allows mcCollisions with weight in between min and max
523524 return ;
524525 }
@@ -861,7 +862,7 @@ struct TrackEfficiency {
861862 return ;
862863 }
863864
864- float pTHat = collision.mcCollision ().ptHard () < 999 . 0f ? collision.mcCollision ().ptHard () : simPtRef / (std::pow (collision.mcCollision ().weight (), 1.0 / pTHatExponent));
865+ float pTHat = collision.mcCollision ().ptHard () < pTHatSettingSentinelValue ? collision.mcCollision ().ptHard () : simPtRef / (std::pow (collision.mcCollision ().weight (), 1.0 / pTHatExponent));
865866 if (pTHat < ptHatMin || pTHat > ptHatMax) { // only allows mcCollisions with weight in between min and max
866867 return ;
867868 }
@@ -890,7 +891,7 @@ struct TrackEfficiency {
890891 return ;
891892 }
892893
893- float pTHat = collision.mcCollision ().ptHard () < 999 . 0f ? collision.mcCollision ().ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
894+ float pTHat = collision.mcCollision ().ptHard () < pTHatSettingSentinelValue ? collision.mcCollision ().ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
894895 if (pTHat < ptHatMin || pTHat > ptHatMax) { // only allows mcCollisions with weight in between min and max
895896 return ;
896897 }
@@ -915,7 +916,7 @@ struct TrackEfficiency {
915916 return ;
916917 }
917918
918- float pTHat = mcCollision.ptHard () < 999 . 0f ? mcCollision.ptHard () : simPtRef / (std::pow (mcCollision.weight (), 1.0 / pTHatExponent));
919+ float pTHat = mcCollision.ptHard () < pTHatSettingSentinelValue ? mcCollision.ptHard () : simPtRef / (std::pow (mcCollision.weight (), 1.0 / pTHatExponent));
919920 if (pTHat < ptHatMin || pTHat > ptHatMax) { // only allows mcCollisions with weight in between min and max
920921 return ;
921922 }
@@ -978,7 +979,7 @@ struct TrackEfficiency {
978979 return ;
979980 }
980981
981- float pTHat = mcCollision.ptHard () < 999 . 0f ? mcCollision.ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
982+ float pTHat = mcCollision.ptHard () < pTHatSettingSentinelValue ? mcCollision.ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
982983 if (pTHat < ptHatMin || pTHat > ptHatMax) { // only allows mcCollisions with weight in between min and max
983984 return ;
984985 }
@@ -1072,7 +1073,7 @@ struct TrackEfficiency {
10721073 registry.fill (HIST (" h_collisions" ), 3.5 );
10731074 registry.fill (HIST (" h2_centrality_collisions" ), centrality, 3.5 );
10741075
1075- float pTHat = collision.mcCollision ().ptHard () < 999 . 0f ? collision.mcCollision ().ptHard () : simPtRef / (std::pow (collision.mcCollision ().weight (), 1.0 / pTHatExponent));
1076+ float pTHat = collision.mcCollision ().ptHard () < pTHatSettingSentinelValue ? collision.mcCollision ().ptHard () : simPtRef / (std::pow (collision.mcCollision ().weight (), 1.0 / pTHatExponent));
10761077 if (pTHat < ptHatMin || pTHat > ptHatMax) { // only allows mcCollisions with weight in between min and max
10771078 return ;
10781079 }
@@ -1109,7 +1110,7 @@ struct TrackEfficiency {
11091110 registry.fill (HIST (" h_collisions" ), 3.5 );
11101111 registry.fill (HIST (" h_collisions_weighted" ), 3.5 , eventWeight);
11111112
1112- float pTHat = collision.mcCollision ().ptHard () < 999 . 0f ? collision.mcCollision ().ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
1113+ float pTHat = collision.mcCollision ().ptHard () < pTHatSettingSentinelValue ? collision.mcCollision ().ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
11131114 if (pTHat < ptHatMin || pTHat > ptHatMax) { // only allows mcCollisions with weight in between min and max
11141115 return ;
11151116 }
@@ -1124,7 +1125,7 @@ struct TrackEfficiency {
11241125 // float centrality = checkCentFT0M ? mcCollision.centFT0M() : mcCollision.centFT0C(); mcCollision.centFT0C() isn't filled at the moment; can be added back when it is
11251126
11261127 float eventWeight = mcCollision.weight ();
1127- float pTHat = mcCollision.ptHard () < 999 . 0f ? mcCollision.ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
1128+ float pTHat = mcCollision.ptHard () < pTHatSettingSentinelValue ? mcCollision.ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
11281129 registry.fill (HIST (" h2_mccollision_pthardfromweight_pthardfromhepmcxsection" ), simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent)), mcCollision.ptHard ());
11291130
11301131 float centrality = -1 ;
@@ -1196,7 +1197,7 @@ struct TrackEfficiency {
11961197 // float centrality = checkCentFT0M ? mcCollision.centFT0M() : mcCollision.centFT0C(); mcCollision.centFT0C() isn't filled at the moment; can be added back when it is
11971198
11981199 float eventWeight = mcCollision.weight ();
1199- float pTHat = mcCollision.ptHard () < 999 . 0f ? mcCollision.ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
1200+ float pTHat = mcCollision.ptHard () < pTHatSettingSentinelValue ? mcCollision.ptHard () : simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent));
12001201 registry.fill (HIST (" h2_mccollision_pthardfromweight_pthardfromhepmcxsection" ), simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent)), mcCollision.ptHard ());
12011202 registry.fill (HIST (" h2_mccollision_pthardfromweight_pthardfromhepmcxsection_weighted" ), simPtRef / (std::pow (eventWeight, 1.0 / pTHatExponent)), mcCollision.ptHard (), eventWeight);
12021203
0 commit comments