MUSUN.fcl
Go to the documentation of this file.
1 BEGIN_PROLOG
2 
3  # No experiment specific configurations because MUSUN is detector agnostic
4 
5 standard_MUSUN:
6 {
7  module_type: "MUSUN"
8 
9  PDG: -13 # PdG codes of particles to make
10  ChargeRatio: 1.38 # Ratio of particle/anti-particle
11  # For one particle type set equal to 0.
12 
13  # Directory where ROOT pdf file is. If not using one of the standard configurations
14  # you will !!NEED!! to change this to your own directory!!!
15  InputDir: "/dune/data/users/warburton/MUSUN/"
16  InputFile1: "muint-dune4850-mr-new.dat" # Table of muon intensities for theta and phi
17  InputFile2: "musp-dune4850-mr-new.dat" # Binary file of energies for theta and depths
18  InputFile3: "depth-dune4850-mr-new.dat" # Table of slant depths for theta and phi
19 
20  CavernAngle: 7 # Angle of the detector from the East to South.
21  RockDensity: 2.70 # Default rock density is 2.70 g cm-3. If this is
22  # changed then the three input files need to be
23  # remade. If there is a desire for this contact
24  # Vitaly Kudryavtsev at V.Kudryavtsev@shef.ac.uk
25 
26  Emin: 1. # Minimum Kinetic Energy GeV
27  Emax: 1000000. # Maximum Kinetic Energy GeV
28 
29  Thetamin: 0. # Minimum angle, must be more than or equal to 0.
30  Thetamax: 90. # Maximum angle, must be less than or equal to 90
31 
32  Phimin: 0. # Minimum angle, must be more than or equal to 0.
33  Phimax: 360. # Maximum angle, must be less than or equal to 360
34 
35  # You have to change all these dimensions if you change the size of the
36  # parallelepiped. The dimensions are given in cm.
37  # Ideally want the surface on which the muons are sampled to be extended by 5 m
38  # in all directions into rock, except top (y) where extension is 7m.
39  # Using /dune/app/users/talion/PrintDetBounds.C and specifiying which GDML file is
40  # desired, the co-ordinates of detector centre in the world volume and the extent
41  # of rock from detector are written out.
42  # Macro is also in dunetpc/dune/Geometry/gdml/PrintDetBounds.C
43  # Below co-ordinates are for dune10kt_v1_nowires.gdml, output is given too...
44  # Units: cm
45  # Coordinates: x makes a RHS y is vertical z is beam direction
46  # World Dimensions: x = 9554.03 y = 9417.7 z = 14321.9
47  # Det Dimensions: x = 1954.03 y = 1817.7 z = 6721.89
48  # Det Center: x = 0 y = 120.75 z = 2778
49  #
50  # Detector Enclosure Boundaries (everything outside is rock):
51  # x: -977.014 to 977.014
52  # y: -788.099 to 1029.6
53  # z: -582.946 to 6138.94
54  #
55  # Amount of rock on each side:
56  # -x: 3800
57  # +x: 3800
58  # -y: 3920.75
59  # +y: 3679.25
60  # -z: 6578
61  # +z: 1022
62 
63  igflag: 1 # If want muons sampled from parallelepiped, not a sphere.
64  Xmin: -1477. # Minumum X position cm
65  Xmax: 1477. # Maximum X position cm
66  Ymin: -1288. # Minumum Y position cm
67  Ymax: 1730. # Maximum Y position cm
68  Zmin: -1083. # Minumum Z position cm
69  Zmax: 6639. # Maximum Z position cm
70 
71 # Vitaly's Fortran code for LBNE uses these positions.
72 # Xmin: -1500. # Minumum X position cm
73 # Xmax: 1500. # Maximum X position cm
74 # Ymin: -1400. # Minumum Y position cm
75 # Ymax: 1600. # Maximum Y position cm
76 # Zmin: -3500. # Minumum Z position cm
77 # Zmax: 3500. # Maximum Z position cm
78 
79  T0: 0. # starting time
80  SigmaT: 0. # variation in the starting time
81  TDist: 0 # 0 - uniform, 1 - gaussian
82 }
83 
84 END_PROLOG