Changeset 60 in svn

Timestamp:

Nov 27, 2008, 8:53:09 PM (16 years ago)

Author:

severine ovyn

Message:

final version

File:

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

trunk/routines/resolutions.C

@@ -82 +82 @@
   if(params[2]==0.000000)sprintf(tempResol,"#frac{#sigma(E_{T}^{rec}/E_{T}^{MC})}{<E_{T}^{rec}/E_{T}^{MC}>} = #frac{%f}{#sqrt{E_{T}^{MC}}} #oplus %f",params[1]/100,params[0]/100);
   else sprintf(tempResol,"#frac{#sigma(E_{T}^{rec}/E_{T}^{MC})}{<E_{T}^{rec}/E_{T}^{MC}>} = #frac{%f}{#sqrt{E_{T}^{MC}}} #oplus #frac{%f}{E_{T}^{MC}} #oplus %f",params[1]/100,params[2]/100,params[0]/100);
-  
+ 
   TPaveText *leg1 = MakeTPave(0.4,0.6,0.8,0.65,tempResol);
   leg1->Draw();
 
+  TPaveText *events = MakeTPave(0.2,0.75,0.35,0.8,"Events: pp #rightarrow gg ");
+  events->Draw();
+
+
   TPaveText *Delphes = MakeTPave(0.2,0.15,0.35,0.2,"MG/ME + Delphes");
   Delphes->Draw();
@@ -102 +106 @@
   TTree *Analyze = (TTree*)f1->Get("Analysis");
   
-  const Int_t numBin=9;
-  double bins[numBin]={0,20,40,60,80,100,150,220,400};
+  const Int_t numBin=6;
+  double bins[numBin]={0,20,40,60,80,140};
   TProfile *ETSminusET = new TProfile("ETSminusET","Electron resolution: E_{T}^{rec}/E_{T}^{mc}",(numBin-1),bins,-10,10);
   
@@ -173 +177 @@
   Delphes->Draw();
   
-  TPaveText *events = MakeTPave(0.2,0.75,0.35,0.8,"t#bar{t} #rightarrow l^{+}l^{-}#nu#bar{#nu}b#bar{b}");
+  TPaveText *events = MakeTPave(0.2,0.75,0.35,0.8,"Events: WHq'#rightarrow W#tau#tauq'#rightarrowjjl#tauq', m_{H}=150 GeV ");
   events->Draw();
   
@@ -188 +192 @@
   TTree *Analyze = (TTree*)f1->Get("Analysis");
   
-  TCanvas *ct = new TCanvas("ct","Tau information",0,0,1000,425);
-  ct->cd(); int frame=0;
-  ct->Divide(2,1);
-  
-  ct->cd(++frame);
+  TCanvas *ct1 = new TCanvas("ct1","Tau information",100,100,600,450);
+  ct1->cd(); 
+  
   TH1F *tauEnergy =MakeNormTH1F(20,0.8,1,Analyze,"TauJetPTResol.EnergieCen>>tauEnergy",1, 0, 1,2,false);
   tauEnergy->Draw();
@@ -204 +206 @@
   TPaveText *ctau = MakeTPave(0.3,0.75,0.45,0.8,"C_{#tau} = #frac{#sum E^{towers} (#DeltaR = 0.15)}{E^{jet}}");
   ctau->Draw();
-  TPaveText *events1 = MakeTPave(0.3,0.65,0.45,0.7,"Events: t#bar{t} #rightarrow l^{+}#nul^{-}#bar{#nu}b#bar{b}");
+  TPaveText *events1 = MakeTPave(0.3,0.65,0.45,0.7,"Events: WHq'#rightarrow W#tau#tauq'#rightarrowjjl#tauq', m_{H}=150 GeV ");
   events1->Draw();
   
   
-  ct->cd(++frame);
+  TCanvas *ct2 = new TCanvas("ct2","Tau information",100,100,600,450);
+  ct2->cd();
   TH1F *NumTrack =MakeNormTH1F(6,0,6,Analyze,"TauJetPTResol.NumTrack>>NumTrack",1, 0, 1,2,false);
   NumTrack->Draw();
@@ -219 +222 @@
   TPaveText *numtracks = MakeTPave(0.6,0.75,0.85,0.8,"#DeltaR < 0.4, p_{T}^{track} > 2 GeV");
   numtracks->Draw();
-  TPaveText *events = MakeTPave(0.6,0.65,0.85,0.7,"Events: t#bar{t} #rightarrow l^{+}#nul^{-}#bar{#nu}b#bar{b}");
+  TPaveText *events = MakeTPave(0.6,0.65,0.85,0.7,"Events: WHq'#rightarrow W#tau#tauq'#rightarrowjjl#tauq', m_{H}=150 GeV ");
   events->Draw();
    
@@ -230 +233 @@
   gROOT->Reset();
   
-  TFile *f1 = new TFile("test2.root","read");
+  TFile *f1 = new TFile("PTMIS.root","read");
   TTree *Analyze = (TTree*)f1->Get("Analysis");
   
   TF1 *fitfun = new TF1("user","[0]*sqrt(x)",0,600);
   const Int_t numBin=7;
-//  double bins[numBin]={10,100,150,200,250,300,350,400,450,500};
   double bins[numBin]={200,250,300,350,400,450,500};
-//  double bins[numBin]={0,10,20,30,40,50,60,70};
   TProfile *ETSoverET = new TProfile("ETSoverET","ETmis resol",(numBin-1),bins,-1000,1000);
-  Analyze->Draw("ETmisResol.ExSmeare:ETmisResol.SEt>>ETSoverET");
-  Analyze->Fit("user","RQ");
-  Analyze->Fit("user","RQ");
  
  
@@ -313 +311 @@
 void General()
 {
-  //JetResol();
-  //ElecResol();
+  JetResol();
+  ElecResol();
   ETmisResol();
-  //TauJetInfo();
-  
-}
-
-
+  TauJetInfo();
+  
+}
+
+