dune_radiological_model_decay0.fcl
Go to the documentation of this file.
1 #include "services_dune.fcl"
2 
3 BEGIN_PROLOG
4 # parameters for a radiological generation generator. Currently, it generates radiologicals within rectangular prisms
5 # with sides parallel to the x, y, and z axes, and within a specified time window. The parameters Nuclide through
6 # T1 are arrays, so that many volumes can be specified to be radioactive with different nuclides.
7 
8 # remove excluded volume inside APA because we will have 39Ar there too!
9 # -JReichenbacher (01/11/2017)
10 # v2.01/JS: use only generator and not also arGen
11 # v2.1-02142017/JR: moved photon paddles by (0.5+0.476+0.001) cm in x
12 # v2.2-02222017/JR: cut out LAr volume inside APA around PDs (-0.477,...,0.477) cm in x
13 # v2.3-03212017/JR: include neutrons emitted from concrete of membrane cryostat structure (temporarily placed at aluminum field cage positions in sim.)
14 # v2.4-05122017/JS: Removed the LAr cut out. Added in the 210po generation near the PDs. Disabled FC (Until Neutrons can be pushed to LArData.)
15 # v2.5-05252017/JS: Added in 85Kr as a background in LAr using the WARP collaboration rates. https://arxiv.org/pdf/astro-ph/0603131.pdf
16 # v2.6-05302017/JS: Added neutronGen back in, using the Concrete_DUNE1 spectrum (by JR). New spectrum will be in dune_pardata.
17 # v2.7-06062017/JS: Added material specific generation for LAr contaminants.
18 # v2.8-06092017/JS: Reconfigured this to be an include for a typical job. This file will now host generators for the various radiologicals, but will not run a simulation itself.
19 
20 #The Generator names provided that should be called from other files are:
21 #dune10kt_1x2x6_39ar #Ar39 in LAr
22 #dune10kt_1x2x6_APA #Co60 in APA frame
23 #dune10kt_1x2x6_neutron #Originally use for K40 in Field Cage. Now repurposed for Neutrons from concrete.
24 #dune10kt_1x2x6_CPA #K40 from CPA
25 #dune10kt_1x2x6_85kr #Kr in LAr
26 #dune10kt_1x2x6_222rn #Rn in LAr
27 #dune10kt_1x2x6_210po #po for 1 cm in front of PDs (approximates effect of Rn daughters on PDs
28 #
29 
30 dune10kt_1x2x6_weird_beta_example_in_LAr:
31 {
32  # module_type: "Decay0Gen"
33  module_type: "SpectrumVolumeGen"
34  isotope: "11"
35  material: "LAr"
36  spectrum_p_min: 0
37  spectrum_p_max: 0.01
38  spectrum: [0,0,1,1,1,1,4,4,1,1,0,0]
39  ## OR
40  # spectrum_p_min: 0
41  # spectrum_p_max: 0.01
42  # nbins:100
43  # function: "[0]*x+[1]*x*x+[2]"
44  # parameters: [100,0,0]
45  ## OR
46  # spectrum: [0,1,4,4,1.5]
47  # bins: [0.,0.02,0.04,0.08,0.15,0.3]
48  ## OR
49  # bins: [0.,0.02,0.04,0.08,0.15,0.3]
50  # function: "[0]*x+[1]*x*x+[2]"
51  # parameters: [100,0,0]
52  BqPercc: 0.00141 # activity -- Becquerels per cc. 0.00141 assumes 1.01 Bq/kg (typical for 39Ar) and a density of 1.396 g/cc for LAr
53  volume_rand: "volCryostat"
54 }
55 
56 dune10kt_1x2x6_85Kr_in_LAr:{
57  module_type: "Decay0Gen"
58  isotope: "Kr85"
59  volume_rand: "volCryostat"
60  material: "LAr"
61  BqPercc: 0.00016 #Rate tied to LAr using rates established by the WARP collaboration. https://arxiv.org/pdf/astro-ph/0603131.pdf
62 }
63 
64 dune10kt_1x2x6_39Ar_in_LAr:{
65  module_type: "Decay0Gen"
66  isotope: "Ar39"
67  volume_rand: "volCryostat"
68  material: "LAr"
69  BqPercc: 0.00141 # activity -- Becquerels per cc. 0.00141 assumes 1.01 Bq/kg (typical for 39Ar) and a density of 1.396 g/cc for LAr
70 }
71 
72 
73 dune10kt_1x2x6_42Ar_in_LAr:{
74  module_type: "Decay0Gen"
75  decay_chain:{
76  isotope_0: "Ar42"
77  isotope_1: "K42"
78  }
79  volume_rand: "volCryostat"
80  material: "LAr"
81  BqPercc: 0.0000001283768
82 }
83 
84 dune10kt_1x2x6_60Co_in_APA:{
85  module_type: "Decay0Gen"
86  material: "STEEL_STAINLESS_Fe7Cr2Ni"
87  volume_rand: "volCryostat"
88  volume_gen:".*APA.*"
89  isotope: "Co60"
90  BqPercc: 0.000082 # activity -- Becquerels per cc. 0.000082 assumes 45.5 mBq/kg (MPIK 2008 high measuremen)
91 }
92 
93 dune10kt_1x2x6_40K_in_CPA:{
94  module_type: "Decay0Gen"
95  isotope: "K40"
96  volume_rand: "volCryostat"
97  volume_gen:".*Cathode.*"
98  material: "STEEL_STAINLESS_Fe7Cr2Ni"
99  BqPercc: 0.0027195 # activity -- Becquerels per cc. 0.0027195 assumes 4.9 Bq/kg (measured for 40K) and a density of 1.85 g/cc for 3 mm thick G10 -JReichenbacher (01/11/2017)
100 }
101 
102 dune10kt_1x2x6_238U_fullchain_in_LAr:{
103  module_type:"Decay0Gen"
104  BqPercc: 0.000001395
105  material: "LAr"
106  volume_rand: "volCryostat"
107  volume_gen: ".*"
108  decay_chain:{ # ALL THESE ARE ASSUMED TO BE @ EQUILIBRIUM!!
109  isotope_0:"U238"
110  isotope_1:"Th234"
111  isotope_2:"Pa234m"
112  isotope_3:"U234"
113  isotope_4:"Th230"
114  isotope_5:"Ra226"
115  isotope_6:"Rn222"
116  isotope_7:"Po218"
117  isotope_8:"Pb214"
118  isotope_9:"Bi214"
119  isotope_10:"Pb210"
120  isotope_11:"Bi210"
121  isotope_12:"Po210"
122  }
123 }
124 
125 dune10kt_1x2x6_neutron_from_rock: #Center region is added to simulate neutrons comming in from the rock at the outermost APAs (on the edge of the cryostat) per J. Reichenbacher's request. This should more accurately represent the worst case scenario in the FD.
126 {
127  module_type: "RadioGen"
128  Nuclide: ["Concrete_DUNE10MeV", "Concrete_DUNE10MeV", "Concrete_DUNE10MeV", "Concrete_DUNE10MeV", "Concrete_DUNE10MeV" ]
129  Material: [".*",".*",".*",".*", ".*"]
130  BqPercc: [ 0.0000076, 0.0000076, 0.0000076, 0.0000076, 0.0000076 ] # activity -- Becquerels per cc. 0.0000076 assumes 10 neutrons/y/g per 10 ppm U-238 in concrete and a density of 2.40 g/cc for 10 cm mean depth in concrete -JReichenbacher (03/21/2017)
131  T0: [ -2246000.,-2246000.,-2246000.,-2246000.,-2246000 ] # ending time in ns
132  T1: [ 2246000., 2246000., 2246000., 2246000., 2246000 ] # ending time in ns
133 
134  X0: [ 0. , 0. , 0., 0., -0.5 ] # in cm in world coordinates, bottom corner of box
135  X1: [ 350. , 350. , 350., 350., 0.5 ] # in cm in world coordinates, top corner of box
136  Y0: [ -600. , -600. , -601., 600., -600. ] # in cm in world coordinates, bottom corner of box
137  Y1: [ 600. , 600. , -600., 601., 600. ] # in cm in world coordinates, top corner of box
138  Z0: [ -0.5, 1394.5, 0., 0., 0. ] # in cm in world coordinates, bottom corner of box
139  Z1: [ 0.5, 1395.5, 1395., 1395., 1395. ] # in cm in world coordinates, top corner of box
140 }
141 
142 # dune10kt_1x2x6_PDs:{
143  # module_type: "Decay0Gen"
144  # Nuclide: [ "222Rn" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
145  # Material: [".*"]
146  # BqPercc: [ 0.000005 ] # activity -- Becquerels per cc. 0.000005 assumes 2x 12.5% PD coverage in APA frame and a surface alpha activity requirement of 0.2 Bq/m^2 -JReichenbacher (02/14/2017)
147  # X0: [ 0.477] # in cm in world coordinates, bottom corner of box (moved it by 0.476 cm)
148  # X1: [ 1.477] # in cm in world coordinates, top corner of box (moved it by 0.476 cm)
149  # Y0: [ -600. ] # in cm in world coordinates, bottom corner of box
150  # Y1: [ 600. ] # in cm in world coordinates, top corner of box
151  # Z0: [ 0 ] # in cm in world coordinates, bottom corner of box
152  # Z1: [ 1395 ] # in cm in world coordinates, top corner of box
153  # T0: [ -2246000 ]
154  # T1: [ 2246000. ] # ending time in ns
155  # }
156 
157 # dune10kt_1x2x6_CPA:{
158  # module_type: "Decay0Gen"
159  # Nuclide: [ "40K" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
160  # Material: [".*"]
161  # BqPercc: [ 0.0027195 ] # activity -- Becquerels per cc. 0.0027195 assumes 4.9 Bq/kg (measured for 40K) and a density of 1.85 g/cc for 3 mm thick G10 -JReichenbacher (01/11/2017)
162  # X0: [ 349.5 ] # in cm in world coordinates, bottom corner of box
163  # X1: [ 350.5 ] # in cm in world coordinates, top corner of box
164  # Y0: [ -600. ] # in cm in world coordinates, bottom corner of box
165  # Y1: [ 600. ] # in cm in world coordinates, top corner of box
166  # Z0: [ 0 ] # in cm in world coordinates, bottom corner of box
167  # Z1: [ 1395 ] # in cm in world coordinates, top corner of box
168  # T0: [ -2246000 ]
169  # T1: [ 2246000. ] # ending time in ns
170  # }
171 
172 # dune10kt_1x2x6_42ar: @local::dune10kt_1x2x6_39ar
173 # dune10kt_1x2x6_42ar.Nuclide: [ "42Ar" ]
174 # dune10kt_1x2x6_42ar.BqPercc: [] #Rate from paper by Barabash in conference proceding http://iopscience.iop.org/article/10.1088/1742-6596/718/6/062004/pdf
175 
176 # dune10kt_1x2x6_85kr: @local::dune10kt_1x2x6_39ar
177 # dune10kt_1x2x6_85kr.Nuclide: [ "85Kr" ]
178 # dune10kt_1x2x6_85kr.BqPercc: [ 0.00016 ] #Rate tied to LAr using rates established by the WARP collaboration. https://arxiv.org/pdf/astro-ph/0603131.pdf
179 
180 # dune10kt_1x2x6_222rn: @local::dune10kt_1x2x6_39ar
181 # dune10kt_1x2x6_222rn.Nuclide: [ "222Rn" ]
182 # dune10kt_1x2x6_222rn.BqPercc: [ 0.000005584 ] # DUNE FD requirement of 1 mBq/kg (-JS) times 4 (to emulate the 4 subsequent alpha decays in secular equilibrium) -JReichenbacher (01/11/2017)
183 
184 
185 # dune10kt_1x2x6_210po: @local::dune10kt_1x2x6_PDs
186 
187 protodunesp_39ar:
188 {
189  module_type: "Decay0Gen"
190  Nuclide: [ "39Ar" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
191  Material: ["LAr"]
192  BqPercc: [ 0.00141 ] # activity -- Becquerels per cc. 0.00141 assumes 1.01 Bq/kg (typical for 39Ar) and a density of 1.396 g/cc for LAr
193  X0: [ -360. ] # in cm in world coordinates, bottom corner of box
194  X1: [ 360. ] # in cm in world coordinates, top corner of box
195  Y0: [ 0. ] # in cm in world coordinates, bottom corner of box
196  Y1: [ 600. ] # in cm in world coordinates, top corner of box
197  Z0: [ 0. ] # in cm in world coordinates, bottom corner of box
198  Z1: [ 700. ] # in cm in world coordinates, top corner of box
199  T0: [ -2500000. ] # max drift time (2.25 ms) + trigger offset (0.25 ms)
200  T1: [ 3000000. ] # ending time in ns (readout window size)
201 }
202 
203 protodunesp_42ar: @local::protodunesp_39ar
204 protodunesp_42ar.Nuclide: [ "42Ar" ]
205 protodunesp_42ar.BqPercc: [ 0.0000001283768 ] #Rate from paper by Barabash in conference proceding http://iopscience.iop.org/article/10.1088/1742-6596/718/6/062004/pdf
206 
207 protodunesp_85kr: @local::protodunesp_39ar
208 protodunesp_85kr.Nuclide: [ "85Kr" ]
209 protodunesp_85kr.BqPercc: [ 0.00016 ] #Rate tied to LAr using rates established by the WARP collaboration. https://arxiv.org/pdf/astro-ph/0603131.pdf
210 
211 protodunesp_222rn: @local::protodunesp_39ar
212 protodunesp_222rn.Nuclide: [ "222Rn" ]
213 protodunesp_222rn.BqPercc: [ 0.000005584 ] # DUNE FD requirement of 10 mBq/kg times 4 (to emulate the 4 subsequent alpha decays in secular equilibrium) -JReichenbacher (01/11/ 2017)
214 
215 END_PROLOG
Generated by   doxygen 1.8.11