source: Sophya/trunk/ArchTOIPipe/Kernel/fitsringwtr.cc@ 4010

Last change on this file since 4010 was 2361, checked in by cecile, 22 years ago

nouveaux tuyaux de ring : theta, phi puis Nsamples mesures
File size: 6.6 KB
Line 
1// ArchTOIPipe (C) CEA/DAPNIA/SPP IN2P3/LAL
2// Eric Aubourg
3// Cecile Renault
4// Christophe Magneville
5// Reza Ansari
6// $Id: fitsringwtr.cc,v 1.7 2003-04-17 13:47:34 cecile Exp $
7
8#include "fitsringwtr.h"
9#include "toimanager.h"
10
11extern void fits_lock();
12extern void fits_unlock();
13
14FITSRingWriter::FITSRingWriter(string fn) {
15 fname = fn;
16 fstatus = 0;
17 // Open file
18 remove(fname.c_str());
19 fits_lock();
20 fits_create_file(&fptr,fname.c_str(),&fstatus);
21 fits_unlock();
22 name = "ringwriter";
23 cout << "FITSRINGWriter:: FileName=" << name << endl;
24 }
25
26FITSRingWriter::~FITSRingWriter() {
27}
28
29void FITSRingWriter::setNSamples(long n) {
30 nSamples = n;
31}
32
33void FITSRingWriter::setNRings(long n) {
34 nRings = n;
35}
36
37void FITSRingWriter::setAperture(double a) {
38 aperture = a;
39}
40
41void FITSRingWriter::init() {
42 declareInput("signal");
43}
44
45void FITSRingWriter::run() {
46 //int snb = getMinIn(); ON IGNORE ICI
47 //int sne = getMaxIn(); ON IGNORE ICI
48
49 char buf[100];
50
51 fits_lock();
52 fits_create_img(fptr,BYTE_IMG,0,NULL,&fstatus);
53 fits_write_date(fptr, &fstatus);
54 fits_write_comment(fptr,"ArchTOIPipe FITSRingWriter",&fstatus);
55 strcpy(buf, "RINGSET");
56 fits_write_key(fptr, TSTRING, "PDMTYPE", &buf, "Planck Data Model Type", &fstatus);
57 fits_write_key(fptr, TLONG, "NCIRCLES", &nRings, "Number of rings", &fstatus);
58 fits_unlock();
59
60 char* ttype = "DATA";
61 char* tform = "1D";
62 char xname[100];
63
64 for (int ring = 0; ring<nRings; ring++) {
65 sprintf(xname, "CRING_%d", ring);
66
67 wontNeedBefore(ring*(nSamples+2));
68 double theta = getData(0, ring*(nSamples+2)+0);
69 double phi = getData(0, ring*(nSamples+2)+1);
70
71 fits_lock();
72 fits_create_tbl(fptr,BINARY_TBL,nSamples,1,&ttype,&tform,NULL,xname,&fstatus);
73 strcpy(buf, "RINGDATA");
74 fits_write_key(fptr, TSTRING, "PDMTYPE", &buf, "Planck Data Model Type", &fstatus);
75 fits_write_key(fptr, TDOUBLE, "CIRTHETA", &theta, "Theta angle in radians", &fstatus);
76 fits_write_key(fptr, TDOUBLE, "CIRPHI", &phi, "Phi angle in radians", &fstatus);
77 fits_write_key(fptr, TDOUBLE, "CIRAPER", &aperture, "Aperture angle in radians", &fstatus);
78 fits_write_key(fptr, TLONG, "NSAMPLES", &nSamples, "Number of samples on the circle", &fstatus);
79 fits_unlock();
80
81 cout << nSamples << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << endl;
82
83 for (int i=0; i<nSamples; i++) {
84 double data = getData(0, ring*(nSamples+2)+2+i);
85 if (i > 1000) cout << i << " " << data << endl;
86 fits_lock();
87 fits_write_col(fptr, TDOUBLE, 1,i+1,1,1, &data, &fstatus);
88 fits_unlock();
89 }
90 }
91
92 fits_lock();
93 fits_close_file(fptr, &fstatus);
94 fits_report_error(stderr, fstatus);
95 fits_unlock();
96 cout << "fitsringwriter done " << endl;
97}
Note: See TracBrowser for help on using the repository browser.