BaselinesSet.h

Go to the documentation of this file.

#if     !defined(_BASELINESSET_H)
 
#include "ASDMEntities.h"
#include "Error.h"
#include "BaselineMetadata.h"
#include "SwitchCyclesList.h"
#include "DataDescriptionsSet.h"

#include <AxisName.h>
#include <DataContent.h>


namespace sdmbin {
  
  class BaselinesSet : public DataDescriptionsSet
  {
  public:
    BaselinesSet();

    BaselinesSet( vector<Tag>  v_antennaIdArray, 
                  vector<int>  v_feedIdArray,
                  vector<int>  v_phasedArrayList,
                  vector<bool> v_antennaUsedArray, 
                  DataDescriptionsSet& dataDescriptionsSet
                  );

    BaselinesSet(const BaselinesSet &);

    ~BaselinesSet();

    unsigned int transferId(unsigned int na,  unsigned int ndd, unsigned int nbin);

    unsigned int transferId(unsigned int na,  unsigned int nfe, unsigned int ndd, unsigned int nbin);
    
    unsigned int transferId(unsigned int na1, unsigned int na2, unsigned int ndd, unsigned int nbin, unsigned int napc);

    unsigned int transferId(unsigned int na1, unsigned int na2, unsigned int nfe, unsigned int ndd,  unsigned int nbin, unsigned int napc);

    unsigned int baselineIndex( Tag antennaId)  ; // throw (Error);

    unsigned int baselineIndex( Tag antennaId1, Tag antennaId2) ; // throw (Error);

    unsigned int baselineIndex( unsigned int na1, unsigned int na2) ; // throw(Error);

    unsigned int antenna1(unsigned int baselineIndex)  ; // throw (Error);

    unsigned int antenna2(unsigned int baselineIndex)  ; // throw (Error);

    unsigned int feedIndex(Tag antennaId, int feedId) ; // throw (Error);

    Tag         getEffAntennaId(unsigned int na) ; // throw (Error);

    unsigned int getNumAntennas();


    int         getFeedId(unsigned int na, unsigned int nfe) ; // throw (Error);

    unsigned int getNumEffAntennas();

    unsigned int getNumBaselines();

    unsigned int getNumEffBaselines();
    
    unsigned int getNumPDTvalue(Enum<DataContent> e_dc, EnumSet<AxisName> es_an, bool effective);

  private:
    vector<Tag>  v_AntennaIdArray_;   
    vector<int>  v_FeedIdArray_;      
    vector<int>  v_PhasedArrayList_;  
    unsigned int numFeeds_;           
    unsigned int numAntennas_;        
    unsigned int numEffAntennas_;     
    unsigned int numBaselines_;       
    unsigned int numEffBaselines_;    
    vector<Tag>  v_effAntennaIdArray_;

  };


  inline unsigned int BaselinesSet::transferId(unsigned int na1,unsigned int na2,unsigned int nfe, unsigned int ndd,unsigned int nbin,unsigned int napc){

    unsigned int baselineidx = baselineIndex(na1,na2);                   //cout << "baselineidx=" << baselineidx << endl;

    unsigned int v_cumulCrossSize_ndd=0; if(v_cumulCrossSize_.size()!=0)v_cumulCrossSize_ndd = v_cumulCrossSize_[ndd];
    unsigned int v_crossSize_ndd=0;      if(v_crossSize_.size()!=0)     v_crossSize_ndd      = v_crossSize_[ndd];

    unsigned int n = 
      baselineidx*numFeeds_*sumCrossSize_ +                           // /baseline
      nfe*sumCrossSize_ +                                             // /baseline/feed     
      v_cumulCrossSize_ndd +                                          // /baseline/feed/datadesc
      nbin*(v_crossSize_ndd/numBin(ndd)) +                            // /baseline/feed/datadesc/bin
      napc*(v_crossSize_ndd/(numBin(ndd)*numApc_));                   // /baseline/feed/datadesc/bin/napc

    return n;
  }

  inline unsigned int DataDescriptionsSet::numBin(unsigned int ndd) { // throw(Error){
    // if(ndd>=v_basebandName_.size())
    //   Error(SERIOUS,"The dataDescription index must be smaller than %d",
    //      v_basebandName_.size());
    return v_numBin_[ndd];  
  }

  inline unsigned int BaselinesSet::baselineIndex( unsigned int na1, unsigned int na2) { // throw(Error){
    unsigned int iIdx = min (na1, na2);
    unsigned int jIdx = max (na1, na2);
    return jIdx * (jIdx - 1) / 2 + iIdx;
  }
}
#define _BASELINESSET_H
#endif