Changeset 39 in svn for trunk/interface/FuncDef.h

Timestamp:

Nov 18, 2008, 10:30:58 AM (16 years ago)

Author:

severine ovyn

Message:

ok for ETmis

File:

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

trunk/interface/FuncDef.h

@@ -79 +79 @@
 
   string nom = histo.erase(0,histo.find(">>")+2);
-  TH1F *h = new TH1F(nom.c_str(),"",200,0,3);
+  TH1F *h = new TH1F(nom.c_str(),"",50,-3,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","gaus",RangMin,RangMax);
+  TF1 *Gauss = new TF1("Gauss","gaus",-3,3);
   Gauss->FixParameter(1,MeanFix);
-  h->Fit("Gauss","R");
-  h->Fit("Gauss","RI");
-  h->Fit("Gauss","RI");
+  h->Fit("Gauss","QR");
+  h->Fit("Gauss","QRI");
+  h->Fit("Gauss","QRI");
   Double_t* params = Gauss->GetParameters();
   rms=params[2];
@@ -102 +99 @@
 }
 
-void GaussValuesETmis(TTree * Analyze,string histo,double &rms, double &mean, string min,string max)
+void GaussValuesETmis(TTree * Analyze,string histo,double &rms, string min,string max)
 {
   string temp = histo;
@@ -109 +106 @@
 
   string nom = histo.erase(0,histo.find(">>")+2);
-  TH1F *h = new TH1F(nom.c_str(),"",50,-300,300);
+  TH1F *h = new TH1F(nom.c_str(),"",20,-100,100);
 
   string all = min + " && " + max;
@@ -117 +114 @@
   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");
+  //TF1 *Gauss = new TF1("Gauss","gaus");
+  h->Fit("Gauss","QR");
+  h->Fit("Gauss","QRI");
+  h->Fit("Gauss","QRI");
   Double_t* params = Gauss->GetParameters();
   rms=params[2];
-  mean=params[1];
+  //rms=h->GetRMS();
+  //mean=params[1];
   h->Draw("P");
   h->GetXaxis()->SetTitle("E_{T}^{rec}-E_{T}^{MC} [GeV]");