In This Package:

FloatingFeePedestalSvc.cc

Go to the documentation of this file.

#include "FloatingFeePedestalSvc.h"
#include "Event/ReadoutPmtCrate.h"
#include "Event/ReadoutPmtChannel.h"
#include "Event/Readout.h"

#include "Conventions/Electronics.h"
#include "Context/TimeStamp.h"

using namespace std;
using namespace DayaBay;

FloatingFeePedestalSvc::FloatingFeePedestalSvc(const std::string& name, ISvcLocator *svc)
  : 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()
{
}

StatusCode FloatingFeePedestalSvc::initialize()
{
  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()
{
  StatusCode sc;
  sc = this->Service::finalize();

  return sc;
}

StatusCode FloatingFeePedestalSvc::queryInterface(const InterfaceID& riid,
                                                  void** ppvInterface)
{
  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;
}

// Record the preAdc value from past events,                                                                              
StatusCode FloatingFeePedestalSvc::cumulate( const DayaBay::ReadoutHeader* rh )
{
  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;
}

// Return the FEE Calibration data
double FloatingFeePedestalSvc::pedestal( const DayaBay::FeeChannelId& channelId,
                                         int gain ) const
{
  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 )
{
  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;
}

---