コード例 #1
0
ファイル: objects.cpp プロジェクト: 5432935/crossbridge
int Table::operator == (Object &tblobj) 
{
	char *a;
	if (strcmp (a=tblobj.Type (), Type ()))						// Type check
		return 0;
	Table *tblptr = (Table *) &tblobj;
	if (tblptr->Nelem () == nelem)								// Count check
	{
	Object *b;
		for (int i = 0; i < nelem; i++)
			if (*(b=tblptr->Get (i)) != *(Get(i)))				// Elementwise check
				return 0;
	}
	else
		return 0;
	return 1;
}
コード例 #2
0
int main (int argc, char *argv[])
{
    loadBDTSignalRegions();
    loadCnCSignalRegions();

    printBoxedMessage("Starting plot generation");

    // ####################
    // ##   Init tools   ##
    // ####################

    // Create a sonic Screwdriver
    SonicScrewdriver screwdriver;

    // #########################################################
    // ##   Create ProcessClasses (and associated datasets)   ##
    // #########################################################

    screwdriver.AddProcessClass("1ltop",    "1l top",                            "background",  kRed-7);
    screwdriver.AddProcessClass("ttbar_2l", "t#bar{t} #rightarrow l^{+}l^{-}",   "background",  kCyan-3);
    screwdriver.AddProcessClass("W+jets",   "W+jets",                            "background",  kOrange-2);
    screwdriver.AddProcessClass("rare",     "rare",                              "background",  kMagenta-5);
    screwdriver.AddProcessClass("data",     "data",                              "data",        COLORPLOT_BLACK);

    // ##########################
    // ##    Create Regions    ##
    // ##########################

    string globalTag = argv[1];
    string globalLabel = argv[2];

    vector<string> signalRegionsTagList;
    vector<string> signalRegionsLabelList;
    for (int i = 3 ; i < argc ; i++)
    {
        signalRegionsTagList.push_back(argv[i]);
        string label = signalRegionLabel(argv[i],"root");
        replace(label, "BDT T2tt-",     "BDT ");
        replace(label, "BDT T2bw025-",  "BDT ");
        replace(label, "BDT T2bw050-",  "BDT ");
        replace(label, "BDT T2bw075-",  "BDT ");
        replace(label, "T2tt, ",        ""    );
        replace(label, "T2bw x=0.25, ", ""    );
        replace(label, "T2bw x=0.50, ", ""    );
        replace(label, "T2bw x=0.75, ", ""    );
        signalRegionsLabelList.push_back(label);
    }

    for (unsigned int i = 0 ; i < signalRegionsTagList.size() ; i++)
    {
        screwdriver.AddRegion(signalRegionsTagList[i],
                              signalRegionsLabelList[i],
                              0);
    }

    // ##########################
    // ##   Create Channels    ##
    // ##########################

    screwdriver.AddChannel("signalRegion_MTpeak",   globalLabel+";Preselection;M_{T} peak",   0);
    screwdriver.AddChannel("signalRegion_MTtail",   globalLabel+";Preselection;M_{T} tail",   0);

    screwdriver.AddChannel("0btag_MTpeak",          globalLabel+";0 b-tag Control Region;M_{T} peak",        0);
    screwdriver.AddChannel("0btag_MTtail",          globalLabel+";0 b-tag Control Region;M_{T} tail",        0);

    screwdriver.AddChannel("reversedVeto_MTpeak",   globalLabel+";Reversed veto Control Region;M_{T} peak",  0);
    screwdriver.AddChannel("reversedVeto_MTtail",   globalLabel+";Reversed veto Control Region;M_{T} tail",  0);

    screwdriver.AddChannel("2leptons",              globalLabel+";2 leptons Control Region;No M_{T} cut",                0);
    screwdriver.AddChannel("2leptons_MTpeak",       globalLabel+";2 leptons Control Region;M_{T} peak",      0);
    screwdriver.AddChannel("2leptons_MTtail",       globalLabel+";2 leptons Control Region;M_{T} tail",      0);

    // ########################
    // ##   Create Figures   ##
    // ########################

    screwdriver.AddFigurePerProcess("yield","Yield", "logY");

    // ########################################
    // ##       Create histograms and        ##
    // ##  schedule type of plots to produce ##
    // ########################################

    // Schedule plots
    screwdriver.SchedulePlots("1DDataMCComparisonFigure");

    // Config plots

    screwdriver.SetGlobalStringOption("1DStackFigurePerProcess",  "includeSignal",                    "stack");
    screwdriver.SetGlobalFloatOption ("1DStackFigurePerProcess",  "factorSignal",                     1.0    );

    screwdriver.SetGlobalStringOption("Plot", "infoTopRight", "CMS Internal");
    screwdriver.SetGlobalStringOption("Plot", "infoTopLeft",  "#sqrt{s} = 8 TeV, L = 19.5 fb^{-1}");

    screwdriver.SetGlobalBoolOption("Plot", "exportPdf", true);
    screwdriver.SetGlobalBoolOption("Plot", "exportEps", false);
    screwdriver.SetGlobalBoolOption("Plot", "exportPng", false);

    // #########################
    // ##   Set the figures   ##
    // #########################

    // Loop on the signal regions and set the figure values...
    for (unsigned int i = 0 ; i < signalRegionsTagList.size() ; i++)
    {
        // Read input tables
        Table rawYields    = Table(string("./results/latest/prediction/rawYieldTables/")+signalRegionsTagList[i]+".tab");
        Table scaleFactors = Table(string("./results/latest/prediction/scaleFactors/")  +signalRegionsTagList[i]+".tab");
        Table tableMTtailCorrection = Table("./results/latest/MTtailCorrection/SF_MTtail.tab");

        Figure SF_pre           = scaleFactors.Get("value","SF_pre");
        Figure SF_post          = scaleFactors.Get("value","SF_post");
        Figure SF_0btag         = scaleFactors.Get("value","SF_0btag");
        Figure SF_vetopeak      = scaleFactors.Get("value","SF_vetopeak");

        string signalRegionTag_ = signalRegionsTagList[i];
        if (findSubstring(signalRegionTag_,"BDT"))
        {
            size_t pos;
            pos = signalRegionTag_.find("_low");
            if (pos != string::npos) signalRegionTag_ = signalRegionTag_.substr(0,pos);
            pos = signalRegionTag_.find("_medium");
            if (pos != string::npos) signalRegionTag_ = signalRegionTag_.substr(0,pos);
            pos = signalRegionTag_.find("_high");
            if (pos != string::npos) signalRegionTag_ = signalRegionTag_.substr(0,pos);
        }
        Figure SF_MTtail_Wjets = tableMTtailCorrection.Get("SFR_Wjets",signalRegionTag_);
        Figure SF_MTtail_1ltop = tableMTtailCorrection.Get("SFR_1ltop",signalRegionTag_);

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "signalRegion_MTpeak", rawYields.Get("signalRegion_MTpeak","1ltop"   ) * SF_post);
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "signalRegion_MTpeak", rawYields.Get("signalRegion_MTpeak","ttbar_2l") * SF_pre );
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "signalRegion_MTpeak", rawYields.Get("signalRegion_MTpeak","W+jets"  ) * SF_post);
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "signalRegion_MTpeak", rawYields.Get("signalRegion_MTpeak","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "signalRegion_MTpeak", rawYields.Get("signalRegion_MTpeak","data"    ));

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "signalRegion_MTtail", rawYields.Get("signalRegion_MTtail","1ltop"   ) * SF_post * SF_MTtail_1ltop);
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "signalRegion_MTtail", rawYields.Get("signalRegion_MTtail","ttbar_2l") * SF_pre                   );
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "signalRegion_MTtail", rawYields.Get("signalRegion_MTtail","W+jets"  ) * SF_post * SF_MTtail_Wjets);
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "signalRegion_MTtail", rawYields.Get("signalRegion_MTtail","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "signalRegion_MTtail", rawYields.Get("signalRegion_MTtail","data"    ));

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "0btag_MTpeak",        rawYields.Get("0btag_MTpeak","1ltop"   ) * SF_0btag);
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "0btag_MTpeak",        rawYields.Get("0btag_MTpeak","ttbar_2l"));
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "0btag_MTpeak",        rawYields.Get("0btag_MTpeak","W+jets"  ) * SF_0btag);
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "0btag_MTpeak",        rawYields.Get("0btag_MTpeak","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "0btag_MTpeak",        rawYields.Get("0btag_MTpeak","data"    ));

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "0btag_MTtail",        rawYields.Get("0btag_MTtail","1ltop"   ) * SF_0btag * SF_MTtail_1ltop);
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "0btag_MTtail",        rawYields.Get("0btag_MTtail","ttbar_2l"));
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "0btag_MTtail",        rawYields.Get("0btag_MTtail","W+jets"  ) * SF_0btag * SF_MTtail_Wjets);
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "0btag_MTtail",        rawYields.Get("0btag_MTtail","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "0btag_MTtail",        rawYields.Get("0btag_MTtail","data"    ));

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "reversedVeto_MTpeak", rawYields.Get("reversedVeto_MTpeak","1ltop"   ) * SF_vetopeak);
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "reversedVeto_MTpeak", rawYields.Get("reversedVeto_MTpeak","ttbar_2l") * SF_pre );
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "reversedVeto_MTpeak", rawYields.Get("reversedVeto_MTpeak","W+jets"  ) * SF_vetopeak);
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "reversedVeto_MTpeak", rawYields.Get("reversedVeto_MTpeak","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "reversedVeto_MTpeak", rawYields.Get("reversedVeto_MTpeak","data"    ));

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "reversedVeto_MTtail", rawYields.Get("reversedVeto_MTtail","1ltop"   ) * SF_vetopeak * SF_MTtail_1ltop);
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "reversedVeto_MTtail", rawYields.Get("reversedVeto_MTtail","ttbar_2l") * SF_pre );
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "reversedVeto_MTtail", rawYields.Get("reversedVeto_MTtail","W+jets"  ) * SF_vetopeak * SF_MTtail_Wjets);
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "reversedVeto_MTtail", rawYields.Get("reversedVeto_MTtail","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "reversedVeto_MTtail", rawYields.Get("reversedVeto_MTtail","data"    ));

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "2leptons_MTpeak",     rawYields.Get("2leptons_MTpeak","1ltop"   ));
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "2leptons_MTpeak",     rawYields.Get("2leptons_MTpeak","ttbar_2l"));
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "2leptons_MTpeak",     rawYields.Get("2leptons_MTpeak","W+jets"  ));
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "2leptons_MTpeak",     rawYields.Get("2leptons_MTpeak","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "2leptons_MTpeak",     rawYields.Get("2leptons_MTpeak","data"    ));

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "2leptons_MTtail",     rawYields.Get("2leptons_MTtail","1ltop"   ));
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "2leptons_MTtail",     rawYields.Get("2leptons_MTtail","ttbar_2l"));
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "2leptons_MTtail",     rawYields.Get("2leptons_MTtail","W+jets"  ));
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "2leptons_MTtail",     rawYields.Get("2leptons_MTtail","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "2leptons_MTtail",     rawYields.Get("2leptons_MTtail","data"    ));

        screwdriver.SetFigure("yield", "1ltop",    signalRegionsTagList[i], "2leptons",            rawYields.Get("2leptons","1ltop"   ));
        screwdriver.SetFigure("yield", "ttbar_2l", signalRegionsTagList[i], "2leptons",            rawYields.Get("2leptons","ttbar_2l"));
        screwdriver.SetFigure("yield", "W+jets",   signalRegionsTagList[i], "2leptons",            rawYields.Get("2leptons","W+jets"  ));
        screwdriver.SetFigure("yield", "rare",     signalRegionsTagList[i], "2leptons",            rawYields.Get("2leptons","rare"    ));
        screwdriver.SetFigure("yield", "data",     signalRegionsTagList[i], "2leptons",            rawYields.Get("2leptons","data"    ));
    }

    // ##############################
    // ##   Make and write plots   ##
    // ##############################

    cout << endl;
    cout << "   > Making plots..." << endl;
    screwdriver.MakePlots();
    cout << "   > Saving plots..." << endl;
    screwdriver.WritePlots("./results/latest/summaryPlots/yields/"+globalTag+"/");

    printBoxedMessage("Plot generation completed");

    // #############################
    // ##   Post-plotting tests   ##
    // #############################

    printBoxedMessage("Program done.");
    return (0);
}
int main (int argc, char *argv[])
{

  printBoxedMessage("Starting plot generation");

  // ####################
  // ##   Init tools   ##
  // ####################

     // Create a sonic Screwdriver
      SonicScrewdriver screwdriver;

     // ##########################
     // ##   Create Variables   ##
     // ##########################

     screwdriver.AddVariable("MT",             "M_{T}",                   "GeV",    40,0,400,       &(myEvent.MT),                    "logY");

     // #########################################################
     // ##   Create ProcessClasses (and associated datasets)   ##
     // #########################################################

        screwdriver.AddProcessClass("1ltop", "1l top",                             "background",kRed-7);
            #ifdef USING_TTBAR_POWHEG
                screwdriver.AddDataset("ttbar_powheg",                "1ltop",  0, 0);
            #endif
            #ifdef USING_TTBAR_MADGRAPH
                screwdriver.AddDataset("ttbar_madgraph_1l",             "1ltop",  0, 0);
            #endif
            //screwdriver.AddDataset("ttbar_madgraph_scaledown",    "1ltop",  0, 0);
            //screwdriver.AddDataset("ttbar_madgraph_scaleup",      "1ltop",  0, 0);
            //screwdriver.AddDataset("ttbar_madgraph_matchingdown", "1ltop",  0, 0);
            //screwdriver.AddDataset("ttbar_madgraph_matchingup",   "1ltop",  0, 0);
            //screwdriver.AddDataset("ttbar_madgraph_mass166-5",    "1ltop",  0, 0);
            //screwdriver.AddDataset("ttbar_madgraph_mass178-5",    "1ltop",  0, 0);

            screwdriver.AddDataset("singleTop_st",                  "1ltop",  0, 0);


     screwdriver.AddProcessClass("ttbar_2l", "t#bar{t} #rightarrow l^{+}l^{-}", "background",kCyan-3);
            #ifdef USING_TTBAR_MADGRAPH
                screwdriver.AddDataset("ttbar_madgraph_2l",   "ttbar_2l",  0, 0);
            #endif


     screwdriver.AddProcessClass("W+jets",   "W+jets",                          "background",kOrange-2);
             screwdriver.AddDataset("W+jets",    "W+jets", 0, 0);

     screwdriver.AddProcessClass("rare",   "rare",                              "background",kMagenta-5);
             screwdriver.AddDataset("rare",   "rare", 0, 0);

     screwdriver.AddProcessClass("data",   "data",                              "data",COLORPLOT_BLACK);
            screwdriver.AddDataset("SingleElec",   "data", 0, 0);
            screwdriver.AddDataset("SingleMuon",   "data", 0, 0);
            screwdriver.AddDataset("DoubleElec",   "data", 0, 0);
            screwdriver.AddDataset("DoubleMuon",   "data", 0, 0);
            screwdriver.AddDataset("MuEl",         "data", 0, 0);

     // ##########################
     // ##    Create Regions    ##
     // ##########################

     screwdriver.AddRegion("0btag_noMTCut",         "Preselection;0 b-tag Control Region, no M_{T} cut",    &goesIn0BtagControlRegion          );

     // ##########################
     // ##   Create Channels    ##
     // ##########################

     screwdriver.AddChannel("singleLepton", "e/#mu-channels",            &goesInSingleLeptonChannel);

  // ########################################
  // ##       Create histograms and        ##
  // ##  schedule type of plots to produce ##
  // ########################################

     // Create histograms
     screwdriver.Create1DHistos();

     // Schedule plots
     screwdriver.SchedulePlots("1DDataMCComparison");

     // Config plots

     screwdriver.SetGlobalStringOption("DataMCComparison",  "includeSignal",                    "stack");
     screwdriver.SetGlobalFloatOption ("DataMCComparison",  "factorSignal",                     1.0    );

     screwdriver.SetGlobalStringOption("Plot", "infoTopRight", "CMS Preliminary");
     screwdriver.SetGlobalStringOption("Plot", "infoTopLeft",  "#sqrt{s} = 8 TeV, L = 19.5 fb^{-1}");

     screwdriver.SetGlobalBoolOption("Plot", "exportPdf", true);
     screwdriver.SetGlobalBoolOption("Plot", "exportEps", false);
     screwdriver.SetGlobalBoolOption("Plot", "exportPng", false);

  // ########################################
  // ##       Run over the datasets        ##
  // ########################################

  vector<string> datasetsList;
  screwdriver.GetDatasetList(&datasetsList);

  cout << "   > Reading datasets... " << endl;
  cout << endl;

    Table scaleFactors = Table("../prediction/scaleFactors/preselection.tab");

    Figure SF_pre           = scaleFactors.Get("value","SF_pre");
    Figure SF_post          = scaleFactors.Get("value","SF_post");
    Figure SF_0btag         = scaleFactors.Get("value","SF_0btag");
    Figure SF_vetopeak      = scaleFactors.Get("value","SF_vetopeak");
    Figure SF_MTtail_1ltop  = scaleFactors.Get("value","SF_MTtail_1ltop");
    Figure SF_MTtail_Wjets  = scaleFactors.Get("value","SF_MTtail_Wjets");

  for (unsigned int d = 0 ; d < datasetsList.size() ; d++)
  {
     string currentDataset = datasetsList[d];
     string currentProcessClass = screwdriver.GetProcessClass(currentDataset);

     sampleName = currentDataset;
     sampleType = screwdriver.GetProcessClassType(currentProcessClass);

     // Open the tree
     TFile f((string(FOLDER_BABYTUPLES)+currentDataset+".root").c_str());
     TTree* theTree = (TTree*) f.Get("babyTuple");

     intermediatePointers pointers;
     InitializeBranchesForReading(theTree,&myEvent,&pointers);

     myEvent.mStop       = -1;
     myEvent.mNeutralino = -1;

  // ########################################
  // ##        Run over the events         ##
  // ########################################

     bool ttbarDatasetToBeSplitted = false;
      if (findSubstring(currentDataset,"ttbar")
      && (currentDataset != "ttbar_madgraph_1l")
      && (currentDataset != "ttbar_madgraph_2l"))
          ttbarDatasetToBeSplitted = true;

      int nEntries = theTree->GetEntries();
      for (int i = 0 ; i < nEntries ; i++)
      {
          if (i % (nEntries / 50) == 0) printProgressBar(i,nEntries,currentDataset);

          // Get the i-th entry
          ReadEvent(theTree,i,&pointers,&myEvent);

          float weight = getWeight();

          // Split 1-lepton ttbar and 2-lepton ttbar
          string currentProcessClass_ = currentProcessClass;
          if (ttbarDatasetToBeSplitted && (myEvent.numberOfGenLepton == 2))
              currentProcessClass_ = "ttbar_2l";

          if (myEvent.MT > 100)
          {
                   if (currentProcessClass_ == "1ltop" ) weight *= SF_MTtail_1ltop.value();
              else if (currentProcessClass_ == "W+jets") weight *= SF_MTtail_Wjets.value();
          }

          screwdriver.AutoFillProcessClass(currentProcessClass_,weight);
      }

      printProgressBar(nEntries,nEntries,currentDataset);
      cout << endl;
      f.Close();

  }

    // #############################
    // ##   Apply scale factors   ##
    // #############################

    screwdriver.ApplyScaleFactor("W+jets",   "0btag_noMTCut",        "singleLepton", SF_0btag);
    screwdriver.ApplyScaleFactor("1ltop",    "0btag_noMTCut",        "singleLepton", SF_0btag);

  // ###################################
  // ##   Make plots and write them   ##
  // ###################################

  cout << endl;
  cout << "   > Making plots..." << endl;
  screwdriver.MakePlots();
  cout << "   > Saving plots..." << endl;
  screwdriver.WritePlots("./plots/fullMTwithSFRinTail/");

  printBoxedMessage("Plot generation completed");

  // #############################
  // ##   Post-plotting tests   ##
  // #############################

  printBoxedMessage("Program done.");
  return (0);
}
コード例 #4
0
backgroundEstimation::backgroundEstimation(string signalRegionLabel_)
{
    signalRegionLabel = signalRegionLabel_;

    // ########################
    // #  Read raw MC yields  #
    // ########################

    rawYieldTable = Table("rawYieldTables/"+signalRegionLabel+".tab");
    secondLeptonInAcceptanceYieldTable = Table("secondLeptonInAcceptance/"+signalRegionLabel+".tab");

    // ###################################################################
    // #  Initialize the table containing the raw mc and the prediction  #
    // ###################################################################

    vector<string> predictionTableColumns = { "raw_mc", "prediction" };
    predictionTable = Table(predictionTableColumns, processesTagList);

    // Raw MC
    Figure N1ltop_mc  = rawYieldTable.Get("signalRegion_MTtail","1ltop"   );
    Figure N2ltop_mc  = rawYieldTable.Get("signalRegion_MTtail","ttbar_2l");
    Figure Nwjets_mc  = rawYieldTable.Get("signalRegion_MTtail","W+jets"  );
    Figure Nrare_mc   = rawYieldTable.Get("signalRegion_MTtail","rare"    );
    Figure NSumBkg_mc = N1ltop_mc+N2ltop_mc+Nwjets_mc+Nrare_mc;

    predictionTable.Set("raw_mc","1ltop",    N1ltop_mc );
    predictionTable.Set("raw_mc","ttbar_2l", N2ltop_mc );
    predictionTable.Set("raw_mc","W+jets",   Nwjets_mc );
    predictionTable.Set("raw_mc","rare",     Nrare_mc  );
    predictionTable.Set("raw_mc","totalSM",  NSumBkg_mc);

    // ######################################################
    // #  Initialize the table containing the scale factors #
    // ######################################################

    vector<string> dummy = { "value" };
    scaleFactorTable = Table(dummy,scaleFactorsTagList);


    // ####################################
    // #  Read results from other modules #
    // ####################################

    #ifdef USING_CR4_AND_CR5_CHECKS_FROM_LOOSENED_SIGNAL_REGIONS
    	// Read the uncertainty from CR4 CR5 modeling checks
        Table temporaryTable;

        if (findSubstring(signalRegionLabel,"BDT"))
    	    temporaryTable = Table("../checkCR4CR5modeling/results/uncertainties_BDT.tab");
        else
    	    temporaryTable = Table("../checkCR4CR5modeling/results/uncertainties_CnC.tab");

        ttll_CR4CR5ModelingUncertainty = temporaryTable.Get("value",signalRegionLabel).error();
    #else
    	ttll_CR4CR5ModelingUncertainty = PLACEHOLDER_CR4_CR5_MODELING_UNCERTAINTY;
    #endif

    ttll_nJetsModelingUncertainty = DILEPTON_TTBAR_JET_MULTIPLICITY_MODELING_UNCERTAINTY;

    // ##################################################################
    // #  Initialize the table containing the systematic uncertainties  #
    // ##################################################################

    vector<string> uncertainties = { "absolute", "relative" };
    systematicsUncertainties = Table(uncertainties, systematicsTagList);

    // Reset the systematics flags
    ResetSystematics();
}