double SigFun(double *x, double *par) { double val=par[0]*log(10.0)*exp(-x[0]*log(10.0)); return val; } //===================================================================== //-------------- Main script ------------------------------------------ //___________________________________________________________________________ // List of plots: // 0:"MetAll"; 1:"MetSig"; 2:"Met"; 3:"Et1"; 4:"Et2"; 5:"Etem"; 6:"Etjet"; 7:"Ht"; 8:"M"; 9:"Mej"; // 10:"Mextra"; 11:"Mjj"; 12:"Nem"; 13:"Njet15"; 14:"Qt"; 15:"dPhi"; 16:"dPhi1"; 17:"dPhi2"; 18:"dPhi3"; // 19:"Njet20"; 20:"Njet25"; 21:"dPhi4"; 22:"dPhi5; 23:"DeltaR" void myAnalysisPlots_gmsbpre(int plot_id=0, int rebincode=1) { gROOT->Reset(); gStyle->SetOptFit(1111); c1 = new TCanvas("c1","Met Model Study",200,10,700,500); c1->GetFrame()->SetBorderSize(12); char fName_wvx[200]; char fName_wvx1[200]; char fName_tri[200]; char fName_tri1[200]; char fName[200]; char fName1[200]; char fObject[200]; char fPlotTitle[200]; char fPlotTitle1[200]; char fPlotTitle2[200]; char fName_GMSB[200]; float scale_F[14]; scale_F[0]=0.0635; // Bauer W(e)+pho scale_F[1]=0.0605; // Bauer W(mu)+pho scale_F[2]=0.0609; // Bauer W(tau)+pho scale_F[3]=0.00393; // Bauer Z(e)+pho scale_F[4]=0.00395; // Bauer Z(mu)+pho scale_F[5]=0.00396; // Bauer Z(tau)+pho scale_F[6]=0.0387; // Pythia Z(tau) scale_F[7]=0.056; // Pythia Z(e) scale_F[8]=0.0876; // Pythia Z(mu) scale_F[9]=0.00234; // Pythia ttbar scale_F[10]=0.210; // Pythia W(e) scale_F[11]=0.296; // Pythia W(tau) scale_F[12]=0.308; // Pythia W(mu) scale_F[13]=0.001035; // Pythia W(mu) char fName_gg[14][200]; double qcd_scale=0.1342; double qcd_scale_er=0.00075; double h_int=0.0; double f_par=0.1; //________ scale for e+pho double w_scale_sig; double w_scale_sig_err; //________ scale factor for B.H. double bh_scale; double bh_scale_er; //________ scale factor for C.R. double cr_scale; double cr_scale_er; //-------- scale factor for GMSB siganl double sig_scale; double sig_scale_er; bh_scale=0.0082; bh_scale_er=0.0082; //________ scale factor for C.R. cr_scale=0.00063; cr_scale_er=0.00063; //________ scale for e+pho (with sig>5.0 cut) // w_scale_sig=0.748; // w_scale_sig_err=0.047; w_scale_sig=1.0; w_scale_sig_err=0.23; //-------- scale for GMSB signal sig_scale=0.02297*2591/133330; sig_scale_er=0.1*sig_scale; switch(plot_id) { case 0: sprintf(fPlotTitle,"MetAll"); sprintf(fPlotTitle1,"WphoMetAll_srv"); sprintf(fPlotTitle2,"MetGenAll"); break; case 1: sprintf(fPlotTitle,"MetSig"); sprintf(fPlotTitle1,"WphoMetSig_srv"); sprintf(fPlotTitle2,"MetGenSig"); break; case 2: sprintf(fPlotTitle,"Met"); sprintf(fPlotTitle1,"WphoMet_srv"); break; case 3: sprintf(fPlotTitle,"Et1"); sprintf(fPlotTitle1,"WphoEt1_srv"); break; case 4: sprintf(fPlotTitle,"Et2"); sprintf(fPlotTitle1,"WphoEt2_srv"); break; case 5: sprintf(fPlotTitle,"Etem"); sprintf(fPlotTitle1,"WphoEtem_srv"); break; case 6: sprintf(fPlotTitle,"Etjet"); sprintf(fPlotTitle1,"WphoEtjet_srv"); break; case 7: sprintf(fPlotTitle,"Ht"); sprintf(fPlotTitle1,"WphoHt_srv"); break; case 8: sprintf(fPlotTitle,"M"); sprintf(fPlotTitle1,"WphoM_srv"); break; case 9: sprintf(fPlotTitle,"Mej"); sprintf(fPlotTitle1,"WphoMej_srv"); break; case 10: sprintf(fPlotTitle,"Mextra"); sprintf(fPlotTitle1,"WphoMextra_srv"); break; case 11: sprintf(fPlotTitle,"Mjj"); sprintf(fPlotTitle1,"WphoMjj_srv"); break; case 12: sprintf(fPlotTitle,"Nem"); sprintf(fPlotTitle1,"WphoNem_srv"); break; case 13: sprintf(fPlotTitle,"Njet15"); sprintf(fPlotTitle1,"WphoNjet15_srv"); break; case 14: sprintf(fPlotTitle,"Qt"); sprintf(fPlotTitle1,"WphoQt_srv"); break; case 15: sprintf(fPlotTitle,"dPhi"); sprintf(fPlotTitle1,"WphoDelPhi_srv"); break; case 16: sprintf(fPlotTitle,"dPhi1"); sprintf(fPlotTitle1,"WphoDelPhi1_srv"); break; case 17: sprintf(fPlotTitle,"dPhi2"); sprintf(fPlotTitle1,"WphoDelPhi2_srv"); break; case 18: sprintf(fPlotTitle,"dPhi3"); sprintf(fPlotTitle1,"WphoDelPhi3_srv"); break; case 19: sprintf(fPlotTitle,"Njet20"); sprintf(fPlotTitle1,"WphoNjet20_srv"); break; case 20: sprintf(fPlotTitle,"Njet25"); sprintf(fPlotTitle1,"WphoNjet25_srv"); break; case 21: sprintf(fPlotTitle,"dPhi4"); sprintf(fPlotTitle1,"WphoDelPhi4_srv"); break; case 22: sprintf(fPlotTitle,"dPhi5"); sprintf(fPlotTitle1,"WphoDelPhi5_srv"); break; case 23: sprintf(fPlotTitle,"DeltaR"); sprintf(fPlotTitle1,"DeltaR_srv"); break; default: sprintf(fPlotTitle,"MetAll"); sprintf(fPlotTitle1,"WphoMetAll_srv"); } //------------ di-pho sample sprintf(fName,"qcd_data/myhisto_QCDdata_Ht000_dphi0.00_sig0_new.root"); // data sprintf(fName1,"qcd_data/myhisto_QCDdata_Ht000_dphi0.00_sig0_new.root"); // data //___________ MC di-pho files sprintf(fName_gg[0],"ewk_mc/myhisto_Bauer_Wepho_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[1],"ewk_mc/myhisto_Bauer_Wmupho_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[2],"ewk_mc/myhisto_Bauer_Wtaupho_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[3],"ewk_mc/myhisto_Bauer_Zeepho_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[4],"ewk_mc/myhisto_Bauer_Zmumupho_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[5],"ewk_mc/myhisto_Bauer_Ztautaupho_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[6],"ewk_mc/myhisto_Pythia_Ztautau_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[7],"ewk_mc/myhisto_Pythia_Zee_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[8],"ewk_mc/myhisto_Pythia_Zmumu_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[9],"ewk_mc/myhisto_Pythia_ttbarIncl_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[10],"ewk_mc/myhisto_Pythia_We_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[11],"ewk_mc/myhisto_Pythia_Wtau_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[12],"ewk_mc/myhisto_Pythia_Wmu_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_gg[13],"ewk_mc/myhisto_Pythia_Znunupho_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_wvx,"wrongvx/myhisto_wrongvx_Ht000_dphi0.00_sig0_new.root"); sprintf(fName_tri,"tripho/myhisto_tripho_Ht000_dphi0.00_sig0_new.root"); //------------ GMSB signal sprintf(fName_GMSB,"signal/myhisto_140_0_Ht000_dphi0.00_sig0_new.root"); TFile* myFile1= new TFile(fName); //connect file, everything stays in memory till "delete" TFolder* myFolder1 = (TFolder* ) myFile1->Get("Ana;1"); // get folder TFile* myFile11= new TFile(fName1); //connect file, everything stays in memory till "delete" TFolder* myFolder11 = (TFolder* ) myFile1->Get("Ana;1"); //------------ Beam Halo sprintf(fName,"beamhalo/myhisto_BH_Ht000_dphi0.00_sig0_new.root"); TFile* myFile4= new TFile(fName); //connect file, everything stays in memory till "delete" TFolder* myFolder4 = (TFolder* ) myFile4->Get("Ana;1"); // get folder //------------ Cosmics sprintf(fName,"cosmics/myhisto_CR_Ht000_dphi0.00_sig0_new.root"); TFile* myFile5= new TFile(fName); //connect file, everything stays in memory till "delete" TFolder* myFolder5 = (TFolder* ) myFile5->Get("Ana;1"); // get folder //----------- GMSB signal TFile* myFile7=new TFile(fName_GMSB); TFolder* myFolder7 = (TFolder* ) myFile7->Get("Ana;1"); //---wrong vx TFile* myFile9=new TFile(fName_wvx); TFolder* myFolder9 = (TFolder* ) myFile9->Get("Ana;1"); //---wrong tripho TFile* myFile10=new TFile(fName_tri); TFolder* myFolder10 = (TFolder* ) myFile10->Get("Ana;1"); //___________ MC di-pho files TFile* myFile_gg0= new TFile(fName_gg[0]); TFile* myFile_gg1= new TFile(fName_gg[1]); TFile* myFile_gg2= new TFile(fName_gg[2]); TFile* myFile_gg3= new TFile(fName_gg[3]); TFile* myFile_gg4= new TFile(fName_gg[4]); TFile* myFile_gg5= new TFile(fName_gg[5]); TFile* myFile_gg6= new TFile(fName_gg[6]); TFile* myFile_gg7= new TFile(fName_gg[7]); TFile* myFile_gg8= new TFile(fName_gg[8]); TFile* myFile_gg9= new TFile(fName_gg[9]); TFile* myFile_gg10= new TFile(fName_gg[10]); TFile* myFile_gg11= new TFile(fName_gg[11]); TFile* myFile_gg12= new TFile(fName_gg[12]); TFile* myFile_gg13= new TFile(fName_gg[13]); TFolder* myFolder_gg0 = (TFolder* ) myFile_gg0->Get("Ana;1"); TFolder* myFolder_gg1 = (TFolder* ) myFile_gg1->Get("Ana;1"); TFolder* myFolder_gg2 = (TFolder* ) myFile_gg2->Get("Ana;1"); TFolder* myFolder_gg3 = (TFolder* ) myFile_gg3->Get("Ana;1"); TFolder* myFolder_gg4 = (TFolder* ) myFile_gg4->Get("Ana;1"); TFolder* myFolder_gg5 = (TFolder* ) myFile_gg5->Get("Ana;1"); TFolder* myFolder_gg6 = (TFolder* ) myFile_gg6->Get("Ana;1"); TFolder* myFolder_gg7 = (TFolder* ) myFile_gg7->Get("Ana;1"); TFolder* myFolder_gg8 = (TFolder* ) myFile_gg8->Get("Ana;1"); TFolder* myFolder_gg9 = (TFolder* ) myFile_gg9->Get("Ana;1"); TFolder* myFolder_gg10 = (TFolder* ) myFile_gg10->Get("Ana;1"); TFolder* myFolder_gg11 = (TFolder* ) myFile_gg11->Get("Ana;1"); TFolder* myFolder_gg12 = (TFolder* ) myFile_gg12->Get("Ana;1"); TFolder* myFolder_gg13 = (TFolder* ) myFile_gg13->Get("Ana;1"); TH1F *h_gg0; TH1F *h_gg1; TH1F *h_gg2; TH1F *h_gg3; TH1F *h_gg4; TH1F *h_gg5; TH1F *h_gg6; TH1F *h_gg7; TH1F *h_gg8; TH1F *h_gg9; TH1F *h_gg10; TH1F *h_gg11; TH1F *h_gg12; TH1F *h_gg13; //_________ getting MC histo for di-pho sprintf(fObject,"MyJetFilter/Hist/Ana_data/%s",fPlotTitle); // get object name h_gg0= (TH1F*) myFolder_gg0->FindObject(fObject); h_gg1= (TH1F*) myFolder_gg1->FindObject(fObject); h_gg2= (TH1F*) myFolder_gg2->FindObject(fObject); h_gg3= (TH1F*) myFolder_gg3->FindObject(fObject); h_gg4= (TH1F*) myFolder_gg4->FindObject(fObject); h_gg5= (TH1F*) myFolder_gg5->FindObject(fObject); h_gg6= (TH1F*) myFolder_gg6->FindObject(fObject); h_gg7= (TH1F*) myFolder_gg7->FindObject(fObject); h_gg8= (TH1F*) myFolder_gg8->FindObject(fObject); h_gg9= (TH1F*) myFolder_gg9->FindObject(fObject); h_gg10= (TH1F*) myFolder_gg10->FindObject(fObject); h_gg11= (TH1F*) myFolder_gg11->FindObject(fObject); h_gg12= (TH1F*) myFolder_gg12->FindObject(fObject); h_gg13= (TH1F*) myFolder_gg13->FindObject(fObject); //------------------------------------------------------------------------ //--------- Accessing data histograms TH1F *h_data; TH1F *h_fake; TH1F *h_bh; TH1F *h_cr; TH1F *h_sig; //GMSB signal TH1F *h_wvx; TH1F *h_wvx1; TH1F *h_tri; TH1F *h_tri1; sprintf(fObject,"MyJetFilter/Hist/Ana_data/%s",fPlotTitle); h_sig= (TH1F*) myFolder7->FindObject(fObject); sprintf(fObject,"MyJetFilter/Hist/Ana_data/%s",fPlotTitle); // get object name h_data= (TH1F*) myFolder1->FindObject(fObject); // get histogram sprintf(fObject,"MyJetFilter/Hist/Ana_bckg/%s",fPlotTitle); // get object name h_fake= (TH1F*) myFolder1->FindObject(fObject); // get histogram sprintf(fObject,"JetModule2/Hist/Ana_data/%s",fPlotTitle); h_wvx= (TH1F*) myFolder9->FindObject(fObject); sprintf(fObject,"JetModule2/Hist/Ana_bckg/%s",fPlotTitle); h_wvx1= (TH1F*) myFolder9->FindObject(fObject); sprintf(fObject,"MyJetFilter/Hist/Ana_data/%s",fPlotTitle); h_tri= (TH1F*) myFolder10->FindObject(fObject); sprintf(fObject,"MyJetFilter/Hist/Ana_bckg/%s",fPlotTitle); // get object name h_tri1= (TH1F*) myFolder10->FindObject(fObject); sprintf(fObject,"MyJetFilter/Hist/Ana_data/%s",fPlotTitle); h_bh= (TH1F*) myFolder4->FindObject(fObject); // get histogram sprintf(fObject,"MyJetFilter/Hist/Ana_data/%s",fPlotTitle); // get object name h_cr= (TH1F*) myFolder5->FindObject(fObject); // get histogram if(rebincode>1) { h_gg0->Rebin(rebincode); h_gg1->Rebin(rebincode); h_gg2->Rebin(rebincode); h_gg3->Rebin(rebincode); h_gg4->Rebin(rebincode); h_gg5->Rebin(rebincode); h_gg6->Rebin(rebincode); h_gg7->Rebin(rebincode); h_gg8->Rebin(rebincode); h_gg9->Rebin(rebincode); h_gg10->Rebin(rebincode); h_gg11->Rebin(rebincode); h_gg12->Rebin(rebincode); h_gg13->Rebin(rebincode); h_data->Rebin(rebincode); h_fake->Rebin(rebincode); h_bh->Rebin(rebincode); h_cr->Rebin(rebincode); h_sig->Rebin(rebincode); h_wvx->Rebin(rebincode); h_wvx1->Rebin(rebincode); h_tri->Rebin(rebincode); h_tri1->Rebin(rebincode); } //* I need to do this in order to clean memory from ROOT prompt TFile* mytmpfile= new TFile("tmp.tmp","recreate"); // create auxiliary file mytmpfile->cd(); //-------------------------------------- Met histograms int nbins=h_data->GetNbinsX(); float binwidth=h_data->GetBinWidth(1); float lowedge=(h_data->GetBinCenter(1))-binwidth/2.0; float highedge=lowedge+nbins*binwidth; float mbins=h_sig->GetNbinsX(); float binwidth1=h_sig->GetBinWidth(1); float lowedge1=(h_sig->GetBinCenter(1))-binwidth/2.0; float highedge1=lowedge1+mbins*binwidth; Char_t title[200]; Char_t name[200]; sprintf(name,"h"); //sprintf(title,"CDF Run II Preliminary, 2.0 fb^{-1}: #gamma#gamma signal sample"); TH1F* myFnl_dt= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_gn= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_tot= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_tot_mG= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_tot_pG= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_bh= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_cr= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_noncol= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_sig= new TH1F(name,title,mbins,lowedge1,highedge1); TH1F* h_MCgg= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_MCgg= new TH1F(name,title,nbins,lowedge,highedge); TH1F* myFnl_wvx=new TH1F(name, title, nbins,lowedge,highedge); TH1F* myFnl_wvx1=new TH1F(name, title, nbins,lowedge,highedge); TH1F* myFnl_tri=new TH1F(name, title, nbins,lowedge,highedge); TH1F* myFnl_tri1=new TH1F(name, title, nbins,lowedge,highedge); myFnl_MCgg->Sumw2(); h_MCgg->Add(h_gg0,h_gg1,scale_F[0],scale_F[1]); h_MCgg->Add(h_gg2,scale_F[2]); h_MCgg->Add(h_gg3,scale_F[3]); h_MCgg->Add(h_gg4,scale_F[4]); h_MCgg->Add(h_gg5,scale_F[5]); h_MCgg->Add(h_gg6,scale_F[6]); h_MCgg->Add(h_gg7,scale_F[7]); h_MCgg->Add(h_gg8,scale_F[8]); h_MCgg->Add(h_gg9,scale_F[9]); h_MCgg->Add(h_gg10,scale_F[10]); h_MCgg->Add(h_gg11,scale_F[11]); h_MCgg->Add(h_gg12,scale_F[12]); h_MCgg->Add(h_gg13,scale_F[13]); myFnl_wvx->Sumw2(); myFnl_wvx1->Sumw2(); myFnl_tri->Sumw2(); myFnl_tri1->Sumw2(); myFnl_dt->Sumw2(); myFnl_gn->Sumw2(); myFnl_tot->Sumw2(); myFnl_bh->Sumw2(); myFnl_cr->Sumw2(); myFnl_noncol->Sumw2(); myFnl_sig->Sumw2(); //GMSB signal float binX; float bin_value; float bin_value1; float bin_error; for(int i=0; iGetBinContent(i+1); bin_error=h_MCgg->GetBinError(i+1); bin_error=(bin_error*bin_error*w_scale_sig*w_scale_sig+bin_value*bin_value*w_scale_sig_err*w_scale_sig_err); myFnl_MCgg->SetBinContent(i+1,bin_value*w_scale_sig); myFnl_MCgg->SetBinError(i+1,sqrt(bin_error)); bin_value=h_fake->GetBinContent(i+1); bin_error=h_fake->GetBinError(i+1); myFnl_gn->SetBinContent(i+1,bin_value); myFnl_gn->SetBinError(i+1,bin_error); bin_value=h_data->GetBinContent(i+1); bin_error=h_data->GetBinError(i+1); myFnl_dt->SetBinContent(i+1,bin_value); myFnl_dt->SetBinError(i+1,bin_error); h_int=h_int+bin_value; bin_value=h_bh->GetBinContent(i+1); bin_error=h_bh->GetBinError(i+1); myFnl_bh->SetBinContent(i+1,bin_value*bh_scale); myFnl_bh->SetBinError(i+1,bin_error*bh_scale_er); bin_value=h_cr->GetBinContent(i+1); bin_error=h_cr->GetBinError(i+1); myFnl_cr->SetBinContent(i+1,bin_value*cr_scale); myFnl_cr->SetBinError(i+1,bin_error*cr_scale_er); bin_value=h_wvx->GetBinContent(i+1); bin_error=h_wvx->GetBinError(i+1); myFnl_wvx->SetBinContent(i+1,bin_value*qcd_scale); myFnl_wvx->SetBinError(i+1,bin_error*qcd_scale_er); bin_value=h_wvx1->GetBinContent(i+1); bin_error=h_wvx1->GetBinError(i+1); myFnl_wvx1->SetBinContent(i+1,bin_value*qcd_scale); myFnl_wvx1->SetBinError(i+1,bin_error*qcd_scale_er); bin_value=h_tri->GetBinContent(i+1); bin_error=h_tri->GetBinError(i+1); myFnl_tri->SetBinContent(i+1,bin_value*qcd_scale); myFnl_tri->SetBinError(i+1,bin_error*qcd_scale_er); bin_value=h_tri1->GetBinContent(i+1); bin_error=h_tri1->GetBinError(i+1); myFnl_tri1->SetBinContent(i+1,bin_value*qcd_scale); myFnl_tri1->SetBinError(i+1,bin_error*qcd_scale_er); //----- GMSB part bin_value=h_sig->GetBinContent(i+1); bin_error=h_sig->GetBinError(i+1); myFnl_sig->SetBinContent(i+1,bin_value*sig_scale); myFnl_sig->SetBinError(i+1,bin_error*sig_scale_er); } myFnl_noncol->Add(myFnl_bh,myFnl_cr,1.0,1.0); // myFnl_tot->Add(myFnl_gn,myFnl_MCgg,1.0,1.0); myFnl_tot->Add(myFnl_MCgg,1.0); myFnl_tot->Add(myFnl_gn,1.0); myFnl_tot->Add(myFnl_wvx,1.0); myFnl_tot->Add(myFnl_tri,1.0); myFnl_tot->Add(myFnl_noncol,1.0); for(int i=0; iGetBinContent(i+1); bin_error=myFnl_tot->GetBinError(i+1); myFnl_tot_pG->SetBinContent(i+1,bin_value+bin_error); if(bin_value-bin_error>0.0) myFnl_tot_mG->SetBinContent(i+1,bin_value-bin_error); else myFnl_tot_mG->SetBinContent(i+1,0.0); myFnl_tot_mG->SetBinError(i+1,0.0); myFnl_tot_pG->SetBinError(i+1,0.0); } std::cout<<"Integral="<GetXaxis()->SetTitle("MetSig"); break; case 2: myFnl_tot->GetXaxis()->SetTitle("#slash{E}_{T}, GeV"); break; case 3: myFnl_tot->GetXaxis()->SetTitle("E_{T}^{#gamma1}, GeV"); break; case 4: myFnl_tot->GetXaxis()->SetTitle("E_{T}^{#gamma2}, GeV"); break; case 5: myFnl_tot->GetXaxis()->SetTitle("E_{T}^{em3}, GeV"); break; case 6: myFnl_tot->GetXaxis()->SetTitle("E_{T}^{jet1}, GeV"); break; case 7: myFnl_tot->GetXaxis()->SetTitle("H_{T} (GeV)"); break; case 8: myFnl_tot->GetXaxis()->SetTitle("M_{#gamma#gamma}, GeV/c^{2}"); break; case 9: myFnl_tot->GetXaxis()->SetTitle("M_{#gammaj}, GeV/c^{2}"); break; case 10: myFnl_tot->GetXaxis()->SetTitle("M_{emL-emE}, GeV/c^{2}"); break; case 11: myFnl_tot->GetXaxis()->SetTitle("M_{jj}, GeV/c^{2}"); break; case 12: myFnl_tot->GetXaxis()->SetTitle("N_{em}"); break; case 13: myFnl_tot->GetXaxis()->SetTitle("N_{jet}(E_{T}>15)"); break; case 14: myFnl_tot->GetXaxis()->SetTitle("Q_{T}, GeV"); break; case 15: myFnl_tot->GetXaxis()->SetTitle("#Delta#phi_{#gamma#gamma} (rad)"); break; case 16: myFnl_tot->GetXaxis()->SetTitle("#Delta#phi_{#slash{E}_{T}-#gamma1}, rad"); break; case 17: myFnl_tot->GetXaxis()->SetTitle("#Delta#phi_{#slash{E}_{T}-#gamma2}, rad"); break; case 18: myFnl_tot->GetXaxis()->SetTitle("#Delta#phi_{#slash{E}_{T}-jet1}, rad"); break; case 19: myFnl_tot->GetXaxis()->SetTitle("N_{jet}(E_{T}>20)"); break; case 20: myFnl_tot->GetXaxis()->SetTitle("N_{jet}(E_{T}>25)"); break; case 21: myFnl_tot->GetXaxis()->SetTitle("#Delta#phi_{#slash{E}_{T}-jet1}, rad"); break; case 22: myFnl_tot->GetXaxis()->SetTitle("#Delta#phi_{#slash{E}_{T}-jet1}, rad"); break; case 23: myFnl_tot->GetXaxis()->SetTitle("#Delta R, two EM objects"); break; default: myFnl_tot->GetXaxis()->SetTitle("#slash{E}_{T}, GeV"); } myFnl_tot->SetStats(0); myFnl_tot->GetYaxis()->SetTitle("Events per 0.5"); myFnl_tot->GetXaxis()->SetLabelSize(0.04); myFnl_tot->GetYaxis()->SetLabelSize(0.04); myFnl_tot->GetYaxis()->SetTitleSize(.055); myFnl_tot->GetXaxis()->SetTitleSize(.055); myFnl_tot->SetFillColor(5); myFnl_tot->SetLineWidth(1); if(plot_id<15) { myFnl_tot->SetMinimum(0.01); myFnl_tot->SetMaximum(1.5*h_int); } if(plot_id==19 || plot_id==20) { myFnl_tot->SetMinimum(0.01); myFnl_tot->SetMaximum(10.0*h_int); } if(plot_id>=15 && plot_id<19) { myFnl_tot->SetMinimum(0.0); myFnl_tot->SetMaximum(3.4); } if(plot_id>=21) { myFnl_tot->SetMinimum(0.0); myFnl_tot->SetMaximum(h_int); } //--- drawing Total Background myFnl_tot->SetFillStyle(1001); myFnl_tot->Draw("hist"); //___ drawing Total Background+Sigma myFnl_tot_pG->SetLineColor(2); // myFnl_tot_pG->Draw("same"); //___ drawing "Ele+pho" myFnl_MCgg->SetFillStyle(1001); myFnl_MCgg->SetFillColor(7); myFnl_MCgg->Draw("hist same"); myFnl_wvx->SetFillStyle(1001); myFnl_wvx->SetFillColor(40); // myFnl_wvx->Draw("hist same"); myFnl_tri->SetFillStyle(1001); myFnl_tri->SetFillColor(20); // myFnl_tri->Draw("hist same"); //___ drawing Non-collision if(myFnl_noncol->GetEntries()>0) { myFnl_noncol->SetFillStyle(1001); myFnl_noncol->SetFillColor(6); myFnl_noncol->Draw("hist same"); } //___ drawing Total Background-Sigma myFnl_tot_mG->SetLineColor(2); // myFnl_tot_mG->Draw("same"); //___ drawing GMSB signal myFnl_sig->SetLineColor(1); myFnl_sig->SetLineStyle(2); myFnl_sig->SetLineWidth(3); myFnl_sig->Add(myFnl_tot); // myFnl_sig->Draw("hist same"); // myFnl_sig->Add(myFnl_tot); //___ drawing Data myFnl_dt->SetMarkerStyle(20); myFnl_dt->SetMarkerSize(1); myFnl_dt->Draw("e p same"); // myFnl_dt->Draw("axissame"); TPaveText *pt = new TPaveText(0.518367,0.916168,0.918143,0.979042,"brNDC"); // TPaveText *pt = new TPaveText(0.568233,0.808383,0.868009,0.871257,"brNDC"); pt->SetName(""); pt->SetFillStyle(0); pt->SetLineColor(0); pt->SetTextFont(42); pt->SetTextSize(0.0389222); TText *text = pt->AddText("CDF Run II Preliminary, 2.6 fb^{-1}"); pt->Draw(); TPaveText *pt = new TPaveText(0.0367,0.916168,0.438143,0.979042,"brNDC"); // TPaveText *pt = new TPaveText(0.568233,0.808383,0.868009,0.871257,"brNDC"); pt->SetName(""); pt->SetFillStyle(0); pt->SetLineColor(0); pt->SetTextFont(42); pt->SetTextSize(0.0389222); TText *text = pt->AddText("#gamma#gamma+#slash{E}_{T} analysis in GMSB"); pt->Draw(); TLegend *leg1; if(plot_id<15) leg1= new TLegend(0.35,0.55,0.85,0.85,""); // if(plot_id==19 || plot_id==20) leg1= new TLegend(0.45,0.60,0.80,0.85,""); // if(plot_id>=15 && plot_id<19) leg1= new TLegend(0.15,0.60,0.50,0.85,""); // if(plot_id>=21) leg1= new TLegend(0.15,0.60,0.50,0.85,""); // leg1= new TLegend(0.45,0.60,0.80,0.85,""); // leg1->SetHeader("#tilde{#chi}^{0}_{1} mass=140 GeV, lifetime=0 ns"); leg1->SetFillColor(0); leg1->SetFillStyle(0); leg1->SetBorderSize(0); leg1->SetTextSize(0.063); leg1->SetTextFont(42); leg1->AddEntry(myFnl_dt,"data","LP"); // leg1->AddEntry(myFnl_tot_pG,"Bckg#pm#sigma_{total}","f"); leg1->AddEntry(myFnl_tot,"QCD with fake #slash{E}_{T}","f"); // leg1->AddEntry(myFnl_wvx,"Wrong Vertex with fake #slash{E}_{T}","f"); // leg1->AddEntry(myFnl_tri,"Tri-Photon with fake #slash{E}_{T}","f"); TF1 *func; if(plot_id==1) { func=new TF1("sig",SigFun,0.0,6.7,1); func->SetParameter(0,f_par); func->SetLineColor(2); // func->Draw("same"); // leg1->AddEntry(func,"Fake #slash{E}_{T} only","L"); } leg1->AddEntry(myFnl_MCgg,"EWK with real #slash{E}_{T}","f"); leg1->AddEntry(myFnl_noncol,"Non-collision","f"); // leg1->AddEntry(myFnl_sig,"GMSB signal","f"); // leg1->AddEntry("#tilde{#chi}^{0}_{1} mass=140 GeV, lifetime=0 ns"); leg1->Draw(); TLatex la; la.SetTextColor(1); la.SetNDC(); la.SetTextAlign(10); la.SetTextSize(0.02); la.DrawLatex(.45,.35,"Overflow Bin"); if(plot_id>1) { TLatex lax2; lax2.SetTextAlign(0); lax2.SetTextSize(0.04); double y_pos=0.35*h_int; double x_pos=0.15*highedge; if(plot_id<15) { y_pos=0.15*h_int; x_pos=0.60*highedge; } if(plot_id==19 || plot_id==20) { y_pos=0.15*h_int; x_pos=0.60*highedge; } // lax2.DrawLatex(x_pos,y_pos,"#slash{E}_{T}-significance>5.0"); } myFnl_dt->Draw("axissame"); // in ROOT prompt do > delete gDirectory->GetFile() }