Changeset 9525 for trunk/MagicSoft/Mars/mhflux

Timestamp:

12/10/09 10:08:55 (15 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc

@@ -598 +598 @@
         *fLog << " - Signal Events           " << fEventsSignal << endl;
         *fLog << " - Background Events       " << fEventsBackground << endl;
+        *fLog << " - E/sqrt(B>=Alpha)        " << fEventsExcess/TMath::Sqrt(TMath::Max(fEventsBackground,fScaleFactor)) << endl;
         *fLog << " - Chi^2/ndf (Signal)      " << fChiSqSignal << endl;
         *fLog << " - Chi^2/ndf (Background)  " << fChiSqBg << endl;
         *fLog << " - Signal integrated up to " << fIntegralMax << "°" << endl;
-        *fLog << " - Scale Factor (Off)      " << fScaleFactor << endl;
+        *fLog << " - Off Scale Alpha (Off)   " << fScaleFactor << endl;
     }
 }
@@ -858 +859 @@
         return GetGausSigma();
     case kWeakSource:
-        return GetEventsBackground()<1 ? -GetEventsExcess() : -GetEventsExcess()/TMath::Sqrt(GetEventsBackground());
+        if (GetEventsExcess()<1)
+            return 0;
+        return -GetEventsExcess()/TMath::Sqrt(TMath::Max(GetEventsBackground(), GetEventsBackground()));
     case kWeakSourceLogExcess:
-        return GetEventsBackground()<1 ? -GetEventsExcess() : -GetEventsExcess()/TMath::Sqrt(GetEventsBackground())*TMath::Log10(GetEventsExcess());
+        if (GetEventsExcess()<1)
+            return 0;
+        return -GetEventsExcess()/TMath::Sqrt(TMath::Max(GetEventsBackground(), GetEventsBackground()))*TMath::Log10(GetEventsExcess());
     }
     return 0;