btagPerfPOGFormulas_nofit::btagPerfPOGFormulas_nofit(const char* inputfile) {
// just a sigmoid. Used to parametrize the efficiencies
std::cout<<"I'm using the hardcoded CSV b-tag SFs"<<std::endl;
esdata_ = TFile::Open(inputfile);
h_eff_csvl_b_brl_ = (TH1F*)esdata_->Get( "CSVL/h_eff_bTagOverGoodJet_ptb_Barrel" );
h_eff_csvl_b_fwd_ = (TH1F*)esdata_->Get( "CSVL/h_eff_bTagOverGoodJet_ptb_Endcaps" );
h_eff_csvl_c_brl_ = (TH1F*)esdata_->Get( "CSVL/h_eff_bTagOverGoodJet_ptc_Barrel" );
h_eff_csvl_c_fwd_ = (TH1F*)esdata_->Get( "CSVL/h_eff_bTagOverGoodJet_ptc_Endcaps" );
h_eff_csvl_l_brl_ = (TH1F*)esdata_->Get( "CSVL/h_eff_bTagOverGoodJet_ptl_Barrel" );
h_eff_csvl_l_fwd_ = (TH1F*)esdata_->Get( "CSVL/h_eff_bTagOverGoodJet_ptl_Endcaps" );
h_eff_csvm_b_brl_ = (TH1F*)esdata_->Get( "CSVM/h_eff_bTagOverGoodJet_ptb_Barrel" );
h_eff_csvm_b_fwd_ = (TH1F*)esdata_->Get( "CSVM/h_eff_bTagOverGoodJet_ptb_Endcaps" );
h_eff_csvm_c_brl_ = (TH1F*)esdata_->Get( "CSVM/h_eff_bTagOverGoodJet_ptc_Barrel" );
h_eff_csvm_c_fwd_ = (TH1F*)esdata_->Get( "CSVM/h_eff_bTagOverGoodJet_ptc_Endcaps" );
h_eff_csvm_l_brl_ = (TH1F*)esdata_->Get( "CSVM/h_eff_bTagOverGoodJet_ptl_Barrel" );
h_eff_csvm_l_fwd_ = (TH1F*)esdata_->Get( "CSVM/h_eff_bTagOverGoodJet_ptl_Endcaps" );
h_eff_csvt_b_brl_ = (TH1F*)esdata_->Get( "CSVT/h_eff_bTagOverGoodJet_ptb_Barrel" );
h_eff_csvt_b_fwd_ = (TH1F*)esdata_->Get( "CSVT/h_eff_bTagOverGoodJet_ptb_Endcaps" );
h_eff_csvt_c_brl_ = (TH1F*)esdata_->Get( "CSVT/h_eff_bTagOverGoodJet_ptc_Barrel" );
h_eff_csvt_c_fwd_ = (TH1F*)esdata_->Get( "CSVT/h_eff_bTagOverGoodJet_ptc_Endcaps" );
h_eff_csvt_l_brl_ = (TH1F*)esdata_->Get( "CSVT/h_eff_bTagOverGoodJet_ptl_Barrel" );
h_eff_csvt_l_fwd_ = (TH1F*)esdata_->Get( "CSVT/h_eff_bTagOverGoodJet_ptl_Endcaps" );
GetSFlmeanCSVL00 = GetSFlmean("CSV", "L", 0.0, 0.5, "ABCD");
GetSFlminCSVL00 = GetSFlmin("CSV", "L", 0.0, 0.5, "ABCD");
GetSFlmaxCSVL00 = GetSFlmax("CSV", "L", 0.0, 0.5, "ABCD");
GetSFlmeanCSVL05 = GetSFlmean("CSV", "L", 0.5, 1.0, "ABCD");
GetSFlminCSVL05 = GetSFlmin("CSV", "L", 0.5, 1.0, "ABCD");
GetSFlmaxCSVL05 = GetSFlmax("CSV", "L", 0.5, 1.0, "ABCD");
GetSFlmeanCSVL10 = GetSFlmean("CSV", "L", 1.0, 1.5, "ABCD");
GetSFlminCSVL10 = GetSFlmin("CSV", "L", 1.0, 1.5, "ABCD");
GetSFlmaxCSVL10 = GetSFlmax("CSV", "L", 1.0, 1.5, "ABCD");
GetSFlmeanCSVL15 = GetSFlmean("CSV", "L", 1.5, 2.4, "ABCD");
GetSFlminCSVL15 = GetSFlmin("CSV", "L", 1.5, 2.4, "ABCD");
GetSFlmaxCSVL15 = GetSFlmax("CSV", "L", 1.5, 2.4, "ABCD");
GetSFlmeanCSVM00 = GetSFlmean("CSV", "M", 0.0, 0.8, "ABCD");
GetSFlminCSVM00 = GetSFlmin("CSV", "M", 0.0, 0.8, "ABCD");
GetSFlmaxCSVM00 = GetSFlmax("CSV", "M", 0.0, 0.8, "ABCD");
GetSFlmeanCSVM08 = GetSFlmean("CSV", "M", 0.8, 1.6, "ABCD");
GetSFlminCSVM08 = GetSFlmin("CSV", "M", 0.8, 1.6, "ABCD");
GetSFlmaxCSVM08 = GetSFlmax("CSV", "M", 0.8, 1.6, "ABCD");
GetSFlmeanCSVM16 = GetSFlmean("CSV", "M", 1.6, 2.4, "ABCD");
GetSFlminCSVM16 = GetSFlmin("CSV", "M", 1.6, 2.4, "ABCD");
GetSFlmaxCSVM16 = GetSFlmax("CSV", "M", 1.6, 2.4, "ABCD");
GetSFlmeanCSVT00 = GetSFlmean("CSV", "T", 0.0, 2.4, "ABCD");
GetSFlminCSVT00 = GetSFlmin("CSV", "T", 0.0, 2.4, "ABCD");
GetSFlmaxCSVT00 = GetSFlmax("CSV", "T", 0.0, 2.4, "ABCD");
}
TF1* plotmean(TString tagger, TString TaggerStrength, float Etamin, float Etamax, TString opt = "" , int col = 1, float lineWidth = 1, int lineStyle = 1)
{
TF1* f = GetSFlmean(tagger,TaggerStrength,Etamin,Etamax);
if( f != NULL )
{
f->SetLineColor(col);
f->SetMinimum(0.4);
f->SetMaximum(1.6);
f->SetLineWidth(lineWidth);
f->SetLineStyle(lineStyle);
f->Draw(opt);
}
//else cout << "NULL pointer returned... Function seems not to exist" << endl;
return f;
}
