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

Timestamp:

01/19/03 14:52:29 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.cc

@@ -158 +158 @@
         return -5.;
 
+    const UInt_t n = geom->GetNumPixels();
+
     Float_t minval = (*this)[0].GetNumPhotons();
 
@@ -164 +166 @@
         const MCerPhotPix &pix = (*this)[i];
 
+        const UInt_t id = pix.GetPixId();
+        if (id>=n)
+            continue;
+
         Float_t testval = pix.GetNumPhotons();
 
         if (geom)
-            testval *= geom->GetPixRatio(pix.GetPixId());
+            testval *= geom->GetPixRatio(id);
 
         if (testval < minval)
@@ -187 +193 @@
         return 50.;
 
+    const UInt_t n = geom->GetNumPixels();
+
     Float_t maxval = (*this)[0].GetNumPhotons();
 
@@ -193 +201 @@
         const MCerPhotPix &pix = (*this)[i];
 
+        const UInt_t id = pix.GetPixId();
+        if (id>=n)
+            continue;
+
         Float_t testval = pix.GetNumPhotons();
-
         if (geom)
-            testval *= geom->GetPixRatio(pix.GetPixId());
+            testval *= geom->GetPixRatio(id);
 
         if (testval > maxval)
@@ -322 +333 @@
     return NULL;
 }
+
+/*
+// --------------------------------------------------------------------------
+//
+// Use this function to sum photons in events together.
+//
+Bool_t MCerPhotEvt::AddEvent(const MCerPhotEvt &evt)
+{
+    if (evt.fNumPixels<=0)
+    {
+        *fLog << "Warning - Event to be added has no pixels." << endl;
+        return kFALSE;
+    }
+    if (fNumPixels<=0)
+    {
+        *fLog << "Warning - Event to add pixels to has no pixels." << endl;
+        return kFALSE;
+    }
+
+    for (UInt_t i=0; i<evt.fNumPixels; i++)
+    {
+        const UInt_t id = evt[i].GetPixId();
+
+        MCerPhotPix *pix2 = GetPixById(id);
+        if (!pix2)
+        {
+            *fLog << "Error - Pixel#" << dec << id << " does not exist in this event!" << endl;
+            return kFALSE;
+        }
+
+        pix2->AddNumPhotons(evt[i].GetNumPhotons());
+    }
+    return kTRUE;
+}
+*/