Changeset 2456 in Sophya for trunk/SophyaExt/FitsIOServer

Timestamp:

Nov 18, 2003, 1:06:24 AM (22 years ago)

Author:

cmv

Message:

amelioration architecture classe cmv 18/11/03

Location:

trunk/SophyaExt/FitsIOServer

Files:

• fabtcolread.cc (modified) (28 diffs)
• fabtcolread.h (modified) (8 diffs)

trunk/SophyaExt/FitsIOServer/fabtcolread.cc

@@ -43 +43 @@
 {
   FitsFN = "";
-  NHdu = -1;
+  NHdu = IHdu = HduType = 0;
+  HasBeenPos = false;
   FitsPtr = NULL;
 }
@@ -52 +53 @@
 FitsOpenFile::FitsOpenFile(FitsOpenFile& fof)
 {
-  Init(fof.GetFileName().c_str());
+  Init(fof.FileName().c_str());
 }
 
@@ -62 +63 @@
   Delete();
   FitsFN = "";
-  NHdu = 0;
+  NHdu = IHdu = HduType = 0;
+  HasBeenPos = false;
 }
 
@@ -80 +82 @@
 void FitsOpenFile::Init(const char* fname)
 {
-   // Parametres Generaux
-   FitsFN = fname;
-   NHdu = 0;
-   FitsPtr = NULL;
-
-   // Ouverture du fichier
-   if(FitsFN.size() <= 0 )
-      throw ParmError("FitsOpenFile::Init: Fits file name error\n");
-
-   int sta = 0;
-   if(fits_open_file(&FitsPtr,FitsFN.c_str(),READONLY,&sta)) {
-     printerror(sta);
-     FitsPtr = NULL;
-     throw NullPtrError("FitsOpenFile::Init: Error opening Fits file\n");
-   }
-
-   // Get number of hdu
-   if(fits_get_num_hdus(FitsPtr,&NHdu,&sta)) {
-     printerror(sta);
-     NHdu = 0;
-     Delete();
-     throw NotAvailableOperation("FitsOpenFile::Init: Error getting NHdu\n");
-   }
-   if(NHdu<=0) {
-     Delete();
-     throw SzMismatchError("FitsOpenFile::Init: Bad NHdu\n");
-   }
-
- }
+  // Parametres Generaux
+  FitsFN = fname;
+  NHdu = IHdu = HduType = 0;
+  HasBeenPos = false;
+  FitsPtr = NULL;
+
+  // Ouverture du fichier
+  if(FitsFN.size() <= 0 )
+     throw ParmError("FitsOpenFile::Init: Fits file name error\n");
+
+  int sta = 0;
+  if(fits_open_file(&FitsPtr,FitsFN.c_str(),READONLY,&sta)) {
+    printerror(sta);
+    FitsPtr = NULL;
+    throw NullPtrError("FitsOpenFile::Init: Error opening Fits file\n");
+  }
+
+  // Get number of hdu
+  if(fits_get_num_hdus(FitsPtr,&NHdu,&sta)) {
+    printerror(sta);
+    NHdu = 0;
+    Delete();
+    throw NotAvailableOperation("FitsOpenFile::Init: Error getting NHdu\n");
+  }
+  if(NHdu<=0) {
+    Delete();
+    throw SzMismatchError("FitsOpenFile::Init: Bad NHdu\n");
+  }
+
+  MoveToHDU(1);
+}
+
+/*! Move to an HDU
+\param ihdu: hdu number to move
+\warning ihdu = [1,nhdu]
+\return 0 if positionning failed, ihdu if success
+*/
+int FitsOpenFile::MoveToHDU(int ihdu)
+{
+ if(FitsPtr==NULL)
+   throw NullPtrError("FitsOpenFile::MoveToHDU: no fits file open FitsPtr==NULL\n");
+ int ihdusave = IHdu;
+ if(ihdu<=0) ihdu=1; if(ihdu>NHdu) ihdu=NHdu;
+ int sta=0;
+ if(fits_movabs_hdu(FitsPtr,ihdu,&HduType,&sta)) {
+   printerror(sta);
+   // On se repositionne ou on etait
+   fits_movabs_hdu(FitsPtr,ihdusave,&HduType,&sta);
+   IHdu = ihdusave;
+ } else IHdu = ihdu;
+ return IHdu;
+}
+
+/*! Move to the first HDU of a certain type
+\param hdutype: type of the hdu
+\param hdudeb: start at that hdu
+\return the type of HDU the file is positionned
+*/
+int FitsOpenFile::MoveToFirst(int hdutype,int ihdudeb)
+{
+ if(ihdudeb<=0) ihdudeb=1; if(ihdudeb>NHdu) ihdudeb=NHdu;
+ int ihdusave = IHdu;
+ for(int ihdu=ihdudeb;ihdu<=NHdu;ihdu++) {
+   MoveToHDU(ihdu);
+   if(HduType==hdutype) break;
+ }
+ // Si echec, on se repositionne ou on etait
+ if(HduType!=hdutype) MoveToHDU(ihdusave);
+ return HduType;
+}
+
+/*! Move to the last HDU of a certain type
+\param hdutype: type of the hdu
+\param hdudeb: stop at that hdu
+\return the type of HDU the file is positionned
+*/
+int FitsOpenFile::MoveToLast(int hdutype,int ihdudeb)
+{
+ if(ihdudeb<=0) ihdudeb=1; if(ihdudeb>NHdu) ihdudeb=NHdu;
+ int ihdusave = IHdu;
+ for(int ihdu=NHdu;ihdu>=ihdudeb;ihdu--) {
+   MoveToHDU(ihdu);
+   if(HduType==hdutype) break;
+ }
+ // Si echec, on se repositionne ou on etait
+ if(HduType!=hdutype) MoveToHDU(ihdusave);
+ return HduType;
+}
+
+/*! Print */
+void FitsOpenFile::Print(void)
+{
+ cout<<"FitsOpenFile::Print: "<<FitsFN
+     <<" hdu="<<IHdu<<"/"<<NHdu<<" type="<<HduType
+     <<" hasbeenpos="<<HasBeenPos<<endl;
+}
 
 //////////////////////////////////////////////////////////////
@@ -121 +190 @@
 double FitsOpenFile::ReadKey(fitsfile *fitsptr,char *keyname)
 {
- if(keyname==NULL) return 0.;
+ if(keyname==NULL || fitsptr==NULL) return 0.;
  int sta=0; double val=0.;
  if(fits_read_key(fitsptr,TDOUBLE,keyname,&val,NULL,&sta))
@@ -136 +205 @@
 long FitsOpenFile::ReadKeyL(fitsfile *fitsptr,char *keyname)
 {
- if(keyname==NULL) return 0;
+ if(keyname==NULL || fitsptr==NULL) return 0;
  int sta=0; long val=0;
  if(fits_read_key(fitsptr,TLONG,keyname,&val,NULL,&sta))
@@ -151 +220 @@
 string FitsOpenFile::ReadKeyS(fitsfile *fitsptr,char *keyname)
 {
- if(keyname==NULL) return (string)"";
+ if(keyname==NULL || fitsptr==NULL) return (string)"";
  int sta=0; char val[FLEN_VALUE];
  if(fits_read_key(fitsptr,TSTRING,keyname,val,NULL,&sta))
@@ -189 +258 @@
   // Select the column to be read
   FitsABTColRd fbt(fof,"BoloMuv_28",0,1000,1,3);
+  FitsABTColRd fbt2(fof,"BoloMuv_29",0,1000,1,3);
   fbt.SetDebug(3);
   fbt.Print(3);
@@ -194 +264 @@
   for(long i=0;i<fbt.GetNbLine();i++) {
     double x = fbt.Read(i);
-    if(i%lpmod==0) cout<<i<<": "<<x<<endl;
+    double y = fbt2.Read(i);
+    if(i%lpmod==0) cout<<i<<": "<<x<<", "<<y<<endl;
   }
   // Read into a vector
@@ -203 +274 @@
   n = fbt.Read(10,-1,data);
     cout<<"Number of values read: "<<n<<endl;
+  TVector<double> data2;
+  fbt2.Read(32,50,data);
   // Close the fits file
   delete fof;
@@ -218 +291 @@
   \param lp : debug level
   \verbatim
-  - if ihdu=0 or ihdu>nhdu first binary or ASCII table is taken
+  - if ihdu<=0 first BINARY or ASCII table is taken
+  - if ihdu>nhdu ihdu is set to nhdu
   - bsens>0    read forward
     bsens<0    read backward
@@ -247 +321 @@
 {
   Init(fbt.GetFitsOpenFile(),fbt.GetColLabel().c_str()
-      ,fbt.GetColNum(),fbt.GetHDU()
+      ,fbt.GetColNum(),fbt.HDU()
       ,fbt.GetBLen(),fbt.GetBSens(),fbt.DbgLevel);
 }
@@ -254 +328 @@
 FitsABTColRd::FitsABTColRd()
 {
- FitsFN = "";
  ColLabel = ""; ColTUnit = ""; ColTForm = "";
  ColNum = -1; ColTypeCode = 0;
- IHdu = 0; NHdu = 0; HduType = 0;
  NBcol = 0; NBline = 0;
  SetNulVal(); SetDebug(0);
@@ -271 +343 @@
 {
  // Initialisation des Parametres Generaux
- FitsFN = "";
- ColLabel = collabel;
- ColTUnit = "";
- ColTForm = "";
- ColNum = colnum;
- ColTypeCode = 0;
- IHdu = ihdu;
- NHdu = 0;
- HduType = 0;
- NBcol = 0;
- NBline = 0;
- SetNulVal();
- SetDebug(lp);
+ ColLabel=collabel; ColTUnit=""; ColTForm=""; ColNum=colnum; ColTypeCode=0;
+ NBcol = 0; NBline = 0;
+ SetNulVal(); SetDebug(lp);
  NFitsRead = 0;
- FitsOF = NULL;
- FitsPtr = NULL;
+ FitsOF = NULL; FitsPtr = NULL;
  LineDeb = LineFin = -1;
  Buffer = NULL;
- ChangeBuffer(blen,bsens);
 
  // Caracteristiques du FitsOpenFile
  FitsOF = fof;
- if(FitsOF==NULL) {
-   Delete();
+ if(FitsOF==NULL)
    throw NullPtrError("FitsABTColRd::Init: FitsOpenFile pointer is NULL\n");
- }
+
  FitsPtr = FitsOF->GetFitsPtr();
- if(FitsPtr==NULL) {
-   Delete();
+ if(FitsPtr==NULL)
    throw NullPtrError("FitsABTColRd::Init: FitsPtr pointer is NULL\n");
- }
- NHdu = FitsOF->GetNHdu();
- if(DbgLevel>1) cout<<"FitsABTColRd::Init  NHdu="<<NHdu<<endl;
- if(NHdu<=0) {
-   Delete();
-   throw SzMismatchError("FitsABTColRd::Init: Bad NHdu\n");
- }
- FitsFN = FitsOF->GetFileName();
- if(FitsFN.size() <= 0 ) {
-   Delete();
-   throw ParmError("FitsABTColRd::Init: Fits file name error\n");
- }
 
  int sta = 0;
+ if(ihdu<0) ihdu=0; if(ihdu>NHDU()) ihdu=NHDU();
 
  // Get HDU for bin/ascii table
- // si IHdu <=0 || >NHdu on cherche la 1ere bin/ascii table
- // sinon on se positionne sur IHdu
- if(IHdu<=0 || IHdu>NHdu)
-   for(int ihdu=1;ihdu<=NHdu;ihdu++) {
-     if(fits_movabs_hdu(FitsPtr,ihdu,&HduType,&sta)) FitsOpenFile::printerror(sta);
-     if(DbgLevel>1) cout<<"...Init ihdu="
-                        <<ihdu<<" HduType="<<HduType<<endl;
-     if(HduType==BINARY_TBL || HduType==ASCII_TBL) {IHdu = ihdu; break;}
+ // ATTENTION: le fichier est ouvert mais non positionne sur un HDU,
+ // une classe utilisant ce fichier doit le positionner sur un HDU.
+ // Par contre, si une autre classe utilise ce meme FitsOpenFile,
+ // elle ne peut le positionner que sur ce meme HDU !
+ if(FitsOF->GetPosStatus()==false) {
+   if(ihdu==0) { // find the first BINARY then the first ASCII
+     int rc = FitsOF->MoveToFirst(BINARY_TBL);
+     if(rc!=BINARY_TBL) FitsOF->MoveToFirst(ASCII_TBL);
+   } else {
+     int rc = FitsOF->MoveToHDU(ihdu);
+     if(rc!=ihdu)
+       throw RangeCheckError("FitsABTColRd::Init: Error moving to requested HDU\n");
    }
- if(IHdu<=0 || IHdu>NHdu) {
-   cout<<"NO BINARY or ASCII hdu found"<<endl;
-   IHdu = 0; Delete();
-   throw TypeMismatchExc("FitsABTColRd::Init: NO BINARY or ASCII hdu found\n");
- }
- if(fits_movabs_hdu(FitsPtr,IHdu,&HduType,&sta)) {
-   FitsOpenFile::printerror(sta); Delete();
-   throw RangeCheckError("FitsABTColRd::Init: Error moving to requested HDU\n");
- }
- if(HduType!=BINARY_TBL && HduType!=ASCII_TBL) {
-   Delete();
+ } else { // Fits file has already been positionned
+   if(ihdu>0 && ihdu!=HDU())
+     throw RangeCheckError("FitsABTColRd::Init: file already posit. at another HDU\n");
+ }
+
+ // Check HDUType and set position status to TRUE
+ if(HDUType()!=BINARY_TBL && HDUType()!=ASCII_TBL)
    throw TypeMismatchExc("FitsABTColRd::Init: HDU not ASCII/BINARY table\n");
- }
+ if(DbgLevel>1) cout<<"...Init ihdu="<<ihdu<<" HduType="<<HDUType()<<endl;
+ FitsOF->SetPosStatus(true);
 
  // Get number of columns
  if(fits_get_num_cols(FitsPtr,&NBcol,&sta)) {
-   FitsOpenFile::printerror(sta); Delete();
+   FitsOpenFile::printerror(sta);
    throw NotAvailableOperation("FitsABTColRd::Init: Error getting number of columns\n");
  }
  if(DbgLevel>1) cout<<"...Init  NBcol="<<NBcol<<endl;
- if(NBcol<1) {
-   Delete();
+ if(NBcol<1)
    throw RangeCheckError("FitsABTColRd::Init: Bad number of colums\n");
- }
 
  // Get number of rows
  if(fits_get_num_rows(FitsPtr,&NBline,&sta)) {
-   FitsOpenFile::printerror(sta); Delete();
+   FitsOpenFile::printerror(sta);
    throw NotAvailableOperation("FitsABTColRd::Init: Error getting number of rows\n");
  }
  if(DbgLevel>1) cout<<"...Init  NBline="<<NBline<<endl;
- if(NBline<1) {
-   Delete();
+ if(NBline<1)
    throw RangeCheckError("FitsABTColRd::Init: Bad number of rows\n");
- }
 
  // Get column number
@@ -367 +411 @@
    strcpy(labelcol,ColLabel.c_str());
    if(fits_get_colnum(FitsPtr,CASESEN,labelcol,&ColNum,&sta)) {
-     FitsOpenFile::printerror(sta); Delete();
+     FitsOpenFile::printerror(sta);
      throw NotAvailableOperation("FitsABTColRd::Init: Error getting column name\n");
    }
@@ -373 +417 @@
  }
  if(DbgLevel>1) cout<<"...Init  ColNum="<<ColNum<<endl;
- if(ColNum<0 || ColNum>=NBcol) {
-   Delete();
+ if(ColNum<0 || ColNum>=NBcol)
    throw RangeCheckError("FitsABTColRd::Init: Bad column number\n");
- }
 
  // Get column type
  if(fits_get_coltype(FitsPtr,ColNum+1,&ColTypeCode,NULL,NULL,&sta)) {
-   FitsOpenFile::printerror(sta); Delete();
+   FitsOpenFile::printerror(sta);
    throw ParmError("FitsABTColRd::Init: Error getting column type\n");
  }
  if(DbgLevel>1) cout<<"...Init ColTypeCode="<<ColTypeCode<<endl;
  if(ColTypeCode==TSTRING || ColTypeCode==TCOMPLEX ||  ColTypeCode==TDBLCOMPLEX
-                         || ColTypeCode<0 ) {
-   Delete();
+                         || ColTypeCode<0 )
    throw ParmError("FitsABTColRd::Init: Selected column is not Numerical\n");
- }
 
  // Get column name back, tunit, tform
@@ -394 +434 @@
  long repeat=0; double tscale=1., tzero=0.;
  int rc=0;
- if(HduType==BINARY_TBL) {
+ if(HDUType()==BINARY_TBL) {
    fits_get_bcolparms(FitsPtr,ColNum+1,labelcol,tunit,tform
                      ,&repeat,&tscale,&tzero,NULL,tdisp,&sta);
@@ -402 +442 @@
  }
  if(rc) {
-   FitsOpenFile::printerror(sta); Delete();
+   FitsOpenFile::printerror(sta);
    throw RangeCheckError("FitsABTColRd::Init: Error getting the column caracteristics\n");
  }
@@ -408 +448 @@
  ColTUnit = tunit;
  ColTForm = tform;
+
+ // Set the buffer for reading
+ ChangeBuffer(blen,bsens);
 
  if(DbgLevel)
@@ -779 +822 @@
    <<" ncols="<<NBcol<<" nrows="<<NBline;
  if(lp>0) os<<" NRead="<<NFitsRead;
- os<<"\n... "<<FitsFN<<"["<<IHdu<<"/"<<NHdu<<"]"
-   <<"\n... Label["<<ColNum<<"]="<<ColLabel
-   <<" TypeCode="<<ColTypeCode
+ os<<"\n... "<<FileName()<<"["<<HDU()<<"/"<<NHDU()<<" type="<<HDUType()<<"]"
+   <<"\n... Label["<<ColNum<<"]="<<ColLabel<<" TypeCode="<<ColTypeCode
    <<" TUnit="<<ColTUnit<<" TForm="<<ColTForm
    <<endl;
@@ -828 +870 @@
   \param lp : debug level
   \verbatim
-  - if ihdu=0 or ihdu>nhdu first binary or ASCII table is taken
+  - if ihdu<=0 first BINARY or ASCII table is taken
+  - if ihdu>nhdu ihdu is set to nhdu
   - bsens>0    read forward
     bsens<0    read backward
@@ -878 +921 @@
  FitsOpenFile* fof = new FitsOpenFile(*fbt.GetFitsOpenFile());
  Init(fof,fbt.GetColLabel().c_str()
-     ,fbt.GetColNum(),fbt.GetHDU()
+     ,fbt.GetColNum(),fbt.HDU()
      ,fbt.GetBLen(),fbt.GetBSens(),fbt.DbgLevel);
 }
@@ -902 +945 @@
   Constructor.
   \param fof : Pointer to the Class for opening the FITS file
-  \param ihdu : number of the HDU where the column is.
+  \param ihdu : number of the HDU where the image is.
   \param lp : debug level
   \verbatim
-  - if ihdu=0 or ihdu>nhdu first binary or ASCII table is taken
+  - if ihdu<=0 first IMAGE hdu is taken
+  - if ihdu>nhdu ihdu is set to nhdu
   \endverbatim
   \warning ihdu = [1,nhdu]
@@ -917 +961 @@
 FitsImg2DRd::FitsImg2DRd(FitsImg2DRd& fbt)
 {
- Init(fbt.GetFitsOpenFile(),fbt.GetHDU(),fbt.DbgLevel);
+ Init(fbt.GetFitsOpenFile(),fbt.HDU(),fbt.DbgLevel);
 }
 
@@ -923 +967 @@
 FitsImg2DRd::FitsImg2DRd()
 {
- FitsFN = "";
- IHdu = 0; NHdu = 0; HduType = 0;
  Naxis[0] = Naxis[1] = 0;
  SetNulVal(); SetDebug(0);
@@ -934 +976 @@
 {
  // Initialisation des Parametres Generaux
- FitsFN = "";
- IHdu = ihdu; NHdu = 0; HduType = 0;
  Naxis[0] = Naxis[1] = 0;
  SetNulVal(); SetDebug(lp);
@@ -944 +984 @@
  if(FitsOF==NULL)
    throw NullPtrError("FitsImg2DRd::Init: FitsOpenFile pointer is NULL\n");
+
  FitsPtr = FitsOF->GetFitsPtr();
  if(FitsPtr==NULL)
    throw NullPtrError("FitsImg2DRd::Init: FitsPtr pointer is NULL\n");
- NHdu = FitsOF->GetNHdu();
- if(DbgLevel>1) cout<<"FitsImg2DRd::Init  NHdu="<<NHdu<<endl;
- if(NHdu<=0)
-   throw SzMismatchError("FitsImg2DRd::Init: Bad NHdu\n");
- FitsFN = FitsOF->GetFileName();
- if(FitsFN.size() <= 0 )
-   throw ParmError("FitsImg2DRd::Init: Fits file name error\n");
 
  int sta = 0;
+ if(ihdu<0) ihdu=0; if(ihdu>NHDU()) ihdu=NHDU(); 
 
  // Get HDU 2D image
- // si IHdu <=0 || >NHdu on cherche la 1ere image
- // sinon on se positionne sur IHdu
- if(IHdu<=0 || IHdu>NHdu)
-   for(int ihdu=1;ihdu<=NHdu;ihdu++) {
-     if(fits_movabs_hdu(FitsPtr,ihdu,&HduType,&sta)) FitsOpenFile::printerror(sta);
-     if(DbgLevel>1)
-       cout<<"...Init ihdu="<<ihdu<<" HduType="<<HduType<<endl;
-     if(HduType==IMAGE_HDU) {IHdu = ihdu; break;}
+ // ATTENTION: ... cf blabla equivalent dans FitsABTColRd::Init()
+ if(FitsOF->GetPosStatus()==false) {
+   if(ihdu==0) { // find the first IMAGE_HDU
+     FitsOF->MoveToFirst(IMAGE_HDU);
+   } else {
+     int rc = FitsOF->MoveToHDU(ihdu);
+     if(rc!=ihdu)
+       throw RangeCheckError("FitsABTColRd::Init: Error moving to requested HDU\n");
    }
- if(IHdu<=0 || IHdu>NHdu) {
-   cout<<"NO IMAGE_HDU hdu found"<<endl;
-   IHdu = 0;
-   throw TypeMismatchExc("FitsImg2DRd::Init: NO IMAGE_HDU hdu found\n");
- }
- if(fits_movabs_hdu(FitsPtr,IHdu,&HduType,&sta)) {
-   FitsOpenFile::printerror(sta);
-   throw RangeCheckError("FitsImg2DRd::Init: Error moving to requested HDU\n");
- }
- if(HduType!=IMAGE_HDU)
+ } else { // Fits file has already been positionned
+   if(ihdu>0 && ihdu!=HDU())
+     throw RangeCheckError("FitsABTColRd::Init: file already posit. at another HDU\n");
+ }
+
+ // Check HDUType and set position status to TRUE
+ if(HDUType()!=IMAGE_HDU)
    throw TypeMismatchExc("FitsImg2DRd::Init: HDU not IMAGE_HDU\n");
+ FitsOF->SetPosStatus(true);
 
  // Get NAXIS 1 et 2
@@ -986 +1019 @@
  }
  if(DbgLevel>1)
-   cout<<"...Init(hdu="<<IHdu<<")  NAXIS1="<<Naxis[0]<<" NAXIS2="<<Naxis[1]<<" (nfound="<<nfound<<")"<<endl;
+   cout<<"...Init(hdu="<<HDU()<<")  NAXIS1="<<Naxis[0]<<" NAXIS2="
+       <<Naxis[1]<<" (nfound="<<nfound<<")"<<endl;
  if(nfound!=2 || Naxis[0]<=0 || Naxis[1]<=0)
    throw NotAvailableOperation("FitsImg2DRd::Init: bad Naxis[0-1] value\n");

trunk/SophyaExt/FitsIOServer/fabtcolread.h

@@ -24 +24 @@
   virtual ~FitsOpenFile();
 
-  inline string GetFileName() {return FitsFN;}
-  inline int GetNHdu() {return NHdu;}
-  inline fitsfile* GetFitsPtr() {return FitsPtr;}
+  inline string FileName() const {return FitsFN;}
+  //! Get the number of the HDU read
+  inline int HDU(void) const {return IHdu;}
+  //! Get the number of the HDU type
+  inline int HDUType(void) const {return HduType;}
+  //! Get the number of HDU in file
+  inline int NHDU() const {return NHdu;}
+  //! Get the CFISTIO fits file pointer
+  inline fitsfile* GetFitsPtr() const {return FitsPtr;}
+  //! Set the positionning status of the file
+  inline void SetPosStatus(bool sta=true) {HasBeenPos = sta;}
+  //! Get the positionning status of the file
+  inline bool GetPosStatus(void) const {return HasBeenPos;}
+
+  int MoveToHDU(int ihdu);
+  int MoveToFirst(int hdutype,int ihdudeb=1);
+  int MoveToLast(int hdutype,int ihdudeb=1);
+  void Print(void);
 
   static double ReadKey(fitsfile *fitsptr,char *keyname);
@@ -38 +53 @@
 
   string FitsFN;
-  int NHdu;
+  int NHdu, IHdu, HduType;
   fitsfile *FitsPtr;
+  bool HasBeenPos;
 };
 
@@ -62 +78 @@
   double Read(long n,bool usebuffer=true);
 
-  long   Read(long n1,long n2,TVector<uint_2>& data);
-  long   Read(long n1,long n2,TVector<int_4>& data);
-  long   Read(long n1,long n2,TVector<int_8>& data);
-  long   Read(long n1,long n2,TVector<float>& data);
-  long   Read(long n1,long n2,TVector<double>& data);
+  long Read(long n1,long n2,TVector<uint_2>& data);
+  long Read(long n1,long n2,TVector<int_4>& data);
+  long Read(long n1,long n2,TVector<int_8>& data);
+  long Read(long n1,long n2,TVector<float>& data);
+  long Read(long n1,long n2,TVector<double>& data);
 
   //! return the value of the first row
@@ -79 +95 @@
 
   //! Set debug level
-  inline void    SetDebug(int lp=0) {DbgLevel = (unsigned short) lp;}
+  inline void SetDebug(int lp=0) {DbgLevel = (unsigned short) lp;}
   //! Set null value to be return when reading null data (0=return the data)
-  inline void    SetNulVal(double nulval=0.) {NulVal = nulval;}
+  inline void SetNulVal(double nulval=0.) {NulVal = nulval;}
   //! Get the FITS file name
-  inline string  GetFileName(void) {return FitsFN;}
+  inline string FileName(void) const
+         {if(FitsOF) return FitsOF->FileName(); else return (string)"";}
   //! Get the pointer to FitsOpenFile
-  inline FitsOpenFile* GetFitsOpenFile(void) {return FitsOF;}
+  inline FitsOpenFile* GetFitsOpenFile(void) const {return FitsOF;}
   //! Get the FITS file pointer (cfistio pointer)
-  inline fitsfile* GetFitsPointer(void) {return FitsPtr;}
+  inline fitsfile* GetFitsPtr(void) const {return FitsPtr;}
   //! Get the number of HDU in the FITS file
-  inline int     GetNHDU(void) {return NHdu;}
+  inline int NHDU(void) const
+         {if(FitsOF) return FitsOF->NHDU(); else return 0;}
   //! Get the number of the HDU read
-  inline int     GetHDU(void) {return IHdu;}
+  inline int HDU(void) const
+         {if(FitsOF) return FitsOF->HDU(); else return 0;}
   //! Get the HDU type
-  inline int     GetHDUType(void) {return HduType;}
+  inline int HDUType(void) const
+         {if(FitsOF) return FitsOF->HDUType(); else return 0;}
   //! Get the number of rows in the FITS HDU to be read
-  inline long    GetNbLine(void) {return NBline;}
+  inline long GetNbLine(void) const {return NBline;}
   //! Get the number of columns in the FITS HDU to be read
-  inline int     GetNbCol(void) {return NBcol;}
+  inline int GetNbCol(void) const {return NBcol;}
   //! Get the columns number that is read
-  inline int     GetColNum(void) {return ColNum;}
+  inline int GetColNum(void) const {return ColNum;}
   //! Get the columns label that is read
-  inline string  GetColLabel(void) {return ColLabel;}
+  inline string GetColLabel(void) const {return ColLabel;}
   //! Get the columns type code that is read
-  inline int     GetColTypeCode(void) {return ColTypeCode;}
+  inline int GetColTypeCode(void) const {return ColTypeCode;}
   //! Get the columns fits tunit that is read
-  inline string  GetColTUnit(void) {return ColTUnit;}
+  inline string GetColTUnit(void) const {return ColTUnit;}
   //! Get the columns fits tform that is read
-  inline string  GetColTForm(void) {return ColTForm;}
+  inline string GetColTForm(void) const {return ColTForm;}
   //! Get the read requested buffer length
-  inline long    GetBLen(void) {return BuffLen;}
+  inline long GetBLen(void) const {return BuffLen;}
   //! Get the read buffer direction
-  inline long    GetBSens(void) {return BuffSens;}
+  inline long GetBSens(void) const {return BuffSens;}
   //! Print to os
-  virtual void   Print(ostream& os,int lp=1) const;
+  virtual void Print(ostream& os,int lp=1) const;
   //! Print to stdout
-  inline  void   Print(int lp=1) const {Print(cout,lp);}
+  inline  void Print(int lp=1) const {Print(cout,lp);}
   //! Get the read effective buffer length
-  inline long    GetNBuffer(void) {return NBuffer;}
+  inline long GetNBuffer(void) const {return NBuffer;}
   //! Get the read bufferpointer
   inline double* GetBuffer(void) {return Buffer;}
@@ -126 +146 @@
   void Delete(void);
 
-  string FitsFN,ColLabel,ColTUnit,ColTForm;
-  int ColNum,ColTypeCode,IHdu,NHdu,HduType,NBcol;
+  string ColLabel,ColTUnit,ColTForm;
+  int ColNum,ColTypeCode,NBcol;
   long NBline;
 
@@ -136 +156 @@
   unsigned long NFitsRead;
   FitsOpenFile* FitsOF;
-  fitsfile* FitsPtr;
+  fitsfile* FitsPtr; // Redite avec FitsOF->FitsPtr, mais utilise partout!
   long LineDeb, LineFin;
   double *Buffer;
@@ -175 +195 @@
   string ReadKeyS(char *keyname);
 
-  long   Read(TMatrix<uint_2>& data);
-  long   Read(TMatrix<int_4>& data);
-  long   Read(TMatrix<int_8>& data);
-  long   Read(TMatrix<float>& data);
-  long   Read(TMatrix<double>& data);
+  long Read(TMatrix<uint_2>& data);
+  long Read(TMatrix<int_4>& data);
+  long Read(TMatrix<int_8>& data);
+  long Read(TMatrix<float>& data);
+  long Read(TMatrix<double>& data);
 
   //! Set debug level
@@ -186 +206 @@
   inline void    SetNulVal(double nulval=0.) {NulVal = nulval;}
   //! Get the pointer to FitsOpenFile
-  inline FitsOpenFile* GetFitsOpenFile(void) {return FitsOF;}
+  inline FitsOpenFile* GetFitsOpenFile(void) const {return FitsOF;}
+  //! Get the FITS file pointer (cfistio pointer)
+  inline fitsfile* GetFitsPtr(void) const {return FitsPtr;}
+  //! Get the number of HDU in the FITS file
+  inline int NHDU(void) const
+         {if(FitsOF) return FitsOF->NHDU(); else return 0;}
   //! Get the number of the HDU read
-  inline int     GetHDU(void) {return IHdu;}
+  inline int HDU(void) const
+         {if(FitsOF) return FitsOF->HDU(); else return 0;}
   //! Get the HDU type
-  inline int     GetHDUType(void) {return HduType;}
+  inline int HDUType(void) const
+         {if(FitsOF) return FitsOF->HDUType(); else return 0;}
   //! Get NAXIS1
-  inline long    GetNaxis1(void) {return Naxis[0];}
+  inline long Naxis1(void) const {return Naxis[0];}
   //! Get NAXIS2
-  inline long    GetNaxis2(void) {return Naxis[1];}
+  inline long Naxis2(void) const {return Naxis[1];}
 
 protected:
   void Init(FitsOpenFile* fof,int ihdu,int lp);
 
-  string FitsFN;
-  int IHdu,NHdu,HduType;
   long Naxis[2];
-
   double NulVal;
   unsigned short DbgLevel;
 
   FitsOpenFile* FitsOF;
-  fitsfile* FitsPtr;
-};
+  fitsfile* FitsPtr; // Redite avec FitsOF->FitsPtr, mais utilise partout!
+};