PulserDriftCalScheme.cxx

Go to the documentation of this file.

// $Id: PulserDriftCalScheme.cxx,v 1.18 2006/06/15 01:57:48 rhatcher Exp $
//
// Calibrator for new PulserCalibration drift correction
//
// pa@hep.ucl.ac.uk
//
// $Log: PulserDriftCalScheme.cxx,v $
// Revision 1.18  2006/06/15 01:57:48  rhatcher
// Initialize all elements of the "Drift" struct after creating it, not all
// paths through the code set each element.
//
// Revision 1.17  2006/04/19 18:50:38  rhatcher
// Change DetectorType:: to Detector:: everywhere.
//
// Revision 1.16  2006/01/26 10:28:13  cpw1
// Fix for new CALPULSERFITS
//
// Revision 1.15  2006/01/25 17:31:45  rhatcher
// In FitIsOK() comment out line using CalPulserFit::GetNumPoints() because
// that method not longer exists for CalPulserFit.
//
// Revision 1.14  2006/01/25 13:48:33  cpw1
// Use zero-corrected values and implement checks on fits
//
// Revision 1.13  2005/09/28 09:30:23  cpw1
// Add option to use CalPulserFits method
//
// Revision 1.12  2005/05/13 18:07:14  tagg
// Make a few warnings MAXMSG from MSG.
//
// Revision 1.11  2005/04/09 10:33:42  hartnell
//
// Get rid of "Only 343 entries for..." messages here too.
//
// Revision 1.10  2005/03/22 12:37:04  tagg
// Set up to use errors system.
//
// Applies a 50% error if no drift data can be found.
// Does NOT apply any error associated with drift point; I'm not yet
// confident enough in Phil's code to figure this out.
//
// Tidy up code so to remove redundant call level.
//
// Revision 1.9  2005/02/01 17:33:17  tagg
// Very cosmetic changes to code to get working with gcc 3.4. Mostly
// extra top-level semicolons (what a drag!) and getting CVSID in the
// right place (i.e. before #including DbiResultPtr.tpl)
//
// Revision 1.8  2005/01/28 20:21:00  cbs
// Removed the 2 hour time lag in GetTemperature() for CalDet.
// Fixed a bug in the CalDet2003Check() functions in Pulser*CalScheme.cxx
// Commented out the CalDet2003Check() and the check that the channel is of
// type VA before applying the temperature calibration in the
// Pulser*CalScheme::ApplyCalib() and Pulser*CalScheme::InverseCalib().
// The VA check needs to go back in if temperature calibrations are ever
// switched on for Near, Far detectors.
//
// Revision 1.7  2004/12/07 15:38:35  cbs
// VA temperature correction now gets applied to QIE hits during
// CalDet 2003 N/F running
//
// Revision 1.6  2004/12/06 18:17:07  cbs
// Fixed a small mistake in the application of the temperature correction
//
// Revision 1.5  2004/12/06 15:27:25  cbs
// Added code to Drift/SigLin CalSchemes to perform a temperature correction
// to the siglin values. A correction factor of 0.22% / degree C is applied
// to all hits from VA channels. The reference temperature is 18 degC.
//
// Modified the existing temperature correction in StripToStripCalScheme.
// A factor of 0.08% / degree C is applied to all hits from all channels to
// correct for residual measured changes in response at CalDet after the
// VA dependence had been removed.
//
// The code uses GetTemperature() to get the current temperature.
//
// Revision 1.4  2004/11/15 18:52:19  phil.adamson
// If we can't calibrate a strip, try a nearby one.
//
// Revision 1.3  2004/11/04 23:21:12  phil.adamson
// Fix stupid error
//
// Revision 1.2  2004/11/04 22:17:37  phil.adamson
// Minor fix for new calib
//
// Revision 1.1  2004/11/04 04:46:53  phil.adamson
// New for new LI calib
//
#include "Calibrator/PulserDriftCalScheme.h"
#include "Calibrator.h"
#include "MessageService/MsgService.h"
#include "Plex/PlexSEIdAltL.h"
#include "Plex/PlexStripEndId.h"
#include "Plex/PlexPinDiodeId.h"
#include "Plex/PlexHandle.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include "Validity/VldRange.h"
#include "Validity/VldTimeStamp.h"
#include "PulserCalibration/PulserConventions.h"
#include "PulserCalibration/PulserPinScale.h"
#include "PulserCalibration/PulserXScale.h"

ClassImp(PulserDriftCalScheme)
CVSID("$Id: PulserDriftCalScheme.cxx,v 1.18 2006/06/15 01:57:48 rhatcher Exp $");


//......................................................................
PulserDriftCalScheme::PulserDriftCalScheme(): fPlex(0)
{
   MSG("Calib",Msg::kVerbose) 
     << "PulserDriftCalScheme::PulserDriftCalScheme" 
     << endl;

   Registry r;
   r.Set("CalMode",0);
   r.Set("MinNumPulses",500);
   r.Set("VATemperatureCalibration","off");
   InitializeConfig(r);

   fRefTemp = 18.0; //reference temperature
   fVATempCorFactor = -0.0022; //frac diff in VA electronics response with temp

}

//......................................................................
void PulserDriftCalScheme::DoReset( const VldContext& vc )
{
   // Reference tables should be for fixed context
   VldTimeStamp ts(2005,6,2,0,0,0);
   VldContext ref(vc.GetDetector(),vc.GetSimFlag(),ts);
   fXScale.NewQuery(ref,Pulser::kNearEnd);
   fFarXScale.NewQuery(ref,Pulser::kFarEnd);
   fPinScale.NewQuery(ref,0);
   fNearLow.NewQuery(ref,Pulser::kNearLow);
   fFarLow.NewQuery(ref,Pulser::kFarLow);
   fNearHigh.NewQuery(ref,Pulser::kNearHigh);
   fFarHigh.NewQuery(ref,Pulser::kFarHigh);
   fNearFar.NewQuery(ref,Pulser::kNearFar);

   // Drift values
   fNearDrift.NewQuery(vc,Pulser::kNearEnd);
   fFarDrift.NewQuery(vc,Pulser::kFarEnd);
   fPinDrift.NewQuery(vc,Pulser::kNearEnd);
   if (fPlex) delete fPlex;
   fPlex = new PlexHandle(vc);
  
   // Don't bother with warning messages, but do increment errors.

  if(fXScale.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fFarXScale.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fNearLow.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fFarLow.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fNearHigh.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fFarHigh.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fNearFar.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fNearDrift.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fFarDrift.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
  if(fPinDrift.GetNumRows()<=0) IncrementErrors(kDriftCalibrator,kMissingTable);
}

//......................................................................
void PulserDriftCalScheme::ConfigModified()
{
  
  const char* str_tempcal;

  bool ok = true;
  ok = ok && GetConfig().Get("CalMode",     fCalMode);
  ok = ok && GetConfig().Get("MinNumPulses",fMinNumPulses);
  ok = ok && GetConfig().Get("VATemperatureCalibration",str_tempcal);
  if(!ok) MSG("Calib",Msg::kWarning) << "PulserDriftCalibrator "
                                          << " Problem when configuring. " <<endl;
  
  // Interpret strings.
  if(strcasecmp(str_tempcal,"on")==0) fDoTempCal = true;
  else fDoTempCal = false;

}
//......................................................................
void PulserDriftCalScheme::PrintConfig( std::ostream& os ) const
{
  os << "  CalMode:      " << fCalMode << endl;
  os << "  MinNumPulses: " << fMinNumPulses << endl;
  os << "  VA Temperature calibration is " << ((fDoTempCal)?"on":"off") << endl;
  if(fDoTempCal) {
    os << "  VA Temp Correction Factor = 1 / [ 1 - " 
       << fVATempCorFactor << "*(" 
       << fRefTemp << " - temp) ]" << endl;
  }
}


//......................................................................
FloatErr  PulserDriftCalScheme::GetDriftCorrected(FloatErr rawcharge, 
                                           const PlexStripEndId& seid) const
{

  //check for VA temperature correction and that channel is VA:
  float tempCor = 1.0;
  if(fPlex && fDoTempCal) {    
    //if(fPlex->GetRawChannelId(seid).GetElecType()==ElecType::kVA ||
    // this->CalDet2003Check())
    tempCor = GetVATemperatureCorrection();
  }

  // Drift correction = "xscale" * driftPIN/driftADC
  // where xscale is our calibration base, and looks like ADC/PIN
  // and comes from the slope of the first gaincurve at the low end.

  Drift d = GetDriftPoint(seid);
  if (d.adc > 0) {
    Float_t value = rawcharge * (float)( d.xscale * d.pin/d.adc * tempCor);
    Float_t error2 = (float) (d.erradc * d.erradc / (d.adc * d.adc));
    if ( d.pin > 0) {error2+= d.errpin*d.errpin/(d.pin*d.pin);}
    if ( d.xscale > 0) {error2 += d.errxscale*d.errxscale/(d.xscale*d.xscale);}
    return FloatErr(value,value*sqrt(error2));
  }
  else {
    // Go through other stripends sharing this pixel spot
    PlexSEIdAltL altlist = fPlex->GetSEIdAltL(fPlex->GetRawChannelId(seid));
    for (PlexSEIdAltL::const_iterator it = altlist.begin();
         it!=altlist.end();it++) {
      Drift d = GetDriftPoint((*it).GetSEId());
      if (d.adc > 0) {
        MAXMSG("Calib",Msg::kWarning,20) 
          << "Using stripend "<<(*it).GetSEId()<<" to calibrate "<<seid<<endl;
      Float_t value = rawcharge * (float)( d.xscale * d.pin/d.adc * tempCor);
      Float_t error2 = (float) (d.erradc * d.erradc / (d.adc * d.adc));
      if ( d.pin > 0) {error2+= d.errpin*d.errpin/(d.pin*d.pin);}
      if ( d.xscale > 0) {error2 += d.errxscale*d.errxscale/(d.xscale*d.xscale);}
      return FloatErr(value,value*sqrt(error2));
      }
    }
  }

  // Try other strips on same plane for Anatael.

  int numstrips = seid.NumberOfStrips();
  PlexStripEndId otherstrip = seid;
  for (int n = 0; n< numstrips; n++) {
    otherstrip.SetStrip(n);
    Drift d = GetDriftPoint(otherstrip);
    if (d.adc > 0) {
      Float_t value = rawcharge * (float)( d.xscale * d.pin/d.adc * tempCor);
      Float_t error2 = (float) (d.erradc * d.erradc / (d.adc * d.adc));
      if ( d.pin > 0) {error2+= d.errpin*d.errpin/(d.pin*d.pin);}
      if ( d.xscale > 0) {error2 += d.errxscale*d.errxscale/(d.xscale*d.xscale);}
      return FloatErr(value,value*sqrt(error2));
    }
  }

  MAXMSG("Calib",Msg::kWarning,20) 
    << "No drift point for  StripEnd " << seid.AsString()
    << " so NOT CALIBRATING THIS CHANNEL\n";
  
  IncrementErrors(kDriftCalibrator,kMissingRow,seid);
  FloatErr result = rawcharge * tempCor;
  result*=FloatErr(1,0.5); // Add another 50% error.
  return result;
}   


FloatErr PulserDriftCalScheme::DecalDrift(FloatErr undrifted, 
                                          const PlexStripEndId& seid) const
{

  //check for temperature correction:
  float tempCor = 1.0;
  if(fPlex && fDoTempCal) {    
    //if(fPlex->GetRawChannelId(seid).GetElecType()==ElecType::kVA ||
    // this->CalDet2003Check())       
    tempCor = GetVATemperatureCorrection();
  }

  // Drift correction = "xscale" * driftPIN/driftADC
  // where xscale is our calibration base, and looks like ADC/PIN
  // and comes from the slope of the first gaincurve at the low end.
  
  Drift d = GetDriftPoint(seid);
  if (d.adc > 0) {
    return undrifted * (float)((1.0/d.xscale) * (d.adc/d.pin) / tempCor); // FIXME: No error returned
  }
  else {
    // Go through other stripends sharing this pixel spot
    PlexSEIdAltL altlist = fPlex->GetSEIdAltL(fPlex->GetRawChannelId(seid));
    for (PlexSEIdAltL::const_iterator it = altlist.begin();
         it!=altlist.end();it++) {
      Drift d = GetDriftPoint((*it).GetSEId());
      if (d.adc > 0) {
        MAXMSG("Calib",Msg::kWarning,20) 
          << "Using stripend "<<(*it).GetSEId()<<" to uncalibrate "<<seid<<endl;
        return undrifted * (float)((1.0/d.xscale) * (d.adc/d.pin) / tempCor); // FIXME: No error returned
      }
      
    }
  }
  
  MAXMSG("Calib",Msg::kWarning,20) 
    << "No drift point for  StripEnd " << seid.AsString()
    << " so NOT UNCALIBRATING THIS CHANNEL\n";
  
  IncrementErrors(kDriftCalibrator,kMissingRow,seid);
  FloatErr result = undrifted / tempCor;
  result *= FloatErr(1,0.5); // Add 50% error
  return result;
}   



FloatErr PulserDriftCalScheme::GetDriftPinDiodeValue(PlexLedId ledid, Int_t highOrLow) const
{

   PlexHandle ph(GetContext());
   std::pair<PlexPinDiodeId,PlexPinDiodeId> myPins = 
      ph.GetPinDiodeIds(ledid); // NB first = low, second = high
   
   if(highOrLow==0) {//high ==0  low ==1
     const PulserDriftPin *highpin =fPinDrift.GetRowByIndex(myPins.second.GetEncoded());
     if(highpin==0) {
       MSG("Calib",Msg::kWarning) 
          << "No Drift Point database row for Pin " << myPins.second 
          << " (indexed as " << myPins.second.GetEncoded() << " there are "
          << fPinDrift.GetNumRows() << " rows\n.";
       IncrementErrors(kDriftCalibrator,kMissingRow,ledid);
       return -1;
     }  
     else {
       MSG("Calib",Msg::kDebug) 
          << "Got Drift Point database row for Pin " << myPins.second 
          << " (indexed as " << myPins.second.GetEncoded() << " )\n";
       if(highpin->GetNumEntries()<fMinNumPulses) {
        MSG("Calib",Msg::kWarning) 
           << "Only " << highpin->GetNumEntries() << " entries for  " 
           << myPins.second 
           << " (indexed as " << myPins.first.GetEncoded() << " )" << endl;
        return -1;
       }
       Float_t value = highpin->ZCMean();
       Float_t error = highpin->ZCError();
       return FloatErr(value,error);
       //       return FloatErr(highpin->ZCMean(),highpin->ZCError());
     }
   }
   else if(highOrLow==1) { //low gain pin 
     const PulserDriftPin *lowpin =fPinDrift.GetRowByIndex(myPins.first.GetEncoded());
     if(lowpin==0) {
       MSG("Calib",Msg::kWarning) 
          << "No Drift Point database row for Pin " << myPins.first
          << " (indexed as " << myPins.first.GetEncoded() << " )\n";
       return -1;
     }  
     else {
       MSG("Calib",Msg::kDebug) 
          << "Got Drift Point database row for Pin " << myPins.first 
          << " (indexed as " << myPins.first.GetEncoded() << " )\n";
       if(lowpin->GetNumEntries()<fMinNumPulses) {
         MSG("Calib",Msg::kWarning) 
            << "Only " << lowpin->GetNumEntries() << " entries for  " 
            << myPins.first
            << " (indexed as " << myPins.first.GetEncoded() << " )" << endl;
         return -1;
       }
       return FloatErr(lowpin->ZCMean(),lowpin->ZCError());
     }
   }           
   
   return -1;
}

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




PulserDriftCalScheme::Drift 
  PulserDriftCalScheme::GetDriftPoint(PlexStripEndId seid) const
{
  Drift drift;
  drift.adc = drift.pin = drift.xscale = drift.erradc = 
    drift.errpin = drift.errxscale = 0;  // initialize values in struct

  // First find drift data for strip
  
  const PulserDrift *neardrift =
    fNearDrift.GetRowByIndex(seid.BuildPlnStripEndKey());
  if(neardrift==0) {
    MAXMSG("Calib",Msg::kWarning,20) 
      << "No Near Drift Point database row for StripEnd " << seid 
      << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
    IncrementErrors(kDriftCalibrator,kMissingRow,seid);
    return GetDriftFar(seid);
  }
  MSG("Calib",Msg::kDebug) 
    << "Got Near Drift row for StripEnd " << seid 
    << " (indexed as " << seid.BuildPlnStripEndKey() << " )"
    << " aggregate. " << neardrift->GetAggregateNo() << endl;


  if(neardrift->GetNumEntries()<fMinNumPulses) {
    MAXMSG("Calib",Msg::kWarning,20) 
      << "Only " << neardrift->GetNumEntries() << " entries for  " 
      << seid
      << " (indexed as " << seid.BuildPlnStripEndKey() << " )" << endl;
    IncrementErrors(kDriftCalibrator,kDataInsufficient,seid);
    return GetDriftFar(seid);
  }

  drift.adc = neardrift->ZCMean();
  drift.erradc = neardrift->ZCError();
  if (drift.adc<1) return GetDriftFar(seid);
  if (drift.adc>8000) {
    // If the near side has a big number, try to use the far side.
    // If the far side doesn't have data, we'll have to use the near side, 
    // even though it's a big number.
       Drift fardrift = GetDriftFar(seid);
          if (fardrift.adc>0) return fardrift;
  }
  
  // Now need to get PIN info

  PlexLedId veryTempLedId = PlexLedId(neardrift->GetAggregateNo());
  PlexLedId tempLedId = PlexLedId(GetContext().GetDetector(),
                                  veryTempLedId.GetPulserBox(),
                                  veryTempLedId.GetLedInBox());
 
  FloatErr highpin=GetDriftPinDiodeValue(tempLedId, 0); // high gain
  FloatErr lowpin=GetDriftPinDiodeValue(tempLedId, 1); // low gain

  // Phil's scheme: use both pins and PulserXscale table
  if (fCalMode == 0 ) {
    const PulserPinScale *pinscale = fPinScale.GetRowByIndex(neardrift->GetAggregateNo());
    if (pinscale==0) {
      MAXMSG("Calib",Msg::kWarning,50) 
       << "No PinScale database row for StripEnd " 
       << seid 
       << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
     IncrementErrors(kDriftCalibrator,kMissingRow,seid);
     return GetDriftFar(seid);
    }  
    MSG("Calib",Msg::kDebug) 
      << "Got PinScale row for " 
      << " aggregate. " << pinscale->GetAggregateNo() << endl;
  
    drift.pin = pinscale->ScalePins(highpin.GetValue(),lowpin.GetValue());
    if (drift.pin < 0) return GetDriftFar(seid);

    const PulserXScale *xscale = 
      fXScale.GetRowByIndex(seid.BuildPlnStripEndKey());
    if(xscale==0) {
      MSG("Calib",Msg::kWarning) 
        << "No XScale database row for StripEnd " 
        << seid 
        << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
      IncrementErrors(kDriftCalibrator,kMissingRow,seid);
      return GetDriftFar(seid);
    }  
    MSG("Calib",Msg::kDebug) 
      << "Got XScale row for StripEnd " << seid 
      << " (indexed as " << seid.BuildPlnStripEndKey() << " )"
      << " aggregate. " << xscale->GetAggregateNo() << endl;

    drift.xscale = xscale->GetXScale();
    return drift; 
  }

  // Pat's scheme: use CalPulserFits table and high gain pin for preference
  else {
    // Check high-gain pin
    if(highpin.GetValue() > 100.) {
      const CalPulserFits *fitptrh = 
         fNearHigh.GetRowByIndex(seid.BuildPlnStripEndKey());
      // Check fit
      if(fitptrh!=0) {
        if(FitIsOK(*fitptrh)) {
          MSG("Calib",Msg::kDebug) 
            << "Good CalPulserFits Near-High row for StripEnd " << seid 
            << " (indexed as " << seid.BuildPlnStripEndKey() << " )"
            << " aggregate " << fitptrh->GetAggregateNo() << endl;
          drift.pin = highpin.GetValue();
          drift.errpin = highpin.GetError();
          drift.xscale = fitptrh->GetSlope();
          drift.errxscale = fitptrh->GetSlopeErr();
          return drift;
        }
      }
      if(fitptrh==0) {
        MSG("Calib",Msg::kWarning) 
          << "No CalPulserFits Near-High database row for StripEnd " 
          << seid 
          << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
        IncrementErrors(kDriftCalibrator,kMissingRow,seid);
      }
    }
    // Near-high no good: use far in preference to low-gain pin
    Drift fardrift = GetDriftFar(seid);
    if(fardrift.adc > 0) {return fardrift;}
    // Far no good, so try low-gain pin
    if(lowpin.GetValue() > 100.) {
      const CalPulserFits *fitptrl = 
         fNearLow.GetRowByIndex(seid.BuildPlnStripEndKey());
      // Check fit
      if(fitptrl!=0) {
        if(FitIsOK(*fitptrl)) {
          MSG("Calib",Msg::kInfo) 
            << "Good CalPulserFits Near-Low row for StripEnd " << seid 
            << " (indexed as " << seid.BuildPlnStripEndKey() << " )"
            << " aggregate " << fitptrl->GetAggregateNo() << endl;
          drift.pin = lowpin.GetValue();
          drift.errpin = lowpin.GetError();
          drift.xscale = fitptrl->GetSlope();
          drift.errxscale = fitptrl->GetSlopeErr();
          return drift;
        }
      }
      if(fitptrl==0) {
        MSG("Calib",Msg::kWarning) 
          << "No CalPulserFits Near-Low database row for StripEnd " 
          << seid 
          << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
        IncrementErrors(kDriftCalibrator,kMissingRow,seid);
      }
    }
    // If we get to here we have failed to find good data
    return DriftZero();
  }

}  


PulserDriftCalScheme::Drift PulserDriftCalScheme::GetDriftFar(PlexStripEndId seid) const
{
  // Return drift point bits using the far end, or zeros if we can't manage it
  Drift drift;
  drift.adc = drift.pin = drift.xscale = drift.erradc = 
    drift.errpin = drift.errxscale = 0;  // initialize values in struct

   const PulserDrift *fardrift =
     fFarDrift.GetRowByIndex(seid.BuildPlnStripEndKey());
   
   if(fardrift==0) {
     MAXMSG("Calib",Msg::kWarning,20) 
       << "No Far Drift Point database row for StripEnd " << seid 
       << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
     IncrementErrors(kDriftCalibrator,kMissingRow,seid);
     return DriftZero();
   }
   MSG("Calib",Msg::kDebug) 
     << "Got Far Drift row for StripEnd " << seid 
     << " (indexed as " << seid.BuildPlnStripEndKey() << " )"
     << " aggregate. " << fardrift->GetAggregateNo() << endl;
   if(fardrift->GetNumEntries()<fMinNumPulses) {
     MAXMSG("Calib",Msg::kWarning,20) 
       << "Only " << fardrift->GetNumEntries() << " entries for  " 
       << seid
       << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
     IncrementErrors(kDriftCalibrator,kDataInsufficient,seid);
     return DriftZero();
   }

   drift.adc = fardrift->ZCMean();
   drift.erradc = fardrift->ZCError();
   if (drift.adc < 1) return DriftZero();

   PlexLedId veryTempLedId = PlexLedId(fardrift->GetAggregateNo());
   PlexLedId tempLedId = PlexLedId(GetContext().GetDetector(),
                                   veryTempLedId.GetPulserBox(),
                                   veryTempLedId.GetLedInBox());

  //Now need to get PIN info.

  FloatErr highpin=GetDriftPinDiodeValue(tempLedId, 0); // high gain
  FloatErr lowpin=GetDriftPinDiodeValue(tempLedId, 1); // low gain

  // Phil's scheme using both pins
  if(fCalMode == 0) {
    const PulserPinScale *pinscale = fPinScale.GetRowByIndex(fardrift->GetAggregateNo());
    if (pinscale==0) {
      MAXMSG("Calib",Msg::kWarning,50) 
       << "No PinScale database row for StripEnd " 
       << seid 
       << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
     IncrementErrors(kDriftCalibrator,kMissingRow,seid);
     return DriftZero();
    }  
    MSG("Calib",Msg::kDebug) 
      << "Got PinScale row for " 
      << " aggregate. " << pinscale->GetAggregateNo() << endl;
  
    drift.pin = pinscale->ScalePins(highpin.GetValue(),lowpin.GetValue());
    if (drift.pin < 0) return DriftZero();

     const PulserXScale *xscale = 
       fFarXScale.GetRowByIndex(seid.BuildPlnStripEndKey());
     if(xscale==0) {
       MSG("Calib",Msg::kWarning) 
         << "No Far XScale database row for StripEnd " 
         << seid 
         << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
       IncrementErrors(kDriftCalibrator,kMissingRow,seid);
       return DriftZero();
     }  
     MSG("Calib",Msg::kDebug) 
       << "Got XScale row for StripEnd " << seid 
       << " (indexed as " << seid.BuildPlnStripEndKey() << " )"
       << " aggregate. " << xscale->GetAggregateNo() << endl;

     drift.xscale = xscale->GetXScale();
     return drift;
  }
  else {
    // Pat's scheme using CalPulserFits
    // Check for good high-gain pin
    if(highpin.GetValue() > 100.) {
      const CalPulserFits *fitptrh = 
         fFarHigh.GetRowByIndex(seid.BuildPlnStripEndKey());
      // Check for good fit
      if(fitptrh!=0) {
        if(FitIsOK(*fitptrh)) {
          MSG("Calib",Msg::kInfo) 
          << "Good CalPulserFits Far-High row for StripEnd " << seid 
          << " (indexed as " << seid.BuildPlnStripEndKey() << " )"
          << " aggregate " << fitptrh->GetAggregateNo() << endl;
          drift.pin = highpin.GetValue();
          drift.errpin = highpin.GetError();
          drift.xscale = fitptrh->GetSlope();
          drift.errxscale = fitptrh->GetSlopeErr();
          return drift;
        }
      }
      if(fitptrh==0) {
         MSG("Calib",Msg::kWarning) 
           << "No CalPulserFits Far-High database row for StripEnd " 
           << seid 
           << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
         IncrementErrors(kDriftCalibrator,kMissingRow,seid);
      }
    }
    // If we get to here, need to use low-gain pin
    if(lowpin.GetValue() > 100.) {
      const CalPulserFits *fitptrl = 
         fFarLow.GetRowByIndex(seid.BuildPlnStripEndKey());
      // Check for good fit
      if(fitptrl!=0) {
        if(FitIsOK(*fitptrl)) {
          MSG("Calib",Msg::kInfo) 
            << "Good CalPulserFits Far-Low row for StripEnd " << seid 
            << " (indexed as " << seid.BuildPlnStripEndKey() << " )"
            << " aggregate " << fitptrl->GetAggregateNo() << endl;
          drift.pin = lowpin.GetValue();
          drift.errpin = lowpin.GetError();
          drift.xscale = fitptrl->GetSlope();
          drift.errxscale = fitptrl->GetSlopeErr();
          return drift;
        }
      }
      if(fitptrl==0) {
        MSG("Calib",Msg::kWarning) 
          << "No CalPulserFits Far-Low database row for StripEnd " 
          << seid 
          << " (indexed as " << seid.BuildPlnStripEndKey() << " )\n";
          IncrementErrors(kDriftCalibrator,kMissingRow,seid);
      }
    }
    // If we get here nothing is good
    return DriftZero();    
  }
}

PulserDriftCalScheme::Drift PulserDriftCalScheme::DriftZero() const
{
  Drift d;
  d.adc = 0;
  d.pin = 0;
  d.xscale = 0;
  d.erradc = 0;
  d.errpin = 0;
  d.errxscale = 0;
  return d;
}

//......................................................................
Float_t PulserDriftCalScheme::GetVATemperatureCorrection() const
{
  //======================================================================
  // Purpose: Do VA temperature correction
  //======================================================================
  
  float temp = Calibrator::Instance().GetTemperature();

  float correction = 1 - fVATempCorFactor * (fRefTemp - temp);
  
  return 1./correction;

  MSG("Calib",Msg::kVerbose) << "Got Temperature Correction of: " 
                             << correction << endl;

}

//......................................................................
Bool_t PulserDriftCalScheme::CalDet2003Check() const
{
  const VldContext &vldc = GetContext();
  //is this CalDet:
  if(vldc.GetDetector()!=Detector::kCalDet) return false;
  //is this 2003:
  VldTimeStamp vldstart = VldTimeStamp(2003,8,22,0,0,0); //T7 Sept 2003 N/F
  VldTimeStamp vldend   = VldTimeStamp(2003,9,26,0,0,0); //Run Nums: 70000's
  if(vldc.GetTimeStamp()<vldstart) return false;
  if(vldc.GetTimeStamp()>=vldend) return false;
  //must be CalDet 2003
  return true;
}
//......................................................................
Bool_t PulserDriftCalScheme::FitIsOK(const CalPulserFits& fit) const
{
  if(fit.GetNPFit() < 3) return false;
  if(fit.GetChisq() < 0 || fit.GetChisq() > 100) return false;
  if(fit.GetSlope() == 0) return false;
  if(fit.GetSlopeErr()/fit.GetSlope() > 0.01) return false;
  return true;
}

---