Changeset 4293 for trunk/MagicSoft/Mars

Timestamp:

06/15/04 10:10:09 (20 years ago)

Author:

aliu

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mifae

Files:

• Changelog (modified) (1 diff)
• library/MIslandClean.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/Changelog

@@ -25 +25 @@
       islands
 
-   * library/MIslands.[h,cc],MIslandCalc.cc
+   * library/MIslands.[h,cc],MIslandCalc.[h,cc], MIslandCleaning.[h,cc]
      - add a new island cleaning which consists in removing all the
       islands except the larger one

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

@@ -26 +26 @@
 //   islands are supposed to be eliminated.     
 //
+// Other cleanings that are allowed in this code are:
+// 
+//    - Resting only with the continent, i.e. the larger island
 //    
 //  Example:
@@ -164 +167 @@
     
   //
-  //eliminates the island with a signal-to-noise 
-  //lower than a given limit 
+  //eliminates the island with a time spread  
+  //higher than a given limit 
   //
   //else if( fIslandCleaningMethod == kTiming ){
@@ -188 +191 @@
       }     
   }
+
+  //
+  //eliminates all the islands except the continent, 
+  //i.e. the larger island in the event
+  //
+  else if( fIslandCleaningMethod == 3 ){
+    Int_t nisl = fIsl->GetIslNum();
+
+    Int_t max = -1000;    
+    Int_t idMax;
+
+    for(Int_t isl = 0; isl<nisl ; isl++)
+      {
+        if (fIsl->GetPixNum(isl)>max)
+         {
+           max = fIsl->GetPixNum(isl);
+           idMax = isl;
+         } 
+      }     
+
+    for(Int_t isl = 0; isl<nisl ; isl++)
+      {
+        if (isl != idMax)
+          {
+            for(Int_t idx = 0;idx<577; idx++)
+              {
+                MCerPhotPix &pix  = (*fEvt)[idx];
+                
+                if (fIsl->GetIslId(idx) == isl)
+                  pix.SetPixelUnused();
+              }
+          }
+      }
+  }
   
   fEvt->SetReadyToSave();