PulserPinScale.cxx

Go to the documentation of this file.

// $Id: PulserPinScale.cxx,v 1.4 2004/11/15 13:03:07 west Exp $
//
// PulserPinScale
//
// Dbi object for raw drift point info
//
// Author: Phil Adamson 7/2004
//
// $Log: PulserPinScale.cxx,v $
// Revision 1.4  2004/11/15 13:03:07  west
// gcc 3.4: Replace DbiReader/DbiWriter.cxx -> .tpl and move after CVSID
//
// Revision 1.3  2004/11/04 22:17:37  phil.adamson
// Minor fix for new calib
//
// Revision 1.2  2004/11/04 09:16:08  west
// Add dummy versions of missing ctors and dtors
//
// Revision 1.1  2004/11/04 04:22:11  phil.adamson
// Extra bits for new calibration
//
//

#include "MessageService/MsgService.h"
#include "PulserCalibration/PulserPinScale.h"
#include "DatabaseInterface/DbiResultSet.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include "TGraphErrors.h"
#include "TF1.h"


ClassImp(PulserPinScale)

  //........

CVSID("$Id: PulserPinScale.cxx,v 1.4 2004/11/15 13:03:07 west Exp $");
  
#include "DatabaseInterface/DbiResultPtr.tpl"
  template class DbiResultPtr<PulserPinScale>;
#include "DatabaseInterface/DbiWriter.tpl"
  template class DbiWriter<PulserPinScale>;

//  To be replaced: dummy ctor and dtor just to allow linking.
PulserPinScale::PulserPinScale() {}
PulserPinScale::~PulserPinScale() {}

PulserPinScale::PulserPinScale(PulserGainPin &high, PulserGainPin &low, 
                               int aggNo)
{
  fAggregateNo = aggNo;
  // Do fit here
  assert(high.GetNumPoints()==low.GetNumPoints());
  const Float_t *x_high=high.GetMean();
  const Float_t *ex_high=high.GetError();
  const Float_t *y_low = low.GetMean();  
  const Float_t *ey_low = low.GetError();  

  TGraphErrors g(high.GetNumPoints(),x_high,y_low,ex_high,ey_low);
  float xmin=300;
  float xmax = 8000;

  for (int i=0;i<high.GetNumPoints();i++) {
    if (y_low[i] < 300 && x_high[i]>300) xmin =x_high[i]+1;
    if (y_low[i] > 8000 && x_high[i] < 8000) {
      xmax = x_high[i]-1;
      break;
    }
  }

  TF1 f("line","[0]*x",xmin,xmax);
  g.Fit("line","QR");
  MSG("Pulser",Msg::kDebug)<<"Aggregate "<<aggNo<<" fit slope "<<f.GetParameter(0)<<endl;

  fHighGainScale = 1;
  fLowGainScale = 1.0/f.GetParameter(0);

}


PulserPinScale::PulserPinScale(float high, float low, int aggNo)
{
  fHighGainScale = high;
  fLowGainScale = low;
  fAggregateNo = aggNo;
}

void PulserPinScale::Fill(DbiResultSet& rs,
                          const DbiValidityRec* /* vrec */)
{
  rs >> fAggregateNo >> fHighGainScale >> fLowGainScale;

} 

void PulserPinScale::Store(DbiOutRowStream& ors,
                          const DbiValidityRec* /* vrec */) const
{
    ors << fAggregateNo 
        << fHighGainScale 
        << fLowGainScale;
}


Float_t PulserPinScale::ScalePins(Float_t high, Float_t low) const
{
  if (fHighGainScale<=0||high<0) {
    return fLowGainScale*low;
  }
  if (fLowGainScale<=0||low<0) {
    return fHighGainScale*high;
  }
  if (high>8000) {
    return fLowGainScale*low;
  }
  if (low<500 || low > 8000) {
    return fHighGainScale*high;
  }
  if (high< 500) {
    return fLowGainScale*low;
  }
  return 0.5*(fLowGainScale*low + fHighGainScale*high);

}

---