//Creates rootFile with a TTree with the ouput of MVA in it.
void MVAOutput(string Cat = "NoEleMatch_woGwoGSF_Barrel"
)
{
// TChain *inputTree = new TChain("TrainTree");
TChain *inputTree = new TChain("TestTree");
TString pathToFile = "tmva/tmvaRoot/";
// inputTree->Add(pathToFile+Form("/TMVA_v14EleVeto_%s.root",Cat.data()));
inputTree->Add(pathToFile+Form("/V18/TMVA_v18EleVeto_%s.root",Cat.data()));
if(inputTree->IsZombie()){
cout<<"ZOOOOMBIEEE"<<endl;
return;
}
int nEntries = inputTree->GetEntries();
// std::string outputFileName = Form("/data_CMS/cms/ivo/AntiEMVANewTraining/MVAOutput/V14/MvaOutput_Test_%s.root",Cat.data());
std::string outputFileName = Form("/data_CMS/cms/ivo/AntiEMVANewTraining/MVAOutput/V18/MvaOutput_All-TestTree_%s.root",Cat.data());
TFile* outputFile = new TFile (outputFileName.data(),"RECREATE");
TTree* mytree = new TTree("tree", "tree");
float t_NoEleMatch_woGwoGSF_Barrel;
float t_NoEleMatch_woGwGSF_Barrel;
float t_NoEleMatch_wGwoGSF_Barrel;
float t_NoEleMatch_wGwGSF_Barrel;
float t_woGwoGSF_Barrel;
float t_woGwGSF_Barrel;
float t_wGwoGSF_Barrel;
float t_wGwGSF_Barrel;
float t_NoEleMatch_woGwoGSF_Endcap;
float t_NoEleMatch_woGwGSF_Endcap;
float t_NoEleMatch_wGwoGSF_Endcap;
float t_NoEleMatch_wGwGSF_Endcap;
float t_woGwoGSF_Endcap;
float t_woGwGSF_Endcap;
float t_wGwoGSF_Endcap;
float t_wGwGSF_Endcap;
int t_classID;
mytree->Branch("NoEleMatch_woGwoGSF_Barrel",&t_NoEleMatch_woGwoGSF_Barrel,"NoEleMatch_woGwoGSF_Barrel/F");
mytree->Branch("NoEleMatch_woGwGSF_Barrel",&t_NoEleMatch_woGwGSF_Barrel,"NoEleMatch_woGwGSF_Barrel/F");
mytree->Branch("NoEleMatch_wGwoGSF_Barrel",&t_NoEleMatch_wGwoGSF_Barrel,"NoEleMatch_wGwoGSF_Barrel/F");
mytree->Branch("NoEleMatch_wGwGSF_Barrel",&t_NoEleMatch_wGwGSF_Barrel,"NoEleMatch_wGwGSF_Barrel/F");
mytree->Branch("woGwoGSF_Barrel",&t_woGwoGSF_Barrel,"woGwoGSF_Barrel/F");
mytree->Branch("woGwGSF_Barrel",&t_woGwGSF_Barrel,"woGwGSF_Barrel/F");
mytree->Branch("wGwoGSF_Barrel",&t_wGwoGSF_Barrel,"wGwoGSF_Barrel/F");
mytree->Branch("wGwGSF_Barrel",&t_wGwGSF_Barrel,"wGwGSF_Barrel/F");
mytree->Branch("NoEleMatch_woGwoGSF_Endcap",&t_NoEleMatch_woGwoGSF_Endcap,"NoEleMatch_woGwoGSF_Endcap/F");
mytree->Branch("NoEleMatch_woGwGSF_Endcap",&t_NoEleMatch_woGwGSF_Endcap,"NoEleMatch_woGwGSF_Endcap/F");
mytree->Branch("NoEleMatch_wGwoGSF_Endcap",&t_NoEleMatch_wGwoGSF_Endcap,"NoEleMatch_wGwoGSF_Endcap/F");
mytree->Branch("NoEleMatch_wGwGSF_Endcap",&t_NoEleMatch_wGwGSF_Endcap,"NoEleMatch_wGwGSF_Endcap/F");
mytree->Branch("woGwoGSF_Endcap",&t_woGwoGSF_Endcap,"woGwoGSF_Endcap/F");
mytree->Branch("woGwGSF_Endcap",&t_woGwGSF_Endcap,"woGwGSF_Endcap/F");
mytree->Branch("wGwoGSF_Endcap",&t_wGwoGSF_Endcap,"wGwoGSF_Endcap/F");
mytree->Branch("wGwGSF_Endcap",&t_wGwGSF_Endcap,"wGwGSF_Endcap/F");
mytree->Branch("classID",&t_classID,"classID/I");
// float Tau_EtaAtEcalEntrance;
// float Tau_Pt;
// float Tau_EmFraction;
// float Tau_HadrHoP;
// float Tau_HadrEoP;
// float Tau_VisMass;
float Tau_dCrackPhi;
float Tau_dCrackEta;
// float Tau_DecayMode;
float Tau_NumHitsVariable;
// //PFTaus variables
inputTree->SetBranchAddress("Tau_dCrackEta", &Tau_dCrackEta);
inputTree->SetBranchAddress("Tau_dCrackPhi", &Tau_dCrackPhi);
inputTree->SetBranchAddress("_Tau_GSFNumHits_M_Tau_KFNumHits__D__Tau_GSFNumHits_P_Tau_KFNumHits_", &Tau_NumHitsVariable);
int Tau_GsfEleMatch;
int Tau_GenEleMatch;
int Tau_GenEleFromZMatch;
int Tau_GenEleFromZTauTauMatch;
int Tau_GenHadMatch;
int Tau_GenJetMatch;
// float Tau_AbsEta;
float Tau_Eta;
float Tau_EtaAtEcalEntrance;
float Tau_PhiAtEcalEntrance;
float Tau_Pt;
float Tau_KFNumHits;
float Tau_HasGsf;
float Tau_GSFChi2;
float Tau_GSFNumHits;
float Tau_GSFTrackResol;
float Tau_GSFTracklnPt;
float Tau_GSFTrackEta;
float Tau_EmFraction;
float Tau_NumChargedCands;
float Tau_NumGammaCands;
float Tau_HadrHoP;
float Tau_HadrEoP;
float Tau_VisMass;
float Tau_GammaEtaMom;
float Tau_GammaPhiMom;
float Tau_GammaEnFrac;
float Tau_HadrMva;
// float Tau_IsoMVALoose;
float Tau_LooseComb3HitsIso;
float Tau_AntiELoose;
float Tau_AntiEMedium;
float Tau_AntiETight;
float Tau_AntiELooseMVA3;
float Tau_AntiEMediumMVA3;
float Tau_AntiETightMVA3;
float Tau_AntiEVTightMVA3;
float Tau_MatchElePassVeto;
float Tau_DecayMode;
float Elec_AbsEta;
float Elec_Pt;
float Elec_HasSC;
float Elec_PFMvaOutput;
float Elec_EtotOverPin;
float Elec_EeOverPout;
float Elec_EgammaOverPdif;
int Elec_EarlyBrem;//
int Elec_LateBrem;//
float Elec_Logsihih;
float Elec_DeltaEta;
float Elec_HoHplusE;
float Elec_Fbrem;
float Elec_HasGSF;
float Elec_Chi2GSF;
float Elec_GSFNumHits;
float Elec_GSFTrackResol;
float Elec_GSFTracklnPt;
float Elec_GSFTrackEta;
float ElecVeto_Pt;
int classID;
//GsfElectron variables
inputTree->SetBranchAddress("Elec_AbsEta", &Elec_AbsEta);
inputTree->SetBranchAddress("Elec_Pt", &Elec_Pt);
inputTree->SetBranchAddress("Elec_HasSC", &Elec_HasSC);
inputTree->SetBranchAddress("Elec_PFMvaOutput", &Elec_PFMvaOutput);
inputTree->SetBranchAddress("Elec_EtotOverPin", &Elec_EtotOverPin);
inputTree->SetBranchAddress("Elec_EeOverPout", &Elec_EeOverPout);
inputTree->SetBranchAddress("Elec_EgammaOverPdif", &Elec_EgammaOverPdif);
inputTree->SetBranchAddress("Elec_EarlyBrem", &Elec_EarlyBrem);
inputTree->SetBranchAddress("Elec_LateBrem", &Elec_LateBrem);
inputTree->SetBranchAddress("Elec_Logsihih", &Elec_Logsihih);
inputTree->SetBranchAddress("Elec_DeltaEta", &Elec_DeltaEta);
inputTree->SetBranchAddress("Elec_HoHplusE", &Elec_HoHplusE);
inputTree->SetBranchAddress("Elec_Fbrem", &Elec_Fbrem);
inputTree->SetBranchAddress("Elec_HasGSF", &Elec_HasGSF);
inputTree->SetBranchAddress("Elec_Chi2GSF", &Elec_Chi2GSF);
inputTree->SetBranchAddress("Elec_GSFNumHits", &Elec_GSFNumHits);
inputTree->SetBranchAddress("Elec_GSFTrackResol", &Elec_GSFTrackResol);
inputTree->SetBranchAddress("Elec_GSFTracklnPt", &Elec_GSFTracklnPt);
inputTree->SetBranchAddress("Elec_GSFTrackEta", &Elec_GSFTrackEta);
//PFTaus variables
inputTree->SetBranchAddress("Tau_GsfEleMatch", &Tau_GsfEleMatch);
inputTree->SetBranchAddress("Tau_GenEleMatch", &Tau_GenEleMatch);
inputTree->SetBranchAddress("Tau_GenEleFromZMatch", &Tau_GenEleFromZMatch);
inputTree->SetBranchAddress("Tau_GenEleFromZTauTauMatch", &Tau_GenEleFromZTauTauMatch);
inputTree->SetBranchAddress("Tau_GenHadMatch", &Tau_GenHadMatch);
inputTree->SetBranchAddress("Tau_GenJetMatch", &Tau_GenJetMatch);
inputTree->SetBranchAddress("Tau_EtaAtEcalEntrance", &Tau_EtaAtEcalEntrance);
inputTree->SetBranchAddress("Tau_PhiAtEcalEntrance", &Tau_PhiAtEcalEntrance);
inputTree->SetBranchAddress("Tau_Pt", &Tau_Pt);
inputTree->SetBranchAddress("Tau_Eta", &Tau_Eta);
inputTree->SetBranchAddress("Tau_HasGsf", &Tau_HasGsf);
inputTree->SetBranchAddress("Tau_GSFChi2", &Tau_GSFChi2);
inputTree->SetBranchAddress("Tau_GSFNumHits", &Tau_GSFNumHits);
inputTree->SetBranchAddress("Tau_GSFTrackResol", &Tau_GSFTrackResol);
inputTree->SetBranchAddress("Tau_GSFTracklnPt", &Tau_GSFTracklnPt);
inputTree->SetBranchAddress("Tau_GSFTrackEta", &Tau_GSFTrackEta);
inputTree->SetBranchAddress("Tau_KFNumHits", &Tau_KFNumHits);
inputTree->SetBranchAddress("Tau_EmFraction", &Tau_EmFraction);
inputTree->SetBranchAddress("Tau_NumChargedCands", &Tau_NumChargedCands);
inputTree->SetBranchAddress("Tau_NumGammaCands", &Tau_NumGammaCands);
inputTree->SetBranchAddress("Tau_HadrHoP", &Tau_HadrHoP);
inputTree->SetBranchAddress("Tau_HadrEoP", &Tau_HadrEoP);
inputTree->SetBranchAddress("Tau_VisMass", &Tau_VisMass);
inputTree->SetBranchAddress("Tau_GammaEtaMom", &Tau_GammaEtaMom);
inputTree->SetBranchAddress("Tau_GammaPhiMom", &Tau_GammaPhiMom);
inputTree->SetBranchAddress("Tau_GammaEnFrac", &Tau_GammaEnFrac);
inputTree->SetBranchAddress("Tau_HadrMva", &Tau_HadrMva);
// inputTree->SetBranchAddress("Tau_IsoMVALoose", &Tau_IsoMVALoose);
inputTree->SetBranchAddress("Tau_LooseComb3HitsIso", &Tau_LooseComb3HitsIso);
inputTree->SetBranchAddress("Tau_AntiELoose", &Tau_AntiELoose);
inputTree->SetBranchAddress("Tau_AntiEMedium", &Tau_AntiEMedium);
inputTree->SetBranchAddress("Tau_AntiETight", &Tau_AntiETight);
inputTree->SetBranchAddress("Tau_AntiELooseMVA3", &Tau_AntiELooseMVA3);
inputTree->SetBranchAddress("Tau_AntiEMediumMVA3", &Tau_AntiEMediumMVA3);
inputTree->SetBranchAddress("Tau_AntiETightMVA3", &Tau_AntiETightMVA3);
inputTree->SetBranchAddress("Tau_AntiEVTightMVA3", &Tau_AntiEVTightMVA3);
inputTree->SetBranchAddress("ElecVeto_Pt",&ElecVeto_Pt);
inputTree->SetBranchAddress("Tau_DecayMode", &Tau_DecayMode);
inputTree->SetBranchAddress("Tau_MatchElePassVeto", &Tau_MatchElePassVeto);
inputTree->SetBranchAddress("classID", &classID);
inputTree->SetBranchStatus("*", 1);
// string Weight_NoEleMatch_woGwoGSF_BL = "tmva/weights/V14/TMVAClassification_v14EleVeto_NoEleMatch_woGwoGSF_Barrel_BDTG.weights.xml";
// string Weight_NoEleMatch_woGwGSF_BL = "tmva/weights/V14/TMVAClassification_v14EleVeto_NoEleMatch_woGwGSF_Barrel_BDTG.weights.xml";
// string Weight_NoEleMatch_wGwoGSF_BL = "tmva/weights/V14/TMVAClassification_v14EleVeto_NoEleMatch_wGwoGSF_Barrel_BDTG.weights.xml";
// string Weight_NoEleMatch_wGwGSF_BL = "tmva/weights/V14/TMVAClassification_v14EleVeto_NoEleMatch_wGwGSF_Barrel_BDTG.weights.xml";
// string Weight_woGwoGSF_BL = "tmva/weights/V14/TMVAClassification_v14EleVeto_woGwoGSF_Barrel_BDTG.weights.xml";
// string Weight_woGwGSF_BL = "tmva/weights/V14/TMVAClassification_v14EleVeto_woGwGSF_Barrel_BDTG.weights.xml";
// string Weight_wGwoGSF_BL = "tmva/weights/V14/TMVAClassification_v14EleVeto_wGwoGSF_Barrel_BDTG.weights.xml";
// string Weight_wGwGSF_BL = "tmva/weights/V14/TMVAClassification_v14EleVeto_wGwGSF_Barrel_BDTG.weights.xml";
// string Weight_NoEleMatch_woGwoGSF_EC = "tmva/weights/V14/TMVAClassification_v14EleVeto_NoEleMatch_woGwoGSF_Endcap_BDTG.weights.xml";
// string Weight_NoEleMatch_woGwGSF_EC = "tmva/weights/V14/TMVAClassification_v14EleVeto_NoEleMatch_woGwGSF_Endcap_BDTG.weights.xml";
// string Weight_NoEleMatch_wGwoGSF_EC = "tmva/weights/V14/TMVAClassification_v14EleVeto_NoEleMatch_wGwoGSF_Endcap_BDTG.weights.xml";
// string Weight_NoEleMatch_wGwGSF_EC = "tmva/weights/V14/TMVAClassification_v14EleVeto_NoEleMatch_wGwGSF_Endcap_BDTG.weights.xml";
// string Weight_woGwoGSF_EC = "tmva/weights/V14/TMVAClassification_v14EleVeto_woGwoGSF_Endcap_BDTG.weights.xml";
// string Weight_woGwGSF_EC = "tmva/weights/V14/TMVAClassification_v14EleVeto_woGwGSF_Endcap_BDTG.weights.xml";
// string Weight_wGwoGSF_EC = "tmva/weights/V14/TMVAClassification_v14EleVeto_wGwoGSF_Endcap_BDTG.weights.xml";
// string Weight_wGwGSF_EC = "tmva/weights/V14/TMVAClassification_v14EleVeto_wGwGSF_Endcap_BDTG.weights.xml";
string Weight_NoEleMatch_woGwoGSF_BL = "tmva/weights/V18/TMVAClassification_v18EleVeto_NoEleMatch_woGwoGSF_Barrel_BDTG.weights.xml";
string Weight_NoEleMatch_woGwGSF_BL = "tmva/weights/V18/TMVAClassification_v18EleVeto_NoEleMatch_woGwGSF_Barrel_BDTG.weights.xml";
string Weight_NoEleMatch_wGwoGSF_BL = "tmva/weights/V18/TMVAClassification_v18EleVeto_NoEleMatch_wGwoGSF_Barrel_BDTG.weights.xml";
string Weight_NoEleMatch_wGwGSF_BL = "tmva/weights/V18/TMVAClassification_v18EleVeto_NoEleMatch_wGwGSF_Barrel_BDTG.weights.xml";
string Weight_woGwoGSF_BL = "tmva/weights/V18/TMVAClassification_v18EleVeto_woGwoGSF_Barrel_BDTG.weights.xml";
string Weight_woGwGSF_BL = "tmva/weights/V18/TMVAClassification_v18EleVeto_woGwGSF_Barrel_BDTG.weights.xml";
string Weight_wGwoGSF_BL = "tmva/weights/V18/TMVAClassification_v18EleVeto_wGwoGSF_Barrel_BDTG.weights.xml";
string Weight_wGwGSF_BL = "tmva/weights/V18/TMVAClassification_v18EleVeto_wGwGSF_Barrel_BDTG.weights.xml";
string Weight_NoEleMatch_woGwoGSF_EC = "tmva/weights/V18/TMVAClassification_v18EleVeto_NoEleMatch_woGwoGSF_Endcap_BDTG.weights.xml";
string Weight_NoEleMatch_woGwGSF_EC = "tmva/weights/V18/TMVAClassification_v18EleVeto_NoEleMatch_woGwGSF_Endcap_BDTG.weights.xml";
string Weight_NoEleMatch_wGwoGSF_EC = "tmva/weights/V18/TMVAClassification_v18EleVeto_NoEleMatch_wGwoGSF_Endcap_BDTG.weights.xml";
string Weight_NoEleMatch_wGwGSF_EC = "tmva/weights/V18/TMVAClassification_v18EleVeto_NoEleMatch_wGwGSF_Endcap_BDTG.weights.xml";
string Weight_woGwoGSF_EC = "tmva/weights/V18/TMVAClassification_v18EleVeto_woGwoGSF_Endcap_BDTG.weights.xml";
string Weight_woGwGSF_EC = "tmva/weights/V18/TMVAClassification_v18EleVeto_woGwGSF_Endcap_BDTG.weights.xml";
string Weight_wGwoGSF_EC = "tmva/weights/V18/TMVAClassification_v18EleVeto_wGwoGSF_Endcap_BDTG.weights.xml";
string Weight_wGwGSF_EC = "tmva/weights/V18/TMVAClassification_v18EleVeto_wGwGSF_Endcap_BDTG.weights.xml";
TMVA::Reader *reader_NoEleMatch_woGwoGSF_BL = new TMVA::Reader( "!Color:!Silent:Error" );
reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_Pt",&Tau_Pt);
reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi);
// reader_NoEleMatch_woGwoGSF_BL->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_NoEleMatch_woGwoGSF_BL->SetVerbose(kTRUE);
reader_NoEleMatch_woGwoGSF_BL->BookMVA("BDTG",Weight_NoEleMatch_woGwoGSF_BL);
TMVA::Reader *reader_NoEleMatch_woGwGSF_BL = new TMVA::Reader( "!Color:!Silent:Error" );
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_Pt",&Tau_Pt);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_HadrMva",&Tau_HadrMva);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_GSFChi2",&Tau_GSFChi2);
reader_NoEleMatch_woGwGSF_BL->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi);
// reader_NoEleMatch_woGwGSF_BL->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_NoEleMatch_woGwGSF_BL->SetVerbose(kTRUE);
reader_NoEleMatch_woGwGSF_BL->BookMVA("BDTG",Weight_NoEleMatch_woGwGSF_BL);
TMVA::Reader *reader_NoEleMatch_wGwoGSF_BL = new TMVA::Reader( "!Color:!Silent:Error" );
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_Pt",&Tau_Pt);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi);
// reader_NoEleMatch_wGwoGSF_BL->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_NoEleMatch_wGwoGSF_BL->SetVerbose(kTRUE);
reader_NoEleMatch_wGwoGSF_BL->BookMVA("BDTG",Weight_NoEleMatch_wGwoGSF_BL);
TMVA::Reader *reader_NoEleMatch_wGwGSF_BL = new TMVA::Reader( "!Color:!Silent:Error" );
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_Pt",&Tau_Pt);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_HadrMva",&Tau_HadrMva);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_GSFChi2",&Tau_GSFChi2);
reader_NoEleMatch_wGwGSF_BL->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi);
// reader_NoEleMatch_wGwGSF_BL->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_NoEleMatch_wGwGSF_BL->SetVerbose(kTRUE);
reader_NoEleMatch_wGwGSF_BL->BookMVA("BDTG",Weight_NoEleMatch_wGwGSF_BL);
TMVA::Reader *reader_woGwoGSF_BL = new TMVA::Reader( "!Color:!Silent:Error" );
reader_woGwoGSF_BL->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin);
reader_woGwoGSF_BL->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif);
reader_woGwoGSF_BL->AddVariable("Elec_Fbrem",&Elec_Fbrem);
reader_woGwoGSF_BL->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF);
reader_woGwoGSF_BL->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits);
reader_woGwoGSF_BL->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol);
reader_woGwoGSF_BL->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt);
reader_woGwoGSF_BL->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta);
reader_woGwoGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_woGwoGSF_BL->AddVariable("Tau_Pt",&Tau_Pt);
reader_woGwoGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_woGwoGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_woGwoGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_woGwoGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_woGwoGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
reader_woGwoGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi);
// reader_woGwoGSF_BL->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_woGwoGSF_BL->SetVerbose(kTRUE);
reader_woGwoGSF_BL->BookMVA("BDTG",Weight_woGwoGSF_BL);
TMVA::Reader *reader_woGwGSF_BL = new TMVA::Reader( "!Color:!Silent:Error" );
reader_woGwGSF_BL->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin);
reader_woGwGSF_BL->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif);
reader_woGwGSF_BL->AddVariable("Elec_Fbrem",&Elec_Fbrem);
reader_woGwGSF_BL->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF);
reader_woGwGSF_BL->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits);
reader_woGwGSF_BL->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol);
reader_woGwGSF_BL->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt);
reader_woGwGSF_BL->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta);
reader_woGwGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_woGwGSF_BL->AddVariable("Tau_Pt",&Tau_Pt);
reader_woGwGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_woGwGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_woGwGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_woGwGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_woGwGSF_BL->AddVariable("Tau_HadrMva",&Tau_HadrMva);
reader_woGwGSF_BL->AddVariable("Tau_GSFChi2",&Tau_GSFChi2);
reader_woGwGSF_BL->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable);
reader_woGwGSF_BL->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol);
reader_woGwGSF_BL->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt);
reader_woGwGSF_BL->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta);
reader_woGwGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
reader_woGwGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi);
// reader_woGwGSF_BL->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_woGwGSF_BL->SetVerbose(kTRUE);
reader_woGwGSF_BL->BookMVA("BDTG",Weight_woGwGSF_BL);
TMVA::Reader *reader_wGwoGSF_BL = new TMVA::Reader( "!Color:!Silent:Error" );
reader_wGwoGSF_BL->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin);
reader_wGwoGSF_BL->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif);
reader_wGwoGSF_BL->AddVariable("Elec_Fbrem",&Elec_Fbrem);
reader_wGwoGSF_BL->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF);
reader_wGwoGSF_BL->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits);
reader_wGwoGSF_BL->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol);
reader_wGwoGSF_BL->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt);
reader_wGwoGSF_BL->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta);
reader_wGwoGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_wGwoGSF_BL->AddVariable("Tau_Pt",&Tau_Pt);
reader_wGwoGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_wGwoGSF_BL->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands);
reader_wGwoGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_wGwoGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_wGwoGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_wGwoGSF_BL->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom);
reader_wGwoGSF_BL->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom);
reader_wGwoGSF_BL->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac);
reader_wGwoGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
reader_wGwoGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi);
// reader_wGwoGSF_BL->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_wGwoGSF_BL->SetVerbose(kTRUE);
reader_wGwoGSF_BL->BookMVA("BDTG",Weight_wGwoGSF_BL);
TMVA::Reader *reader_wGwGSF_BL = new TMVA::Reader( "!Color:!Silent:Error" );
reader_wGwGSF_BL->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin);
reader_wGwGSF_BL->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif);
reader_wGwGSF_BL->AddVariable("Elec_Fbrem",&Elec_Fbrem);
reader_wGwGSF_BL->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF);
reader_wGwGSF_BL->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits);
reader_wGwGSF_BL->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol);
reader_wGwGSF_BL->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt);
reader_wGwGSF_BL->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta);
reader_wGwGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_wGwGSF_BL->AddVariable("Tau_Pt",&Tau_Pt);
reader_wGwGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_wGwGSF_BL->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands);
reader_wGwGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_wGwGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_wGwGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_wGwGSF_BL->AddVariable("Tau_HadrMva",&Tau_HadrMva);
reader_wGwGSF_BL->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom);
reader_wGwGSF_BL->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom);
reader_wGwGSF_BL->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac);
reader_wGwGSF_BL->AddVariable("Tau_GSFChi2",&Tau_GSFChi2);
reader_wGwGSF_BL->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable);
reader_wGwGSF_BL->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol);
reader_wGwGSF_BL->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt);
reader_wGwGSF_BL->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta);
reader_wGwGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
reader_wGwGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi);
// reader_wGwGSF_BL->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_wGwGSF_BL->SetVerbose(kTRUE);
reader_wGwGSF_BL->BookMVA("BDTG",Weight_wGwGSF_BL);
////////////////////////
TMVA::Reader *reader_NoEleMatch_woGwoGSF_EC = new TMVA::Reader( "!Color:!Silent:Error" );
reader_NoEleMatch_woGwoGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_NoEleMatch_woGwoGSF_EC->AddVariable("Tau_Pt",&Tau_Pt);
reader_NoEleMatch_woGwoGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_NoEleMatch_woGwoGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_NoEleMatch_woGwoGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_NoEleMatch_woGwoGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_NoEleMatch_woGwoGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
// reader_NoEleMatch_woGwoGSF_EC->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_NoEleMatch_woGwoGSF_EC->SetVerbose(kTRUE);
reader_NoEleMatch_woGwoGSF_EC->BookMVA("BDTG",Weight_NoEleMatch_woGwoGSF_EC);
TMVA::Reader *reader_NoEleMatch_woGwGSF_EC = new TMVA::Reader( "!Color:!Silent:Error" );
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_Pt",&Tau_Pt);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_HadrMva",&Tau_HadrMva);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_GSFChi2",&Tau_GSFChi2);
reader_NoEleMatch_woGwGSF_EC->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta);
reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
// reader_NoEleMatch_woGwGSF_EC->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_NoEleMatch_woGwGSF_EC->SetVerbose(kTRUE);
reader_NoEleMatch_woGwGSF_EC->BookMVA("BDTG",Weight_NoEleMatch_woGwGSF_EC);
TMVA::Reader *reader_NoEleMatch_wGwoGSF_EC = new TMVA::Reader( "!Color:!Silent:Error" );
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_Pt",&Tau_Pt);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac);
reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
// reader_NoEleMatch_wGwoGSF_EC->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_NoEleMatch_wGwoGSF_EC->SetVerbose(kTRUE);
reader_NoEleMatch_wGwoGSF_EC->BookMVA("BDTG",Weight_NoEleMatch_wGwoGSF_EC);
TMVA::Reader *reader_NoEleMatch_wGwGSF_EC = new TMVA::Reader( "!Color:!Silent:Error" );
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_Pt",&Tau_Pt);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_HadrMva",&Tau_HadrMva);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_GSFChi2",&Tau_GSFChi2);
reader_NoEleMatch_wGwGSF_EC->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta);
reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
// reader_NoEleMatch_wGwGSF_EC->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_NoEleMatch_wGwGSF_EC->SetVerbose(kTRUE);
reader_NoEleMatch_wGwGSF_EC->BookMVA("BDTG",Weight_NoEleMatch_wGwGSF_EC);
TMVA::Reader *reader_woGwoGSF_EC = new TMVA::Reader( "!Color:!Silent:Error" );
reader_woGwoGSF_EC->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin);
reader_woGwoGSF_EC->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif);
reader_woGwoGSF_EC->AddVariable("Elec_Fbrem",&Elec_Fbrem);
reader_woGwoGSF_EC->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF);
reader_woGwoGSF_EC->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits);
reader_woGwoGSF_EC->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol);
reader_woGwoGSF_EC->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt);
reader_woGwoGSF_EC->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta);
reader_woGwoGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_woGwoGSF_EC->AddVariable("Tau_Pt",&Tau_Pt);
reader_woGwoGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_woGwoGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_woGwoGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_woGwoGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_woGwoGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
// reader_woGwoGSF_EC->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_woGwoGSF_EC->SetVerbose(kTRUE);
reader_woGwoGSF_EC->BookMVA("BDTG",Weight_woGwoGSF_EC);
TMVA::Reader *reader_woGwGSF_EC = new TMVA::Reader( "!Color:!Silent:Error" );
reader_woGwGSF_EC->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin);
reader_woGwGSF_EC->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif);
reader_woGwGSF_EC->AddVariable("Elec_Fbrem",&Elec_Fbrem);
reader_woGwGSF_EC->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF);
reader_woGwGSF_EC->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits);
reader_woGwGSF_EC->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol);
reader_woGwGSF_EC->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt);
reader_woGwGSF_EC->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta);
reader_woGwGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_woGwGSF_EC->AddVariable("Tau_Pt",&Tau_Pt);
reader_woGwGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_woGwGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_woGwGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_woGwGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_woGwGSF_EC->AddVariable("Tau_HadrMva",&Tau_HadrMva);
reader_woGwGSF_EC->AddVariable("Tau_GSFChi2",&Tau_GSFChi2);
reader_woGwGSF_EC->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable);
reader_woGwGSF_EC->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol);
reader_woGwGSF_EC->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt);
reader_woGwGSF_EC->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta);
reader_woGwGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
// reader_woGwGSF_EC->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_woGwGSF_EC->SetVerbose(kTRUE);
reader_woGwGSF_EC->BookMVA("BDTG",Weight_woGwGSF_EC);
TMVA::Reader *reader_wGwoGSF_EC = new TMVA::Reader( "!Color:!Silent:Error" );
reader_wGwoGSF_EC->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin);
reader_wGwoGSF_EC->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif);
reader_wGwoGSF_EC->AddVariable("Elec_Fbrem",&Elec_Fbrem);
reader_wGwoGSF_EC->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF);
reader_wGwoGSF_EC->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits);
reader_wGwoGSF_EC->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol);
reader_wGwoGSF_EC->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt);
reader_wGwoGSF_EC->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta);
reader_wGwoGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_wGwoGSF_EC->AddVariable("Tau_Pt",&Tau_Pt);
reader_wGwoGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_wGwoGSF_EC->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands);
reader_wGwoGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_wGwoGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_wGwoGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_wGwoGSF_EC->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom);
reader_wGwoGSF_EC->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom);
reader_wGwoGSF_EC->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac);
reader_wGwoGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
// reader_wGwoGSF_EC->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_wGwoGSF_EC->SetVerbose(kTRUE);
reader_wGwoGSF_EC->BookMVA("BDTG",Weight_wGwoGSF_EC);
TMVA::Reader *reader_wGwGSF_EC = new TMVA::Reader( "!Color:!Silent:Error" );
reader_wGwGSF_EC->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin);
reader_wGwGSF_EC->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif);
reader_wGwGSF_EC->AddVariable("Elec_Fbrem",&Elec_Fbrem);
reader_wGwGSF_EC->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF);
reader_wGwGSF_EC->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits);
reader_wGwGSF_EC->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol);
reader_wGwGSF_EC->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt);
reader_wGwGSF_EC->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta);
reader_wGwGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance);
reader_wGwGSF_EC->AddVariable("Tau_Pt",&Tau_Pt);
reader_wGwGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction);
reader_wGwGSF_EC->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands);
reader_wGwGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP);
reader_wGwGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP);
reader_wGwGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass);
reader_wGwGSF_EC->AddVariable("Tau_HadrMva",&Tau_HadrMva);
reader_wGwGSF_EC->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom);
reader_wGwGSF_EC->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom);
reader_wGwGSF_EC->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac);
reader_wGwGSF_EC->AddVariable("Tau_GSFChi2",&Tau_GSFChi2);
reader_wGwGSF_EC->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable);
reader_wGwGSF_EC->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol);
reader_wGwGSF_EC->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt);
reader_wGwGSF_EC->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta);
reader_wGwGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta);
// reader_wGwGSF_EC->AddVariable("Tau_DecayMode",&Tau_DecayMode);
reader_wGwGSF_EC->SetVerbose(kTRUE);
reader_wGwGSF_EC->BookMVA("BDTG",Weight_wGwGSF_EC);
cout<< "Number of entries : "<<nEntries<<endl;
for (int iEntry = 0; iEntry<nEntries ; iEntry++){
if(iEntry%100000==0) cout << iEntry <<"/"<<nEntries << endl;
inputTree->GetEntry(iEntry);
t_NoEleMatch_woGwoGSF_Barrel = -99;
t_NoEleMatch_woGwGSF_Barrel = -99;
t_NoEleMatch_wGwoGSF_Barrel = -99;
t_NoEleMatch_wGwGSF_Barrel = -99;
t_woGwoGSF_Barrel = -99;
t_woGwGSF_Barrel = -99;
t_wGwoGSF_Barrel = -99;
t_wGwGSF_Barrel = -99;
t_NoEleMatch_woGwoGSF_Endcap = -99;
t_NoEleMatch_woGwGSF_Endcap = -99;
t_NoEleMatch_wGwoGSF_Endcap = -99;
t_NoEleMatch_wGwGSF_Endcap = -99;
t_woGwoGSF_Endcap = -99;
t_woGwGSF_Endcap = -99;
t_wGwoGSF_Endcap = -99;
t_wGwGSF_Endcap = -99;
// cout<<"classID :"<<classID<<endl;
t_classID = classID;
if (Cat.find("NoEleMatch_woGwoGSF_Barrel")!=std::string::npos) t_NoEleMatch_woGwoGSF_Barrel = reader_NoEleMatch_woGwoGSF_BL->EvaluateMVA("BDTG");
if (Cat.find("NoEleMatch_woGwGSF_Barrel")!=std::string::npos) t_NoEleMatch_woGwGSF_Barrel = reader_NoEleMatch_woGwGSF_BL->EvaluateMVA("BDTG");
if (Cat.find("NoEleMatch_wGwoGSF_Barrel")!=std::string::npos) t_NoEleMatch_wGwoGSF_Barrel = reader_NoEleMatch_wGwoGSF_BL->EvaluateMVA("BDTG");
if (Cat.find("NoEleMatch_wGwGSF_Barrel")!=std::string::npos) t_NoEleMatch_wGwGSF_Barrel = reader_NoEleMatch_wGwGSF_BL->EvaluateMVA("BDTG");
if (Cat.find("woGwoGSF_Barrel")!=std::string::npos) t_woGwoGSF_Barrel = reader_woGwoGSF_BL->EvaluateMVA("BDTG");
if (Cat.find("woGwGSF_Barrel")!=std::string::npos) t_woGwGSF_Barrel = reader_woGwGSF_BL->EvaluateMVA("BDTG");
if (Cat.find("wGwoGSF_Barrel")!=std::string::npos) t_wGwoGSF_Barrel = reader_wGwoGSF_BL->EvaluateMVA("BDTG");
if (Cat.find("wGwGSF_Barrel")!=std::string::npos) t_wGwGSF_Barrel = reader_wGwGSF_BL->EvaluateMVA("BDTG");
if (Cat.find("NoEleMatch_woGwoGSF_Endcap")!=std::string::npos) t_NoEleMatch_woGwoGSF_Endcap = reader_NoEleMatch_woGwoGSF_EC->EvaluateMVA("BDTG");
if (Cat.find("NoEleMatch_woGwGSF_Endcap")!=std::string::npos) t_NoEleMatch_woGwGSF_Endcap = reader_NoEleMatch_woGwGSF_EC->EvaluateMVA("BDTG");
if (Cat.find("NoEleMatch_wGwoGSF_Endcap")!=std::string::npos) t_NoEleMatch_wGwoGSF_Endcap = reader_NoEleMatch_wGwoGSF_EC->EvaluateMVA("BDTG");
if (Cat.find("NoEleMatch_wGwGSF_Endcap")!=std::string::npos) t_NoEleMatch_wGwGSF_Endcap = reader_NoEleMatch_wGwGSF_EC->EvaluateMVA("BDTG");
if (Cat.find("woGwoGSF_Endcap")!=std::string::npos) t_woGwoGSF_Endcap = reader_woGwoGSF_EC->EvaluateMVA("BDTG");
if (Cat.find("woGwGSF_Endcap")!=std::string::npos) t_woGwGSF_Endcap = reader_woGwGSF_EC->EvaluateMVA("BDTG");
if (Cat.find("wGwoGSF_Endcap")!=std::string::npos) t_wGwoGSF_Endcap = reader_wGwoGSF_EC->EvaluateMVA("BDTG");
if (Cat.find("wGwGSF_Endcap")!=std::string::npos) t_wGwGSF_Endcap = reader_wGwGSF_EC->EvaluateMVA("BDTG");
mytree->Fill();
}
mytree->Write();
cout<<"Creating file : "<<outputFileName.data()<<endl;
// inputFile->Close();
outputFile->Close();
return;
}
