Changeset 2355 in Sophya for trunk/ArchTOIPipe/ProcWSophya/toi2ring.cc

Timestamp:

Mar 28, 2003, 5:06:40 PM (23 years ago)

Author:

cecile

Message:

ok

File:

• trunk/ArchTOIPipe/ProcWSophya/toi2ring.cc (modified) (5 diffs)

trunk/ArchTOIPipe/ProcWSophya/toi2ring.cc

@@ -3 +3 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: toi2ring.cc,v 1.4 2003-02-24 14:47:10 cecile Exp $
+// $Id: toi2ring.cc,v 1.5 2003-03-28 16:06:17 cecile Exp $
 
 #include "toimanager.h"
@@ -61 +61 @@
  int_8 Nrings=0,i;
  
- double phi_prec=-1.;
+ double phi_prec=-1.,ctheta_prec,cphi_prec;
+ double phi,val,ctheta,cphi;
+
  double* r = new double[Nsamples];
  double* wr = new double[Nsamples];
@@ -69 +71 @@
  }
 
- 
  cout << "loop " << snb << " " << sne << endl;
  
@@ -75 +76 @@
 
   uint_8 flg = 0;
-  double ctheta = getData(0,k);
-  double cphi = getData(1,k);
-  double phi = getData(2,k);
-  double val = getData(3,k);
+  ctheta = getData(0,k);
+  cphi = getData(1,k);
+  phi = getData(2,k);
+  val = getData(3,k);
   mSnRead++;
 
@@ -85 +86 @@
   wr[i]++;
 
-  cout << k << " " << phi << " " << phi_prec << " " << Nrings <<endl;
   
   if(phi < phi_prec || k == sne){
+    
+    //  if (k==sne) { wfits.setNRings(Nrings); wwfits.setNRings(Nrings);}
     
     for(i=0;i<Nsamples;i++) {
       if(wr[i] == 0) r[i] = -999999.;
       else r[i] /= wr[i];
-      cout << "PUT " << i << " " << Nrings << " " << Nsamples << " " << ctheta << " " << cphi << " " << r[i]<< endl;
-      putData(0,Nrings*Nsamples+i,ctheta);
-      putData(1,Nrings*Nsamples+i,cphi);
+      
+      putData(0,Nrings*Nsamples+i,ctheta_prec);
+      putData(1,Nrings*Nsamples+i,cphi_prec);
       putData(2,Nrings*Nsamples+i,r[i]);
-      putData(3,Nrings*Nsamples+i,ctheta);
-      putData(4,Nrings*Nsamples+i,cphi);
+      putData(3,Nrings*Nsamples+i,ctheta_prec);
+      putData(4,Nrings*Nsamples+i,cphi_prec);
       putData(5,Nrings*Nsamples+i,wr[i]);
       wr[i]=0;
       r[i]=0.;
+    
     }
     
-    Nrings++;
-    return;
-    
+    Nrings++;    
+    cout <<" Ring "<< Nrings << " written " << endl;
   }
-  
-
   phi_prec=phi;
+  ctheta_prec=ctheta;
+  cphi_prec=cphi;
   totnscount++;
-}
+} 
 
 cout<<"RING2Toi::run: Samples Read "<<mSnRead