Changes between Version 4 and Version 5 of HowToUseNUISCOMP-GENIE

Timestamp:

Jun 12, 2017, 6:34:00 PM (7 years ago)

Author:

Patrick Stowell

Comment:

--

HowToUseNUISCOMP-GENIE

@@ -346 +346 @@
 }}}
 
-With this new card file we can then run our comparisons again, 
-
+With this new card file we can then run our comparisons again, but this time save the output to a different file
+{{{
+$ nuiscomp -c genie_tutorial.card -o genie_samples_reweighted.root
+
+[LOG Fitter]: Starting nuiscomp.exe
+
+...
+
+[LOG Fitter]: Setting up nuiscomp
+[LOG Fitter]: Number of parameters :  2
+[LOG Fitter]: Read genie_parameter : MaCCRES = 1 :
+[LOG Fitter]: Read genie_parameter : NormCCRES = -1 :
+
+...
+
+[LOG Fitter]: Setting up FitWeight Engine
+[LOG Fitter]: Registed Dial Enum : MaCCRES 5 5011
+[LOG Fitter]: Setting up GENIE RW : genierw
+
+...
+
+[LOG Reconf]:--- Starting Reconfigure iter. 0
+[LOG Reconf]:--- Event Manager Reconfigure
+[LOG Reconf]:--- MiniBooNE_CC1pip_XSec_1DTpi_nu : Processed 0 events. [M, W] = [33, 1]
+[LOG Reconf]:--- MiniBooNE_CC1pip_XSec_1DTpi_nu : Processed 500000 events. [M, W] = [13, 1]
+[LOG Reconf]:--- MiniBooNE_CC1pip_XSec_1DTpi_nu : Processed 1000000 events. [M, W] = [13, 1]
+[LOG Reconf]:--- MiniBooNE_CC1pip_XSec_1DTpi_nu : Processed 1500000 events. [M, W] = [11, 1.17986]
+[LOG Reconf]:--- MiniBooNE_CC1pip_XSec_1DTpi_nu : Processed 2000000 events. [M, W] = [51, 1]
+
+...
+
+[LOG Minmzr]:- Writing each of the data classes...
+[LOG Sample]:-- Written Histograms: MiniBooNE_CC1pip_XSec_1DTpi_nu
+[LOG Sample]:-- Written Histograms: MiniBooNE_CC1pip_XSec_1DTu_nu
+[LOG Fitter]: ------------------------------------ -
+[LOG Fitter]: Comparison Complete.
+[LOG Fitter]: ------------------------------------ -
+}}}
+
+You might have noticed that now during the event processing stage, the values [M, W] correspond to the [Mode, Weight] for the event processed at that point, and that now we have added reweight dials, the weights are not always equal to 1.0.
+
+Now lets open up both our comparison files and compare the outputs
+{{{
+$ root genie_samples.root genie_samples_reweighted.root
+root [2] TBrowser b
+}}}
+
+We now have two files with the exact same structure, the only differences between these will be the values given in each of the MC histograms.
+
+
+
+If we compare the MiniBooNE_CC1pip_XSec_1DTpi_nu_MC histograms on a single plot against the data we can see that the normalisation of the reweighed prediction has shifted upwards, and the likelihood for that prediction given the data has gotten worse as a result.
+
+
+
+
+
+