Merge pull request #77 from NNPDF/CMS_SINGLETOP_TCH_R_8TEV

CMS_SINGLETOP_TCH_R_8TEV

filters/CMS_SINGLETOP_TCH_R_8TEV.cc

@@ -0,0 +1,63 @@
+/*
+Measurement of the ratio of the single top quark and single antitop quark cross sections in the t-channel @LHC CMS 8 TeV
+There is a single point
+
+LHC-CMS 8 TeV
+---------------
+
+emu events with b-tagged jets (L=19.7 1/fb)
+Archived as: 1403.7366v2
+Published in: JHEP 06 (2014) 090
+R_{t-ch.} = 1.95 ± 0.10 (stat) ± 0.19 (syst)
+
+The systematic uncertainty is composed of the following 9 uncertainties (see Table 4 in the paper for the full breakdown and for more details):
+JES, JER, MET, and pileup
+b-tagging and mis-tag
+Leptop reconstruction/trig.
+QCD multijet estimation
+W+jets, ttbar estimation
+Other backgrounds ratio
+Signal modeling
+PDF uncertainty
+Simulation sample size
+*/
+
+#include "CMS_SINGLETOP_TCH_R_8TEV.h"
+
+void CMS_SINGLETOP_TCH_R_8TEVFilter::ReadData()
+{
+  // Opening file
+  fstream f1;
+  stringstream datafile("");
+  datafile << dataPath() 
+	   << "rawdata/" << fSetName << "/" << fSetName << ".data";
+  f1.open(datafile.str().c_str(), ios::in);
+
+  if (f1.fail()) 
+    {
+      cerr << "Error opening data file " << datafile.str() << endl;
+      exit(-1);
+    }
+
+  // Starting filter
+
+  string line;
+  int idum;
+  double cme;
+
+  getline(f1,line);
+  istringstream lstream(line);
+  lstream >> idum >> cme;
+
+  fKin1[0] = 0.;
+  fKin2[0] = Mt*Mt;          // Top mass
+  fKin3[0] = cme*1000;       // Sqrt(s)
+
+  lstream >> fData[0];       // Central value
+  lstream >> fStat[0];       // Absolute statistical uncertainty
+
+  lstream >> fSys[0][0].add; // Absolute total systematic uncertainty
+  fSys[0][0].mult = fSys[0][0].add*100/fData[0]; // Multiplicative total systematic uncertainty
+  fSys[0][0].type = MULT;
+  fSys[0][0].name = "UNCORR";      
+}

inc/CMS_SINGLETOP_TCH_R_8TEV.h

@@ -0,0 +1,15 @@
+/**
+ *  \class CMS_SINGLETOP_TCH_R_8TEV
+ *  \brief CMS_SINGLETOP_TCH_R_8TEV CommonData processor
+ */
+
+#pragma once
+
+#include "buildmaster_utils.h"
+
+class CMS_SINGLETOP_TCH_R_8TEVFilter: public CommonData {
+public: CMS_SINGLETOP_TCH_R_8TEVFilter():
+  CommonData("CMS_SINGLETOP_TCH_R_8TEV") { ReadData(); }
+private:
+  void ReadData();
+};

meta/CMS_SINGLETOP_TCH_R_8TEV.yaml

@@ -0,0 +1,4 @@
+ndata:  	  1
+nsys:   	  1
+setname:	  "CMS_SINGLETOP_TCH_R_8TEV"
+proctype:	  "INC"

rawdata/CMS_SINGLETOP_TCH_R_8TEV/CMS_SINGLETOP_TCH_R_8TEV.data

@@ -0,0 +1 @@
+1   8.0   1.95   0.10   0.19

src/buildmaster.cc

@@ -66,6 +66,7 @@
 #include "CMSWMU8TEV.h"
 #include "FutureColliders.h"
 #include "CMS_SINGLETOP_TCH_TOT_7TEV.h"
+#include "CMS_SINGLETOP_TCH_R_8TEV.h"
 #include "CMS_SINGLETOP_TCH_R_13TEV.h"
 
 /**
@@ -244,6 +245,7 @@ void InitCommonData(vector<CommonData*>& target)
   target.push_back(new CMSTTBARTOT13TEVFilter());
   target.push_back(new CMSTTBARTOTFilter());
   target.push_back(new CMS_SINGLETOP_TCH_TOT_7TEVFilter());
+  target.push_back(new CMS_SINGLETOP_TCH_R_8TEVFilter());
   target.push_back(new CMS_SINGLETOP_TCH_R_13TEVFilter());
 
   // ************************* ZEUS ******************************