Changeset 504 in svn for trunk

Timestamp:

Jul 22, 2009, 7:20:03 PM (15 years ago)

Author:

Xavier Rouby

Message:

added ATLAS expectation band for etmis

File:

• trunk/routines/resolutions_atlas.C (modified) (6 diffs)

trunk/routines/resolutions_atlas.C

@@ -39 +39 @@
 #include "TF1.h"
 #include "TGraph.h"
+#include "TMultiGraph.h"
+#include "TGraphErrors.h"
 #include "TLegend.h"
 
@@ -395 +397 @@
   TProfile *ETSoverET = new TProfile("ETSoverET","ETmis resol",numBin-1,bins,-1000,1000);
   TF1 *fitfun = new TF1("user","[0]*sqrt(x)",0,max);
-  fitfun->SetLineColor(kBlue);
+  fitfun->SetLineColor(kBlack);
+  fitfun->SetLineStyle(2);
  
  
@@ -435 +438 @@
   TCanvas *c6 = new TCanvas("c6","ETmis resolution",100,100,600,450);
   c6->cd();
+
+  TF1 *fdo = new TF1("fdo","0.53*sqrt(x)",1,900);
+  TF1 *fup = new TF1("fup","0.57*sqrt(x)",1,900);
+  int N = 100;
+  float xx[N], yy[N], eex[N], eey[N];
+  for (int i=0; i<N; i++) {
+     xx[i]  = 9*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;
@@ -441 +458 @@
   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]");
@@ -448 +464 @@
   gr11->Fit("user","RQ");
   gr11->Fit("user","RQ");
-  gr11->GetYaxis()->SetRangeUser(0,60);
-  gr11->GetXaxis()->SetRangeUser(1,max);
-  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->GetXaxis()->SetRangeUser(1,900);
+  mgr->SetMaximum(16);
 
   Double_t* params = fitfun->GetParameters(); 
@@ -468 +490 @@
   TPaveText *atlas = MakeTPave(0.65,0.25,0.8,0.30,"ATLAS-like detector");
   atlas->Draw();
+
+  TLegend *legend = new TLegend(0.69,0.6,0.94,0.74,NULL,"NDC");
+  legend->AddEntry(gr,"ATLAS","f");
+  legend->AddEntry(gr11,"Delphes","l");
+  legend->SetFillColor(10);
+  legend->SetBorderSize(0);
+  legend->Draw();
+
   
   delete fitfun;