Changeset 1965 for trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.cc

Timestamp:

04/19/03 18:39:05 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.cc (modified) (10 diffs)

trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.cc

@@ -27 +27 @@
 
 #include <math.h>
+#include <limits.h>
 #include <fstream.h>
 
@@ -120 +121 @@
     {
         const MCerPhotPix &pix = (*this)[i];
+        if (!pix.IsPixelUsed())
+            continue;
 
         if (id == pix.GetPixId() && pix.IsPixelUsed())
@@ -139 +142 @@
     {
         const MCerPhotPix &pix = (*this)[i];
+        if (!pix.IsPixelUsed())
+            continue;
 
         if (id == pix.GetPixId() && pix.IsPixelCore())
@@ -160 +165 @@
     const UInt_t n = geom->GetNumPixels();
 
-    Float_t minval = (*this)[0].GetNumPhotons();
-
-    for (UInt_t i=1; i<fNumPixels; i++)
-    {
-        const MCerPhotPix &pix = (*this)[i];
+    Float_t minval = FLT_MAX;
+
+    for (UInt_t i=1; i<fNumPixels; i++)
+    {
+        const MCerPhotPix &pix = (*this)[i];
+        if (!pix.IsPixelUsed())
+            continue;
 
         const UInt_t id = pix.GetPixId();
@@ -195 +202 @@
     const UInt_t n = geom->GetNumPixels();
 
-    Float_t maxval = (*this)[0].GetNumPhotons();
-
-    for (UInt_t i=1; i<fNumPixels; i++)
-    {
-        const MCerPhotPix &pix = (*this)[i];
+    Float_t maxval = -FLT_MAX;
+
+    for (UInt_t i=1; i<fNumPixels; i++)
+    {
+        const MCerPhotPix &pix = (*this)[i];
+        if (!pix.IsPixelUsed())
+            continue;
 
         const UInt_t id = pix.GetPixId();
@@ -229 +238 @@
     {
         const MCerPhotPix &pix = (*this)[i];
+        if (!pix.IsPixelUsed())
+            continue;
 
         Float_t testval = pix.GetNumPhotons()/pix.GetErrorPhot();
@@ -247 +258 @@
         return -5.;
 
-    Float_t maxval = (*this)[0].GetNumPhotons()/(*this)[0].GetErrorPhot();
-
-    for (UInt_t i=1; i<fNumPixels; i++)
-    {
-        const MCerPhotPix &pix = (*this)[i];
+    Float_t maxval = -FLT_MAX;
+
+    for (UInt_t i=1; i<fNumPixels; i++)
+    {
+        const MCerPhotPix &pix = (*this)[i];
+        if (!pix.IsPixelUsed())
+            continue;
 
         Float_t testval = pix.GetNumPhotons()/pix.GetErrorPhot();
@@ -272 +285 @@
         return 50.;
 
-    Float_t minval = (*this)[0].GetErrorPhot();
-
-    for (UInt_t i=1; i<fNumPixels; i++)
-    {
-        const MCerPhotPix &pix = (*this)[i];
+    Float_t minval = FLT_MAX;
+
+    for (UInt_t i=1; i<fNumPixels; i++)
+    {
+        const MCerPhotPix &pix = (*this)[i];
+        if (!pix.IsPixelUsed())
+            continue;
 
         Float_t testval = pix.GetErrorPhot();
@@ -300 +315 @@
         return 50.;
 
-    Float_t maxval = (*this)[0].GetErrorPhot();
-
-    for (UInt_t i=1; i<fNumPixels; i++)
-    {
-        const MCerPhotPix &pix = (*this)[i];
+    Float_t maxval = -FLT_MAX;
+
+    for (UInt_t i=1; i<fNumPixels; i++)
+    {
+        const MCerPhotPix &pix = (*this)[i];
+        if (!pix.IsPixelUsed())
+            continue;
 
         Float_t testval = pix.GetErrorPhot();
@@ -368 +385 @@
 }
 */
+
+void MCerPhotEvt::Scale(Double_t f)
+{
+    fPixels->ForEach(MCerPhotPix, Scale)(f);
+}