[2328] | 1 | // ArchTOIPipe (C) CEA/DAPNIA/SPP IN2P3/LAL
|
---|
| 2 | // Eric Aubourg
|
---|
| 3 | // Christophe Magneville
|
---|
| 4 | // Reza Ansari
|
---|
[2397] | 5 | // $Id: piotoirdr.cc,v 1.9 2003-06-05 02:18:06 aubourg Exp $
|
---|
[2328] | 6 |
|
---|
| 7 | #include "piotoirdr.h"
|
---|
| 8 |
|
---|
| 9 | extern void fits_lock();
|
---|
| 10 | extern void fits_unlock();
|
---|
| 11 |
|
---|
[2396] | 12 | //char** environ; // For darwin dylibs with piolib
|
---|
[2392] | 13 |
|
---|
[2328] | 14 | PIOTOIReader::PIOTOIReader(string grp, string obj, string flg) {
|
---|
| 15 | group = grp;
|
---|
| 16 | object = obj;
|
---|
| 17 | flagdef = flg;
|
---|
| 18 | pioGroup = NULL;
|
---|
| 19 |
|
---|
| 20 | }
|
---|
| 21 |
|
---|
| 22 | PIOTOIReader::~PIOTOIReader() {
|
---|
| 23 | if (pioGroup != NULL) {
|
---|
| 24 | // ??
|
---|
| 25 | }
|
---|
| 26 | }
|
---|
| 27 |
|
---|
| 28 | void PIOTOIReader::init() {
|
---|
| 29 | fits_lock();
|
---|
[2359] | 30 | pioGroup = PIOOpenTOI(const_cast<char*>(group.c_str()), "r");
|
---|
[2328] | 31 | fits_unlock();
|
---|
| 32 | readBounds();
|
---|
| 33 | bufferSize = 1000;
|
---|
| 34 | cout << "PIOTOIReader : opened " << group
|
---|
| 35 | << " " << snBegin << " - " << snEnd << endl;
|
---|
| 36 | if (flagdef != "") {
|
---|
| 37 | cout << " with flagdef = " << flagdef << endl;
|
---|
| 38 | }
|
---|
| 39 |
|
---|
| 40 | declareOutput(object);
|
---|
| 41 | }
|
---|
| 42 |
|
---|
| 43 | void PIOTOIReader::readBounds() {
|
---|
| 44 | PIOSTRING* flgName;
|
---|
| 45 | PIOINT nbFlg;
|
---|
| 46 | PIOSTRING* toiType;
|
---|
| 47 | PIOSTRING* toiName;
|
---|
| 48 | PIOLONG* beginIndex;
|
---|
| 49 | PIOLONG* endIndex;
|
---|
| 50 | PIOINT nbTOI;
|
---|
| 51 | PIOSTRING ROIGroup;
|
---|
| 52 |
|
---|
| 53 | fits_lock();
|
---|
| 54 | PIOInfoTOI(&flgName,
|
---|
| 55 | &nbFlg,
|
---|
| 56 | &toiType,
|
---|
| 57 | &toiName,
|
---|
| 58 | &beginIndex,
|
---|
| 59 | &endIndex,
|
---|
| 60 | &nbTOI,
|
---|
| 61 | ROIGroup,
|
---|
| 62 | pioGroup);
|
---|
| 63 |
|
---|
[2369] | 64 | // On recupere un tableau de beginIndex, il faut retrouver celui qui nous concerne
|
---|
| 65 | for (int i=0; i<nbTOI; i++) {
|
---|
| 66 | cout << toiName[i] << endl;
|
---|
| 67 | if (object == toiName[i]) {
|
---|
| 68 | snBegin = beginIndex[i];
|
---|
| 69 | snEnd = endIndex[i];
|
---|
| 70 | cout << snBegin << " " << snEnd << endl;
|
---|
| 71 | break;
|
---|
| 72 | }
|
---|
| 73 | }
|
---|
[2328] | 74 |
|
---|
| 75 | PIOFreeInfoTOI(flgName,
|
---|
| 76 | toiType,
|
---|
| 77 | toiName,
|
---|
| 78 | beginIndex,
|
---|
| 79 | endIndex);
|
---|
| 80 | fits_unlock();
|
---|
| 81 | }
|
---|
| 82 |
|
---|
| 83 | void PIOTOIReader::run() {
|
---|
| 84 | cout << "PIOTOIReader run " << group << "[" << object << "] "
|
---|
| 85 | << snBegin << " " << snEnd << endl;
|
---|
| 86 |
|
---|
| 87 | for (int snb = snBegin; snb <= snEnd; snb += bufferSize) {
|
---|
| 88 | int sne = snb + bufferSize - 1;
|
---|
| 89 | if (sne > snEnd) sne=snEnd;
|
---|
| 90 | // TBD : check if PIODOUBLE != double;
|
---|
| 91 | char command[80];
|
---|
| 92 | sprintf(command, "Begin=%d; End=%d", snb, sne);
|
---|
[2397] | 93 | double* data = NULL;
|
---|
[2328] | 94 | uint_8* flags = NULL;
|
---|
[2396] | 95 | PIOBYTE* pioflags = NULL;
|
---|
[2328] | 96 | fits_lock();
|
---|
| 97 | PIOLONG n = PIOReadTOI((void**) &data,
|
---|
| 98 | const_cast<char*>(object.c_str()),
|
---|
| 99 | "PIODOUBLE",
|
---|
| 100 | command,
|
---|
| 101 | pioGroup);
|
---|
[2397] | 102 | if ( n<0) {
|
---|
| 103 | cout << " PIOReadTOI : " << n << endl;
|
---|
| 104 | abort();
|
---|
| 105 | }
|
---|
[2328] | 106 | if (flagdef != "") {
|
---|
[2396] | 107 | PIOLONG nf = PIOReadFLGObjectMask(&pioflags,
|
---|
[2382] | 108 | const_cast<char*>(flagdef.c_str()),
|
---|
| 109 | command,
|
---|
| 110 | pioGroup);
|
---|
[2328] | 111 | if (nf != n) {
|
---|
[2382] | 112 | cerr << "*** PIO Error, inconsistent read, data " << n << " flags "
|
---|
| 113 | << nf << endl;
|
---|
[2328] | 114 | abort();
|
---|
| 115 | }
|
---|
[2382] | 116 | flags = new uint_8[nf];
|
---|
[2386] | 117 | for (int i=0; i<nf; i++) flags[i] = (uint_8)pioflags[i];
|
---|
[2382] | 118 | PIODeleteFLG(pioflags, pioGroup);
|
---|
[2328] | 119 | }
|
---|
| 120 | fits_unlock();
|
---|
| 121 |
|
---|
| 122 |
|
---|
| 123 | if (flags != NULL) {
|
---|
| 124 | putData(0, snb, sne-snb+1, data, flags);
|
---|
| 125 | } else {
|
---|
| 126 | putData(0, snb, sne-snb+1, data);
|
---|
| 127 | }
|
---|
| 128 |
|
---|
| 129 | fits_lock();
|
---|
| 130 | PIODeleteTOI(data, pioGroup);
|
---|
[2369] | 131 | // if (flags != NULL) PIODeleteFLG(flags, pioGroup); // PB !
|
---|
[2382] | 132 | if (flags != NULL) delete[] flags;
|
---|
[2328] | 133 | fits_unlock();
|
---|
| 134 | }
|
---|
[2332] | 135 | PIOCloseTOI(pioGroup);
|
---|
[2328] | 136 | }
|
---|