Changeset 5080 for trunk/MagicSoft/Mars/mhflux

Timestamp:

09/16/04 17:00:50 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MAlphaFitter.h (modified) (3 diffs)
• MHFalseSource.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhflux/MAlphaFitter.h

@@ -42 +42 @@
     MAlphaFitter() : fFunc(new TF1("", "gaus(0) + pol1(3)", 0, 90)), fSigInt(10), fSigMax(75), fBgMin(45), fBgMax(85), fPolynomOrder(1), fCoefficients(3+fPolynomOrder+1)
     {
+        gROOT->GetListOfFunctions()->Remove(fFunc);
     }
 
@@ -67 +68 @@
         TF1 *fcn = f.fFunc;
         f.fFunc = new TF1(*fFunc);
+        gROOT->GetListOfFunctions()->Remove(f.fFunc);
         delete fcn;
     }
@@ -74 +76 @@
     void SetBackgroundFitMin(Float_t s)    { fBgMin        = s; }
     void SetBackgroundFitMax(Float_t s)    { fBgMax        = s; }
-    void SetPolynomOrder(Int_t s)          { fPolynomOrder = s; delete fFunc; fFunc=new TF1 ("", Form("gaus(0) + pol%d(3)", s)); fCoefficients.Set(3+s+1); fCoefficients.Reset(); }
+    void SetPolynomOrder(Int_t s)          { fPolynomOrder = s; delete fFunc; fFunc=new TF1 ("", Form("gaus(0) + pol%d(3)", s));
+        gROOT->GetListOfFunctions()->Remove(fFunc);
+        fCoefficients.Set(3+s+1); fCoefficients.Reset(); }
 
     Double_t GetEventsExcess() const       { return fEventsExcess; }

trunk/MagicSoft/Mars/mhflux/MHFalseSource.cc

@@ -115 +115 @@
 
 #include <TF1.h>
+#include <TF2.h>
 #include <TH2.h>
 #include <TGraph.h>
@@ -729 +730 @@
     gPad->Modified();
     gPad->cd();
+}
+
+Double_t FcnGauss2d(Double_t *x, Double_t *par)
+{
+    TVector2 v = TVector2(x[0], x[1]).Rotate(par[5]*TMath::DegToRad());
+
+    const Double_t g0 = TMath::Gaus(v.X(), par[1], par[2]);
+    const Double_t g1 = TMath::Gaus(v.Y(), par[3], par[4]);
+
+    //Gaus(Double_t x, Double_t mean=0, Double_t sigma=1, Bool_t norm=kFALSE);
+    return par[0]*(g0+g1);
 }
 
@@ -897 +909 @@
         }
 
+    *fLog << "Done." << endl;
+
     h0a.GetXaxis()->SetRangeUser(0, maxalpha0*1.5);
     h0b.GetXaxis()->SetRangeUser(0, maxalpha0*1.5);
@@ -926 +940 @@
     hist->Draw("colz");
     hist->SetBit(kCanDelete);
+
+
+    TF2 f2d("", FcnGauss2d, -1.5, 1.5, -1.5, 1.5, 6);
+    f2d.SetParName(0, "Max   sigma");
+    f2d.SetParName(1, "Mean_1  deg");
+    f2d.SetParName(2, "Sigma_1 deg");
+    f2d.SetParName(3, "Mean_2  deg");
+    f2d.SetParName(4, "Sigma_2 deg");
+    f2d.SetParName(5, "Phi     deg");
+    f2d.SetParLimits(1, -1, 1);    // mu_1
+    f2d.SetParLimits(3, -1, 1);    // mu_2
+    f2d.SetParLimits(2, 0, 1);     // sigma_1
+    f2d.SetParLimits(4, 0, 1);     // sigma_2
+    f2d.SetParLimits(5, 0, 90);    // phi
+    f2d.SetParameter(0, maxs); // A
+    f2d.SetParameter(1, 0);    // mu_1
+    f2d.SetParameter(2, 0.2);  // sigma_1
+    f2d.SetParameter(3, 0);    // mu_2
+    f2d.SetParameter(4, 0.2);  // sigma_2
+    f2d.SetParameter(5, 0);    // phi
+    hist->Fit(&f2d, "IN0");
+    //f2d.DrawClone("cont2same")->SetBit(kCanDelete);
+
+
     catalog->Draw("mirror same *");
     c->cd(3);