Changeset 4218 for trunk/MagicSoft/Mars/mtemp/mifae/library/MIslandCalc.cc

Timestamp:

05/27/04 16:38:48 (20 years ago)

Author:

aliu

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mifae/library/MIslandCalc.cc (modified) (16 diffs)

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

@@ -3 +3 @@
 !
 !   Author(s): Ester Aliu, 2/2004 <aliu@ifae.es>
-!  
+!   
+!   Last Update: 5/2004
 !
 \* ======================================================================== */
@@ -18 +19 @@
 //   - fPixNum[20];        //  number of pixels in the island
 //   - fSigToNoise[20];    //  signal to noise of the island
-//   - fTimeHi[20][577]    //  mean of the arrival time in the hi gain  
-//   - fTimeLo[20][577]    //  mean of the arrival time in the lo gain
-//   - fTimeSpreadHi[20];  //  mean arrival time spread of the island  
-//   - fTimeSpreadLo[20];  //  mean arrival time spread of the island
+//   - fTime[20][577]      //  mean of the arrival time  
+//   - fTimeSpread[20];    //  mean arrival time spread of the island
 //
 //
@@ -30 +29 @@
 //   MArrivalTimeCam
 //
-//  Output Containers:
+// Output Containers:
 //   MIslands
 //
@@ -109 +108 @@
       {
         fIsl = (MIslands*)pList->FindCreateObj("MIslands", AddSerialNumber(fIslName));
-        cout << "kk1" << endl;
+        //cout << "kk1" << endl;
       }
     else
       {
       fIsl = (MIslands*)pList->FindCreateObj(AddSerialNumber("MIslands"));
-      cout << "kk2" << endl;
+      //cout << "kk2" << endl;
       }
     if (!fIsl)
@@ -130 +129 @@
 
   Int_t    npix = 577;
+
   Int_t    sflag;
   Int_t    control;
@@ -148 +148 @@
   
 
+  MCerPhotPix *pix;
+
   //loop over all pixels
-  for(Int_t idx=0 ; idx<npix ; idx++)
-    { 
+  MCerPhotEvtIter Next(fEvt, kFALSE);
+  
+  while ((pix=static_cast<MCerPhotPix*>(Next())))
+    {
+      const Int_t idx = pix->GetPixId();
+
       const MGeomPix &gpix  = (*fCam)[idx];
       const Int_t    nnmax  = gpix.GetNumNeighbors();
-      
-      if(idx<0 || idx==npix)
-        {
-          cout << "Pixel (software) index out of range [0-576]. Available number of pixels=" << npix << endl;
-          return 1;
-        }
-      
-      
+
       if( fEvt->IsPixelUsed(idx)) 
         {
@@ -193 +192 @@
   fIslNum = nvect;
   
-
+  
   // Repeated Chain Corrections
-
+  
   for(Int_t i = 1; i <= nvect; i++)
     {
@@ -203 +202 @@
           for(Int_t k = 0; k < npix; k++)
             {
-        
+              
               if (vect[i][k] == 1 && vect[j][k] == 1)
                 {
@@ -225 +224 @@
     }
   
- 
+  
   Int_t l = 1;
-
+  
   for(Int_t i = 1;  i<= nvect ; i++)
     {
@@ -239 +238 @@
         }
     }
-   
+  
   
   //set the number of islands in one event
@@ -249 +248 @@
   Float_t time[577];
   Float_t timeVariance[fIslNum];
-  //Float_t timeHi[577], timeLo[577];  
-  //Float_t timeVarianceHi[fIslNum];
-  //Float_t timeVarianceLo[fIslNum];
- 
+  
   //reset the "sets" functions
   if (fIslNum <1)
@@ -262 +258 @@
       fIsl->SetSigToNoise(i,-1);
       fIsl->SetTimeSpread(i,-1);
-      //      fIsl->SetTimeSpreadHi(i,-1);
-      // fIsl->SetTimeSpreadLo(i,-1);
-
+    
       for(Int_t idx = 0; idx<npix; idx++)
         {
           fIsl->SetIslId(idx, -1);
           fIsl->SetArrivalTime(i, idx, -1 );
-          // fIsl->SetArrivalTimeHiGain(i, idx, -1 );
-          // fIsl->SetArrivalTimeLoGain(i, idx, -1);
         }
     }
  
-    
-  for(Int_t i = 1; i<=fIslNum ; i++)
+   
+   for(Int_t i = 1; i<=fIslNum ; i++)
     {
       Int_t n = 0;
       Int_t ncore = 0;
-
-      Float_t MINhi = 10000;
-      Float_t MINlo = 10000;
-      Float_t MAXhi = 0;
-      Float_t MAXlo = 0;
+      
       Float_t MIN = 10000;
       Float_t MAX = 0;
-
-      //cout << "Isl #" << i << endl;
 
       signal = 0;
@@ -293 +279 @@
       fPixNum[i-1] = 0;
       timeVariance[i-1] = 0;
-      // timeVarianceHi[i-1] = 0;
-      // timeVarianceLo[i-1] = 0;
-
+      
       for(Int_t idx=0 ; idx<npix ; idx++)
         {
-          const MCerPhotPix &pix = (*fEvt)[idx];
+          
+          MCerPhotPix *pix = fEvt->GetPixById(idx);
           const MPedestalPix &ped  = (*fPed)[idx];
           const MArrivalTimePix &timepix = (*fTime)[idx];
-          
+
+          if (pix == NULL) break;
+            
           if (vect[i][idx]==1){
-    
+            
             fPixNum[i-1]++;
-            signal += pix.GetNumPhotons() * (fCam->GetPixRatio(idx));
+            signal += pix->GetNumPhotons() * (fCam->GetPixRatio(idx));
             noise += pow(ped.GetPedestalRms(),2);
-
+            
             time[n] = timepix.IsLoGainUsed() ? timepix.GetArrivalTimeLoGain() : timepix.GetArrivalTimeHiGain();
-
-            // timeHi[n] = timepix.GetArrivalTimeHiGain(); 
-            // cout << "timeHi" << timeHi[n] << endl;       
-
-            //timeLo[n] = timepix.GetArrivalTimeLoGain(); 
-            //cout << "timeHi" << timeLo[n] << endl;        
-        
+            
+            
             if (fEvt->IsPixelCore(idx)){ 
               
-              /*  if (timeHi[n] < MINhi)
-                MINhi = timeHi[n];
-              if (timeLo[n] < MINlo)
-                MINlo = timeLo[n];
-
-              if (timeHi[n] > MAXhi)
-                MAXhi = timeHi[n];
-              if (timeLo[n] > MAXlo)
-                MAXlo = timeLo[n];
-
-              */
-
               if (time[n] > MAX)
                 MAX = time[n];
@@ -335 +305 @@
                 MIN = time[n];
               
-              //control2[n] = 1;   
-              //timeVarianceHi[i-1] += timeHi[n];
-              //timeVarianceLo[i-1] += timeLo[n];
               ncore++;
             }
-
+            
             fIsl->SetIslId(idx, i-1);
             fIsl->SetArrivalTime(i-1, idx, time[n]);
-            // fIsl->SetArrivalTimeHiGain(i-1, idx, timeHi[n]);
-            // fIsl->SetArrivalTimeLoGain(i-1, idx, timeLo[n]);
-        
+                
             n++;
           }
@@ -351 +316 @@
         }
 
-      Float_t mean = timeVariance[i-1]/ncore;
-      // Float_t meanHi = timeVarianceHi[i-1]/ncore;
-      // Float_t meanLo = timeVarianceLo[i-1]/ncore;
-      
+      // Float_t mean = timeVariance[i-1]/ncore;
+           
       timeVariance[i-1] = 0;
-      //timeVarianceHi[i-1] = 0;
-      //timeVarianceLo[i-1] = 0;
-
-      /*
-      //old
-        for (Int_t k = 0; k <n ; k++)
-        { 
-        if (control2[k] == 1){
-        
-        timeVarianceHi[i-1] += pow(timeHi[k]- meanHi,2); 
-        timeVarianceLo[i-1] += pow(timeLo[k]- meanLo,2); 
-        }
-        }
-        timeVarianceHi[i-1] = sqrt(timeVarianceHi[i-1]/(ncore-1)); 
-        timeVarianceLo[i-1] = sqrt(timeVarianceLo[i-1]/(ncore-1)); */
-      
-      //timeVarianceHi[i-1] = (MAXhi - MINhi)/ncore; 
-      //timeVarianceLo[i-1] = (MAXlo - MINlo)/ncore; 
-
+     
       timeVariance[i-1] = (MAX - MIN)/ncore; 
       timeVariance[i-1] = (MAX - MIN)/ncore; 
@@ -383 +328 @@
       fIsl->SetSigToNoise(i-1,fSigToNoise[i-1]);
       fIsl->SetTimeSpread(i-1, timeVariance[i-1]);
-      //fIsl->SetTimeSpreadHi(i-1, timeVarianceHi[i-1]);
-      //fIsl->SetTimeSpreadLo(i-1, timeVarianceLo[i-1]);
-
-      //cout << " TimeHi Spread: " << timeVarianceHi[i-1] << endl;
-      //cout << " # core pixels: " << ncore << endl;
-      //cout << " # used pixels: " << n << endl;
-      //cout << " SigToNoise: " << fSigToNoise[i-1] << endl << endl;
       
     }