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

Timestamp:

Jul 14, 2009, 3:23:30 PM (15 years ago)

Author:

Xavier Rouby

Message:

updates needed for atlas jet/met resolution

File:

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

trunk/interface/FuncDef.h

@@ -83 +83 @@
   Analyze->Draw(temp.c_str(),all.c_str());
   h->SetMarkerSize(0.6);
-  double MeanFix = h->GetMean();
+  //double MeanFix = h->GetMean();
   double RangMin = (h->GetMean()-1.5*h->GetRMS());
   double RangMax = (h->GetMean()+1.5*h->GetRMS());
@@ -101 +101 @@
 
 
-void GaussValuesAsymmetry(TTree * Analyze,string histo, 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(),"",500,-10,40);
-  string all = min + " && " + max;
+void GaussValuesAsymmetry(TTree * Analyze,string histo, double &mean, string min,string max,string cut="1==1")
+{
+  string temp = histo;
+  string mintemp = min;
+  string maxtemp = max;
+
+  string nom = histo.erase(0,histo.find(">>")+2);
+  TH1F *h = new TH1F(nom.c_str(),nom.c_str(),50,-3,3);
+  //string all = min + " && " + max + " && abs(JetPTResol.Eta) < 0.5";
+  string all = min + " && " + max + " && " + cut;
+  Analyze->Draw(temp.c_str(),all.c_str());
+  h->SetMarkerSize(0.6);
+//  mean = h->GetMean();
+    TF1 *Gauss = new TF1("Gauss","gaus",-0.4,0.4);
+    //Gauss->FixParameter(1,MeanFix);
+  h->Fit("Gauss","QR");
+  h->Fit("Gauss","QRI");
+  h->Fit("Gauss","QRI");
+  Double_t* params = Gauss->GetParameters();
+  mean=params[1];
+  h->Draw("P");
+  h->GetXaxis()->SetTitle("E_{T}^{rec}/E_{T}^{MC} [GeV]");
+  Gauss->Delete();
+//  double RangMin = (h->GetMean()-2*h->GetRMS());
+//  double RangMax = (h->GetMean()+2*h->GetRMS());
+}
+
+void GaussValuesAsymmetry2(TTree * Analyze,string histo, double &mean, string min,string max,string cut="1==1")
+{
+  string temp = histo;
+  string mintemp = min;
+  string maxtemp = max;
+
+  string nom = histo.erase(0,histo.find(">>")+2);
+  TH1F *h = new TH1F(nom.c_str(),nom.c_str(),50,-3,3);
+  string all = min + " && " + max + " && " + cut;
   Analyze->Draw(temp.c_str(),all.c_str());
   h->SetMarkerSize(0.6);
   mean = h->GetMean();
-  double RangMin = (h->GetMean()-2*h->GetRMS());
-  double RangMax = (h->GetMean()+2*h->GetRMS());
-}
-
+}
 
 
@@ -160 +184 @@
   Analyze->Draw(temp.c_str(),all.c_str());
   h->SetMarkerSize(0.6);
-  double RangMin = (h->GetMean()-1.9*h->GetRMS());
-  double RangMax = (h->GetMean()+1.9*h->GetRMS());
+  double RangMin = (h->GetMean()-2.5*h->GetRMS());
+  double RangMax = (h->GetMean()+2.5*h->GetRMS());
   //TF1 *Gauss = new TF1("Gauss","gaus");
   TF1 *Gauss = new TF1("Gauss","gaus",RangMin,RangMax);