Changeset 153 in svn for trunk/interface

Timestamp:

Jan 7, 2009, 12:24:34 PM (16 years ago)

Author:

severine ovyn

Message:

used for publi

File:

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

trunk/interface/FuncDef.h

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