1 #include "services_dune.fcl"
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.
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.
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
33 module_type: "RadioGen"
34 Nuclide: [ "39Ar" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
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
44 T1: [ 2246000. ] # ending time in ns
48 module_type: "RadioGen"
49 Nuclide: [ "60Co" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
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
59 T1: [ 2246000. ] # ending time in ns
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.
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
79 module_type: "RadioGen"
80 Nuclide: [ "222Rn" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
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
90 T1: [ 2246000. ] # ending time in ns
94 module_type: "RadioGen"
95 Nuclide: [ "40K" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
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
105 T1: [ 2246000. ] # ending time in ns
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
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
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)
121 dune10kt_1x2x6_210po: @local::dune10kt_1x2x6_PDs
125 module_type: "RadioGen"
126 Nuclide: [ "39Ar" ] # list of nuclides to simulate, supported so far: 39Ar, 60Co, 85Kr, 40K, 232Th, 238U, 222Rn
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)
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
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
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)
1.8.11
