CalTempCalibration.cxx

Go to the documentation of this file.

// $Id: CalTempCalibration.cxx,v 1.4 2005/03/11 16:17:48 west Exp $
//
// CalTempCalibration
//
// Package: Dbi (Database Interface).

// Begin_Html<img src="../../pedestrians.gif" align=center>
// <a href="../source_warning.html">Warning for beginners</a>.<br>
// Also see <a href="../../root_crib/index.html">The ROOT Crib</a> and
// <a href="../index.html">The MINOS Class User Guide</a>End_Html
//
//rjn@hep.ucl.ac.uk
#include "Calibrator/CalTempCalibration.h"
#include "MessageService/MsgService.h"
#include "DatabaseInterface/DbiOutRowStream.h"
#include "DatabaseInterface/DbiResultSet.h"
#include "DatabaseInterface/DbiValidityRec.h"

#define REFTEMP 18.0          // reference temperature
#define TCONVERT_M -0.016     // conversion factors: temperature->light output
#define TCONVERT_ICEPT 4.342  // from linear fit

ClassImp(CalTempCalibration)

//   Definition of static data members
//   *********************************

CVSID("$Id: CalTempCalibration.cxx,v 1.4 2005/03/11 16:17:48 west Exp $\n  \
      CVSID_DBIRESULTPTR ");

//  Instantiate associated Result Pointer class.
//  *******************************************

#include "DatabaseInterface/DbiResultPtr.tpl"
template class  DbiResultPtr<CalTempCalibration>;

#include "DatabaseInterface/DbiWriter.tpl"
template class  DbiWriter<CalTempCalibration>;

// Definition of member functions (alphabetical order)
// ***************************************************


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

void CalTempCalibration::Fill(DbiResultSet& rs, 
                              const DbiValidityRec* /* vrec */) {


//
//  Purpose:  Fill object from Result Set
//
//  Arguments: 
//    rs           in    Result Set used to fill object
//    vrec         in    Associated validity record (or 0 if filling
//                                                    DbiValidityRec)
//  Return:    
//
//  Contact:   N. West
//
//  Specification:-
//  =============
//
//  o Fill object from current row of Result Set.

//  Program Notes:-
//  =============

//  This method demonstrates both the "dumb" fill method (just
//  load the data as it comes) and the smart method (check column
//  name and load according to column order). 

   if ( rs.TableName() == "CALTEMPCALIBRATION" ) {
      // Dumb method.
     rs >> fTemp >> fTempErr;
   }
   else {
      
      // Smart method
      Int_t numCol = rs.NumCols();
      //  The first column (SeqNo) has already been processed.
      for (Int_t curCol = rs.HasRowCounter() ? 3 : 2; curCol <= numCol; ++curCol) {
         string colName = rs.CurColName();
         if ( colName == "TEMP" )       rs >> fTemp;
         else if ( colName == "TEMPERR" )    rs >> fTempErr;
         else {
            MSG("Dbi",Msg::kDebug) << "Ignoring column " << curCol 
                                   << "(" << colName << ")"
                                   << "; not part of CalTempCalibration" 
                                   << endl;
         rs.IncrementCurCol();
         }
      }
   }
   
   
}


//.....................................................................
void CalTempCalibration::Store(DbiOutRowStream& ors,
                       const DbiValidityRec* /* vrec */) const {
   //
//
//  Purpose:  Stream object to output row stream
//
//  Arguments: 
//    ors          in     Output row stream.
//    vrec         in    Associated validity record (or 0 if filling
//                                                    DbiValidityRec)
//
//  Return:    
//
//  Contact:   N. West
//
//  Specification:-
//  =============
//
//  o  Stream object to output row stream.

//  Program Notes:-
//  =============

//  None.

  ors  << fTemp << fTempErr;
    
}

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

Float_t CalTempCalibration::GetCorrection() const {

  return (REFTEMP*TCONVERT_M + TCONVERT_ICEPT)
    /(fTemp*TCONVERT_M + TCONVERT_ICEPT);

}

Float_t CalTempCalibration::GetCorrected(const Float_t rawcharge) const {
//
//
//  Purpose: To apply temperature calibration to charge,NPEs,MIPs,etc
//
//  Arguments: 
//    xxxxxxxxx    in    yyyyyy
//
//  Return:    
//
//  Contact:   R.Nichol
//
//  Specification:-
//  =============
//
//  o 

//  Program Notes:-
//  =============

//  None.
  Float_t cor = GetCorrection();

   MSG("Calib",Msg::kVerbose) << " Input " << rawcharge
                              << " Temp " << fTemp
                              << " Corrected " 
                              << cor*rawcharge << "\n";
   return (cor*rawcharge);
}

---