Changeset 603 for trunk/MagicSoft/Mars/manalysis

Timestamp:

02/21/01 14:41:54 (24 years ago)

Author:

harald

Message:

implemented the Image Cleaning a la CT1 to the class MCerPhotEvt

changed the readCT1.C file to show the effects of the image cleaning
a la CT1

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MCerPhotEvt.cc (modified) (8 diffs)
• MCerPhotEvt.h (modified) (3 diffs)

trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.cc

@@ -18 +18 @@
   fPixId    = pix ; 
   fIsUsed   = kTRUE ; 
+  fIsCore   = kFALSE ; 
   fPhot     = phot ; 
   fErrPhot  = errphot ; 
@@ -26 +27 @@
   fPixId    = pix ; 
   fIsUsed   = kTRUE ; 
+  fIsUsed   = kFALSE ; 
   fPhot     = phot ; 
   fErrPhot  = errphot ; 
@@ -41 +43 @@
       cout << " UnUsed " ;
 
+  if ( fIsCore == kTRUE )
+    cout << "   Core " ; 
+  else 
+    if ( fIsCore == kFALSE )
+      cout << "        " ;
+
   cout << "  Nphot= " << fPhot
        << "  Error(Nphot) = " << fErrPhot
@@ -127 +135 @@
     {
       entry = ((MCerPhotPix *) fPixels->At(il))->GetPhotons() ; 
-      noise =  3 * ((MCerPhotPix *) fPixels->At(il))->GetErrorPhot() ; 
+      noise = ((MCerPhotPix *) fPixels->At(il))->GetErrorPhot() ; 
       
       if (entry < 3 * noise ) 
@@ -154 +162 @@
             id2 = fNN->GetNN(id, in ) ; 
             
-            if ( PixelIsUsed(id2) == kTRUE ) 
-              cout << " hulibu " << id << "/" << id2 << endl ;  
-            
-            
+            if (id2 >=0 ) { 
+              if ( PixelIsUsed(id2) == kTRUE ) 
+                itest++ ; 
+            }   
+          } 
+          
+          // 
+          //   check if no neighbor, then set unused 
+          //
+          if ( itest == 0 ) 
+            { 
+              ((MCerPhotPix *) fPixels->At(il))->SetPixelUnused() ; 
+            }
+
+        } 
+    } 
+  
+  // now we declare all pixels that survive as CorePixels
+  
+  for (Int_t il=0; il<fPixels->GetEntries(); il++ ) 
+    {
+      if ( ((MCerPhotPix *) fPixels->At(il))->IsPixelUsed() == kTRUE )
+        { 
+          ((MCerPhotPix *) fPixels->At(il))->SetCorePixel() ; 
+        }
+    }
+
+} 
+
+void MCerPhotEvt::CleanLevel3()
+{
+  //   Look for the boundary pixels around the core pixels
+  //   if a pixel has more than 2.5 sigma, and a core neigbor
+  //   it is declared as used. 
+  
+  Float_t entry, noise ; 
+  Int_t   id, id2 ; 
+  for (Int_t il=0; il<fPixels->GetEntries(); il++ ) 
+    {
+      if ( ((MCerPhotPix *) fPixels->At(il))->IsCorePixel() == kFALSE )
+        {  
+          entry = ((MCerPhotPix *) fPixels->At(il))->GetPhotons() ; 
+          noise = ((MCerPhotPix *) fPixels->At(il))->GetErrorPhot() ; 
+      
+          if (entry > 2.5 * noise ) { 
+            id = ((MCerPhotPix *) fPixels->At(il))->GetPixId()  ; 
+            for (Int_t in=0 ; in < 6 ; in++ )
+              { 
+                id2 = fNN->GetNN(id, in ) ; 
+                if (id2 >=0 ) 
+                  { 
+                    if ( PixelIsCore(id2) == kTRUE ) { 
+                      ((MCerPhotPix *) fPixels->At(il))->SetPixelUsed() ; 
+                      break ; 
+                    } 
+                  } 
+              } 
           } 
         } 
@@ -163 +224 @@
 } 
 
+
+
+
 Bool_t MCerPhotEvt::PixelExist(Int_t id ) 
 { 
@@ -171 +235 @@
       if ( id == ((MCerPhotPix *) fPixels->At(il))->GetPixId() ) { 
         
-        cout << " PixelExist " << il ; 
+        // cout << " PixelExist " << il ; 
         return kTRUE ; 
       } 
@@ -189 +253 @@
            ((MCerPhotPix *) fPixels->At(il))->IsPixelUsed() == kTRUE ) { 
         
-        cout << " PixelIsUsed  " << il ; 
+        // cout << " PixelIsUsed  " << il ; 
+        return kTRUE ; 
+      } 
+    } 
+
+  return kFALSE ; 
+
+} 
+
+Bool_t MCerPhotEvt::PixelIsCore(Int_t id ) 
+{ 
+  //   Checks if in the pixel list is an entry with pixel id
+  
+  for (Int_t il=0; il<fPixels->GetEntries(); il++ ) 
+    {
+      if ( id == ((MCerPhotPix *) fPixels->At(il))->GetPixId() && 
+           ((MCerPhotPix *) fPixels->At(il))->IsCorePixel() == kTRUE ) { 
+        
         return kTRUE ; 
       } 

trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.h

@@ -18 +18 @@
   Int_t    fPixId     ;  //   the pixel Id
   Bool_t   fIsUsed    ;  //   the pixel is used for calculations --> kTRUE
+  Bool_t   fIsCore    ;  //   the pixel is a Core pixel          --> kTRUE
   Float_t  fPhot      ;  //   The number of Cerenkov photons 
   Float_t  fErrPhot   ;  //   the error of fPhot
@@ -59 +60 @@
     }  
   
+  void SetCorePixel() 
+    { 
+      fIsCore = kTRUE ; 
+    } 
+
+  Bool_t IsCorePixel() 
+    { 
+      return fIsCore ; 
+    }  
   
   
@@ -96 +106 @@
   void CleanLevel1() ; 
   void CleanLevel2() ; 
+  void CleanLevel3() ; 
   
   Bool_t PixelExist( Int_t id ) ; 
   Bool_t PixelIsUsed( Int_t id ) ; 
+  Bool_t PixelIsCore( Int_t id ) ; 
  
   Int_t GetPixelId(Int_t i ) ;