Project

General

Profile

GNuMIFlux.xml

Marco Del Tutto, 12/17/2015 04:18 AM

 
1
<?xml version="1.0" encoding="ISO-8859-1"?>
2

    
3
<!-- Configurations for GNuMI flux  
4
     Each is grouped into a <param_set> that are parsed and applied to 
5
     GNuMIFlux.  Elements should appear in the following order:
6

7
     <verbose> : how noisy to be while parsing <param_set>
8
     <units>   : sets user coordinate units (e.g. "cm", "meter", etc)
9
                 set this first and only once
10
     <beamdir> : defines the relative orientation of the beam in user coord
11
                 three types are currently supported:
12
           "series"    : a sequence of rotations around one axis at a time
13
                         each <rotation> needs value and attributes:
14
                             axis ('x','y','z') units ('deg' 'rad')
15
           "thetaphi3" : 3 pairs of (theta,phi) ala GEANT3
16
                         needs 6 values and attribute: units
17
           "newxyz"    : 3 direction vectors, ala Geant4
18
                         needs 9 values (will normalize each vector)
19
     <beampos> :  define the beam origin in the user coordinate system
20
                  3 values define the point 
21
                  if 6 values treat as ( user xyz ) = ( beam xyz )
22
                  but must have set <beamdir> first
23
     <window> :   define the flux window via three <point> elements
24
                  each <point> is a user coord (x,y,z)
25
     <enumax>:    set limits on the maximum enu energy:  lowlimit, fudgefactor
26
                  set limits on wgt: fudgefactor, # points to use
27
                  final reported energy is max(lowlimit,fudgefactor*enumxscan)
28
                  where is 'enumxscan' is the highest energy seen when 
29
                  scanning for x-y weights.
30
     <reuse>:     set # of times an entry is sequentially reused
31
     <upstreamz>: user coord z to push neutrino orgin to
32
                  if abs(z) > 1e30 then leave on the flux window 
33
                              
34
     The units, beamdir, beampos elements define the user coordinate system 
35
        user_dir =  rot * beam_dir
36
        user_pos = (rot * beam_pos ) + user_pos_beamorigin
37
     If more than one <beamdir> or <beampos> is specified the last
38
     one of each is used.  Results of rotations from <beamdir> elements
39
     are not cummulative.
40
            
41
     Parsing of element contents w/ multiple values is liberal about
42
     spacing, line breaks, grouping and separating characters.
43

44
  -->
45

    
46
<gnumi_config>
47
  <param_set name="microboone-numi">
48
    <!-- verbosity of XML parsing -->
49
    <verbose> 3 </verbose>
50
    <!-- setting user units should be first -->
51
    <units> m </units>
52
    <!-- beamdir must come before beam zero position -->
53
    <!-- direction of the beam in the user coordinate system -->
54
    <!-- beam goes up 3 degrees in the coordinate system of the FD -->
55
    <!-- rotation matrix created by sequence of rotations -->
56
    <!-- ** Single dimension at a time transformation ** -->
57
    <!-- <beamdir type="series">
58
       <rotation axis="x" units="rad"> -0.0582977560 </rotation>
59
       <rotation axis="y" units="rad">  0.4038688478 </rotation>
60
    </beamdir> -->
61
    <!-- ** GEANT4 Style Rotation Matrix ** -->
62
    <beamdir type="newxyz">
63
      [   0.921228671    0.0226872648   0.388359401  ]
64
      [   0.00136256111  0.998103714   -0.061539578  ]
65
      [  -0.389019125    0.0572211871   0.919450845  ]
66
    </beamdir>
67
    <!-- position of beam -->
68
    <!-- this is a translation from ( user coords ) to ( beam coords ) -->
69
    <!-- basically a statement of the origin in the detector coordinate system
70
         has the corresponding value in the beam coordinate system -->
71
    <!-- XML config must have already set the rotation matrix -->
72
    <!-- uboone with numi beam - this puts uboone origin at miniboone origin
73
         close enough for now -->
74
    <beampos> ( 0.0, 0.0, 0.0 ) = ( 54.499, 74.461,  677.611) </beampos>
75
    <!-- define flux window in detector coords -->
76
    <!-- need 3 positions to define a rectoid -->
77
    <!-- not yet support beam coords, 
78
         but if/when must have already specified transform -->
79
    <!-- make this window big enough to cover the detector, also offset it from 
80
         the origin towards the beam in x and y -->
81
    
82
    <window>
83
      <point coord="det"> 5, -5, -35 </point>
84
      <point coord="det"> -5, 2, -35 </point>
85
      <point coord="det"> 5, -5, 20  </point>
86
    </window>
87
    
88
    <enumax> 15. 1.05 </enumax> 
89
    <!--first number is estimated max energy, second is
90
        scale factor to multiply max scan energy by -->
91
  </param_set>
92
  
93
</gnumi_config>