Changeset 7147 for trunk/MagicSoft/Mars/mhflux

Timestamp:

06/13/05 09:41:57 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MHThetaSq.cc (modified) (4 diffs)
• MHThetaSq.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mhflux/MHThetaSq.cc

@@ -34 +34 @@
 // For more detailes see MHAlpha.
 //
+// Version 2:
+// ---------
+//  + UInt_t fNumBinsSignal;
+//  + UInt_t fNumBinsTotal;
+//
 //////////////////////////////////////////////////////////////////////////////
 #include "MHThetaSq.h"
@@ -57 +62 @@
 //
 MHThetaSq::MHThetaSq(const char *name, const char *title)
-    : MHAlpha(name, title), fThetaSq(0)
+    : MHAlpha(name, title), fThetaSq(0), fNumBinsSignal(3), fNumBinsTotal(45)
 {
     //
@@ -120 +125 @@
     // Calculate bining which fits alpha-cut
     const Double_t intmax = fit->GetSignalIntegralMax();
-    const UInt_t   nbins  = 75;
-    const UInt_t   nsig   =  5;
+    const UInt_t   nbins  = fNumBinsTotal;
+    const UInt_t   nsig   = fNumBinsSignal;
 
     MBinning binsa(nbins, 0, nbins*intmax/nsig);
@@ -201 +206 @@
    //     fMap[4] = fMatrix->AddColumn("MTime.GetAxisTime");
 }
+
+Int_t MHThetaSq::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
+{
+    Int_t rc = MHAlpha::ReadEnv(env, prefix, print);
+    if (rc==kERROR)
+        return kERROR;
+
+    if (IsEnvDefined(env, prefix, "NumBinsSignal", print))
+    {
+        SetNumBinsSignal(GetEnvValue(env, prefix, "NumBinsSignal", (Int_t)fNumBinsSignal));
+        rc = kTRUE;
+    }
+    if (IsEnvDefined(env, prefix, "NumBinsTotal", print))
+    {
+        SetNumBinsTotal(GetEnvValue(env, prefix, "NumBinsTotal", (Int_t)fNumBinsTotal));
+        rc = kTRUE;
+    }
+    return rc;
+}

trunk/MagicSoft/Mars/mhflux/MHThetaSq.h

@@ -13 +13 @@
     MParameterD  *fThetaSq; //!
 
+    UInt_t fNumBinsSignal;
+    UInt_t fNumBinsTotal;
+
     Bool_t      GetParameter(const MParList &pl);
     Double_t    GetVal() const;
@@ -23 +26 @@
     void InitMapping(MHMatrix *mat, Int_t type=0);
 
+    // MParContainer
+    Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
+
 public:
     MHThetaSq(const char *name=NULL, const char *title=NULL);
 
-    ClassDef(MHThetaSq, 1) // Theta-Plot which is fitted online
+    void SetNumBinsSignal(UInt_t n) { fNumBinsSignal=TMath::Max(n, 1U); }
+    void SetNumBinsTotal(UInt_t n)  { fNumBinsTotal =TMath::Max(n, 1U); }
+
+    ClassDef(MHThetaSq, 2) // Theta-Plot which is fitted online
 };