Skip to content

Commit 03e06ef

Browse files
ddobrigkddobrigk
committed
Debug histograms added
1 parent dd986e1 commit 03e06ef

File tree

1 file changed

+18
-0
lines changed

1 file changed

+18
-0
lines changed

‎PWGLF/TableProducer/Strangeness/strangenesstofpid.cxx‎

Lines changed: 18 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -522,6 +522,16 @@ struct strangenesstofpid {
522522
histos.add("hDeltaTimeMethodsVsP_negOmPi", "hDeltaTimeMethodsVsP_negOmPi", kTH3F, {axisP, axisEta, axisDeltaTime});
523523
histos.add("hDeltaTimeMethodsVsP_negOmPr", "hDeltaTimeMethodsVsP_negOmPr", kTH3F, {axisP, axisEta, axisDeltaTime});
524524
histos.add("hDeltaTimeMethodsVsP_bachOmKa", "hDeltaTimeMethodsVsP_bachOmKa", kTH3F, {axisP, axisEta, axisDeltaTime});
525+
526+
histos.add("hMethodComparison_posLaPi", "hMethodComparison_posLaPi", kTH2F, {axisTime, axisTime});
527+
histos.add("hMethodComparison_posLaPr", "hMethodComparison_posLaPr", kTH2F, {axisTime, axisTime});
528+
histos.add("hMethodComparison_negLaPi", "hMethodComparison_negLaPi", kTH2F, {axisTime, axisTime});
529+
histos.add("hMethodComparison_negLaPr", "hMethodComparison_negLaPr", kTH2F, {axisTime, axisTime});
530+
531+
histos.add("hMethodComparison_posXiPi", "hMethodComparison_posXiPi", kTH2F, {axisTime, axisTime});
532+
histos.add("hMethodComparison_posXiPr", "hMethodComparison_posXiPr", kTH2F, {axisTime, axisTime});
533+
histos.add("hMethodComparison_negXiPi", "hMethodComparison_negXiPi", kTH2F, {axisTime, axisTime});
534+
histos.add("hMethodComparison_negXiPr", "hMethodComparison_negXiPr", kTH2F, {axisTime, axisTime});
525535
}
526536

527537
// list memory consumption at start if running in modes with more output
@@ -863,6 +873,7 @@ struct strangenesstofpid {
863873
(timeLambda + timePositivePr) / (pTra.tofSignal() - pTra.tofEvTime()),
864874
positiveP, v0.positiveeta());
865875
if(calculationMethod.value==2 && std::abs(timePositivePr_Method0-o2::aod::v0data::kNoTOFValue)>o2::aod::v0data::kEpsilon && std::abs(timePositivePr_Method1-o2::aod::v0data::kNoTOFValue)>o2::aod::v0data::kEpsilon){
876+
histos.fill(HIST("hMethodComparison_posLaPr"), timePositivePr_Method0, timePositivePr_Method1);
866877
histos.fill(HIST("hDeltaTimeMethodsVsP_posLaPr"), positiveP, v0.positiveeta(), (timePositivePr_Method0 - timePositivePr_Method1)*positiveCosine);
867878
}
868879
if (doQANSigma)
@@ -871,6 +882,7 @@ struct strangenesstofpid {
871882
if (std::abs(v0.mAntiLambda() - 1.115683) < v0Group.qaMassWindow && fabs(pTra.tpcNSigmaPi()) < v0Group.qaTPCNSigma && fabs(nTra.tpcNSigmaPr()) < v0Group.qaTPCNSigma) {
872883
histos.fill(HIST("h2dDeltaTimePositiveLambdaPi"), v0.p(), v0.eta(), deltaTimePositiveLambdaPi);
873884
if(calculationMethod.value==2 && std::abs(timePositivePi_Method0-o2::aod::v0data::kNoTOFValue)>o2::aod::v0data::kEpsilon && std::abs(timePositivePi_Method1-o2::aod::v0data::kNoTOFValue)>o2::aod::v0data::kEpsilon){
885+
histos.fill(HIST("hMethodComparison_posLaPi"), timePositivePi_Method0, timePositivePi_Method1);
874886
histos.fill(HIST("hDeltaTimeMethodsVsP_posLaPi"), positiveP, v0.positiveeta(), (timePositivePi_Method0 - timePositivePi_Method1)*positiveCosine);
875887
}
876888
if (doQANSigma)
@@ -895,6 +907,7 @@ struct strangenesstofpid {
895907
(timeLambda + timeNegativePi) / (nTra.tofSignal() - nTra.tofEvTime()),
896908
negativeP, v0.negativeeta());
897909
if(calculationMethod.value==2 && std::abs(timeNegativePi_Method0-o2::aod::v0data::kNoTOFValue)>o2::aod::v0data::kEpsilon && std::abs(timeNegativePi_Method1-o2::aod::v0data::kNoTOFValue)>o2::aod::v0data::kEpsilon){
910+
histos.fill(HIST("hMethodComparison_negLaPi"), timeNegativePi_Method0, timeNegativePi_Method1);
898911
histos.fill(HIST("hDeltaTimeMethodsVsP_negLaPi"), negativeP, v0.negativeeta(), (timeNegativePi_Method0 - timeNegativePi_Method1)*negativeCosine);
899912
}
900913
if (doQANSigma)
@@ -903,6 +916,7 @@ struct strangenesstofpid {
903916
if (std::abs(v0.mAntiLambda() - 1.115683) < v0Group.qaMassWindow && fabs(pTra.tpcNSigmaPi()) < v0Group.qaTPCNSigma && fabs(nTra.tpcNSigmaPr()) < v0Group.qaTPCNSigma) {
904917
histos.fill(HIST("h2dDeltaTimeNegativeLambdaPr"), v0.p(), v0.eta(), deltaTimeNegativeLambdaPr);
905918
if(calculationMethod.value==2 && std::abs(timeNegativePr_Method0-o2::aod::v0data::kNoTOFValue)>o2::aod::v0data::kEpsilon && std::abs(timeNegativePr_Method1-o2::aod::v0data::kNoTOFValue)>o2::aod::v0data::kEpsilon){
919+
histos.fill(HIST("hMethodComparison_negLaPr"), timeNegativePr_Method0, timeNegativePr_Method1);
906920
histos.fill(HIST("hDeltaTimeMethodsVsP_negLaPr"), negativeP, v0.negativeeta(), (timeNegativePr_Method0 - timeNegativePr_Method1)*negativeCosine);
907921
}
908922
if (doQANSigma)
@@ -1154,9 +1168,11 @@ struct strangenesstofpid {
11541168
histos.fill(HIST("h2dbachDeltaTimeAsXiPi"), cascade.p(), cascade.eta(), bachDeltaTimeAsXiPi);
11551169
if(calculationMethod.value==2){
11561170
if(std::abs(posFlightPr_Method0)>o2::aod::cascdata::kEpsilon && std::abs(posFlightPr_Method1)>o2::aod::cascdata::kEpsilon){
1171+
histos.fill(HIST("hMethodComparison_posXiPr"), posFlightPr_Method0, posFlightPr_Method1);
11571172
histos.fill(HIST("hDeltaTimeMethodsVsP_posXiPr"), positiveP, cascade.positiveeta(), (posFlightPr_Method0 - posFlightPr_Method1)*positiveCosine);
11581173
}
11591174
if(std::abs(negFlightPi_Method0)>o2::aod::cascdata::kEpsilon && std::abs(negFlightPi_Method1)>o2::aod::cascdata::kEpsilon){
1175+
histos.fill(HIST("hMethodComparison_negXiPi"), negFlightPi_Method0, negFlightPi_Method1);
11601176
histos.fill(HIST("hDeltaTimeMethodsVsP_negXiPi"), negativeP, cascade.negativeeta(), (negFlightPi_Method0 - negFlightPi_Method1)*negativeCosine);
11611177
}
11621178
if(std::abs(bachFlightPi_Method0)>o2::aod::cascdata::kEpsilon && std::abs(bachFlightPi_Method1)>o2::aod::cascdata::kEpsilon){
@@ -1197,9 +1213,11 @@ struct strangenesstofpid {
11971213
histos.fill(HIST("h2dbachDeltaTimeAsXiPi"), cascade.p(), cascade.eta(), bachDeltaTimeAsXiPi);
11981214
if(calculationMethod.value==2){
11991215
if(std::abs(posFlightPi_Method0-o2::aod::cascdata::kNoTOFValue)>o2::aod::cascdata::kEpsilon && std::abs(posFlightPi_Method1-o2::aod::cascdata::kNoTOFValue)>o2::aod::cascdata::kEpsilon){
1216+
histos.fill(HIST("hMethodComparison_posXiPr"), posFlightPi_Method0, posFlightPi_Method1);
12001217
histos.fill(HIST("hDeltaTimeMethodsVsP_posXiPi"), positiveP, cascade.positiveeta(), (posFlightPi_Method0 - posFlightPi_Method1)*positiveCosine);
12011218
}
12021219
if(std::abs(negFlightPr_Method0-o2::aod::cascdata::kNoTOFValue)>o2::aod::cascdata::kEpsilon && std::abs(negFlightPr_Method1-o2::aod::cascdata::kNoTOFValue)>o2::aod::cascdata::kEpsilon){
1220+
histos.fill(HIST("hMethodComparison_negXiPr"), negFlightPr_Method0, negFlightPr_Method1);
12031221
histos.fill(HIST("hDeltaTimeMethodsVsP_negXiPr"), negativeP, cascade.negativeeta(), (negFlightPr_Method0 - negFlightPr_Method1)*negativeCosine);
12041222
}
12051223
if(std::abs(bachFlightPi_Method0-o2::aod::cascdata::kNoTOFValue)>o2::aod::cascdata::kEpsilon && std::abs(bachFlightPi_Method1-o2::aod::cascdata::kNoTOFValue)>o2::aod::cascdata::kEpsilon){

0 commit comments

Comments
 (0)