Changeset 7151 for trunk/MagicSoft/Mars/mhflux/MHAlpha.cc

Timestamp:

06/13/05 10:20:58 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MHAlpha.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhflux/MHAlpha.cc

@@ -292 +292 @@
     if (!fResult)
         return kFALSE;
+    fSigma = (MParameterD*)const_cast<MParList*>(pl)->FindCreateObj("MParameterD", "GaussSigma");
+    if (!fSigma)
+        return kFALSE;
 
     //fExcess = (MParameterD*)const_cast<MParList*>(pl)->FindCreateObj("MParameterD", "MExcess");
@@ -469 +472 @@
             size = fHillas->GetSize();
         energy = fEnergy   ? fEnergy->GetVal() : (fHillas?fHillas->GetSize():1000);
-        theta  = fPointPos ? fPointPos->GetZd()   : 0;
-    }
-
-    //if (size>0)
-    //    alpha /= (2.4 + 1.13*(log10((energy-14)/0.37)-5)*(log10((energy-14)/0.37)-5))/15;
+        theta  = fPointPos ? fPointPos->GetZd() : 0;
+    }
 
     // enhance histogram if necessary
@@ -480 +480 @@
     // Fill histograms
     fHist.Fill(theta, energy, TMath::Abs(alpha), w);
-
-    // Check cuts
-    /*
-    if ( (fEnergyMin>=0 && energy<fEnergyMin) ||
-         (fEnergyMax>=0 && energy>fEnergyMax) ||
-         (fSizeMin>=0   && size  <fSizeMin)   ||
-         (fSizeMax>=0   && size  >fSizeMin) )
-         return kTRUE;
-         */
 
     if (!fSkipHistTime)
@@ -851 +842 @@
 
     fResult->SetVal(fFit.GetMinimizationValue());
+    fSigma->SetVal(fFit.GetGausSigma());
 
     if (!fSkipHistEnergy)