ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/Development/Plotting/Modules/ZbTools.C
(Generate patch)

Comparing UserCode/cbrown/Development/Plotting/Modules/ZbTools.C (file contents):
Revision 1.9 by buchmann, Wed Nov 14 10:46:45 2012 UTC vs.
Revision 1.10 by buchmann, Wed Nov 14 14:52:54 2012 UTC

# Line 23 | Line 23 | using namespace PlottingSetup;
23   namespace ZbData {
24   vector<float> data_over_mc;
25  
26 < TCut ZplusBsel("pt1>20&&pt2>20&&mll>60&&mll<120");
26 > TCut ZplusBsel("pt1>20&&pt2>20&&mll>60&&mll<120&&id1==id2");
27   TCut LeadingB("Zb3010_bTagProbCSVBP[0]>0.898");
28   TCut EtaB("abs(Zb3010_pfJetGoodEta[0])<1.3");
29 < TCut PhiZcut("abs(pfJetDphiZ[0])>2.7");
29 > TCut PhiZcut("abs(Zb3010_pfJetDphiZ[0])>2.7");
30    
31   }
32  
# Line 220 | Line 220 | void draw_Zb_kin_vars() {
220     draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>125&&pt<1000",40,0,2,"#alpha","Official/alpha_pt_125_to_1000",1);
221      draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000",40,0,2,"#alpha","Official/alpha_pt_30_to_1000",1);
222    
223 <  draw_kin_variable("pfBJetDphiZ[0]",ZplusBsel&&LeadingB&&"pt>10&&pt<1000&&pfBJetDphiZ[0]>0",30,0,3.2,"#delta#phi (Z,b lead)","DeltaPhiZBlead",0);
223 >  draw_kin_variable("Zb3010_pfBJetDphiZ[0]",ZplusBsel&&LeadingB&&"pt>10&&pt<1000&&pfBJetDphiZ[0]>0",30,0,3.2,"#delta#phi (Z,b lead)","DeltaPhiZBlead",0);
224    
225   /*
226    draw_kin_variable("pt1",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,0,200,"p_{T}^{l1}","Lep/pt1",1);
# Line 259 | Line 259 | Value get_Zb_data_over_mc(string variabl
259    TH1F *hdata  = allsamples.Draw("hdata",variable,1000,0,30, "#alpha [GeV]", "events", cut,data,luminosity);
260    TH1F *hmc    = allsamples.Draw("hmc" , variable,1000,0,30, "#alpha [GeV]", "events", cut,mc,  luminosity);
261    
262 <  hdata->Fit("gaus","QE");
263 <  hmc->Fit("gaus","QE");
264 <  TF1 *fdata = (TF1*) hdata->GetFunction("gaus");
265 <  TF1 *fmc   = (TF1*) hmc->GetFunction("gaus");
266 <    
267 < /*  float a=hdata->GetMean();
262 >  float a=hdata->GetMean();
263    float b=hmc->GetMean();
264    float da=hdata->GetMeanError();
265 <  float db=hmc->GetMeanError();*/
271 <  float a = fdata->GetParameter(1);
272 <  float da = fdata->GetParError(1);
273 <  float b = fmc->GetParameter(1);
274 <  float db = fmc->GetParError(1);
265 >  float db=hmc->GetMeanError();
266    float factor = a / b;
267    float error = TMath::Sqrt(  (1/(b*b)) * (da*da) + ((a*a)/(b*b))*db*db);
268    delete hdata;
# Line 894 | Line 885 | void do_basic_ZB_analysis() {
885   //     compare_selection("Zb40");
886  
887    
888 < //  print_all_b_yields();
889 < //  draw_mpf_vars();
888 >  print_all_b_yields();
889 >  draw_mpf_vars();
890   //  DrawEvilCutFlow();
891    
892 < //  draw_Zb_kin_vars();
892 >  draw_Zb_kin_vars();
893    
894    new_data_mc_agreement_2d();
895    

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines