StTpcDbElectronics.cc

/*****************************************************************
 *
 * $Id: StTpcDbElectronics.cc,v 1.7 2008/08/02 14:33:09 fisyak Exp $
 *
 * Author: Manuel Calderon de la Barca Sanchez & Brian Lasiuk Sept 13, 1999
 *
 *****************************************************************
 * Description: Electronics parameters for TRS taken from DB for
 *              the STAR Main TPC              
 *
 *****************************************************************
 *
 * $Log: StTpcDbElectronics.cc,v $
 * Revision 1.7  2008/08/02 14:33:09  fisyak
 * new interface to tpcT0
 *
 * Revision 1.6  2008/06/20 15:01:12  fisyak
 * move from StTrsData to StTpcRawData
 *
 * Revision 1.5  2000/03/15 17:39:48  calderon
 * Remove beeps
 *
 * Revision 1.4  2000/02/10 01:21:49  calderon
 * Switch to use StTpcDb.
 * Coordinates checked for consistency.
 * Fixed problems with StTrsIstream & StTrsOstream.
 *
 * Revision 1.3  2000/01/10 23:14:29  lasiuk
 * Include MACROS for compatiblity with SUN CC5
 *
 * Revision 1.2  1999/12/08 02:10:41  calderon
 * Modified to eliminate warnings on Linux.
 *
 * Revision 1.1  1999/10/11 23:55:20  calderon
 * Version with Database Access and persistent file.
 * Not fully tested due to problems with cons, it
 * doesn't find the local files at compile time.
 * Yuri suggests forcing commit to work directly with
 * files in repository.
 *
 ******************************************************************/
#include "SystemOfUnits.h"
#include "StTpcDbElectronics.hh"
//#include "StUtilities/StMessMgr.h"
#include "StTpcDb/StTpcDb.h"

#ifndef ST_NO_EXCEPTIONS
#   include <stdexcept>
#   if !defined(ST_NO_NAMESPACES)
        using std::invalid_argument;
#   endif
#endif

StTpcElectronics* StTpcDbElectronics::mInstance = 0; // static data member

StTpcDbElectronics::StTpcDbElectronics() { /* nopt */ }

StTpcDbElectronics::StTpcDbElectronics(StTpcDb* globalDbPointer)
{
    gTpcDbPtr            = globalDbPointer;
    mNominalGain         = gTpcDbPtr->Electronics()->nominalGain();
    mSamplingFrequency   = gTpcDbPtr->Electronics()->samplingFrequency();
    mTZero               = gTpcDbPtr->Electronics()->tZero();
    mAdcConversion       = gTpcDbPtr->Electronics()->adcConversion();    
    mAdcConversionCharge = gTpcDbPtr->Electronics()->adcCharge();    
    mNumberOfTimeBins    = gTpcDbPtr->Electronics()->numberOfTimeBins();    
    mAveragePedestal     = static_cast<int>(gTpcDbPtr->Electronics()->averagePedestal());
    mShapingTime         = gTpcDbPtr->Electronics()->shapingTime();
    mTau                 = gTpcDbPtr->Electronics()->tau();

#ifndef ST_NO_NAMESPACES
    using namespace units;
#endif

    //Units Integrity
    mNominalGain         *= (millivolt)/(coulomb*1.e-15);  // mV/fC
    mSamplingFrequency   *= MHz;
    mTZero               *= microsecond;
    mAdcConversion       *= (millivolt);
    mAdcConversionCharge *= (coulomb*1.e-15);
    mShapingTime         *= nanosecond;
    mTau                 *= nanosecond;
}

StTpcElectronics*
StTpcDbElectronics::instance()
{
    if (!mInstance) {
#ifndef ST_NO_EXCEPTIONS
        throw invalid_argument("StTpcDbElectronics::getInstance(): Argument Missing!");
#else
        std::cerr << "StTpcDbElectronics::getInstance(): Argument Missing!" << endl;
        std::cerr << "No arguments for instantiantion" << endl;
        std::cerr << "Exiting..." << endl;
#endif
    }
    return mInstance;
}

StTpcElectronics*
StTpcDbElectronics::instance(StTpcDb* gTpcDbPtr)
{
    if (!mInstance) {
        mInstance = new StTpcDbElectronics(gTpcDbPtr);
    }
    else {
        std::cerr << "StTpcDbElectronics::instance()"  << endl;
        std::cerr << "\tWARNING:" << endl;
        std::cerr << "\tSingleton class is already instantiated" << endl;
        std::cerr << "\tArgument  ignored!!" << endl;
        std::cerr << "\tContinuing..." << endl;
    }
    return mInstance;
}

double StTpcDbElectronics::channelGain(int sector, int row, int pad) const
{
    // This should be where channel by channel gains are looked up
    // Note: the DB has getGain(), getOnlineGain(), getNominalGain(), and
    //       getRelativeGain().  We use Nominal gain I believe.
  return 1.;//gTpcDbPtr->Gain(sector)->getNominalGain(row,pad);
}

double StTpcDbElectronics::channelGain(StTpcPadCoordinate& coord) const
{
  return channelGain(coord.sector(), coord.row(), (Int_t) coord.pad());
}

int StTpcDbElectronics::pedestal(int sector, int row, int pad, int timeB) const
{
    return averagePedestal();
}

int StTpcDbElectronics::pedestal(StTpcPadCoordinate& coord) const
{
  return pedestal(coord.sector(), coord.row(), (Int_t) coord.pad(), (Int_t) coord.timeBucket());
}
double StTpcDbElectronics::tZero(int sector, int row, int pad) const
{
  return gTpcDbPtr->tpcT0()->T0(sector,row,pad);
}
double StTpcDbElectronics::tZero(StTpcPadCoordinate& coord) const
{
  return tZero(coord.sector(), coord.row(), (Int_t) coord.pad());
}
void StTpcDbElectronics::print(ostream& os) const
{
#ifndef ST_NO_NAMESPACES
    using namespace units;
#endif
    os << "Electronics Data Base Parameters"                         << endl;
    os << "======================================="                         << endl;
    os << "Analog:"                                                         << endl;
    os << "nominalGain:         " << mNominalGain/((volt*.001)/(coulomb*1.e-15))  << " mV/fC"       << endl;
    os << "samplingFrequency:   " << mSamplingFrequency/MHz   << " MHz"         << endl;
    os << "tZero:               " << mTZero/microsecond       << " us"          << endl;
    os << "shapingTime:         " << mShapingTime/nanosecond  << " ns"          << endl;
    os << "shapingTime2:        " << mTau/nanosecond          << " ns"          << endl;

    os << "\nDigital:"                                                      << endl;
    os << "adcConversion:       " << mAdcConversion/(volt*.001) << " mV/channel" << endl;
    os << "adcConversionCharge: " << mAdcConversionCharge/(coulomb*1.e-15)  << " mV/fC"      << endl;    
    os << "numberOfTimeBins:    " << mNumberOfTimeBins                      << endl;    
    os << "averagePedestal:     " << mAveragePedestal      << " channels"   << endl;
    os << endl;
//     os << "Example of Gain per pad for Sector 1, Row 1" << endl;
//     os << " Row      Pad      Gain for this Pad" << endl;
//     os << "==========================================" << endl;
//     for(int i=0; i<gTpcDbPtr->PadPlaneGeometry()->numberOfPadsAtRow(1); i++) {
//      os.width(3);
//      os.setf(ios::right,ios::adjustfield);
//      os << 1;
//      os.width(9);
//      os << i+1;
//      os.width(15);
//      os << channelGain(1,1,i+1) << endl;
//     }
//     os << endl;
}