-
Notifications
You must be signed in to change notification settings - Fork 45
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge pull request #192 from SBNSoftware/mdeltutt_issue147
Cosmic InTime Workflow with Refactored LArG4
- Loading branch information
Showing
7 changed files
with
185 additions
and
25 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
127 changes: 116 additions & 11 deletions
127
sbndcode/JobConfigurations/standard/g4/g4_simphotontime_filter.fcl
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1,16 +1,121 @@ | ||
| # File: g4_simphotontime_filter.fcl | ||
| # Purpose: A geant4 fcl, supposed to run after prodcorsika_proton_intime_filter.fcl, or similar, | ||
| # for making cosmics in time samples. | ||
| # | ||
| # This runs the new, refactored, LArG4 simulation. | ||
|
|
||
| #include "mergesimsources_sbnd.fcl" | ||
|
|
||
| #include "standard_g4_sbnd.fcl" | ||
|
|
||
| physics.producers.larg4outtime: @local::physics.producers.largeant | ||
| physics.producers.larg4outtime.InputLabels: [ "GenInTimeSorter:outtime" ] | ||
| # Add process to finish the electron drift simulation for the intimes | ||
| physics.producers.simdriftintime: @local::sbnd_simdrift | ||
|
|
||
| # Add a geant4 process, that will run only on the outtime cosmics | ||
| physics.producers.loader: { module_type: "PhysListLoader" } | ||
| physics.producers.larg4outtime: @local::sbnd_larg4 | ||
| physics.producers.ionandscintouttime: @local::sbnd_ionandscint | ||
| physics.producers.pdfastsimouttime: @local::sbnd_pdfastsim_par | ||
| physics.producers.simdriftouttime: @local::sbnd_simdrift | ||
|
|
||
| # Set the appropriate input labels, to run geant4 only on the outtime cosmics | ||
| physics.producers.simdriftintime.SimulationLabel: "ionandscintintime:priorSCE" | ||
| physics.producers.larg4outtime.inputCollections: [ "GenInTimeSorter:outtime" ] | ||
| physics.producers.ionandscintouttime.InputModuleLabels: ["larg4outtime"] | ||
| physics.producers.pdfastsimouttime.SimulationLabel: "ionandscintouttime:priorSCE" | ||
| physics.producers.simdriftouttime.SimulationLabel: "ionandscintouttime:priorSCE" | ||
|
|
||
| # Add processes for light simulation outside the active volume (AV) for the intimes | ||
| physics.producers.ionandscintoutintime: @local::sbnd_ionandscint_out | ||
| physics.producers.pdfastsimoutintime: @local::sbnd_pdfastsim_pvs | ||
| physics.producers.ionandscintoutintime.InputModuleLabels: ["larg4intime"] | ||
| physics.producers.pdfastsimoutintime.SimulationLabel: "ionandscintoutintime" | ||
|
|
||
| # Add processes for light simulation outside the active volume (AV) for the outtimes | ||
| physics.producers.ionandscintoutouttime: @local::sbnd_ionandscint_out | ||
| physics.producers.pdfastsimoutouttime: @local::sbnd_pdfastsim_pvs | ||
| physics.producers.ionandscintoutouttime.InputModuleLabels: ["larg4outtime"] | ||
| physics.producers.pdfastsimoutouttime.SimulationLabel: "ionandscintoutouttime" | ||
|
|
||
| # Add a process that merges the MCParticles | ||
| physics.producers.largeant: @local::sbnd_merge_sim_sources | ||
| physics.producers.largeant.FillMCParticles: true | ||
| physics.producers.largeant.InputSourcesLabels: [ "larg4intime", "larg4outtime"] | ||
|
|
||
| # Add a process that merges the SimEnergyDeposits | ||
| physics.producers.ionandscint: @local::sbnd_merge_sim_sources | ||
| physics.producers.ionandscint.FillSimEnergyDeposits: true | ||
| physics.producers.ionandscint.InputSourcesLabels: [ "ionandscintintime", "ionandscintouttime"] | ||
|
|
||
| # Add a process that merges the SimChannels | ||
| physics.producers.simdrift: @local::sbnd_merge_sim_sources | ||
| physics.producers.simdrift.FillSimChannels: true | ||
| physics.producers.simdrift.InputSourcesLabels: [ "simdriftintime", "simdriftouttime"] | ||
|
|
||
| # Add a process that merges the AuxDetSimChannels TODO | ||
| physics.producers.crtsimch: @local::sbnd_merge_sim_sources | ||
| physics.producers.crtsimch.FillAuxDetSimChannels: true | ||
| physics.producers.crtsimch.InputSourcesLabels: [ "?", "?"] | ||
|
|
||
| # Add a process that merges the SimPhotons inside the AV | ||
| physics.producers.pdfastsim: @local::sbnd_merge_sim_sources | ||
| physics.producers.pdfastsim.FillSimPhotons: true | ||
| physics.producers.pdfastsim.InputSourcesLabels: [ "pdfastsimintime", "pdfastsimouttime"] | ||
|
|
||
| # Add a process that merges the SimPhotons outside the AV | ||
| physics.producers.pdfastsimout: @local::sbnd_merge_sim_sources | ||
| physics.producers.pdfastsimout.FillSimPhotons: true | ||
| physics.producers.pdfastsimout.InputSourcesLabels: [ "pdfastsimoutintime", "pdfastsimoutouttime"] | ||
|
|
||
| # Add all these new modules to the simulate path | ||
| physics.simulate: [ rns | ||
| ### Complete intime drift simulation | ||
| , simdriftintime | ||
| ### Do full Geant4 simulation for the outtimes | ||
| , loader | ||
| , larg4outtime | ||
| , ionandscintouttime | ||
| , pdfastsimouttime | ||
| , simdriftouttime | ||
| ### Simulate the light outside the AV | ||
| , ionandscintoutintime | ||
| , pdfastsimoutintime | ||
| , ionandscintoutouttime | ||
| , pdfastsimoutouttime | ||
| ### Merge the intime and outtime paths | ||
| , largeant | ||
| , ionandscint | ||
| , simdrift | ||
| , pdfastsim | ||
| , pdfastsimout | ||
| ### Do truth-level reconstruction | ||
| , mcreco | ||
| ] | ||
|
|
||
| physics.producers.largeant: { | ||
| module_type: "MergeSimSources" | ||
| InputSourcesLabels: [ "larg4intime","larg4outtime"] | ||
| TrackIDOffsets: [ 10000000,20000000 ] | ||
| } | ||
| services.ParticleListAction.keepGenTrajectories: ["GenInTimeSorter"] | ||
|
|
||
| physics.simulate: [ rns, larg4outtime, largeant, mcreco ] | ||
| # Drop the intime and outtime collections, which have now been | ||
| # been merged into a 'largeant' collection | ||
| outputs.out1.outputCommands: [ "keep *_*_*_*" | ||
| # Drop G4 | ||
| , "drop *_larg4intime_*_*" | ||
| , "drop *_larg4outtime_*_*" | ||
| # Drop IonAndScint Inside AV | ||
| , "drop *_ionandscintintime_*_*" | ||
| , "drop *_ionandscintouttime_*_*" | ||
| # Drop PDFastSim Inside AV | ||
| , "drop *_pdfastsimintime_*_*" | ||
| , "drop *_pdfastsimouttime_*_*" | ||
| # Drop SimDrift Inside AV | ||
| , "drop *_simdriftintime_*_*" | ||
| , "drop *_simdriftouttime_*_*" | ||
| # Drop IonAndScint Outside AV | ||
| , "drop *_ionandscintoutintime_*_*" | ||
| , "drop *_ionandscintoutouttime_*_*" | ||
| # Drop PDFastSim Ouside AV | ||
| , "drop *_pdfastsimoutintime_*_*" | ||
| , "drop *_pdfastsimoutouttime_*_*" | ||
| ] | ||
|
|
||
| outputs.out1.outputCommands: [ "keep *_*_*_*", | ||
| "drop *_larg4intime_*_*", | ||
| "drop *_larg4outtime_*_*"] | ||
| # Remove unnecesary processes | ||
| physics.producers.ionandscintout: @erase |
58 changes: 48 additions & 10 deletions
58
sbndcode/JobConfigurations/standard/gen/corsika/prodcorsika_proton_intime_filter.fcl
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1,28 +1,66 @@ | ||
| # File: prodcorsika_proton_intime_filter.fcl | ||
| # Purpose: Generates CORSIKA events filtering out events with no cosmics intime | ||
| # | ||
| # This runs the new, refactored, LArG4 simulation. | ||
| # | ||
| # Ported from uBooNE gen-in-time fhicl for use by SBND | ||
| # by Gray Putnam <[email protected]> | ||
|
|
||
| #include "filtersgenintime_sbnd.fcl" | ||
| #include "filterssimphotonlitetime_sbnd.fcl" | ||
| #include "largeantmodules_sbnd.fcl" | ||
|
|
||
| #include "larg4_sbnd.fcl" | ||
| #include "ionandscint_sbnd.fcl" | ||
| #include "PDFastSim_sbnd.fcl" | ||
| #include "simdrift_sbnd.fcl" | ||
|
|
||
| #include "prodcorsika_cosmics_proton.fcl" | ||
|
|
||
| # Ported from uBooNE gen-in-time fhicl for use by SBND | ||
| # by Gray Putnam <[email protected]> | ||
|
|
||
| # Add the g4 services | ||
| services: { | ||
| @table::services | ||
| @table::sbnd_g4_services | ||
| } | ||
|
|
||
| physics.producers.larg4intime: @local::sbnd_largeant | ||
| # Rename generator to corsika | ||
| physics.producers.corsika: @local::physics.producers.generator | ||
|
|
||
| # Add a generation filter, that splits in intime and outtime cosmics | ||
| physics.filters.GenInTimeSorter: @local::sbnd_filtergenintime | ||
|
|
||
| # Add a geant4 process, only on the intime cosmics | ||
| physics.producers.loader: { module_type: "PhysListLoader" } | ||
| physics.producers.larg4intime: @local::sbnd_larg4 | ||
| physics.producers.ionandscintintime: @local::sbnd_ionandscint | ||
| physics.producers.pdfastsimintime: @local::sbnd_pdfastsim_par | ||
|
|
||
| # Add a filter on the geant4 intime output, based on sim photons | ||
| physics.filters.timefilter: @local::sbnd_timefilterssimphotonlitetime | ||
|
|
||
| physics.producers.corsika: @local::physics.producers.generator | ||
| # Add all these new modules to the simulate path | ||
| physics.simulate: [ corsika | ||
| , GenInTimeSorter | ||
| , loader | ||
| , larg4intime | ||
| , ionandscintintime | ||
| , pdfastsimintime | ||
| , timefilter | ||
| , rns | ||
| ] | ||
|
|
||
| # Remove unnecesary processes | ||
| physics.producers.generator: @erase | ||
|
|
||
| # Set the appropriate input labels | ||
| physics.producers.larg4intime.inputCollections: ["GenInTimeSorter:intime"] | ||
| physics.producers.ionandscintintime.InputModuleLabels: ["larg4intime"] | ||
| physics.producers.pdfastsimintime.SimulationLabel: "ionandscintintime:priorSCE" | ||
| physics.filters.timefilter.SimPhotonsLiteCollectionLabel: "pdfastsimintime" | ||
|
|
||
| services.ParticleListAction.keepGenTrajectories: ["GenInTimeSorter"] | ||
|
|
||
| physics.simulate: [ corsika, GenInTimeSorter, larg4intime, timefilter, rns ] | ||
| outputs.out1.SelectEvents: [ "simulate" ] | ||
|
|
||
| physics.producers.larg4intime.KeepParticlesInVolumes: ["volCryostat", "volTaggerTopHigh", "volTaggerTopLow", "volTaggerSideLeft", "volTaggerSideRight", "volTaggerFaceFront", "volTaggerFaceBack", "volTaggerBot"] | ||
| physics.producers.larg4intime.InputLabels: [ "GenInTimeSorter:intime" ] | ||
| #physics.producers.larg4intime.SparsifyTrajectory: "true" | ||
| process_name: CosmicsCorsikaProtonGenAndG4InTime | ||
|
|
||
| process_name: CosmicsCorsikaCMCGenAndG4InTime |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,17 @@ | ||
|
|
||
| BEGIN_PROLOG | ||
|
|
||
| sbnd_merge_sim_sources : { | ||
| module_type: "MergeSimSources" | ||
| FillMCParticles: false | ||
| FillSimPhotons: false | ||
| FillSimChannels: false | ||
| FillAuxDetSimChannels: false | ||
| FillSimEnergyDeposits: false | ||
| InputSourcesLabels: [ "larg4intime", "larg4outtime"] | ||
| TrackIDOffsets: [ 10000000, 20000000 ] | ||
| StoreReflected: true | ||
| EnergyDepositInstanceLabels: [ "priorSCE" ] | ||
| } | ||
|
|
||
| END_PROLOG |