void ZbbMCFMgraphCT10(){ TFile* RootFile= new TFile("ZbbMCFMgraphCT10.root","CREATE"); int NbinedgesZPt = 6; int NbinedgesbbdR = 6; int NbinedgesbbM = 6; Float_t binedgesZPt[7] = {0.,25., 40., 55., 70., 95., 250. }; double axZPt[] = {0,0,0,0,0,0};double aexZPt[] = {0,0,0,0,0,0}; for(int i=0;i!=NbinedgesZPt;i++)axZPt[i]=0.5*binedgesZPt[i]+0.5*binedgesZPt[i+1]; Float_t binedgesbbdR[7] = {0.4, 1.15, 1.9, 2.4, 2.8, 3.2, 5.0 }; double axbbdR[] = {0,0,0,0,0,0};double aexbbdR[] = {0,0,0,0,0,0}; for(int i=0;i!=NbinedgesbbdR;i++)axbbdR[i]=0.5*binedgesbbdR[i]+0.5*binedgesbbdR[i+1]; Float_t binedgesbbM[7] = {10., 45., 65., 85., 115., 165, 350.}; double axbbM[] = {0,0,0,0,0,0};double aexbbM[] = {0,0,0,0,0,0}; for(int i=0;i!=NbinedgesbbM;i++)axbbM[i]=0.5*binedgesbbM[i]+0.5*binedgesbbM[i+1]; double ayZPt []={88.7463,92.6488,97.8508,77.5556,74.1817,69.9009}; double ayZPt_stat []={1.62048,3.29034,2.43893,1.81598,1.2376,0.865823}; double ayZPt_scale_up []={10.9597,14.9038,12.4894,12.0177,11.0146,12.3058}; double ayZPt_scale_down[]={7.17699,14.8259,11.2784,9.0269,9.0075,10.971}; double ayZPt_pdf_up []={2.13183,2.10539,2.17018,1.68261,1.59592,1.74342}; double ayZPt_pdf_down []={2.22468,2.19859,2.23978,1.75667,1.71806,1.88209}; double ayZPt_as_up []={2.76763,3.68487,1.84017,0.628658,0.195874,2.58713}; double ayZPt_as_down []={0.280963,0.0447232,1.78487,3.17411,0.195024,0.807662}; double aybbdR []={61.9975,76.0985,73.727,80.1596,98.232,111.626}; double aybbdR_stat []={1.36115,1.85372,1.89096,1.95081,2.54774,3.31728}; double aybbdR_scale_up []={9.564,8.9552,8.9407,18.1453,17.7108,15.637}; double aybbdR_scale_down[]={7.8039,9.57249,9.2397,5.87869,11.2326,17.3537}; double aybbdR_pdf_up []={0.927618,1.45526,1.53733,1.65204,2.17738,3.02964}; double aybbdR_pdf_down []={0.911554,1.48947,1.60371,1.71635,2.25379,3.34459}; double aybbdR_as_up []={0.699927,0.149583,0.197352,1.08186,2.11717,0.930498}; double aybbdR_as_down []={0.687098,3.08153,0.196504,1.05752,2.15915,4.92739}; double aybbM []={74.0819,73.1399,74.8065,86.1594,87.2404,90.71}; double aybbM_stat []={1.96512,2.9745,1.04721,1.55451,1.7063,3.58649}; double aybbM_scale_up []={10.5601,9.6872,10.0264,17.0619,11.3371,13.3118}; double aybbM_scale_down[]={9.09441,6.8225,9.0886,5.2576,19.6346,15.638}; double aybbM_pdf_up []={1.26459,1.57268,1.67106,1.9169,2.29729,2.51105}; double aybbM_pdf_down []={1.26535,1.60286,1.72272,1.99047,2.46724,2.85416}; double aybbM_as_up []={1.24689,2.02558,2.01832,0.125368,0.449048,0.69553}; double aybbM_as_down []={2.8363,1.93925,2.02858,0.125044,3.17571,0.686989}; for(int i=0;i!=NbinedgesZPt;i++)ayZPt[i] =ayZPt[i] /(-binedgesZPt[i]+binedgesZPt[i+1]); for(int i=0;i!=NbinedgesZPt;i++)ayZPt_stat[i] =ayZPt_stat[i] /(-binedgesZPt[i]+binedgesZPt[i+1]); for(int i=0;i!=NbinedgesZPt;i++)ayZPt_scale_up[i] =ayZPt_scale_up[i] /(-binedgesZPt[i]+binedgesZPt[i+1]); for(int i=0;i!=NbinedgesZPt;i++)ayZPt_scale_down[i]=ayZPt_scale_down[i]/(-binedgesZPt[i]+binedgesZPt[i+1]); for(int i=0;i!=NbinedgesZPt;i++)ayZPt_pdf_up[i] =ayZPt_pdf_up[i] /(-binedgesZPt[i]+binedgesZPt[i+1]); for(int i=0;i!=NbinedgesZPt;i++)ayZPt_pdf_down[i] =ayZPt_pdf_down[i] /(-binedgesZPt[i]+binedgesZPt[i+1]); for(int i=0;i!=NbinedgesZPt;i++)ayZPt_as_up[i] =ayZPt_as_up[i] /(-binedgesZPt[i]+binedgesZPt[i+1]); for(int i=0;i!=NbinedgesZPt;i++)ayZPt_as_down[i] =ayZPt_as_down[i] /(-binedgesZPt[i]+binedgesZPt[i+1]); for(int i=0;i!=NbinedgesbbdR;i++)aybbdR[i] =aybbdR[i] /(-binedgesbbdR[i]+binedgesbbdR[i+1]); for(int i=0;i!=NbinedgesbbdR;i++)aybbdR_stat[i] =aybbdR_stat[i] /(-binedgesbbdR[i]+binedgesbbdR[i+1]); for(int i=0;i!=NbinedgesbbdR;i++)aybbdR_scale_up[i] =aybbdR_scale_up[i] /(-binedgesbbdR[i]+binedgesbbdR[i+1]); for(int i=0;i!=NbinedgesbbdR;i++)aybbdR_scale_down[i]=aybbdR_scale_down[i]/(-binedgesbbdR[i]+binedgesbbdR[i+1]); for(int i=0;i!=NbinedgesbbdR;i++)aybbdR_pdf_up[i] =aybbdR_pdf_up[i] /(-binedgesbbdR[i]+binedgesbbdR[i+1]); for(int i=0;i!=NbinedgesbbdR;i++)aybbdR_pdf_down[i] =aybbdR_pdf_down[i] /(-binedgesbbdR[i]+binedgesbbdR[i+1]); for(int i=0;i!=NbinedgesbbdR;i++)aybbdR_as_up[i] =aybbdR_as_up[i] /(-binedgesbbdR[i]+binedgesbbdR[i+1]); for(int i=0;i!=NbinedgesbbdR;i++)aybbdR_as_down[i] =aybbdR_as_down[i] /(-binedgesbbdR[i]+binedgesbbdR[i+1]); for(int i=0;i!=NbinedgesbbM;i++)aybbM[i] =aybbM[i] /(-binedgesbbM[i]+binedgesbbM[i+1]); for(int i=0;i!=NbinedgesbbM;i++)aybbM_stat[i] =aybbM_stat[i] /(-binedgesbbM[i]+binedgesbbM[i+1]); for(int i=0;i!=NbinedgesbbM;i++)aybbM_scale_up[i] =aybbM_scale_up[i] /(-binedgesbbM[i]+binedgesbbM[i+1]); for(int i=0;i!=NbinedgesbbM;i++)aybbM_scale_down[i]=aybbM_scale_down[i]/(-binedgesbbM[i]+binedgesbbM[i+1]); for(int i=0;i!=NbinedgesbbM;i++)aybbM_pdf_up[i] =aybbM_pdf_up[i] /(-binedgesbbM[i]+binedgesbbM[i+1]); for(int i=0;i!=NbinedgesbbM;i++)aybbM_pdf_down[i] =aybbM_pdf_down[i] /(-binedgesbbM[i]+binedgesbbM[i+1]); for(int i=0;i!=NbinedgesbbM;i++)aybbM_as_up[i] =aybbM_as_up[i] /(-binedgesbbM[i]+binedgesbbM[i+1]); for(int i=0;i!=NbinedgesbbM;i++)aybbM_as_down[i] =aybbM_as_down[i] /(-binedgesbbM[i]+binedgesbbM[i+1]); TGraphAsymmErrors* ZPT_gae_stat = new TGraphAsymmErrors(NbinedgesZPt, axZPt, ayZPt, aexZPt, aexZPt, ayZPt_stat , ayZPt_stat ); ZPT_gae_stat ->Draw("AP"); TGraphAsymmErrors* ZPT_gae_scale = new TGraphAsymmErrors(NbinedgesZPt, axZPt, ayZPt, aexZPt, aexZPt, ayZPt_scale_down, ayZPt_scale_up); ZPT_gae_scale->Draw("AP"); TGraphAsymmErrors* ZPT_gae_pdf = new TGraphAsymmErrors(NbinedgesZPt, axZPt, ayZPt, aexZPt, aexZPt, ayZPt_pdf_down , ayZPt_pdf_up ); ZPT_gae_pdf ->Draw("AP"); TGraphAsymmErrors* ZPT_gae_as = new TGraphAsymmErrors(NbinedgesZPt, axZPt, ayZPt, aexZPt, aexZPt, ayZPt_as_down , ayZPt_as_up ); ZPT_gae_as ->Draw("AP"); TGraphAsymmErrors* bbdR_gae_stat = new TGraphAsymmErrors(NbinedgesbbdR, axbbdR, aybbdR, aexbbdR, aexbbdR, aybbdR_stat , aybbdR_stat ); bbdR_gae_stat ->Draw("AP"); TGraphAsymmErrors* bbdR_gae_scale = new TGraphAsymmErrors(NbinedgesbbdR, axbbdR, aybbdR, aexbbdR, aexbbdR, aybbdR_scale_down, aybbdR_scale_up); bbdR_gae_scale->Draw("AP"); TGraphAsymmErrors* bbdR_gae_pdf = new TGraphAsymmErrors(NbinedgesbbdR, axbbdR, aybbdR, aexbbdR, aexbbdR, aybbdR_pdf_down , aybbdR_pdf_up ); bbdR_gae_pdf ->Draw("AP"); TGraphAsymmErrors* bbdR_gae_as = new TGraphAsymmErrors(NbinedgesbbdR, axbbdR, aybbdR, aexbbdR, aexbbdR, aybbdR_as_down , aybbdR_as_up ); bbdR_gae_as ->Draw("AP"); TGraphAsymmErrors* bbM_gae_stat = new TGraphAsymmErrors(NbinedgesbbM, axbbM, aybbM, aexbbM, aexbbM, aybbM_stat , aybbM_stat ); bbM_gae_stat ->Draw("AP"); TGraphAsymmErrors* bbM_gae_scale = new TGraphAsymmErrors(NbinedgesbbM, axbbM, aybbM, aexbbM, aexbbM, aybbM_scale_down, aybbM_scale_up); bbM_gae_scale->Draw("AP"); TGraphAsymmErrors* bbM_gae_pdf = new TGraphAsymmErrors(NbinedgesbbM, axbbM, aybbM, aexbbM, aexbbM, aybbM_pdf_down , aybbM_pdf_up ); bbM_gae_pdf ->Draw("AP"); TGraphAsymmErrors* bbM_gae_as = new TGraphAsymmErrors(NbinedgesbbM, axbbM, aybbM, aexbbM, aexbbM, aybbM_as_down , aybbM_as_up ); bbM_gae_as ->Draw("AP"); ZPT_gae_stat ->SetName("ZPT_gae_stat"); ZPT_gae_scale ->SetName("ZPT_gae_scale"); ZPT_gae_pdf ->SetName("ZPT_gae_pdf"); ZPT_gae_as ->SetName("ZPT_gae_as"); bbdR_gae_stat ->SetName("bbdR_gae_stat"); bbdR_gae_scale ->SetName("bbdR_gae_scale"); bbdR_gae_pdf ->SetName("bbdR_gae_pdf"); bbdR_gae_as ->SetName("bbdR_gae_as"); bbM_gae_stat ->SetName("bbM_gae_stat"); bbM_gae_scale ->SetName("bbM_gae_scale"); bbM_gae_pdf ->SetName("bbM_gae_pdf"); bbM_gae_as ->SetName("bbM_gae_as"); ZPT_gae_stat ->Write(); ZPT_gae_scale ->Write(); ZPT_gae_pdf ->Write(); ZPT_gae_as ->Write(); bbdR_gae_stat ->Write(); bbdR_gae_scale ->Write(); bbdR_gae_pdf ->Write(); bbdR_gae_as ->Write(); bbM_gae_stat ->Write(); bbM_gae_scale ->Write(); bbM_gae_pdf ->Write(); bbM_gae_as ->Write(); RootFile->Close(); return ; }