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

Comparing UserCode/cbrown/Development/Plotting/Modules/Plotting_Functions.C (file contents):
Revision 1.52 by buchmann, Thu Aug 16 14:30:32 2012 UTC vs.
Revision 1.75 by buchmann, Tue Dec 4 12:31:52 2012 UTC

# Line 33 | Line 33 | void todo() {
33  
34  
35  
36 < void find_one_peak_combination(TCut specialcut, float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, float &MCSigma, float &MCSigmaError, float &DataSigma, float& DataSigmaError, stringstream &result, bool doPUreweighting = true, string saveas="")
36 > void find_one_peak_combination(TCut specialcut, bool SwitchOffNJetsCut, float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, float &MCSigma, float &MCSigmaError, float &DataSigma, float& DataSigmaError, stringstream &result, bool doPUreweighting = true, string saveas="")
37   {
38    // Temporarily switch off PU reweighting, if asked
39    TCut weightbackup=cutWeight;
# Line 42 | Line 42 | void find_one_peak_combination(TCut spec
42    int nbins=100;
43    if(PlottingSetup::DoBTag) nbins=25;
44    
45 +  TCut nJetsCut(cutnJets);
46 +  if(SwitchOffNJetsCut) nJetsCut=specialcut;
47 +  
48    TCanvas *tempcan = new TCanvas("tempcan","Temporary canvas for peak finding preparations");
49 <  TH1F *rawJZBeemmMC      = allsamples.Draw("rawJZBeemmMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&specialcut,mc, luminosity);
50 <  TH1F *rawJZBeemmData    = allsamples.Draw("rawJZBeemmData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&specialcut,data, luminosity);
51 <  TH1F *rawJZBemMC      = allsamples.Draw("rawJZBemMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&specialcut,mc, luminosity);
52 <  TH1F *rawJZBemData    = allsamples.Draw("rawJZBemData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&specialcut,data, luminosity);
49 >  TH1F *rawJZBeemmMC      = allsamples.Draw("rawJZBeemmMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&nJetsCut&&specialcut,mc, luminosity);
50 >  TH1F *rawJZBeemmData    = allsamples.Draw("rawJZBeemmData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&nJetsCut&&specialcut,data, luminosity);
51 >  TH1F *rawJZBemMC      = allsamples.Draw("rawJZBemMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&nJetsCut&&specialcut,mc, luminosity);
52 >  TH1F *rawJZBemData    = allsamples.Draw("rawJZBemData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&nJetsCut&&specialcut,data, luminosity);
53    TH1F *rawttbarjzbeemmMC;
54  
55    if(method==doKM) {
56      //we only need this histo for the KM fitting...
57 <    rawttbarjzbeemmMC = allsamples.Draw("rawttbarjzbeemmMC",jzbvariablemc,nbins, -50,50, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&specialcut,mc,luminosity,allsamples.FindSample("TTJet"));
57 >    rawttbarjzbeemmMC = allsamples.Draw("rawttbarjzbeemmMC",jzbvariablemc,nbins, -50,50, "JZB [GeV]", "events",cutmass&&cutOSSF&&nJetsCut&&specialcut,mc,luminosity,allsamples.FindSample("TTJet"));
58      MCPeak=find_peak(rawJZBeemmMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,MCSigmaError,method,saveas);
59      DataPeak=find_peak(rawJZBeemmData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,DataSigmaError,method,saveas);
60      delete rawttbarjzbeemmMC;
# Line 84 | Line 87 | void find_one_peak_combination(TCut spec
87    delete tempcan;
88   }
89  
90 < void find_peaks(float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, stringstream &result, bool doPUreweighting, stringstream &datajzb, stringstream &mcjzb)
90 > void find_peaks(float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, stringstream &result, bool doPUreweighting, stringstream &datajzb, stringstream &mcjzb, string sSpecialCut="", bool SwitchOffNJetsCut=false)
91   {
92 +  bool overunderflowstatus=addoverunderflowbins;
93 +  switch_overunderflow(false);
94 +  
95 +  TCut SpecialCut("mll>=0");
96 +  if(sSpecialCut!="") SpecialCut=TCut(sSpecialCut.c_str());
97    
98    bool DoInvidualeemmPeaks=false;
99    
# Line 101 | Line 109 | void find_peaks(float &MCPeak,float &MCP
109    float mcSigma,mcSigmaError, dataSigma, dataSigmaError;
110    
111    dout << "Finding global peak : " << endl;
112 <  find_one_peak_combination(TCut(""),mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,"");
112 >  find_one_peak_combination(SpecialCut,SwitchOffNJetsCut,mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,"");
113    
114    if(DoInvidualeemmPeaks) {
115      dout << "Finding peak for electrons : " << endl;
116 <    find_one_peak_combination(TCut("id1==0"),mceepeak,mceepeakerr,dataeepeak,dataeepeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_ele");
116 >    find_one_peak_combination(SpecialCut&&TCut("id1==0"),SwitchOffNJetsCut,mceepeak,mceepeakerr,dataeepeak,dataeepeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_ele");
117      dout << "Finding peak for muons : " << endl;
118 <    find_one_peak_combination(TCut("id1==1"),mcmmpeak,mcmmpeakerr,datammpeak,datammpeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_mu");
118 >    find_one_peak_combination(SpecialCut&&TCut("id1==1"),SwitchOffNJetsCut,mcmmpeak,mcmmpeakerr,datammpeak,datammpeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_mu");
119      
120      datajzb << "(" << jzbvariabledata;
121      mcjzb << "(" << jzbvariablemc;
# Line 133 | Line 141 | void find_peaks(float &MCPeak,float &MCP
141      else mcjzb << "+ (id1!=id2)*" << TMath::Abs(mcpeak) << " ";
142      
143      mcjzb << ")";
144 +    
145    } else {
146      datajzb << "(" << jzbvariabledata;
147        mcjzb << "(" << jzbvariablemc;
# Line 148 | Line 157 | void find_peaks(float &MCPeak,float &MCP
157      mcjzb << ")";
158    }
159    
160 <  
160 >    MCPeak=mcpeak;
161 >    MCPeakError=mcpeakerr;
162 >    DataPeak=datapeak;
163 >    DataPeakError=datapeakerr;
164 >    switch_overunderflow(overunderflowstatus);
165   }
166  
167   void make_special_obs_pred_mll_plot(string datajzb, string mcjzb, float jzbthreshold, float binWidth = 5.0) {
# Line 292 | Line 305 | void make_special_mll_plot(int nbins, fl
305    datahistoOSSF->SetMarkerSize(DataMarkerSize);
306    datahistoOSSF->Draw();
307  
308 <  mcstackOSSF.Draw("same");
308 >  mcstackOSSF.Draw("histo,same");
309    datahistoOSSF->Draw("same");
310  
311    datahistoOSOF->SetMarkerColor(TColor::GetColor("#FE642E"));
# Line 357 | Line 370 | void draw_ratio_plot(TH1* hdata, THStack
370    oneline->SetLineWidth(1);
371    oneline->Draw("same");
372   }
373 +
374   float make_one_OFSF_plot(string variable, string addcut, string legendTitle, int nbins, float min, float max, float ymax, bool logscale,
375 <                        string xlabel, string filename, bool plotratio=true, bool loadlastminmax=false, float legendPosition=0.55) {
375 >                        string xlabel, string filename, float legendPosition=0.55) {
376  
377    TCut ibasiccut=basiccut;
378    bool draw_separation_lines=false;
# Line 372 | Line 386 | float make_one_OFSF_plot(string variable
386    cutSF = cutOSSF&&cutnJets&&ibasiccut;
387  
388    TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
389 <  ckin->SetLogy(logscale);
389 >  TPad* rcan = new TPad("rcan","rcan",0,0,1,1);
390 >  rcan->SetLogy(logscale);
391 >  rcan->cd();
392 >
393 >  std::cout << "OF/SF comparison: variable = " << variable << ", cut = " << cutSF.GetTitle() << std::endl;
394    TH1F *datahistoSF = allsamples.Draw("datahistoSF",variable,nbins,min,max, xlabel, "events",cutSF,data,luminosity);
395    TH1F *datahistoOF = allsamples.Draw("datahistoOF",variable,nbins,min,max, xlabel, "events",cutOF,data,luminosity);
396 <  string signal("LM3");
397 <  TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max);
398 <  int idx = signalsamples.FindSample(signal)[0];
399 <  (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cutSF);
400 <  signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
401 <  signalhisto->SetLineColor((signalsamples.collection)[idx].samplecolor);
402 <  signalhisto->SetLineStyle(2);
396 > //   string signal("LM3");
397 > //   TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max);
398 > //   int idx = signalsamples.FindSample(signal)[0];
399 > //   (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cutSF);
400 > //   signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
401 > //   signalhisto->SetLineColor((signalsamples.collection)[idx].samplecolor);
402 > //   signalhisto->SetLineStyle(2);
403    datahistoSF->SetMarkerSize(DataMarkerSize);
404    datahistoOF->SetLineColor(kRed);
405  
# Line 392 | Line 410 | float make_one_OFSF_plot(string variable
410    }
411    if (ymax<0) {
412      if ( logscale ) datahistoSF->SetMaximum(5.3*datahistoSF->GetMaximum());
413 <    else datahistoSF->SetMaximum(1.5*datahistoSF->GetMaximum());
413 >    else datahistoSF->SetMaximum(1.8*datahistoSF->GetMaximum());
414    } else {
415      datahistoSF->SetMaximum(ymax);
416    }
# Line 411 | Line 429 | float make_one_OFSF_plot(string variable
429    } else {
430      mleg->AddEntry((TObject*)0, "", "");
431    }
432 <  mleg->AddEntry(signalhisto, "LM3", "L");
432 >  //mleg->AddEntry(signalhisto, "LM3", "L");
433  
434    datahistoSF->Draw("E1");
435    if (datahistoOF->Integral()>0) datahistoOF->Draw("HIST,SAMES");
436 <  signalhisto->Draw("HIST,SAMES");
436 >  //signalhisto->Draw("HIST,SAMES");
437    mleg->Draw();
438    DrawPrelim();
439 <  CompleteSave(ckin, "SFOF/" + filename);
439 >  if (datahistoOF->Integral()>0) {
440 >    save_with_ratio( datahistoSF, datahistoOF, rcan, "SFOF/" + filename, false, false, "SF/OF" );
441 >  } else {
442 >    CompleteSave(rcan, "SFOF/" + filename);
443 >    delete rcan;
444 >  }
445  
446    datahistoSF->Delete();
447    datahistoOF->Delete();
448 <  signalhisto->Delete();
448 >  //signalhisto->Delete();
449    delete mleg;
450    delete ckin;
451  
# Line 430 | Line 453 | float make_one_OFSF_plot(string variable
453  
454   }
455  
456 < void make_OFSF_plots(string variable, string addcut, int nbins, float min, float max, bool logscale, string xlabel, string filename,
457 <                bool plotratio=true, bool loadlastminmax=false, float legendPosition=0.55) {
456 > // Compare data to data
457 > float make_data_comparison_plot(string variable, TCut cut, int nbins, float min, float max, float ymax, bool logscale,
458 >                                string xlabel, string filename, float legendPosition=0.55) {
459 >
460 >  TCut ibasiccut=basiccut&&cut;
461 >
462 >  TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
463 >  TPad* rcan = new TPad("rcan","rcan",0,0,1,1);
464 >  rcan->SetLogy(logscale);
465 >  rcan->cd();
466 >
467 >  std::cout << "Data comparison: variable = " << variable << ", cut = " << ibasiccut.GetTitle() << std::endl;
468 >
469 >  TH1F *data1 = allsamples.Draw("data1",variable,nbins,min,max, xlabel, "events",ibasiccut,data,luminosity);
470 >  TH1F *data2 = comparesamples.Draw("data2",variable,nbins,min,max, xlabel, "events",ibasiccut,data,luminosity);
471 >  data1->Scale(5.0/9.2);
472 >
473 >  data1->SetMarkerSize(DataMarkerSize);
474 >  data2->SetLineColor(kRed);
475  
476 <  string mllcuts[] = { "mll>20","mll>20&&mll<70", "mll>70&&mll<110", "mll>110" };
477 <  string mllcutname[] = { "m_{ll} > 20 GeV", "20 < m_{ll} < 70 GeV", "70 < m_{ll} < 110 GeV", "m_{ll} > 110 GeV" };
476 >  if ( !logscale ) {
477 >    data1->SetMinimum(0); // Defaults
478 >  } else {
479 >    data1->SetMinimum(0.5);
480 >  }
481 >  if (ymax<0) {
482 >    if ( logscale ) data1->SetMaximum(5.3*data1->GetMaximum());
483 >    else data1->SetMaximum(1.5*data1->GetMaximum());
484 >  } else {
485 >    data1->SetMaximum(ymax);
486 >  }
487 >
488 >  float ymaxSet = data1->GetMaximum();
489 >
490 >  data1->GetXaxis()->SetTitle(xlabel.c_str());
491 >  data1->GetYaxis()->SetTitle("Events");
492 >  data1->GetXaxis()->CenterTitle();
493 >  data1->GetYaxis()->CenterTitle();
494 >
495 >  TLegend *mleg = make_legend("",legendPosition,0.7,false,legendPosition+0.2);
496 >  mleg->AddEntry(data1, "2012 data (scaled)", "PL");
497 >  mleg->AddEntry(data2, "2011 data", "L");
498 >
499 >  data1->Draw("E1");
500 >  data2->Draw("HIST,SAMES");
501 >  mleg->Draw();
502 >  DrawPrelim();
503 >  save_with_ratio( data1, data2, rcan, "compareData/" + filename, false, false, "new/old" );
504 >
505 >  data1->Delete();
506 >  data2->Delete();
507 >  delete mleg;
508 >  delete ckin;
509 >
510 > }
511 >
512 > void make_OFSF_plots(string variable, string addcut, int nbins, float min, float max, bool logscale,
513 >                     string xlabel, string filename, float legendPosition=0.55) {
514 >
515 >  string mllcuts[] = { "mll>20","mll>15&&mll<70", "mll>75&&mll<105", "mll>120" };
516 >  string mllcutname[] = { "m_{ll} > 20 GeV", "20 < m_{ll} < 70 GeV", "70 < m_{ll} < 110 GeV", "m_{ll} > 120 GeV" };
517    string plotname[] = {"_all","_low","_peak","_high"};
518    float ymax;
519 <  for ( int i=0; i<4; ++i ) {
519 >
520 >  int start = 0;
521 >  if ( !PlottingSetup::openBox ) start = 3;
522 >
523 >  for ( int i=start; i<4; ++i ) {
524      if ( addcut != "" ) mllcuts[i] += "&&"+addcut;
525 <    if ( i==0 ) {
526 <      ymax = make_one_OFSF_plot(variable, mllcuts[i],            mllcutname[i], nbins, min, max, -1, logscale, xlabel,
527 <                         filename+plotname[i], plotratio, loadlastminmax, legendPosition );
525 >    if ( i==start ) {
526 >      ymax = make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, -1, logscale, xlabel,
527 >                         filename+plotname[i], legendPosition );
528      } else {
529 <      make_one_OFSF_plot(variable, mllcuts[i],            mllcutname[i], nbins, min, max, ymax, logscale, xlabel,
530 <                         filename+plotname[i], plotratio, loadlastminmax, legendPosition );
529 >      make_one_OFSF_plot(variable, mllcuts[i],        mllcutname[i], nbins, min, max, ymax, logscale, xlabel,
530 >                         filename+plotname[i], legendPosition );
531      }
532      make_one_OFSF_plot(variable, "id1==1&&id1==id2&&"+mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel,
533 <                       filename+plotname[i]+"_mm", plotratio, loadlastminmax, legendPosition );
533 >                       filename+plotname[i]+"_mm", legendPosition );
534      make_one_OFSF_plot(variable, "id1==0&&id1==id2&&"+mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel,
535 <                       filename+plotname[i]+"_ee", plotratio, loadlastminmax, legendPosition );
535 >                       filename+plotname[i]+"_ee", legendPosition );
536    }
537  
538   }
# Line 463 | Line 546 | void make_kin_plot(string variable, stri
546   //  TCut basiccut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])&&(passed_triggers||!is_data)");
547    TCut ibasiccut=basiccut;
548    bool draw_separation_lines=false;
549 +  bool drawsignal = false;
550    
551    if(isPF) ibasiccut=basiccut&&"pfjzb[0]>-998";
552  
# Line 484 | Line 568 | void make_kin_plot(string variable, stri
568    //Step 2: Refine the cut
569    TCut cut;
570    cut=cutmass&&cutOSSF&&cutnJets&&ibasiccut;
571 <  if(filename=="nJets") cut=cutmass&&cutOSSF&&ibasiccut;
572 <  if(filename=="nJets_osof") cut=cutmass&&cutOSOF&&ibasiccut;
571 >  if(filename=="nJets" || filename=="nJets_inclusive" || filename=="nJets_met100" || filename=="nJets_inclusive_met100") cut=cutmass&&cutOSSF&&ibasiccut;
572 >  if(filename=="nJets_osof" || filename=="nJets_osof_inclusive" || filename=="nJets_osof_met100" || filename=="nJets_osof_inclusive_met100") cut=cutmass&&cutOSOF&&ibasiccut;
573    if(filename=="nJets_nocuts_except_mll_ossf") cut=cutmass&&cutOSSF;
574 <  if(filename=="mll") {
574 >  if(filename=="mll"||filename=="mll_met100") {
575          cut=cutOSSF&&cutnJets&&ibasiccut;
576          draw_separation_lines=true;
577    }
578 <  if(filename=="mll_ee") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==0";
579 <  if(filename=="mll_osof") {
578 >  if(filename=="mll_ee"||filename=="mll_ee_met100") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==0";
579 >  if(filename=="mll_mm"||filename=="mll_mm_met100") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==1";
580 >  if(filename=="mll_osof"||filename=="mll_osof_met100") {
581          cut=cutOSOF&&cutnJets&&ibasiccut;
582          draw_separation_lines=true;
583    }
499 –  if(filename=="mll_mm") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==1";
584    if(Contains(filename,"aboveJZB")) cut=cutOSSF&&cutnJets&&ibasiccut;
585    if(Contains(filename,"mll_ee_above")) cut=cut&&"id1==0";
586    if(Contains(filename,"mll_mm_above")) cut=cut&&"id1==1";
# Line 507 | Line 591 | void make_kin_plot(string variable, stri
591    if(filename=="mll_inclusive_mm") cut=cutOSSF&&"id1==1";
592    if(filename=="pfJetGoodEta_0") cut=cutOSSF&&cutmass&&ibasiccut&&cutnJets;
593    if(filename=="pfJetGoodPt_0") cut=cutOSSF&&cutmass&&ibasiccut&&cutnJets;
594 +  if(filename=="numVtx") cut=cutmass&&ibasiccut;
595  
596    TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
597    ckin->SetLogy(logscale);
# Line 524 | Line 609 | void make_kin_plot(string variable, stri
609    if (logscale) datahisto->SetMaximum(5.3*datahisto->GetMaximum());
610    else datahisto->SetMaximum(1.3*datahisto->GetMaximum());
611  
612 <  cout << "******** Cut used : " << (const char*) cut << endl;
612 >  cout << "******** Cut used : " << (const char*) cut << " for plot " << filename << endl;
613    if(loadlastminmax) {
614          datahisto->SetMinimum(lastrange_min);
615          datahisto->SetMaximum(lastrange_max);
# Line 538 | Line 623 | void make_kin_plot(string variable, stri
623    string signal("LM3");
624    TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max);
625    int idx = signalsamples.FindSample(signal)[0];
626 <  (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cut);
627 <  signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
628 <  signalhisto->SetLineColor(kOrange);
626 >  if(drawsignal) (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cut);
627 >  if(drawsignal) signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
628 >  if(drawsignal) signalhisto->SetLineColor(kOrange);
629  
630    THStack mcstack   = allsamples.DrawStack("mcstack",  variable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
631    datahisto->Draw("e1");
632    ckin->Update();
633 <  mcstack.Draw("same");
633 >  mcstack.Draw("histo,same");
634    
635    datahisto->Draw("same,e1");
636    TLegend *kinleg = allsamples.allbglegend();
# Line 588 | Line 673 | void make_kin_plot(string variable, stri
673          kinpad->cd();
674          kinpad->SetLogy(logscale);
675          datahisto->Draw("e1");
676 <        mcstack.Draw("same");
677 <        signalhisto->Draw("same");
676 >        mcstack.Draw("histo,same");
677 >        if(drawsignal) signalhisto->Draw("same");
678          datahisto->Draw("same,e1");
679          datahisto->Draw("same,axis");
680          if(RestrictToMassPeak&&draw_separation_lines) {
# Line 597 | Line 682 | void make_kin_plot(string variable, stri
682            upperboundary->Draw("same");
683          }
684            
685 <        kinleg->AddEntry("signalhisto",signal.c_str(),"l");
685 >        if(drawsignal) kinleg->AddEntry("signalhisto",signal.c_str(),"l");
686          kinleg->Draw();
687          write_cut->Draw();
688          DrawPrelim();
# Line 612 | Line 697 | void make_kin_plot(string variable, stri
697          else CompleteSave(ckin,"kin/"+filename);
698    }
699    datahisto->Delete();
700 +  delete signalhisto;
701 +  delete ckin;
702 + }
703 +
704 +
705 + void make_plain_kin_plot(string variable, string addcut, int nbins, float min, float max, bool logscale,
706 +                   string xlabel, string filename, bool isPF=true, bool plotratio=true, bool loadlastminmax=false ) {
707 +  TCut cut=TCut(addcut.c_str())&&basiccut;
708 +
709 +  TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
710 +  ckin->SetLogy(logscale);
711 +  TH1F *datahisto = allsamples.Draw("datahisto",variable,nbins,min,max, xlabel, "events",cut,data,luminosity);
712 +  datahisto->SetMarkerSize(DataMarkerSize);
713 +  if ( !logscale ) datahisto->SetMinimum(0); // Defaults
714 +  else datahisto->SetMinimum(0.5);
715 +  // Custom max.
716 +  if (logscale) datahisto->SetMaximum(5.3*datahisto->GetMaximum());
717 +  else datahisto->SetMaximum(1.3*datahisto->GetMaximum());
718 +
719 +  cout << "******** Cut used : " << (const char*) cut << " for plot " << filename << endl;
720 +  if(loadlastminmax) {
721 +        datahisto->SetMinimum(lastrange_min);
722 +        datahisto->SetMaximum(lastrange_max);
723 +        if(logscale) {
724 +                datahisto->SetMinimum(pow(10,lastrange_min));
725 +                datahisto->SetMaximum(pow(10,lastrange_max));
726 +        }              
727 +  }
728 +
729 +  THStack mcstack   = allsamples.DrawStack("mcstack",  variable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
730 +  datahisto->Draw("e1");
731 +  ckin->Update();
732 +  mcstack.Draw("histo,same");
733 +  
734 +  datahisto->Draw("same,e1");
735 +  TLegend *kinleg = allsamples.allbglegend();
736 +  kinleg->Draw();
737 +
738 +  lastrange_min=ckin->GetUymin();
739 +  lastrange_max=ckin->GetUymax();
740 +
741 +  if ( plotratio ) {
742 +        TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
743 +        kinpad->cd();
744 +        kinpad->SetLogy(logscale);
745 +        datahisto->Draw("e1");
746 +        mcstack.Draw("histo,same");
747 +        datahisto->Draw("same,e1");
748 +        datahisto->Draw("same,axis");
749 +          
750 +        kinleg->Draw();
751 +        DrawPrelim();
752 +        string saveas="kin/"+filename;
753 +        if(isPF) saveas="kin/"+filename+"__PF";
754 +        save_with_ratio(datahisto,mcstack,kinpad->cd(),saveas);
755 +  } else {
756 +        if(isPF) CompleteSave(ckin,"kin/"+filename+"__PF");
757 +        else CompleteSave(ckin,"kin/"+filename);
758 +  }
759 +  datahisto->Delete();
760    delete ckin;
761   }
762  
763 +
764   void make_JES_plot(TCut cut, string name) {
765  
766    int nbins=10;
# Line 648 | Line 794 | void make_JES_plot(TCut cut, string name
794    JESunc->SetFillColor(TColor::GetColor("#00ADE1"));
795    JESunc->SetFillStyle(3002);
796    datahisto->Draw("e1");
797 <  mcstack.Draw("same");
797 >  mcstack.Draw("histo,same");
798    JESunc->Draw("2");
799    datahisto->Draw("same,e1");
800    TLegend *kinleg = allsamples.allbglegend();
# Line 660 | Line 806 | void make_JES_plot(TCut cut, string name
806  
807   }
808  
809 + string Cut2Str(TCut cut) {
810 +  return ((string)(const char*)cut);
811 + }
812 +
813 + void MakeElegantTwoThreeComparisons() {
814 +    
815 +    TCut signal("met[4]>100&&met[4]<150&&mll>20"&&PlottingSetup::basicqualitycut);
816 +    TCut TwoJets("pfJetGoodNum40==2");
817 +    TCut ThreeJets("pfJetGoodNum40>=3");
818 +    
819 +    float min=20;
820 +    float max=300;
821 +    int nbins=int((max-min)/5.0);
822 +    string xlabel="m_{ll} [GeV]";
823 +    
824 +    TCanvas *can2 = new TCanvas("can2","can2");
825 +    TH1F *TwoOF = allsamples.Draw("TwoOF","mll",nbins,min,max, xlabel, "events",cutOSOF&&signal&&TwoJets,mc,luminosity,allsamples.FindSample("TTJets_"));
826 +    TH1F *TwoSF = allsamples.Draw("TwoSF","mll",nbins,min,max, xlabel, "events",cutOSSF&&signal&&TwoJets,mc,luminosity,allsamples.FindSample("TTJets_"));
827 +    TH1F *ThreeOF = allsamples.Draw("ThreeOF","mll",nbins,min,max, xlabel, "events",cutOSOF&&signal&&ThreeJets,mc,luminosity,allsamples.FindSample("TTJets_"));
828 +    TH1F *ThreeSF = allsamples.Draw("ThreeSF","mll",nbins,min,max, xlabel, "events",cutOSSF&&signal&&ThreeJets,mc,luminosity,allsamples.FindSample("TTJets_"));
829 +    
830 +    can2->cd(1);
831 +    
832 +    TwoSF->SetFillColor(kWhite);
833 +    TwoOF->SetFillColor(kWhite);
834 +    ThreeSF->SetFillColor(kWhite);
835 +    ThreeOF->SetFillColor(kWhite);
836 +    
837 +    TPad *kin3 = new TPad("kin3","kin3",0,0,1,1);
838 +    kin3->cd();
839 +    TwoOF->Draw("histo");
840 +    TwoSF->Draw("e1,same");
841 +    TLegend *leg_SFOF_Two = make_legend();
842 +    leg_SFOF_Two->SetHeader("t#bar{t} MC:");
843 +    leg_SFOF_Two->AddEntry(TwoOF,"OF, nJets==2","l");
844 +    leg_SFOF_Two->AddEntry(TwoSF,"SF, nJets==2","lp");
845 +    leg_SFOF_Two->Draw();
846 +    save_with_ratio(TwoSF,TwoOF,kin3->cd(),"DoubleRatios/SFOF_Two");
847 +    
848 +    TPad *kin4 = new TPad("kin4","kin4",0,0,1,1);
849 +    kin4->cd();
850 +    ThreeOF->Draw("histo");
851 +    ThreeSF->Draw("e1,same");
852 +    TLegend *leg_SFOF_Three = make_legend();
853 +    leg_SFOF_Three->SetHeader("t#bar{t} MC:");
854 +    leg_SFOF_Three->AddEntry(TwoOF,"OF, nJet#geq3","l");
855 +    leg_SFOF_Three->AddEntry(TwoSF,"SF, nJets#geq3","lp");
856 +    leg_SFOF_Three->Draw();
857 +    save_with_ratio(ThreeSF,ThreeOF,kin4->cd(),"DoubleRatios/SFOF_Three");
858 +    
859 +    
860 +    TwoSF->Scale(ThreeSF->Integral()/TwoSF->Integral());
861 +    TPad *kin = new TPad("kin","kin",0,0,1,1);
862 +    kin->cd();
863 +    TwoSF->Draw("histo");
864 +    ThreeSF->Draw("e1,same");
865 +    TLegend *leg_SF_Two_Three = make_legend();
866 +    leg_SF_Two_Three->SetHeader("t#bar{t} MC:");
867 +    leg_SF_Two_Three->AddEntry(TwoSF,"#splitline{SF, nJets==2}{ (scaled)}","l");
868 +    leg_SF_Two_Three->AddEntry(ThreeSF,"SF, nJets#geq3","lp");
869 +    leg_SF_Two_Three->Draw();
870 +    save_with_ratio(TwoSF,ThreeSF,kin->cd(),"DoubleRatios/SF_Two_Three");
871 +    
872 +    TwoOF->Scale(ThreeOF->Integral()/TwoOF->Integral());
873 +    TPad *kin2 = new TPad("kin2","kin2",0,0,1,1);
874 +    kin2->cd();
875 +    TwoOF->Draw("histo");
876 +    ThreeOF->Draw("e1,same");
877 +    TLegend *leg_OF_Two_Three = make_legend();
878 +    leg_OF_Two_Three->SetHeader("t#bar{t} MC:");
879 +    leg_OF_Two_Three->AddEntry(TwoSF,"#splitline{OF, nJets==2}{ (scaled)}","l");
880 +    leg_OF_Two_Three->AddEntry(ThreeSF,"OF, nJets#geq3","lp");
881 +    leg_OF_Two_Three->Draw();
882 +    save_with_ratio(TwoOF,ThreeOF,kin2->cd(),"DoubleRatios/OF_Two_Three");
883 +    
884 + }
885 +
886 +
887   void do_kinematic_plots(string mcjzb, string datajzb, bool doPF=false)
888   {
889 + //  switch_overunderflow(true);
890    bool dolog=true;
891    bool nolog=false;
892 +
893 +  bool doOFSF = false;
894 +  bool doKin  = true;
895 +  bool doDataComp = false;
896 +  bool MakeTwoThreeJetComparison = false;
897 +  
898 +
899    if(doPF) write_warning(__FUNCTION__,"Please use caution when trying to produce PF plots; not all versions of the JZB trees have these variables!");
900    float mll_low=50;
901    float mll_hi=160;
902    if(!PlottingSetup::RestrictToMassPeak) {
903          mll_low=20;
904 <        mll_hi=300;
904 >        mll_hi=320;
905    }
906 +  
907 +  if(MakeTwoThreeJetComparison) MakeElegantTwoThreeComparisons();
908 +
909 +  if ( doOFSF ) {
910 +     make_OFSF_plots("mll", "met[4]>100", 60, 20., 320., false, "m_{ll}", "mll");
911 +
912 + //      make_OFSF_plots("pfJetGoodNum40", "met[4]>100",                        7, 3, 10, true, "#(jets)", "njets");
913 + //      make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", "njets_btagVeto");
914 + //      make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30>0",  7, 3, 10, true, "#(jets)", "njets_AtLeastOneBJet30");
915 +
916 + //      make_OFSF_plots("pfJetGoodNumBtag30", "met[4]>100", 5, 0, 5, true, "#(b-jets)", "nbjets");
917 + //      make_OFSF_plots("pfJetGoodPtBtag[0]", "met[4]>100&&pfJetGoodNumBtag30>0", 20, 0, 400, true, "p_{T}(leading b-jet)", "ptb1");
918 +
919 + //     make_OFSF_plots("iso1", "met[4]>100", 20, 0, 0.3, true, "lepton 1 isolation", "iso1");
920 + //     make_OFSF_plots("iso2", "met[4]>100", 20, 0, 0.3, true, "lepton 2 isolation", "iso2");
921 +    //  make_OFSF_plots("pt1",  "met[4]>100", 30, 0., 300., true, "p_{T,1}", "pt1");
922 +    //  make_OFSF_plots("pt2",  "met[4]>100", 22, 0., 220., true, "p_{T,2}", "pt2");
923 + //     make_OFSF_plots("eta1", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{1}", "eta1", 0.15);
924 + //     make_OFSF_plots("eta2", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{2}", "eta2", 0.15);
925 +    //  make_OFSF_plots("phi1", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", "phi1", 0.2);
926 +    //  make_OFSF_plots("phi2", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", "phi2", 0.2);
927 + //     make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", "met[4]>100", 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", "jpt1pt2", 0.2);
928 +    make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", "dphij1met", 0.2);
929 + //    make_OFSF_plots("TMath::Abs(dphi)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", "dphi", 0.2);
930 + //    make_OFSF_plots("TMath::Abs(dphiMet1)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", "dphiMet1", 0.2);
931 + //    make_OFSF_plots("TMath::Abs(dphiMet2)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", "dphiMet2", 0.2);
932 + //    make_OFSF_plots("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(l,MET)|)", "dphilc");
933 + //    make_OFSF_plots("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(jet,MET)|)", "dphijc");
934 + //    make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(l,MET)|)", "dphilco");
935 + //    make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(pfJetDphiMet[0])), TMath::Min( (TMath::Pi()-TMath::Abs(pfJetDphiMet[1])), (TMath::Pi() - TMath::Abs(pfJetDphiMet[2]))))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(jet,MET)|)", "dphijco");
936 +  }
937 +
938 +  if ( doDataComp) {
939 +    TCut mllCut("");
940 +    float massmin = 15.;
941 +    float massmax = 315;
942 +    int massnbins = 60;
943 + /*    if ( !PlottingSetup::openBox ) {
944 +      mllCut = "mll>70";
945 +      massmin = 120;
946 +      massmax = 360;
947 +      massnbins = 14;
948 +    }*/
949 +
950 +    TCut cutSignal = cutmass&&cutnJets&&"met[4]>100";
951 +    
952 + /*    make_data_comparison_plot("numVtx", "",40, -0.5, 39.5,-1., true, "N(Vertices)", "numVtx");
953 +    make_data_comparison_plot("pfJetGoodEta[0]", "",50, -3.0, 3.0,-1., true, "N(Jets)", "pfJetGoodEta0");
954 +    
955 +    
956 +    make_data_comparison_plot("pfJetGoodNum40", cutOSOF,10, -0.5, 9.5,-1., true, "N(Jets)", "njets_OF");
957 +    make_data_comparison_plot("met[4]", cutOSSF,60, 0., 300.,-1., true, "pfMET", "met_SF_inclusive");
958 +    make_data_comparison_plot("met[1]", cutOSSF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_SF_inclusive");
959 +    make_data_comparison_plot("met[4]", cutOSOF,60, 0., 300.,-1., true, "pfMET", "met_OF_inclusive");
960 +    make_data_comparison_plot("met[1]", cutOSOF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_OF_inclusive");
961 +    
962 +    make_data_comparison_plot("met[4]", cutnJets&&cutOSSF,60, 0., 300.,-1., true, "pfMET", "met_SF_exclusive");
963 + //    make_data_comparison_plot("met[1]", cutnJets&&cutOSSF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_SF_exclusive");
964 +    make_data_comparison_plot("met[4]", cutnJets&&cutOSOF,60, 0., 300.,-1., true, "pfMET", "met_OF_exclusive");
965 + //    make_data_comparison_plot("met[1]", cutnJets&&cutOSOF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_OF_exclusive");
966 + */
967 +  
968 +    string basiccutsaved = (const char*)basiccut;
969 +    
970 +    string newbasiccut=basiccutsaved;
971 +    
972 +    size_t found;
973 +    found = newbasiccut.find( "pt2>20" );
974 +    while (found!=string::npos){
975 +      newbasiccut.replace( found, string( "pt2>20" ).length(), "pt2>10" );
976 +      found = newbasiccut.find( "pt2>20" );
977 +    }
978 + //    basiccut=TCut(newbasiccut.c_str());
979 + //    make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[4]>150")&&cutOSSF,60, 0., 300.,-1., false, "mll", "mll_SF_Aachen_pt2010_met");
980 + //    make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[4]>150")&&cutOSOF,60, 0., 300.,-1., false, "mll", "mll_OF_Aachen_pt2010_met");
981 + //    make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[1]>150")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_Aachen_pt2010_t1cpfmet");
982 + //    make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[1]>150")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_Aachen_pt2010_t1cpfmet");
983 + //    basiccut=TCut(basiccutsaved.c_str());
984 +    
985 + //     make_data_comparison_plot("pfJetGoodNum40", cutOSSF, 8, 0., 8.,-1., true, "#(jets)", "njets_SF");
986 + //     make_data_comparison_plot("pfJetGoodNum40", cutOSOF, 8, 0., 8.,-1., true, "#(jets)", "njets_OF");
987 +    
988 + //     make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[4]>100")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_ETH_met");
989 + //     make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[4]>100")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_ETH_met");
990 + //     make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[1]>100")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_ETH_t1cpfmet");
991 + //     make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[1]>100")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_ETH_t1cpfmet");
992 +
993 +    
994 + /*     make_data_comparison_plot("mll", cutOSSF,60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive");
995 +     make_data_comparison_plot("mll", cutOSSF&&TCut("id1==0"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive_ee");
996 +     make_data_comparison_plot("mll", cutOSSF&&TCut("id1==1"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive_mm");
997 +     make_data_comparison_plot("mll", cutOSOF,60, 15., 315.,-1., true, "m_{ll}", "mll_OF_inclusive");
998 +    make_data_comparison_plot("mll", cutnJets&&cutOSSF&&TCut("id1==0"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive_ee");
999 +    make_data_comparison_plot("mll", cutnJets&&cutOSSF&&TCut("id1==1"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive_mm");
1000 +    make_data_comparison_plot("mll", cutnJets&&cutOSSF,60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive");
1001 +    make_data_comparison_plot("mll", cutnJets&&cutOSOF,60, 15., 315.,-1., true, "m_{ll}", "mll_OF_exclusive");
1002 + */
1003 +    make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut,           massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig");
1004 +    make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut,           massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig");
1005 +    
1006 +    make_data_comparison_plot("mll", TCut("abs(eta1)<1.4&&abs(eta2)<1.4") && cutOSSF&&cutSignal&&mllCut,           massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_central");
1007 +    make_data_comparison_plot("mll", TCut("abs(eta1)<1.4&&abs(eta2)<1.4") && cutOSOF&&cutSignal&&mllCut,           massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_central");
1008 +    
1009 +    make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_btagVeto");
1010 +    make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_btagVeto");
1011 +    
1012 +    make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0",  massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_AtLeastOneBJet");
1013 +    make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0",  massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_AtLeastOneBJet");
1014 +
1015   /*
1016 <  make_OFSF_plots("mll", "met[4]>100", 60, 20., 320., false, "m_{ll}", "mll", false, false);
1017 <  make_OFSF_plots("pfJetGoodNum40", "met[4]>100", 7, 3, 10, true, "#(jets)", "njets", false, false);
1018 <  make_OFSF_plots("pt1",  "met[4]>100", 30, 0., 300., true, "p_{T,1}", "pt1", false, false);
1019 <  make_OFSF_plots("pt2",  "met[4]>100", 22, 0., 220., true, "p_{T,2}", "pt2", false, false);
1020 <  make_OFSF_plots("eta1", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{1}", "eta1", false, false, 0.15);
1021 <  make_OFSF_plots("eta2", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{2}", "eta2", false, false, 0.15);
1022 <  make_OFSF_plots("phi1", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", "phi1", false, false, 0.2);
1023 <  make_OFSF_plots("phi2", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", "phi2", false, false, 0.2);
683 <  make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", "met[4]>100", 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", "jpt1pt2", false, false, 0.2);
684 <  make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", "dphij1met", false, false, 0.0);
685 <  make_OFSF_plots("TMath::Abs(dphi)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", "dphi", false, false, 0.2);
686 <  make_OFSF_plots("TMath::Abs(dphiMet1)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", "dphiMet1", false, false, 0.2);
687 <  make_OFSF_plots("TMath::Abs(dphiMet2)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", "dphiMet2", false, false, 0.2);
688 <  make_OFSF_plots("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(l,MET)|)", "dphilc", false, false, 0.2);
689 <  make_OFSF_plots("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(jet,MET)|)", "dphijc", false, false, 0.2);
690 <  make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(l,MET)|)", "dphilco", false, false, 0.2);
691 <  make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(pfJetDphiMet[0])), TMath::Min( (TMath::Pi()-TMath::Abs(pfJetDphiMet[1])), (TMath::Pi() - TMath::Abs(pfJetDphiMet[2]))))", "met[4]>100", 16, 0, 3.2, false, "Min(#Pi - |#Delta#phi(jet,MET)|)", "dphijco", false, false, 0.2);
692 <
693 <  make_kin_plot("met[4]","",70,0,350,dolog,"MET [GeV]","met",doPF,true);
694 <  make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll",doPF,true);
695 <  make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof",doPF,true,true);
696 <  make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee",doPF,true);
697 <  make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm",doPF,true);
698 <  make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive",doPF,true);
699 <  make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_ee",doPF,true);
700 <  make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_mm",doPF,true);
701 <  make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_osof",doPF,true);
702 <  make_kin_plot("mll","",(int)((350-mll_low))/5,mll_low,350,dolog,"m_{ll} [GeV]","mll_inclusive_highrange",doPF);
703 <  if(!doPF) make_special_mll_plot((int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]");
704 <
705 <  make_kin_plot("pfJetGoodPt[0]/pfJetGoodPt[1]","",45,1,10,dolog,"pt_{j}^{1}/pt_{j}^{2}","j1j2ratio",doPF,true);
706 <  make_kin_plot("TMath::Abs(pfJetDphiMet[0])","",32,0,3.2,nolog,"|#Delta#phi(jet1,MET)|","dphiJ1MET",doPF,true);
707 <
708 <  make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets",doPF);
709 <  make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof",doPF);
710 <  make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_nocuts_except_mll_ossf",doPF);
711 <
712 <  make_kin_plot("numVtx","",(int)(30.5-(-0.5)),-0.5,30.5,nolog,"N(Vtx)","numVtx",doPF);
713 < //  make_kin_plot("jetpt[0]","",40,0,200,dolog,"leading jet p_{T} [GeV]","pfJetGoodPt_0",doPF);
714 < //  make_kin_plot("jeteta[0]","",40,-5,5,nolog,"leading jet #eta","pfJetGoodEta_0",doPF);
1016 >    make_data_comparison_plot("mll", mllCut&&cutOSOF&&cutSignal, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig");
1017 >    make_data_comparison_plot("mll", cutmass&&cutOSSF&&"met[4]>100&&met[4]<150&&pfJetGoodNum40==2", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_CR");
1018 >    make_data_comparison_plot("mll", cutmass&&cutOSOF&&"met[4]>100&&met[4]<150&&pfJetGoodNum40==2", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_CR");
1019 >
1020 >    make_data_comparison_plot("pfJetGoodNum40", cutOSSF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(jets)", "njets_SF_sig");
1021 >    make_data_comparison_plot("pfJetGoodNum40", cutOSOF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(jets)", "njets_OF_sig");
1022 >    make_data_comparison_plot("pfJetGoodNumBtag30", cutOSSF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(b-jets)", "nbjets_SF_sig");
1023 >    make_data_comparison_plot("pfJetGoodNumBtag30", cutOSOF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(b-jets)", "nbjets_OF_sig");
1024    */
1025 <  make_kin_plot("pt","",50,0,500,dolog,"Z p_{T} [GeV]","Zpt",doPF);
1026 <  make_kin_plot("pt1","",50,0,200,nolog,"p_{T} [GeV]","pt1",doPF);
1027 <  make_kin_plot("pt2","",50,0,200,nolog,"p_{T} [GeV]","pt2",doPF);
1028 <  make_kin_plot("eta1","",40,-3,3,nolog,"#eta_{l}","eta",doPF);
1029 <  make_kin_plot("jzb[1]","",100,-150,200,dolog,"JZB [GeV]","jzb_ossf",doPF);
1030 <  stringstream jzbcut;
1031 <  jzbcut << "((is_data&&("<<datajzb<<")>100)||(!is_data&&("<<mcjzb<<")>100))";
1032 <  make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB100",doPF,true);
1033 <  make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB100",doPF,true);
1034 <  make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB100",doPF,true);
1035 <  make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB100",doPF,true);
1036 <  stringstream jzbcut2;
1037 <  jzbcut2 << "((is_data&&("<<datajzb<<")>150)||(!is_data&&("<<mcjzb<<")>150))";
1038 <  make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB150",doPF,true);
1039 <  make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB150",doPF,true);
1040 <  make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB150",doPF,true);
1041 <  make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB150",doPF,true);
1042 <  stringstream jzbcut3;
1043 <  jzbcut3 << "((is_data&&("<<datajzb<<")>50)||(!is_data&&("<<mcjzb<<")>50))";
1044 <  make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB50",doPF,true);  
1045 <  make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB50",doPF,true,true);
1046 <  make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB50",doPF,true);  
1047 <  make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB50",doPF,true);  
1048 <
1049 <  make_kin_plot("mll","met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV] (MET>100GeV)","mll_met100_ll",doPF,true);
1050 <  //make_kin_plot("mll","met[4]>150&&id1==0",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ee} [GeV] (MET>150GeV)","mll_met150_ee",doPF,true);
1051 <  //make_kin_plot("mll","met[4]>150&&id1==1",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{#mu#mu} [GeV] (MET>150GeV)","mll_met150_mm",doPF,true);
1052 <
1053 <  make_special_obs_pred_mll_plot(datajzb,mcjzb,0);
1054 <  make_special_obs_pred_mll_plot(datajzb,mcjzb,50);
1055 <  make_special_obs_pred_mll_plot(datajzb,mcjzb,80);
1056 <  make_special_obs_pred_mll_plot(datajzb,mcjzb,100);
1025 >  }
1026 >
1027 >
1028 >  if ( doKin ) {
1029 >    string mllCut("");
1030 > /*  
1031 >    make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_0j",doPF,true);
1032 >    make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_1j",doPF,true);
1033 >    make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_2j",doPF,true);
1034 >    make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_3j",doPF,true);
1035 >
1036 >    make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_0j",doPF,true);
1037 >    make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_1j",doPF,true);
1038 >    make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_2j",doPF,true);
1039 >    make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_3j",doPF,true);
1040 > */
1041 >
1042 >    make_plain_kin_plot("pfTightHT",Cut2Str(cutOSOF&&cutnJets),120,0,1200,dolog,"H_{T}","HT_3jets_OF",doPF,true);
1043 >    make_plain_kin_plot("pfTightHT",Cut2Str(cutOSSF&&cutnJets) ,120,0,1200,dolog,"H_{T}","HT_3jets_SF",doPF,true);
1044 >    make_plain_kin_plot("pfTightHT",Cut2Str(cutOSOF&&cutnJets&&TCut("met[4]>100")),120,0,1200,dolog,"H_{T}","HT_3jets_OF_MET100",doPF,true);
1045 >    make_plain_kin_plot("pfTightHT",Cut2Str(cutOSSF&&cutnJets&&TCut("met[4]>100")) ,120,0,1200,dolog,"H_{T}","HT_3jets_SF_MET100",doPF,true);
1046 >
1047 >    
1048 >    make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_0j",doPF,true);
1049 >    make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_1j",doPF,true);
1050 >    make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_2j",doPF,true);
1051 >    make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_3j",doPF,true);
1052 >
1053 >    make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_0j",doPF,true);
1054 >    make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_1j",doPF,true);
1055 >    make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_2j",doPF,true);
1056 >    make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_3j",doPF,true);
1057 >
1058 >     /*make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&cutnJets&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_SF_20to70",doPF,true);
1059 >     make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&cutnJets&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_OF_20to70",doPF,true);
1060 >    
1061 >     make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&cutnJets&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_SF_above120",doPF,true);
1062 >     make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&cutnJets&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_OF_above120",doPF,true);*/
1063 >    
1064 > /*     make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&TCut("pfJetGoodNum40==2")&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_SF_20to70___2jetsonly",doPF,true);
1065 >     make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("pfJetGoodNum40==2")&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_OF_20to70___2jetsonly",doPF,true);
1066 >    
1067 >     make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&TCut("pfJetGoodNum40==2")&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_SF_above120___2jetsonly",doPF,true);
1068 >     make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("pfJetGoodNum40==2")&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_OF_above120___2jetsonly",doPF,true);
1069 > */
1070 >    write_info(__FUNCTION__,"Last plot");
1071 >    
1072 >    
1073 >    
1074 >    
1075 >    
1076 >    
1077 >    
1078 >    
1079 >     return;
1080 >    
1081 >    
1082 >    
1083 >     make_kin_plot("pfJetGoodPt[0]/pfJetGoodPt[1]","",45,1,10,dolog,"pt_{j}^{1}/pt_{j}^{2}","j1j2ratio",doPF,true);
1084 >     make_kin_plot("TMath::Abs(pfJetDphiMet[0])","",32,0,3.2,nolog,"|#Delta#phi(jet1,MET)|","dphiJ1MET",doPF,true);
1085 >    // Plots in signal region
1086 > //     make_kin_plot("met[4]","",70,0,350,dolog,"MET [GeV]","met",doPF,true);
1087 >     make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll",doPF,true);
1088 >     make_kin_plot("mll","mll>20",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof",doPF,true,true);
1089 >     make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee",doPF,true);
1090 >     make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm",doPF,true);
1091 >    
1092 >     make_kin_plot("pfJetGoodNum40",mllCut,9,-0.5,8.5,dolog,"nJets","nJets",doPF);
1093 >     make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof",doPF);
1094 >    
1095 > //     make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_met100",doPF,true);
1096 > //     make_kin_plot("mll","mll>20&&met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof_met100",doPF,true,true);
1097 > //     make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee_met100",doPF,true);
1098 > //     make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm_met100",doPF,true);
1099 >    
1100 >    
1101 > //     make_kin_plot("pfJetGoodNum40","met[4]>100"+mllCut,9,-0.5,8.5,dolog,"nJets","nJets_met100",doPF);
1102 > //     make_kin_plot("pfJetGoodNum40","met[4]>100",9,-0.5,8.5,dolog,"nJets","nJets_osof_met100",doPF);
1103 >
1104 >    // Further inclusive invariant mass plots
1105 >    make_kin_plot("mll",mllCut,(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive",doPF,true);
1106 > //     make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_ee",doPF,true);
1107 > //     make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_mm",doPF,true);
1108 >    make_kin_plot("mll",mllCut,(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_osof",doPF,true);
1109 >    
1110 >    //make_kin_plot("mll","",(int)((350-mll_low))/5,mll_low,350,dolog,"m_{ll} [GeV]","mll_inclusive_highrange",doPF);
1111 >    //if(!doPF) make_special_mll_plot((int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]");
1112 >    
1113 >    
1114 >    // Number of jets
1115 >    make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_inclusive",doPF);
1116 >    make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof_inclusive",doPF);
1117 >    //make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_nocuts_except_mll_ossf",doPF);
1118 >    
1119 >    // Others
1120 >    make_kin_plot("numVtx","",(int)(30.5-(-0.5)),-0.5,30.5,nolog,"N(Vtx)","numVtx",doPF);
1121 >    //  make_kin_plot("jetpt[0]","",40,0,200,dolog,"leading jet p_{T} [GeV]","pfJetGoodPt_0",doPF);
1122 >    //  make_kin_plot("jeteta[0]","",40,-5,5,nolog,"leading jet #eta","pfJetGoodEta_0",doPF);
1123 >    make_kin_plot("pt","",50,0,500,dolog,"Z p_{T} [GeV]","Zpt",doPF);
1124 >    make_kin_plot("pt1","",50,0,200,nolog,"p_{T} [GeV]","pt1",doPF);
1125 >    make_kin_plot("pt2","",50,0,200,nolog,"p_{T} [GeV]","pt2",doPF);
1126 >    make_kin_plot("eta1","",40,-3,3,nolog,"#eta_{l}","eta",doPF);
1127 >    make_kin_plot("jzb[1]","",100,-150,200,dolog,"JZB [GeV]","jzb_ossf",doPF);
1128 >    //   stringstream jzbcut;
1129 >    //   jzbcut << "((is_data&&("<<datajzb<<")>100)||(!is_data&&("<<mcjzb<<")>100))";
1130 >    //   make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB100",doPF,true);
1131 >    //   make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB100",doPF,true);
1132 >    //   make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB100",doPF,true);
1133 >    //   make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB100",doPF,true);
1134 >    //   stringstream jzbcut2;
1135 >    //   jzbcut2 << "((is_data&&("<<datajzb<<")>150)||(!is_data&&("<<mcjzb<<")>150))";
1136 >    //   make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB150",doPF,true);
1137 >    //   make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB150",doPF,true);
1138 >    //   make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB150",doPF,true);
1139 >    //   make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB150",doPF,true);
1140 >    //   stringstream jzbcut3;
1141 >    //   jzbcut3 << "((is_data&&("<<datajzb<<")>50)||(!is_data&&("<<mcjzb<<")>50))";
1142 >    //   make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB50",doPF,true);  
1143 >    //   make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB50",doPF,true,true);
1144 >    //   make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB50",doPF,true);  
1145 >    //   make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB50",doPF,true);  
1146 >    
1147 >    //   make_kin_plot("mll","met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV] (MET>100GeV)","mll_met100_ll",doPF,true);
1148 >    //make_kin_plot("mll","met[4]>150&&id1==0",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ee} [GeV] (MET>150GeV)","mll_met150_ee",doPF,true);
1149 >    //make_kin_plot("mll","met[4]>150&&id1==1",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{#mu#mu} [GeV] (MET>150GeV)","mll_met150_mm",doPF,true);
1150 >  }
1151 >    
1152 > //   make_special_obs_pred_mll_plot(datajzb,mcjzb,0);
1153 > //   make_special_obs_pred_mll_plot(datajzb,mcjzb,50);
1154 > //   make_special_obs_pred_mll_plot(datajzb,mcjzb,80);
1155 >   make_special_obs_pred_mll_plot(datajzb,mcjzb,100);
1156   //   make_special_obs_pred_mll_plot(datajzb,mcjzb,150);
1157   //   make_special_obs_pred_mll_plot(datajzb,mcjzb,200);
1158   //   make_special_obs_pred_mll_plot(datajzb,mcjzb,250);
1159  
1160 <  make_JES_plot(cutmass&&cutOSSF&&basiccut,"_ossf");
1161 <  make_JES_plot(cutmass&&cutOSOF&&basiccut,"_osof");
754 <
755 <  string regioncut[4] =    {"mll>0&&","mll>20&&mll<70&&","mll>70&&mll<110&&","mll>110&&"};
756 <  string regionname[4] = {"global","lowmass","Zregion","highmass"};
757 <  bool StoreRatioPlots=false;
758 <  for(int iregion=0;iregion<4;iregion++) {
759 <    cout << regionname[iregion] << " : " << regioncut[iregion] << endl;
760 <
761 <    //Z properties
762 <    make_kin_plot("pt","met[4]>100&&id1==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutOSSF&&cutnJets),40,0,200,dolog,"p_{T}^{ee} [GeV] (MET>100GeV)","Z/ZPT_ee_met100_"+regionname[iregion],doPF,StoreRatioPlots);
763 <    make_kin_plot("pt","met[4]>100&&id1==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutOSSF&&cutnJets),40,0,200,dolog,"p_{T}^{#mu#mu} [GeV] (MET>100GeV)","Z/ZPT_mm_met100_"+regionname[iregion],doPF,StoreRatioPlots);
764 <    
765 <    make_kin_plot("phi","met[4]>100&&id1==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutOSSF&&cutnJets),32,-3.2,3.2,nolog,"#phi^{ee} (MET>100GeV)","Z/ZPhi_ee_met100_"+regionname[iregion],doPF,StoreRatioPlots);
766 <    make_kin_plot("phi","met[4]>100&&id1==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutOSSF&&cutnJets),32,-3.2,3.2,nolog,"#phi^{#mu#mu} (MET>100GeV)","Z/ZPhi_mm_met100_"+regionname[iregion],doPF,StoreRatioPlots);
767 <    
768 <    make_kin_plot("eta","met[4]>100&&id1==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutOSSF&&cutnJets),30,-3.0,3.0,nolog,"#eta^{ee} (MET>100GeV)","Z/ZEta_ee_met100_"+regionname[iregion],doPF,StoreRatioPlots);
769 <    make_kin_plot("eta","met[4]>100&&id1==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutOSSF&&cutnJets),30,-3.0,3.0,nolog,"#eta^{#mu#mu} (MET>100GeV)","Z/ZEta_mm_met100_"+regionname[iregion],doPF,StoreRatioPlots);
770 <    
771 <    // lepton properties
772 <    make_kin_plot("pt1","met[4]>100&&id1==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),40,0,200,dolog,"p_{T}^{e1} [GeV] (MET>100GeV)","Lepton/LeadingElectronPt_met100_"+regionname[iregion],doPF,StoreRatioPlots);
773 <    make_kin_plot("pt1","met[4]>100&&id1==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),40,0,200,dolog,"p_{T}^{#mu1} [GeV] (MET>100GeV)","Lepton/LeadingMuonPt_met100_"+regionname[iregion],doPF,StoreRatioPlots);
774 <    make_kin_plot("pt2","met[4]>100&&id2==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),40,0,200,dolog,"p_{T}^{e2} [GeV] (MET>100GeV)","Lepton/SubLeadingElectronPt_met100_"+regionname[iregion],doPF,StoreRatioPlots);
775 <    make_kin_plot("pt2","met[4]>100&&id2==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),40,0,200,dolog,"p_{T}^{#mu2} [GeV] (MET>100GeV)","Lepton/SubLeadingMuonPt_met100_"+regionname[iregion],doPF,StoreRatioPlots);
776 <    
777 <    make_kin_plot("phi1","met[4]>100&&id1==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),32,-3.2,3.2,nolog,"#phi^{e1} (MET>100GeV)","Lepton/LeadingElectronPhi_met100_"+regionname[iregion],doPF,StoreRatioPlots);
778 <    make_kin_plot("phi1","met[4]>100&&id1==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),32,-3.2,3.2,nolog,"#phi^{#mu1} (MET>100GeV)","Lepton/LeadingMuonPhi_met100_"+regionname[iregion],doPF,StoreRatioPlots);
779 <    make_kin_plot("phi2","met[4]>100&&id2==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),32,-3.2,3.2,nolog,"#phi^{e2} (MET>100GeV)","Lepton/SubLeadingElectronPhi_met100_"+regionname[iregion],doPF,StoreRatioPlots);
780 <    make_kin_plot("phi2","met[4]>100&&id2==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),32,-3.2,3.2,nolog,"#phi^{#mu2} (MET>100GeV)","Lepton/SubLeadingMuonPhi_met100_"+regionname[iregion],doPF,StoreRatioPlots);
781 <    
782 <    make_kin_plot("eta1","met[4]>100&&id1==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),30,-3.0,3.0,nolog,"#eta^{e1} (MET>100GeV)","Lepton/LeadingElectronEta_met100_"+regionname[iregion],doPF,StoreRatioPlots);
783 <    make_kin_plot("eta1","met[4]>100&&id1==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),30,-3.0,3.0,nolog,"#eta^{#mu1} (MET>100GeV)","Lepton/LeadingMuonEta_met100_"+regionname[iregion],doPF,StoreRatioPlots);
784 <    make_kin_plot("eta2","met[4]>100&&id2==0&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),30,-3.0,3.0,nolog,"#eta^{e2} (MET>100GeV)","Lepton/SubLeadingElectronEta_met100_"+regionname[iregion],doPF,StoreRatioPlots);
785 <    make_kin_plot("eta2","met[4]>100&&id2==1&&"+regioncut[iregion]+(const char*)TCut(cutmass&&cutnJets),30,-3.0,3.0,nolog,"#eta^{#mu2} (MET>100GeV)","Lepton/SubLeadingMuonEta_met100_"+regionname[iregion],doPF,StoreRatioPlots);
786 <    
787 <    
788 <  }
1160 > //   make_JES_plot(cutmass&&cutOSSF&&basiccut,"_ossf");
1161 > //   make_JES_plot(cutmass&&cutOSOF&&basiccut,"_osof");
1162  
1163 +  switch_overunderflow(false);
1164   }
1165  
1166   void make_comp_plot( string var, string xlabel, string filename, float jzbcut, string mcjzb, string datajzb,
# Line 826 | Line 1200 | void make_comp_plot( string var, string
1200        if (ymax>ymin) datahisto->SetMaximum(ymax);
1201        lm3histo->SetLineStyle(2);
1202        datahisto->Draw("e1");
1203 <      mcstack.Draw("same");
1203 >      mcstack.Draw("histo,same");
1204        datahisto->Draw("same,e1");
1205        lm3histo->Draw("hist,same");
1206        TLegend *kinleg = allsamples.allbglegend((sRegions[iregion]+(ijzb?"neg":"pos")).c_str());
# Line 856 | Line 1230 | void region_comparison_plots(string mcjz
1230  
1231   //  switch_overunderflow(true);  // switching overflow/underflow bins on
1232      
1233 <  
1233 >  switch_overunderflow(true);
1234      flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- the arguments changed
1235    for(int ijzb=0;ijzb<(int)jzb_cuts.size();ijzb++) {
1236      float jzbcut=jzb_cuts[ijzb]; // Comparison plots are done for this JZB cut
# Line 1173 | Line 1547 | if you want to start from scratch (witho
1547  
1548   void do_prediction_plot(string jzb, TCanvas *globalcanvas, float high, int use_data, bool overlay_signal = false,string subdir="" )
1549   {
1550 < //  switch_overunderflow(true);
1550 >  
1551    bool is_data=false;
1552    bool use_signal=false;
1553    if(use_data==1) is_data=true;
1554    if(use_data==2) use_signal=true;
1555 <  int nbins=50;//100;
1556 <  if(is_data) nbins=50;
1555 >  int nbins=int(high/10);//100;
1556 >  if(is_data) nbins=int(high/10);
1557    float low=0;
1558 <  float hi=500;
1558 >  float hi=high;
1559 >  
1560 >  stringstream cutpositiveS;
1561 >  cutpositiveS << "(" << jzb << ">0)";
1562 >  TCut cutpositive(cutpositiveS.str().c_str());
1563 >  stringstream cutnegativeS;
1564 >  cutnegativeS << "(" << jzb << "<0)";
1565 >  TCut cutnegative(cutnegativeS.str().c_str());
1566 >  
1567    
1568    TH1F *blankback = new TH1F("blankback","blankback",int(high/10),0,high);
1569 <  TH1F *RcorrJZBeemm   = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1570 <  TH1F *LcorrJZBeemm   = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1571 <  TH1F *RcorrJZBem     = allsamples.Draw("RcorrJZBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1572 <  TH1F *LcorrJZBem     = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1569 >  TH1F *RcorrJZBeemm   = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1570 >  TH1F *LcorrJZBeemm   = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1571 >  TH1F *RcorrJZBem     = allsamples.Draw("RcorrJZBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1572 >  TH1F *LcorrJZBem     = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1573  
1574    blankback->GetXaxis()->SetTitle(RcorrJZBeemm->GetXaxis()->GetTitle());
1575    blankback->GetYaxis()->SetTitle(RcorrJZBeemm->GetYaxis()->GetTitle());
# Line 1202 | Line 1584 | void do_prediction_plot(string jzb, TCan
1584    
1585    TH1F *RcorrJZBeemmNoS;
1586  
1587 <    //these are for the ratio
1588 <  
1589 <  TH1F *JRcorrJZBeemm   = allsamples.Draw("JRcorrJZBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1590 <  TH1F *JLcorrJZBeemm   = allsamples.Draw("JLcorrJZBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1591 <  TH1F *JRcorrJZBem     = allsamples.Draw("JRcorrJZBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1210 <  TH1F *JLcorrJZBem     = allsamples.Draw("JLcorrJZBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1587 >  //these are for the ratio
1588 >  TH1F *JRcorrJZBeemm   = allsamples.Draw("JRcorrJZBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1589 >  TH1F *JLcorrJZBeemm   = allsamples.Draw("JLcorrJZBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1590 >  TH1F *JRcorrJZBem     = allsamples.Draw("JRcorrJZBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1591 >  TH1F *JLcorrJZBem     = allsamples.Draw("JLcorrJZBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1592    
1593    TH1F *JRcorrJZBSBem;
1594    TH1F *JLcorrJZBSBem;
1595    TH1F *JRcorrJZBSBeemm;
1596    TH1F *JLcorrJZBSBeemm;
1597    
1598 <  if(use_data==2 || overlay_signal) RcorrJZBeemmNoS = allsamples.Draw("RcorrJZBeemmNoS",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,false);
1598 >  if(use_data==2 || overlay_signal) RcorrJZBeemmNoS = allsamples.Draw("RcorrJZBeemmNoS",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,false);
1599  
1600    
1601    if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
1602 <    RcorrJZBSBem   = allsamples.Draw("RcorrJZBSBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1603 <    LcorrJZBSBem   = allsamples.Draw("LcorrJZBSBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1604 <    RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1605 <    LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1602 >    RcorrJZBSBem   = allsamples.Draw("RcorrJZBSBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1603 >    LcorrJZBSBem   = allsamples.Draw("LcorrJZBSBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1604 >    RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1605 >    LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1606      
1607      //these are for the ratio
1608 <    JRcorrJZBSBem   = allsamples.Draw("JRcorrJZBSBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1609 <    JLcorrJZBSBem   = allsamples.Draw("JLcorrJZBSBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1610 <    JRcorrJZBSBeemm = allsamples.Draw("JRcorrJZBSBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1611 <    JLcorrJZBSBeemm = allsamples.Draw("JLcorrJZBSBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1608 >    JRcorrJZBSBem   = allsamples.Draw("JRcorrJZBSBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1609 >    JLcorrJZBSBem   = allsamples.Draw("JLcorrJZBSBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1610 >    JRcorrJZBSBeemm = allsamples.Draw("JRcorrJZBSBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1611 >    JLcorrJZBSBeemm = allsamples.Draw("JLcorrJZBSBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
1612    }
1613    
1614    TH1F *lm4RcorrJZBeemm;
1615 <  if(overlay_signal || use_data == 2 || use_data == 1) lm4RcorrJZBeemm = allsamples.Draw("lm4RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,allsamples.FindSample("LM"));
1615 >  if(overlay_signal || use_data == 2 || use_data == 1) lm4RcorrJZBeemm = allsamples.Draw("lm4RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,allsamples.FindSample("LM"));
1616    
1617    flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- prediction changed.
1618    
# Line 1295 | Line 1676 | void do_prediction_plot(string jzb, TCan
1676    SQRT(BpredSys);
1677    BpredSys->Divide(JBpred);
1678  
1298 –  
1679    flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- prediction changed
1680    TH1F *Tpred = (TH1F*)RcorrJZBem->Clone("Bpred");
1681    Tpred->Add(LcorrJZBem,-1.0);
# Line 1348 | Line 1728 | void do_prediction_plot(string jzb, TCan
1728      Bpred->Draw("hist,same");
1729      //analytical_function[0]->Draw("same"); analytical_function[1]->Draw("same");analytical_function[2]->Draw("same");
1730      RcorrJZBeemm->Draw("e1x0,same");//HAVE IT ON TOP!
1731 <    lm4RcorrJZBeemm->Draw("hist,same");
1731 >    //lm4RcorrJZBeemm->Draw("hist,same");
1732      legBpred->AddEntry(RcorrJZBeemm,"observed","p");
1733      legBpred->AddEntry(Bpred,"predicted","l");
1734   //    legBpred->AddEntry(analytical_function[1],"predicted fit","l");
1735   //    legBpred->AddEntry(analytical_function[2],"stat. uncert.","l");
1736 <    legBpred->AddEntry(lm4RcorrJZBeemm,(allsamples.collection[allsamples.FindSample("LM")[0]].samplename).c_str(),"l");
1736 > //     legBpred->AddEntry(lm4RcorrJZBeemm,(allsamples.collection[allsamples.FindSample("LM")[0]].samplename).c_str(),"l");
1737      legBpred->Draw();
1738      DrawPrelim();
1739  
# Line 1363 | Line 1743 | void do_prediction_plot(string jzb, TCan
1743      Bpred->SetLineWidth(2);
1744      predcomppad->cd();
1745      predcomppad->SetLogy(1);
1366 –    
1746      TH1F *jzbnegative = (TH1F*)LcorrJZBeemm->Clone("jzbnegative");
1747      TH1F *sidebandsemu = (TH1F*)Bpred->Clone("sidebandsemu");
1748      sidebandsemu->Add(jzbnegative,-1);
# Line 1385 | Line 1764 | void do_prediction_plot(string jzb, TCan
1764   //    lm4RcorrJZBeemm->SetLineColor(kOrange+1);
1765      lm4RcorrJZBeemm->SetLineWidth(2);
1766      //lm4RcorrJZBeemm->SetLineWidth(2); // paper style. overruled.
1767 <    lm4RcorrJZBeemm->Draw("histo,same");
1767 > //     lm4RcorrJZBeemm->Draw("histo,same");
1768      DrawPrelim();
1769      TLegend *speciallegBpred = make_legend("",0.45,0.55);
1770      //TLegend *speciallegBpred = make_legend("",0.35,0.55); // paper style. overruled.
# Line 1395 | Line 1774 | void do_prediction_plot(string jzb, TCan
1774      if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) speciallegBpred->AddEntry(sidebandsemu,"Sidebands/e#mu (data)","f");
1775      else speciallegBpred->AddEntry(sidebandsemu,"e#mu (data)","f");
1776   //    speciallegBpred->AddEntry(lm4RcorrJZBeemmC,"LM4","l");
1777 <    speciallegBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
1777 > //     speciallegBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
1778      speciallegBpred->Draw();
1779      save_with_ratio(JRcorrJZBeemm,JBpred,predcomppad,subdir+"Bpred_Data_____PredictionComposition",true,true,"data/pred",BpredSys);
1780      
# Line 1404 | Line 1783 | void do_prediction_plot(string jzb, TCan
1783   //    THStack kostack = allsamples.DrawStack("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,!is_data, luminosity,use_signal);
1784      blankback->Draw();
1785   //    kostack.Draw("same");
1786 <    predcomposition.Draw();
1786 >    predcomposition.Draw("hist");
1787      Bpred->Draw("hist,same");
1788      //analytical_function[0]->Draw("same"); analytical_function[1]->Draw("same");analytical_function[2]->Draw("same");
1789      RcorrJZBeemm->Draw("e1x0,same");//HAVE IT ON TOP!
# Line 1431 | Line 1810 | void do_prediction_plot(string jzb, TCan
1810      legBpred->AddEntry(Bpred,"MC predicted","l");    
1811      if(versok) legBpred->AddEntry((TObject*)0,"",""); // Just for alignment // causes seg fault on root v5.18
1812      if(versok) legBpred->AddEntry((TObject*)0,"",""); // causes seg fault on root v5.18
1813 <    if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
1813 > //     if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
1814      legBpred->Draw();
1815      DrawMCPrelim();
1816      Bpredsaveas="Bpred_MC";
# Line 1463 | Line 1842 | void do_prediction_plot(string jzb, TCan
1842        
1843        RcorrJZBeemmNoS->SetLineStyle(2);
1844        legBpred2->AddEntry(RcorrJZBeemmNoS,"MC B","l");
1845 <      legBpred2->AddEntry(lm4RcorrJZBeemm,"MC S","l");
1845 > //       legBpred2->AddEntry(lm4RcorrJZBeemm,"MC S","l");
1846        legBpred2->Draw();
1847        RcorrJZBeemmNoS->SetLineColor(TColor::GetColor("#61210B"));
1848        RcorrJZBeemmNoS->Draw("histo,same");
# Line 1482 | Line 1861 | void do_prediction_plot(string jzb, TCan
1861    //save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,Bpredsaveas,true,use_data!=1,ytitle);
1862    save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,subdir+Bpredsaveas,true,false,ytitle,BpredSys);//not extending the y range anymore up to 4
1863  
1864 +
1865    
1866    flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
1867    // The part below is meaningless for the offpeak analysis (it's a comparison of the different estimates but there is but one estimate!)
# Line 1540 | Line 1920 | void do_prediction_plot(string jzb, TCan
1920        legBpredc->AddEntry(Bpredem,"MC OFZP","l");
1921        legBpredc->AddEntry(BpredSBem,"MC OFSB","l");
1922        legBpredc->AddEntry(BpredSBeemm,"MC SFSB","l");
1923 <      if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
1923 > //       if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
1924        legBpredc->Draw();
1925        CompleteSave(globalcanvas,subdir+"Bpred_MCwithS_comparison");
1926      }
# Line 1576 | Line 1956 | void do_prediction_plot(string jzb, TCan
1956      delete JLcorrJZBSBeemm;
1957    }
1958    if(overlay_signal || use_data==2) delete lm4RcorrJZBeemm;
1579 –  switch_overunderflow(false);
1959   }
1960  
1961   void do_prediction_plots(string mcjzb, string datajzb, float DataSigma, float MCSigma, bool overlay_signal ) {
1962 +  switch_overunderflow(true);
1963    TCanvas *globalcanvas = new TCanvas("globalcanvas","Prediction Canvas");
1964    do_prediction_plot(datajzb,globalcanvas,jzbHigh ,data,overlay_signal);
1965    if ( !PlottingSetup::Approved ) {
# Line 1588 | Line 1968 | void do_prediction_plots(string mcjzb, s
1968    } else {
1969      write_info(__FUNCTION__,"You set approved to true, therefore not producing prediction/observation plots for MC with and without signal.");
1970    }
1971 < }
1972 <
1593 < void do_ratio_plot(int is_data,vector<float> binning, string jzb, TCanvas *can, float high=-9999) {
1594 <  bool do_data=0;
1595 <  bool dosignal=0;
1596 <  if(is_data==1) do_data=1;
1597 <  if(is_data==2) dosignal=1;
1598 <  TH1F *RcorrJZBeemm   = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,do_data, luminosity,dosignal);
1599 <  TH1F *LcorrJZBeemm   = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,do_data, luminosity,dosignal);
1600 <  TH1F *RcorrJZBem     = allsamples.Draw("RcorrJZBem",jzb.c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,do_data, luminosity,dosignal);
1601 <  TH1F *LcorrJZBem     = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,do_data, luminosity,dosignal);
1602 <  
1603 <  TH1F *RcorrJZBSBem;
1604 <  TH1F *LcorrJZBSBem;
1605 <  TH1F *RcorrJZBSBeemm;
1606 <  TH1F *LcorrJZBSbeemm;
1607 <  
1608 <  if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
1609 <    RcorrJZBSBem   = allsamples.Draw("RcorrJZBSbem",jzb.c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,do_data, luminosity,dosignal);
1610 <    LcorrJZBSBem   = allsamples.Draw("LcorrJZBSbem",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,do_data, luminosity,dosignal);
1611 <    RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSbeemm",jzb.c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,do_data, luminosity,dosignal);
1612 <    LcorrJZBSbeemm = allsamples.Draw("LcorrJZBSbeemm",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,do_data, luminosity,dosignal);
1613 <  }
1614 <  
1615 <
1616 <
1617 <
1618 <  flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
1619 <  TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred");
1620 <  if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
1621 <    Bpred->Add(RcorrJZBem,1.0/3);
1622 <    Bpred->Add(LcorrJZBem,-1.0/3);
1623 <    Bpred->Add(RcorrJZBSBem,1.0/3);
1624 <    Bpred->Add(LcorrJZBSBem,-1.0/3);
1625 <    Bpred->Add(RcorrJZBSBeemm,1.0/3);
1626 <    Bpred->Add(LcorrJZBSbeemm,-1.0/3);
1627 <  } else {
1628 <    Bpred->Add(RcorrJZBem,1.0);
1629 <    Bpred->Add(LcorrJZBem,-1.0);
1630 <  }
1631 <
1632 <  can->cd();
1633 <  can->SetLogy(0);
1634 <  Bpred->SetLineColor(kRed);
1635 <  Bpred->SetStats(0);
1636 <  if(high>0) Bpred->GetXaxis()->SetRangeUser(0,high);
1637 <  TH1F *JZBratioforfitting=(TH1F*)RcorrJZBeemm->Clone("JZBratioforfitting");
1638 <  JZBratioforfitting->Divide(Bpred);
1639 <  TGraphAsymmErrors *JZBratio = histRatio(RcorrJZBeemm,Bpred,is_data,binning,false);
1640 <  
1641 <  
1642 <  JZBratio->SetTitle("");
1643 <  JZBratio->GetYaxis()->SetRangeUser(0.0,9.0);
1644 < //  if(is_data==1) JZBratio->GetXaxis()->SetRangeUser(0,jzbHigh );
1645 <  
1646 <  TF1 *pol0 = new TF1("pol0","[0]",0,1000);
1647 <  TF1 *pol0d = new TF1("pol0","[0]",0,1000);
1648 < //  straightline_fit->SetParameter(0,1);
1649 <  JZBratioforfitting->Fit(pol0,"Q0R","",0,30);
1650 <  pol0d->SetParameter(0,pol0->GetParameter(0));
1651 <  
1652 <  JZBratio->GetYaxis()->SetTitle("Observed/Predicted");
1653 <  JZBratio->GetXaxis()->SetTitle("JZB [GeV]");
1654 <  if ( high>0 ) JZBratio->GetXaxis()->SetRangeUser(0.0,high);
1655 <  JZBratio->GetYaxis()->SetNdivisions(519);
1656 <  JZBratio->GetYaxis()->SetRangeUser(0.0,4.0);
1657 <  JZBratio->GetYaxis()->CenterTitle();
1658 <  JZBratio->GetXaxis()->CenterTitle();
1659 <  JZBratio->SetMarkerSize(DataMarkerSize);
1660 <  JZBratio->Draw("AP");
1661 <  /////----------------------------
1662 <  TPaveText *writeresult = new TPaveText(0.15,0.78,0.49,0.91,"blNDC");
1663 <  writeresult->SetFillStyle(4000);
1664 <  writeresult->SetFillColor(kWhite);
1665 <  writeresult->SetTextFont(42);
1666 <  writeresult->SetTextSize(0.03);
1667 <  writeresult->SetTextAlign(12);
1668 <  ostringstream jzb_agreement_data_text;
1669 <  jzb_agreement_data_text<< setprecision(2) << "mean =" << pol0->GetParameter(0) << " #pm " << setprecision(1) <<  pol0->GetParError(0);
1670 <  if(is_data==1) fitresultconstdata=pol0->GetParameter(0);// data
1671 <  if(is_data==0) fitresultconstmc=pol0->GetParameter(0); // monte carlo, no signal
1672 < /*  if(is_data) writeresult->AddText("Data closure test");
1673 <  else writeresult->AddText("MC closure test");
1674 <  */
1675 <  writeresult->AddText(jzb_agreement_data_text.str().c_str());
1676 < //  writeresult->Draw("same");
1677 < //  pol0d->Draw("same");
1678 <  TF1 *topline =  new TF1("","1.5",0,1000);
1679 <  TF1 *bottomline =  new TF1("","0.5",0,1000);
1680 <  topline->SetLineColor(kBlue);
1681 <  topline->SetLineStyle(2);
1682 <  bottomline->SetLineColor(kBlue);
1683 <  bottomline->SetLineStyle(2);
1684 < //  topline->Draw("same");
1685 < //  bottomline->Draw("same");
1686 <  TF1 *oneline = new TF1("","1.0",0,1000);
1687 <  oneline->SetLineColor(kBlue);
1688 <  oneline->SetLineStyle(1);
1689 <  oneline->Draw("same");
1690 <  TLegend *phony_leg = make_legend("ratio",0.6,0.55,false);//this line is just to have the default CMS Preliminary (...) on the canvas as well.
1691 <  if(is_data==1) DrawPrelim();
1692 <  else DrawMCPrelim();
1693 <  TLegend *leg = new TLegend(0.55,0.75,0.89,0.89);
1694 <  leg->SetTextFont(42);
1695 <  leg->SetTextSize(0.04);
1696 < //   if(is_data==1) leg->SetHeader("Ratio (data)");
1697 < //   else leg->SetHeader("Ratio (MC)");
1698 <
1699 <  TString MCtitle("MC ");
1700 <  if (is_data==1) MCtitle = "";
1701 <
1702 <  leg->SetFillStyle(4000);
1703 <  leg->SetFillColor(kWhite);
1704 <  leg->SetTextFont(42);
1705 < //  leg->AddEntry(topline,"+20\% sys envelope","l");
1706 <  leg->AddEntry(JZBratio,MCtitle+"obs / "+MCtitle+"pred","p");
1707 <  leg->AddEntry(oneline,"ratio = 1","l");
1708 < //  leg->AddEntry(pol0d,"fit in [0,30] GeV","l");
1709 < //  leg->AddEntry(bottomline,"#pm50% envelope","l");
1710 <
1711 <
1712 <  //leg->Draw("same"); // no longer drawing legend
1713 <
1714 <  if(is_data==1) CompleteSave(can, "jzb_ratio_data");
1715 <  if(is_data==0) CompleteSave(can, "jzb_ratio_mc");
1716 <  if(is_data==2) CompleteSave(can, "jzb_ratio_mc_BandS");//special case, MC with signal!
1717 <  
1718 <  delete RcorrJZBeemm;
1719 <  delete LcorrJZBeemm;
1720 <  delete RcorrJZBem;
1721 <  delete LcorrJZBem;
1722 <  
1723 <  delete RcorrJZBSBem;
1724 <  delete LcorrJZBSBem;
1725 <  delete RcorrJZBSBeemm;
1726 <  delete LcorrJZBSbeemm;
1727 < }
1728 <
1729 < void do_ratio_plots(string mcjzb,string datajzb,vector<float> ratio_binning) {
1730 <  TCanvas *globalc = new TCanvas("globalc","Ratio Plot Canvas");
1731 <  globalc->SetLogy(0);
1732 <  
1733 <  do_ratio_plot(mc,ratio_binning,mcjzb,globalc, jzbHigh  );
1734 <  do_ratio_plot(data,ratio_binning,datajzb,globalc, jzbHigh  );
1735 <  do_ratio_plot(mcwithsignal,ratio_binning,mcjzb,globalc, jzbHigh  );
1971 >  delete globalcanvas;
1972 >  switch_overunderflow(false);
1973   }
1974  
1975   string give_jzb_expression(float peak, int type) {
# Line 1920 | Line 2157 | void draw_pure_jzb_histo(TCut cut,string
2157    datahisto->SetMinimum(0.1);
2158    datahisto->SetMarkerSize(DataMarkerSize);
2159    datahisto->Draw("e1");
2160 <  mcstack.Draw("same");
2160 >  mcstack.Draw("histo,same");
2161    datahisto->Draw("same,e1");
2162    
2163    TLegend *leg;
# Line 1951 | Line 2188 | void draw_pure_jzb_histo(TCut cut,string
2188    datahisto->SetMinimum(0.1);
2189    datahisto->SetMarkerSize(DataMarkerSize);
2190    datahisto->Draw("e1");
2191 <  mcstack.Draw("same");
2191 >  mcstack.Draw("histo,same");
2192    datahisto->Draw("same,e1");
2193    leg->SetHeader("");
2194    leg->Draw();
# Line 2000 | Line 2237 | void produce_stretched_jzb_plots(string
2237    
2238  
2239   void diboson_plots(string mcjzb, string datajzb,vector<float> ratio_binning) {
2240 +  switch_overunderflow(true);
2241    vector<int> SamplesToBeModified = allsamples.FindSampleBySampleName("WW/WZ/ZZ");
2242    
2243    if(SamplesToBeModified.size()==0 || SamplesToBeModified[0]==-1) {
# Line 2039 | Line 2277 | void diboson_plots(string mcjzb, string
2277      dout << "     Reset xs for sample " << (allsamples.collection)[SamplesToBeModified[i]].samplename << " from " << Upxs << " to " << (allsamples.collection)[SamplesToBeModified[i]].xs << " (by a factor of " << stretchfactor << ") and reset the correct name (from " << Upname << ")" << endl;
2278      
2279    }
2280 <  
2280 > //   switch_overunderflow(false);
2281   }
2282  
2283  
# Line 2074 | Line 2312 | void draw_normalized_data_vs_data_histo(
2312  
2313  
2314   void jzb_plots(string mcjzb, string datajzb,vector<float> ratio_binning) {
2315 +  switch_overunderflow(true);
2316    TCanvas *can = new TCanvas("can","JZB Plots Canvas");
2317    float max=jzbHigh ;
2318    float min=-120;
# Line 2110 | Line 2349 | void jzb_plots(string mcjzb, string data
2349      if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_OFSB_coarse",can,coarse_binning);
2350    }
2351    delete can;
2352 +  switch_overunderflow(false);
2353   }
2354  
2355  
# Line 2650 | Line 2890 | void ttbar_sidebands_comparison(string m
2890    //in the case of the on peak analysis, we compare the 3 control regions to the real value
2891    //in the case of the OFF peak analysis, we compare our control region to the real value
2892    TCut weightbackup=cutWeight;
2893 +  switch_overunderflow(true);
2894    
2895    bool doPURW=false;
2896    
# Line 2671 | Line 2912 | void ttbar_sidebands_comparison(string m
2912        
2913    float simulatedlumi = luminosity; //in pb please - adjust to your likings
2914  
2915 <  TH1F *TZem  = systsamples.Draw("TZem",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2916 <  TH1F *nTZem = systsamples.Draw("nTZem","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2915 >  TH1F *TZem  = systsamples.Draw("TZem",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT"));
2916 >  TH1F *nTZem = systsamples.Draw("nTZem","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT"));
2917    TH1F *TSem;
2918    TH1F *nTSem;
2919 <  TH1F *TZeemm  = systsamples.Draw("TZeemm",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2920 <  TH1F *nTZeemm = systsamples.Draw("nTZeemm","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2919 >  TH1F *TZeemm  = systsamples.Draw("TZeemm",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT"));
2920 >  TH1F *nTZeemm = systsamples.Draw("nTZeemm","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT"));
2921    TH1F *TSeemm;
2922    TH1F *nTSeemm;
2923    
2924    if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
2925 <    TSem    = systsamples.Draw("TSem",       mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2926 <    nTSem   = systsamples.Draw("nTSem",  "-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2927 <    TSeemm  = systsamples.Draw("TSeemm",     mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2928 <    nTSeemm = systsamples.Draw("nTSeemm","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2925 >    TSem    = systsamples.Draw("TSem",       mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT"));
2926 >    nTSem   = systsamples.Draw("nTSem",  "-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT"));
2927 >    TSeemm  = systsamples.Draw("TSeemm",     mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT"));
2928 >    nTSeemm = systsamples.Draw("nTSeemm","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT"));
2929    }
2930  
2931    TCanvas *tcan = new TCanvas("tcan","tcan");
# Line 2716 | Line 2957 | void ttbar_sidebands_comparison(string m
2957      histos.push_back(TSem);
2958      histos.push_back(TSeemm);
2959    }
2960 <  draw_all_ttbar_histos(tcan,histos,"histo",1);
2960 >  draw_all_ttbar_histos(tcan,histos,"histo",8);
2961    
2962    TLegend *leg = make_legend("MC t#bar{t}",0.6,0.65,false);
2963    leg->AddEntry(TZeemm,"SFZP","l");
# Line 2728 | Line 2969 | void ttbar_sidebands_comparison(string m
2969    leg->Draw("same");
2970    DrawMCPrelim(simulatedlumi);
2971    CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_shape_comparison");
2972 +  
2973    TH1F *TZemcopy = (TH1F*)TZem->Clone("TZemcopy");
2974    TH1F *TZeemmcopy = (TH1F*)TZeemm->Clone("TZeemmcopy");
2975    TH1F *TSeemmcopy;
# Line 2793 | Line 3035 | void ttbar_sidebands_comparison(string m
3035    DrawMCPrelim(simulatedlumi);
3036    CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_shape_comparison_ratio");
3037    
3038 +  
3039    if (0) { // Turn this off: we don't need this
3040  
3041      ///-------------- second part: only look at the quantity we actually care about!
# Line 2898 | Line 3141 | void ttbar_sidebands_comparison(string m
3141  
3142    delete tcan;
3143    cutWeight=weightbackup;
3144 +  switch_overunderflow(false);
3145   }
3146  
3147   void ttbar_sidebands_comparison(string mcjzb, vector<float> jzb_binning) {
# Line 2932 | Line 3176 | void ttbar_sidebands_comparison(string m
3176    nicer_binning.push_back(150);
3177    //nicer_binning.push_back(175);
3178    nicer_binning.push_back(200);
3179 <  nicer_binning.push_back(250);
3180 <  nicer_binning.push_back(300);
3181 <  nicer_binning.push_back(400);
3179 > //  nicer_binning.push_back(250);
3180 > //  nicer_binning.push_back(300);
3181 > //  nicer_binning.push_back(400);
3182    
3183    ttbar_sidebands_comparison(mcjzb,nicer_binning, "ttbar/");
3184   }
3185  
3186  
3187 < void zjets_prediction_comparison(string mcjzbWithPUa) {
3187 > void zjets_prediction_comparison(string mcjzbWithPUa, TCut massregioncut, string massregionname) {
3188 >  cout << "****************************************************************************************************************************************************************" << endl;
3189    TCanvas *zcan = new TCanvas("zcan","zcan");  
3190   //  zcan->SetLogy(1);
3191    TCut weightbackup=cutWeight;
3192    
3193 <  bool UsePURW=false;
3193 >  bool UsePURW=true;
3194    
3195    
3196    string mcjzb;
3197    if(UsePURW) {
3198      mcjzb=mcjzbWithPUa;
3199 +    cout << "Using PURW peak positions" << endl;
3200    } else {
3201      // Do it without PU re-weighting
3202 +    cutWeight="1.0";
3203      float MCPeakNoPU=0,MCPeakErrorNoPU=0,DataPeakNoPU=0,DataPeakErrorNoPU=0,MCSigma=0,DataSigma=0;
3204      stringstream resultsNoPU;
3205      stringstream noPUdatajzb;
# Line 2963 | Line 3210 | void zjets_prediction_comparison(string
3210      
3211      mcjzb = noPUmcjzb.str();
3212      
2966 –    cutWeight="1.0";
3213    }
3214  
3215    
# Line 2987 | Line 3233 | void zjets_prediction_comparison(string
3233    string var( "abs("+mcjzb+")" );
3234    
3235    TCut notTau("abs(genMID1)!=15");
3236 <  TCut ONLYTau("mll>20");
3236 >  TCut ee_mm_tautau("mll>0");
3237    
3238  
3239 <  TH1F *hJZBpos = systsamples.Draw("hJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3240 <  TH1F *hJZBneg = systsamples.Draw("hJZBneg",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3239 >  TH1F *hJZBpos = systsamples.Draw("hJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3240 >  TH1F *hJZBneg = systsamples.Draw("hJZBneg",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3241    
3242    hJZBpos->SetLineColor(kBlack);
3243    hJZBneg->SetLineColor(kRed);
# Line 3006 | Line 3252 | void zjets_prediction_comparison(string
3252    leg->AddEntry(hJZBneg,"Predicted","l");
3253    leg->Draw("same");
3254    DrawMCPrelim(simulatedlumi);
3255 <  CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction");
3255 >  CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction"+any2string(massregionname));
3256    
3257    TH1F* hratio = (TH1F*)hJZBpos->Clone("hratio");
3258    hratio->Divide(hJZBneg);
# Line 3015 | Line 3261 | void zjets_prediction_comparison(string
3261      cout << "Positive: " << hJZBpos->GetBinContent(i) << "   vs    Negative : " << hJZBneg->GetBinContent(i) << "   (ratio : " << hJZBpos->GetBinContent(i) / hJZBneg->GetBinContent(i) << endl;
3262    }
3263    
3264 <  zcan->SetLogy(0);
3265 <  hratio->GetYaxis()->SetRangeUser(0,2.0);
3264 > //  zcan->SetLogy(0);
3265 >  hratio->GetYaxis()->SetRangeUser(0,2.5);
3266    hratio->GetYaxis()->SetTitle("Observed/Predicted");
3267    hratio->Draw("e1");
3268    
# Line 3025 | Line 3271 | void zjets_prediction_comparison(string
3271    TLine *bottom = new TLine(binning[0],0.75,binning[binning.size()-1],0.75);
3272    
3273  
3274 +  cout << __LINE__ << endl;
3275    top->SetLineColor(kBlue);top->SetLineStyle(2);
3276    bottom->SetLineColor(kBlue);bottom->SetLineStyle(2);
3277    center->SetLineColor(kBlue);
# Line 3038 | Line 3285 | void zjets_prediction_comparison(string
3285    leg2->AddEntry(bottom,"syst. envelope","l");
3286    leg2->Draw("same");
3287    DrawMCPrelim(simulatedlumi);
3288 <  CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction_ratio");
3288 >  CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction_ratio"+any2string(massregionname));
3289    
3290    TCut reducedNJets(cutnJets);
3291    
3292 <  TH1F *TAUhJZBpos = systsamples.Draw("TAUhJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3293 <  TH1F *LcorrJZBeemm     = systsamples.Draw("LcorrJZBeemm",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3294 <  TH1F *RcorrJZBem       = systsamples.Draw("RcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3295 <  TH1F *LcorrJZBem       = systsamples.Draw("LcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3292 >  TH1F *TAUhJZBpos       = systsamples.Draw("TAUhJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3293 >  TH1F *LcorrJZBeemm     = systsamples.Draw("LcorrJZBeemm",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3294 >  TH1F *RcorrJZBem       = systsamples.Draw("RcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3295 >  TH1F *LcorrJZBem       = systsamples.Draw("LcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3296  
3297 +  cout << __LINE__ << endl;
3298    TH1F *RcorrJZBSBem;
3299    TH1F *LcorrJZBSBem;
3300    TH1F *RcorrJZBSBeemm;
3301    TH1F *LcorrJZBSBeemm;
3302  
3303 +  cout << __LINE__ << endl;
3304    if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
3305 <    RcorrJZBSBem   = systsamples.Draw("RcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3306 <    LcorrJZBSBem   = systsamples.Draw("LcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3307 <    RcorrJZBSBeemm = systsamples.Draw("RcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3308 <    LcorrJZBSBeemm = systsamples.Draw("LcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3305 >    RcorrJZBSBem   = systsamples.Draw("RcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3306 >    LcorrJZBSBem   = systsamples.Draw("LcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3307 >    RcorrJZBSBeemm = systsamples.Draw("RcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3308 >    LcorrJZBSBeemm = systsamples.Draw("LcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY"));
3309    }
3310    
3311    TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred");
# Line 3074 | Line 3323 | void zjets_prediction_comparison(string
3323    
3324    Bpred->SetLineColor(kRed);
3325  
3326 +  cout << __LINE__ << endl;
3327    TAUhJZBpos->SetLineColor(kBlack);
3328    Bpred->SetLineColor(kRed);
3329    
# Line 3087 | Line 3337 | void zjets_prediction_comparison(string
3337    TAUleg->AddEntry(Bpred,"Predicted","l");
3338    TAUleg->Draw("same");
3339    DrawMCPrelim(simulatedlumi);
3340 <  CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction");
3340 >  CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction__"+any2string(massregionname));
3341    
3342    TH1F* TAUhratio = (TH1F*)TAUhJZBpos->Clone("TAUhratio");
3343    TAUhratio->Divide(Bpred);
# Line 3105 | Line 3355 | void zjets_prediction_comparison(string
3355    center->Draw("same");
3356    bottom->Draw("same");
3357    
3358 <  TLegend *TAUleg2 = make_legend("MC Z+jets #rightarrow #tau#tau",0.25,0.75,false);
3358 >  TLegend *TAUleg2 = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.25,0.75,false);
3359    TAUleg2->AddEntry(TAUhratio,"obs / pred","pe");
3360    TAUleg2->AddEntry(bottom,"syst. envelope","l");
3361    TAUleg2->Draw("same");
3362    DrawMCPrelim(simulatedlumi);
3363 <  CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction_ratio");
3363 >  CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction_ratio"+any2string(massregionname));
3364    
3365    delete Bpred;
3366    delete TAUhJZBpos;
3367    delete LcorrJZBeemm;
3368    delete RcorrJZBem;
3369    delete LcorrJZBem;
3370 +  delete hJZBpos;
3371 +  delete hJZBneg;
3372    
3373    if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
3374      delete RcorrJZBSBem;
# Line 3129 | Line 3381 | void zjets_prediction_comparison(string
3381    delete zcan;
3382    cutWeight=weightbackup;
3383  
3384 +  cout << __LINE__ << endl;
3385   }
3386  
3387  
3388 <
3388 > void zjets_prediction_comparison(string mcjzbWithPUa) {
3389 >  zjets_prediction_comparison(mcjzbWithPUa, TCut(""), "nomasscut");
3390 >  zjets_prediction_comparison(mcjzbWithPUa, TCut("abs(mll-91)<20"), "Zwindow_20");
3391 >  zjets_prediction_comparison(mcjzbWithPUa, TCut("mll>20&&mll<70"), "LowMassRegion2070");
3392 >  zjets_prediction_comparison(mcjzbWithPUa, TCut("mll>110"), "HighMassRegion110");
3393 > }
3394 >  
3395   void sideband_assessment(string datajzb, float min=30.0, float max=50.0) {
3396    tout << endl << endl;
3397    stringstream bordercut;

Diff Legend

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