close
Warning:
Can't synchronize with repository "(default)" ("(default)" is not readable or not a Git repository.). Look in the Trac log for more information.
- Timestamp:
-
Jan 8, 2017, 4:13:29 PM (8 years ago)
- Author:
-
Clarence Wret
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v19
|
v20
|
|
81 | 81 | The output of `neutroot2`, contains a ROOT file containing `nevents` `neutvect` objects which has all the information about each event in it. This is read in by NUISANCE to generate predictions. `neutroot2` also saves the predicted flux (`flux_numu`) and event rate (`evtrt_numu`) predictions automatically into the output file so NUISANCE already has everything it needs. |
82 | 82 | |
83 | | To read NEUT inputs the NEUT type just needs to be specified in the NUISANCE input card file, e.g. for MiniBooNE CCQE: |
84 | | `sample MiniBooNE_CCQE_XSec_1DQ2_nu NEUT:/path/to/neut/input/file.root` |
| 83 | To read NEUT inputs the NEUT type just needs to be specified in the NUISANCE input card file, e.g. for MiniBooNE CCQE |
| 84 | |
| 85 | {{{ |
| 86 | sample MiniBooNE_CCQE_XSec_1DQ2_nu NEUT:/path/to/neut/input/file.root |
| 87 | }}} |
85 | 88 | |
86 | 89 | |
… |
… |
|
109 | 112 | The user specifies the parameter(s) and the range which should be varied as: |
110 | 113 | |
111 | | `neut_parameter MaCCQE 0.0 -4.0 4.0 1.0 FREE` |
| 114 | {{{ |
| 115 | neut_parameter MaCCQE 0.0 -4.0 4.0 1.0 FREE |
| 116 | }}} |
112 | 117 | |
113 | 118 | where `MaCCQE` is the parameter to be varied, `0.0` is the starting values, `-4.0` is the minimum value, `4.0` is the maximum value, `1.0` is the step-size to take, and `FREE` means the parameter is free to vary in the minimisation. |
… |
… |
|
117 | 122 | The whole list of parameters can be found in `$NEUT_ROOT/src/reweight/NSyst.h` and the `return`ed `std::string`s (e.g. `MaCCQE` and `CA5RES`). |
118 | 123 | |
| 124 | {{{ |
| 125 | #!div class="important" |
119 | 126 | The starting, minimum, maximum and step values are all in units of 1-sigma uncertainties. So `-4.0` for the minimum means -4*1-sigma ranges, and `0.0` for the starting value means the nominal value nominal. |
120 | 127 | |
… |
… |
|
126 | 133 | |
127 | 134 | The output of the fit is in the units of `par_var`, not `par_value`. So an output for `MaCCQE` of `0.7` would equate to `1.21*(1+1*0.165) = 1.35` GeV. |
128 | | |
| 135 | }}} |
129 | 136 | |
130 | 137 | = Joint Flux Inputs = |