dune_radiological_model.fcl
Go to the documentation of this file.
1 #include "services_dune.fcl"
2 
3 BEGIN_PROLOG
4 
5 # parameters for a radiological generation generator. Currently, it generates radiologicals within rectangular prisms
6 # with sides parallel to the x, y, and z axes, and within a specified time window. The parameters Nuclide through
7 # T1 are arrays, so that many volumes can be specified to be radioactive with different nuclides.
8 
9 # remove excluded volume inside APA because we will have 39Ar there too!
10 # -JReichenbacher (01/11/2017)
11 # v2.01/JS: use only generator and not also arGen
12 # v2.1-02142017/JR: moved photon paddles by (0.5+0.476+0.001) cm in x
13 # v2.2-02222017/JR: cut out LAr volume inside APA around PDs (-0.477,...,0.477) cm in x
14 # v2.3-03212017/JR: include neutrons emitted from concrete of membrane cryostat structure (temporarily placed at aluminum field cage positions in sim.)
15 # 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.)
16 # 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
17 # v2.6-05302017/JS: Added neutronGen back in, using the Concrete_DUNE1 spectrum (by JR). New spectrum will be in dune_pardata.
18 # v2.7-06062017/JS: Added material specific generation for LAr contaminants.
19 # 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.
20 
21 #The Generator names provided that should be called from other files are:
22 #dune10kt_1x2x6_39ar #Ar39 in LAr
23 #dune10kt_1x2x6_APA #Co60 in APA frame
24 #dune10kt_1x2x6_neutron #Originally use for K40 in Field Cage. Now repurposed for Neutrons from concrete.
25 #dune10kt_1x2x6_CPA #K40 from CPA
26 #dune10kt_1x2x6_85kr #Kr in LAr
27 #dune10kt_1x2x6_222rn #Rn in LAr
28 #dune10kt_1x2x6_210po #po for 1 cm in front of PDs (approximates effect of Rn daughters on PDs
29 #
30 
31 dune10kt_1x2x6_39ar:
32 {
33  module_type: "RadioGen"
34  Nuclide: [ "39Ar" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
35  Material: ["LAr"]
36  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
37  X0: [ -475. ] # in cm in world coordinates, bottom corner of box
38  X1: [ 475. ] # in cm in world coordinates, top corner of box
39  Y0: [ -750. ] # in cm in world coordinates, bottom corner of box
40  Y1: [ 800. ] # in cm in world coordinates, top corner of box
41  Z0: [ -55. ] # in cm in world coordinates, bottom corner of box
42  Z1: [ 1500. ] # in cm in world coordinates, top corner of box
43  T0: [ -2246000 ]
44  T1: [ 2246000. ] # ending time in ns
45 }
46 
47 dune10kt_1x2x6_APA:{
48  module_type: "RadioGen"
49  Nuclide: [ "60Co" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
50  Material: [".*"]
51  BqPercc: [ 0.000082 ] # activity -- Becquerels per cc. 0.000082 assumes 45.5 mBq/kg (MPIK 2008 high measurement for 60Co) and virtual steel density inside APA frame of 0.0018 kg/cc = 260kg/(6.3m x 2.3m x 1cm) -JReichenbacher (01/11/2017)
52  X0: [ -0.5 ] # in cm in world coordinates, bottom corner of box
53  X1: [ 0.5 ] # in cm in world coordinates, top corner of box
54  Y0: [ -600. ] # in cm in world coordinates, bottom corner of box
55  Y1: [ 600. ] # in cm in world coordinates, top corner of box
56  Z0: [ 0 ] # in cm in world coordinates, bottom corner of box
57  Z1: [ 1395 ] # in cm in world coordinates, top corner of box
58  T0: [ -2246000 ]
59  T1: [ 2246000. ] # ending time in ns
60 }
61 
62 dune10kt_1x2x6_neutron: #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.
63 {
64  module_type: "RadioGen"
65  Nuclide: ["Concrete_DUNE10MeV", "Concrete_DUNE10MeV", "Concrete_DUNE10MeV", "Concrete_DUNE10MeV", "Concrete_DUNE10MeV" ]
66  Material: [".*",".*",".*",".*", ".*"]
67  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)
68  X0: [ 0.0, 0.0, 0.0, 0.0 , -0.5 ] # in cm in world coordinates, bottom corner of box
69  X1: [ 350., 350., 350., 350., 0.5 ] # in cm in world coordinates, top corner of box
70  Y0: [ -600., -600, -601., 600., -600 ] # in cm in world coordinates, bottom corner of box
71  Y1: [ 600. , 600, -600., 601., 600 ] # in cm in world coordinates, top corner of box
72  Z0: [ -0.5, 1394.5, 0, 0, 0 ] # in cm in world coordinates, bottom corner of box
73  Z1: [ 0.5, 1395.5, 1395, 1395, 1395] # in cm in world coordinates, top corner of box
74  T0: [ -2246000.,-2246000.,-2246000.,-2246000.,-2246000 ] # ending time in ns
75  T1: [ 2246000., 2246000., 2246000., 2246000., 2246000 ] # ending time in ns
76 }
77 
78 dune10kt_1x2x6_PDs:{
79  module_type: "RadioGen"
80  Nuclide: [ "222Rn" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
81  Material: [".*"]
82  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)
83  X0: [ 0.477] # in cm in world coordinates, bottom corner of box (moved it by 0.476 cm)
84  X1: [ 1.477] # in cm in world coordinates, top corner of box (moved it by 0.476 cm)
85  Y0: [ -600. ] # in cm in world coordinates, bottom corner of box
86  Y1: [ 600. ] # in cm in world coordinates, top corner of box
87  Z0: [ 0 ] # in cm in world coordinates, bottom corner of box
88  Z1: [ 1395 ] # in cm in world coordinates, top corner of box
89  T0: [ -2246000 ]
90  T1: [ 2246000. ] # ending time in ns
91 }
92 
93 dune10kt_1x2x6_CPA:{
94  module_type: "RadioGen"
95  Nuclide: [ "40K" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
96  Material: [".*"]
97  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)
98  X0: [ 349.5 ] # in cm in world coordinates, bottom corner of box
99  X1: [ 350.5 ] # in cm in world coordinates, top corner of box
100  Y0: [ -600. ] # in cm in world coordinates, bottom corner of box
101  Y1: [ 600. ] # in cm in world coordinates, top corner of box
102  Z0: [ 0 ] # in cm in world coordinates, bottom corner of box
103  Z1: [ 1395 ] # in cm in world coordinates, top corner of box
104  T0: [ -2246000 ]
105  T1: [ 2246000. ] # ending time in ns
106 }
107 
108 dune10kt_1x2x6_42ar: @local::dune10kt_1x2x6_39ar
109 dune10kt_1x2x6_42ar.Nuclide: [ "42Ar" ]
110 dune10kt_1x2x6_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
111 
112 dune10kt_1x2x6_85kr: @local::dune10kt_1x2x6_39ar
113 dune10kt_1x2x6_85kr.Nuclide: [ "85Kr" ]
114 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
115 
116 dune10kt_1x2x6_222rn: @local::dune10kt_1x2x6_39ar
117 dune10kt_1x2x6_222rn.Nuclide: [ "222Rn" ]
118 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)
119 
120 
121 dune10kt_1x2x6_210po: @local::dune10kt_1x2x6_PDs
122 
123 protodunesp_39ar:
124 {
125  module_type: "RadioGen"
126  Nuclide: [ "39Ar" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
127  Material: ["LAr"]
128  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
129  X0: [ -360. ] # in cm in world coordinates, bottom corner of box
130  X1: [ 360. ] # in cm in world coordinates, top corner of box
131  Y0: [ 0. ] # in cm in world coordinates, bottom corner of box
132  Y1: [ 600. ] # in cm in world coordinates, top corner of box
133  Z0: [ 0. ] # in cm in world coordinates, bottom corner of box
134  Z1: [ 700. ] # in cm in world coordinates, top corner of box
135  T0: [ -2500000. ] # max drift time (2.25 ms) + trigger offset (0.25 ms)
136  T1: [ 3000000. ] # ending time in ns (readout window size)
137 }
138 
139 protodunesp_42ar: @local::protodunesp_39ar
140 protodunesp_42ar.Nuclide: [ "42Ar" ]
141 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
142 
143 protodunesp_85kr: @local::protodunesp_39ar
144 protodunesp_85kr.Nuclide: [ "85Kr" ]
145 protodunesp_85kr.BqPercc: [ 0.00016 ] #Rate tied to LAr using rates established by the WARP collaboration. https://arxiv.org/pdf/astro-ph/0603131.pdf
146 
147 protodunesp_222rn: @local::protodunesp_39ar
148 protodunesp_222rn.Nuclide: [ "222Rn" ]
149 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)
150 
151 END_PROLOG