InterferometryCache.h

#ifndef INTERFEROMETRY_CACHE_H
#define INTERFEROMETRY_CACHE_H

#include "AnitaConventions.h"


#include <vector>


namespace Acclaim
{


  class AnalysisReco;
  class CrossCorrelator;
  
  class InterferometryCache {

    friend class InterferometricMap;

  public:
    InterferometryCache();

    InterferometryCache(CrossCorrelator* cc, const AnalysisReco* reco);

    void init(CrossCorrelator* cc, const AnalysisReco* reco, bool forceCacheRecalculation = false);
  
    inline int coarseIndex(int pol, int combo, int phiBin, int thetaBin){  
      return ((pol*fNumCombos + combo)*fNCoarseBinsPhi + phiBin)*fNCoarseBinsTheta + thetaBin;
    }

    inline double coarseDt(int pol, int combo, int phiBin, int thetaBin){
      return fDeltaTs[coarseIndex(pol, combo, phiBin, thetaBin)];
    }

    inline int partBAsIndex(int pol, int combo, int fineThetaBin){
      return (pol*fNumCombos + combo)*fNFineBinsTheta + fineThetaBin;
    }

    inline int zoomedCosPartIndex(int pol, int ant, int finePhiBin){
      return (pol*NUM_SEAVEYS + ant)*fNFineBinsPhi + finePhiBin;
    }

    inline int part21sIndex(int pol, int combo, int finePhiBin){
      return (pol*fNumCombos + combo)*fNFineBinsPhi + finePhiBin;
    }

    void populateCache(CrossCorrelator* cc, const AnalysisReco* reco);

    void populateFineCache(CrossCorrelator* cc, const AnalysisReco* reco);


    int fUseOffAxisDelay; 
    bool fInitialized; 
    int fNCoarseBinsPhi; 
    int fNCoarseBinsTheta; 
    int fNumCombos; 
    std::vector<double> fDeltaTs; 

    int fNFineBinsPhi; 
    int fNFineBinsTheta; 

    //[AnitaPol::kNotAPol][NUM_COMBOS][NUM_BINS_THETA_ZOOM_TOTAL]
    std::vector<double> fPartBAsZoom; 

    //[AnitaPol::kNotAPol][NUM_COMBOS][NUM_BINS_PHI_ZOOM_TOTAL]
    std::vector<double> fPart21sZoom; 

    // NUM_BINS_PHI_ZOOM_TOTAL
    std::vector<double> fZoomedThetaWaves; 

    // NUM_BINS_THETA_ZOOM_TOTAL
    std::vector<double> fZoomedTanThetaWaves; 

    // NUM_BINS_THETA_ZOOM_TOTAL
    std::vector<double> fZoomedCosThetaWaves; 

    // NUM_BINS_THETA_ZOOM_TOTAL
    std::vector<double> fDtFactors; 

    // NUM_BINS_PHI_ZOOM_TOTAL
    std::vector<double> fZoomedPhiWaveLookup; 

    // NUM_BINS_PHI_ZOOM_TOTAL
    std::vector<double> fZoomedCosPartLookup; 

    // [AnitaPol::kNotAPol][NUM_COMBOS][NUM_BINS_PHI_ZOOM_TOTAL]
    std::vector<double> fOffAxisDelays; 

    // [AnitaPol::kNotAPol][NUM_COMBOS][NUM_BINS_PHI_ZOOM_TOTAL]    
    std::vector<double> fOffAxisDelaysDivided; 
  };
}

#endif