//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // // MVA // //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ void MVA(TString signal = "ttDM0001scalar0500", TString systematic = "LepElepTup", bool doMVATrain = false, bool doMVARead = true) { if (!doMVATrain && !doMVARead) return; gInterpreter->ExecuteMacro("../test/PaperStyle.C"); gSystem->mkdir(trainingdir, kTRUE); gSystem->mkdir(applicationdir, kTRUE); _systematic = systematic; _inputdir = "../minitrees/" + _systematic + "/TTDM/"; (TMVA::gConfig().GetIONames()).fWeightFileDir = weightsdir; // Training //---------------------------------------------------------------------------- if (doMVATrain) MVATrain(signal); // Reading //---------------------------------------------------------------------------- if (doMVARead) { MVARead(signal, signal); if(_systematic == "nominal") MVARead(signal, "01_Data"); MVARead(signal, "14_HZ"); MVARead(signal, "10_HWW"); MVARead(signal, "06_WW"); MVARead(signal, "02_WZTo3LNu"); MVARead(signal, "03_ZZ"); MVARead(signal, "11_Wg"); MVARead(signal, "07_ZJets"); MVARead(signal, "09_TTV"); MVARead(signal, "04_TTTo2L2Nu"); MVARead(signal, "05_ST"); if(_systematic == "nominal") MVARead(signal, "00_Fakes"); } }
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // // MVA // //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ void MVA(TString signal = "monoH_2HDM_MZp-600_MA0-400", bool doMVATrain = true, bool doMVARead = false) { if (!doMVATrain && !doMVARead) return; gInterpreter->ExecuteMacro("../test/PaperStyle.C"); gSystem->mkdir(trainingdir, kTRUE); gSystem->mkdir(applicationdir, kTRUE); (TMVA::gConfig().GetIONames()).fWeightFileDir = weightsdir; // Training //---------------------------------------------------------------------------- if (doMVATrain) MVATrain(signal); // Reading //---------------------------------------------------------------------------- if (doMVARead) { MVARead(signal, signal); MVARead(signal, "01_Data"); MVARead(signal, "14_HZ"); MVARead(signal, "10_HWW"); MVARead(signal, "06_WW"); MVARead(signal, "02_WZTo3LNu"); MVARead(signal, "03_ZZ"); MVARead(signal, "11_Wg"); MVARead(signal, "07_ZJets"); MVARead(signal, "09_TTV"); MVARead(signal, "04_TTTo2L2Nu"); MVARead(signal, "05_ST"); MVARead(signal, "00_Fakes"); } }
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // // MVA // //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ void MVA(float metPfType1_cut = 0., float mt2ll_cut = 0., TString signal = "ttDM0001scalar00010", // bool doMVATrain = 1, bool doMVARead = 1) { if (!doMVATrain && !doMVARead) return; gInterpreter->ExecuteMacro("../test/PaperStyle.C"); gSystem->mkdir(trainingdir, kTRUE); gSystem->mkdir(applicationdir, kTRUE); (TMVA::gConfig().GetIONames()).fWeightFileDir = weightsdir; // Training //---------------------------------------------------------------------------- if (doMVATrain) MVATrain(metPfType1_cut,mt2ll_cut,signal); // Reading //---------------------------------------------------------------------------- if (doMVARead) { //TString MVA_id = Form("met%.0f_mt2ll%.0f", metPfType1_cut, mt2ll_cut); TString MVA_id = Form("mt2ll%.0f", mt2ll_cut); MVARead(MVA_id, signal, "00_Fakes_reduced_1outof6"); MVARead(MVA_id, signal, "01_Data_reduced_1outof6"); MVARead(MVA_id, signal, "02_WZTo3LNu"); MVARead(MVA_id, signal, "03_VZ"); MVARead(MVA_id, signal, "04_TTTo2L2Nu"); MVARead(MVA_id, signal, "05_ST"); MVARead(MVA_id, signal, "06_WW"); MVARead(MVA_id, signal, "07_ZJets"); MVARead(MVA_id, signal, "09_TTV"); MVARead(MVA_id, signal, "10_HWW"); MVARead(MVA_id, signal, "11_Wg"); MVARead(MVA_id, signal, "12_Zg"); MVARead(MVA_id, signal, "13_VVV"); MVARead(MVA_id, signal, "14_HZ"); MVARead(MVA_id, signal, "15_WgStar"); MVARead(MVA_id, signal, signal); //MVARead(MVA_id, signal, "ttDM0001scalar00010"); //MVARead(MVA_id, signal, "ttDM0001scalar00020"); //MVARead(MVA_id, signal, "ttDM0001scalar00050"); } }