In This Package:

SimpPmtCalibTool Class Reference

#include <SimpPmtCalibTool.h>

Inheritance diagram for SimpPmtCalibTool:

[legend]Collaboration diagram for SimpPmtCalibTool:[legend]List of all members.

Public Types
	SUCCESS
	NO_INTERFACE
	VERSMISMATCH
	LAST_ERROR
	SUCCESS
	NO_INTERFACE
	VERSMISMATCH
	LAST_ERROR
enum	Status
Public Member Functions
	SimpPmtCalibTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~SimpPmtCalibTool ()
virtual StatusCode	calibration (const DayaBay::DaqCrate *, DayaBay::CalibReadout *)
	This is the extension. Sub classes must provide it.
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
INTupleSvc *	ntupleSvc () const
INTupleSvc *	evtColSvc () const
IDataProviderSvc *	detSvc () const
IDataProviderSvc *	evtSvc () const
IIncidentSvc *	incSvc () const
IChronoStatSvc *	chronoSvc () const
IHistogramSvc *	histoSvc () const
IAlgContextSvc *	contextSvc () const
DataObject *	put (IDataProviderSvc *svc, DataObject *object, const std::string &address, const bool useRootInTES=true) const
DataObject *	put (DataObject *object, const std::string &address, const bool useRootInTES=true) const
Gaudi::Utils::GetData< TYPE >::return_type	get (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
Gaudi::Utils::GetData< TYPE >::return_type	get (const std::string &location, const bool useRootInTES=true) const
TYPE *	getDet (IDataProviderSvc *svc, const std::string &location) const
TYPE *	getDet (const std::string &location) const
bool	exist (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
bool	exist (const std::string &location, const bool useRootInTES=true) const
bool	existDet (IDataProviderSvc *svc, const std::string &location) const
bool	existDet (const std::string &location) const
TYPE *	getOrCreate (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
TYPE *	getOrCreate (const std::string &location, const bool useRootInTES=true) const
TOOL *	tool (const std::string &type, const std::string &name, const IInterface *parent=0, bool create=true) const
TOOL *	tool (const std::string &type, const IInterface *parent=0, bool create=true) const
SERVICE *	svc (const std::string &name, const bool create=true) const
IUpdateManagerSvc *	updMgrSvc () const
IDataProviderSvc *	fastContainersSvc () const
StatusCode	Error (const std::string &msg, const StatusCode st=StatusCode::FAILURE, const size_t mx=10) const
StatusCode	Warning (const std::string &msg, const StatusCode st=StatusCode::FAILURE, const size_t mx=10) const
StatusCode	Print (const std::string &msg, const StatusCode st=StatusCode::SUCCESS, const MSG::Level lev=MSG::INFO) const
StatusCode	Assert (const bool ok, const std::string &message="", const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
StatusCode	Assert (const bool ok, const char *message, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
StatusCode	Exception (const std::string &msg, const GaudiException &exc, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
StatusCode	Exception (const std::string &msg, const std::exception &exc, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
StatusCode	Exception (const std::string &msg="no message", const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
MsgStream &	msgStream (const MSG::Level level) const
MsgStream &	always () const
MsgStream &	fatal () const
MsgStream &	err () const
MsgStream &	error () const
MsgStream &	warning () const
MsgStream &	info () const
MsgStream &	debug () const
MsgStream &	verbose () const
MsgStream &	msg () const
const Statistics &	counters () const
StatEntity &	counter (const std::string &tag) const
MSG::Level	msgLevel () const
bool	msgLevel (const MSG::Level level) const
void	resetMsgStream () const
bool	typePrint () const
bool	propsPrint () const
bool	statPrint () const
bool	errorsPrint () const
long	printStat (const MSG::Level level=MSG::ALWAYS) const
long	printErrors (const MSG::Level level=MSG::ALWAYS) const
long	printProps (const MSG::Level level=MSG::ALWAYS) const
void	registerCondition (const std::string &condition, StatusCode(CallerClass::*mf)()=NULL)
void	registerCondition (const std::string &condition, CondType *&condPtrDest, StatusCode(CallerClass::*mf)()=NULL)
void	registerCondition (char *condition, StatusCode(CallerClass::*mf)()=NULL)
void	registerCondition (TargetClass *condition, StatusCode(CallerClass::*mf)()=NULL)
StatusCode	runUpdate ()
TransientFastContainer< T > *	getFastContainer (const std::string &location, typename TransientFastContainer< T >::size_type initial=0)
StatusCode	release (const IInterface *interface) const
virtual unsigned long	release ()
const std::string &	context () const
const std::string &	rootInTES () const
double	globalTimeOffset () const
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
virtual unsigned long	addRef ()
virtual const std::string &	name () const
virtual const std::string &	type () const
virtual const IInterface *	parent () const
virtual StatusCode	configure ()
virtual StatusCode	start ()
virtual StatusCode	stop ()
virtual StatusCode	terminate ()
virtual StatusCode	reinitialize ()
virtual StatusCode	restart ()
virtual Gaudi::StateMachine::State	FSMState () const
virtual Gaudi::StateMachine::State	targetFSMState () const
virtual StatusCode	sysInitialize ()
virtual StatusCode	sysStart ()
virtual StatusCode	sysStop ()
virtual StatusCode	sysFinalize ()
virtual StatusCode	sysReinitialize ()
virtual StatusCode	sysRestart ()
virtual StatusCode	setProperty (const Property &p)
virtual StatusCode	setProperty (const std::string &s)
virtual StatusCode	setProperty (const std::string &n, const std::string &v)
StatusCode	setProperty (const std::string &name, const TYPE &value)
virtual StatusCode	getProperty (Property *p) const
virtual const Property &	getProperty (const std::string &name) const
virtual StatusCode	getProperty (const std::string &n, std::string &v) const
virtual const std::vector< Property * > &	getProperties () const
PropertyMgr *	getPropertyMgr ()
ISvcLocator *	serviceLocator () const
ISvcLocator *	svcLoc () const
IMessageSvc *	msgSvc () const
IToolSvc *	toolSvc () const
StatusCode	setProperties ()
StatusCode	service (const std::string &name, T *&svc, bool createIf=true) const
StatusCode	service (const std::string &type, const std::string &name, T *&svc) const
void	declInterface (const InterfaceID &, void *)
Property *	declareProperty (const std::string &name, T &property, const std::string &doc="none") const
Property *	declareRemoteProperty (const std::string &name, IProperty *rsvc, const std::string &rname="") const
IAuditorSvc *	auditorSvc () const
IMonitorSvc *	monitorSvc () const
void	declareInfo (const std::string &name, const T &var, const std::string &desc) const
void	declareInfo (const std::string &name, const std::string &format, const void *var, int size, const std::string &desc) const
virtual const std::string &	type () const =0
virtual const IInterface *	parent () const =0
virtual StatusCode	configure ()=0
virtual StatusCode	start ()=0
virtual StatusCode	stop ()=0
virtual StatusCode	terminate ()=0
virtual StatusCode	reinitialize ()=0
virtual StatusCode	restart ()=0
virtual Gaudi::StateMachine::State	FSMState () const =0
virtual StatusCode	sysInitialize ()=0
virtual StatusCode	sysStart ()=0
virtual StatusCode	sysStop ()=0
virtual StatusCode	sysFinalize ()=0
virtual StatusCode	sysReinitialize ()=0
virtual StatusCode	sysRestart ()=0
virtual unsigned long	refCount () const =0
virtual const std::string &	name () const =0
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvInterface)=0
virtual unsigned long	addRef ()=0
virtual unsigned long	release ()=0
Static Public Member Functions
static const InterfaceID &	interfaceID ()
static const InterfaceID &	interfaceID ()
static const InterfaceID &	interfaceID ()
static const InterfaceID &	interfaceID ()
	Retrieve interface ID.
Public Attributes
	SUCCESS
	NO_INTERFACE
	VERSMISMATCH
	LAST_ERROR
Protected Types
typedef std::map< std::string, StatEntity >	Statistics
typedef std::map< std::string, unsigned int >	Counter
typedef std::vector< IAlgTool * >	AlgTools
typedef std::pair< IInterface *, std::string >	ServiceEntry
typedef std::vector< ServiceEntry >	Services
Protected Member Functions
StatusCode	releaseTool (const IAlgTool *tool) const
StatusCode	releaseSvc (const IInterface *svc) const
int	outputLevel () const
virtual unsigned long	refCount () const
IntegerProperty &	outputLevelProperty ()
void	initOutputLevel (Property &prop)
Static Protected Attributes
static const bool	IgnoreRootInTES
static const bool	UseRootInTES
Private Attributes
std::string	m_calibDataSvcName
ICalibDataSvc *	m_calibDataSvc
std::string	m_cableSvcName
bool	m_useDynamicPedestal
bool	m_useUncalibratedCharge
bool	m_useUncalibratedTime
ICableSvc *	m_cableSvc

---

Detailed Description

Definition at line 26 of file SimpPmtCalibTool.h.

---

Constructor & Destructor Documentation

SimpPmtCalibTool::SimpPmtCalibTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 31 of file SimpPmtCalibTool.cc.

  :  GaudiTool(type, name, parent)
  ,  m_calibDataSvc(0)
  ,  m_cableSvc(0)
{
  declareInterface< ICalibTool >(this);
  declareProperty("CalibDataSvcName", m_calibDataSvcName = "StaticCalibDataSvc",
                  "Name of Pmt calibration data Service");
  declareProperty("CableSvcName", m_cableSvcName = "StaticCableSvc",
                  "Name of Cable Service");
  declareProperty("UseDynamicPedestal", m_useDynamicPedestal=true, 
                  "Use pedestal value embedded with each hit"); 
  declareProperty("UseUncalibratedCharge", m_useUncalibratedCharge=false, 
                  "Use approximate charge calibration"); 
  declareProperty("UseUncalibratedTime", m_useUncalibratedTime=false, 
                  "Use approximate time calibration"); 
}

SimpPmtCalibTool::~SimpPmtCalibTool

(

)

[virtual]

Definition at line 51 of file SimpPmtCalibTool.cc.

{}

---

Member Function Documentation

StatusCode SimpPmtCalibTool::calibration	(	const DayaBay::DaqCrate *	,
		DayaBay::CalibReadout *
	)			[virtual]

This is the extension. Sub classes must provide it.

Implements ICalibTool.

Definition at line 72 of file SimpPmtCalibTool.cc.

{
  debug() << "running calibration() " << endreq;
  StatusCode sc;
  
  // Convert to PMT crate readout
  const DayaBay::DaqPmtCrate* pmtReadout = readout->asPmtCrate();
  if(!pmtReadout){
    error() << "calibration():  no pmt readouts." << endreq;
    return StatusCode::FAILURE;
  }
  
  DayaBay::CalibReadoutPmtCrate* calibReadoutCrate 
    = dynamic_cast<DayaBay::CalibReadoutPmtCrate*>(calibReadout);
  // Check CalibReadoutCrate
  if(calibReadoutCrate == NULL) {
    error() << "calibration():  no calib pmt readouts." << endreq;
    return StatusCode::FAILURE;
  }
  
  const std::vector<DayaBay::DaqPmtChannel*>& channels 
    = pmtReadout->channelReadouts();
  std::vector<DayaBay::DaqPmtChannel*>::const_iterator channelIter, 
    channelEnd = channels.end();
  
  int task = 0;
  Context context = calibReadoutCrate->header()->context();
  ServiceMode svcMode(context, task);
  
  // Simple calibration channel by channel
  for(channelIter = channels.begin();
      channelIter!=channelEnd; 
      channelIter++){ 
    const DayaBay::DaqPmtChannel* channel = *channelIter; 
    const DayaBay::FeeChannelId& channelId = channel->channelId();
    
    // Get pmt calib data from CalibDataSvc
    const DayaBay::DetectorSensor& sensDetId 
      = m_cableSvc->sensor(channelId, svcMode);
    
    calibReadoutCrate->addSensor(sensDetId);
    DayaBay::CalibReadoutPmtChannel& cro_ch 
      = *(calibReadoutCrate->sensor(sensDetId)); 
    
    const DayaBay::PmtCalibData* pmtCalib = 0;
    if(!m_useUncalibratedCharge || !m_useUncalibratedTime){
      pmtCalib = m_calibDataSvc->pmtCalibData(sensDetId, svcMode);
      if(!pmtCalib) {
        error() << "No calibration properties for PMT: " << sensDetId 
                << endreq;
        return StatusCode::FAILURE;
      }
    }
    
    const DayaBay::FeeCalibData* feeCalib = 0; 
    if(!m_useDynamicPedestal){ 
      feeCalib = m_calibDataSvc->feeCalibData(channelId, svcMode); 
      if(!feeCalib) {
        error() << "No calibration properties for FEE channel: " 
                << channelId << endreq;
        return StatusCode::FAILURE;
      }
    }
    
    // simple tdc calibration
    std::vector<double> calib_time;
    std::vector<double> calib_charge;
    double tdcToTime = -(1.0/DayaBay::TdcFrequencyHz) * Units::second;
    for(unsigned int hitIdx=0; hitIdx<channel->hitCount(); hitIdx++){
      // Convert TDC to nanoseconds, and remove PMT time offset
      double calibTime = channel->tdc(hitIdx) * tdcToTime;
      if(!m_useUncalibratedTime) calibTime -= pmtCalib->m_timeOffset;
      calib_time.push_back( calibTime ); 
      
      double adcBaseline = 0.;
      double adcToSpe = 0.;
      if(channel->adc(hitIdx)==0) {
        // Catch invalid ADC values, set to zero charge
        calib_charge.push_back(0.);
        continue;
      }
      // Find ADC pedestal for this channel 
      if(m_useDynamicPedestal){
        adcBaseline = channel->preAdcAvg(hitIdx); 
      }else{
        if( channel->isHighGainAdc(hitIdx) ){
          adcBaseline = feeCalib->m_adcBaselineHigh;
        }else{
          adcBaseline = feeCalib->m_adcBaselineLow;
        }
      }
      
      // Find ADC to charge conversion
      if(m_useUncalibratedCharge){
        // Use mean charge conversion, instead of PMT calibration parameters
        adcToSpe = 1. / 20.;
        if(! channel->isHighGainAdc(hitIdx) ) adcToSpe = 1.0;
      }else{
        // Use PMT calibration parameters
        if(channel->isHighGainAdc(hitIdx)) {  //high
          adcToSpe = 1.0 / pmtCalib->m_speHigh;
        }else{
          adcToSpe = 1.0 / pmtCalib->m_speLow;
        }
      }
      calib_charge.push_back( (channel->adc(hitIdx) - adcBaseline) * adcToSpe);
      
    }
    cro_ch.setTime(calib_time);
    cro_ch.setCharge(calib_charge);
  }
  
  return StatusCode::SUCCESS;
}

StatusCode SimpPmtCalibTool::initialize

(

)

[virtual]

Reimplemented from GaudiTool.

Definition at line 53 of file SimpPmtCalibTool.cc.

{
  // Get CalibData Service
  m_calibDataSvc = svc<ICalibDataSvc>(m_calibDataSvcName, true);

  // Get CableData Service
  m_cableSvc = svc<ICableSvc>(m_cableSvcName, true);

  debug() << "successfully get CableSvc and CalibSvc" << endreq;

  return StatusCode::SUCCESS;
}

StatusCode SimpPmtCalibTool::finalize

(

)

[virtual]

Reimplemented from GaudiTool.

Definition at line 67 of file SimpPmtCalibTool.cc.

{
  return StatusCode::SUCCESS;
}

const InterfaceID & ICalibTool::interfaceID

(

)

[static, inherited]

Retrieve interface ID.

Reimplemented from IAlgTool.

Definition at line 8 of file ICalibTool.cc.

{ 
    return IID_ICalibTool; 
}

---

Member Data Documentation

std::string SimpPmtCalibTool::m_calibDataSvcName [private]

Definition at line 44 of file SimpPmtCalibTool.h.

ICalibDataSvc* SimpPmtCalibTool::m_calibDataSvc [private]

Definition at line 47 of file SimpPmtCalibTool.h.

std::string SimpPmtCalibTool::m_cableSvcName [private]

Definition at line 50 of file SimpPmtCalibTool.h.

bool SimpPmtCalibTool::m_useDynamicPedestal [private]

Definition at line 53 of file SimpPmtCalibTool.h.

bool SimpPmtCalibTool::m_useUncalibratedCharge [private]

Definition at line 56 of file SimpPmtCalibTool.h.

bool SimpPmtCalibTool::m_useUncalibratedTime [private]

Definition at line 59 of file SimpPmtCalibTool.h.

ICableSvc* SimpPmtCalibTool::m_cableSvc [private]

Definition at line 62 of file SimpPmtCalibTool.h.

---

The documentation for this class was generated from the following files:

• /home/dayabay/installs/trunk_2011_04_11_opt/NuWa-trunk/dybgaudi/Calibration/CalibAlg/src/components/SimpPmtCalibTool.h
• /home/dayabay/installs/trunk_2011_04_11_opt/NuWa-trunk/dybgaudi/Calibration/CalibAlg/src/components/SimpPmtCalibTool.cc

---