PulserGainLog.cxx

Go to the documentation of this file.

// $Id: PulserGainLog.cxx,v 1.4 2007/03/01 17:06:40 rhatcher Exp $
//
// PulserGainLog
//
// Dbi object
//
// Author: Phil Adamson 7/2003
//
// $Log: PulserGainLog.cxx,v $
// Revision 1.4  2007/03/01 17:06:40  rhatcher
// fix a bunch of CVSID lines.  The quoted string must be "$Id: $" (especially
// common is to leave off the trailing $).
//
// Revision 1.3  2004/10/21 10:20:16  west
// gcc 3.4: Replace DbiReader/DbiWriter.cxx -> .tpl and move after CVSID
//
// Revision 1.2  2004/10/05 21:39:48  murgia
// Fixed code to write points with low light levels. Added call to EndFile.
//
// Revision 1.1  2004/09/08 21:50:51  pa
// Checkin of new PulserCalibration version. MIght even work
//
//

#include "PulserCalibration/PulserGainLog.h"
#include "MessageService/MsgService.h"
#include "DatabaseInterface/DbiResultSet.h"
#include "DatabaseInterface/DbiValidityRec.h"
ClassImp(PulserGainLog)

CVSID("$Id: PulserGainLog.cxx,v 1.4 2007/03/01 17:06:40 rhatcher Exp $");

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


PulserGainLog::PulserGainLog() : fTimeStart(0,0) , fTimeEnd(0,0) 
{
  fCodeVer = 0;
  fAggregateNo = 0;
  fFirst = fLast = fNumPoints = 0;
  for (int i=0;i<40;i++) fHeight[i] = fWidth[i] = fNum[i] = fTriggers[i] = 0;
}

PulserGainLog::PulserGainLog(Int_t codever) : fCodeVer(codever), 
  fTimeStart(0,0) , fTimeEnd(0,0) 
{
  fFirst = fLast = fNumPoints = fAggregateNo = 0;
  for (int i=0;i<40;i++) fHeight[i] = fWidth[i] = fNum[i] = fTriggers[i] = 0;
}

 
void PulserGainLog::Fill(DbiResultSet& rs,
                          const DbiValidityRec* /* vrec */)
{
  rs >> fAggregateNo;
  rs >> fCodeVer;
  rs >> fIdent;
  rs >> fNumPoints;
  rs >> fFirst >> fLast;
  rs >> fTimeStart >> fTimeEnd;
  for (int i=0;i<40;i++) rs >> fHeight[i];
  for (int i=0;i<40;i++) rs >> fWidth[i];
  for (int i=0;i<40;i++) rs >> fNum[i];
  for (int i=0;i<40;i++) rs >> fTriggers[i];
} 

void PulserGainLog::Store(DbiOutRowStream& ors,
                           const DbiValidityRec* /* vrec */) const
{
  ors << fAggregateNo;
  ors << fCodeVer;
  ors << fIdent;
  ors << fNumPoints;
  ors << fFirst << fLast;
  ors << fTimeStart << fTimeEnd;
  for (int i=0;i<40;i++) ors << fHeight[i];
  for (int i=0;i<40;i++) ors << fWidth[i];
  for (int i=0;i<40;i++) ors << fNum[i];
  for (int i=0;i<40;i++) ors << fTriggers[i];
}


bool PulserGainLog::IsComplete() const
{
  MSG("Pulser",Msg::kDebug)<<"fFirst, fNumPoints "<< fFirst << ", " << fNumPoints <<endl;
    return ((fFirst==1)&&(fLast==fNumPoints));
  //  return ((fFirst==1)&&(fLast==20));
}


bool PulserGainLog::Add(const PulserGainLog &pgl)
{
  if (fIdent!=pgl.GetIdent()) return false; 
// Don't add if our current gc has different (led, num points etc.) from pgl. 


  const Float_t *height = pgl.GetHeight();
  const Float_t *width = pgl.GetWidth();
  const Float_t *numpulses = pgl.GetNumPulses();
  const Float_t *numtriggers = pgl.GetNumTriggers();


  // Allow Add() to add to either end of the curve.
  if (pgl.GetFirstPoint()==1+fLast) {
    // add to end
    fLast = pgl.GetLastPoint();
    fTimeEnd = pgl.GetTimeEnd();
  }
  else if(pgl.GetLastPoint()==fFirst-1) {
    // add to start
    fFirst = pgl.GetFirstPoint();
    fTimeStart = pgl.GetTimeStart();
  }
  else {
    // Fragments don't join up
    return false;
  }

  for (int i = pgl.GetFirstPoint();i<=pgl.GetLastPoint();i++) {
    fHeight[i-1] = height[i-1];
    fWidth[i-1] = width[i-1];
    fNum[i-1] = numpulses[i-1];
    fTriggers[i-1] = numtriggers[i-1];
    //    ++fNumPoints;
  }

  return true;
}

bool PulserGainLog::Add(const PulserSummaryList* psl)
{


    MSG("Pulser",Msg::kDebug)<<"In GainLog Add(psl). Aggregate Number: " << fAggregateNo << endl; 
    MSG("Pulser",Msg::kDebug)<<"In GainLog Add(psl). fFirst, fLast, fIdent: " << fFirst << ", " << fLast << ", " << fIdent << endl;
    MSG("Pulser",Msg::kDebug)<<"In GainLog Add(psl). Trying to add point: " << psl->GetPointIndex()  << ", " << psl->GetIdent() << endl;

  if (fAggregateNo != 0) {
    // We already have some data
    if (fAggregateNo != psl->GetAggregateNo()) return false;
    if (fIdent != psl->GetIdent()) return false;
    if (fLast != (psl->GetPointIndex()-1)) return false;
  } 
  else {
    fAggregateNo = psl->GetAggregateNo();
    fIdent = psl->GetIdent();
    fFirst = psl->GetPointIndex();
    fTimeStart = psl->GetVldContext().GetTimeStamp();
    fNumPoints = psl->GetExpectedPoints();
  }
  fLast = psl->GetPointIndex();
  fTimeEnd = psl->GetVldContext().GetTimeStamp();
  fHeight[fLast-1] = psl->GetHeight();
  fWidth[fLast-1] = psl->GetWidth();
  fNum[fLast-1] = psl->GetNumPulses();
  fTriggers[fLast-1] = psl->GetTriggers();

    MSG("Pulser",Msg::kDebug)<<"In GainLog Add(psl). Aggregate Number: " << fAggregateNo << endl; 
    MSG("Pulser",Msg::kDebug)<<"In GainLog Add(psl). fFirst, fLast, fIdent: " << fFirst << ", " << fLast << ", " << fIdent << endl;
    MSG("Pulser",Msg::kDebug)<<"In GainLog Add(psl). Trying to add point: " << psl->GetPointIndex()  << ", " << psl->GetIdent() << endl;
 
    //  ++fNumPoints;
  return true;
}

---