MIPCalibrator.cxx

Go to the documentation of this file.

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

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


MIPCalibrator::MIPCalibrator(Int_t calMode)
   : fResPtr(),fCalMode(calMode)
{
 MSG("MIPCalibrator",Msg::kVerbose) << "MIPCalibrator::MIPCalibrator" 
                                   << endl;

}
MIPCalibrator::MIPCalibrator(VldContext vc, Int_t calMode)
   : fResPtr(vc),fCalMode(calMode)
{
 //assert(fResPtr.GetNumRows()!=0);

  if(calMode>0){
    if(fResPtr.GetNumRows()==0) {
      MSG("Calibrator",Msg::kWarning) 
        << "No rows in database table with validity context "
        << vc.AsString() << " will assert() here" << ".\n";
      MSG("Calibrator",Msg::kWarning)
        << "TimeStamp " << vc.GetTimeStamp() << "\n";
      MSG("Calibrator",Msg::kWarning)
        << "SimFlag " << (int)vc.GetSimFlag() << "\n";
      MSG("Calibrator",Msg::kWarning)
        << "Detector " << (int)vc.GetDetector() << "\n";
      
      
      if ( ! ((int)vc.GetDetector() & 127) ) {
        MSG("Calibrator",Msg::kWarning)
          << "Detector " << (int)vc.GetDetector() 
          << " is the problem.\n";
      }
      if ( ! ((int)vc.GetSimFlag()  & 127) ) {
        MSG("Calibrator",Msg::kWarning)
          << "Simflag " << (int)vc.GetSimFlag() 
          << " is the problem.\n";
      }
      
      
      
      const DbiValidityRec* fred = fResPtr.GetValidityRec();
      const VldRange freda=fred->GetVldRange();
      MSG("Calibrator",Msg::kWarning) 
        << "bugger "
        << freda.AsString() << ".\n";
    }
  }
  else{
      MSG("Calibrator",Msg::kInfo) << "No MIP table loaded - default used.\n";

  }

}

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

Float_t MIPCalibrator::Calibrate(Float_t siglin,
                                         PlexStripEndId seid) const
{
//======================================================================
// Purpose: Calibrate method
//======================================================================
   return ApplyCalib(siglin,seid);
   
}

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

Float_t  MIPCalibrator::ApplyCalib(Float_t siglin, 
                                           PlexStripEndId seid) const
{
//======================================================================
// Purpose: The actual numerical calibration, taking Raw Charge and
// converting to Pulser Units.
//======================================================================
 MSG("Calibrator",Msg::kVerbose) 
      << "MIPCalibrator on " << siglin << " siglin in seid " 
      << seid << " (indexed as " << seid.BuildPlnStripEndKey() <<" )\n";  
  // Now need to get the row which corresponds to the stripendnum.
 const CalMIPCalibration* mipcal = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
 if(mipcal ==0) {
    static int msglimit = 20;  // no more than 20 complaints
    if (msglimit) {
      MSG("Calibrator",Msg::kWarning) 
        << "No database row for StripEnd " << seid << " (indexed as " 
        << seid.BuildPlnStripEndKey() << " )\n";
      if (--msglimit == 0) 
        MSG("Calibrator",Msg::kWarning)
          << " .. last warning of this type" << endl;
    }
    // in case where db lookup fails and warning given, apply default
    // calibration (1MIP = 639.6 ADU)
    return siglin/639.6;
 }
 Float_t MIP = mipcal->GetMIP(siglin);
 
 return MIP; 
}

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

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

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


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



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

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

}

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

---