1 | | = [ page being built ] = |
2 | | This page aims to outline how to use NuWro to create a sample of interactions on your favourite target. We'll also outline how to feed this generated sample through the NUISANCE framework and get good results. |
| 1 | [[PageOutline(1-3,,pullout,unnumbered)]] |
| 2 | |
| 3 | = !NuWro Inputs = |
| 4 | This page aims to outline how to use !NuWro to create a sample of interactions on your favourite target. We'll also outline how to feed this generated sample through the NUISANCE framework and get good results. |
| 5 | |
| 6 | The latest version of the !NuWro generator is setup to automatically save the predicted event rate and flux distributions for the generated events. This is everything NUISANCE needs to make properly normalised cross-section predictions with the standard !NuWro output. |
| 7 | |
| 8 | The latest version of !NuWro with recently developed reweight engine can be obtained from here: |
| 9 | |
| 10 | Note: If you would like to use events generated with an older version of !NuWro this is still possible. Please see: |
| 11 | |
| 12 | |
| 13 | = Generate !NuWro Events for NUISANCE = |
| 14 | |
| 15 | |
| 16 | = Loading !NuWro v12+ events into NUISANCE = |
| 17 | NUISANCE's input handlers are setup to automatically load in the standard output of !NuWro v12 onwards. |
| 18 | Once you have generated events the MC sample can be loaded in by specifying the "NUWRO" type in the NUISANCE sample input field. An example of this is shown in the following card file. |
| 19 | |
| 20 | sample MiniBooNE_CCQE_XSec_1DQ2_nu NUWRO:/path/to/nuwro/events.root |
| 21 | |
| 22 | |
| 23 | = Loading !NuWro <v12 events into NUISANCE = |
| 24 | For versions of !NuWro less than v12 there are no flux and event rate histograms saved into the generator output. It is still possible to use these events but NUISANCE must reconstruct the histograms at runtime (a slow process). |
| 25 | |
| 26 | This must be done by using the flag input.regen_nuwro_plots=1 at runtime. When this is set NUISANCE will extract the flux histogram spectrum from the !NuWro parameters class. It will then bin all the events as a function of Enu and use the !NuWro event weight to normalise this event spectrum to a properly normalised predicted event rate for the given input flux. |
| 27 | |
| 28 | The following cardfile is an example of how to setup this type of input for the MiniBooNE CCQE dataset. |
| 29 | |
| 30 | = Using !NuWro reweight for comparisons = |
| 31 | Support for the !NuWro reweight package is built into NUISANCE. It has immediate access to any reweight dial available in the version of !NuWro you built NUISANCE against. An example of a nuiscomp cardfile that will generate predictions with alternate model predictions is shown below: |
| 32 | |
| 33 | |
| 34 | |
| 35 | = Using !NuWro reweight for dial tuning = |
| 36 | As for all other reweight engines any dial can be treated as a free parameter for dial tuning studies in nuismin. The following cardfile when used with nuismin will vary the CCQE and Resonant axial mass parameters until a best fit is found to the MiniBooNE CCQELike dataset. |
| 37 | |
| 38 | |
| 39 | |
| 40 | |
| 41 | |
| 42 | |
| 43 | |
| 44 | |
| 45 | |
| 46 | |
| 47 | |
| 48 | |
| 49 | |