|
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.
- Timestamp:
-
Jan 8, 2017, 9:55:34 AM (8 years ago)
- Author:
-
Clarence Wret
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v73
|
v74
|
|
396 | 396 | The 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]. |
397 | 397 | |
398 | | Reading the [https://arxiv.org/pdf/1605.07964v2.pdf paper], we see the requirement: |
| 398 | Reading 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 | |
| 402 | and |
399 | 403 | |
400 | 404 | "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." |
401 | 405 | |
402 | | so our signal definition needs to extract these quantities from the `FitEvent` object. |
| 406 | so our signal definition needs to extract these quantities from the `FitEvent` object and perform the selection upon them. |
403 | 407 | |
404 | 408 | We 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. |
405 | 409 | |
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: |
| 410 | 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)`, which looks promising. This helper function searches for: |
407 | 411 | |
408 | 412 | * 1 single lepton of type `nuPDG-1` |
… |
… |
|
410 | 414 | * makes sure the neutrino energy is between `EnuMin` and `EnuMax` |
411 | 415 | |
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`. |
| 416 | We 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`. |
413 | 417 | |
414 | 418 | In summary, we get for `src/T2K/T2K_SignalDef.cxx`: |
|