StripToStripCalScheme.cxx

Go to the documentation of this file.

// $Id: StripToStripCalScheme.cxx,v 1.11 2007/03/29 19:56:13 tagg Exp $
//
// The 'new' strip-to-strip calibration scheme.
//
// Applies the numbers from the CALSTRIPTOSTRIP table.
// These numbers are such that they take out ALL variations between
// strips, normalized to the center of the strip.
//
// i.e. for an event in the center of the strip, this calibration should
// remove all strip-to-strip effects.
//
// In other words, this calibration takes out:
// 1) Gain
// 2) All fibre attenuation from pixel to middle of strip
// 3) Scintillator light output variation strip-to-strip.
//
// This constrasts with the MuonSigCorrCalScheme,
// which takes out:
// 1) Gain
// 2) Corrections to mapper data
//
// n.tagg1@physics.ox.ac.uk
#include "StripToStripCalScheme.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(StripToStripCalScheme)
CVSID("$Id: StripToStripCalScheme.cxx,v 1.11 2007/03/29 19:56:13 tagg Exp $");

//......................................................................
StripToStripCalScheme::StripToStripCalScheme()
{
  MSG("Calib",Msg::kVerbose) << "StripToStripCalScheme::StripToStripCalScheme" 
                                  << endl;
  Registry r;
  r.Set("TemperatureCalibration","off");
  r.Set("MuonSource","Cosmics");
  r.Set("SpectrometerFudge",1.0);
  InitializeConfig(r);

  // Hard code these for now.
  fRefTemp = 18.0;       // reference temperature
  fTempCorFactor = -0.0008;  // frac diff in light output with temperature

}

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

  MSG("Calib",Msg::kVerbose) << "StripToStripCalScheme::DoReset(VldContext) Task: " << fMuonTask << endl;

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

//......................................................................
void StripToStripCalScheme::ConfigModified()
{

  const char* str_tempcal;
  const char* str_source;


  bool ok = true;
  ok = ok && GetConfig().Get("TemperatureCalibration",str_tempcal);
  ok = ok && GetConfig().Get("MuonSource",str_source);
  ok = ok && GetConfig().Get("SpectrometerFudge",fSpectrometerFudge);
  if(!ok) MSG("Calib",Msg::kError) << "Problem setting up StripToStripCalScheme config." << endl;
  
  // Interpret strings.
  if(strcasecmp(str_tempcal,"on")==0) fDoTempCal = true;
  else fDoTempCal = false;
  
  if(strncasecmp(str_source,"cosmic",6)==0) fMuonTask = 0;
  else if(strncasecmp(str_source,"paul",4)==0) fMuonTask = 1;
  else if(strncasecmp(str_source,"ps",2)==0) fMuonTask = 2;
  else if(strncasecmp(str_source,"alt",3)==0) fMuonTask = 3;
  else {
    MSG("Calib",Msg::kError) << "I don't understand the MuonSource setting " << str_source << endl;
    MSG("Calib",Msg::kError) << "Valid settings are: Cosmics, Paul's Cosmics, PS muons, and Alternative PS muons"<< endl;
    MSG("Calib",Msg::kError) << "Defaulting to cosmics"<< endl;
    fMuonTask=0;
  }

  // Ensure that the DB has been changed for this event.
  Reset(GetContext(),true); 
}

//......................................................................
void StripToStripCalScheme::PrintConfig( std::ostream& os ) const
{
  os << " (NEW Strip-to-strip calibration scheme) " << endl;
  os << "  Muon Source = " << fMuonTask << " = " 
     << GetConfig().GetCharString("MuonSource") << endl;
  os << " Temperature calibration is " << ((fDoTempCal)?"on":"off") << endl;
  if(fDoTempCal) {
    os << "  Temp Correction Factor = 1 / [ 1 - " 
       << fTempCorFactor << "*(" 
       << fRefTemp << " - temp) ]" << endl;
  }
  if(fSpectrometerFudge!=1.0) {
      os << "  Spectrometer fudged by ratio: " << fSpectrometerFudge << endl;
  }
}


//......................................................................
FloatErr StripToStripCalScheme::GetStripToStripCorrected(FloatErr siglin, 
                                                       const PlexStripEndId& seid) const
{
  

  MSG("Calib",Msg::kVerbose) << "StripToStripCalScheme on " << siglin 
                             << " siglin in seid " << seid.BuildPlnStripEndKey() <<"\n";
  // Now need to get the row which corresponds to the stripendnum.

  FloatErr correction = 1.0;
  if(fDoTempCal)
    correction = GetTemperatureCorrection();

  // Don't complain about veto shield entries; just do it.
  if(seid.IsVetoShield()) return  siglin * correction;

  if(seid.IsNearSpect()) correction *= fSpectrometerFudge;

  const CalStripToStrip* muoncal = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
  
  if(muoncal ==0) {                             
    if(fResPtr.GetNumRows()>0) {
      MAXMSG("Calib",Msg::kWarning,10) 
       << "StripToStripCalScheme: No database row for StripEnd " 
       << seid.AsString("c") << " key " << seid.BuildPlnStripEndKey() << "\n";
     
     IncrementErrors(kStripCalibrator,kMissingRow,seid);
   }

    return siglin * correction * FloatErr(1,0.5); // Add 50% error 
 }

  
  return correction * siglin / FloatErr(muoncal->GetResponse(), muoncal->GetResponseErr());
}

//......................................................................
FloatErr StripToStripCalScheme::DecalStripToStrip(FloatErr sigcorr,
                                                const PlexStripEndId& seid) const
{

  MSG("Calib",Msg::kVerbose) << "StripToStripCalScheme decal on " << sigcorr 
                             << " siglin in seid " << seid.AsString() <<"\n";
  
  FloatErr correction = 1.0;
  if(fDoTempCal)
    correction = GetTemperatureCorrection();

  // Don't complain about veto shield entries; just do it.
  if(seid.IsVetoShield()) return  sigcorr / correction;

  if(seid.IsNearSpect()) correction *= fSpectrometerFudge;

  // Now need to get the row which corresponds to the stripendnum.
  const CalStripToStrip* muoncal = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
  
  if(muoncal ==0) {
    if(fResPtr.GetNumRows()>0) {
      MAXMSG("Calib",Msg::kWarning,10) 
        << "StripToStripCalScheme: No database row for StripEnd " << seid.BuildPlnStripEndKey() << "\n";
      
     IncrementErrors(kStripCalibrator,kMissingRow,seid);
    }
    
    return sigcorr / correction * FloatErr(1,0.5); // Add 50% for unknown.
  }
  
  return sigcorr / correction * FloatErr(muoncal->GetResponse(), muoncal->GetResponseErr());
}

//......................................................................
FloatErr StripToStripCalScheme::GetTemperatureCorrection() const
{
  //======================================================================
  // Purpose: Do temperature correction
  //======================================================================
  
  float temp = Calibrator::Instance().GetTemperature();
  
  float correction = 1 - fTempCorFactor * (fRefTemp - temp);
  
  return 1./correction;

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

}

---