void plotHistos_Wen() {
gROOT->LoadMacro("tdrstyle.C");
gROOT->ProcessLine("setTDRStyle()");
TH1::SetDefaultSumw2(1);
//gROOT->SetBatch(1);
TString plotdir = "plots/";
if (gSystem->AccessPathName(plotdir))
gSystem->mkdir(plotdir);
////////////////////////////////////////////////////////////////////////////
// Task 1 (a) //
// - please comment out other tasks //
////////////////////////////////////////////////////////////////////////////
/*
// Read from ntuples
Events * ev = new Events();
TCut cutmc_all = "Vtype==0"; // change this to your channel
TCut cutdata_all = "Vtype==0"; // change this to your channel
TString process = "VH"; // can try other processes: VH, WJ, ZJ, TT
ev->read(cutmc_all, cutdata_all, process);
TString var = "HCSV_mass"; // the variable to plot
TCut cut = "V_pt>120"; // the selection cut
TString title = ";m(jj) [GeV]"; // the title of the histogram
TString plotname = process + "_Hmass"; // the name of the image file
int nbinsx = 15; // number of bins
double xlow = 30.0; // the low edge of x-axis
double xup = 255.0; // the upper edge of x-axis
TString options = "!plotLog:!plotNorm"; // use "plotLog" to plot on log-y scale,
// use "!plotLog" to plot on linear scale;
// use "plotNorm" to plot normalized plots,
// use "!plotNorm" otherwise.
// Use "ev->ZH" for ZH, or "ev->WjLF", "ev->WjHF", "ev->ZjLF", "ev->ZjHF", "ev->TT" for other processes
MakePlot(ev->ZH, var, cut, title, nbinsx, xlow, xup, plotname, plotdir, options);
// You can put in the parameters directly as in the following commented out line:
//MakePlot(ev->ZH, "HCSV_pt", cut, "; p_{T}(jj) [GeV]", 16, 0, 240., process + "_Hpt", plotdir, options);
*/
////////////////////////////////////////////////////////////////////////////
// Task 1 (b) //
// - please comment out other tasks //
////////////////////////////////////////////////////////////////////////////
/*
// Read from ntuples
Events * ev = new Events();
TCut cutmc_all = "Vtype==0"; // change this to your channel
TCut cutdata_all = "Vtype==0"; // change this to your channel
ev->read(cutmc_all, cutdata_all, "VH:ZJ"); // read both VH and ZJ processes
TString var = "HCSV_mass";
TCut cut = "V_pt>120";
//TCut cut = "V_pt>120 && Jet_pt[hJCidx[0]]>30 && Jet_pt[hJCidx[1]]>30"; // for Wln, Znn, change to tighter cut
TString title = ";m(jj) [GeV]";
TString plotname = "ZH_vs_ZJ_Hmass";
int nbinsx = 15;
double xlow = 30.0;
double xup = 255.0;
TString options = "!plotLog:plotNorm";
// Using "ev->ZH" for ZH and "ev->ZjHF" for Z+HF
MakePlot2(ev->ZH, ev->ZjHF, var, cut, title, nbinsx, xlow, xup, plotname, plotdir, options);
*/
////////////////////////////////////////////////////////////////////////////
// Task 2 //
// - please comment out other tasks //
////////////////////////////////////////////////////////////////////////////
// Zmm______________________________________________________________________
//TString channel = "Zmm";
//Loose cuts for this channel
//TCut cutmc_all = "Vtype==0 && V_pt>100 && met_pt<60"; // V related
//cutmc_all += "75<V_mass && V_mass<105"; //Z mass window
//cutmc_all += "Jet_pt[hJCidx[0]]>20 && Jet_pt[hJCidx[1]]>20 && abs(Jet_eta[hJCidx[0]])<2.4 && abs(Jet_eta[hJCidx[1]])<2.4";// H jets
//cutmc_all += "min(Jet_btagCSV[hJCidx[0]], Jet_btagCSV[hJCidx[1]])>0.4"; //loose b-tag
//cutmc_all += "Sum$(Jet_pt > 20 & abs(Jet_eta) < 2.4) < 3"; //no additional jets
//TCut met_cleaning = "Flag_hbheFilterNew &&Flag_hbheIsoFilter && Flag_goodVertices &&Flag_eeBadScFilter &&Flag_CSCTightHaloFilter";
//TCut mc_trigger = "HLT_ZmmHbbAll==1";
//TCut data_trigger = "HLT_BIT_HLT_IsoMu18_v==1";
//TCut cutdata_all = cutmc_all;
//cutdata_all += data_trigger;
//cutmc_all += mc_trigger;
// Scale factors in order of: WjLF, WjHF, ZjLF, ZjHF, TT
// NOTE: WjLF, WjHF are not needed for Zll
//double scalefactors[5] = {1.00, 1.00, 1.00, 1.00, 1.00};
// Zee______________________________________________________________________
//TString channel = "Zee";
//Loose cuts for this channel
//TCut cutmc_all = "Vtype==1 && V_pt>100 && met_pt<60"; // V related
//cutmc_all += "75<V_mass && V_mass<105"; //Z mass window
//cutmc_all += "Jet_pt[hJCidx[0]]>20 && Jet_pt[hJCidx[1]]>20 && abs(Jet_eta[hJCidx[0]])<2.4 && abs(Jet_eta[hJCidx[1]])<2.4";// H jets
//cutmc_all += "min(Jet_btagCSV[hJCidx[0]], Jet_btagCSV[hJCidx[1]])>0.4"; //loose b-tag
//cutmc_all += "Sum$(Jet_pt > 20 & abs(Jet_eta) < 2.4) < 3"; //no additional jets
//TCut met_cleaning = "Flag_hbheFilterNew &&Flag_hbheIsoFilter && Flag_goodVertices &&Flag_eeBadScFilter &&Flag_CSCTightHaloFilter";
//TCut mc_trigger = "HLT_ZeeHbbAll==1";
//TCut data_trigger = "HLT_BIT_HLT_Ele27_eta2p1_WPLoose_Gsf_v==1";
//TCut cutdata_all = cutmc_all;
//cutdata_all += data_trigger;
//cutmc_all += mc_trigger;
// Scale factors in order of: WjLF, WjHF, ZjLF, ZjHF, TT
// NOTE: WjLF, WjHF are not needed for Zll
//double scalefactors[5] = {1.00, 1.00, 1.00, 1.00, 1.00};
// Wmn______________________________________________________________________
//TString channel = "Wmn";
//Loose cuts for this channel
//TCut cutmc_all = "Vtype==2 && V_pt>100 && met_pt>45"; // V related
//cutmc_all += "vLeptons_pt[0]>30 && abs(vLeptons_eta[0])<2.4 && naLeptons==0";//one lepton
//cutmc_all += "Jet_pt[hJCidx[0]]>20 && Jet_pt[hJCidx[1]]>20 && abs(Jet_eta[hJCidx[0]])<2.4 && abs(Jet_eta[hJCidx[1]])<2.4";// H jets
//cutmc_all += "min(Jet_btagCSV[hJCidx[0]], Jet_btagCSV[hJCidx[1]])>0.4"; //loose b-tag
//cutmc_all += "Sum$(Jet_pt > 20 & abs(Jet_eta) < 2.4) < 3"; //no additional jets
//TCut met_cleaning = "Flag_hbheFilterNew &&Flag_hbheIsoFilter && Flag_goodVertices &&Flag_eeBadScFilter &&Flag_CSCTightHaloFilter";
//TCut mc_trigger = "";
//TCut data_trigger = "";
//TCut cutdata_all = cutmc_all;
//cutdata_all += data_trigger;
//cutmc_all += mc_trigger;
// Scale factors in order of: WjLF, WjHF, ZjLF, ZjHF, TT
// NOTE: ZjLF, ZjHF are not needed for Wln
//double scalefactors[5] = {1.00, 1.00, 1.00, 1.00, 1.00};
// Wen______________________________________________________________________
TString channel = "Wen";
//Loose cuts for this channel
TCut cutmc_all = "Vtype==3 && V_pt>100 && met_pt>45"; // V related
cutmc_all += "vLeptons_pt[0]>30 && abs(vLeptons_eta[0])<2.4 && naLeptons==0";//one lepton
cutmc_all += "Jet_pt[hJCidx[0]]>20 && Jet_pt[hJCidx[1]]>20 && abs(Jet_eta[hJCidx[0]])<2.4 && abs(Jet_eta[hJCidx[1]])<2.4";// H jets
cutmc_all += "min(Jet_btagCSV[hJCidx[0]], Jet_btagCSV[hJCidx[1]])>0.4"; //loose b-tag
cutmc_all += "Sum$(Jet_pt > 20 & abs(Jet_eta) < 2.4) < 3"; //no additional jets
TCut met_cleaning = "Flag_hbheFilterNew &&Flag_hbheIsoFilter && Flag_goodVertices &&Flag_eeBadScFilter &&Flag_CSCTightHaloFilter";
TCut mc_trigger = "";
TCut data_trigger = "";
TCut cutdata_all = cutmc_all;
cutdata_all += data_trigger;
cutmc_all += mc_trigger;
// Scale factors in order of: WjLF, WjHF, ZjLF, ZjHF, TT
// NOTE: ZjLF, ZjHF are not needed for Wln
double scalefactors[5] = {1.00, 1.00, 1.00, 1.00, 1.00};
// Znn______________________________________________________________________
//TString channel = "Znn";
// These are loose cuts for all plots in this particular channel
//TCut cutmc_all = "Vtype==4 && HCSV_pt>130 && Jet_pt[hJCidx[0]]>80 && Jet_pt[hJCidx[1]]>30 && abs(Jet_eta[hJCidx[0]])<2.5 && abs(Jet_eta[hJCidx[1]])<2.5 && naLeptons==0 && Sum$(Jet_pt > 20 & abs(Jet_eta) < 2.4)<5";
//&& min(Min$(abs(deltaPhi(met_phi,hJet_phi))),Min$(abs(deltaPhiMETjets(met_phi,aJet_phi,aJet_pt,aJet_eta)))+999*(Sum$(aJet_pt>25 && abs(aJet_eta)<4.5 && aJet_id==1 && aJet_puJetIdL>0)==0) )>0.5";//needs to be updated to 13 TeV
//cutmc_all += "min(Jet_btagCSV[hJCidx[0]], Jet_btagCSV[hJCidx[1]])>0.4 && min(abs(HCSV_phi-V_phi),(2*TMath::Pi())-abs(HCSV_phi-V_phi))>2.0"; // tighter cut
//TCut cutdata_all = cutmc_all;
//cutmc_all *= "(triggerFlags[42]==1 || triggerFlags[39]==1 || triggerFlags[41]==1)"; // apply trigger bits and MET cleaning for MC (they are already applied on data)
//cutmc_all *= "triggercorrMET(met_pt)";
// Scale factors in order of: WjLF, WjHF, ZjLF, ZjHF, TT
//double scalefactors[5] = {1.00, 1.00, 1.00, 1.00, 1.00};
// All channels_____________________________________________________________
// Read from ntuples
Events * ev = new Events();
ev->read(cutmc_all, cutdata_all);
// Set the scale factors
ev->set_sf(scalefactors);
// Optimize these five variables (default: recommendations for Zll)
//double vpt = 150.;
//double hpt = 0.;
//double maxcsv = 0.679;
//double mincsv = 0.5;
//double dPhi = 0.;
// Optimize these five variables (default: recommendations for Wln)
double vpt = 150.;
double hpt = 100.;
double maxcsv = 0.898;
double mincsv = 0.5;
double dPhi = 2.95;
// Optimize these five variables (default: recommendations for Znn)
//double vpt = 170.; // for Znn, pT(V) = MET
//double hpt = 170.;
//double maxcsv = 0.898;
//double mincsv = 0.5;
//double dPhi = 2.95;
// If doing cut and count analysis, cut on HCSV_mass by changing the values of minhmass and maxhmass
//double minhmass = 0.;
//double maxhmass = 9999.;
double minhmass = 110.;
double maxhmass = 140.;
// These are tight cuts for this particular plot
//TCut cutmc = Form("V_pt>%.2f && HCSV_pt>%.2f && max(Jet_btagCSV[hJCidx[0]],Jet_btagCSV[hJCidx[1]])>%.3f && min(Jet_btagCSV[hJCidx[0]], Jet_btagCSV[hJCidx[1]])>%.3f && HVdPhi>%.2f && %.2f<HCSV_mass && HCSV_mass<%.2f", vpt, hpt, maxcsv, mincsv, dPhi, minhmass, maxhmass);
TCut cutmc = "";
TCut cutdata = cutmc;
TString var = "HCSV_mass";
TString title = ";m(jj) [GeV]";
TString plotname = channel + "_Hmass";
int nbinsx = 15;
double xlow = 30.0;
double xup = 255.0;
TString options = "printStat:plotSig:plotData:!plotLog";
MakePlots(ev, var, cutmc, cutdata, title, nbinsx, xlow, xup, plotname, plotdir, options);
// Or, just put in them directly as in the following commented out line:
//MakePlots(ev, "HCSV_mass", cutmc, cutdata, "m(jj) [GeV]", 14, 30.0, 255.0, channel+"_Hmass", plotdir, "printStat:plotSig:!plotData:!plotLog");
////////////////////////////////////////////////////////////////////////////
// Task 3 //
// - please comment out other tasks, but keep Task 2 //
////////////////////////////////////////////////////////////////////////////
/*
TString dcname = Form("vhbb_%s_13TeV.txt", channel.Data()); // the datacard name
TString wsname = plotdir + plotname +".root"; // the workspace name
bool useshapes = false;
TString options1 = "!unblind:SplusB";
// For cut-and-count analysis, apply HCSV_mass cut before calling MakeDatacard(...)
MakeDatacard(channel, dcname, wsname, useshapes, options1);
// For shape analysis, remove HCSV_mass cut before calling MakeDatacard(...)
//cutmc = Form("V_pt>%.2f && HCSV_pt>%.2f && max(Jet_btagCSV[hJCidx[0]],Jet_btagCSV[hJCidx[1]])>%.3f && min(Jet_btagCSV[hJCidx[0]], Jet_btagCSV[hJCidx[1]])>%.3f && min(abs(HCSV_phi-V_phi),(2*TMath::Pi())-abs(HCSV_phi-V_phi))>%.2f", vpt, hpt, maxcsv, mincsv, dPhi);
//cutdata = cutmc;
//plotname = channel + "_Hmass_shapes";
//MakePlots(ev, var, cutmc, cutdata, title, nbinsx, xlow, xup, plotname, plotdir, options);
//dcname = Form("vhbb_shapes_%s_13TeV.txt", channel.Data()); // the datacard name
//wsname = plotdir + plotname +".root"; // the workspace name
//useshapes = true;
//options1 = "unblind:SplusB";
//MakeDatacard(channel, dcname, wsname, useshapes, options1);
*/
delete ev;
}
