ChannelSubsetMatch.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2008 by Paul Demorest
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __Pulsar_ChannelSubsetMatch_h
#define __Pulsar_ChannelSubsetMatch_h
 
#include "Pulsar/Config.h"
#include "Pulsar/Database.h"
 
namespace Pulsar {
 
  class Archive;

  class ChannelSubsetMatch {
 
  public:

    ChannelSubsetMatch ();

    ~ChannelSubsetMatch ();

    bool match (const Archive* super, const Archive* sub);

    bool match (const Database::StaticEntry* super,
                const Database::StaticEntry* sub);

    int super_channel (const Archive* super, const Archive* sub, int subchan);

    int sub_channel (const Archive* super, const Archive* sub, int superchan);

    std::string get_reason() const { return reason; }
 
    template<class Container>
    unsigned match_channel (const Container* container, double frequency,
        double tolerance=0.0);
 
  protected:

    std::string reason;

    double freq_tol;

    double bw_tol;
  };
 
}
 
template<class Container>
unsigned Pulsar::ChannelSubsetMatch::match_channel (const Container* container,
                                                    double frequency,
                                                    double tolerance)
{
  if (tolerance==0.0)
  {
    // Use default value for class
    tolerance = freq_tol * frequency;
  }
  for (unsigned j=0; j<container->get_nchan(); j++)
  {
    double cfreq = container->get_centre_frequency(j);
    if (fabs(frequency-cfreq) < tolerance)
      return j;
  }
 
  throw Error (InvalidParam, "Pulsar::ChannelSubsetMatch::match_channel",
               "no match for frequency=%lf", frequency);
}
 
#endif