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 73 and Version 74 of HowToAddSample


Ignore:
Timestamp:
Jan 8, 2017, 9:55:34 AM (7 years ago)
Author:
Clarence Wret
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • HowToAddSample

    v73 v74  
    396396The base class function `MeasurementBase::Reconfigure()` specifies `isSignal(FitEvent *)` which is what we will implement in our class. Since we plan on re-using the signal definition for multiple distributions, we should include it in the `src/T2K/T2K_SignalDef.cxx` implementation, as mentioned [#point_signal earlier].
    397397
    398 Reading the [https://arxiv.org/pdf/1605.07964v2.pdf paper], we see the requirement:
     398Reading the [https://arxiv.org/pdf/1605.07964v2.pdf paper], we see the two requirements:
     399
     400 "The CC1π^+^cross section is described by the particles leaving the nucleus, i.e. one muon, one positive pion and any number of nucleons" (
     401
     402and
    399403
    400404 "The analysis presented restricts the kinematic phase-space to the region defined by p,,µ,, > 200 MeV/c, p,,π,,^+^ > 200 MeV/c, cos(θ,,µ,,) > 0.3 and cos(θ,,π+,, ) > 0.3."
    401405
    402 so our signal definition needs to extract these quantities from the `FitEvent` object.
     406so our signal definition needs to extract these quantities from the `FitEvent` object and perform the selection upon them.
    403407
    404408We prefer a clear unambiguous naming scheme, so let's name the namespace function `isCC1pip_T2K_H2O`. The function will need to take the `FitEvent` object along with the `EnuMin` and `EnuMax` defined in the experiment signal definition.
    405409
    406 In the signal definition we can use already declared functions of the namespace `SignalDef`. These live in `src/Utils/SignalDef.cxx`. We can see a `SignalDef::isCC1pi(FitEvent *, int nuPDG, int piPDG, double EnuMin, double EnuMax)`. This helper function searches for:
     410We can use already declared functions of the namespace `SignalDef`: these live in `src/Utils/SignalDef.cxx`. We can see a `SignalDef::isCC1pi(FitEvent *, int nuPDG, int piPDG, double EnuMin, double EnuMax)`, which looks promising. This helper function searches for:
    407411
    408412 * 1 single lepton of type `nuPDG-1`
     
    410414 * makes sure the neutrino energy is between `EnuMin` and `EnuMax`
    411415
    412 We can then add our kinematic requirements on the muon and pion in our `SignalDef::isCC1pip_T2K_H2O(FitEvent *event, double EnuMin, double EnuMax` implementation in `src/T2K/T2K_SignalDef.cxx`.
     416We can then add our kinematic phase-space requirements on the muon and pion in our `SignalDef::isCC1pip_T2K_H2O(FitEvent *event, double EnuMin, double EnuMax` implementation in `src/T2K/T2K_SignalDef.cxx`.
    413417
    414418In summary, we get for `src/T2K/T2K_SignalDef.cxx`: