PEGainCalibrator.cxx

Go to the documentation of this file.

// $Id: PEGainCalibrator.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/PEGainCalibrator.h"
#include "Plex/PlexSEIdAltL.h"
#include "Plex/PlexStripEndId.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include "Validity/VldRange.h"

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

PEGainCalibrator::PEGainCalibrator(int calMode)
   : fResPtr(), fCalMode(calMode)
{
   MSG("Calibrator",Msg::kVerbose) << "PEGainCalibrator::PEGainCalibrator" 
                                   << endl;
}
PEGainCalibrator::PEGainCalibrator(const VldContext vc, int calMode)
   : fResPtr(vc),fCalMode(calMode)
{
   //assert(fResPtr.GetNumRows()!=0);
   //MSG("PE Calibrator",Msg::kVerbose) 
   // << "Number of valid rows in PE gain dB table is "
   // << fResPtr.GetNumRows()  <<  ".\n";
  if(fCalMode>0){
    if(fResPtr.GetNumRows()==0) {
      static int msglimit = 20;  // no more than 20 buggered messages
      if (msglimit) {
        MSG("PE Calibrator",Msg::kWarning) 
          << "No rows in database table with validity context "
          << vc.AsString() << " will assert() here" << ".\n";
        MSG("PE Calibrator",Msg::kWarning)
          << "TimeStamp " << vc.GetTimeStamp() << "\n";
        MSG("PE Calibrator",Msg::kWarning)
          << "SimFlag " << (int)vc.GetSimFlag() << "\n";
        MSG("PE Calibrator",Msg::kWarning)
          << "Detector " << (int)vc.GetDetector() << "\n";
        
        if ( ! ((int)vc.GetDetector() & 127) ) {
          MSG("PE Calibrator",Msg::kWarning)
            << "Detector " << (int)vc.GetDetector() 
            << " is the problem.\n";
        }
        if ( ! ((int)vc.GetSimFlag()  & 127) ) {
          MSG("PE Calibrator",Msg::kWarning)
            << "Simflag " << (int)vc.GetSimFlag() 
            << " is the problem.\n";
        }
        
        const DbiValidityRec* fred = fResPtr.GetValidityRec();
        const VldRange freda=fred->GetVldRange();
        MSG("PE Calibrator",Msg::kWarning) 
          << "bugger "
          << freda.AsString() << ".\n";
        
        if (--msglimit == 0) 
          MSG("PE Calibrator",Msg::kWarning)
            << " ... last warning of this type" << endl;
      }
    }
  }
  else{
          MSG("PE Calibrator",Msg::kInfo) << " No PE table loaded - defaults used.\n";

  }
}

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

Float_t PEGainCalibrator::Calibrate(Float_t rawcharge,
                                         PlexStripEndId seid) const
{
//======================================================================
// Purpose: Calibrate method
//======================================================================
   if (fCalMode) return ApplyCalib(rawcharge,seid);
   else          return rawcharge/60.0;   
}

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

Float_t  PEGainCalibrator::ApplyCalib(Float_t rawcharge, 
                                           PlexStripEndId seid) const
{
//======================================================================
// Purpose: The actual numerical calibration, taking Raw Charge and
// converting to PEs.
//======================================================================

   Int_t psekey = seid.BuildPlnStripEndKey();
   //MSG("Calibrator",Msg::kVerbose) 
   //   << "PEGainCalibrator on " << rawcharge << " raw charge in seid " 
   //   << seid << " (indexed as " << psekey <<" )\n";

   // Now need to get the row which corresponds to the stripendnum.
   const CalADCtoPE* dpgc = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
   if(dpgc ==0) {
     static int msglimit = 20;  // no more than 20 missing row messages
     if (msglimit) {
       MSG("Calibrator",Msg::kWarning) 
         << "No database row for StripEnd " << seid << " (indexed as " 
         << psekey << " )\n";
       if (--msglimit == 0)
         MSG("Calibrator",Msg::kWarning)
           << " ... last warning of this type" << endl;
     }
     return rawcharge;
   }

 return dpgc->GetPE(rawcharge);
}

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



void PEGainCalibrator::ReInit(VldContext vc) {

  // if calibration is off, just return
  if (fCalMode == 0) return;

  // MSG("Calibrator",Msg::kVerbose) 
  //    << "PEGainCalibrator::ReInit(VldContext)" << endl;

  fResPtr.NewQuery(vc,0);

  MSG("PE Calibrator",Msg::kVerbose) 
    << "PEGainCalibrator::ReInit(VldContext)  mode " << fCalMode << endl
    << "Number of valid rows in PE gain dB table is "
    << fResPtr.GetNumRows()  <<  ".\n";

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


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

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

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

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

---