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:
-
Nov 15, 2016, 12:43:23 PM (8 years ago)
- Author:
-
Patrick Stowell
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v1
|
v2
|
|
5 | 5 | |
6 | 6 | = How it works = |
7 | | For each entry in a given list of input files NUISANCE will temporarily create an InputHandler and use it to extract the total rate, flux, and number of entries. |
| 7 | For each entry in a given list of input files NUISANCE will temporarily create an !InputHandler and use it to extract the total rate, flux, and number of entries. |
8 | 8 | |
9 | | A TChain is then formed for the provided generator type and each input is added into the list. NUISANCE keeps track of how many entries were in each TTree and uses this to calculate a separate scaling factor for each set of events. When a JOINT input is being used this scale factor is assigned as the FitEvent.InputWeight value to allow it to be propagated to the event predictions. |
| 9 | A TChain is then formed for the provided generator type and each input is added into the list. NUISANCE keeps track of how many entries were in each TTree and uses this to calculate a separate scaling factor for each set of events. When a JOINT input is being used this scale factor is assigned as the !FitEvent.!InputWeight value to allow it to be propagated to the event predictions. |
10 | 10 | |
11 | 11 | To get the best results and ensure that the statistical error on each sample is correct we suggest that the number of events in each sample is proportional to the relative flux integral of each component. For example if generating 1 million events with a numu + numubar flux with integrals of 150 and 500 respectively (arbritrary norm), then the number of events in each sample should be numu=230,769 and numubar=769,230. |
… |
… |
|
28 | 28 | }}} |
29 | 29 | |
| 30 | = Using a JOINT input = |
| 31 | A joint input file can then be read into NUISANCE as its own JOINT type as follows |
| 32 | |
| 33 | '''Reading JOINT NEUT MB rhc''' |
| 34 | {{{ |
| 35 | sample MiniBooNE_CCQELike_XSec_1DQ2_antinu JOINT:/data/events/neut/miniboone/neut_miniboone_rhc_numu_numubar.txt |
| 36 | }}} |
| 37 | |
| 38 | '''Reading JOINT GENIE MINERvA fhc''' |
| 39 | {{{ |
| 40 | sample MINERvA_CCQE_XSec_1DQ2_nu JOINT:/data/events/genie/minerva/genie_minerva_fhc.joint |
| 41 | }}} |