DeltaRM.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2005 - 2017 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/Polarimetry/Pulsar/DeltaRM.h
 
#ifndef __Pulsar_DeltaRM_h
#define __Pulsar_DeltaRM_h
 
#include "Pulsar/PhaseWeight.h"
#include "Pulsar/FrequencyIntegrate.h"
 
#include "Estimate.h"
 
namespace Pulsar {
 
  class Archive;

  class DeltaRM  {
 
  public:

    DeltaRM ();

    ~DeltaRM ();

    void set_rotation_measure (const Estimate<double>& rm)
    { rotation_measure = rm; }

    Estimate<double> get_rotation_measure () const
    { return rotation_measure; }

    unsigned get_used_bins () { return used_bins; }

    void set_data (const Archive*);

    void set_subint (unsigned idx) { isubint = idx; }
    unsigned get_subint () const { return isubint; }

    void refine ();

    void set_threshold (float t) { threshold = t; }
    float get_threshold () const { return threshold; }

    void set_onpulse (PhaseWeight* w) { onpulse_weights = w; }
    
    void set_include (const std::vector<unsigned>& bins);

    void set_exclude (const std::vector<unsigned>& bins);

    void set_policy (FrequencyIntegrate::Divided*);
 
  protected:
 
    float threshold;
    unsigned used_bins;

    std::vector<unsigned> include_bins;

    std::vector<unsigned> exclude_bins;

    Estimate<double> rotation_measure;

    Reference::To<const Archive> data;

    unsigned isubint = 0;

    Reference::To<PhaseWeight> onpulse_weights;

    Reference::To<FrequencyIntegrate::Divided> policy;
  };
 
}
 
#endif