Changeset 473 in svn

Timestamp:

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

Author:

Xavier Rouby

Message:

update des labels / securites dans le code

Location:

trunk/routines

Files:

• resolutions.C (modified) (6 diffs)
• resolutions_atlas.C (modified) (9 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;
 }

trunk/routines/resolutions_atlas.C

@@ -55 +55 @@
   //TFile *f1 = new TFile("JET2_atlas_midpoint_newCaloRes.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");
   
@@ -124 +125 @@
   Double_t* params = fitfun->GetParameters();
 
-  fitfun->SetLineColor(1);
+  fitfun->SetLineColor(kBlue);
   gr11->Fit("user","QRN");
   gr11->Fit("user","QRN");
@@ -136 +137 @@
 
   TF1 *fitfunDelphes = new TF1("userDelphes",tempResol2,7,1000);
-  fitfunDelphes->SetLineColor(596);
+  fitfunDelphes->SetLineColor(kBlue);
   fitfunDelphes->SetLineStyle(7);
   fitfunDelphes->Draw("same");
 
-  fitfunATLAS->SetLineColor(1);
+  fitfunATLAS->SetLineColor(kRed);
   fitfunATLAS->SetLineWidth(2);
   fitfunATLAS->Draw("same");
@@ -168 +169 @@
   TFile *f1 = new TFile("ETMIS2_ATLAS.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");
   
@@ -256 +258 @@
   TFile *f1 = new TFile("TAUJET2_ATLAS.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");
   
@@ -299 +302 @@
   gROOT->Reset();
   
-  TFile *f1 = new TFile("JET2_ATLAS.root","read");
+  //TFile *f1 = new TFile("JET2_ATLAS.root","read");
+  TFile *f1 = new TFile("JET2_atlas_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};
@@ -333 +339 @@
       char tempName[500];
       sprintf(tempName,"ETmisResol.ExSmeare>>hETSoverET%d",i);
+      //sprintf(tempName,"ETmisResol.EtSmeare>>hETSoverET%d",i); // non-gaussian, so it does not work properly
       string mystringName(tempName);
       
@@ -347 +354 @@
   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");
@@ -357 +366 @@
   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,"ATLAS-like detector");
+  atlas->Draw();
   
   delete fitfun;