SigMappedCalibrator.cxx

Go to the documentation of this file.

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

CVSID("$Id: SigMappedCalibrator.cxx,v 1.1 2004/07/26 16:52:20 tagg Exp $");
SigMappedCalibrator::SigMappedCalibrator(Int_t calMode)
   : fMuonPtr(), fMapperFitsPtr(), fMapperDataPtr(), fCalMode(calMode)
{
   MSG("Calibrator",Msg::kVerbose) << "SigMappedCalibrator::SigMappedCalibrator" 
                                   << endl;
}

SigMappedCalibrator::SigMappedCalibrator(const VldContext vc, Int_t calMode)
   : fMuonPtr(vc),fMapperFitsPtr(vc),fMapperDataPtr(vc),fCalMode(calMode)
{
   //assert(fMuonPtr.GetNumRows()!=0);

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

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

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

Float_t  SigMappedCalibrator::ApplyCalib(Float_t siglin, Float_t Y,
                                           PlexStripEndId seid) const
{
//======================================================================
// Purpose: The actual numerical calibration.
  // depending upon mode variable, uses mapper fit, muon fit, or
  // 'raw' mapper data
//======================================================================

  Float_t MapCorrected;
  Float_t MuonCorrected;
  Float_t MuonPlusMapCorrected;
   
 MSG("Calibrator",Msg::kVerbose) << "SigMappedCalibrator on " << siglin 
                                 << " siglin in seid " << seid 
                                 << " (indexed as " << seid.BuildPlnStripEndKey() 
                                 <<" )\n";
 
 switch (fCalMode){
 case 0:
    // in this case forgo calibration
    
    break;
 case 1:
    {
    // default:  Use fits to mapper only
       const CalMapperFits* mappercal = fMapperFitsPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
       if(mappercal ==0) {
          static int msglimit = 20;  // no more than 20 complaints
          if (msglimit) {
            MSG("MapperCalibrator",Msg::kWarning) 
              << "No database row for StripEnd " << seid 
              << " (indexed as " << seid.BuildPlnStripEndKey() << ")" << endl;
            if (--msglimit == 0)
              MSG("MapperCalibrator",Msg::kWarning)
                << " ... last warning of this type" << endl;
          }
          return siglin;
       }
       
       MapCorrected = mappercal->AttenCorrected(siglin,Y);
       return MapCorrected;
       
       break;
    }
 case 2:
    {
       // use full mapper table only
       
       break;
    }
 case 3:
    {
       // use muon fits 
       const CalMuonCalibration* muoncal = fMuonPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
       if(muoncal ==0) {
       MSG("MuonCalibrator",Msg::kWarning) 
          << "No database row for StripEnd " << seid.BuildPlnStripEndKey() << "\n";
       return siglin;
       }
       MuonCorrected = muoncal->AttenCorrected(siglin,Y);
       return MuonCorrected;
       break;
    }
 case 4:
    {
       // use muon fits, then mapper table as tweaks
       // use muon fits 
       const CalMuonCalibration* muoncal = fMuonPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
       const CalMapperData* mapperdata = fMapperDataPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
       if(muoncal ==0) {
          MSG("MuonCalibrator",Msg::kWarning) 
             << "No database row for StripEnd " << seid.BuildPlnStripEndKey() << "\n";
          return siglin;
       }
       
       MuonCorrected = muoncal->AttenCorrected(siglin,Y); 
       MuonPlusMapCorrected = mapperdata->ResponseCorrected(MuonCorrected,Y);
       return MuonPlusMapCorrected;
       break;
    }

 }


 return siglin; 

}

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

void SigMappedCalibrator::ReInit(VldContext vc) {
   MSG("Calibrator",Msg::kVerbose) << "SigMappedCalibrator::ReInit(VldContext)" << endl;
   switch (fCalMode){
   case 0:
      // in this case forgo calibration
      break;
   case 1:
      // default:  Use fits to mapper only
      fMapperFitsPtr.NewQuery(vc,0);
      break;
   case 2:
      // use full mapper table only
      break;
   case 3:
      // use muon fits 
      fMuonPtr.NewQuery(vc,0);
      break;
   case 4:
      // use muon fits, then mapper table as tweaks
      fMuonPtr.NewQuery(vc,0);
      fMapperDataPtr.NewQuery(vc,0);
      break;
   default:
      //Don't know what to do.
      break;
   }
}

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

const CalMapperFits* SigMappedCalibrator::GetCalMapperFits(PlexStripEndId seid) const
{
 return fMapperFitsPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
}

---