Changes between Version 2 and Version 3 of HowToJoint

Timestamp:

Jun 13, 2017, 2:18:57 PM (7 years ago)

Author:

Patrick Stowell

Comment:

--

HowToJoint

@@ -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.
 
-= Writing a JOINT input =
-A joint input is simply a list of the normal NUISANCE inputs saved into a text format. Examples are below.
-
-'''neut_miniboone_rhc_numu_numubar.txt'''
+'''Reading JOINT GENIE MINERvA fhc'''
+In this example both a FHC numu and FHC numubar event sample will be added together to create a prediction for the MiniBooNE CCQELike data.
 {{{
-NEUT:/data/events/neut/miniboone/rhc_numubar.root
-NEUT:/data/events/neut/miniboone/rhc_numu.root
+sample  MiniBooNE_CCQELike_XSec_1DQ2_nu  GENIE:(/data/events/genie/minerva/genie_minerva_fhc.numuroot,
 }}}
 
-'''genie_minerva_fhc.joint'''
+It is also possible to add different event modes together. For example say you have generated events in GENIE with the Default model, and you would like to add MEC events too. You could read them in as
 {{{
-GENIE:/data/events/genie/minerva/fhc_numu.root
-GENIE:/data/events/genie/minerva/fhc_numubar.root
-GENIE:/data/events/genie/minerva/fhc_nue.root
-GENIE:/data/events/genie/minerva/fhc_nuebar.root
+sample  MiniBooNE_CCQE_XSec_1DQ2_nu  GENIE:(/data/events/genie/minerva/genie_minerva_fhc_numu_default.root,/data/events/genie/minerva/genie_minerva_fhc_numu_mec.root)
 }}}
 
-= Using a JOINT input =
-A joint input file can then be read into NUISANCE as its own JOINT type as follows
 
-'''Reading JOINT NEUT MB rhc'''
 {{{
-sample  MiniBooNE_CCQELike_XSec_1DQ2_antinu  JOINT:/data/events/neut/miniboone/neut_miniboone_rhc_numu_numubar.txt
+#!div class="important" 
+**WARNING** : When reading in joint inputs the maximum number of events defined by the configuration option "MAXEVENTS" must be set to "-1" otherwise NUISANCE will try and run over the first portion of the event file, and will skip out events in the second input file, completely messing up the relative contribution from each event file. 
+
+In NUISANCE v2r0 there is a check for this, but on older versions there is not. You have been warned!
 }}}
-
-'''Reading JOINT GENIE MINERvA fhc'''
-{{{
-sample  MINERvA_CCQE_XSec_1DQ2_nu  JOINT:/data/events/genie/minerva/genie_minerva_fhc.joint
-}}}