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; }