Ignore:
Timestamp:
Feb 24, 2003, 3:14:51 PM (23 years ago)
Author:
cecile
Message:

* empty log message *

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/ArchTOIPipe/ProcWSophya/toi2ring.cc

    r2319 r2329  
    33//                               Christophe Magneville
    44//                               Reza Ansari
    5 // $Id: toi2ring.cc,v 1.2 2003-01-30 17:19:25 cecile Exp $
     5// $Id: toi2ring.cc,v 1.3 2003-02-24 14:14:51 cecile Exp $
    66
    77#include "toimanager.h"
     
    6868
    6969 
     70 cout << "loop " << snb << " " << sne << endl;
    7071 
    7172for(int k=snb;k<=sne;k++) {
     
    8283  wr[i]++;
    8384
    84   cout << phi << " " << phi_prec << " " << Nrings <<endl;
     85  cout << k << " " << phi << " " << phi_prec << " " << Nrings <<endl;
    8586 
    86   if(phi < phi_prec){
     87  if(phi < phi_prec || k == sne){
    8788   
    88     putData(0,Nrings,ctheta);
    89     putData(1,Nrings,cphi);
    90 
    9189    for(i=0;i<Nsamples;i++) {
    9290      if(wr[i] == 0) r[i] = -999999.;
    93       cout << i << " " << Nrings << " " << Nsamples << " " << ctheta << " " << cphi << " " << r[i]<< endl;
     91      else r[i] /= wr[i];
     92      cout << "PUT " << i << " " << Nrings << " " << Nsamples << " " << ctheta << " " << cphi << " " << r[i]<< endl;
     93      putData(0,Nrings*Nsamples+i,ctheta);
     94      putData(1,Nrings*Nsamples+i,cphi);
    9495      putData(2,Nrings*Nsamples+i,r[i]);
    95       //      putData(1,Nrings*Nsamples+i,wr[i]);
    9696      wr[i]=0;
    9797      r[i]=0.;
     
    9999   
    100100    Nrings++;
     101    return;
     102   
    101103  }
    102104 
Note: See TracChangeset for help on using the changeset viewer.