Changeset 2311 in Sophya for trunk/ArchTOIPipe/ProcWSophya/toi2map.cc

Timestamp:

Jan 10, 2003, 9:04:47 AM (23 years ago)

Author:

aubourg

Message:

toi2map.cc

File:

• trunk/ArchTOIPipe/ProcWSophya/toi2map.cc (modified) (4 diffs)

trunk/ArchTOIPipe/ProcWSophya/toi2map.cc

@@ -3 +3 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: toi2map.cc,v 1.31 2002-10-24 10:30:50 aubourg Exp $
+// $Id: toi2map.cc,v 1.32 2003-01-10 08:04:47 aubourg Exp $
 
 #include "machdefs.h"
@@ -158 +158 @@
 cout<<" (getData() not bufferized)"<<endl;
 #endif
-
+ cout << "calib = " << mCalibFactor << endl;
 // Remplissage des spheres
 for(int s=snb;s<=sne;s++) {
@@ -217 +217 @@
   int_4 ipix = mMap->PixIndexSph(theta,phi);
   if ((ipix < 0) || (ipix >= mMap->NbPixels()) ) continue;
+  /*
+    if (mNSnFill % 1000 == 0) {
+    cout << ipix << " " << (*mMap)(ipix) << " + " << bolo*mCalibFactor <<
+    " " << ((*mWMap)(ipix)) << endl;
+      }
+  */
   (*mMap)(ipix) += bolo*mCalibFactor;    
   ((*mWMap)(ipix)) += 1;
@@ -228 +234 @@
  for(int_4 i=0;i<mMap->NbPixels();i++) {
    r_8 wf = (*mWMap)(i);
+   /*
+     if (i%100 == 0) {
+     cout << "pix " << i << " w = " << wf << " map = " <<
+     (*mMap)(i);
+     }
+   */
    if(wf>0.) {mNpixFill++; (*mMap)(i) /= wf;}
+   /*
+     if (i%100 == 0) {
+     cout << " -> " << (*mMap)(i);
+     }
+   */
    int_4 nf = int_4(wf);
    if(nf>=NFILL) nf=NFILL-1; NFill[nf]++;