GainVoltPmtParameters.h

#ifndef GainVoltPmtParameters_H_INCLUDED
#define GainVoltPmtParameters_H_INCLUDED
#include <Stiostream.h>
#include <vector>
using std::vector;
#include "PowerLawFit.h"
#include "PmtIdentifier.h"

class GainVoltPmtParameters 
{
public:
  GainVoltPmtParameters();
  GainVoltPmtParameters(double refVoltage, double gains[5]);
  GainVoltPmtParameters(int n, double volts[], double gains[]);
  GainVoltPmtParameters(const GainVoltPmtParameters&parameters);
  ~GainVoltPmtParameters();
  GainVoltPmtParameters & operator=(const GainVoltPmtParameters&parameters);

  void set(double refVolt, double gains[5]);
  void set(int n, double volts[], double gains[]);

  friend ostream& operator<<(ostream& os, GainVoltPmtParameters& object);
  friend istream& operator>>(istream& is, GainVoltPmtParameters& object);
  void fit();
  double getMultConstant() const;
  double getExponent() const;
  double getGain(double voltage) const;
  double getVoltage(double gain) const;
  int getNPoints() const;
  void print();
  PmtIdentifier & getPmtIdentifier();

  static void setDefaults(double multCoefficient, double exponent);

  bool isValid() const;

  static int ioMode;
protected:
  PmtIdentifier _id;
  double _a;
  double _b;
  map<double,double> _data;
  PowerLawFit<double> _fit;

  static double _defaultA;
  static double _defaultB;

};


#endif