void L1UCTEfficiency::Loop(const char * var, const char * binning, const char *name, const char *opt)
{
if (fChain == 0) return;
Long64_t nentries = fChain->GetEntriesFast();
std::cout << "L1UCTEfficiency> data nentries: " << nentries << std::endl;
m_histograms1D["denominator"] = MakePlot( var , m_DSel.str().c_str(), binning, "Var", name,
(TString("hd") + TString(opt)).Data() );
//m_histograms1D["denominator"]->Draw();
m_histograms1D["l1num"] = MakePlot( var , m_l1Sel.str().c_str(),binning , "Var",name,
(TString("l1num") + TString(opt)).Data() );
//m_histograms1D["l1num"]->Draw();
m_histograms1D["l1gnum"] = MakePlot( var, m_l1gSel.str().c_str(),binning, "Var",name,
(TString("l1gnum") + TString(opt)).Data() );
//m_histograms1D["l1gnum"]->Draw();
MakeL1Efficiency( m_histograms1D["l1num"], m_histograms1D["denominator"], "l1Eff");
MakeL1GEfficiency( m_histograms1D["l1gnum"], m_histograms1D["denominator"], "l1gEff");
m_nbins = m_histograms1D["denominator"]->GetXaxis()->GetNbins();
m_xmin = m_histograms1D["denominator"]->GetXaxis()->GetXmin();
m_xmax = m_histograms1D["denominator"]->GetXaxis()->GetXmax();
}
void draw_scatter_plot() {
MakePlot("mht30", 150, 0., 750., "true_mu_ptW", 150, 0., 750.,"p_{T}^{W} [GeV];H_{T}^{miss} [GeV]","true_mu_pt>0&&num_true_mus==1&&num_true_els+num_true_had_taus==0", "reduced_trees/13TeV/TTJets_MSDecaysCKM_central_Tune4C_13TeV-madgraph-tauola_Phys14DR-PU20bx25_PHYS14_25_V1-v1_MINIAODSIM_UCSB2406_v78/*.root", "ttbar_1mu", true);
MakePlot("ht30", 150, 0., 1500., "true_mu_ptW", 150, 0., 750.,"p_{T}^{W} [GeV];H_{T} [GeV]","true_mu_pt>0&&num_true_mus==1&&num_true_els+num_true_had_taus==0", "reduced_trees/13TeV/TTJets_MSDecaysCKM_central_Tune4C_13TeV-madgraph-tauola_Phys14DR-PU20bx25_PHYS14_25_V1-v1_MINIAODSIM_UCSB2406_v78/*.root", "ttbar_1mu", true);
MakePlot("true_mu_pt", 150, 0., 750., "true_mu_ptW", 150, 0., 750.,"p_{T}^{W} [GeV];p_{T}^{#mu} [GeV]","true_mu_pt>0&&num_true_mus==1&&num_true_els+num_true_had_taus==0", "reduced_trees/13TeV/TTJets_MSDecaysCKM_central_Tune4C_13TeV-madgraph-tauola_Phys14DR-PU20bx25_PHYS14_25_V1-v1_MINIAODSIM_UCSB2406_v78/*.root", "ttbar_1mu", true);
MakePlot("mht30", 150, 0., 750., "0.5*(1-cos(true_mu_dTT))", 150, 0., 1.,"#frac{1}{2}(1-cos#Delta#theta_{T});H_{T}^{miss} [GeV]","true_mu_pt>0&&num_true_mus==1&&num_true_els+num_true_had_taus==0", "reduced_trees/13TeV/TTJets_MSDecaysCKM_central_Tune4C_13TeV-madgraph-tauola_Phys14DR-PU20bx25_PHYS14_25_V1-v1_MINIAODSIM_UCSB2406_v78/*.root", "ttbar_1mu", true, "true_mu_dTT");
MakePlot("true_mu_ptW", 150, 0., 750., "0.5*(1-cos(true_mu_dTT))", 150, 0., 1.,"#frac{1}{2}(1-cos#Delta#theta_{T});p_{T}^{W} [GeV]","true_mu_pt>0&&num_true_mus==1&&num_true_els+num_true_had_taus==0", "reduced_trees/13TeV/TTJets_MSDecaysCKM_central_Tune4C_13TeV-madgraph-tauola_Phys14DR-PU20bx25_PHYS14_25_V1-v1_MINIAODSIM_UCSB2406_v78/*.root", "ttbar_1mu", true, "true_mu_dTT");
MakePlot("true_mu_pt", 100, 0., 200., "0.5*(1-cos(true_mu_dTT))", 150, 0., 1.,"#frac{1}{2}(1-cos#Delta#theta_{T});p_{T}^{#mu} [GeV]","true_mu_pt>0&&num_true_mus==1&&num_true_els+num_true_had_taus==0", "reduced_trees/13TeV/TTJets_MSDecaysCKM_central_Tune4C_13TeV-madgraph-tauola_Phys14DR-PU20bx25_PHYS14_25_V1-v1_MINIAODSIM_UCSB2406_v78/*.root", "ttbar_1mu", true, "true_mu_dTT");
}
void UCTRateTree::Loop(const char * var, const char * selec, const char * binning, const char *name)
{
if (fChain == 0) return;
m_debug = true;
Long64_t nentries = fChain->GetEntriesFast();
std::cout << "UCTRateTree> data nentries: " << nentries << std::endl;
m_histograms1D[name] = MakePlot( var, selec, binning, "pt", "L1 trigger studies",name);
m_histograms1D[name]->Scale( L1_CALIB_FACTOR );
TString l1gname(name);
m_histograms1D[l1gname.Data()] = Make_Th_Integral( m_histograms1D[name] );
}
