source: Sophya/trunk/ArchTOIPipe/ProcWSophya/toi2ring.cc@ 2331

// ArchTOIPipe           (C)     CEA/DAPNIA/SPP IN2P3/LAL
//                               Eric Aubourg
//                               Christophe Magneville
//                               Reza Ansari
// $Id: toi2ring.cc,v 1.4 2003-02-24 14:47:10 cecile Exp $

#include "toimanager.h"
#include "pexceptions.h"
#include "ctimer.h"
#include "toi2ring.h"
// La valeur "Pi" doit etre celle de smathconst.h a cause du test sur theta
#include "smathconst.h"

////////////////////////////////////////////////////////////////////////
Toi2RING::Toi2RING(long n)
  :totnscount(0), nRings(0)
{
  setNSamples(n);
}

Toi2RING::~Toi2RING()
{
}


////////////////////////////////////////////////////////////////////////
void Toi2RING::Print(::ostream & os)
{
}

////////////////////////////////////////////////////////////////////////
void Toi2RING::init() {
  cout << "Toi2RING::init" << endl;
  declareInput("Theta");     // input  index 0
  declareInput("Phi");     // input  index 1
  declareInput("Phase");     // input  index 2
  declareInput("Value");     // input  index 3
  declareOutput("CircTheta");   // output index 0
  declareOutput("CircPhi");   // output index 1
  declareOutput("Rvalue");   // output index 2
  declareOutput("WCircTheta");   // output index 3
  declareOutput("WCircPhi");   // output index 4
  declareOutput("WRvalue");   // output index 5
}

////////////////////////////////////////////////////////////////////////
void Toi2RING::run() {

long snb = getMinIn();
long sne = getMaxIn();
if(snb>sne) {
  cout<<"Toi2RING::run() - Bad sample interval"<<snb<<" , "<<sne<<endl;
  throw ParmError("Toi2RING::run() - Bad sample interval");
}

//---------------------------------------------------------
try {


uint_4 mSnRead=0, mSnFilled=0, BadCoorRange=0;
 int_8 Nrings=0,i;
 
 double phi_prec=-1.;
 double* r = new double[Nsamples];
 double* wr = new double[Nsamples];
 for(i=0;i<Nsamples;i++) {
   wr[i]=0;
   r[i]=0;
 }

 
 cout << "loop " << snb << " " << sne << endl;
 
for(int k=snb;k<=sne;k++) {

  uint_8 flg = 0;
  double ctheta = getData(0,k);
  double cphi = getData(1,k);
  double phi = getData(2,k);
  double val = getData(3,k);
  mSnRead++;

  i=(int)(phi*Nsamples/360.);
  r[i]+= val;
  wr[i]++;

  cout << k << " " << phi << " " << phi_prec << " " << Nrings <<endl;
  
  if(phi < phi_prec || k == sne){
    
    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(2,Nrings*Nsamples+i,r[i]);
      putData(3,Nrings*Nsamples+i,ctheta);
      putData(4,Nrings*Nsamples+i,cphi);
      putData(5,Nrings*Nsamples+i,wr[i]);
      wr[i]=0;
      r[i]=0.;
    }
    
    Nrings++;
    return;
    
  }
  

  phi_prec=phi;
  totnscount++;
}

cout<<"RING2Toi::run: Samples Read "<<mSnRead
    <<" Filled "<<mSnFilled
    <<" BadCoorRange="<<BadCoorRange<<endl;

//---------------------------------------------------------
} catch (PException & exc) {
  cout<<"Toi2RING: Catched Exception "<<(string)typeid(exc).name()
      <<"\n .... Msg= "<<exc.Msg()<<endl;
}

return;                                                                            
}