From ee33deccd65962a8f512cd78c55c0a00e0c363f3 Mon Sep 17 00:00:00 2001 From: Sean Dobbs Date: Wed, 7 Oct 2026 15:52:01 -0400 Subject: [PATCH 1/2] Implement energy cutoff flags in BCAL pi0 skim Declared but never used --- .../Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/plugins/Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc b/src/plugins/Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc index 2adf10c021..889abdca20 100644 --- a/src/plugins/Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc +++ b/src/plugins/Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc @@ -59,8 +59,8 @@ void JEventProcessor_pi0bcalskim::Init() { auto app = GetApplication(); - MIN_SH1_E = 0.2; - MIN_SH2_E = 0.2; + MIN_SH1_E = 0.4; + MIN_SH2_E = 0.4; WRITE_EVIO = 1; WRITE_HDDM = 0; @@ -228,8 +228,8 @@ void JEventProcessor_pi0bcalskim::Process(const std::shared_ptr& e TLorentzVector sh2_p(sh2_E*sh2_x/sh2_R,sh2_E*sh2_y/sh2_R,sh2_E*sh2_z/sh2_R,sh2_E); TLorentzVector ptot = sh1_p+sh2_p; inv_mass = ptot.M(); - Candidate |= ( (sh2_E>0.4) && (inv_mass<0.25) && (inv_mass>0.05)); - if((sh2_E>0.4) && (inv_mass<0.25) && (inv_mass>0.05)) { + Candidate |= ( (sh2_E>MIN_SH1_E) && (inv_mass<0.25) && (inv_mass>0.05)); + if((sh2_E>MIN_SH2_E) && (inv_mass<0.25) && (inv_mass>0.05)) { if(find(locObjectsToSave.begin(), locObjectsToSave.end(), locBCALShowers[i]) == locObjectsToSave.end()) locObjectsToSave.push_back(static_cast(locBCALShowers[i])); if(find(locObjectsToSave.begin(), locObjectsToSave.end(), locBCALShowers[j]) == locObjectsToSave.end()) From 902db900316bb40ff272cd891beaf4d7025f87e8 Mon Sep 17 00:00:00 2001 From: Sean Dobbs Date: Wed, 7 Oct 2026 18:19:39 -0400 Subject: [PATCH 2/2] bug fixes --- .../Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/plugins/Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc b/src/plugins/Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc index 889abdca20..797407d685 100644 --- a/src/plugins/Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc +++ b/src/plugins/Utilities/pi0bcalskim/JEventProcessor_pi0bcalskim.cc @@ -228,8 +228,8 @@ void JEventProcessor_pi0bcalskim::Process(const std::shared_ptr& e TLorentzVector sh2_p(sh2_E*sh2_x/sh2_R,sh2_E*sh2_y/sh2_R,sh2_E*sh2_z/sh2_R,sh2_E); TLorentzVector ptot = sh1_p+sh2_p; inv_mass = ptot.M(); - Candidate |= ( (sh2_E>MIN_SH1_E) && (inv_mass<0.25) && (inv_mass>0.05)); - if((sh2_E>MIN_SH2_E) && (inv_mass<0.25) && (inv_mass>0.05)) { + Candidate |= ( (sh1_E>MIN_SH1_E) && (sh2_E>MIN_SH2_E) && (inv_mass<0.25) && (inv_mass>0.05)); + if( (sh1_E>MIN_SH1_E) && (sh2_E>MIN_SH2_E) && (inv_mass<0.25) && (inv_mass>0.05)) { if(find(locObjectsToSave.begin(), locObjectsToSave.end(), locBCALShowers[i]) == locObjectsToSave.end()) locObjectsToSave.push_back(static_cast(locBCALShowers[i])); if(find(locObjectsToSave.begin(), locObjectsToSave.end(), locBCALShowers[j]) == locObjectsToSave.end())