Exemplos de MVATrain em C++ (Cpp)

Exemplo n.º 1

Arquivo: MVA.C Projeto: cedricpri/AnalysisCMS

//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
// 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");

    }
}

Exemplo n.º 2

Arquivo: MVA.C Projeto: jgferrero/multivariateDM

//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
// 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");
    }
}

Exemplo n.º 3

Arquivo: MVA.C Projeto: amanjong/AnalysisCMS

//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
// 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");
    }
}