From 22ed77c98f6fb9dcc1aa5ac0069d920671d97090 Mon Sep 17 00:00:00 2001
From: Jihee Kim <jihee.kim@anl.gov>
Date: Fri, 13 Aug 2021 14:52:10 +0000
Subject: [PATCH] Resolve "imaging_ecal_pion0 failing Job"
---
benchmarks/B0/scripts/gen_far_forward_protons.cxx | 6 +++---
benchmarks/imaging_ecal/options/imaging_2dcluster.py | 4 ++--
2 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/benchmarks/B0/scripts/gen_far_forward_protons.cxx b/benchmarks/B0/scripts/gen_far_forward_protons.cxx
index 33d5f2e4..1b363ca1 100644
--- a/benchmarks/B0/scripts/gen_far_forward_protons.cxx
+++ b/benchmarks/B0/scripts/gen_far_forward_protons.cxx
@@ -53,12 +53,12 @@ void gen_far_forward_protons(int n_events = 100,
//std::cout << std::sqrt(px*px + py*py + pz*pz) - p << " is zero? \n";
// type 1 is final state
- // pdgid 11 - proton 0.510 MeV/c^2
+ // pdgid 11 - proton 938.272 MeV/c^2
GenParticlePtr p3 = std::make_shared<GenParticle>(
FourVector(
px, py, pz,
- sqrt(p*p + (0.000511 * 0.000511))),
- 11, 1);
+ sqrt(p*p + (0.938272 * 0.938272))),
+ 2212, 1);
GenVertexPtr v1 = std::make_shared<GenVertex>();
v1->add_particle_in(p1);
diff --git a/benchmarks/imaging_ecal/options/imaging_2dcluster.py b/benchmarks/imaging_ecal/options/imaging_2dcluster.py
index 45669e17..6802e07d 100644
--- a/benchmarks/imaging_ecal/options/imaging_2dcluster.py
+++ b/benchmarks/imaging_ecal/options/imaging_2dcluster.py
@@ -17,7 +17,7 @@ with open('config/emcal_barrel_calibration.json') as f:
calib_data = json.load(f)['electron']
kwargs = dict()
-kwargs['sf'] = float(calib_data['sampling_fraction_scfi'])
+kwargs['img_sf'] = float(calib_data['sampling_fraction_img'])
# input arguments through environment variables
kwargs['input'] = os.environ.get('CB_EMCAL_SIM_FILE', '../topside/barrel_pion0_5GeV.root')
@@ -80,7 +80,7 @@ imcal_barrel_clreco = RecoCoG('imcal_barrel_clreco',
outputClusterCollection='EcalBarrelClusters',
outputInfoCollection='EcalBarrelClustersInfo',
logWeightBase=6.2,
- samplingFraction=kwargs['sf'])
+ samplingFraction=kwargs['img_sf'])
podout.outputCommands = ['keep *']
--
GitLab