Changeset 13393

Timestamp:

04/20/12 11:04:35 (13 years ago)

Author:

tbretz

Message:

Implemented the possibility to fit the projection per default with a gaussian.

File:

• trunk/Mars/mhist/MHCamEvent.cc (modified) (5 diffs)

trunk/Mars/mhist/MHCamEvent.cc

@@ -59 +59 @@
 #include "MHCamEvent.h"
 
+#include <TF1.h>
+#include <TPaveStats.h>
+
 #include "MLog.h"
 #include "MLogManip.h"
@@ -323 +326 @@
 }
 
-void MHCamEvent::Paint(Option_t *)
+void MHCamEvent::Paint(Option_t *opt)
 {
     TVirtualPad *pad = gPad;
@@ -333 +336 @@
         if (h->GetMaximum()>0)
             gPad->SetLogy();
+
+        if (TString(opt).Contains("gaus", TString::kIgnoreCase))
+        {
+            // Create a copy of the function if it is not yet there
+            if (!h->GetFunction("gaus"))
+                h->Fit("gaus", "q");
+
+            TF1 *f = h->GetFunction("gaus");
+            if (f)
+            {
+                f->SetParameters(h->GetMean(), h->GetRMS(), h->GetMaximum());
+                f->SetLineWidth(1);
+                h->Fit("gaus", "q");
+            }
+
+            TPaveStats *s1 = dynamic_cast<TPaveStats*>(h->GetListOfFunctions()->FindObject("stats"));
+            if (s1)
+            {
+                if (s1->GetOptStat()==1111)
+                    s1->SetOptStat(1000000001);
+                if (s1->GetOptFit()==1)
+                    s1->SetOptFit(1111);
+            }
+
+        }
     }
 
@@ -350 +378 @@
 }
 
-void MHCamEvent::Draw(Option_t *)
+void MHCamEvent::Draw(Option_t *opt)
 {
     TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this);
@@ -356 +384 @@
     pad->SetBorderMode(0);
 
-    AppendPad();
+    AppendPad(opt);
 
     TString name = MString::Format("%s_1", pad->GetName());