prodsingle_libraryana.fcl
Go to the documentation of this file.
1 #include "services_microboone.fcl"
2 #include "backtrackerservice.fcl"
3 #include "particleinventoryservice.fcl"
4 #include "caldata.fcl"
5 #include "hitfindermodules.fcl"
6 #include "singles.fcl"
7 #include "largeantmodules.fcl"
8 #include "detsimmodules.fcl"
9 #include "mccheatermodules.fcl"
10 #include "opticaldetectormodules.fcl"
11 #include "photpropmodules.fcl"
12 
13 
14 process_name: SinglesGen
15 
16 services:
17 {
18  # Load the service that manages root files for histograms.
19  TFileService: { fileName: "fastOptical_reconstruct_pmt15_3events_qe100_hist_uboone.root" }
20  TimeTracker: {}
21  RandomNumberGenerator: {} #ART native random number generator
22  @table::microboone_simulation_services
23  BackTrackerService: @local::microboone_backtrackerservice # defined in backtrackerservice.fcl
24  ParticleInventoryService: @local::standard_particleinventoryservice
25 }
26 services.PhotonVisibilityService: @local::standard_photonvisibilityservice_new
27 
28 #Start each new event with an empty event.
29 source:
30 {
31  module_type: EmptyEvent
32  maxEvents: 3 # Number of events to create
33  firstRun: 1 # Run number to use for this file
34  firstEvent: 1 # number of first event in the file
35 }
36 
37 # Define and configure some modules to do work on each event.
38 # First modules are defined; they are scheduled later.
39 # Modules are grouped by type.
40 physics:
41 {
42 
43  producers:
44  {
45  generator: @local::microboone_singlep
46  largeant: @local::microboone_largeant
47  daq: @local::microboone_simwire
48  backtrack: @local::standard_backtrackerloader
49  ffthit: @local::microboone_hitfinder
50  caldata: @local::microboone_calwire
51  opdigi: @local::microboone_opdigi
52  ophit: @local::microboone_ophitfinder
53  #ophitfinder: @local::standard_oplowintensityhitfinder
54  opflashfinder: @local::microboone_opflashfinder
55  opflashnew: @local::microboone_opflashnew
56  }
57 
58  analyzers:
59  {
60  libana: @local::microboone_photonlibraryanalyzer
61  largana: @local::microboone_largeantana
62  pmtresponse: @local::microboone_simphotoncounter
63  opdigiana: @local::microboone_opdigiana
64  ophitana: @local::microboone_ophitana
65  opflashana: @local:: microboone_opflashana
66  }
67 
68  #define the producer and filter modules for this path, order matters,
69  simulate: [ generator ]
70  analyzeIt: [ libana ]
71 
72  #define the output stream, there could be more than one if using filters
73  stream1: [ out1 ]
74 
75  #trigger_paths is a keyword and contains the paths that modify the art::event,
76  #ie filters and producers
77  trigger_paths: [simulate]
78 
79  #end_paths is a keyword and contains the paths that do not modify the art::Event,
80  #ie analyzers and output streams. these all run simultaneously
81  end_paths: [analyzeIt, stream1]
82 }
83 
84 #block to define where the output goes. if you defined a filter in the physics
85 #block and put it in the trigger_paths then you need to put a SelectEvents: {SelectEvents: [XXX]}
86 #entry in the output stream you want those to go to, where XXX is the label of the filter module(s)
87 outputs:
88 {
89  out1:
90  {
91  module_type: RootOutput
92  fileName: "fastOptical_reconstruct_pmt15_3events_qe100_gen_uboone.root" #default file name, can override from command line with -o or --output
93  }
94 }
95 
96 #manually setting parameters
97 services.DetectorPropertiesService.ScintYield: 24000
98 #(default in Utilities/larproperties.fcl number of photons per MeV of deposited energy)
99 
100 services.DetectorPropertiesService.ScintYieldRatio: 0.3
101 #(default in Utilities/larproperties.fcl)
102 
103 services.DetectorPropertiesService.ScintByParticleType: true
104 #(default in Utilities/larproperties.fcl, only used with fast optical simulation)
105 
106 services.DetectorPropertiesService.EnableCerenkovLight: false
107 #(default is true in Utilities/larproperties.fcl doesn't work with fast optical simulation
108 
109 services.DetectorPropertiesService.RayleighEnergies: [ 4.95, 5.69, 6.51, 7.53, 8.50, 9.67, 9.95, 10.57 ]
110 #(default in Utilities/larproperties.fcl, energy in eV)
111 
112 services.DetectorPropertiesService.RayleighSpectrum: [ 3790, 1850, 930, 410, 192, 90, 38, 10 ]
113 #(default in Utilities/larproperties.fcl, length in cm)
114 
115 services.DetectorPropertiesService.AbsLengthEnergies: [ 4, 5, 6, 7, 8, 9, 10, 11 ]
116 #(default in Utilities/larproperties.fcl, energy in eV)
117 
118 services.DetectorPropertiesService.AbsLengthSpectrum: [ 2000., 2000., 2000., 2000., 2000., 2000., 2000., 2000. ]
119 #(default in Utilities/larproperties.fcl, length in cm)
120 
121 physics.analyzers.pmtresponse.QuantumEfficiency: 1.00 #0.03
122 #(default is 0.015 in OpticalDetector/opticaldetectormodules.fcl)
123 
124 physics.analyzers.pmtresponse.WavelengthCutLow: 0
125 #(default is 0 in OpticalDetector/opticaldetectormodules.fcl)
126 
127 physics.analyzers.pmtresponse.WavelengthCutHigh: 1000
128 #(default is 1000 in OpticalDetector/opticaldetectormodules.fcl)
129 
130 # enable optical physics in LArG4
131 services.LArG4Parameters.EnabledPhysics: [ "Em",
132  "FastOptical",
133  "SynchrotronAndGN",
134  "Ion",
135  "Hadron",
136  "Decay",
137  "HadronElastic",
138  "Stopping",
139  "NeutronTrackingCut" ]
140 
141 # enable this custom physics list
142 services.LArG4Parameters.UseCustomPhysics: true
143 
144 # disable cerenkov light
145 services.DetectorPropertiesService.EnableCerenkovLight: false
146 
147 
148 #17.5 < x < 238.8 cm (110.65 center)
149 #-99 < y < 99.1 cm
150 #17.5 < z < 1019.3 cm
151 
152 #in front on PMT15
153 physics.producers.generator.T0: [500]
154 physics.producers.generator.X0: [20.]
155 physics.producers.generator.Y0: [60.]
156 physics.producers.generator.Z0: [532.]
157 
158 physics.producers.generator.SigmaX: [0.]
159 physics.producers.generator.SigmaY: [0.]
160 physics.producers.generator.SigmaZ: [0.]
161 physics.producers.generator.PosDist: 0 #0 uniform, 1 gaussian
162 
163 physics.producers.generator.PDG: [2212]
164 physics.producers.generator.P0: [0.040]
165 physics.producers.generator.SigmaP: [0.]
166 physics.producers.generator.PDist: 1 #0 uniform, 1 gaussian
167 physics.producers.generator.SigmaThetaXZ: [ 0 ]
168 physics.producers.generator.SigmaThetaYZ: [ 0 ]
169 
170 services.PhotonVisibilityService.LibraryFile: "PhotonPropagation/LibraryData/lib8984855.root"