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 (Archive*);
 
    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<Archive> data;
 
    Reference::To<PhaseWeight> onpulse_weights;
 
    Reference::To<FrequencyIntegrate::Divided> policy;
  };
 
}
 
#endif