/////////////////////////////////////////////////////////////////////////////// // /////////////////////////////////////////////////////////////////////////////// double gX1 = 0.; double gX2 = 0.; double gTextFontSize = 1; //_____________________________________________________________________________ void plot_init() { //----------------------------------------------------------------------------- // step 1: make sure that all the needed scripts are loaded //----------------------------------------------------------------------------- char macro[200]; const char* script[] = { "plot_common.C", 0 }; const char* work_dir = gSystem->Getenv("SRT_LOCAL"); for (int i=0; script[i] != 0; i++) { sprintf(macro,"%s/Electroweak/scripts/%s",work_dir,script[i]); if (! gInterpreter->IsLoaded(macro)) gInterpreter->LoadMacro(macro); } } //_____________________________________________________________________________ void plot_w_width_direct(int mode = 0) { //----------------------------------------------------------------------------- // experimental results // // res +stat -stat +sys -sys lum Y title size value size color //----------------------------------------------------------------------------- Res_t sm_2004 = { 2.0921, 0.0025, 0.0025, 0. , 0. , 0., 0.25 , "SM prediction " , 0.04, "2.0921 #pm 0.0025", 0.04}; Res_t cdf_run1 = { 2.04 , 0.11 , 0.11 , 0.09, 0.09, 0., 0.85 , "CDF Run I " , 0.04, "2.05 #pm 0.10 #pm 0.08 ",0.04}; Res_t d0_run1 = { 2.23 , 0.15 , 0.14 , 0.1 , 0.1 , 0., 0.75 , "D0 Run I " , 0.04, "2.23 #pm^{0.15}_{0.14} #pm 0.1 ", 0.04}; Res_t tev_run1 = { 2.102 , 0.106 , 0.106 , 0.0 , 0.0 , 0., 0.65 , "Tevatron Run I " , 0.04, "2.102 #pm 0.106 " , 0.04}; Res_t d0_2004 = { 2.011 , 0.142 , 0.142 , 0. , 0. , 0., 0.40 , "D0 2004 prelim " , 0.04, "2.011 #pm 0.142 " , 0.04}; Res_t lep_2004 = { 2.15 , 0.091 , 0.091 , 0. , 0. , 0., 0.20 , "LEP prelim " , 0.04, "2.15 #pm 0.091 ", 0.04}; Res_t pdg_2004 = { 2.124 , 0.041 , 0.041 , 0. , 0. , 0., 0.15 , "PDG 2004 " , 0.04, "2.124 #pm 0.041 " , 0.04}; //----------------------------------------------------------------------------- gStyle->SetOptStat(0); plot_init(); //----------------------------------------------------------------------------- // canvas first //----------------------------------------------------------------------------- TCanvas* v = new TCanvas("w_width_direct","w_width_direct",0,0,680,440); //----------------------------------------------------------------------------- // draw box, axis and global title //----------------------------------------------------------------------------- TH2F* h = new TH2F("h_direct_width","",1,1.8,3.,1,0,1); h->SetNdivisions(502,"X"); h->SetNdivisions(500,"Y"); h->SetTickLength(0.01,"X"); h->SetTickLength(0.01,"Y"); h->GetXaxis()->SetLabelFont(42); h->GetXaxis()->SetTitle("G(W), GeV"); h->Draw(); TLatex* title = new TLatex(2.0,1.03,"Direct W width measurements"); title->SetTextFont(42); title->Draw(); //----------------------------------------------------------------------------- // defaults X's //----------------------------------------------------------------------------- gX1 = 2.45; gX2 = 2.7; gTextFontSize = 0.04; //----------------------------------------------------------------------------- // draw world average band //----------------------------------------------------------------------------- draw_band(sm_2004,5); //----------------------------------------------------------------------------- // experimental points //----------------------------------------------------------------------------- draw_point(cdf_run1,1); draw_point(d0_run1 ,1); draw_point(tev_run1,1); draw_point(lep_2004,4); draw_point(d0_2004 ,2); draw_point(pdg_2004,1); // if (mode != 0) draw_point(cdf_03); h->Draw("same"); } //_____________________________________________________________________________ void plot_w_width_indirect(int mode = 0) { //----------------------------------------------------------------------------- // experimental results //----------------------------------------------------------------------------- Res_t sm_2004 = { 2.0921, 0.0025, 0.0025, 0., 0., 0., 0.25 , "SM prediction " , 0.04, "2.0921 #pm 0.0025", 0.04}; Res_t cdf_r1a = { 2.16 , 0.17 , 0.17 , 0., 0., 0., 0.70 , "CDF Run IA(e+#mu)" , 0.04, "2.16 #pm 0.17 " , 0.04}; Res_t cdf_2004 = { 2.079 , 0.041 , 0.041 , 0., 0., 0., 0.60 , "CDF 2004 " , 0.04, "2.079 #pm 0.041" , 0.04}; Res_t pdg_2004 = { 2.124 , 0.041 , 0.041 , 0., 0., 0., 0.40 , "PDG 2004 " , 0.04, "2.124 #pm 0.041 ", 0.04}; gStyle->SetOptStat(0); plot_init(); //----------------------------------------------------------------------------- // canvas first //----------------------------------------------------------------------------- TCanvas* v = new TCanvas("w_indirect_width","w_indirect_width",0,0,700,600); //----------------------------------------------------------------------------- // draw box, axis and global title //----------------------------------------------------------------------------- TH2F* h = new TH2F("h_indirect_width","",1,1.5,3.,1,0,1); h->SetNdivisions(502,"X"); h->SetNdivisions(500,"Y"); h->SetTickLength(0.01,"X"); h->SetTickLength(0.01,"Y"); h->GetXaxis()->SetLabelFont(42); h->GetXaxis()->SetTitle("G(W), GeV"); h->Draw(); TLatex* title = new TLatex(2.,1.03,"W width measurements"); title->SetTextFont(42); title->Draw(); //----------------------------------------------------------------------------- // draw world average band //----------------------------------------------------------------------------- draw_band(sm_2004,5); //----------------------------------------------------------------------------- // experimental points //----------------------------------------------------------------------------- draw_point(cdf_2004,1); draw_point(pdg_2004,1); // if (mode != 0) draw_point(cdf_03,2); h->Draw("same"); } //_____________________________________________________________________________ void plot_sigma_w(int mode = 0) { gStyle->SetOptStat(0); plot_init(); //----------------------------------------------------------------------------- // canvas first, then draw box, axis and global title //----------------------------------------------------------------------------- TCanvas* v = new TCanvas("sigma_w","sigma_w",0,0,600,700); TH2F* h = new TH2F("h_sigma_w","",1,1000.,6000.,1,0,1); h->SetNdivisions(502,"X"); h->SetNdivisions(500,"Y"); h->SetTickLength(0.01,"X"); h->SetTickLength(0.01,"Y"); h->GetXaxis()->SetLabelFont(42); h->GetXaxis()->SetTitle("#sigma #times B, pb"); h->Draw(); TLatex* title = new TLatex(1150.,1.02,"Tevatron W #rightarrow l #nu cross section measurements"); title->SetTextFont(42); title->SetTextSize(0.04); title->Draw(); //----------------------------------------------------------------------------- // experimental results //----------------------------------------------------------------------------- Res_t thr_2004 = { 2687. , 54. , 54. , 0., 0., 0., 0.40 , "NNLO " , 0.03, "2687 #pm 54 pb ", 0.03}; Res_t cdf_wmnu_2004 = { 2768, 16, 16, 64, 64, 166, 0.75 , "CDF'04 (#mu)" , 0.029, "2768 #pm 16 #pm 64 #pm 166 pb", 0.029}; Res_t cdf_wenu_2004 = { 2780, 14, 14, 60, 60, 167, 0.70 , "CDF'04 (e)" , 0.029, "2780 #pm 14 #pm 60 #pm 167 pb", 0.029}; Res_t cdf_wlnu_2004 = { 2775, 10, 10, 53, 53, 167, 0.65 , "CDF'04 (e+#mu)" , 0.029, "2775 #pm 10 #pm 53 #pm 167 pb", 0.029}; Res_t cdf_wenu_2004_plug = { 2874, 34, 34, 167, 167, 172, 0.60 , "CDF'04 (e, plug)" , 0.029, "2874 #pm 34 #pm 167 #pm 172 pb", 0.029}; Res_t cdf_wenu_2004_plug_ = { -1.e10, 0,0,0,0, 0, 0.57 , "(preliminary)" , 0.02, "", 0.02}; Res_t cdf_wtnu_2003 = { 2620, 70, 70, 210, 210, 160, 0.53 , "CDF'03(#tau)" , 0.029, "2620 #pm 70 #pm 210 #pm 160 pb", 0.029}; Res_t cdf_wtnu_2003_ = { -1.e10, 0,0,0,0, 0, 0.50 , "(preliminary)" , 0.02, "", 0.02}; Res_t d0_wenu_2004 = { 2865.2, 8.3, 8.3, 74.7, 74.7, 186.2, 0.25 , "D0'04 (e)" , 0.029, "2865 #pm 8 #pm 75 #pm 186 pb", 0.029}; Res_t d0_wenu_2004_ = { -1.e10, 0,0,0,0, 0, 0.23 , "(preliminary)" , 0.02, "", 0.02}; Res_t d0_wmnu_2003 = { 3226, 128, 128, 100, 100, 322, 0.15 , "D0'03 (#mu)" , 0.029, "3226#pm128#pm100#pm322 pb", 0.029}; Res_t d0_wmnu_2003_ = { -1.e10, 8.3, 8.3, 74.7, 74.7, 186.2, 0.13 , "(preliminary)" , 0.02, "", 0.02}; gX1 = 1100; gX2 = 3800; //----------------------------------------------------------------------------- // draw world average band //----------------------------------------------------------------------------- draw_band(thr_2004,5); //----------------------------------------------------------------------------- // experimental points //----------------------------------------------------------------------------- draw_point(cdf_wmnu_2004,1); draw_point(cdf_wenu_2004,1); draw_point(cdf_wlnu_2004,4); draw_point(cdf_wenu_2004_plug,1); draw_point(cdf_wenu_2004_plug_,1); draw_point(cdf_wtnu_2003,1); draw_point(cdf_wtnu_2003_,1); draw_point(d0_wenu_2004,1); draw_point(d0_wenu_2004_,1); draw_point(d0_wmnu_2003,1); draw_point(d0_wmnu_2003_,1); h->Draw("same"); } //_____________________________________________________________________________ void plot_sigma_z(int mode = 0) { //----------------------------------------------------------------------------- // experimental results //----------------------------------------------------------------------------- Res_t thr_2004 = { 253.1 , 5. , 5. , 0., 0., 0., 0.40 , "NNLO " , 0.03, "253.1 #pm 0.005 ", 0.03}; Res_t cdf_zee_2004 = { 255.8, 3.9, 3.9, 5.5, 5.5, 15.4, 0.75 , "CDF'04 (e)" , 0.03, "255.8 #pm 3.9 #pm 5.5 #pm 15.4 pb", 0.028}; Res_t cdf_zmm_2004 = { 248.0, 5.9, 5.9, 7.6, 7.6, 14.9, 0.70 , "CDF'04 (#mu)" , 0.03, "248.0 #pm 5.9 #pm 7.6 #pm 14.9 pb", 0.028}; Res_t cdf_zll_2004 = { 254.9, 3.3, 3.3, 4.6, 4.6, 15.2, 0.65 , "CDF'04 (e+#mu)" , 0.03, "254.9 #pm 3.3 #pm 4.6 #pm 15.2 pb", 0.028}; Res_t cdf_ztt_2004 = { 242., 48.0, 48.0, 26., 26., 15, 0.60 , "CDF'04 (#tau)" , 0.03, "242 #pm 48.0 #pm 26 #pm 15 pb", 0.028}; Res_t cdf_ztt_2004_ = { -1.e10, 0,0,0,0, 0, 0.57 , " (preliminary)" , 0.02, "", 0.02}; Res_t d0_zee_2004 = { 264.9, 3.9, 3.9, 9.9, 9.9, 17.2, 0.40 , "D0'04 (e)" , 0.03, "264.9 #pm 3.9 #pm 9.9 #pm 17.2 pb", 0.028}; Res_t d0_zee_2004_ = { -1.e10, 0,0,0,0, 0, 0.37 , " (preliminary)" , 0.02, "", 0.02}; Res_t d0_zmm_2004 = { 291.3, 3.0, 3.0, 6.9, 6.9, 18.8, 0.30 , "D0'04 (#mu)" , 0.03, "291.3 #pm 3.0 #pm 6.9 #pm 18.8 pb", 0.028}; Res_t d0_zmm_2004_ = { -1.e10, 0,0,0,0, 0, 0.27 , " (preliminary)" , 0.02, "", 0.02}; Res_t d0_ztt_2004 = { 256., 16.0, 16.0, 17., 17., 16, 0.20 , "D0'04 (#tau)" , 0.03, "256 #pm 16 #pm 17 #pm 16 pb", 0.028}; Res_t d0_ztt_2004_ = { -1.e10, 0,0,0,0, 0, 0.17 , " (preliminary)" , 0.02, "", 0.02}; //----------------------------------------------------------------------------- gStyle->SetOptStat(0); plot_init(); //----------------------------------------------------------------------------- // canvas first //----------------------------------------------------------------------------- TCanvas* v = new TCanvas("sigma_z","sigma_z",0,0,600,700); //----------------------------------------------------------------------------- // draw box, axis and global title //----------------------------------------------------------------------------- TH2F* h = new TH2F("h_sigma_z","",1,0.,600.,1,0,1); h->SetNdivisions(502,"X"); h->SetNdivisions(500,"Y"); h->SetTickLength(0.01,"X"); h->SetTickLength(0.01,"Y"); h->GetXaxis()->SetLabelFont(42); h->GetXaxis()->SetTitle("#sigma #times B, pb"); h->Draw(); TLatex* title = new TLatex(30.,1.02,"Tevatron Z #rightarrow l^{+} l^{-} cross section measurements"); title->SetTextFont(42); title->SetTextSize(0.04); title->Draw(); gX1 = 10.; gX2 = 330; //----------------------------------------------------------------------------- // draw world average band //----------------------------------------------------------------------------- draw_band(thr_2004,5); //----------------------------------------------------------------------------- // experimental points //----------------------------------------------------------------------------- draw_point(cdf_zmm_2004,1); draw_point(cdf_zee_2004,1); draw_point(cdf_zll_2004,4); draw_point(cdf_ztt_2004 ,1); draw_point(cdf_ztt_2004_,1); draw_point(d0_zmm_2004 ,1); draw_point(d0_zmm_2004_,1); draw_point(d0_zee_2004 ,1); draw_point(d0_zee_2004_,1); draw_point(d0_ztt_2004 ,1); draw_point(d0_ztt_2004_,1); h->Draw("same"); } //_____________________________________________________________________________ void plot_sigma_ww(int mode = 0) { //----------------------------------------------------------------------------- // experimental results //----------------------------------------------------------------------------- Res_t mc_nlo_2004 = { 12.5 , 0.8 , .8 , 0., 0., 0., 0.40 , "NNLO " , 0.03, "253.1 #pm 0.005 ", 0.03 } ; Res_t cdf_ww_dl_2004 = { 14.3, 5.6, 4.9, 1.6, 1.6, 0.9, 0.80 , "CDF'04 (e)" , 0.03, "14.8 #pm ^{5.6}_{4.9} #pm 1.6 #pm 0.9 pb", 0.028}; Res_t cdf_ww_dl_2004_ = { -1.e10, 0,0,0,0, 0, 0.77 , "(preliminary)" , 0.02, "", 0.02}; Res_t cdf_ww_lt_2004 = { 19.4, 5.1, 5.1, 3.5, 3.5, 1.2, 0.70 , "CDF'04 (e)" , 0.03, "19.4 #pm 5.1 #pm 3.5 #pm 1.2 pb", 0.028}; Res_t cdf_ww_lt_2004_ = { -1.e10, 0,0,0,0, 0, 0.67 , "(preliminary)" , 0.02, "", 0.02}; Res_t d0_ww_ll_2004 = { 13.8, 4.3, 3.8, 1.0, 0.8, 0.9, 0.20 , "D0'04(combined)" , 0.03, "13.8 #pm ^{4.3}_{3.8} #pm ^{1.0}_{0.8} #pm 0.9 pb", 0.028}; Res_t d0_ww_ll_2004_ = { -1.e10, 0,0,0,0, 0, 0.17 , "(preliminary)" , 0.02, "", 0.02}; Res_t d0_ww_ee_2004 = { 14.6, 11.0, 8.4, 1.3, 1.2, 0.95, 0.27 , "D0'04(ee)" , 0.03, "14.6 #pm ^{11.0}_{8.4} #pm ^{1.3}_{1.2} #pm 0.95 pb", 0.028 }; Res_t d0_ww_ee_2004_ = { -1.e10, 0,0,0,0, 0, 0.24 , "(preliminary)" , 0.02, "", 0.02}; Res_t d0_ww_em_2004 = { 13.6, 5.1, 4.3, 1.1, 0.9, 0.9, 0.34 , "D0'04(e#mu)" , 0.03, "13.6 #pm ^{5.1}_{4.3} #pm ^{1.1}_{0.9} #pm 0.9 pb", 0.028}; Res_t d0_ww_em_2004_ = { -1.e10, 0,0,0,0, 0, 0.31 , "(preliminary)" , 0.02, "", 0.02}; Res_t d0_ww_mm_2004 = { 13.2, 14.8, 10.7, 1.8, 1.1, 0.9, 0.41 , "D0'04(#mu#mu)" , 0.03, "13.2 #pm ^{14.8}_{10.7} #pm ^{1.8}_{1.1} #pm 0.9 pb", 0.028}; Res_t d0_ww_mm_2004_ = { -1.e10, 0,0,0,0, 0, 0.38 , "(preliminary)" , 0.02, "", 0.02}; gStyle->SetOptStat(0); plot_init(); //----------------------------------------------------------------------------- // canvas first //----------------------------------------------------------------------------- TCanvas* v = new TCanvas("sigma_ww","sigma_ww",0,0,600,700); //----------------------------------------------------------------------------- // draw box, axis and global title //----------------------------------------------------------------------------- TH2F* h = new TH2F("h_sigma_ww","",1,-35.,65.,1,0,1); h->SetNdivisions(502,"X"); h->SetNdivisions(500,"Y"); h->SetTickLength(0.01,"X"); h->SetTickLength(0.01,"Y"); h->GetXaxis()->SetLabelFont(42); h->GetXaxis()->SetTitle("#sigma (p#bar{p} #rightarrow WW), pb"); h->Draw(); TLatex* title = new TLatex(-28.,1.02, "Tevatron p#bar{p} #rightarrow WW #rightarrow l^{+} l^{-} cross section measurements"); title->SetTextFont(42); title->SetTextSize(0.03); title->Draw(); //----------------------------------------------------------------------------- gX1 = -34.; gX2 = 30; //----------------------------------------------------------------------------- // draw MC@NLO band //----------------------------------------------------------------------------- draw_band(mc_nlo_2004,5); //----------------------------------------------------------------------------- // experimental points //----------------------------------------------------------------------------- draw_point(cdf_ww_dl_2004 ,1); draw_point(cdf_ww_dl_2004_,1); draw_point(cdf_ww_lt_2004 ,1); draw_point(cdf_ww_lt_2004_,1); draw_point(d0_ww_ee_2004 ,1); draw_point(d0_ww_ee_2004_,1); draw_point(d0_ww_em_2004 ,1); draw_point(d0_ww_em_2004_,1); draw_point(d0_ww_mm_2004 ,1); draw_point(d0_ww_mm_2004_,1); draw_point(d0_ww_ll_2004 ,4); draw_point(d0_ww_ll_2004_,4); h->Draw("same"); }