Changeset 27 in svn for trunk/routines

Timestamp:

Nov 13, 2008, 9:43:35 AM (16 years ago)

Author:

severine ovyn

Message:

Jet resol OK

File:

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

trunk/routines/resolutions.C

@@ -9 +9 @@
 #include "interface/FuncDef.h"
 
-void JetResol(TTree *Analyze)
+void JetResol()
 {
-   TCanvas *c1 = new TCanvas("c1","JET resol",100,100,600,450);
-   c1->cd(); 
+
+   setTDRStyle();
+   gROOT->Reset();
+
+   TFile *f1 = new TFile("JETResol.root","read");
+   TTree *Analyze = (TTree*)f1->Get("Analysis");
   
    const Int_t numBin=7;
-   double bins[numBin]={0,20,60,100,220,420,860};
+   double bins[numBin]={0,20,40,60,80,100,220};
    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];
 
-   TCanvas *c2 = new TCanvas("c2","JET resol",0,0,1000,650);
-   c2->cd(); int frame=0;
-   c2->Divide(6,2);
+   TCanvas *c1 = new TCanvas("c1","JET resol",0,0,1000,650);
+   c1->cd(); int frame=0;
+   c1->Divide(6,2);
 
    float x[numBin-1];
    float y[numBin-1];
+   float ex[numBin-1];
+   float ey[numBin-1];
+
    float finval=0;//valeur a remplir puis a fitter
 
@@ -38 +44 @@
         cout<<"binMin "<<binMin<<" binMax "<<binMax<<" bin center "<<binCenter<<endl;
         char tempMin[500];
-        sprintf(tempMin,"JetPTResol.NonSmearePT>%d",binMin);
+        if(i==0)binMin=5;
+        sprintf(tempMin,"JetPTResol.NonSmearePT > %d",binMin);
         string mystringMin(tempMin);
         char tempMax[500];
-        sprintf(tempMax,"JetPTResol.NonSmearePT<%d",binMax);
+        sprintf(tempMax,"JetPTResol.NonSmearePT < %d",binMax);
         string mystringMax(tempMax);
         char tempName[500];
-        sprintf(tempName,"JetPTResol.SmearePT>>hETSoverET%d",i);
+        sprintf(tempName,"(JetPTResol.SmearePT)>>hETSoverET%d",i);
         string mystringName(tempName);
 
-        c2->cd(++frame);
+        c1->cd(++frame);
         GaussValues(Analyze,tempName,rms[i],mean[i],mystringMin,mystringMax);
         x[i]=binCenter;
         finval=rms[i]/mean[i];
-        y[i]=(finval*2);
+        y[i]=(finval*100);
+        ex[i]=0;
+        ey[i]=0;
         cout<<"finval "<<finval<<" mean val "<<mean[i]<<" rms value "<<rms[i]<<endl;
-
       }
 
-   TCanvas *c3 = new TCanvas("c3","JET resol",100,100,600,450);
-   c3->cd();
+   TCanvas *c2 = new TCanvas("c2","JET resol",100,100,600,450);
+   c2->cd();
 
-   TF1 *fitfun = new TF1("user","sqrt(pow([0],2)+pow([1]/sqrt(x),2)+pow([2]/x,2))",1,3000);
+   TF1 *fitfun = new TF1("user","sqrt(pow([0],2)+pow([1]/sqrt(x),2)+pow([2]/x,2))",1,170);
 
-   TGraph *gr11 = new TGraph((numBin-1),x,y);
+   TGraphErrors *gr11 = new TGraphErrors((numBin-1),x,y,ex,ey);
    gr11->Draw("AP");
+   gr11->SetTitle("");
    gr11->GetXaxis()->SetTitle("E_{T}^{MC} [GeV]");
-   gr11->GetYaxis()->SetTitle("#sigma(E_{T}^{rec}/E_{T}^{MC})/<E_{T}^{rec}/E_{T}^{MC}>");
+   gr11->GetYaxis()->SetRangeUser(0,50);
+   gr11->GetYaxis()->SetTitle("#sigma(E_{T}^{rec}/E_{T}^{MC})_{fit}/<E_{T}^{rec}/E_{T}^{MC}>_{fit}");
 
-   //fitfun->FixParameter(0,0.0541930);
-   //fitfun->FixParameter(1,0.634908);
-   gr11->Fit("user","RQ");
-   gr11->Fit("user","RQ");
-   gr11->Fit("user","RQ");
+   Double_t* params = fitfun->GetParameters();
+
    gr11->Fit("user","R");
+   gr11->Fit("user","RI");
+   gr11->Fit("user","RI");
+   gr11->Fit("user","RI");
+   char tempResol[500];
+   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();
 
-
- delete fitfun;
-// delete Gauss;
+   delete fitfun;
 }
 
 
-void ETmisResol(TTree *Analyze)
+void ETmisResol()
 {
 
-   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");
+   setTDRStyle();
+   gROOT->Reset();
 
-   double rms[numBin];
-   double mean[numBin];
+   TFile *f1 = new TFile("PTMIS.root","read");
+   TTree *Analyze = (TTree*)f1->Get("Analysis");
+
+   const Int_t numBin=6;
+   double bins[numBin]={10,20,30,50,80,100};
+   TProfile *ETSoverET = new TProfile("ETSoverET","ETmis resol",(numBin-1),bins,-1000,1000);
+
+   double rms[numBin-1];
+   double mean[numBin-1];
 
    TCanvas *c5 = new TCanvas("c5","gerrors2",0,0,1000,650);
@@ -98 +112 @@
    c5->Divide(6,2);
 
-   float x[numBin-1];
-   float y[numBin-1];
+   double x[numBin-1];
+   double y[numBin-1];
 
    for ( int i=0; i<(numBin-1); i++)  // premiÚre bin : i ==1 et pas i == 0
@@ -109 +123 @@
         cout<<"binMin "<<binMin<<" binMax "<<binMax<<" bin center "<<binCenter<<endl;
         char tempMin[500];
-        sprintf(tempMin,"(ETmisResol.Et)>%d",binMin);
+        sprintf(tempMin,"ETmisResol.Et>%d",binMin);
         string mystringMin(tempMin);
         char tempMax[500];
-        sprintf(tempMax,"(ETmisResol.Et)<%d",binMax);
+        sprintf(tempMax,"ETmisResol.Et<%d",binMax);
         string mystringMax(tempMax);
+
         char tempName[500];
-        sprintf(tempName,"(ETmisResol.EtSmeare-ETmisResol.Et)>>hETSoverET%d",i);
+        sprintf(tempName,"ETmisResol.EtSmeare>>hETSoverET%d",i);
         string mystringName(tempName);
 
         c5->cd(++frame);
-        GaussValues(Analyze,tempName,rms[i],mean[i],mystringMin,mystringMax);
+        GaussValuesETmis(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;
@@ -130 +144 @@
    c6->cd();
 
-   TF1 *fitfun = new TF1("user","sqrt(pow([0],2)+pow([1]/sqrt(x),2)+pow([2]/x,2))",1,3000);
+   TF1 *fitfun = new TF1("user","sqrt(pow([0],2)+pow([1]/sqrt(x),2)+pow([2]/x,2))",0,200);
 
-   TGraph *gr11 = new TGraph((numBin-1),x,y);
+   TGraph *gr11 = new TGraph((numBin-1),x,rms);
    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->GetXaxis()->SetTitle("MC^{MET} [GeV]");
+   gr11->GetYaxis()->SetTitle("#sigma(REC^{MET}-MC^{MET})");
    gr11->Fit("user","RQ");
    gr11->Fit("user","RQ");
@@ -144 +155 @@
    gr11->Fit("user","R");
 
+   TCanvas *c7 = new TCanvas("c7","ETmis resolution",100,100,600,450);
+   c7->cd();
+
+   TGraph *gr12 = new TGraph((numBin-1),x,mean);
+   gr12->Draw("AP");
+   gr12->GetXaxis()->SetTitle("MC^{MET} [GeV]");
+   gr12->GetYaxis()->SetTitle("<(REC^{MET}-MC^{MET})>");
 
 
  delete fitfun;
-// delete Gauss;
 }
 
 void General()
 {
-  setTDRStyle();
-  gROOT->Reset();
+  JetResol();
+  ETmisResol();
+}
 
-  TFile *f1 = new TFile("Smearing.root","read");
-  TTree *Analyze = (TTree*)f1->Get("Analysis");
-  JetResol(Analyze);
 
-}