Changeset 19158

Timestamp:

08/13/18 22:53:29 (6 years ago)

Author:

Daniela Dorner

Message:

ignore negative values

Location:

trunk/Mars/mimage

Files:

• MHillas.cc (modified) (2 diffs)
• MHillasExt.cc (modified) (2 diffs)
• MNewImagePar.cc (modified) (2 diffs)

trunk/Mars/mimage/MHillas.cc

@@ -291 +291 @@
             continue;
 
+        const Float_t nphot = pix.GetNumPhotons();   // [#phot]
+        if (nphot<0)
+            continue;
+
         const MGeom &gpix = geom[i];
-
-        const Float_t nphot = pix.GetNumPhotons();
 
         fSize  += nphot;                             // [counter]
@@ -339 +341 @@
             continue;
 
+        const Float_t nphot = pix.GetNumPhotons();   // [#phot]
+        if (nphot<0)
+            continue;
+
         const MGeom &gpix = geom[i];
 
         const Float_t dx = gpix.GetX() - fMeanX;     // [mm]
         const Float_t dy = gpix.GetY() - fMeanY;     // [mm]
-
-        const Float_t nphot = pix.GetNumPhotons();   // [#phot]
 
         corrxx += nphot * dx*dx;                     // [mm^2]

trunk/Mars/mimage/MHillasExt.cc

@@ -191 +191 @@
             continue;
 
+        Double_t nphot = pix.GetNumPhotons();   // [#phot]
+        if (nphot<0)
+            continue;
+
         const MGeom &gpix = geom[i];
 
@@ -196 +200 @@
         const Double_t y = gpix.GetY();
         const Double_t t = pix.GetArrivalTime();
-
-        Double_t nphot = pix.GetNumPhotons();        // [1]
 
         // --- time slope ----

trunk/Mars/mimage/MNewImagePar.cc

@@ -159 +159 @@
             continue;
 
+        // signal in pixel
+        Double_t nphot = pix.GetNumPhotons();
+        if (nphot<0)
+            continue;
+
         // Get geometry of pixel
         const MGeom &gpix = geom[i];
@@ -211 +216 @@
         fNumUsedPixels++;
         fUsedArea += gpix.GetA();
-
-        // signal in pixel
-        Double_t nphot = pix.GetNumPhotons();
 
         //