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

Timestamp:

Nov 7, 2008, 6:37:38 PM (16 years ago)

Author:

severine ovyn

Message:

Working on the resolution

File:

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

trunk/routines/resolutions.C

@@ -9 +9 @@
 #include "interface/FuncDef.h"
 
-void Plot()
+void JetResol(TTree *Analyze)
 {
-
-  setTDRStyle();
-  gROOT->Reset();
-
-  TFile *f1 = new TFile("Smearing.root","read");
-  TTree *Analyze = (TTree*)f1->Get("Analysis");
-
-
-//$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
-//$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
-
-   TCanvas *c1 = new TCanvas("c1","gerrors2",100,100,600,450);
+   TCanvas *c1 = new TCanvas("c1","JET resol",100,100,600,450);
    c1->cd(); 
   
-   const Int_t numBin=5;
-//   double bins[numBin]={5,10,25,55,85,115,205,400};
-   double bins[5]={10,30,50,80,120};
-   TProfile *ETSoverET = new TProfile("ETSoverET","Resolution des electrons",(numBin-1),bins,-10,10);
+   const Int_t numBin=7;
+   double bins[numBin]={0,20,60,100,220,420,860};
+   TProfile *ETSoverET = new TProfile("ETSoverET","Jet resolution: E_{T}^{rec}/E_{T}^{mc}",(numBin-1),bins,-10,10);
    Analyze->Draw("JetPTResol.SmearePT:JetPTResol.NonSmearePT>>ETSoverET","JetPTResol.NonSmearePT>1");
 
    double rms[numBin];
    double mean[numBin];
-   //TF1 *Gauss = new TF1("Gauss","gaus",0,4);
 
-   TCanvas *c2 = new TCanvas("c2","gerrors2",0,0,1000,650);
+   TCanvas *c2 = new TCanvas("c2","JET resol",0,0,1000,650);
    c2->cd(); int frame=0;
    c2->Divide(6,2);
-//   c2->Divide(6,2);
 
-   float x[numBin];
-   float y[numBin];
+   float x[numBin-1];
+   float y[numBin-1];
    float finval=0;//valeur a remplir puis a fitter
 
-   for ( int i=0; i<numBin; i++)  // premiÚre bin : i ==1 et pas i == 0
+   for ( int i=0; i<(numBin-1); i++)  // premiÚre bin : i ==1 et pas i == 0
       {
         TAxis *xaxis = ETSoverET->GetXaxis();
@@ -64 +50 @@
         GaussValues(Analyze,tempName,rms[i],mean[i],mystringMin,mystringMax);
         x[i]=binCenter;
-        mean[i]=ETSoverET->GetBinContent(i+1);
-   //     rms[i]=ETSoverET->GetBinError(i+1);
         finval=rms[i]/mean[i];
-        y[i]=finval;
+        y[i]=(finval*2);
         cout<<"finval "<<finval<<" mean val "<<mean[i]<<" rms value "<<rms[i]<<endl;
 
       }
 
-y[0]=0.2;
-
-   TCanvas *c3 = new TCanvas("c3","gerrors2",100,100,600,450);
+   TCanvas *c3 = new TCanvas("c3","JET resol",100,100,600,450);
    c3->cd();
 
-   TF1 *fitfun = new TF1("user","sqrt(pow([0],2)+pow([1]/sqrt(x),2)+pow([2]/x,2))",10,500);
+   TF1 *fitfun = new TF1("user","sqrt(pow([0],2)+pow([1]/sqrt(x),2)+pow([2]/x,2))",1,3000);
 
-   TGraph *gr11 = new TGraph((numBin),x,y);
+   TGraph *gr11 = new TGraph((numBin-1),x,y);
    gr11->Draw("AP");
    gr11->GetXaxis()->SetTitle("E_{T}^{MC} [GeV]");
@@ -96 +78 @@
 // delete Gauss;
 }
+
+
+void ETmisResol(TTree *Analyze)
+{
+
+   TCanvas *c4 = new TCanvas("c4","ETmis resol",100,100,600,450);
+   c4->cd(); 
+  
+   const Int_t numBin=10;
+   double bins[numBin]={0,10,20,40,80,160,320,640,1080,2160};
+   TProfile *ETSoverET = new TProfile("ETSoverET","ETmis resol",(numBin-1),bins,-10,10);
+   Analyze->Draw("(ETmisResol.EtSmeare-ETmisResol.Et):(ETmisResol.Et)>>ETSoverET","(ETmisResol.Et)>1");
+
+   double rms[numBin];
+   double mean[numBin];
+
+   TCanvas *c5 = new TCanvas("c5","gerrors2",0,0,1000,650);
+   c5->cd(); int frame=0;
+   c5->Divide(6,2);
+
+   float x[numBin-1];
+   float y[numBin-1];
+
+   for ( int i=0; i<(numBin-1); i++)  // premiÚre bin : i ==1 et pas i == 0
+      {
+        TAxis *xaxis = ETSoverET->GetXaxis();
+        float binCenter = xaxis->GetBinCenter(i+1);
+        int binMin = (int)xaxis->GetBinLowEdge(i+1);
+        int binMax = (int)xaxis->GetBinUpEdge(i+1);
+        cout<<"binMin "<<binMin<<" binMax "<<binMax<<" bin center "<<binCenter<<endl;
+        char tempMin[500];
+        sprintf(tempMin,"(ETmisResol.Et)>%d",binMin);
+        string mystringMin(tempMin);
+        char tempMax[500];
+        sprintf(tempMax,"(ETmisResol.Et)<%d",binMax);
+        string mystringMax(tempMax);
+        char tempName[500];
+        sprintf(tempName,"(ETmisResol.EtSmeare-ETmisResol.Et)>>hETSoverET%d",i);
+        string mystringName(tempName);
+
+        c5->cd(++frame);
+        GaussValues(Analyze,tempName,rms[i],mean[i],mystringMin,mystringMax);
+        x[i]=binCenter;
+        mean[i]=ETSoverET->GetBinContent(i+1);
+        y[i]=rms[i];
+        cout<<" mean val "<<mean[i]<<" rms value "<<rms[i]<<endl;
+
+      }
+
+   TCanvas *c6 = new TCanvas("c6","ETmis resolution",100,100,600,450);
+   c6->cd();
+
+   TF1 *fitfun = new TF1("user","sqrt(pow([0],2)+pow([1]/sqrt(x),2)+pow([2]/x,2))",1,3000);
+
+   TGraph *gr11 = new TGraph((numBin-1),x,y);
+   gr11->Draw("AP");
+   gr11->GetXaxis()->SetTitle("E_{T}^{MC} [GeV]");
+   gr11->GetYaxis()->SetTitle("#sigma(E_{T}^{rec}>");
+
+   //fitfun->FixParameter(0,0.0541930);
+   //fitfun->FixParameter(1,0.634908);
+   gr11->Fit("user","RQ");
+   gr11->Fit("user","RQ");
+   gr11->Fit("user","RQ");
+   gr11->Fit("user","R");
+
+
+
+ delete fitfun;
+// delete Gauss;
+}
+
+void General()
+{
+  setTDRStyle();
+  gROOT->Reset();
+
+  TFile *f1 = new TFile("Smearing.root","read");
+  TTree *Analyze = (TTree*)f1->Get("Analysis");
+  JetResol(Analyze);
+
+}