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 HowToAddSample


Ignore:
Timestamp:
Jan 5, 2017, 2:02:06 PM (8 years ago)
Author:
Clarence Wret
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • HowToAddSample

    v5 v6  
    55This is an informal step-by-step guide on how to add samples into the NUISANCE framework.
    66
    7 '''Versions:''' `NUISANCE v1r0`, `NEUT 5.3.6`
     7For this tutorial I'll be adding the T2K CC1pi+ H2O data. The data comes from [http://t2k-experiment.org/results/nd280data-numu-cc1pi-xs-on-h2o-2015/ the T2K site] and [https://arxiv.org/abs/1605.07964 arxiv]. The data is supplied both in a ROOT file and a number of .csv files; I'll be using the ROOT file here.
    88
    9 For this tutorial I'll be adding the T2K CC1pi+ H2O data. The data comes from [http://t2k-experiment.org/results/nd280data-numu-cc1pi-xs-on-h2o-2015/ the T2K site] and [https://arxiv.org/abs/1605.07964 arxiv]. The data is supplied both in a ROOT file and a number of .csv files; I'll be using the ROOT file here.
     9'''Versions:''' `NUISANCE v1r0`, `NEUT 5.3.6`, `arxiv 1605.07964v2.pdf`
    1010
    1111== Examining the data and choosing distributions ==
     
    1414=== Finding the neutrino flux and generating events ===
    1515
    16 The first issue at hand is to find the flux for the experiment. If we don't have this (or its very difficult to find), there is simply no way to reproduce the results (unless the measurement is a total cross-section without any phase space cuts).
     16The first issue at hand is to find the flux for the experiment. If we don't have this (or its very difficult to find), we cant make a generator prediction---unless the measurement is a total cross-section without any phase space cuts.
    1717
    18 A quick search through the document finds it is documented in reference ![12]. You could also find it in our [wiki:ExperimentFlux flux list].
     18A quick search through the document finds it is documented in reference ![12]. It is also in our [wiki:ExperimentFlux flux list].
    1919
    20 I then generate NEUT 5.3.6 events with a suitable card-file (see [wiki:HowToNeut our NEUT guide]), making sure I've chosen the correct target (water) and flux, and that I've made the model selections I want (e.g. Rein-Sehgal or Berger-Sehgal coherent model).
     20I then generate events in NEUT 5.3.6 with a suitable card-file, see our [wiki:HowToNeut NEUT guide] on how to do this. We need the correct target (water) and flux, and perform the model selections we want (e.g. Rein-Sehgal or Berger-Sehgal coherent model). The procedure is very similar for other generators too.
    2121
    22 I normally generate around 1M events with all interaction modes turned on. This way I'll make sure to get all interaction modes into the topological cross-section (e.g. I'll get a small amount of CCQE events which excite a pion from a secondary nucleon interaction, leading to a CCQE+1pi ~ CC1pi+ final state, which is signal)
     22Aim to generate around 1M events with all interaction modes turned on. This way we make sure to get all interaction modes into the topologically defined cross-section. We get a small amount of CCQE events which excite a pion from the outgoing nucleon interacting in the nucleus to kick out a pion, leading to a CCQE+1pi ~ CC1pi+ final state, which is signal for this particular sample.
    2323
    2424
    2525=== Choosing a distribution ===
    2626
    27 The T2K CC1pi+ H2O data release contains various distributions in FIG 4. In this tutorial I'll look at adding one kinematic distribution and one derived distribution: p_mu and E^{rec}_nu.
     27The T2K CC1pi+ H2O data release contains various distributions in FIG 4. In this tutorial I'll look at adding one kinematic distribution and one derived distribution: p_mu and E^{rec}_nu respectively. INSERT IMAGES
    2828
    29 Generally speaking, we try to add all available distributions in a publication. Some distributions will have detector effects "folded" into them (i.e. they will be raw detector-level data). We can only use these if there is some method which bring detector-level variables (e.g. p_mu seen in the detector) to truth-level variable (e.g. p_mu seen after correcting for the detector effects).
     29In NUISANCE we try to add all available distributions from a publication. '''However''', some distributions will have detector effects "folded" into them (i.e. they will be raw detector-level data). We can only use these if there is some method which bring detector-level variables (e.g. p_mu seen in the detector) to truth-level variable (e.g. p_mu seen after correcting for the detector effects).