In This Package:

FloatingFeePedestalSvc Namespace Reference

Namespaces
namespace	__init__
Functions
	FloatingFeePedestalSvc (const std::string &name, ISvcLocator *svc)
	~FloatingFeePedestalSvc ()
StatusCode	initialize ()
StatusCode	finalize ()
StatusCode	queryInterface (const InterfaceID &riid, void **ppvInterface)
StatusCode	cumulate (const DayaBay::ReadoutHeader *rh)
double	pedestal (const DayaBay::FeeChannelId &channelId, int gain) const
double	updatePedestal (list< unsigned int > &PreAdcs, unsigned int newPreAdc)

---

Function Documentation

FloatingFeePedestalSvc::FloatingFeePedestalSvc	(	const std::string &	name,
		ISvcLocator *	svc
	)

Definition at line 17 of file FloatingFeePedestalSvc.cc.

  : Service(name,svc)
{
  declareProperty("Limit",m_limit=500,
                  "Number of preAdcs to document");

  declareProperty("TimeCutInSec",m_timeCutInSec=0.0001,  // 100 us
                  "Time cut in second to remove afterpulses");
}

FloatingFeePedestalSvc::~FloatingFeePedestalSvc

(

)

Definition at line 27 of file FloatingFeePedestalSvc.cc.

{
}

StatusCode FloatingFeePedestalSvc::initialize

(

)

Definition at line 31 of file FloatingFeePedestalSvc.cc.

{
  StatusCode sc;
  sc = this->Service::initialize();

  m_finePreAdc.clear();
  m_crsePreAdc.clear();

  m_finePedes.clear();
  m_crsePedes.clear();

  m_lastTrig = TimeStamp::GetBOT();

  return sc;
}

StatusCode FloatingFeePedestalSvc::finalize

(

)

Definition at line 47 of file FloatingFeePedestalSvc.cc.

{
  StatusCode sc;
  sc = this->Service::finalize();

  return sc;
}

StatusCode FloatingFeePedestalSvc::queryInterface	(	const InterfaceID &	riid,
		void **	ppvInterface
	)

Definition at line 55 of file FloatingFeePedestalSvc.cc.

{
  StatusCode sc = StatusCode::FAILURE;
  if (ppvInterface) {
    *ppvInterface = 0;
    
    if (IFloatingFeePedestalSvc::interfaceID().versionMatch(riid)) {
      *ppvInterface = static_cast<IFloatingFeePedestalSvc*>(this);
      sc = StatusCode::SUCCESS;
      addRef();
    }
    else sc = Service::queryInterface( riid, ppvInterface );
  }
  return sc;
}

StatusCode FloatingFeePedestalSvc::cumulate

(

const DayaBay::ReadoutHeader *

rh

)

Definition at line 73 of file FloatingFeePedestalSvc.cc.

{
  MsgStream log(msgSvc(), "FloatingFeePedestalSvc");

  const Readout* readout = rh->readout();
  if(!readout) {
    log << MSG::ERROR <<"Failed to get readout from header"<<endreq;
    return StatusCode::FAILURE;
  }

  if( readout->detector().detectorId() != DetectorId::kAD1 &&
      readout->detector().detectorId() != DetectorId::kAD2 &&
      readout->detector().detectorId() != DetectorId::kAD3 &&
      readout->detector().detectorId() != DetectorId::kAD4 &&
      readout->detector().detectorId() != DetectorId::kIWS &&
      readout->detector().detectorId() != DetectorId::kOWS ) {
    return StatusCode::SUCCESS;
  }

  TimeStamp trigTime = readout->triggerTime();
  
  if( (m_lastTrig == TimeStamp::GetBOT()) ||
      ((trigTime-m_lastTrig).GetSeconds() > m_timeCutInSec) ) {
    
    const ReadoutPmtCrate* pPmtCrate=0;
    pPmtCrate = dynamic_cast< const DayaBay::ReadoutPmtCrate* > ( readout );
    if(!pPmtCrate) {
      log << MSG::ERROR <<"Failed to get ReadoutPmtCrate"<<endreq;
      return StatusCode::FAILURE;
    }
    
    ReadoutPmtCrate::PmtChannelReadouts channels = pPmtCrate->channelReadout();
    ReadoutPmtCrate::PmtChannelReadouts::iterator ci, ci_end = channels.end();

    for(ci = channels.begin(); ci!=ci_end; ci++) {

      ReadoutPmtChannel channel = ci->second;
      FeeChannelId channelId = channel.channelId();
      
      int preAdc = channel.pedestal()[0];
      int range  = channel.adcRange()[0];

      if( range == FeeGain::kHigh ) {
        list<unsigned int>& preAdcs = m_finePreAdc[ channelId ];
        m_finePedes[ channelId ] = updatePedestal( preAdcs, preAdc );
      }
      if( range == FeeGain::kLow ) {
        list<unsigned int>& preAdcs = m_crsePreAdc[ channelId ];
        m_crsePedes[ channelId ] = updatePedestal( preAdcs, preAdc );
      }
    }
    
    m_lastTrig = trigTime;
  }

  return StatusCode::SUCCESS;
}

double FloatingFeePedestalSvc::pedestal	(	const DayaBay::FeeChannelId &	channelId,
		int	gain
	)			const

Definition at line 140 of file FloatingFeePedestalSvc.cc.

{
  MsgStream log(msgSvc(), "FloatingFeePedestalSvc");

  if( gain==FeeGain::kHigh ) {
    PedestalArray::const_iterator ci;
    ci = m_finePedes.find( channelId );
    if( ci != m_finePedes.end() ) {
      return ci->second;
    } else {
      log << MSG::WARNING <<"Incomplete fine range pedestal map. Try preADC."<<endreq;
      return -1000;
    }
  }

  if( gain==FeeGain::kLow ) {
    PedestalArray::const_iterator ci;
    ci = m_crsePedes.find( channelId );
    if( ci != m_crsePedes.end() ) {
      return ci->second;
    } else {
      log << MSG::WARNING <<"Incomplete coarse range pedestal map. Try preADC."<<endreq;
      return -1000;
    }
  }

  log << MSG::ERROR <<"Query with wrong gain parameter"<<endreq;
  return -1000;
}

double FloatingFeePedestalSvc::updatePedestal	(	list< unsigned int > &	PreAdcs,
		unsigned int	newPreAdc
	)

Definition at line 171 of file FloatingFeePedestalSvc.cc.

{
  PreAdcs.push_front( newPreAdc );
  if( PreAdcs.size()>m_limit ) PreAdcs.pop_back();

  double sum=0;
  int count = PreAdcs.size();
  list<unsigned int>::iterator li,liend=PreAdcs.end();
  for( li=PreAdcs.begin(); li!=liend; li++ ) {
    sum+= (*li);
  }
  
  sum = sum/count;
  
  return sum;
}

---