CalPulserFits.cxx

Go to the documentation of this file.

// $Id: CalPulserFits.cxx,v 1.4 2006/03/29 12:41:38 cpw1 Exp $
//
// CalPulserFits
//
// Dbi object for results of linear fits to pulser gain data.
//
// Author: Pat Ward 13/09/2005
//
// $Log: CalPulserFits.cxx,v $
// Revision 1.4  2006/03/29 12:41:38  cpw1
// Mods for 60-point gain curve
//
// Revision 1.3  2006/01/25 12:47:03  cpw1
// New definition of database table with more columns
//
// Revision 1.2  2005/10/27 09:21:43  cpw1
// Add number of fit points
//
// Revision 1.1  2005/09/28 09:33:11  cpw1
// Replacement for PulserXScale
//
//

#include "MessageService/MsgService.h"
#include "Calibrator/CalPulserFits.h"

ClassImp(CalPulserFits)

CVSID("$Id: CalPulserFits.cxx,v 1.4 2006/03/29 12:41:38 cpw1 Exp $\n  \
      CVSID_DBIRESULTPTR ");

#include "DatabaseInterface/DbiResultPtr.tpl"
  template class DbiResultPtr<CalPulserFits>;
#include "DatabaseInterface/DbiWriter.tpl"
  template class DbiWriter<CalPulserFits>;

// Default constructor and destructor
CalPulserFits::CalPulserFits() {}
CalPulserFits::~CalPulserFits() {}

// Old constructor which fills row with values
CalPulserFits::CalPulserFits(Int_t aggNo, Int_t seid, Int_t npfit, Double_t slope, Double_t slopeerr, Double_t xoffset, Double_t chisq) : 
  fAggregateNo(aggNo),
  fStripEnd(seid),
  fNPFit(npfit),
  fSlope(slope),
  fSlopeErr(slopeerr), 
  fXOffset(xoffset),
  fChisq(chisq) 

{
  // Zero variables not used with this constructor
  fFitType = 0;
  fPar1 = 0;
  fPar2 = 0;
  fPar3 = 0;
  fPar4 = 0;
  fMeanRes = 0;
  fMaxRes = 0;
  fAdcMax = 0;
  fXtra1 = 0;
  fXtra2 = 0;
}

// New constructor which fills row with values
CalPulserFits::CalPulserFits(Int_t aggNo, Int_t seid, Int_t ftype, TF1* fit, Double_t meanres, Double_t maxres, Double_t adcmax) : 
  fAggregateNo(aggNo),
  fStripEnd(seid),
  fFitType(ftype), 
  fMeanRes(meanres), 
  fMaxRes(maxres),
  fAdcMax(adcmax)

{
  // No fit
  if(ftype < -25 || fit == 0) {
    fNPFit = 0;
    fSlope = 0;
    fSlopeErr = 0;
    fXOffset = 0;
    fPar1 = 0;
    fPar2 = 0;
    fPar3 = 0;
    fPar4 = 0;
    fChisq = 0;
  }

  // Line fit only
  if(abs(ftype) == 1) {
    fNPFit = fit->GetNumberFitPoints();
    fSlope = fit->GetParameter(0);
    fSlopeErr = fit->GetParError(0);
    fXOffset = fit->GetParameter(1);
    fPar1 = 0;
    fPar2 = 0;
    fPar3 = 0;
    fPar4 = 0;
    fChisq = fit->GetChisquare();
  }

  if(abs(ftype) == 5 || abs(ftype) == 25) {
    fNPFit = fit->GetNumberFitPoints();
    fSlope = fit->GetParameter(7);
    fSlopeErr = fit->GetParError(7);
    if(fSlope != 0) {
      fSlope = 1/fSlope;
      fSlopeErr = fSlope*fSlope*fSlopeErr;
    }
    fXOffset = fit->GetParameter(6);
    fPar1 = fit->GetParameter(1);
    fPar2 = fit->GetParameter(2);
    if(abs(ftype) == 25) {
      fPar3 = fit->GetParameter(3);
      fPar4 = fit->GetParameter(4);
    }
    else {
      fPar3 = 0;
      fPar4 = 0;
    }
    fChisq = fit->GetChisquare();
  }

  fXtra1 = 0;
  fXtra2 = 0;
}

void CalPulserFits::Fill(DbiResultSet& rs,
                          const DbiValidityRec* /* vrec */)
{
  rs >> fAggregateNo 
     >> fStripEnd 
     >> fFitType
     >> fNPFit
     >> fSlope 
     >> fSlopeErr 
     >> fXOffset 
     >> fPar1
     >> fPar2
     >> fPar3
     >> fPar4
     >> fChisq
     >> fMeanRes
     >> fMaxRes
     >> fAdcMax
     >> fXtra1 
     >> fXtra2;

} 

void CalPulserFits::Store(DbiOutRowStream& ors,
                          const DbiValidityRec* /* vrec */) const
{
    ors << fAggregateNo 
        << fStripEnd
        << fFitType
        << fNPFit
        << fSlope  
        << fSlopeErr 
        << fXOffset 
        << fPar1
        << fPar2
        << fPar3
        << fPar4
        << fChisq
        << fMeanRes
        << fMaxRes
        << fAdcMax
        << fXtra1 
        << fXtra2;
}

---