TemperatureCalScheme.cxx

Go to the documentation of this file.

// $Id: TemperatureCalScheme.cxx,v 1.10 2006/04/19 18:50:38 rhatcher Exp $
//
// Temperature
// n.tagg1@physics.ox.ac.uk
#include "TemperatureCalScheme.h"

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

ClassImp(TemperatureCalScheme)
CVSID("$Id: TemperatureCalScheme.cxx,v 1.10 2006/04/19 18:50:38 rhatcher Exp $");

//......................................................................
TemperatureCalScheme::TemperatureCalScheme()
{
  MSG("Calib",Msg::kVerbose) << "TemperatureCalScheme::TemperatureCalScheme" 
                             << endl;
  Registry r;
  r.Set("DefaultTemp",18.0);
  InitializeConfig(r);
}

//......................................................................
Float_t TemperatureCalScheme::GetTemperature(Int_t whichTemp) const
{

  // MC
  // The MC temperature should track the real temp.
  // But we are going to opt for the simple 'ideal detector'
  if(GetContext().GetSimFlag()==SimFlag::kMC)
    return fDefaultTemp;


  // Far

  Detector::Detector_t det = GetContext().GetDetector();

  if(det == Detector::kFar) {
    const Dcs_Env_Far* row = fDcsFarPtr.GetRow(0);
    if(!row) {
      // This shouldn't ever happen unless the table is empty.
      IncrementErrors(kThermometer,kMissingRow);
      return fDefaultTemp;
    }

    float farenheit = 0;
    switch(whichTemp) {
    default: 
      { MAXMSG("Calib",Msg::kWarning,1) << "There is no thermometer at position " << whichTemp << ". Using 1." << endl; }
    case 1: 
      farenheit = row->GetTemp1(); break;
    case 2: 
      farenheit = row->GetTemp2(); break;
    }

    // Finally, convert to celcius. (Silly Americans.)
    return (farenheit-32.0)/1.8;    
  }

  // Near
  else if(det==Detector::kNear) {
    const Dcs_Env_Near* row = fDcsNearPtr.GetRow(0);
    if(!row) {
      // This shouldn't ever happen unless the table is empty.
      IncrementErrors(kThermometer,kMissingRow);
      return fDefaultTemp;
    }

    float farenheit = row->GetTemp();

    // Finally, convert to celcius. (Silly Americans.)
    return (farenheit-32.0)/1.8;    
  }

  // CalDet
  else if(det==Detector::kCalDet) {
    if(fTempResPtr.GetNumRows()==0) {
      // We've already warned the user if this has happened, so just return.
      return fDefaultTemp;
    }
    
    const CalTempCalibration* tempcal = fTempResPtr.GetRow(0);
    if(tempcal ==0 ) {
      MAXMSG("Calib",Msg::kWarning,10) 
        << "No Temperature database row " << endl;
      IncrementErrors(kThermometer,kMissingRow);
      return fDefaultTemp;
    }
    
    MSG("Calib",Msg::kVerbose) << "Got Temperature Correction of: " 
                               << tempcal->GetTemp() << endl;
    return tempcal->GetTemp();
  }

  return fDefaultTemp;
}


//......................................................................
void TemperatureCalScheme::DoReset( const VldContext& vc )
{
  // Monte Carlo:
  if(GetContext().GetSimFlag()==SimFlag::kMC)
    return;

  // Caldet:

  if(vc.GetDetector() == Detector::kCalDet) {
    //VldTimeStamp v2hour(-7200,0);
    //v2hour.Add(vc.GetTimeStamp());
    //cout << vc.AsString() << "\t" << v2hour.AsString() << endl;
    VldContext newvc(vc.GetDetector(),vc.GetSimFlag(),vc.GetTimeStamp());
    fTempResPtr.NewQuery(newvc,0);
    if(fTempResPtr.GetNumRows()==0) {
      MAXMSG("Calib",Msg::kWarning,1) 
        << "No rows in temperature database table CALTEMPCALIBRATION"
        << vc.AsString() << ".\n";         
      IncrementErrors(kThermometer,kMissingTable);
    }  
  }

  // FD:
  else if(vc.GetDetector() == Detector::kFar) {
    fDcsFarPtr.NewQuery(vc);
    if(fDcsFarPtr.GetNumRows()==0) {
      MAXMSG("Calib",Msg::kWarning,1)
        << "No rows in DCS far detector environment database DCS_ENV_FAR. Temperature cal not avalible."
        << vc.AsString() << endl;
      IncrementErrors(kThermometer,kMissingTable);
    }
  }

  // Near:
  else if(vc.GetDetector() == Detector::kNear) {
    fDcsNearPtr.NewQuery(vc);
    if(fDcsNearPtr.GetNumRows()==0) {
      MAXMSG("Calib",Msg::kWarning,1)
        << "No rows in DCS far detector environment database DCS_ENV_NEAR. Temperature cal not avalible."
        << vc.AsString() << endl;
      IncrementErrors(kThermometer,kMissingTable);
    }
  }
}

//......................................................................
void TemperatureCalScheme::ConfigModified()
{
  Bool_t ok = true;
  ok = ok && GetConfig().Get("DefaultTemp",fDefaultTemp);
  if(!ok) MSG("Calib",Msg::kError) << "TemperatureCalScheme had a problem configuring itself!" << endl;  
}

//......................................................................
void TemperatureCalScheme::PrintConfig( std::ostream& os ) const
{
  os << "  Default Temperature = " << fDefaultTemp << endl;
}

---