Changeset 19917 for trunk/Mars/mbase

Timestamp:

12/17/19 17:26:49 (5 years ago)

Author:

tbretz

Message:

Added a possibility to produce a pre-defined wavelength distirbution.

Location:

trunk/Mars/mbase

Files:

• MSpline3.cc (modified) (1 diff)
• MSpline3.h (modified) (1 diff)

trunk/Mars/mbase/MSpline3.cc

@@ -234 +234 @@
     }
 
-    return TMath::TwoPi()*MSpline3(x.GetArray(), y.GetArray(), n).Integral();
-}
+    return MSpline3(x.GetArray(), y.GetArray(), n).Integral();
+}
+
+MSpline3 MSpline3::GetIntegralSpline() const
+{
+    const Int_t n = GetNp();
+
+    MArrayD x(n+2);
+    MArrayD y(n+2);
+
+    GetKnot(0, x[0], y[0]);
+
+    y[0] = 0;
+
+    GetKnot(0, x[1], y[1]);
+
+    y[1] = 0;
+
+    for (int i=0; i<n-1; i++)
+    {
+        GetKnot(i+1, x[i+2], y[i+2]);
+        y[i+2] = y[i+1]+IntegralBin(i, x[i+2]);
+    }
+
+    for (int i=2; i<=n; i++)
+        y[i] /= y[n];
+
+    y[n+1] = 1;
+
+    x[0]   = x[1]  -(x[2]-x[1]);
+    x[n+1] = x[n]  +(x[n]-x[n-1]);
+
+    return MSpline3(x.GetArray(), y.GetArray(), n+2, "b1 e1");
+}
+
 
 

trunk/Mars/mbase/MSpline3.h

@@ -64 +64 @@
     Double_t IntegralSolidAngle() const;
 
+    MSpline3 GetIntegralSpline() const;
+
     void Scale(double=1);