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


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

--

Legend:

Unmodified
Added
Removed
Modified
  • HowToJoint

    v1 v2  
    55
    66= How it works =
    7 For 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.
     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.
    88
    9 A 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.
     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.
    1010
    1111To 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.
     
    2828}}}
    2929
     30= Using a JOINT input =
     31A joint input file can then be read into NUISANCE as its own JOINT type as follows
     32
     33'''Reading JOINT NEUT MB rhc'''
     34{{{
     35sample  MiniBooNE_CCQELike_XSec_1DQ2_antinu  JOINT:/data/events/neut/miniboone/neut_miniboone_rhc_numu_numubar.txt
     36}}}
     37
     38'''Reading JOINT GENIE MINERvA fhc'''
     39{{{
     40sample  MINERvA_CCQE_XSec_1DQ2_nu  JOINT:/data/events/genie/minerva/genie_minerva_fhc.joint
     41}}}