MIPCalScheme.cxx

Go to the documentation of this file.

// $Id: MIPCalScheme.cxx,v 1.9 2005/03/22 12:29:52 tagg Exp $
//
// Prototype MIP Calibrator class
//
// rjn@hep.ucl.ac.uk
// Rebuild n.tagg1@physics.ox.ac.uk 07/04
//
// Jeff Hartnell (Oct/04):
// Added new MIPCalMethod config scheme to allow better usage 
// of task in database 
#include "MIPCalScheme.h"
#include "MessageService/MsgService.h"
#include "Plex/PlexSEIdAltL.h"
#include "Plex/PlexStripEndId.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include "Validity/VldRange.h"

ClassImp(MIPCalScheme)
CVSID("$Id: MIPCalScheme.cxx,v 1.9 2005/03/22 12:29:52 tagg Exp $");

//......................................................................
MIPCalScheme::MIPCalScheme()
{
 MSG("Calib",Msg::kVerbose) << "MIPCalScheme::MIPCalScheme" 
                                   << endl;
 Registry r;
 r.Set("defaultSigCorPerMip",565.0);
 r.Set("MIPCalMethod","TrackWindow");
 InitializeConfig(r);
}

//......................................................................
FloatErr  MIPCalScheme::GetMIP(FloatErr sigcor,
                           const PlexStripEndId& seid) const
{
  
  // If there is more than one row in the table, look up the seid
  // Otherwise, assume that the first row applies to the whole detector.

  const CalMIPCalibration* mipcal;
  

  if(fResPtr.GetNumRows()==1) {
    // There is only one row. Use it.
    mipcal = fResPtr.GetRow(0);
  } else { 
    if(!seid.IsValid()) {
      // We've been given no particular strip end. Use the first row
      mipcal = fResPtr.GetRow(0);
    } else {
      // Use the row for this particular strip end.
      mipcal = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
    }
  }

  if(mipcal ==0) {
    if(fResPtr.GetNumRows()>0) {
      MAXMSG("Calib",Msg::kWarning,10) 
        << "MIPCalScheme: No database row for StripEnd " << seid.AsString() << "\n";
      
      IncrementErrors(kMIPCalibrator,kMissingRow,seid);
    }
    
    return sigcor / (float)fDefaultSigCorPerMip;
  }
  
  return sigcor / mipcal->GetScale();
}


//......................................................................
FloatErr MIPCalScheme::DecalMIP(FloatErr mip,   
                            const PlexStripEndId& seid) const
{
  
  // If there is more than one row in the table, look up the seid
  // Otherwise, assume that the first row applies to the whole detector.

  const CalMIPCalibration* mipcal;
  
  if(fResPtr.GetNumRows()==1) {                
    // Return the only row
    mipcal = fResPtr.GetRow(0);
  } else {
    if(!seid.IsValid()) {
      // We've been handed 'no row in particular'. Use the first one.
      mipcal = fResPtr.GetRow(0);
    } else {
      // Use the row given.
      mipcal = fResPtr.GetRowByIndex(seid.BuildPlnStripEndKey());
    }
  }

  if(mipcal ==0) {
    if(fResPtr.GetNumRows()>0) {
      MAXMSG("Calib",Msg::kWarning,10) 
        << "MIPCalScheme: No database row for StripEnd " << seid.AsString() << "\n";
      
      IncrementErrors(kMIPCalibrator,kMissingRow,seid);
    }

    return mip * (float)fDefaultSigCorPerMip;
  }
  
  
  return mip * mipcal->GetScale();
}

//.....................................................................
void MIPCalScheme::DoReset( const VldContext& context )
{
  MSG("Calib",Msg::kVerbose) << "MIPCalScheme::DoReset()" << endl;
  
  fResPtr.NewQuery(context,fTask);
  
  if(fResPtr.GetNumRows()==0) {
    MAXMSG("Calib",Msg::kWarning,10) 
      << "No rows in database table with validity context "
      << context.AsString() << "  No calibration will be applied." << endl;
    IncrementErrors(kMIPCalibrator,kMissingTable);
  }  
}

//.....................................................................
void MIPCalScheme::ConfigModified()
{                                               
  bool ok = true;
  ok=ok && GetConfig().Get("defaultSigCorPerMip",fDefaultSigCorPerMip);

  const char* str_source;
  ok=ok && GetConfig().Get("MIPCalMethod",str_source);
  if(strncasecmp(str_source,"cosmic",6)==0) fTask=0;
  else if(strncasecmp(str_source,"paul",4)==0) fTask=1;
  else if(strncasecmp(str_source,"ps",2)==0) fTask=2;
  else if(strncasecmp(str_source,"alt",3)==0) fTask=3;
  else if(strncasecmp(str_source,"TrackWindow",6)==0) fTask=10;
  else {
    MSG("Calib",Msg::kError) 
      <<"I don't understand the MIPCalMethod setting: "
      <<str_source<<endl;
    MSG("Calib",Msg::kError)
      <<"Valid settings are: TrackWindow, Cosmics, Paul's Cosmics, PS muons, and Alternative PS muons"<<endl;
    MSG("Calib",Msg::kError) 
      <<"Defaulting to MIPCalMethod=window"<<endl;
    fTask=10;
  }

  if(!ok) MSG("Calib",Msg::kError)
    <<"Problem configuring MIPCalScheme."<<endl;

  // Force a re-query of the DB for this event.
  this->Reset(GetContext(),true); 
}

//.....................................................................
void MIPCalScheme::PrintConfig( std::ostream& os ) const
{
  os<<"  Using MIPCalMethod="<<GetConfig().GetCharString("MIPCalMethod")
    <<" (DB task="<<fTask<<")"<<endl
    <<"  Default SigCor/SigMap per MIP="<<fDefaultSigCorPerMip<<endl;
}

---