Skip to content
Snippets Groups Projects
Commit 22ed77c9 authored by Jihee Kim's avatar Jihee Kim Committed by Sylvester Joosten
Browse files

Resolve "imaging_ecal_pion0 failing Job"

parent 3a26c550
Branches
No related tags found
1 merge request!160Resolve "imaging_ecal_pion0 failing Job"
...@@ -53,12 +53,12 @@ void gen_far_forward_protons(int n_events = 100, ...@@ -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"; //std::cout << std::sqrt(px*px + py*py + pz*pz) - p << " is zero? \n";
// type 1 is final state // 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>( GenParticlePtr p3 = std::make_shared<GenParticle>(
FourVector( FourVector(
px, py, pz, px, py, pz,
sqrt(p*p + (0.000511 * 0.000511))), sqrt(p*p + (0.938272 * 0.938272))),
11, 1); 2212, 1);
GenVertexPtr v1 = std::make_shared<GenVertex>(); GenVertexPtr v1 = std::make_shared<GenVertex>();
v1->add_particle_in(p1); v1->add_particle_in(p1);
......
...@@ -17,7 +17,7 @@ with open('config/emcal_barrel_calibration.json') as f: ...@@ -17,7 +17,7 @@ with open('config/emcal_barrel_calibration.json') as f:
calib_data = json.load(f)['electron'] calib_data = json.load(f)['electron']
kwargs = dict() kwargs = dict()
kwargs['sf'] = float(calib_data['sampling_fraction_scfi']) kwargs['img_sf'] = float(calib_data['sampling_fraction_img'])
# input arguments through environment variables # input arguments through environment variables
kwargs['input'] = os.environ.get('CB_EMCAL_SIM_FILE', '../topside/barrel_pion0_5GeV.root') 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', ...@@ -80,7 +80,7 @@ imcal_barrel_clreco = RecoCoG('imcal_barrel_clreco',
outputClusterCollection='EcalBarrelClusters', outputClusterCollection='EcalBarrelClusters',
outputInfoCollection='EcalBarrelClustersInfo', outputInfoCollection='EcalBarrelClustersInfo',
logWeightBase=6.2, logWeightBase=6.2,
samplingFraction=kwargs['sf']) samplingFraction=kwargs['img_sf'])
podout.outputCommands = ['keep *'] podout.outputCommands = ['keep *']
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment