mcc12_gen_protoDune_beam_cosmics.fcl

Go to the documentation of this file.

#include "services_dune.fcl"
#include "singles_dune.fcl"
#include "corsika_protodune.fcl"
#include "protodunebeam.fcl"
#include "dune_radiological_model.fcl"

process_name: SinglesGen

services:
{
  # Load the service that manages root files for histograms.
  TFileService: { fileName: "gensingle_protoDUNE_hist.root" }
  TimeTracker:       @local::dune_time_tracker
  MemoryTracker:     @local::dune_memory_tracker
  RandomNumberGenerator: {} #ART native random number generator
  FileCatalogMetadata:  @local::art_file_catalog_mc
  @table::protodune_simulation_services
}

#Start each new event with an empty event.
source:
{
  module_type: EmptyEvent
  timestampPlugin: { plugin_type: "GeneratedEventTimestamp" }
  maxEvents:   1000000
  firstRun:    1           # Run number to use for this file
  firstEvent:  1           # number of first event in the file
}

# 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:
 {
   generator: @local::protodune_beam
   cosmicgenerator: @local::protodune_corsika_cmc  
   ar39:  @local::protodunesp_39ar
   ar42:  @local::protodunesp_42ar
   kr85:  @local::protodunesp_85kr
   rn222: @local::protodunesp_222rn
 }


 #define the producer and filter modules for this path, order matters, 
 #filters reject all following items.  see lines starting physics.producers below
 simulate: [ generator, cosmicgenerator, ar39, ar42, kr85, rn222 ]
 
 #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:    "gensingle_protoDUNE.root" #default file name, can override from command line with -o or --output
   dataTier: "generated"
   compressionLevel: 1
 }
}

#Set generator parameters
#Corresponds to beam window at center of left TPC

# Specific to the GEANT H4 Beam simulation by Nikos and paths at CERN
physics.producers.generator.GoodParticleTreeName: "NTuples/GoodParticle"
#physics.producers.generator.AllParticlesTreeName: "Detector/NP04front"
physics.producers.generator.AllParticlesTreeName: "Detector/NP04front" # Background particles should fire from here
physics.producers.generator.StartEvent: 0
# Set the beam position and direction to fire the particles in the desired way
#physics.producers.generator.BeamX: 7.966 # In cm, taken from protoDUNE_gensingle.fcl
#physics.producers.generator.BeamX: 5.55 # In cm, based on studies of the beam pipe position
#physics.producers.generator.BeamY: 460.84 
#physics.producers.generator.BeamZ: -191.60
physics.producers.generator.BeamX: 8.905247 # In cm
physics.producers.generator.BeamY: 461.749224 
physics.producers.generator.BeamZ: -196.11
physics.producers.generator.BeamThetaShift: 0.0114276
physics.producers.generator.BeamPhiShift: -0.0761654
# Specific rotations for the good particles
physics.producers.generator.RotateXZ: 0.0
physics.producers.generator.RotateYZ: 0.0 
# Rotations for the final beam monitors and background particles
#physics.producers.generator.RotateMonitorsXZ: -10.3 # [Degrees]
#physics.producers.generator.RotateMonitorsYZ: 11.7  # [Degrees]
# Jake Calcutt: Replacing these with angles for opposite order rotations
physics.producers.generator.RotateMonitorsXZ: -10.51372 # [Degrees]
physics.producers.generator.RotateMonitorsYZ: 11.50886  # [Degrees]

# Jake Calcutt: Extra parameters for beam interface info
#physics.producers.generator.L1: 1.98
#physics.producers.generator.L2: 1.69472
#physics.producers.generator.L3: 2.11666
#physics.producers.generator.BeamBend: .12003
#
#physics.producers.generator.LMag: 2.1645
#physics.producers.generator.NominalP: 1.
#physics.producers.generator.B: 1.294817


# Overlay info
physics.producers.generator.ReadoutWindow: 4.0

# We will set intensity and file name in the specific .fcl files for each momentum bin.