Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Scoring planes gone missing? #1213

Open
bryngemark opened this issue Sep 27, 2023 · 0 comments
Open

Scoring planes gone missing? #1213

bryngemark opened this issue Sep 27, 2023 · 0 comments

Comments

@bryngemark
Copy link
Contributor

Describe the bug
I don't see scoring planes in my recent inclusive electron sample. Actually this might be related to the recoil electron/EcalVeto problems I see and describe in #1212

To Reproduce
Follow the recipe in #1212; I attach the config here too for convenience

#!/bin/python                                                                                                                                       

import os
import sys
from LDMX.Framework import ldmxcfg
p=ldmxcfg.Process('sim')                                                                                   
p.run = 1
p.maxEvents = 10000
p.maxTriesPerEvent = 1

from LDMX.SimCore import generators as gen
from LDMX.SimCore import simulator
                                                                            
sim = simulator.simulator( "inclusive_single_e" )
sim.setDetector( 'ldmx-det-v14', True )  #True to include scoring planes                                                                    
sim.description = "One 4 GeV electron shot from far upstream."
sim.generators.append(gen.single_4gev_e_upstream_tagger())
sim.beamSpotSmear = [0.,0.,0.]
sim.description = '4 GeV inclusive pencil beam simulation'

import LDMX.Ecal.EcalGeometry
import LDMX.Ecal.ecal_hardcoded_conditions
import LDMX.Hcal.HcalGeometry
import LDMX.Hcal.hcal_hardcoded_conditions

from LDMX.Ecal import digi as eDigi
from LDMX.Ecal import vetos
                                                              
ecalReco   =eDigi.EcalRecProducer('ecalRecon')
ecalDigi = eDigi.EcalDigiProducer()
ecalReco.v14()

ecalVeto   =vetos.EcalVetoProcessor('ecalVetoBDT')
p.sequence=[ sim, ecalDigi, ecalReco, ecalVeto ] 

p.keep = [ "drop MagnetScoringPlaneHits", "drop TrackerScoringPlaneHits", "drop HcalScoringPlaneHits"]
p.outputFiles=["simoutput.root"]

Desired behavior
I'd like the scoring plane collections (except the ones I drop) to show up, and I thought the True in setDetector still did exactly this?

Environment:
Production image of v3.2.10 on LDCS

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
2 participants