Changeset 473 in svn for trunk/routines/resolutions.C

Timestamp:

Jul 13, 2009, 5:26:01 PM (15 years ago)

Author:

Xavier Rouby

Message:

update des labels / securites dans le code

File:

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

trunk/routines/resolutions.C

@@ -49 +49 @@
   
   TFile *f1 = new TFile("JET2.root","read");
+  if(!f1->IsOpen()) { cout << "could not open "<< f1->GetName() << ". Exiting..." << endl; return;}
+  if(!f1->FindKey("Analysis")) { cout << "Bad input file, could not find the \"Analysis\" tree. Exiting..." << endl; return;}
   TTree *Analyze = (TTree*)f1->Get("Analysis");
   
@@ -153 +155 @@
   
   TFile *f1 = new TFile("ETMIS2.root","read");
+  if(!f1->IsOpen()) { cout << "could not open "<< f1->GetName() << ". Exiting..." << endl; return;}
+  if(!f1->FindKey("Analysis")) { cout << "Bad input file, could not find the \"Analysis\" tree. Exiting..." << endl; return;}
   TTree *Analyze = (TTree*)f1->Get("Analysis");
   
@@ -240 +244 @@
   
   TFile *f1 = new TFile("TAUJET2.root","read");
+  if(!f1->IsOpen()) { cout << "could not open "<< f1->GetName() << ". Exiting..." << endl; return;}
+  if(!f1->FindKey("Analysis")) { cout << "Bad input file, could not find the \"Analysis\" tree. Exiting..." << endl; return;}
   TTree *Analyze = (TTree*)f1->Get("Analysis");
   
@@ -283 +289 @@
   gROOT->Reset();
   
-  TFile *f1 = new TFile("JET2.root","read");
+  //TFile *f1 = new TFile("JET2.root","read");
+  TFile *f1 = new TFile("JET2_cms_midpoint.root","read");
+  
+  if(!f1->IsOpen()) { cout << "could not open "<< f1->GetName() << ". Exiting..." << endl; return;}
+  if(!f1->FindKey("Analysis")) { cout << "Bad input file, could not find the \"Analysis\" tree. Exiting..." << endl; return;}
   TTree *Analyze = (TTree*)f1->Get("Analysis");
   
   TF1 *fitfun = new TF1("user","[0]*sqrt(x)",0,600);
+  fitfun->SetLineColor(kBlue);
   const Int_t numBin=6;
   double bins[numBin]={180,260,340,420,500,580};
@@ -330 +341 @@
   x[numBin]=0;
   y[numBin]=0;
-  TGraph *gr11 = new TGraph((numBin),x,rms);
+  TGraph *gr11 = new TGraph(numBin,x,rms);
+  gr11->SetTitle("");
+  gr11->SetMarkerStyle(3);
   gr11->Draw("AP");
-  gr11->GetXaxis()->SetTitle("Offline sum of E_{T} [GeV]");
+  gr11->GetXaxis()->SetTitle("#Sigma E_{T} [GeV]");
   gr11->GetYaxis()->SetTitle("Resolution of x-component  of MET [GeV]");
   gr11->Fit("user","RQ");
@@ -340 +353 @@
   gr11->GetYaxis()->SetRangeUser(0,60);
   gr11->GetXaxis()->SetRangeUser(1,600);
+  gr11->SetMaximum(16);
 
   Double_t* params = fitfun->GetParameters(); 
  
   char tempResol[500];
-  sprintf(tempResol,"%f * #sqrt{E_{T}}",params[0]);
-
-  TPaveText *leg1 = MakeTPave(0.4,0.6,0.8,0.65,tempResol);
+  sprintf(tempResol,"E_{x}^{mis} resolution: %.2f #times #sqrt{E_{T}}",params[0]);
+
+  TPaveText *leg1 = MakeTPave(0.26,0.84,0.43,0.89,tempResol);
   leg1->Draw();
 
-  TPaveText *leg2 = MakeTPave(0.2,0.8,0.8,0.85,"WHq'#rightarrow WWWq'#rightarrow lllq', m_{H}=150 GeV ");
+  TPaveText *leg2 = MakeTPave(0.15,0.75,0.45,0.8,"Events: pp #rightarrow ggX");
   leg2->Draw();
 
-  TPaveText *Delphes = MakeTPave(0.2,0.15,0.35,0.2,"MG/ME + Delphes");
+  TPaveText *Delphes = MakeTPave(0.65,0.19,0.8,0.24,"MG/ME + Pythia + Delphes");
   Delphes->Draw();
-  
-  
+ 
+  TPaveText *atlas = MakeTPave(0.65,0.25,0.8,0.30,"CMS-like detector");
+  atlas->Draw();
+ 
   delete fitfun;
 }