Changeset 2361 in Sophya for trunk/ArchTOIPipe/Kernel


Ignore:
Timestamp:
Apr 17, 2003, 3:47:34 PM (22 years ago)
Author:
cecile
Message:

nouveaux tuyaux de ring : theta, phi puis Nsamples mesures

Location:
trunk/ArchTOIPipe/Kernel
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/ArchTOIPipe/Kernel/fitsringrdr.cc

    r2360 r2361  
    22//                               Christophe Magneville
    33//                               Reza Ansari
    4 // $Id: fitsringrdr.cc,v 1.4 2003-04-04 14:17:03 cecile Exp $
     4// $Id: fitsringrdr.cc,v 1.5 2003-04-17 13:47:34 cecile Exp $
    55
    66#include "fitsringrdr.h"
    … …  
    7474void FITSRingReader::init() {
    7575
    76   declareOutput("CircTheta");
    77   declareOutput("CircPhi");
    7876  declareOutput("signal");
    79   openFile(fname);
     77  openFile(fname); 
    8078}
    8179
    … …  
    121119
    122120void FITSRingReader::run() {
     121    cout << "!!!!!!! RING " <<  endl;
    123122
    124123  fits_lock();
    … …  
    142141    fits_read_key(fptr, TDOUBLE, "CIRPHI", &phi, commt, &fstatus);
    143142    fits_read_key(fptr, TDOUBLE, "CIRAPER", &aperture, commt, &fstatus);
    144     fits_read_key(fptr, TLONG, "NSAMPLES", &nSamples, commt, &fstatus);
     143    fits_read_key(fptr, TLONG,   "NSAMPLES", &nSamples, commt, &fstatus);
    145144    fits_unlock();
    146     wontNeedBefore(ring*nSamples);
     145    wontNeedBefore(ring*(nSamples+2));
    147146    cout << "!!!!!!! RING " << ring << endl;
    148147   
     148    double data;
     149    fits_lock();
     150    putData(0, ring*(nSamples+2)+0,theta);
     151    putData(0, ring*(nSamples+2)+1,phi);
     152    fits_unlock();     
    149153    for (int i=0; i<nSamples; i++) {
    150       double data;
    151154      fits_lock();
    152155      fits_read_col_dbl(fptr,1,i+1,1,1,0,&data,&anyNul,&fstatus);
    153       putData(0, ring*nSamples+i,theta);
    154       putData(1, ring*nSamples+i,phi);
    155       putData(2, ring*nSamples+i,data);
     156      putData(0, ring*(nSamples+2)+2+i,data);
    156157      fits_unlock();     
    157158    }

  • trunk/ArchTOIPipe/Kernel/fitsringwtr.cc

    r2355 r2361  
    44//                               Christophe Magneville
    55//                               Reza Ansari
    6 // $Id: fitsringwtr.cc,v 1.6 2003-03-28 16:06:40 cecile Exp $
     6// $Id: fitsringwtr.cc,v 1.7 2003-04-17 13:47:34 cecile Exp $
    77
    88#include "fitsringwtr.h"
    … …  
    4040
    4141void FITSRingWriter::init() {
    42   declareInput("CircTheta");
    43   declareInput("CircPhi");
    4442  declareInput("signal");
    4543}
    … …  
    5957  fits_write_key(fptr, TLONG,   "NCIRCLES", &nRings, "Number of rings", &fstatus);
    6058  fits_unlock();
    61  
     59                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
    6260  char* ttype = "DATA";
    6361  char* tform = "1D";
    … …  
    6765    sprintf(xname, "CRING_%d", ring);
    6866
    69     double theta = getData(0, ring*nSamples);
    70     double phi = getData(1, ring*nSamples);
     67    wontNeedBefore(ring*(nSamples+2));
     68    double theta = getData(0, ring*(nSamples+2)+0);
     69    double phi   = getData(0, ring*(nSamples+2)+1);
    7170
    7271    fits_lock();
    … …  
    7978    fits_write_key(fptr, TLONG, "NSAMPLES", &nSamples, "Number of samples on the circle", &fstatus);
    8079    fits_unlock();
    81     wontNeedBefore(ring*nSamples);
    8280
     81    cout << nSamples << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << endl;
     82   
    8383    for (int i=0; i<nSamples; i++) {
    84       double data = getData(2, ring*nSamples+i);
     84      double data = getData(0, ring*(nSamples+2)+2+i);
     85      if (i > 1000) cout << i << " "  << data << endl;
    8586      fits_lock();
    8687      fits_write_col(fptr, TDOUBLE, 1,i+1,1,1, &data, &fstatus);
Note: See TracChangeset for help on using the changeset viewer.