standard_g4_icarus_21apr21_laronly.fcl

Go to the documentation of this file.

#include "services_icarus_simulation.fcl"
#include "largeantmodules_icarus.fcl"
#include "mcreco.fcl"

process_name: G4

services:
{
  @table::icarus_g4_services
}

#source is now a root file
source:
{
  module_type: RootInput
  saveMemoryObjectThreshold: 0
}

# Define and configure some modules to do work on each event.
# First modules are defined; they are scheduled later.
# Modules are grouped by type.
physics:
{

 producers:
 {
   largeant: @local::icarus_largeant
   ionization: {
                 module_type: "larsim/ElectronDrift/ShiftEdepSCE"
                 EDepTag:     "largeant:TPCActive"
                 MakeAnaTree: false
               }
   mcreco:   @local::standard_mcreco
   rns:        { module_type: "RandomNumberSaver" }

 }

 #define the producer and filter modules for this path, order matters,
 #filters reject all following items.  see lines starting physics.producers below
 simulate: [ rns, largeant, ionization, mcreco ]

 #define the output stream, there could be more than one if using filters
 stream1:  [ out1 ]

 #trigger_paths is a keyword and contains the paths that modify the art::event,
 #ie filters and producers
 trigger_paths: [simulate]

 #end_paths is a keyword and contains the paths that do not modify the art::Event,
 #ie analyzers and output streams.  these all run simultaneously
 end_paths:     [stream1]
}

#block to define where the output goes.  if you defined a filter in the physics
#block and put it in the trigger_paths then you need to put a SelectEvents: {SelectEvents: [XXX]}
#entry in the output stream you want those to go to, where XXX is the label of the filter module(s)
outputs:
{
 out1:
 {
   module_type: RootOutput
   fileName:    "%ifb_%tc-%p.root"
   dataTier:    "simulated"
   saveMemoryObjectThreshold: 0
   compressionLevel: 1
   fastCloning: false
 }
}

# drop the sim::SimPhotons
#outputs.out1.outputCommands: ["keep *", "drop sim::SimPhotonss_*_*_*", "drop sim::SimChannels_*_*_*"]

#physics.producers.largeant.KeepParticlesInVolumes: ["volCryostat"] #only keep MCParticles that cross the cryostat

# turn on the SimEnergyDeposit output
services.LArG4Parameters.FillSimEnergyDeposits:       true

physics.producers.largeant.KeepParticlesInVolumes: []
physics.producers.mcreco.MCShowerRecoAlg.MinShowerEnergy: 1 # MeV
physics.producers.mcreco.MCShowerRecoAlg.MinNumDaughters: 0

services.LArG4Parameters.EnabledPhysics: [
         "Em",
         "FastOptical",
         "SynchrotronAndGN",
         "Ion",
         "Hadron",
         "Decay",
         "HadronElastic",
         "Stopping"
      ]
services.LArG4Parameters.ParticleKineticEnergyCut: 1e-6 # 1 keV
physics.producers.largeant.SparsifyTrajectories: true

# Geometry with LAr only
services.Geometry.GDML: "icarus_complete_laronly_20210511.gdml"
services.Geometry.ROOT: "icarus_complete_laronly_20210511.gdml"

services.GeometryConfigurationWriter: {}

services.message.destinations :
{
  STDCOUT:
  {
     type:      "cout"      #tells the message service to output this destination to cout
     threshold: "INFO"      #tells the message service that this destination applies to WARNING and higher level messages
     categories:
     {
       default:
       {
         limit:       1000  #don't print anything at the infomsg level except the explicitly named categories
         reportEvery: 1
       }
     }
  }
}