nuisance is hosted by Hepforge, IPPP Durham
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.

Changes between Version 5 and Version 6 of HowToNuWro


Ignore:
Timestamp:
Nov 6, 2016, 10:35:55 AM (8 years ago)
Author:
Patrick Stowell
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • HowToNuWro

    v5 v6  
    3535}}}
    3636
     37The output file (specified with '-o') will be a ROOT file containing a series of useful histograms that compare the NuWro predictions to the MiniBooNE CCQE data.
     38{{{
     39$ root nuisance_comparison_nuwro_v12.root
     40$ root[0]: _file0->ls();
     41
     42}}}
     43
     44'''MiniBooNE_CCQE_XSec_1DQ2_data:'''
     45
     46
     47'''MiniBooNE_CCQE_XSec_1DQ2_MC:'''
     48
     49
     50'''MiniBooNE_CCQE_XSec_1DQ2_MC_PDG:'''
     51
     52
     53 
    3754
    3855= Loading !NuWro <v12 events into NUISANCE =
     
    83100}}}
    84101
    85 A list of available dial strings can be found in "$NUWRO/src/reweight/NuWroSyst.h".
     102A list of available dial strings can be found in "$NUWRO/src/reweight/NuWroSyst.h". Comparing the plots below with those shown earlier you can see the difference between the NuWro model predictions with MA = 0.0 and +1.0 sigma.
     103
     104'''MiniBooNE_CCQE_XSec_1DQ2_MC (reweighted):'''
     105
     106
    86107
    87108= Using !NuWro reweight for dial tuning =
    88109As 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.
    89110 
     111
     112'''nuwro_reweighted_minimizer.card:'''
     113{{{
     114# par_type      dial_name    value
     115###########################
     116nuwro_parameter  MaCCQE  +1.0  -4.0  4.0  1.0  FREE   
     117
     118# sample   SAMPLE_ID    TYPE:INPUT_TYPE
     119########################################
     120sample    MiniBooNE_CCQE_XSec_1DQ2_nu    NUWRO:/path/to/nuwro/events.root
     121}}}
     122
     123
     124'''nuwro_reweighted_minimizer.card(without comments):'''
     125{{{
     126nuwro_parameter  MaCCQE  +1.0  -4.0  4.0  1.0  FREE   
     127sample    MiniBooNE_CCQE_XSec_1DQ2_nu    NUWRO:/path/to/nuwro/events.root
     128}}}
     129
     130
     131This can be ran using the NUISANCE minimizer app "nuismin":
     132{{{
     133$ nuismin -c nuwro_reweighted_minimizer.card -o nuwro_minimizer_results.root
     134}}}
     135
     136The contents of the nuismin outputfile will be the same as that in the comparison app, but with the predictions now shown for the best fit result. Also included are a series of histograms related to the tuning study (covariance, fit_results, etc)
     137{{{
     138$ root nuwro_minimizer_results.root
     139$ root[0]: _file0->ls();
     140}}}
    90141
    91142