Bug #16389
Backtracker fails on corsika in-time MC
0%
Description
Marco reported an issue with running backtracker on MicroBooNE corsika MC.
setup uboonecode v06_34_00 -q e14:prof
lar -c standard_ana_uboone.fcl /pnfs/uboone/mc/uboone/reconstructed/prodcosmics_corsika_cmc_uboone_intime_reco/reco/prod_intime_fix_v06_26_01/prodcosmics_corsika_cmc_uboone_intime_0_20170410T222014_gen2_13def821-c925-4deb-b6db-1eedbf20cb3e_20170419T053704_reco1_20170419T061659_reco2.root %MSG-s ArtException: PostProcessPath end_path 02-May-2017 20:01:06 CDT PostProcessEvent cet::exception caught in art ---- EventProcessorFailure BEGIN An exception occurred during current event processing ---- EventProcessorFailure BEGIN An exception occurred during current event processing ---- ScheduleExecutionFailure BEGIN ProcessingStopped. ---- ProductNotFound BEGIN A request to resolve an art::Ptr to a product containing items of type: simb::MCTruth with ProductID 1:9 cannot be satisfied because the product cannot be found. The productGetter was not set -- are you trying to dereference a Ptr during mixing? cet::exception going through module ---- ProductNotFound END Exception going through path end_path ---- ScheduleExecutionFailure END ---- EventProcessorFailure END ---- EventProcessorFailure END %MSG %MSG-s ArtException: PostProcessPath end_path 02-May-2017 20:01:06 CDT PostProcessEvent cet::exception caught in art ---- EventProcessorFailure BEGIN An exception occurred during current event processing ---- EventProcessorFailure BEGIN An exception occurred during current event processing ---- ScheduleExecutionFailure BEGIN ProcessingStopped. ---- ProductNotFound BEGIN A request to resolve an art::Ptr to a product containing items of type: simb::MCTruth with ProductID 1:9 cannot be satisfied because the product cannot be found. The productGetter was not set -- are you trying to dereference a Ptr during mixing? cet::exception going through module ---- ProductNotFound END Exception going through path end_path ---- ScheduleExecutionFailure END ---- EventProcessorFailure END ---- EventProcessorFailure END %MSG Art has completed and will exit with status 18.
Allinea shows the problem is at line 128 of BackTracker_service.cc:
for(size_t mc = 0; mc < fMCTruthList.size(); ++mc) LOG_DEBUG("BackTracker") << *(fMCTruthList.at(mc).get());
History
#1 Updated by Lynn Garren almost 4 years ago
- Status changed from New to Assigned
- Assignee set to Brian Rebel
Brian, would you have a quick look? There is also a possibility that there is an upstream bug in CORSIKAGen.
#2 Updated by Lynn Garren almost 4 years ago
- Assignee changed from Brian Rebel to Tingjun Yang
From Brian:
I suggest having whoever is trying to back track those files do the debugging.
Based on the error reported, I would guess that a MCTruth object is simply never made/filled for the Corsika events which would make them impossible to back track. I would have Tingjun or whoever is trying to use it verify that the MCTruth object is made and filled.
#3 Updated by Marco Del Tutto almost 4 years ago
- Subject changed from Backtracker fails on corsika MC to Backtracker fails on corsika in-time MC
Adding a comment:
I also see some track ID without an MCParticle. Talking to Wes, this may be due to the fact that Corsika in-time has two separate paths of geant4 simulation. I am not sure if this is related, but this may lead to a break in the geant-trackID to MCParticle relation and to MCTruth objects not correctly associated to the MCParticles in the event.
#4 Updated by Stefano Soleti over 3 years ago
I have a similar issue with v06_26_01_07
When I do
auto const &generator_handle =
evt.getValidHandle<std::vector<simb::MCTruth>>(_mctruthLabel);
I get
%MSG-s ArtException: PostPathEndRun end_path 11-Oct-2017 17:28:39 UTC PostEndRun
cet::exception caught in art
---- EventProcessorFailure BEGIN
An exception occurred during current event processing
---- EventProcessorFailure BEGIN
An exception occurred during current event processing
---- ScheduleExecutionFailure BEGIN
ProcessingStopped.
---- ProductNotFound BEGIN
A request to resolve an art::Ptr to a product containing items of type: simb::MCTruth with ProductID 1:9
cannot be satisfied because the product cannot be found.
The productGetter was not set -- are you trying to dereference a Ptr during mixing?
cet::exception going through module
---- ProductNotFound END
Exception going through path end_path
---- ScheduleExecutionFailure END
---- EventProcessorFailure END
cet::exception caught in EventProcessor and rethrown
---- EventProcessorFailure END
%MSG
Art has completed and will exit with status 65.
Would be nice to try to fix this, I am happy to help if you could give me guidance
#5 Updated by Marco Del Tutto over 3 years ago
This also affects calls to BuildMCParticleHitMaps in larpandora/LArPandoraInterface/LArPandoraHelper.cxx, as this function throws an exception if it founds a track ID without an MCParticle.
#6 Updated by Gianluca Petrillo over 3 years ago
- complete setup (
larsoft
anduboonecode
version) - configuration file being executed
- pointer to an input file to test on
- any special instruction if not just
lar -c config.fcl -s input.root
#7 Updated by Tingjun Yang over 3 years ago
Gianluca Petrillo wrote:
Could you publish here how to reproduce the problem?
- complete setup (
larsoft
anduboonecode
version)- configuration file being executed
- pointer to an input file to test on
- any special instruction if not just
lar -c config.fcl -s input.root
Hi Gianluca,
Marco's example can still be used to reproduce this problem:
setup uboonecode v06_34_00 -q e14:prof lar -c standard_ana_uboone.fcl /pnfs/uboone/mc/uboone/reconstructed/prodcosmics_corsika_cmc_uboone_intime_reco/reco/prod_intime_fix_v06_26_01/prodcosmics_corsika_cmc_uboone_intime_0_20170410T222014_gen2_13def821-c925-4deb-b6db-1eedbf20cb3e_20170419T053704_reco1_20170419T061659_reco2.root
Shall I reassign this ticket to you?
Thanks,
Tingjun
#8 Updated by Tingjun Yang over 3 years ago
- Assignee changed from Tingjun Yang to Wesley Ketchum
Reassign to Wes.