Changes between Version 47 and Version 48 of HowToAddSample

Timestamp:

Jan 6, 2017, 4:08:35 PM (8 years ago)

Author:

Clarence Wret

Comment:

--

HowToAddSample

@@ -215 +215 @@
 The data and covariance matrix are available in a ROOT file for the sample, so I used `Measurement1D::SetDataFromFile` to set up the data and `Measurement1D::SetCovarFromDataFile` to set up the covariances.
 
+`fName` and the other `std::string`s were set to reasonable names and `fScaleFactor` was 
+
 The `constructor` for the class now looks like:
 
@@ -229 +231 @@
   fPlotTitles = "; p_{#mu} (GeV/c); d#sigma/dp_{#mu} (cm^{2}/(GeV/c)/nucleon)";
   EnuMin = 0.; 
-  EnuMax = 10.;
+  EnuMax = 100.;
   Measurement1D::SetupMeasurement(inputfile, type, rw, fakeDataFile);
 
@@ -325 +327 @@
 }}}
 
-We should now write a helper function in the `FitUtils` namespace to get `E^rec^,,nu,,` from the muon, pion and neutrino `TLorentzVector`s. [arxiv 1605.07964v2.pdf The paper] specifies the exact form in equation 1. The `FitUtils` implementation lives in `src/Utils/FitUtils.cxx`.
+We should now write a helper function in the `FitUtils` namespace to get `E^rec^,,nu,,` from the muon, pion and neutrino `TLorentzVector`s. [https://arxiv.org/pdf/1605.07964v2.pdf The paper] specifies the exact form in equation 1. The `FitUtils` implementation lives in `src/Utils/FitUtils.cxx`.
 
 {{{
@@ -373 +375 @@
 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].
 
-Reading the [arxiv 1605.07964v2.pdf paper], we see the requirement:
+Reading the [https://arxiv.org/pdf/1605.07964v2.pdf paper], we see the requirement:
 
  "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."
 
-
+so our signal definition needs to extract these quantities from the `FitEvent` object.
+
+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.
+
+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:
+
+ * 1 single lepton of type `nuPDG-1` 
+ * 1 single pion of type `piPDG`
+ * makes sure the neutrino energy is between `EnuMin` and `EnuMax`
+
+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`.
+
+In summary, we get for `src/T2K/T2K_SignalDef.cxx`:
+
+{{{
+namespace SignalDef {
+
+  // T2K H2O signal definition
+  bool isCC1pip_T2K_H2O(FitEvent *event, double EnuMin,
+                                 double EnuMax) {
+
+    if (!isCC1pi(event, 14, 211, EnuMin, EnuMax)) return false;
+
+    TLorentzVector Pnu = event->GetHMISParticle(14)->fP; // Can also use FitEvent::GetNeutrinoIn()
+    TLorentzVector Pmu = event->GetHMFSParticle(13)->fP;
+    TLorentzVector Ppip = event->GetHMFSParticle(211)->fP;
+
+    double p_mu = FitUtils::p(Pmu) * 1000;
+    double p_pi = FitUtils::p(Ppip) * 1000;
+    double cos_th_mu = cos(FitUtils::th(Pnu, Pmu));
+    double cos_th_pi = cos(FitUtils::th(Pnu, Ppip));
+
+    if (p_mu <= 200 || p_pi <= 200 || cos_th_mu <= 0.3 || cos_th_pi <= 0.3) {      
+      return false;
+    }
+
+    return true;
+  };
+
+}
+
+}}}
+
+Then back in `src/T2K/T2K_CC1pip_H2O_XSec_1Dpmu_nu.cxx` we simply do:
+
+{{{
+//******************************************************************** 
+// Beware: The H2O analysis has different signal definition to the CH analysis!
+bool T2K_CC1pip_H2O_XSec_1Dpmu_nu::isSignal(FitEvent *event) {
+//******************************************************************** 
+  return SignalDef::isCC1pip_T2K_H2O(event, EnuMin, EnuMax);
+}
+}}}
+
+making sure our [#point_header header file] has `#include "T2K_SignalDef.h"`.
 
 
 = The final implementation =
+
+Here's a summary of the implementation for reference:
+
 `T2K_CC1pip_H2O_XSec_1Dpmu_nu.cxx`:
 
@@ -460 +519 @@
 }}}
 
+
+
 = Make NUISANCE aware of the sample =
 FCN/SampleList