Changeset 5379 for trunk/MagicSoft/Mars/mtemp/mifae/library/MIslandsClean.cc

Timestamp:

11/11/04 09:55:25 (20 years ago)

Author:

aliu

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mifae/library/MIslandsClean.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/library/MIslandsClean.cc

@@ -144 +144 @@
 Int_t MIslandsClean::Process()
 {
-  //
-  //eliminates the island with a signal-to-noise 
-  //lower than a given limit 
-  //
-  //if ( fIslandCleaningMethod == kNoTiming ){
   
   MImgIsland *imgIsl = new MImgIsland;
@@ -156 +151 @@
   Int_t idPix = -1;
 
-  if ( fIslandCleaningMethod == 1 ) {
+  ////////////////////////////////////////////////////
+  //
+  //       TIME SPREAD ISLAND CLEANING
+  //  eliminates the island with a time spread 
+  //         higher than a given limit 
+  /////////////////////////////////////////////////// 
+  if( fIslandCleaningMethod == 0 ){
+    while ((imgIsl=(MImgIsland*)Next())) {
+      
+      //fIslCleanThreshold has different values, FIXME, put two variables
+      if(imgIsl->GetTimeSpread() > fIslCleanThres)
+        {
+          pixNum = imgIsl->GetPixNum();
+          
+          for(Int_t k = 0; k<pixNum; k++)
+            {
+              idPix = imgIsl->GetPixList(k);
+              MCerPhotPix &pix  = (*fEvt)[idPix];
+              pix.SetPixelUnused();
+            }
+        }
+    }
+    
+  }
+  
+
+  ////////////////////////////////////////////////////
+  //
+  //      SIGNAL TO NOISE ISLAND CLEANING
+  //  eliminates the island with a signal-to-noise 
+  //         lower than a given limit 
+  /////////////////////////////////////////////////// 
+  else if ( fIslandCleaningMethod == 1 ) {
     while ((imgIsl=(MImgIsland*)Next())) {
       
@@ -173 +200 @@
   }  
   
-  //
-  //eliminates the island with a time spread  
-  //higher than a given limit 
-  //
-  //else if( fIslandCleaningMethod == kTiming ){
-  else if( fIslandCleaningMethod == 0 ){
-    while ((imgIsl=(MImgIsland*)Next())) {
-
-      //fIslCleanThreshold has different values, FIXME, put two variables
-      if(imgIsl->GetTimeSpread() > fIslCleanThres)
+
+  ////////////////////////////////////////////////////
+  //
+  //      ISLAND SIZE OVER EVENT SIZE ISLAND CLEANING
+  //  eliminates the island with an island size over event size  
+  //         lower than a given limit 
+  /////////////////////////////////////////////////// 
+  else if ( fIslandCleaningMethod == 2 ) {
+    Float_t size = 0;
+    while ((imgIsl=(MImgIsland*)Next())) {
+      size += imgIsl->GetSizeIsl();
+    }
+    
+    Next.Reset();
+    while ((imgIsl=(MImgIsland*)Next())) {
+      
+      if(imgIsl->GetSizeIsl()/size < fIslCleanThres)
         {
           pixNum = imgIsl->GetPixNum();
@@ -193 +227 @@
             }
         }
-    }
-    
-  }
-  
-  //
-  //eliminates all the islands except the continent, 
-  //i.e. the larger island in the event
-  //
+    }   
+  }  
+  
+  
+  ////////////////////////////////////////////////////
+  //
+  //       CONTINENT ISLAND CLEANING
+  //  eliminates all the islands except the continent 
+  //      i.e. the larger island in the event
+  //     according the number of pixels
+  /////////////////////////////////////////////////// 
   else if( fIslandCleaningMethod == 3 ){
     Int_t i = 0;
@@ -219 +256 @@
   }
   
+
+  ////////////////////////////////////////////////////
+  //
+  //       LARGER and SECOND LARGER ISLAND CLEANING
+  // eliminates all the islands except the two biggest 
+  //               ones according size
+  //                
+  /////////////////////////////////////////////////// 
+  else if( fIslandCleaningMethod == 4 ){
+
+    Int_t i = 0;
+    Int_t islnum = fIsl->GetIslNum();
+    Float_t islSize[islnum]; 
+    Int_t islIdx[islnum]; 
+
+    for (Int_t j = 0; j<islnum ; j++){
+      islSize[j] = -1;
+      islIdx[j] = -1;
+    }
+    
+    while ((imgIsl=(MImgIsland*)Next())) {
+      
+      islSize[i] = imgIsl->GetSizeIsl();
+      i++;
+    }
+ 
+    
+    TMath::Sort(islnum, islSize, islIdx, kTRUE);
+    
+    i = 0;
+    Next.Reset();
+    while ((imgIsl=(MImgIsland*)Next())) {
+      if (islnum > 1 && islIdx[0]!=i && islIdx[1]!=i){
+        
+        //cout <<  "removed " << i << " isl 0" << islIdx[0] << " isl 1" << islIdx[1] << endl;
+          
+        pixNum = imgIsl->GetPixNum();
+        
+        for(Int_t k = 0; k<pixNum; k++)
+          {
+            idPix = imgIsl->GetPixList(k);
+            MCerPhotPix &pix  = (*fEvt)[idPix];
+            pix.SetPixelUnused();
+          }
+      }
+      i++;   
+    }  
+  }
+
+
+
+  ///////////////////////////////////////////////////////
+  //
+  //       LARGER and SECOND LARGER ISLAND CLEANING II
+  // eliminates all the islands except the two biggest 
+  // ones according size, if the second one almost has 
+  // the 80% of the size of the biggest one
+  //
+  //                
+  ////////////////////////////////////////////////////// 
+  else if( fIslandCleaningMethod == 5 ){
+
+    Int_t i = 0;
+    Int_t islnum = fIsl->GetIslNum();
+    Float_t islSize[islnum]; 
+    Int_t islIdx[islnum]; 
+
+    for (Int_t j = 0; j<islnum ; j++){
+      islSize[j] = -1;
+      islIdx[j] = -1;
+    }
+    
+    while ((imgIsl=(MImgIsland*)Next())) {
+      
+      islSize[i] = imgIsl->GetSizeIsl();
+      i++;
+    }
+ 
+    
+    TMath::Sort(islnum, islSize, islIdx, kTRUE);
+    
+    i = 0;
+    Next.Reset();
+    while ((imgIsl=(MImgIsland*)Next())) {
+
+      if (islnum > 1 && islIdx[0]!=i && islIdx[1]!=i){
+          
+        pixNum = imgIsl->GetPixNum();
+        
+        for(Int_t k = 0; k<pixNum; k++)
+          {
+            idPix = imgIsl->GetPixList(k);
+            MCerPhotPix &pix  = (*fEvt)[idPix];
+            pix.SetPixelUnused();
+          }
+      }
+      else if(islnum>1 && islSize[islIdx[1]]<0.6*islSize[islIdx[0]]){
+        
+        pixNum = imgIsl->GetPixNum();
+      
+        for(Int_t k = 0; k<pixNum; k++)
+          {
+            idPix = imgIsl->GetPixList(k);
+            MCerPhotPix &pix  = (*fEvt)[idPix];
+            pix.SetPixelUnused();
+          }
+      }
+      i++;   
+    }  
+  }
+  
   fEvt->SetReadyToSave();