fix efficiency interpolation

Lucia Anna Tarasovicova · Lucia Anna Tarasovicova · commit 0fb4fefcb5d5 · 2026-03-02T07:32:16.000+01:00
diff --git a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx
@@ -419,6 +419,8 @@ struct HStrangeCorrelation {
     hPurityUncertaintyHadronMult = static_cast<TH2F*>(listEfficiencies->FindObject("hPurityUncertaintyHadronMult"));
     if (efficiencyFlags.applyEfficiencyPropagation && !efficiencyFlags.applyEffAsFunctionOfMultAndPhi && !hEfficiencyUncertaintyTrigger)
       LOG(fatal) << "Problem getting hEfficiencyUncertaintyTrigger!";
+    if (efficiencyFlags.applyEffAsFunctionOfMult && !hEfficiencyTriggerMult)
+      LOG(fatal) << "Problem getting hEfficiencyTriggerMult!";
     LOG(info) << "Efficiencies now loaded for " << mRunNumber;
   }
 
@@ -1173,7 +1175,7 @@ struct HStrangeCorrelation {
             } else {
               efficiency = hEfficiencyCascade[Index]->Interpolate(ptassoc, assoc.eta());
               if (efficiencyFlags.applyEfficiencyPropagation)
-                efficiencyError = hEfficiencyUncertaintyCascade[Index]->Interpolate(ptassoc, assoc.eta(), assoc.phi());
+                efficiencyError = hEfficiencyUncertaintyCascade[Index]->Interpolate(ptassoc, assoc.eta());
             }
           }
           if (efficiency == 0) { // check for zero efficiency, do not apply if the case

Original file line number	Diff line number	Diff line change
`@@ -419,6 +419,8 @@ struct HStrangeCorrelation {`
`419`	`419`	`hPurityUncertaintyHadronMult = static_cast<TH2F*>(listEfficiencies->FindObject("hPurityUncertaintyHadronMult"));`
`420`	`420`	`if (efficiencyFlags.applyEfficiencyPropagation && !efficiencyFlags.applyEffAsFunctionOfMultAndPhi && !hEfficiencyUncertaintyTrigger)`
`421`	`421`	`LOG(fatal) << "Problem getting hEfficiencyUncertaintyTrigger!";`
	`422`	`+ if (efficiencyFlags.applyEffAsFunctionOfMult && !hEfficiencyTriggerMult)`
	`423`	`+ LOG(fatal) << "Problem getting hEfficiencyTriggerMult!";`
`422`	`424`	`LOG(info) << "Efficiencies now loaded for " << mRunNumber;`
`423`	`425`	`}`
`424`	`426`
`@@ -1173,7 +1175,7 @@ struct HStrangeCorrelation {`
`1173`	`1175`	`} else {`
`1174`	`1176`	`efficiency = hEfficiencyCascade[Index]->Interpolate(ptassoc, assoc.eta());`
`1175`	`1177`	`if (efficiencyFlags.applyEfficiencyPropagation)`
`1176`		`- efficiencyError = hEfficiencyUncertaintyCascade[Index]->Interpolate(ptassoc, assoc.eta(), assoc.phi());`
	`1178`	`+ efficiencyError = hEfficiencyUncertaintyCascade[Index]->Interpolate(ptassoc, assoc.eta());`
`1177`	`1179`	`}`
`1178`	`1180`	`}`
`1179`	`1181`	`if (efficiency == 0) { // check for zero efficiency, do not apply if the case`