Changeset 503 in svn for trunk

Timestamp:

Jul 22, 2009, 7:19:56 PM (15 years ago)

Author:

Xavier Rouby

Message:

added CMS expectation band for etmis

File:

• trunk/routines/resolutions.C (modified) (7 diffs)

trunk/routines/resolutions.C

@@ -39 +39 @@
 #include "TF1.h"
 #include "TGraph.h"
+#include "TMultiGraph.h"
+#include "TGraphErrors.h"
 #include "TLegend.h"
 
@@ -68 +70 @@
   float ex[numBin-1];
   float ey[numBin-1];
+  float erec_over_etruth[numBin-1];
   
   float finval=0;//valeur a remplir puis a fitter
@@ -95 +98 @@
       ex[i]=0;
       ey[i]=0;
+      erec_over_etruth[i]=mean[i]+1;
     }
-  
+ 
+
+  TCanvas *c3 = new TCanvas("c3","JET linearity",100,100,600,450);
+  c3->cd();
+  //for(int i=0; i<numBin-1; i++) cout << x[i] << "\t\t" << erec_over_etruth[i] << endl;
+  TGraph *linearity = new TGraph(numBin-1,x,erec_over_etruth);
+  linearity->Draw("AP");
+  linearity->SetTitle("");
+  linearity->GetXaxis()->SetTitle("E^{Truth} [GeV]");
+  linearity->GetYaxis()->SetTitle("<E^{Reco}/E^{Truth}>");
+  linearity->GetXaxis()->SetRangeUser(30,1000);
+  linearity->SetMinimum(0.85);
+  linearity->SetMaximum(1.20);
+  gPad->SetLogx();
+
+ 
   TCanvas *c2 = new TCanvas("c2","JET resol",100,100,600,450);
   c2->cd();
@@ -297 +316 @@
   
   TF1 *fitfun = new TF1("user","[0]*sqrt(x)",0,600);
-  fitfun->SetLineColor(kBlue);
+  fitfun->SetLineColor(kBlack);
+  fitfun->SetLineStyle(2);
   const Int_t numBin=6;
   double bins[numBin]={180,260,340,420,500,580};
@@ -338 +358 @@
   TCanvas *c6 = new TCanvas("c6","ETmis resolution",100,100,600,450);
   c6->cd();
+
+  TF1 *fdo = new TF1("fdo","0.6*sqrt(x)",1,600);
+  TF1 *fup = new TF1("fup","0.7*sqrt(x)",1,600);
+  int N = 100;
+  float xx[N], yy[N], eex[N], eey[N];
+  for (int i=0; i<N; i++) {
+     xx[i]  = 6*i;
+     yy[i]  = (fup->Eval(xx[i]) + fdo->Eval(xx[i]))/2.;
+     eex[i] = 0;
+     eey[i] = (fup->Eval(xx[i]) - fdo->Eval(xx[i]))/2.;
+  }
+  TGraphErrors *gr = new TGraphErrors(N,xx,yy,eex,eey);
+  gr->SetFillColor(18);
+
   
   x[numBin]=0;
@@ -344 +378 @@
   gr11->SetTitle("");
   gr11->SetMarkerStyle(3);
-  gr11->Draw("AP");
-  gr11->GetXaxis()->SetTitle("#Sigma E_{T} [GeV]");
-  gr11->GetYaxis()->SetTitle("Resolution of x-component  of MET [GeV]");
+  gr11->SetLineStyle(2);
   gr11->Fit("user","RQ");
   gr11->Fit("user","RQ");
   gr11->Fit("user","RQ");
   gr11->Fit("user","RQ");
-  gr11->GetYaxis()->SetRangeUser(0,60);
-  gr11->GetXaxis()->SetRangeUser(1,600);
-  gr11->SetMaximum(16);
+
+
+  TMultiGraph *mgr = new TMultiGraph("mgr","");
+  mgr->Add(gr,"3");
+  mgr->Add(gr11,"P");
+  mgr->Draw("A");
+  mgr->GetXaxis()->SetTitle("#Sigma E_{T} [GeV]");
+  mgr->GetYaxis()->SetTitle("Resolution of x-component  of MET [GeV]");
+  mgr->GetYaxis()->SetRangeUser(0,60);
+  mgr->GetXaxis()->SetRangeUser(1,600);
+  mgr->SetMaximum(16);
 
   Double_t* params = fitfun->GetParameters(); 
- 
+
   char tempResol[500];
   sprintf(tempResol,"E_{x}^{mis} resolution: %.2f #times #sqrt{E_{T}}",params[0]);
@@ -369 +409 @@
   Delphes->Draw();
  
-  TPaveText *atlas = MakeTPave(0.65,0.25,0.8,0.30,"CMS-like detector");
-  atlas->Draw();
+  TPaveText *cms = MakeTPave(0.65,0.25,0.8,0.30,"CMS-like detector");
+  cms->Draw();
+
+  TLegend *legend = new TLegend(0.69,0.6,0.94,0.74,NULL,"NDC");
+  legend->AddEntry(gr,"CMS","f");
+  legend->AddEntry(gr11,"Delphes","l");
+  legend->SetFillColor(10);
+  legend->SetBorderSize(0);
+  legend->Draw();
+
  
   delete fitfun;