TimeCalScheme.cxx

Go to the documentation of this file.

// $Id: TimeCalScheme.cxx,v 1.7 2005/03/22 12:53:04 tagg Exp $
//
// Prototype PE gain Calibrator class
//
// rjn@hep.ucl.ac.uk
#include "TimeCalScheme.h"

#include "MessageService/MsgService.h"
#include "Plex/PlexSEIdAltL.h"
#include "Plex/PlexStripEndId.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include "Validity/VldRange.h"
#include "Conventions/Munits.h"

CVSID("$Id: TimeCalScheme.cxx,v 1.7 2005/03/22 12:53:04 tagg Exp $");
ClassImp(TimeCalScheme)

TimeCalScheme::TimeCalScheme()
{
  Registry r;
  // No configurable options.
  InitializeConfig(r); 
}


//.....................................................................
void TimeCalScheme::DoReset( const VldContext& context )
{
  MSG("Calib",Msg::kVerbose) << "TimeCalScheme::DoReset()" << endl;
  
  fResPtr.NewQuery(context,0);
  
  if(fResPtr.GetNumRows()==0) {
    MAXMSG("Calib",Msg::kWarning,10) 
      << "No rows in database table with validity context "
      << context.AsString() << "  No calibration will be applied." << endl;
    IncrementErrors(kTimeCalibrator,kMissingTable);
  }  
}

//.....................................................................
void TimeCalScheme::ConfigModified()
{                                               
  // No configurable options.
}

//.....................................................................
void TimeCalScheme::PrintConfig( std::ostream& /*os*/ ) const
{
  // Nothing to print.
}

 
//.....................................................................
DoubleErr TimeCalScheme::GetCalibratedTime(DoubleErr rawtime,  
                                          FloatErr,
                                          const PlexStripEndId& seid) const
{

  // Now need to get the row which corresponds to the stripendnum.
  const CalTimeCalibration* dpgc = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());

  if(dpgc==0) {
    if(fResPtr.GetNumRows()>0) {
      MAXMSG("Calib",Msg::kWarning,10) 
        << "TimeCalScheme: No database row for StripEnd " << seid.AsString() << "\n";
      
      IncrementErrors(kTimeCalibrator,kMissingRow,seid);
    }
    return rawtime + FloatErr(0,50*Munits::ns); // Add error of 50ns, since we really don't know.
  }
  return dpgc->GetTime(rawtime);
}


//.....................................................................
DoubleErr TimeCalScheme::DecalTime(DoubleErr caltime,    
                                  FloatErr,
                                  const PlexStripEndId& seid) const
{

  // Now need to get the row which corresponds to the stripendnum.
  const CalTimeCalibration* dpgc = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());

  if(dpgc==0) {
    // Print warning message / log errors if-and-only-if the table is not empty:
    if(fResPtr.GetNumRows()>0) { 
      MAXMSG("Calib",Msg::kWarning,10) 
        << "TimeCalScheme: No database row for StripEnd " << seid.BuildPlnStripEndKey() << "\n";
      
      IncrementErrors(kTimeCalibrator,kMissingRow,seid);
    }

    return caltime + FloatErr(0,50*Munits::ns); // Add error of 50ns, since we really don't know.;
  }

  // No inverse function exists in the row function, so I'll kludge one here.
  // NJT 7/04
  return (caltime +  (double)dpgc->GetOffset())*(double)dpgc->GetScale();
}

---