void ssdl2015card_test() {
md.refresh();
//general parameters ********************* general parameters
string dir="SSDL2015";
string fileName="ssdl2fb_oct"; // not needed for statistics but provides the normalization
string fileList="ssdl2fb_oct"; // put command line that gives all files as in a "ls" command
bool mcOnly = false;
md.anConf.configureNames( dir, fileName, fileList );
md.anConf.configureData(false, 0, mcOnly);
//Lumis( or XSections ) pb-1 & KFactors ************************************
float lumi=1; //pb-1 19470
float energy=13; //TeV
bool useXS=false;
//md.anConf.loadXSDB("XSectionsSpring15.db");
map<string,float> LumisXS;
//via XSect
map<string,float> KFactors;
md.anConf.configureLumi( LumisXS, KFactors, lumi, useXS );
//===============================================================
md.addDataCardSigSample("HHSR1:_standard_prompt_ttW","ttW");
md.addDataCardSample("HHSR1:_standard_prompt_ttZH","rares");
md.addDataCardSample("HHSR1:_standard_prompt_Rares","rares");
md.addDataCardSample("HHSR1:_standard_prompt_WW", "WW");
md.addDataCardSample("HHSR1:_standard_prompt_WZ", "WZ");
md.addDataCardSample("HHSR1:_standard_prompt_TG", "TG");
md.addDataCardSample("HHSR1:_standard_prompt_VG", "VG");
md.addDataCardSample("HHSR1:_flipsappl_data", "flip");
md.addDataCardSample("HHSR1:_fakesappl_data", "fake");
md.addDataCardSample("HHSR1:_standard_prompt_ttZH","pseudodata");
md.addDataCardSample("HHSR1:_standard_prompt_Rares","pseudodata");
md.addDataCardSample("HHSR1:_standard_prompt_WW", "pseudodata");
md.addDataCardSample("HHSR1:_standard_prompt_WZ", "pseudodata");
md.addDataCardSample("HHSR1:_standard_prompt_TG", "pseudodata");
md.addDataCardSample("HHSR1:_standard_prompt_VG", "pseudodata");
md.addDataCardSample("HHSR1:_flipsappl_data", "pseudodata");
md.addDataCardSample("HHSR1:_fakesappl_data", "pseudodata");
//experimental uncertainties
md.addNuisanceParameter("jes","ttW:rares","lnN","1.08:1.07");
//Data-driven methods
md.addNuisanceParameter("chMID","flip","lnN","1.30");
md.addNuisanceParameter("fratio","fake","lnN","1.35");
//syst from stat uncertainties == automatic
//md.addNuisanceParameter("wzstat","wz","lnN",""); //--> already given by the 8/15% on Xsect
// md.addNuisanceParameter("tthstat","tth","lnN","");
// md.addNuisanceParameter("ttwstat","ttw","lnN","");
// md.addNuisanceParameter("ttzstat","ttz","lnN","");
// md.addNuisanceParameter("fakestat","fake","lnN","");
// md.addNuisanceParameter("chMIdstat","chMId","lnN","");
// md.addNuisanceParameter("SIGTAGstat","SIGTAG","lnN","");
// string srs[66]={
// "SR1A", "SR2A", "SR3A", "SR4A", "SR5A", "SR6A", "SR7A", "SR8A",
// "SR9A", "SR10A", "SR11A", "SR12A", "SR13A", "SR14A", "SR15A", "SR16A",
// "SR17A", "SR18A", "SR19A", "SR20A", "SR21A", "SR22A", "SR23A", "SR24A",
// "SR25A", "SR26A", "SR27A", "SR28A", "SR29A", "SR30A", "SR31A", "SR32A",
// "SR1B", "SR2B", "SR3B", "SR4B", "SR5B", "SR6B", "SR7B", "SR8B",
// "SR9B", "SR10B", "SR11B", "SR12B", "SR13B", "SR14B", "SR15B", "SR16B",
// "SR17B", "SR18B", "SR19B", "SR20B", "SR21B", "SR22B", "SR23B", "SR24B",
// "SR25B","SR26B",
// "SR1C", "SR2C", "SR3C", "SR4C", "SR5C", "SR6C", "SR7C", "SR8C"
// };
// for(int isr=0;isr<66;isr++) {
// if(srs[isr]=="CATTAG") {
// md.addNuisanceParameter("tth"+srs[isr]+"stat","tth","lnN","");
// md.addNuisanceParameter("ttw"+srs[isr]+"stat","ttw","lnN","");
// md.addNuisanceParameter("ttz"+srs[isr]+"stat","ttz","lnN","");
// md.addNuisanceParameter("fake"+srs[isr]+"stat","fake","lnN","");
// md.addNuisanceParameter("chMId"+srs[isr]+"stat","chMId","lnN","");
// md.addNuisanceParameter("SIGTAG"+srs[isr]+"stat","SIGTAG","lnN","");
// }
// else {
// md.addNuisanceParameter("tth"+srs[isr]+"stat","tth","lnN","-");
// md.addNuisanceParameter("ttw"+srs[isr]+"stat","ttw","lnN","-");
// md.addNuisanceParameter("ttz"+srs[isr]+"stat","ttz","lnN","-");
// md.addNuisanceParameter("fake"+srs[isr]+"stat","fake","lnN","-");
// md.addNuisanceParameter("chMId"+srs[isr]+"stat","chMId","lnN","-");
// md.addNuisanceParameter("SIGTAG"+srs[isr]+"stat","SIGTAG","lnN","-");
// }
// }//signal regions
// md.addNuisanceParameter("SIGTAGStat","SIGTAG","lnN","");
// md.addNuisanceParameter("ttwStat","ttw","lnN","");
// md.addNuisanceParameter("ttzStat","ttz","lnN","");
// md.addNuisanceParameter("wzStat","wz","lnN","");
//md.addNuisanceParameter("sigNorm","T1t415","lnN","1.10");
// }
//===============================================================
//*********************************************************************²
//Execution macro ******************************************************
//plotting ================
md.dp.setLumiAndEnergy( lumi, energy );
md.prepareDisplay();
// int nCateg=66; //47
// vector<string> _categs(nCateg,"");
// string srs[66]={
// "SR1A", "SR2A", "SR3A", "SR4A", "SR5A", "SR6A", "SR7A", "SR8A",
// "SR9A", "SR10A", "SR11A", "SR12A", "SR13A", "SR14A", "SR15A", "SR16A",
// "SR17A", "SR18A", "SR19A", "SR20A", "SR21A", "SR22A", "SR23A", "SR24A",
// "SR25A", "SR26A", "SR27A", "SR28A", "SR29A", "SR30A", "SR31A", "SR32A",
// "SR1B", "SR2B", "SR3B", "SR4B", "SR5B", "SR6B", "SR7B", "SR8B",
// "SR9B", "SR10B", "SR11B", "SR12B", "SR13B", "SR14B", "SR15B", "SR16B",
// "SR17B", "SR18B", "SR19B", "SR20B", "SR21B", "SR22B", "SR23B", "SR24B",
// "SR25B","SR26B",
// "SR1C", "SR2C", "SR3C", "SR4C", "SR5C", "SR6C", "SR7C", "SR8C"
// }
// _categs.assign(srs, srs+nCateg);
// for(size_t ic=0;ic< _categs.size();ic++) {
// string cat=_categs[ic];
md.makeSingleDataCard("ttW", "global_HHSR1", "selected", "ssdl2015_ttW");
// }
gROOT->ProcessLine(".q");
}
void dataCardProd_HHSR9_T14t12() {
md.refresh();
//general parameters ********************* general parameters
string dir="SSDL2015";
string fileName="ssdl2fb_oct"; // not needed for statistics but provides the normalization
string fileList="ssdl2fb_oct"; // put command line that gives all files as in a "ls" command
bool mcOnly = false;
md.anConf.configureNames( dir, fileName, fileList );
md.anConf.configureData(false, 0, mcOnly);
//Lumis( or XSections ) pb-1 & KFactors ************************************
float lumi=1; //pb-1 19470
float energy=13; //TeV
bool useXS=false;
//md.anConf.loadXSDB("XSectionsSpring15.db");
map<string,float> LumisXS;
//via XSect
map<string,float> KFactors;
md.anConf.configureLumi( LumisXS, KFactors, lumi, useXS );
//===============================================================
md.addDataCardSigSample("HHSR9:_sig_T1t_1200_800","T14t12");
md.addDataCardSample("HHSR9:_standard_prompt_WZ", "WZ");
md.addDataCardSample("HHSR9:_standard_prompt_ttW","ttW");
md.addDataCardSample("HHSR9:_standard_prompt_ttZH","ttHZ");
md.addDataCardSample("HHSR9:_standard_prompt_WW", "WW");
md.addDataCardSample("HHSR9:_standard_prompt_TG", "XG");
md.addDataCardSample("HHSR9:_standard_prompt_VG", "XG");
md.addDataCardSample("HHSR9:_standard_prompt_Rares","rares");
md.addDataCardSample("HHSR9:_flipsappl_data", "flip");
md.addDataCardSample("HHSR9:_fakesappl_data", "fake");
md.addDataCardSample("HHSR9:data","data");
//lumi
md.addNuisanceParameter("lumi","WZ:ttW:ttHZ:WW:XG:rares:T14t12","lnN","1.08:1.08:1.08:1.08:1.08:1.08:1.08");
//experimental uncertainties
md.addNuisanceParameter("jes","ttW:ttZH:WW:XG:rares:T14t12","lnN","1.00:1.01:1.01:1.03:1.03:1.10");
md.addNuisanceParameter("btag","ttW:ttZH:WW:XG:rares:T14t12","lnN","1.11:1.11:1.11:1.09:1.11:1.01");
md.addNuisanceParameter("lEff","ttW:ttZH:WW:XG:rares:T14t12","lnN","1.04:1.04:1.04:1.04:1.04:1.04");
md.addNuisanceParameter("tEff","ttW:ttZH:WW:XG:rares:T14t12","lnN","1.02:1.02:1.02:1.02:1.02:1.02");
md.addNuisanceParameter("tHTE","ttW:ttZH:WW:XG:rares:T14t12","lnN","-:-:-:-:-:-");
//test
md.addNuisanceParameter("BTag","ttW","lnN","");
//Data-driven methods
md.addNuisanceParameter("chMID","flip","lnN","1.30");
md.addNuisanceParameter("fratio","fake","lnN","1.30");
md.addNuisanceParameter("frEwk","fake","lnN","1.10");
//theoretical uncertainties
md.addNuisanceParameter("wzTh","WZ","lnN","1.30");
md.addNuisanceParameter("WWTh","ttW","lnN","1.50");
md.addNuisanceParameter("TGTh","XG","lnN","1.50");
md.addNuisanceParameter("rareTh","rares","lnN","1.50");
//md.addNuisanceParameter("sigTh","T14t12","lnN","SIGTHTAG");
md.addNuisanceParameter("ttWPdf","ttW","lnN","1.04");
md.addNuisanceParameter("ttHZPdf","ttHZ","lnN","1.04");
md.addNuisanceParameter("ttWXs","ttW","lnN","1.13");
md.addNuisanceParameter("ttHZXs","ttHZ","lnN","1.11");
md.addNuisanceParameter("ttWAcc","ttW","lnN","1.03");
md.addNuisanceParameter("ttHZAcc","ttHZ","lnN","1.03");
string srs[66]={
"HHSR1", "HHSR2", "HHSR3", "HHSR4", "HHSR5", "HHSR6", "HHSR7", "HHSR8",
"HHSR9", "HHSR10", "HHSR11", "HHSR12", "HHSR13", "HHSR14", "HHSR15", "HHSR16",
"HHSR17", "HHSR18", "HHSR19", "HHSR20", "HHSR21", "HHSR22", "HHSR23", "HHSR24",
"HHSR25", "HHSR26", "HHSR27", "HHSR28", "HHSR29", "HHSR30", "HHSR31", "HHSR32",
"HLSR1", "HLSR2", "HLSR3", "HLSR4", "HLSR5", "HLSR6", "HLSR7", "HLSR8",
"HLSR9", "HLSR10", "HLSR11", "HLSR12", "HLSR13", "HLSR14", "HLSR15", "HLSR16",
"HLSR17", "HLSR18", "HLSR19", "HLSR20", "HLSR21", "HLSR22", "HLSR23", "HLSR24",
"HLSR25","HLSR26",
"LLSR1", "LLSR2", "LLSR3", "LLSR4", "LLSR5", "LLSR6", "LLSR7", "LLSR8"
};
//statistical uncertainties
for(int isr=0;isr<66;isr++) {
if(srs[isr]=="HHSR9") {
md.addNuisanceParameter("ttW"+srs[isr]+"stat","ttW","gmN","");
md.addNuisanceParameter("ttHZ"+srs[isr]+"stat","ttHZ","lnN","");
md.addNuisanceParameter("WW"+srs[isr]+"stat","WW","lnN","");
md.addNuisanceParameter("XG"+srs[isr]+"stat","XG","lnN","1.18");
md.addNuisanceParameter("rares"+srs[isr]+"stat","rares","lnN","");
md.addNuisanceParameter("fake"+srs[isr]+"stat","fake","lnN","1.07");
md.addNuisanceParameter("flip"+srs[isr]+"stat","flip","lnN","");
md.addNuisanceParameter("T14t12"+srs[isr]+"stat","T14t12","lnN","");
}
else {
md.addNuisanceParameter("ttW"+srs[isr]+"stat","ttW","gmN","-");
md.addNuisanceParameter("ttHZ"+srs[isr]+"stat","ttHZ","lnN","-");
md.addNuisanceParameter("WW"+srs[isr]+"stat","WW","lnN","-");
md.addNuisanceParameter("XG"+srs[isr]+"stat","XG","lnN","-");
md.addNuisanceParameter("rares"+srs[isr]+"stat","rares","lnN","-");
md.addNuisanceParameter("fake"+srs[isr]+"stat","fake","lnN","-");
md.addNuisanceParameter("flip"+srs[isr]+"stat","flip","lnN","-");
md.addNuisanceParameter("T14t12"+srs[isr]+"stat","T14t12","lnN","-");
}
}//signal regions
//md.overwriteNuisanceParameter("ttWHHSR9stat","ttW","2.00");
//===============================================================
//*********************************************************************²
//Execution macro ******************************************************
//plotting ================
md.dp.setLumiAndEnergy( lumi, energy );
md.prepareDisplay();
md.makeSingleDataCard("T14t12", "global_HHSR9", "selected", "ssdl2015_HHSR9_T14t12");
gROOT->ProcessLine(".q");
}
void susy3l_datacard() {
md.refresh();
//general parameters ********************* general parameters
string dir="SUSY3L";
string fileName="limits_160126_OnAndOffZ"; // not needed for statistics but provides the normalization
string fileList="limits_160126_OnAndOffZ"; // put command line that gives all files as in a "ls" command
bool mcOnly = false;
md.anConf.configureNames( dir, fileName, fileList );
md.anConf.configureData(false, 0, mcOnly);
//Lumis( or XSections ) pb-1 & KFactors ************************************
float lumi=2260; //pb-1 19470
float energy=13; //TeV
bool useXS=false;
md.anConf.loadXSDB("XSectionsSpring15.db");
map<string,float> LumisXS;
md.anConf.loadKFDB("kFactorsSpring15.db");
//via XSect
map<string,float> KFactors;
md.anConf.configureLumi( LumisXS, KFactors, lumi, useXS );
//===============================================================
string signal = "T1t4_1200"; md.addDataCardSigSample("T1tttt_mGo1200_mChi800", signal);
//string signal = "T1t4_1500"; md.addDataCardSigSample("T1tttt_mGo1500_mChi100", signal);
//string signal = "T54q_deg"; md.addDataCardSigSample("T5ttttDeg_mGo1000_mStop300_mCh285_mChi280", signal);
//string signal = "T6t2W2_600"; md.addDataCardSigSample("T6ttWW_mSbot600_mCh425_mChi50", signal);
//string signal = "T6t2W2_650"; md.addDataCardSigSample("T6ttWW_mSbot650_mCh150_mChi50", signal);
//string signal = "T54qWZ_315"; md.addDataCardSigSample("T5qqqqWZDeg_mGo1000_mCh315_mChi300_dilep", signal);
//string signal = "T54qWZ_325"; md.addDataCardSigSample("T5qqqqWZDeg_mGo1000_mCh325_mChi300_dilep", signal);
//string signal = "T54qWZ_800"; md.addDataCardSigSample("T5qqqqWZ_mGo1500_mCh800_mChi100_lep", signal);
//string signal = "T54qWZ_1000"; md.addDataCardSigSample("T5qqqqWZ_mGo1200_mCh1000_mChi800_lep", signal);
md.addDataCardSample("GGHZZ4L", "rare");
md.addDataCardSample("VHToNonbb", "rare");
md.addDataCardSample("ZZTo4L", "rare");
md.addDataCardSample("WWZ", "rare");
md.addDataCardSample("WZZ","rare");
md.addDataCardSample("ZZZ", "rare");
md.addDataCardSample("TTTT", "rare");
md.addDataCardSample("tZq_ll", "rare");
md.addDataCardSample("WZTo3LNu","WZ");
md.addDataCardSample("TGJets", "XG");
md.addDataCardSample("TTGJets", "XG");
md.addDataCardSample("WGToLNuG", "XG");
md.addDataCardSample("ZGTo2LG", "XG");
md.addDataCardSample("TTZToLLNuNu","ttZH");
md.addDataCardSample("TTHnobb", "ttZH");
md.addDataCardSample("TTLLJets_m1to10", "ttZH");
md.addDataCardSample("TTWToLNu", "ttW");
md.addDataCardSample("TT_pow", "fake");
md.addDataCardSample("DYJetsToLL_M50", "fake");
md.addDataCardSample("DYJetsToLL_M10to50", "fake");
md.addDataCardSample("TbarToLeptons_tch", "fake");
md.addDataCardSample("TBar_tWch", "fake");
md.addDataCardSample("TToLeptons_sch_amcatnlo", "fake");
md.addDataCardSample("TToLeptons_tch", "fake");
md.addDataCardSample("T_tWch", "fake");
md.addDataCardSample("WJetsToLNu", "fake");
//md.addDataCardSample("data","data");
//shape uncertainties
string nuis = "rare:WZ:XG:ttZH:ttW:fake:";
nuis += signal;
md.addNuisanceParameter("JES",nuis,"shape","");
/*md.addNuisanceParameter("BTAG","ttW:ttZH:WZ:WW:XG:rares:T1ttttBENCH","shape","");
md.addNuisanceParameter("tHTE","ttW:ttZH:WZ:WW:XG:rares:T1ttttBENCH","shape","");
md.addNuisanceParameter("Ewk","fake","shape","");
md.addNuisanceParameter("ttWAcc","ttW","shape","");
md.addNuisanceParameter("ttHZAcc","ttHZ","shape","");
md.addNuisanceParameter("BTAGFS","T1ttttBENCH","shape","");
md.addNuisanceParameter("LepEffFS","T1ttttBENCH","shape","");
md.addNuisanceParameter("ISR","T1ttttBENCH","shape","");
md.addNuisanceParameter("QCDScale","T1ttttBENCH","shape","");
*/
//vector<string> cats(
string cats[30]={
"OnZSR001","OnZSR002", "OnZSR003", "OnZSR004", "OnZSR005", "OnZSR006", "OnZSR007", "OnZSR008", "OnZSR009", "OnZSR010", "OnZSR011", "OnZSR012", "OnZSR013", "OnZSR014", "OnZSR015",
"OffZSR001","OffZSR002", "OffZSR003", "OffZSR004", "OffZSR005", "OffZSR006", "OffZSR007", "OffZSR008", "OffZSR009", "OffZSR010", "OffZSR011", "OffZSR012", "OffZSR013", "OffZSR014", "OffZSR015" };
//vector<string> dss;
string dss[7]={"rare","WZ","XG","ttZH","ttW","fake",signal};
//for(int isr=0;isr<66/*cats.size()*/;isr++) {
// md.addNuisanceParameter("pdfS"+cats[isr],"T1ttttBENCH","shape","");
//md.addNuisanceParameter("flip"+cats[isr]+"stat","flip","shape","");
//md.addNuisanceParameter("fake"+cats[isr]+"stat","fake","shape","");
//for(size_t id=0;id<7/*dss.size()*/;id++) {
// md.addNuisanceParameter(dss[id]+cats[isr]+"stat",dss[id],"shape","");
//}
//}
/*
//Flat uncertanties =================================
//lumi
md.addNuisanceParameter("lumi","ttW:ttHZ:WW:XG:rares:T1ttttBENCH","lnN","1.046:1.046:1.046:1.046:1.046:1.046:1.046");
//md.addNuisanceParameter("lumi","ttW:ttHZ:WW:XG:rares:T1ttttBENCH","lnN","1.12:1.12:1.12:1.12:1.12:1.12:1.12");
//experimental uncertainties
md.addNuisanceParameter("lEff","ttW:ttHZ:WW:XG:rares:T1ttttBENCH","lnN","1.04:1.04:1.04:1.04:1.04:1.04");
md.addNuisanceParameter("tEff","ttW:ttHZ:WW:XG:rares:T1ttttBENCH","lnN","1.02:1.02:1.02:1.02:1.02:1.02");
//Data-driven methods
md.addNuisanceParameter("chMID","flip","lnN","1.30");
md.addNuisanceParameter("fratio","fake","lnN","1.30");
//theoretical uncertainties
md.addNuisanceParameter("wzTh","WZ","lnN","1.30");
md.addNuisanceParameter("WWTh","ttW","lnN","1.50");
md.addNuisanceParameter("TGTh","XG","lnN","1.50");
md.addNuisanceParameter("rareTh","rares","lnN","1.50");
md.addNuisanceParameter("ttWPdf","ttW","lnN","1.04");
md.addNuisanceParameter("ttHZPdf","ttHZ","lnN","1.04");
md.addNuisanceParameter("ttWXs","ttW","lnN","1.13");
md.addNuisanceParameter("ttHZXs","ttHZ","lnN","1.11");
md.addNuisanceParameter("fsHLT","T1ttttBENCH","lnN","1.05");
*/
md.addNuisanceParameter("fakeUnc","fake","lnN","1.30");
md.addNuisanceParameter("sigUnc",signal,"lnN","1.10");
md.addNuisanceParameter("promptUnc","rare:WZ:XG:ttZH:ttW","lnN","1.20:1.20:1.20:1.20:1.20");
//===============================================================
//*********************************************************************²
//Execution macro ******************************************************
//plotting ================
md.dp.setLumiAndEnergy( lumi, energy );
md.prepareDisplay();
int ncategs=30;
// vector<string> categs(
string categs[30]={"global_OnZSR001","global_OnZSR002","global_OnZSR003","global_OnZSR004",
"global_OnZSR005","global_OnZSR006","global_OnZSR007","global_OnZSR008","global_OnZSR009",
"global_OnZSR010","global_OnZSR011","global_OnZSR012","global_OnZSR013","global_OnZSR014","global_OnZSR015",
"global_OffZSR001","global_OffZSR002","global_OffZSR003","global_OffZSR004",
"global_OffZSR005","global_OffZSR006","global_OffZSR007","global_OffZSR008","global_OffZSR009",
"global_OffZSR010","global_OffZSR011","global_OffZSR012","global_OffZSR013","global_OffZSR014","global_OffZSR015"
};
vector<string> vcategs;
for(int i=0;i<30;i++) {
vcategs.push_back( categs[i] );
}
/*
//external uncertainties ===================================
//vector<float> AccTTW
float AccTTW[66]={
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.03,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.03,0.08,
0.03,0.08,0.08,0.08,0.03,0.08,0.03,0.08,
0.08,0.08,0.03,0.08,0.03,0.08,0.03,0.08,
0.08,0.08,
0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,
};
//vector<float> AccTTHZ
float AccTTHZ[66]={
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.03,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.03,0.08,
0.03,0.08,0.08,0.08,0.03,0.08,0.03,0.08,
0.08,0.08,0.03,0.08,0.03,0.08,0.03,0.08,
0.08,0.08,
0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,
};
// vector<float> EWKFake({
// 0.10,0.06,0.10,0.03,0.09,0.10,0.05,0.10,
// 0.09,0.11,0.11,0.09,0.09,0.02,0.10,0.10,
// 0.20,0.03,0.10,0.05,0.05,0.10,0.10,0.10,
// 0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,
// 0.03,0.01,0.10,0.02,0.01,0.10,0.02,0.01,
// 0.01,0.01,0.01,0.01,0.02,0.01,0.10,0.10,
// 0.01,0.10,0.10,0.10,0.10,0.10,0.10,0.10,
// 0.10,0.10,
// 0.10,0.01,0.10,0.10,0.10,0.10,0.10,0.10,
// });
//vector<float> HLTEff
float HLTEff[66]={
0,0.02,0,0.02,0.02,0.02,0.02,0.02,
0,0.02,0,0.02,0.02,0.02,0.02,0.02,
0,0.02,0,0.02,0.02,0.02,0.02,0.02,
0,0.02,0,0.02,0,0.02,0.02,0.02,
0,0.02,0,0.02,0.02,0.02,0,0.02,
0,0.02,0.02,0.02,0,0.02,0,0.02,
0.02,0.02,0,0.02,0,0.02,0,0.02,
0.02,0.02,
0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02
};
// int nCateg=66; //47
// vector<string> _categs(nCateg,"");
//qcd scale for signal
TFile* iQCDFile=new TFile("/shome/mmarionn/MPAF/workdir/QCDScales/scaleUncMASS-.root","read");
TH1F* hQCDUp=(TH1F*)iQCDFile->Get("MASS-ScaleUp");
TH1F* hQCDDown=(TH1F*)iQCDFile->Get("MASS-ScaleDown");
*/
//for(size_t ic=0;ic<66/*categs.size()*/;ic++) {
// for(int isr=0;isr<66/*cats.size()*/;isr++) {
// if(categs[ic].substr(7, categs[ic].size()-7)==cats[isr])
// md.addExternalSystUnc("T1ttttBENCH","pdfS"+cats[isr],0.1, -0.1, categs[ic], "selected");
// else
//md.addExternalSystUnc("T1ttttBENCH","pdfS"+cats[isr],0., 0., categs[ic], "selected");
//}
//md.addExternalSystUnc("T1ttttBENCH","QCDScale",hQCDUp->GetBinContent(ic+1), hQCDDown->GetBinContent(ic+1), categs[ic], "selected");
//md.addExternalSystUnc("fake","frEwk",EWKFake[ic], -1*EWKFake[ic], categs[ic], "selected");
//md.addExternalSystUnc("ttW","ttWAcc",AccTTW[ic], -1*AccTTW[ic], categs[ic], "selected");
//md.addExternalSystUnc("ttHZ","ttHZAcc",AccTTHZ[ic], -1*AccTTHZ[ic], categs[ic], "selected");
//for(size_t id=0;id<7/*dss.size()*/;id++) {
// md.addExternalSystUnc(dss[id],"tHTE",HLTEff[ic], -1*HLTEff[ic], categs[ic], "selected");
//}
//}
for(size_t ic=0;ic<30;ic++) {
for(size_t id=0;id<7;id++) {
md.addExternalSystUnc(dss[id],"JES",0.05,-0.05,vcategs[ic],"selected");
}
}
md.makeMultiDataCard(signal, vcategs, "selected", signal);
md.getStatistics("global");
gROOT->ProcessLine(".q");
}
void ssdlScan_BENCH() {
md.refresh();
//general parameters ********************* general parameters
string dir="SSDL2015";
string fileName="merged_2fb"; // not needed for statistics but provides the normalization
string fileList="merged_2fb_T1tttt-MASS-"; // put command line that gives all files as in a "ls" command
bool mcOnly = false;
md.anConf.configureNames( dir, fileName, fileList );
md.anConf.configureData(false, 0, mcOnly);
//Lumis( or XSections ) pb-1 & KFactors ************************************
float lumi=2160; //pb-1 19470
float energy=13; //TeV
bool useXS=true;
md.anConf.loadXSDB("XSectionsSpring15.db");
map<string,float> LumisXS;
//via XSect
map<string,float> KFactors;
KFactors[ "_standard_prompt_ttW" ] = 1.0213;
KFactors[ "_standard_prompt_ttZH" ] = 1.0213;
KFactors[ "_standard_prompt_WWss" ] = 1.0213;
KFactors[ "_standard_prompt_XG" ] = 1.0213;
KFactors[ "_standard_prompt_Rares" ] = 1.0213;
KFactors[ "T1tttt-MASS-" ] = 1.0213;
// marco put this line here:
// _standard_prompt_WZ: WZTo3LNu: 1.22*xsec ; FillColor=ROOT.kOrange, FakeRate="susy-multilepton/standard_sel.txt\,susy-multilepton/load_is_2_prompt_loop.txt", Label="WZ (prompt)", NormSystematic=0.3
KFactors[ "_standard_prompt_WZ" ] = 1.22;
md.anConf.configureLumi( LumisXS, KFactors, lumi, useXS );
//===============================================================
md.addDataCardSigSample("T1tttt-MASS-","T1tttt_BENCH");
md.addDataCardSample("WZTo3LNu" , "WZ" );
md.addDataCardSample("TTWToLNu" , "ttW" );
md.addDataCardSample("TTZToLLNuNu" , "ttHZ" );
md.addDataCardSample("TTLLJets_m1to10" , "ttHZ" );
md.addDataCardSample("TTHnobb" , "ttHZ" );
md.addDataCardSample("WpWpJJ" , "WW" );
md.addDataCardSample("TTGJets" , "XG" );
md.addDataCardSample("TGJets" , "XG" );
md.addDataCardSample("WGToLNuG" , "XG" );
md.addDataCardSample("ZGTo2LG" , "XG" );
md.addDataCardSample("WWDouble" , "rares" );
md.addDataCardSample("ZZTo4L" , "rares" );
md.addDataCardSample("tZq_ll" , "rares" );
md.addDataCardSample("GGHZZ4L" , "rares" );
md.addDataCardSample("VHToNonbb" , "rares" );
md.addDataCardSample("WWZ" , "rares" );
md.addDataCardSample("ZZZ" , "rares" );
md.addDataCardSample("WZZ" , "rares" );
md.addDataCardSample("TTTT" , "rares" );
md.addDataCardSample("DoubleEG_Run2015C_Oct05_runs_254231_254914" , "data" );
md.addDataCardSample("DoubleEG_Run2015D_Oct05_runs_256630_258158" , "data" );
md.addDataCardSample("DoubleEG_Run2015D_PromptV4_runs_258159_260627" , "data" );
md.addDataCardSample("DoubleMuon_Run2015C_Oct05_runs_254231_254914" , "data" );
md.addDataCardSample("DoubleMuon_Run2015D_Oct05_runs_256630_258158" , "data" );
md.addDataCardSample("DoubleMuon_Run2015D_PromptV4_runs_258159_260627" , "data" );
md.addDataCardSample("MuonEG_Run2015C_Oct05_runs_254231_254914" , "data" );
md.addDataCardSample("MuonEG_Run2015D_Oct05_runs_256630_258158" , "data" );
md.addDataCardSample("MuonEG_Run2015D_PromptV4_runs_258159_260627" , "data" );
md.addDataCardSample("mId:DYJetsToLL_M10to50" , "flip" );
md.addDataCardSample("mId:DYJetsToLL_M50" , "flip" );
//md.addDataCardSample("mId:TTWToLNu" , "flip", -1);
//md.addDataCardSample("mId:TTZToLLNuNu" , "flip", -1);
//md.addDataCardSample("mId:TTLLJets_m1to10" , "flip", -1);
//md.addDataCardSample("mId:TTHnobb" , "flip", -1);
//md.addDataCardSample("mId:WZTo3LNu" , "flip", -1);
//md.addDataCardSample("mId:WpWpJJ" , "flip", -1);
//md.addDataCardSample("mId:TTGJets" , "flip", -1);
//md.addDataCardSample("mId:TGJets" , "flip", -1);
//md.addDataCardSample("mId:WGToLNuG" , "flip", -1);
//md.addDataCardSample("mId:ZGTo2LG" , "flip", -1);
//md.addDataCardSample("mId:WWDouble" , "flip", -1);
//md.addDataCardSample("mId:ZZTo4L" , "flip", -1);
//md.addDataCardSample("mId:tZq_ll" , "flip", -1);
//md.addDataCardSample("mId:GGHZZ4L" , "flip", -1);
//md.addDataCardSample("mId:VHToNonbb" , "flip", -1);
//md.addDataCardSample("mId:WWZ" , "flip", -1);
//md.addDataCardSample("mId:ZZZ" , "flip", -1);
//md.addDataCardSample("mId:WZZ" , "flip", -1);
//md.addDataCardSample("mId:TTTT" , "flip", -1);
md.addDataCardSample("data:Fake:DoubleEG_Run2015C_Oct05_runs_254231_254914" , "fake" );
md.addDataCardSample("data:Fake:DoubleEG_Run2015D_Oct05_runs_256630_258158" , "fake" );
md.addDataCardSample("data:Fake:DoubleEG_Run2015D_PromptV4_runs_258159_260627" , "fake" );
md.addDataCardSample("data:Fake:DoubleMuon_Run2015C_Oct05_runs_254231_254914" , "fake" );
md.addDataCardSample("data:Fake:DoubleMuon_Run2015D_Oct05_runs_256630_258158" , "fake" );
md.addDataCardSample("data:Fake:DoubleMuon_Run2015D_PromptV4_runs_258159_260627", "fake" );
md.addDataCardSample("data:Fake:MuonEG_Run2015C_Oct05_runs_254231_254914" , "fake" );
md.addDataCardSample("data:Fake:MuonEG_Run2015D_Oct05_runs_256630_258158" , "fake" );
md.addDataCardSample("data:Fake:MuonEG_Run2015D_PromptV4_runs_258159_260627" , "fake" );
//md.addDataCardSample("Fake:DYJetsToLL_M10to50" , "fake", -1);
//md.addDataCardSample("Fake:DYJetsToLL_M50" , "fake", -1);
md.addDataCardSample("Fake:TTWToLNu" , "fake", -1);
md.addDataCardSample("Fake:TTZToLLNuNu" , "fake", -1);
md.addDataCardSample("Fake:TTLLJets_m1to10" , "fake", -1);
md.addDataCardSample("Fake:TTHnobb" , "fake", -1);
md.addDataCardSample("Fake:WZTo3LNu" , "fake", -1);
md.addDataCardSample("Fake:WpWpJJ" , "fake", -1);
md.addDataCardSample("Fake:TTGJets" , "fake", -1);
md.addDataCardSample("Fake:TGJets" , "fake", -1);
md.addDataCardSample("Fake:WGToLNuG" , "fake", -1);
md.addDataCardSample("Fake:ZGTo2LG" , "fake", -1);
md.addDataCardSample("Fake:WWDouble" , "fake", -1);
md.addDataCardSample("Fake:ZZTo4L" , "fake", -1);
md.addDataCardSample("Fake:tZq_ll" , "fake", -1);
md.addDataCardSample("Fake:GGHZZ4L" , "fake", -1);
md.addDataCardSample("Fake:VHToNonbb" , "fake", -1);
md.addDataCardSample("Fake:WWZ" , "fake", -1);
md.addDataCardSample("Fake:ZZZ" , "fake", -1);
md.addDataCardSample("Fake:WZZ" , "fake", -1);
md.addDataCardSample("Fake:TTTT" , "fake", -1);
//shape uncertainties
md.addNuisanceParameter("JES","ttW:ttZH:WZ:WW:XG:rares:T1tttt_BENCH","shape","");
md.addNuisanceParameter("BTAG","ttW:ttZH:WZ:WW:XG:rares:T1tttt_BENCH","shape","");
md.addNuisanceParameter("tHTE","ttW:ttZH:WZ:WW:XG:rares:T1tttt_BENCH","shape","");
md.addNuisanceParameter("EWKFR","fake","shape","");
md.addNuisanceParameter("ttWAcc","ttW","shape","");
md.addNuisanceParameter("ttHZAcc","ttHZ","shape","");
md.addNuisanceParameter("BTAGFS","T1tttt_BENCH","shape","");
md.addNuisanceParameter("LepEffFS","T1tttt_BENCH","shape","");
md.addNuisanceParameter("ISR","T1tttt_BENCH","shape","");
md.addNuisanceParameter("QCDScale","T1tttt_BENCH","shape","");
//vector<string> cats(
string cats[66]={
"HHSR1","HHSR2","HHSR3","HHSR4","HHSR5","HHSR6","HHSR7","HHSR8",
"HHSR9","HHSR10","HHSR11","HHSR12","HHSR13","HHSR14","HHSR15","HHSR16",
"HHSR17","HHSR18","HHSR19","HHSR20","HHSR21","HHSR22","HHSR23","HHSR24",
"HHSR25","HHSR26","HHSR27","HHSR28","HHSR29","HHSR30","HHSR31","HHSR32",
"HLSR1","HLSR2","HLSR3","HLSR4","HLSR5","HLSR6","HLSR7","HLSR8",
"HLSR9","HLSR10","HLSR11","HLSR12","HLSR13","HLSR14","HLSR15","HLSR16",
"HLSR17","HLSR18","HLSR19","HLSR20","HLSR21","HLSR22","HLSR23","HLSR24",
"HLSR25","HLSR26",
"LLSR1","LLSR2","LLSR3","LLSR4","LLSR5","LLSR6","LLSR7","LLSR8"
};
//vector<string> dss;
string dss[7]={"ttW","ttHZ","WZ","WW","XG","rares","T1tttt_BENCH"};
for(int isr=0;isr<66/*cats.size()*/;isr++) {
md.addNuisanceParameter("pdfS"+cats[isr],"T1tttt_BENCH","shape","");
md.addNuisanceParameter("flip"+cats[isr]+"stat","flip","shape","");
md.addNuisanceParameter("fake"+cats[isr]+"stat","fake","shape","");
for(size_t id=0;id<7/*dss.size()*/;id++) {
md.addNuisanceParameter(dss[id]+cats[isr]+"stat",dss[id],"shape","");
}
}
//Flat uncertanties =================================
//lumi
md.addNuisanceParameter("lumi","ttW:ttHZ:WW:XG:rares:T1tttt_BENCH","lnN","1.046:1.046:1.046:1.046:1.046:1.046:1.046");
//md.addNuisanceParameter("lumi","ttW:ttHZ:WW:XG:rares:T1tttt_BENCH","lnN","1.12:1.12:1.12:1.12:1.12:1.12:1.12");
//experimental uncertainties
md.addNuisanceParameter("lEff","ttW:ttHZ:WW:XG:rares:T1tttt_BENCH","lnN","1.04:1.04:1.04:1.04:1.04:1.04");
md.addNuisanceParameter("tEff","ttW:ttHZ:WW:XG:rares:T1tttt_BENCH","lnN","1.02:1.02:1.02:1.02:1.02:1.02");
//Data-driven methods
md.addNuisanceParameter("chMID","flip","lnN","1.30");
md.addNuisanceParameter("fratio","fake","lnN","1.30");
//theoretical uncertainties
md.addNuisanceParameter("wzTh","WZ","lnN","1.30");
md.addNuisanceParameter("WWTh","ttW","lnN","1.50");
md.addNuisanceParameter("TGTh","XG","lnN","1.50");
md.addNuisanceParameter("rareTh","rares","lnN","1.50");
md.addNuisanceParameter("ttWPdf","ttW","lnN","1.04");
md.addNuisanceParameter("ttHZPdf","ttHZ","lnN","1.04");
md.addNuisanceParameter("ttWXs","ttW","lnN","1.13");
md.addNuisanceParameter("ttHZXs","ttHZ","lnN","1.11");
md.addNuisanceParameter("fsHLT","T1tttt_BENCH","lnN","1.05");
//===============================================================
//*********************************************************************²
//Execution macro ******************************************************
//plotting ================
md.dp.setLumiAndEnergy( lumi, energy );
md.prepareDisplay();
int ncategs=66;
// vector<string> categs(
string categs[66]={"global_HHSR1","global_HHSR2","global_HHSR3","global_HHSR4",
"global_HHSR5","global_HHSR6","global_HHSR7","global_HHSR8",
"global_HHSR9","global_HHSR10","global_HHSR11","global_HHSR12",
"global_HHSR13","global_HHSR14","global_HHSR15","global_HHSR16",
"global_HHSR17","global_HHSR18","global_HHSR19","global_HHSR20",
"global_HHSR21","global_HHSR22","global_HHSR23","global_HHSR24",
"global_HHSR25","global_HHSR26","global_HHSR27","global_HHSR28",
"global_HHSR29","global_HHSR30","global_HHSR31","global_HHSR32",
"global_HLSR1","global_HLSR2","global_HLSR3","global_HLSR4",
"global_HLSR5","global_HLSR6","global_HLSR7","global_HLSR8",
"global_HLSR9","global_HLSR10","global_HLSR11","global_HLSR12",
"global_HLSR13","global_HLSR14","global_HLSR15","global_HLSR16",
"global_HLSR17","global_HLSR18","global_HLSR19","global_HLSR20",
"global_HLSR21","global_HLSR22","global_HLSR23","global_HLSR24",
"global_HLSR25","global_HLSR26",
"global_LLSR1","global_LLSR2","global_LLSR3","global_LLSR4",
"global_LLSR5","global_LLSR6","global_LLSR7","global_LLSR8",
};
vector<string> vcategs;
for(int i=0;i<66;i++) {
vcategs.push_back( categs[i] );
}
//external uncertainties ===================================
//vector<float> AccTTW
float AccTTW[66]={
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.03,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.03,0.08,
0.03,0.08,0.08,0.08,0.03,0.08,0.03,0.08,
0.08,0.08,0.03,0.08,0.03,0.08,0.03,0.08,
0.08,0.08,
0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,
};
//vector<float> AccTTHZ
float AccTTHZ[66]={
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.03,0.08,0.08,0.08,
0.03,0.08,0.03,0.08,0.08,0.08,0.03,0.08,
0.03,0.08,0.08,0.08,0.03,0.08,0.03,0.08,
0.08,0.08,0.03,0.08,0.03,0.08,0.03,0.08,
0.08,0.08,
0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,
};
// vector<float> EWKFake({
// 0.10,0.06,0.10,0.03,0.09,0.10,0.05,0.10,
// 0.09,0.11,0.11,0.09,0.09,0.02,0.10,0.10,
// 0.20,0.03,0.10,0.05,0.05,0.10,0.10,0.10,
// 0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,
// 0.03,0.01,0.10,0.02,0.01,0.10,0.02,0.01,
// 0.01,0.01,0.01,0.01,0.02,0.01,0.10,0.10,
// 0.01,0.10,0.10,0.10,0.10,0.10,0.10,0.10,
// 0.10,0.10,
// 0.10,0.01,0.10,0.10,0.10,0.10,0.10,0.10,
// });
//vector<float> HLTEff
float HLTEff[66]={
0,0.02,0,0.02,0.02,0.02,0.02,0.02,
0,0.02,0,0.02,0.02,0.02,0.02,0.02,
0,0.02,0,0.02,0.02,0.02,0.02,0.02,
0,0.02,0,0.02,0,0.02,0.02,0.02,
0,0.02,0,0.02,0.02,0.02,0,0.02,
0,0.02,0.02,0.02,0,0.02,0,0.02,
0.02,0.02,0,0.02,0,0.02,0,0.02,
0.02,0.02,
0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02
};
// int nCateg=66; //47
// vector<string> _categs(nCateg,"");
//qcd scale for signal
TFile* iQCDFile=new TFile("/shome/cheidegg/t/runScanStuff/step5/QCDScales/scaleUnc-MASS-.root","read");
TH1F* hQCDUp=(TH1F*)iQCDFile->Get("-MASS-ScaleUp");
TH1F* hQCDDown=(TH1F*)iQCDFile->Get("-MASS-ScaleDown");
for(size_t ic=0;ic<66/*categs.size()*/;ic++) {
for(int isr=0;isr<66/*cats.size()*/;isr++) {
if(categs[ic].substr(7, categs[ic].size()-7)==cats[isr])
md.addExternalSystUnc("T1tttt_BENCH","pdfS"+cats[isr],0.1, -0.1, categs[ic], "selected");
else
md.addExternalSystUnc("T1tttt_BENCH","pdfS"+cats[isr],0., 0., categs[ic], "selected");
}
md.addExternalSystUnc("T1tttt_BENCH","QCDScale",hQCDUp->GetBinContent(ic+1), hQCDDown->GetBinContent(ic+1), categs[ic], "selected");
//md.addExternalSystUnc("fake","frEwk",EWKFake[ic], -1*EWKFake[ic], categs[ic], "selected");
md.addExternalSystUnc("ttW","ttWAcc",AccTTW[ic], -1*AccTTW[ic], categs[ic], "selected");
md.addExternalSystUnc("ttHZ","ttHZAcc",AccTTHZ[ic], -1*AccTTHZ[ic], categs[ic], "selected");
for(size_t id=0;id<7/*dss.size()*/;id++) {
md.addExternalSystUnc(dss[id],"tHTE",HLTEff[ic], -1*HLTEff[ic], categs[ic], "selected");
}
}
md.makeMultiDataCard("T1tttt_BENCH", vcategs, "selected", "ssdlScan_2fb_T1tttt_BENCH");
gROOT->ProcessLine(".q");
}
void dataCardProd_CATTAG_SIGTAG_LUMTAG() {
md.refresh();
//general parameters ********************* general parameters
string dir="SSDL2015";
string fileName="fakeEstim_SR"; // not needed for statistics but provides the normalization
string fileList="fakeEstim_SR"; // put command line that gives all files as in a "ls" command
bool mcOnly = false;
md.anConf.configureNames( dir, fileName, fileList );
md.anConf.configureData(false, 0, mcOnly);
//Lumis( or XSections ) pb-1 & KFactors ************************************
float lumi=LUMTAG; //pb-1 19470
float energy=13; //TeV
bool useXS=false;
md.anConf.loadXSDB("XSectionsPhys14.db");
map<string,float> LumisXS;
//via XSect
map<string,float> KFactors;
md.anConf.configureLumi( LumisXS, KFactors, lumi, useXS );
//===============================================================
md.addDataCardSigSample("CATTAG:SIGNAME","SIGTAG");
md.addDataCardSample("CATTAG:WZJetsTo3LNu", "wz");
md.addDataCardSample("CATTAG:TTH", "tth");
md.addDataCardSample("CATTAG:TTZJets", "ttz");
md.addDataCardSample("CATTAG:TTWJets", "ttw");
md.addDataCardSample("CATTAG:TTJetscharge", "chMId");
md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT100to200charge", "chMId");
md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT200to400charge", "chMId");
md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT400to600charge", "chMId");
md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT600toInfcharge", "chMId");
md.addDataCardSample("CATTAG_Fake:TTJets", "fake");
md.addDataCardSample("CATTAG_Fake:WJetsToLNu_HT100to200", "fake");
md.addDataCardSample("CATTAG_Fake:WJetsToLNu_HT200to400", "fake");
md.addDataCardSample("CATTAG_Fake:WJetsToLNu_HT400to600", "fake");
md.addDataCardSample("CATTAG_Fake:WJetsToLNu_HT600toInf", "fake");
md.addDataCardSample("CATTAG_Fake:DYJetsToLL_M50_HT100to200", "fake");
md.addDataCardSample("CATTAG_Fake:DYJetsToLL_M50_HT200to400", "fake");
md.addDataCardSample("CATTAG_Fake:DYJetsToLL_M50_HT400to600", "fake");
md.addDataCardSample("CATTAG_Fake:DYJetsToLL_M50_HT600toInf", "fake");
// md.addDataCardSample("CATTAG:TTJetsfake", "fake");
// md.addDataCardSample("CATTAG:WJetsToLNu_HT100to200fake", "fake");
// md.addDataCardSample("CATTAG:WJetsToLNu_HT200to400fake", "fake");
// md.addDataCardSample("CATTAG:WJetsToLNu_HT400to600fake", "fake");
// md.addDataCardSample("CATTAG:WJetsToLNu_HT600toInffake", "fake");
// md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT100to200fake", "fake");
// md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT200to400fake", "fake");
// md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT400to600fake", "fake");
// md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT600toInffake", "fake");
md.addDataCardSample("CATTAG:WZJetsTo3LNu", "pseudodata");
md.addDataCardSample("CATTAG:TTH", "pseudodata");
md.addDataCardSample("CATTAG:TTZJets", "pseudodata");
md.addDataCardSample("CATTAG:TTWJets", "pseudodata");
md.addDataCardSample("CATTAG:TTJetscharge", "pseudodata");
md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT100to200charge", "pseudodata");
md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT200to400charge", "pseudodata");
md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT400to600charge", "pseudodata");
md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT600toInfcharge", "pseudodata");
md.addDataCardSample("CATTAG_Fake:TTJets", "pseudodata");
md.addDataCardSample("CATTAG_Fake:WJetsToLNu_HT100to200", "pseudodata");
md.addDataCardSample("CATTAG_Fake:WJetsToLNu_HT200to400", "pseudodata");
md.addDataCardSample("CATTAG_Fake:WJetsToLNu_HT400to600", "pseudodata");
md.addDataCardSample("CATTAG_Fake:WJetsToLNu_HT600toInf", "pseudodata");
md.addDataCardSample("CATTAG_Fake:DYJetsToLL_M50_HT100to200", "pseudodata");
md.addDataCardSample("CATTAG_Fake:DYJetsToLL_M50_HT200to400", "pseudodata");
md.addDataCardSample("CATTAG_Fake:DYJetsToLL_M50_HT400to600", "pseudodata");
md.addDataCardSample("CATTAG_Fake:DYJetsToLL_M50_HT600toInf", "pseudodata");
// md.addDataCardSample("CATTAG:TTJetsfake", "pseudodata");
// md.addDataCardSample("CATTAG:WJetsToLNu_HT100to200fake", "pseudodata");
// md.addDataCardSample("CATTAG:WJetsToLNu_HT200to400fake", "pseudodata");
// md.addDataCardSample("CATTAG:WJetsToLNu_HT400to600fake", "pseudodata");
// md.addDataCardSample("CATTAG:WJetsToLNu_HT600toInffake", "pseudodata");
// md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT100to200fake", "pseudodata");
// md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT200to400fake", "pseudodata");
// md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT400to600fake", "pseudodata");
// md.addDataCardSample("CATTAG:DYJetsToLL_M50_HT600toInffake", "pseudodata");
//experimental uncertainties
md.addNuisanceParameter("jes","ttw:ttz:SIGTAG","lnN","JESTAG");
//extrapolation uncertainties
md.addNuisanceParameter("ttwTail","ttw","lnN","TTWTAILTAG");
md.addNuisanceParameter("ttzTail","ttz","lnN","TTZTAILTAG");
md.addNuisanceParameter("wzB","wz","lnN","WZBTAG");
//flat uncertainties -> cross sections
// md.addNuisanceParameter("ttwNorm","ttw","lnN","1.1"); //->merged with the "tail" source
// md.addNuisanceParameter("ttwNorm","ttz","lnN","1.1");
md.addNuisanceParameter("wzNorm","wz","lnN","WZLUMITAG");
//Data-driven methods
md.addNuisanceParameter("chMID","chMId","lnN","1.50");
md.addNuisanceParameter("fratio","fake","lnN","1.40");
//syst from stat uncertainties == automatic
//md.addNuisanceParameter("wzstat","wz","lnN",""); //--> already given by the 8/15% on Xsect
// md.addNuisanceParameter("tthstat","tth","lnN","");
// md.addNuisanceParameter("ttwstat","ttw","lnN","");
// md.addNuisanceParameter("ttzstat","ttz","lnN","");
// md.addNuisanceParameter("fakestat","fake","lnN","");
// md.addNuisanceParameter("chMIdstat","chMId","lnN","");
// md.addNuisanceParameter("SIGTAGstat","SIGTAG","lnN","");
string srs[66]={
"SR1A", "SR2A", "SR3A", "SR4A", "SR5A", "SR6A", "SR7A", "SR8A",
"SR9A", "SR10A", "SR11A", "SR12A", "SR13A", "SR14A", "SR15A", "SR16A",
"SR17A", "SR18A", "SR19A", "SR20A", "SR21A", "SR22A", "SR23A", "SR24A",
"SR25A", "SR26A", "SR27A", "SR28A", "SR29A", "SR30A", "SR31A", "SR32A",
"SR1B", "SR2B", "SR3B", "SR4B", "SR5B", "SR6B", "SR7B", "SR8B",
"SR9B", "SR10B", "SR11B", "SR12B", "SR13B", "SR14B", "SR15B", "SR16B",
"SR17B", "SR18B", "SR19B", "SR20B", "SR21B", "SR22B", "SR23B", "SR24B",
"SR25B","SR26B",
"SR1C", "SR2C", "SR3C", "SR4C", "SR5C", "SR6C", "SR7C", "SR8C"
};
for(int isr=0;isr<66;isr++) {
if(srs[isr]=="CATTAG") {
md.addNuisanceParameter("tth"+srs[isr]+"stat","tth","lnN","");
md.addNuisanceParameter("ttw"+srs[isr]+"stat","ttw","lnN","");
md.addNuisanceParameter("ttz"+srs[isr]+"stat","ttz","lnN","");
md.addNuisanceParameter("fake"+srs[isr]+"stat","fake","lnN","");
md.addNuisanceParameter("chMId"+srs[isr]+"stat","chMId","lnN","");
md.addNuisanceParameter("SIGTAG"+srs[isr]+"stat","SIGTAG","lnN","");
}
else {
md.addNuisanceParameter("tth"+srs[isr]+"stat","tth","lnN","-");
md.addNuisanceParameter("ttw"+srs[isr]+"stat","ttw","lnN","-");
md.addNuisanceParameter("ttz"+srs[isr]+"stat","ttz","lnN","-");
md.addNuisanceParameter("fake"+srs[isr]+"stat","fake","lnN","-");
md.addNuisanceParameter("chMId"+srs[isr]+"stat","chMId","lnN","-");
md.addNuisanceParameter("SIGTAG"+srs[isr]+"stat","SIGTAG","lnN","-");
}
}//signal regions
// md.addNuisanceParameter("SIGTAGStat","SIGTAG","lnN","");
// md.addNuisanceParameter("ttwStat","ttw","lnN","");
// md.addNuisanceParameter("ttzStat","ttz","lnN","");
// md.addNuisanceParameter("wzStat","wz","lnN","");
//md.addNuisanceParameter("sigNorm","T1t415","lnN","1.10");
// }
//===============================================================
//*********************************************************************²
//Execution macro ******************************************************
//plotting ================
md.dp.setLumiAndEnergy( lumi, energy );
md.prepareDisplay();
// int nCateg=66; //47
// vector<string> _categs(nCateg,"");
// string srs[66]={
// "SR1A", "SR2A", "SR3A", "SR4A", "SR5A", "SR6A", "SR7A", "SR8A",
// "SR9A", "SR10A", "SR11A", "SR12A", "SR13A", "SR14A", "SR15A", "SR16A",
// "SR17A", "SR18A", "SR19A", "SR20A", "SR21A", "SR22A", "SR23A", "SR24A",
// "SR25A", "SR26A", "SR27A", "SR28A", "SR29A", "SR30A", "SR31A", "SR32A",
// "SR1B", "SR2B", "SR3B", "SR4B", "SR5B", "SR6B", "SR7B", "SR8B",
// "SR9B", "SR10B", "SR11B", "SR12B", "SR13B", "SR14B", "SR15B", "SR16B",
// "SR17B", "SR18B", "SR19B", "SR20B", "SR21B", "SR22B", "SR23B", "SR24B",
// "SR25B","SR26B",
// "SR1C", "SR2C", "SR3C", "SR4C", "SR5C", "SR6C", "SR7C", "SR8C"
// }
// _categs.assign(srs, srs+nCateg);
// for(size_t ic=0;ic< _categs.size();ic++) {
// string cat=_categs[ic];
md.makeSingleDataCard("SIGTAG", "nominal", "selected", "ssdl2015_CATTAG_SIGTAG_LUMTAG");
// }
gROOT->ProcessLine(".q");
}
