VaLinearityCalScheme.cxx

Go to the documentation of this file.

// $Id: VaLinearityCalScheme.cxx,v 1.4 2005/08/19 03:03:58 tagg Exp $
//
// Scheme to linearize PIN channels on VA (far detector)
// Interface to G. Barr's CalVaLinearity table.
//
// Can in principle be applied to any VA channel, but usually
// only PIN-diodes are linearized this way.
//
//
// n.tagg1@physics.ox.ac.uk
#include "VaLinearityCalScheme.h"

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

ClassImp(VaLinearityCalScheme)
CVSID("$Id: VaLinearityCalScheme.cxx,v 1.4 2005/08/19 03:03:58 tagg Exp $");

//......................................................................
VaLinearityCalScheme::VaLinearityCalScheme()
{
  MSG("Calib",Msg::kVerbose) << "VaLinearityCalScheme::VaLinearityCalScheme" 
                             << endl;
  Registry r;
  r.Set("Task","PinsOnly");
  InitializeConfig(r);
}

//......................................................................
void VaLinearityCalScheme::DoReset( const VldContext& vc )
{

  MSG("Calib",Msg::kVerbose) << "VaLinearityCalScheme::DoReset()" << endl;

  fResPtr.NewQuery(vc,fTask);
  if(fResPtr.GetNumRows()==0) {
    MAXMSG("Calib",Msg::kWarning,10) 
      << "VaLinearity Scheme: No rows in CALVALINEARITY database table"
      << vc.AsString() << ".\n";         
    IncrementErrors(kVALinCalibrator,kMissingTable);
  }
}

//......................................................................
void VaLinearityCalScheme::ConfigModified()
{
  bool ok = true;
  const char* taskname;
  ok = ok && GetConfig().Get("Task",taskname);

  if     (0==strcasecmp(taskname,"PinsOnly" )) fTask = 0;
  else if(0==strcasecmp(taskname,"All"))       fTask = 1;
  else ok = false;

  if(!ok) MSG("Calib",Msg::kError) << "Problem setting up VaLinearityCalScheme config." << endl;
  
  // Ensure that the DB has been changed for this event.
  Reset(GetContext(),true); 
}


//......................................................................
void VaLinearityCalScheme::PrintConfig( std::ostream& os ) const
{
  os << " VaLinearityCalScheme " << endl;
  os << " Task: " << GetConfig().GetCharString("Task");
}




//......................................................................
FloatErr  VaLinearityCalScheme::GetLinearizedVA(FloatErr rawcharge, 
                                                const RawChannelId& rcid) const
{

  // Reject stupid inputs:
  if(rcid.IsNull()) return rawcharge;
  if(rcid.GetElecType()!=ElecType::kVA) return rawcharge;

  // Get Row ID:
  UInt_t idx = CalVaLinearity::Rcid2Index(rcid);
  
  const CalVaLinearity* row = fResPtr.GetRowByIndex(idx);

  if(row==0) {
    if(fResPtr.GetNumRows()>0) {
      MAXMSG("Calib",Msg::kWarning,10) 
        << "VaLinearityCalScheme: No database row for RawChannelId" << rcid.AsString() 
        << "  Idx:" << Form("%08x",idx) << std::endl;      
      IncrementErrors(kVALinCalibrator,kMissingRow,rcid);
    }
    return rawcharge * FloatErr(1,0.5);
  }
  
  return row->Linearize(rawcharge);
}


//......................................................................
FloatErr VaLinearityCalScheme::DecalVALinearity(FloatErr lincharge, 
                                                const RawChannelId& rcid) const
{

  // Reject stupid inputs:
  if(rcid.IsNull()) return lincharge;
  if(rcid.GetElecType()!=ElecType::kVA) return lincharge;
  

  // Get Row ID:
  UInt_t idx = CalVaLinearity::Rcid2Index(rcid);
  
  const CalVaLinearity* row = fResPtr.GetRowByIndex(idx);

  if(row==0) {
    if(fResPtr.GetNumRows()>0) {
      MAXMSG("Calib",Msg::kWarning,10) 
        << "VaLinearityCalScheme: No database row for RawChannelId" << rcid.AsString() 
        << "  Idx:" << Form("%08x",idx) << std::endl;      
      IncrementErrors(kVALinCalibrator,kMissingRow,rcid);
    }
    return lincharge * FloatErr(1,0.5);
  }
  
  return row->UnLinearize(lincharge);
}

---