Project

General

Profile

Bug #19996

Segfault in TCShower.cxx

Added by Tingjun Yang almost 3 years ago.

Status:
New
Priority:
Normal
Assignee:
Tingjun Yang
Category:
-
Start date:
05/20/2018
Due date:
% Done:

0%

Estimated time:
Duration:

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.

Also available in: Atom PDF