Changes between Version 25 and Version 26 of HowToAddSample

Timestamp:

Jan 5, 2017, 5:22:53 PM (7 years ago)

Author:

Clarence Wret

Comment:

--

HowToAddSample

@@ -8 +8 @@
 
 After this we'll be able to make data/MC comparisons to the T2K CC1pi+ H,,2,,O data-set for NEUT, GENIE, GiBUU and NuWro; all in the same framework. We'll get MC predictions, chi2 values, interaction mode contributions, and more using a consistent signal definition, guaranteed the same across the generators.
+
+'''NUISANCE is currently fully written in C++03 and we would appreciate future contributions to adhere to this standard.'''
 
 {{{
@@ -108 +110 @@
 == Writing the header ==
 
-In the case of our T2K CC1pi+ H,,2,,O data, we're dealing with 1D distributions. They should therefore inherit from the `Measurement1D` class, as mentioned  [#point_base earlier]. The `Measurement1D` class is implemented in `src/FitBase/Measurement1D.cxx`.
-
-
-We will also require a `constructor` and `destructor` for the class, and we'll need to overload `MeasurementBase` methods which define the dependent variable(s) (p,,mu,, in our case) and what our signal is (CC interaction with one muon and one positive pion with no other pions or mesons and any number of nucleons in our case). The `MeasurementBase` functions which we need to overload are `isSignal(FitEvent *event)` and `FillEventVariables(FitEvent *event)`. 
+In the case of our T2K CC1pi+ H,,2,,O data, we're dealing with 1D distributions. They should therefore inherit from the `Measurement1D` base class, as mentioned  [#point_base earlier]. The `Measurement1D` class is implemented in `src/FitBase/Measurement1D.cxx`.
+
+
+We will also require a `constructor` and `destructor` for the class, and we'll need to overload `MeasurementBase` methods which define the dependent variable(s) (p,,mu,, in our case) and what our signal is (CC interaction with one muon and one positive pion with no other pions or mesons and any number of nucleons in our case). The `MeasurementBase` functions which we need to overload are `MeasurementBase::isSignal(FitEvent *event)` and `MeasurementBase::FillEventVariables(FitEvent *event)`. 
 
 The `FitEvent` class is an object which contains information about one single event: all the particles, their kinematics and their status after the interaction, the interaction channel which produced the final state, possible secondary interactions, the interaction target, and so on.  The `FitEvent` class is implemented in `src/FitBase/FitEvent.cxx` and essentially the generator-independent common event format which NUISANCE uses for the generator events.
@@ -138 +140 @@
 }}}
 
+
 == Making the constructor ==
 
-== Specifying the event-level dependent variable(s) ==
-
-== Specifying a signal definition ==
-
-
-= The final implementation =
-`T2K_CC1pip_H2O_XSec_1Dpmu_nu.cxx`:
-
-{{{
-#!cpp
+NUISANCE loads samples through `FCN/SampleList.cxx` in the function `SampleUtils::LoadSample`, which makes a `std::list` of `MeasurementBase` pointers which it later loops over.
+
+Most of the automated setting up is done in `Measurement1D::SetupMeasurement(std::string inputfile, std::string type, FitWeight *rw, std::string fakeDataFile)`, so it makes sense to have a class constructor which passes all of these on to `SetupMeasurement`.
+
+ * inputfile: the location of the input file containing generated events
+ * type: the type of comparison we want to make with this sample
+ * rw: the reweighting engine which may or may not be used to reweight the sample
+ * fakeDataFile: the location of a fake-data input file which is to be used as data instead of the actual data for fake-data studies
+
+
+
+
+
+{{{
 #include "T2K_CC1pip_H2O_XSec_1Dpmu_nu.h"
 
@@ -160 +167 @@
   fName = "T2K_CC1pip_H2O_XSec_1Dpmu_nu";
   fPlotTitles = "; p_{#mu} (GeV/c); d#sigma/dp_{#mu} (cm^{2}/(GeV/c)/nucleon)";
-  EnuMin = 0.;
+  EnuMin = 0.; 
   EnuMax = 10.;
   Measurement1D::SetupMeasurement(inputfile, type, rw, fakeDataFile);
@@ -176 +183 @@
   fScaleFactor = (fEventHist->Integral("width")*1E-38)/double(fNEvents)/TotalIntegratedFlux("width");
 };
+}}}
+
+== Specifying the event-level dependent variable(s) ==
+
+== Specifying a signal definition ==
+
+
+= The final implementation =
+`T2K_CC1pip_H2O_XSec_1Dpmu_nu.cxx`:
+
+{{{
+#!cpp
+#include "T2K_CC1pip_H2O_XSec_1Dpmu_nu.h"
+
+// The muon momentum
+
+//******************************************************************** 
+T2K_CC1pip_H2O_XSec_1Dpmu_nu::T2K_CC1pip_H2O_XSec_1Dpmu_nu(std::string inputfile, FitWeight *rw, std::string  type, std::string fakeDataFile){
+//******************************************************************** 
+
+  fName = "T2K_CC1pip_H2O_XSec_1Dpmu_nu";
+  fPlotTitles = "; p_{#mu} (GeV/c); d#sigma/dp_{#mu} (cm^{2}/(GeV/c)/nucleon)";
+  EnuMin = 0.;
+  EnuMax = 10.;
+  Measurement1D::SetupMeasurement(inputfile, type, rw, fakeDataFile);
+
+  // Data comes in ROOT file
+  // hResultTot is cross-section with all errors
+  // hResultStat is cross-section with stats-only errors
+  // hTruthNEUT is the NEUT cross-section given by experimenter
+  // hTruthGENIE is the GENIE cross-section given by experimenter
+  SetDataFromFile(GeneralUtils::GetTopLevelDir()+"/data/T2K/CC1pip/H2O/nd280data-numu-cc1pi-xs-on-h2o-2015.root","MuMom/hResultTot");
+  SetCovarFromDataFile(GeneralUtils::GetTopLevelDir()+"/data/T2K/CC1pip/H2O/nd280data-numu-cc1pi-xs-on-h2o-2015.root", "MuMom/TotalCovariance", true);
+
+  SetupDefaultHist();
+
+  fScaleFactor = (fEventHist->Integral("width")*1E-38)/double(fNEvents)/TotalIntegratedFlux("width");
+};
 
 //********************************************************************