Changeset 27 in svn for trunk/interface

Timestamp:

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

Author:

severine ovyn

Message:

Jet resol OK

File:

• trunk/interface/FuncDef.h (modified) (3 diffs)

trunk/interface/FuncDef.h

@@ -72 +72 @@
 }
 
-void GaussValues(TTree * Analyze,string histo,double &rms, double &mean, string min, string max)
+void GaussValues(TTree * Analyze,string histo,double &rms, double &mean, string min,string max)
 {
   string temp = histo;
@@ -79 +79 @@
 
   string nom = histo.erase(0,histo.find(">>")+2);
-  TH1F *h = new TH1F(nom.c_str(),"",100,0,2);
-  Analyze->Draw(temp.c_str(),mintemp.c_str(),maxtemp.c_str());
+  TH1F *h = new TH1F(nom.c_str(),"",200,0,3);
+
+  string all = min + " && " + max;
+  Analyze->Draw(temp.c_str(),all.c_str());
   h->SetMarkerSize(0.6);
+//  double MeanFix = ;
   double MeanFix = h->GetMean();
   double RangMin = h->GetMean()-h->GetRMS();
   double RangMax = h->GetMean()+h->GetRMS();
-  TF1 *Gauss = new TF1("Gauss","gausn",RangMin,RangMax);
+  TF1 *Gauss = new TF1("Gauss","gaus",RangMin,RangMax);
   Gauss->FixParameter(1,MeanFix);
   h->Fit("Gauss","R");
@@ -95 +98 @@
   h->Draw("P");
   h->GetXaxis()->SetTitle("E_{T}^{rec}/E_{T}^{MC} [GeV]");
-}
+Gauss->Delete();
+  
+}
+
+void GaussValuesETmis(TTree * Analyze,string histo,double &rms, double &mean, string min,string max)
+{
+  string temp = histo;
+  string mintemp = min;
+  string maxtemp = max;
+
+  string nom = histo.erase(0,histo.find(">>")+2);
+  TH1F *h = new TH1F(nom.c_str(),"",50,-300,300);
+
+  string all = min + " && " + max;
+  Analyze->Draw(temp.c_str(),all.c_str());
+  h->SetMarkerSize(0.6);
+  double RangMin = h->GetMean()-h->GetRMS();
+  double RangMax = h->GetMean()+h->GetRMS();
+  TF1 *Gauss = new TF1("Gauss","gaus",RangMin,RangMax);
+  h->Fit("Gauss","R");
+  h->Fit("Gauss","RI");
+  h->Fit("Gauss","RI");
+  Double_t* params = Gauss->GetParameters();
+  rms=params[2];
+  mean=params[1];
+  h->Draw("P");
+  h->GetXaxis()->SetTitle("E_{T}^{rec}-E_{T}^{MC} [GeV]");
+}
+
 
 TH1F * HiggsMakeNormTH1F(Float_t Scale, Int_t numBin,Float_t minX,Float_t maxX, TTree * Analyze,string histo, Int_t Lcolor, Int_t Fcolor, Int_t Lstyle,Int_t width=2,bool Log=false)