TimeCalibrator.cxx

Go to the documentation of this file.

// $Id: TimeCalibrator.cxx,v 1.1 2004/07/26 16:52:20 tagg Exp $
//
// Prototype PE gain Calibrator class
//
// rjn@hep.ucl.ac.uk
#include "MessageService/MsgService.h"
#include "CandDigit/CandDigitHandle.h"
#include "Calibrator/TimeCalibrator.h"
#include "Plex/PlexSEIdAltL.h"
#include "Plex/PlexStripEndId.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include "Validity/VldRange.h"

CVSID("$Id: TimeCalibrator.cxx,v 1.1 2004/07/26 16:52:20 tagg Exp $");

TimeCalibrator::TimeCalibrator(int calMode)
   :fResPtr(), fCalMode(calMode)
{
   MSG("Calibrator",Msg::kVerbose) << "TimeCalibrator::TimeCalibrator" 
                                   << endl;
}

TimeCalibrator::TimeCalibrator(const VldContext vc, int calMode)
   : fResPtr(vc),fCalMode(calMode)
{
   //assert(fResPtr.GetNumRows()!=0);

  if(fCalMode>0){
    if(fResPtr.GetNumRows()==0) {
        MSG("Time Calibrator",Msg::kWarning) 
          << "No rows in database table with validity context "
          << vc.AsString() << " will assert() here" << ".\n";
        MSG("Time Calibrator",Msg::kWarning)
          << "TimeStamp " << vc.GetTimeStamp() << "\n";
        MSG("Time Calibrator",Msg::kWarning)
          << "SimFlag " << (int)vc.GetSimFlag() << "\n";
        MSG("Time Calibrator",Msg::kWarning)
          << "Detector " << (int)vc.GetDetector() << "\n";
        
        
        if ( ! ((int)vc.GetDetector() & 127) ) {
          MSG("Time Calibrator",Msg::kWarning)
            << "Detector " << (int)vc.GetDetector() 
            << " is the problem.\n";
        }
        if ( ! ((int)vc.GetSimFlag()  & 127) ) {
          MSG("Time Calibrator",Msg::kWarning)
            << "Simflag " << (int)vc.GetSimFlag() 
            << " is the problem.\n";
        }
        
        
        
        const DbiValidityRec* fred = fResPtr.GetValidityRec();
        const VldRange freda=fred->GetVldRange();
        MSG("Time Calibrator",Msg::kWarning) 
          << "bugger "
          << freda.AsString() << ".\n";
    }
  }
  else{
    static int msglimit = 20; // no more than 20 complaints
    if (msglimit) {
        MSG("Time Calibrator",Msg::kInfo   ) << "No Time Cal Table Loaded.\n";
        if (--msglimit == 0) 
          MSG("Time Calibrator",Msg::kInfo   ) 
            << "... last warning of this type." << endl;
    }

  }
}

//......................................................................

Double_t TimeCalibrator::Calibrate(Double_t rawtime,
                                         PlexStripEndId seid) const
{
//======================================================================
// Purpose: Calibrate method
//======================================================================
   return ApplyCalib(rawtime,seid);
   
}

//......................................................................

Double_t  TimeCalibrator::ApplyCalib(Double_t rawtime, 
                                           PlexStripEndId seid) const
{
//======================================================================
// Purpose: The actual numerical calibration, taking Raw Charge and
// converting to Pulser Units.
//======================================================================
   
 MSG("Calibrator",Msg::kVerbose) << "TimeCalibrator on " << rawtime 
                                 << " seid " << seid.BuildPlnStripEndKey() <<"\n";

 // Now need to get the row which corresponds to the stripendnum.
 const CalTimeCalibration* dpgc = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
 if(dpgc ==0) {
   //    MSG("Calibrator",Msg::kWarning) 
   //    << "No database row for StripEnd " << seid.BuildPlnStripEndKey() << "\n";
    return rawtime;
 }
 return dpgc->GetTime(rawtime);
}


//......................................................................

const CalTimeCalibration* TimeCalibrator::GetCalTimeCalibration(PlexStripEndId seid) const
{
  return fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
}

//......................................................................


void TimeCalibrator::ReInit(VldContext vc) {
   MSG("Calibrator",Msg::kVerbose) << "TimeCalibrator::ReInit(VldContext)" << endl;

   fResPtr.NewQuery(vc,0);

   //assert(fResPtr.GetNumRows()!=0);
   if(fResPtr.GetNumRows()==0) {
     static int msglimit = 20;  // no more than 20 complaints
     if (msglimit) {
        MSG("Time Calibrator",Msg::kWarning) 
           << "No rows in database table with validity context "
           << vc.AsString() << " will assert() here" << ".\n";
        MSG("Time Calibrator",Msg::kWarning)
           << "TimeStamp " << vc.GetTimeStamp() << "\n";
        MSG("Time Calibrator",Msg::kWarning)
           << "SimFlag " << (int)vc.GetSimFlag() << "\n";
        MSG("Time Calibrator",Msg::kWarning)
           << "Detector " << (int)vc.GetDetector() << "\n";


        if ( ! ((int)vc.GetDetector() & 127) ) {
           MSG("Time Calibrator",Msg::kWarning)
           << "Detector " << (int)vc.GetDetector() 
           << " is the problem.\n";
        }
        if ( ! ((int)vc.GetSimFlag()  & 127) ) {
           MSG("Time Calibrator",Msg::kWarning)
           << "Simflag " << (int)vc.GetSimFlag() 
           << " is the problem.\n";
        }



        const DbiValidityRec* fred = fResPtr.GetValidityRec();
        const VldRange freda=fred->GetVldRange();
        MSG("Time Calibrator",Msg::kWarning) 
           << "bugger "
           << freda.AsString() << ".\n";

        if (--msglimit == 0)
          MSG("Calibrator",Msg::kWarning)
            << " ... last warning of this type" << endl;
     }
   }
}

//......................................................................

---