icaruscode/fcl/detsim/detsim_1d_icarus.fcl at 3ac2499f389532dac9a9f1e0d1b67ca8d387542e · SBNSoftware/icaruscode · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
#
# Purpose: Runs detector digitisation simulation on G4 output.
#
#
# Optical detector noise is currently disabled.
#
#
# Input
# ------
#
# * `sim::SimChannel` collection (`largeant`): all TPC channels
# * `sim::SimPhoton` collection (`largeant`): all PMT channels
# * `sim::AuxDetSimChannel` collection (`largeant`): all CRT channels
#
#
# Output
# -------
#
# * `daq0:PHYSCRATEDATATPCEE`, `daq1:PHYSCRATEDATATPCEW`,
#   `daq2:PHYSCRATEDATATPCWE`, `daq3:PHYSCRATEDATATPCWW`:
#     TPC digitization waveforms (`raw::RawDigit` collections), one per TPC
# * `opdaq`: optical detector waveforms (`raw::OpDetWaveform` collection)
# * `crtdaq`: cosmic ray tagger hits (`icarus::crt::CRTData` collection)
#
#


# ------------------------------------------------------------------------------
#include "services_icarus_simulation.fcl"
#include "detsimmodules_ICARUS.fcl"
#include "opdetsim_pmt_icarus.fcl"
#include "crtsimmodules_icarus.fcl"
#include "rootoutput_icarus.fcl"
#include "channelmapping_icarus.fcl"


# ------------------------------------------------------------------------------
process_name: DetSim


# ------------------------------------------------------------------------------
services: {
    IICARUSChannelMap: @local::icarus_channelmappinggservice
    @table::icarus_detsim_services
}


# ------------------------------------------------------------------------------
physics:
{
  producers:
  {
    # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    crtdaq: @local::icarus_crtsim


    # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    opdaq:  @local::icarus_simpmt


    # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    daq0:   @local::icarus_simreadoutboard
    daq1:   @local::icarus_simreadoutboard
    daq2:   @local::icarus_simreadoutboard
    daq3:   @local::icarus_simreadoutboard

    # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    rns:    { module_type: "RandomNumberSaver" }

    # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

  } # physics

  simulate: [ rns, opdaq, daq0, daq1, daq2, daq3, crtdaq ]

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

} # physics

# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# point each of the SimWire instances to a different TPC set
#
physics.producers.daq0.OutputInstanceLabel: "PHYSCRATEDATATPCEE"
physics.producers.daq0.TPCVec:              [ [0, 0], [0, 1] ]
physics.producers.daq1.OutputInstanceLabel: "PHYSCRATEDATATPCEW"
physics.producers.daq1.TPCVec:              [ [0, 2], [0, 3] ]
physics.producers.daq2.OutputInstanceLabel: "PHYSCRATEDATATPCWE"
physics.producers.daq2.TPCVec:              [ [1, 0], [1, 1] ]
physics.producers.daq3.OutputInstanceLabel: "PHYSCRATEDATATPCWW"
physics.producers.daq3.TPCVec:              [ [1, 2], [1, 3] ]


# ------------------------------------------------------------------------------
outputs.rootoutput: @local::icarus_rootoutput

#legacy G4 configs
services.Geometry.GDML: "icarus_complete_20220518_overburden.gdml"
services.Geometry.ROOT: "icarus_complete_20220518_overburden.gdml"
physics.producers.crtdaq.G4ModuleLabel: "largeant"
physics.producers.opdaq.InputModule: "largeant"
physics.producers.opdaq.ApplyTimingDelays: false  # keep false for non-overlay MC