Skip to content

Commit 718319b

Browse files
author
Lucia Anna Tarasovicova
committed
fix the trigger eff for mixing
1 parent 32e844a commit 718319b

File tree

1 file changed

+61
-61
lines changed

1 file changed

+61
-61
lines changed

‎PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx‎

Lines changed: 61 additions & 61 deletions
Original file line numberDiff line numberDiff line change
@@ -612,21 +612,22 @@ struct HStrangeCorrelation {
612612
float efficiencyTriggError = 0.0f;
613613
float purityTrigg = 1.0f;
614614
float purityTriggErr = 0.0;
615-
if (!mixing) {
616-
if (efficiencyFlags.applyEfficiencyForTrigger) {
617-
efficiencyTrigg = hEfficiencyTrigger->Interpolate(trigg.pt(), trigg.eta());
615+
if (efficiencyFlags.applyEfficiencyForTrigger) {
616+
efficiencyTrigg = hEfficiencyTrigger->Interpolate(trigg.pt(), trigg.eta());
617+
if (efficiencyFlags.applyPurityTrigger)
618+
purityTrigg = hPurityHadron->Interpolate(trigg.pt());
619+
if (efficiencyFlags.applyEfficiencyPropagation) {
620+
efficiencyTriggError = hEfficiencyUncertaintyTrigger->Interpolate(trigg.pt(), trigg.eta());
618621
if (efficiencyFlags.applyPurityTrigger)
619-
purityTrigg = hPurityHadron->Interpolate(trigg.pt());
620-
if (efficiencyFlags.applyEfficiencyPropagation) {
621-
efficiencyTriggError = hEfficiencyUncertaintyTrigger->Interpolate(trigg.pt(), trigg.eta());
622-
if (efficiencyFlags.applyPurityTrigger)
623-
purityTriggErr = hPurityHadron->Interpolate(trigg.pt());
624-
}
625-
if (efficiencyTrigg == 0) { // check for zero efficiency, do not apply if the case
626-
efficiencyTrigg = 1;
627-
efficiencyTriggError = 0;
628-
}
622+
purityTriggErr = hPurityHadron->Interpolate(trigg.pt());
623+
}
624+
if (efficiencyTrigg == 0) { // check for zero efficiency, do not apply if the case
625+
efficiencyTrigg = 1;
626+
efficiencyTriggError = 0;
629627
}
628+
}
629+
if (!mixing) {
630+
630631
fillTriggerHistogram(histos.get<TH2>(HIST("sameEvent/TriggerParticlesV0")), trigg.pt(), mult, efficiencyTrigg, efficiencyTriggError, purityTrigg, purityTriggErr);
631632
}
632633

@@ -816,36 +817,36 @@ struct HStrangeCorrelation {
816817
float efficiencyTriggError = 0.0f;
817818
float purityTrigg = 1.0f;
818819
float purityTriggErr = 0.0f;
819-
if (!mixing) {
820-
if (efficiencyFlags.applyEfficiencyForTrigger) {
821-
if (efficiencyFlags.applyEffAsFunctionOfMult) {
822-
efficiencyTrigg = hEfficiencyTriggerMult->Interpolate(trigg.pt(), trigg.eta(), mult);
823-
} else {
824-
efficiencyTrigg = hEfficiencyTrigger->Interpolate(trigg.pt(), trigg.eta());
825-
}
820+
if (efficiencyFlags.applyEfficiencyForTrigger) {
821+
if (efficiencyFlags.applyEffAsFunctionOfMult) {
822+
efficiencyTrigg = hEfficiencyTriggerMult->Interpolate(trigg.pt(), trigg.eta(), mult);
823+
} else {
824+
efficiencyTrigg = hEfficiencyTrigger->Interpolate(trigg.pt(), trigg.eta());
825+
}
826+
if (efficiencyFlags.applyPurityTrigger) {
827+
if (efficiencyFlags.applyEffAsFunctionOfMult)
828+
purityTrigg = hPurityHadronMult->Interpolate(trigg.pt(), mult);
829+
else
830+
purityTrigg = hPurityHadron->Interpolate(trigg.pt());
831+
}
832+
if (efficiencyFlags.applyEfficiencyPropagation) {
833+
if (efficiencyFlags.applyEffAsFunctionOfMult)
834+
efficiencyTriggError = hEfficiencyUncertaintyTriggerMult->Interpolate(trigg.pt(), trigg.eta(), mult);
835+
else
836+
efficiencyTriggError = hEfficiencyUncertaintyTrigger->Interpolate(trigg.pt(), trigg.eta());
826837
if (efficiencyFlags.applyPurityTrigger) {
827838
if (efficiencyFlags.applyEffAsFunctionOfMult)
828-
purityTrigg = hPurityHadronMult->Interpolate(trigg.pt(), mult);
839+
purityTriggErr = hPurityUncertaintyHadronMult->Interpolate(trigg.pt(), mult);
829840
else
830-
purityTrigg = hPurityHadron->Interpolate(trigg.pt());
831-
}
832-
if (efficiencyFlags.applyEfficiencyPropagation) {
833-
if (efficiencyFlags.applyEffAsFunctionOfMult)
834-
efficiencyTriggError = hEfficiencyUncertaintyTriggerMult->Interpolate(trigg.pt(), trigg.eta(), mult);
835-
else
836-
efficiencyTriggError = hEfficiencyUncertaintyTrigger->Interpolate(trigg.pt(), trigg.eta());
837-
if (efficiencyFlags.applyPurityTrigger) {
838-
if (efficiencyFlags.applyEffAsFunctionOfMult)
839-
purityTriggErr = hPurityUncertaintyHadronMult->Interpolate(trigg.pt(), mult);
840-
else
841-
purityTriggErr = hPurityUncertaintyHadron->Interpolate(trigg.pt());
842-
}
843-
}
844-
if (efficiencyTrigg == 0) { // check for zero efficiency, do not apply if the case
845-
efficiencyTrigg = 1;
846-
efficiencyTriggError = 0;
841+
purityTriggErr = hPurityUncertaintyHadron->Interpolate(trigg.pt());
847842
}
848843
}
844+
if (efficiencyTrigg == 0) { // check for zero efficiency, do not apply if the case
845+
efficiencyTrigg = 1;
846+
efficiencyTriggError = 0;
847+
}
848+
}
849+
if (!mixing) {
849850
fillTriggerHistogram(histos.get<TH2>(HIST("sameEvent/TriggerParticlesCascade")), trigg.pt(), mult, efficiencyTrigg, efficiencyTriggError, purityTrigg, purityTriggErr);
850851
}
851852
double triggSign = trigg.sign();
@@ -1008,36 +1009,35 @@ struct HStrangeCorrelation {
10081009
float efficiencyTriggerError = 0.0f;
10091010
float purityTrigger = 1.0f;
10101011
float purityTriggerError = 0.0f;
1011-
if (!mixing) {
1012-
1013-
if (efficiencyFlags.applyEfficiencyForTrigger) {
1012+
if (efficiencyFlags.applyEfficiencyForTrigger) {
1013+
if (efficiencyFlags.applyEffAsFunctionOfMult)
1014+
efficiencyTrigger = hEfficiencyTriggerMult->Interpolate(trigg.pt(), trigg.eta(), mult);
1015+
else
1016+
efficiencyTrigger = hEfficiencyTrigger->Interpolate(trigg.pt(), trigg.eta());
1017+
if (efficiencyFlags.applyPurityTrigger) {
10141018
if (efficiencyFlags.applyEffAsFunctionOfMult)
1015-
efficiencyTrigger = hEfficiencyTriggerMult->Interpolate(trigg.pt(), trigg.eta(), mult);
1019+
purityTrigger = hPurityHadronMult->Interpolate(trigg.pt(), mult);
10161020
else
1017-
efficiencyTrigger = hEfficiencyTrigger->Interpolate(trigg.pt(), trigg.eta());
1021+
purityTrigger = hPurityHadron->Interpolate(trigg.pt());
1022+
}
1023+
if (efficiencyFlags.applyEfficiencyPropagation) {
1024+
if (efficiencyFlags.applyEffAsFunctionOfMult)
1025+
efficiencyTriggerError = hEfficiencyUncertaintyTriggerMult->Interpolate(trigg.pt(), trigg.eta(), mult);
1026+
else
1027+
efficiencyTriggerError = hEfficiencyUncertaintyTrigger->Interpolate(trigg.pt(), trigg.eta());
10181028
if (efficiencyFlags.applyPurityTrigger) {
10191029
if (efficiencyFlags.applyEffAsFunctionOfMult)
1020-
purityTrigger = hPurityHadronMult->Interpolate(trigg.pt(), mult);
1030+
purityTriggerError = hPurityUncertaintyHadronMult->Interpolate(trigg.pt(), mult);
10211031
else
1022-
purityTrigger = hPurityHadron->Interpolate(trigg.pt());
1023-
}
1024-
if (efficiencyFlags.applyEfficiencyPropagation) {
1025-
if (efficiencyFlags.applyEffAsFunctionOfMult)
1026-
efficiencyTriggerError = hEfficiencyUncertaintyTriggerMult->Interpolate(trigg.pt(), trigg.eta(), mult);
1027-
else
1028-
efficiencyTriggerError = hEfficiencyUncertaintyTrigger->Interpolate(trigg.pt(), trigg.eta());
1029-
if (efficiencyFlags.applyPurityTrigger) {
1030-
if (efficiencyFlags.applyEffAsFunctionOfMult)
1031-
purityTriggerError = hPurityUncertaintyHadronMult->Interpolate(trigg.pt(), mult);
1032-
else
1033-
purityTriggerError = hPurityUncertaintyHadron->Interpolate(trigg.pt());
1034-
}
1035-
}
1036-
if (efficiencyTrigger == 0) { // check for zero efficiency, do not apply if the case
1037-
efficiencyTrigger = 1;
1038-
efficiencyTriggerError = 0;
1032+
purityTriggerError = hPurityUncertaintyHadron->Interpolate(trigg.pt());
10391033
}
10401034
}
1035+
if (efficiencyTrigger == 0) { // check for zero efficiency, do not apply if the case
1036+
efficiencyTrigger = 1;
1037+
efficiencyTriggerError = 0;
1038+
}
1039+
}
1040+
if (!mixing) {
10411041
if constexpr (requires { triggerTrack.extra(); })
10421042
fillTriggerHistogram(histos.get<TH2>(HIST("sameEvent/TriggerParticlesPion")), trigg.pt(), mult, efficiencyTrigger, efficiencyTriggerError, purityTrigger, purityTriggerError);
10431043
else

0 commit comments

Comments
 (0)