Changeset 11567 for trunk/Mars/mhist/MHCamera.cc

Timestamp:

07/24/11 18:29:30 (13 years ago)

Author:

tbretz

Message:

Do not propagate non finite numbers.

File:

• trunk/Mars/mhist/MHCamera.cc (modified) (5 diffs)

trunk/Mars/mhist/MHCamera.cc

@@ -502 +502 @@
 
         const Double_t val = TestBit(kProfile) ? fArray[i+1]/fBinEntries[i+1] : fArray[i+1];
-        if (MatchSector(i, sector, aidx) && (ball || IsUsed(i)) && val<minimum)
+        if (MatchSector(i, sector, aidx) && (ball || IsUsed(i)) && val<minimum && TMath::Finite(val))
             minimum = val;
     }
@@ -530 +530 @@
 
         const Double_t val = TestBit(kProfile) ? fArray[i+1]/fBinEntries[i+1] : fArray[i+1];
-        if (MatchSector(i, sector, aidx) && (ball || IsUsed(i)) && val>maximum)
+        if (MatchSector(i, sector, aidx) && (ball || IsUsed(i)) && val>maximum && TMath::Finite(val))
             maximum = val;
     }
@@ -773 +773 @@
     // Fill the projected histogram
     for (Int_t idx=0; idx<fNcells-2; idx++)
-        if (IsUsed(idx) && MatchSector(idx, sector, aidx))
+        if (IsUsed(idx) && MatchSector(idx, sector, aidx) && TMath::Finite(GetBinContent(idx+1)))
             h1->Fill(GetBinContent(idx+1));
 
@@ -864 +864 @@
     // Fill the projected histogram
     for (Int_t idx=0; idx<fNcells-2; idx++)
-        if (IsUsed(idx) && MatchSector(idx, sector, aidx))
+        if (IsUsed(idx) && MatchSector(idx, sector, aidx) && TMath::Finite(GetBinContent(idx+1)))
             h1->Fill(TMath::Hypot((*fGeomCam)[idx].GetX(),(*fGeomCam)[idx].GetY())*m2d,
                      GetBinContent(idx+1));
@@ -943 +943 @@
     for (Int_t idx=0; idx<fNcells-2; idx++)
     {
-        if (IsUsed(idx) && MatchSector(idx, TArrayI(), aidx))
+        if (IsUsed(idx) && MatchSector(idx, TArrayI(), aidx) && TMath::Finite(GetPixContent(idx)))
             h1->Fill(TMath::ATan2((*fGeomCam)[idx].GetY(),(*fGeomCam)[idx].GetX())*TMath::RadToDeg()+180,
                      GetPixContent(idx));