All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
simulationservices_icarus.fcl
Go to the documentation of this file.
1 #include "simulationservices.fcl"
2 #include "backtrackerservice.fcl"
4 
5 BEGIN_PROLOG
6 
7 #
8 # Configuration for LArG4Parameters service
9 #
10 # Use as:
11 #
12 # service.LArG4Parmeters: @local::icarus_largeantparameters
13 #
14 #
16 
17  OpticalSimVerbosity: 0 #verbosity of optical simulation, soon to be depricated
19  StoreTrajectories: true
20  VisualizationEnergyCut: 10.e-3 #depricated, in GeV
21  VisualizeNeutrals: false #depricated
22  ModifyProtonCut: false #Whether to modify the default proton cut
23  NewProtonCut: 0.0 #new ProtonCut value, ModifyProtonCut must be set to set new value
24  KeepEMShowerDaughters: false #Don't save secondary, tertiary, etc particles in EM showers
25  LongitudinalDiffusion: 4.0e-9 #in cm^2/ns
26  TransverseDiffusion: 8.8e-9 #in cm^2/ns
27  ElectronClusterSize: 20.0 #number of ionization electrons to drift in a unit
28  MinNumberOfElCluster: 0 #minimum number of electron clusters
29  UseCustomPhysics: true # use the `EnabledPhysics` list of processes:
30  EnabledPhysics: [ "Em", "SynchrotronAndGN", "Ion", "Hadron",
31  "Decay", "HadronElastic", "Stopping"
32  , "FastOptical" # either: use photon library
33  # , OpFastScintillation # or: let Geant4 do the work (for photon visibility library generation)
34  # , "NeutronTrackingCut" # if NeutronTrackingCut is enabled, neutrons will be cut
35  ]
36  CosmogenicK0Bias: 0 # 0 is off. N is the number of secondaries to produce.
37  CosmogenicXSMNBiasOn: 0 # 0 is off. 1 works. 2 still in development.
38  CosmogenicXSMNBiasFactor: 1 # Not more than 5-ish cuz of numerical instabilities.
39  DisableWireplanes: false #if set true, charge drift simulation does not run - used for optical sim jobs OR just when you don't wanna drift the e's.
40  SkipWireSignalInTPCs: [] # put here TPC id's which should not receive ionization electrons - used to simulate TPC geom volumes which are actually dead LAr volumes in protoDUNE
41  FillSimEnergyDeposits: true
42  UseModBoxRecomb: true # use Modified Box recombination instead of Birks
43  UseModLarqlRecomb: true # use LArQL recombination corrections (dependence on EF)
44 
45  #* Recombination factor coefficients come from Nucl.Instrum.Meth.A523:275-286,2004
46  #* * @f$ dE/dx @f$ is given by the voxel energy deposition, but have to convert it to MeV/cm from GeV/voxel width
47  #* * electric field: @f$ E @f$ in kV/cm
48  #* * @f$ R = A/(1 + (dE/dx)*k/E) @f$
49  #* * @f$ A = 0.800 \pm 0.003 @f$
50  #* * @f$ k = 0.0486 @f$ needs to be scaled with Electric field
51  RecombA: 0.800 #< _A_ constant.
52  Recombk: 0.0486 #< _k_ constant, in g/(MeV cm&sup2;)*kV/cm.
53 
54  #* Recombination factor coefficients come from Nucl.Instrum.Meth.A523:275-286,2004
55  #* * @f$ dE/dx @f$ is given by the voxel energy deposition, but have to convert it to MeV/cm from GeV/voxel width
56  #* * electric field: @f$ E @f$ in kV/cm
57  #* * `kModBoxB` needs to be scaled with the electric field.
58  ModBoxA: 0.930
59  ModBoxB: 0.212
60 
61  #* Recombination factor coefficients for LArQL can be found in https://cdcvs.fnal.gov/redmine/projects/larsoft/wiki/LArQL_algorithm
62  #* * @f$ dE/dx @f$ is given by the energy deposition in MeV/cm
63  #* * electric field: @f$ E @f$ in kV/cm
64  LarqlChi0A: 0.00338427
65  LarqlChi0B: -6.57037
66  LarqlChi0C: 1.88418
67  LarqlChi0D: 0.000129379
68  LarqlAlpha: 0.0372
69  LarqlBeta: 0.0124
70 
71  #* ion+excitation work function in eV
72  #* https://doi.org/10.1016/0168-9002(90)90011-T
73  Wph: 19.5
74 
75  # adopt the "correlated" model for ionization and scintillation (SBN DocDB 17964)
76  IonAndScintCalculator: "Correlated"
77  # this enables some improvements on top of the correlated model
78  # (LArQL, using the parameters above; see SBN DocDB 21270)
79  # UseModLarqlRecomb: true
80 
81  # The following parameters specify details of wireplanes or similar
82  # areas with optically parameterized transmissions (Ben J 2013)
83  # volume names to be associated with an optical wireplane model
84  OpticalParamVolumes: ["volTPCPlaneVert_PV"]
85 
86  # specification of which model to use for each volume
87  OpticalParamModels: ["OverlaidWireplanes"]
88 
89  # orientation of each wireplane set
90  # 0 = Xdrift, 1 = Ydrift, 2 = Zdrift
91  OpticalParamOrientations: [0]
92 
93  # This a set of floats which is specific to the particular model used.
94  # For overlaid wireplanes, should be a vector of vectors of
95  # [plane angle, pitch/mm, wire diameter/mm]
96  #
97  # This format is chosen to allow for future extensions to the model
98  # for, eg, DUNE wireplane development.
99  #
100  OpticalParamParameters: [ [ [ 60, 3, 0.15],
101  [-60, 3, 0.15],
102  [ 0, 3, 0.15] ] ]
103  UseLitePhotons: false
104 
105 
106 }
107 
108 
109 #
110 # ICARUS LArVoxelCalculator service configuration.
111 #
112 # Use as:
113 #
114 # services.LArVoxelCalculator: @local::icarus_larvoxelcalculator
115 #
116 #
117 icarus_larvoxelcalculator: {
118  VoxelSizeX: 0.03 #in cm
119  VoxelSizeY: 0.03 #in cm
120  VoxelSizeZ: 0.03 #in cm
121  VoxelSizeT: 5000.0 #in ns
122  VoxelOffsetX: 0.0 #in cm
123  VoxelOffsetY: 0.0 #in cm
124  VoxelOffsetZ: 0.0 #in cm
125  VoxelOffsetT: -2500.0 #in ns
126  VoxelEnergyCut: 1.e-6 #in GeV
127 }
128 
129 
130 
131 #
132 # Backtracking services
133 #
134 # All serviced needed for service-driven backtracking
135 # Use as:
136 #
137 # services: {
138 # ...
139 # @table::icarus_backtracking_services
140 # ...
141 # }
142 #
143 icarus_backtracking_services: {
144  BackTrackerService: @local::standard_backtrackerservice # from `backtrackerservice.fcl` (`larsim`)
145  PhotonBackTrackerService: @local::standard_photonbacktrackerservice # from `photonbacktrackerservice.fcl` (`larsim`)
146  ParticleInventoryService: @local::standard_particleinventoryservice
147 } # icarus_backtracking_services
148 
149 
150 END_PROLOG
151 
Planes which measure V.
Definition: geo_types.h:130
esac done echo Signal files are
Definition: TrainMVA.sh:25
#define the
BEGIN_PROLOG soon to be depricated in GeV VisualizeNeutrals
process_name E
BEGIN_PROLOG g
BEGIN_PROLOG TPC
BEGIN_PROLOG or simple_flux see https
Definition: genie.fcl:8
constexpr double kModBoxB
Modified Box Beta in g/(MeV cm)*kV/cm.
A value measured in the specified unit.
Definition: quantities.h:566
then echo echo For and will not be changed by echo further linking echo echo B echo The symbol is in the uninitialized data multiple common symbols may appear with the echo same name If the symbol is defined the common echo symbols are treated as undefined references For more echo details on common see the discussion of warn common echo in *Note Linker see the discussion of warn common echo in *Note Linker such as a global int variable echo as opposed to a large global array echo echo I echo The symbol is an indirect reference to another symbol This echo is a GNU extension to the a out object file format which is echo rarely used echo echo N echo The symbol is a debugging symbol echo echo R echo The symbol is in a read only data section echo echo S echo The symbol is in an uninitialized data section for small echo objects echo echo T echo The symbol is in the the normal defined echo symbol is used with no error When a weak undefined symbol echo is linked and the symbol is not the value of the echo weak symbol becomes zero with no error echo echo W echo The symbol is a weak symbol that has not been specifically echo tagged as a weak object symbol When a weak defined symbol echo is linked with a normal defined the normal defined echo symbol is used with no error When a weak undefined symbol echo is linked and the symbol is not the value of the echo weak symbol becomes zero with no error echo echo echo The symbol is a stabs symbol in an a out object file In echo this the next values printed are the stabs other field
def recombination
Definition: dedx.py:128
return match has_match and(match.match_pdg==11 or match.match_pdg==-11)
BEGIN_PROLOG soon to be depricated in GeV ModifyProtonCut must be set to set new value KeepEMShowerDaughters
physics simulate
if &&[-z"$BASH_VERSION"] then echo Attempting to switch to bash bash shellSwitch exit fi &&["$1"= 'shellSwitch'] shift declare a IncludeDirectives for Dir in
then echo File list $list not found else cat $list while read file do echo $file sed s
Definition: file_to_url.sh:60
BEGIN_PROLOG soon to be depricated ParticleKineticEnergyCut
do i e
temporary value
process_name larg4outtime ionization
BEGIN_PROLOG icarus_largeantparameters
TimeTrackTreeStorage::TriggerInputSpec_t convert(TimeTrackTreeStorage::Config::TriggerSpecConfig const &config)
services LArG4Parameters UseLitePhotons