Bug #19996
Segfault in TCShower.cxx
Start date:
05/20/2018
Due date:
% Done:
0%
Estimated time:
Description
In TCShower.cxx:
// dE/dx, energy, etc for(auto cid : ss3.CotIDs) { auto& ss = tjs.cots[cid]; auto& stj = tjs.allTraj[ss.ShowerTjID - 1]; unsigned short plane = DecodeCTP(ss.CTP).Plane; ss3.Energy[plane] = ss.Energy; // TODO: calculate the errors in some better way ss3.EnergyErr[plane] = 0.3 * ss.Energy; // TODO: what does MIPEnergy mean anyway? ss3.MIPEnergy[plane] = ss3.EnergyErr[plane]; ss3.MIPEnergyErr[plane] = ss.Energy; std::cout<<plane<<" "<<ss3.dEdx.size()<<" "<<ss.ShowerTjID-1<<" "<<tjs.allTraj.size()<<std::endl; //Added by TY for debugging ss3.dEdx[plane] = stj.dEdx[0]; ss3.dEdxErr[plane] = 0.3 * stj.dEdx[0]; } // ci
For one event, tjs.allTraj.size() = 756, ss.ShowerTjID-1 = 1071706962
lar -c standard_reco_dune10kt_nu_1x2x6.fcl /pnfs/dune/scratch/users/vito/ci_validation/tracking_efficiency_validation/v06_75_00/sim/prodgenie_nu_dune10kt_1x2x6_detsim_021.root --nskip 33 -n 1
I am going to investigate a little.