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 Initial Version and Version 1 of HowToJoint


Ignore:
Timestamp:
Nov 15, 2016, 12:40:24 PM (8 years ago)
Author:
Patrick Stowell
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • HowToJoint

    v1 v1  
     1= JOINT Inputs =
     2
     3Some generators do not support combined beams or targets. The JOINT input type allows NUISANCE to combined different beam files at run time producing a composite input.
     4NOTE: Inputs from different generators cannot currently be merged, you cannot mix a NEUT and GENIE MC file into a single JOINT input.
     5
     6= How it works =
     7For 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
     9A 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
     11To 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.
     12
     13= Writing a JOINT input =
     14A joint input is simply a list of the normal NUISANCE inputs saved into a text format. Examples are below.
     15
     16'''neut_miniboone_rhc_numu_numubar.txt'''
     17{{{
     18NEUT:/data/events/neut/miniboone/rhc_numubar.root
     19NEUT:/data/events/neut/miniboone/rhc_numu.root
     20}}}
     21
     22'''genie_minerva_fhc.joint'''
     23{{{
     24GENIE:/data/events/genie/minerva/fhc_numu.root
     25GENIE:/data/events/genie/minerva/fhc_numubar.root
     26GENIE:/data/events/genie/minerva/fhc_nue.root
     27GENIE:/data/events/genie/minerva/fhc_nuebar.root
     28}}}
     29