ROsFeeTdcTool.cc

Go to the documentation of this file.

#include "ROsFeeTdcTool.h"

#include "Conventions/Detectors.h"
#include "Conventions/Electronics.h"
#include "Conventions/Trigger.h"
#include "Conventions/Site.h"

#include "Event/SimReadoutHeader.h"
#include "Event/ReadoutPmtCrate.h"

#include "Event/ElecCrateHeader.h"
#include "Event/ElecFeeCrate.h"
#include "Event/ElecFeeChannel.h"

#include "Event/SimTrigCommand.h"
#include "Event/SimTrigCommandHeader.h"
#include "Event/SimTrigCommandCollection.h"

#include <set>
#include <vector>

#include "TF1.h"

ROsFeeTdcTool::ROsFeeTdcTool(const std::string& type,
             const std::string& name, 
             const IInterface* parent)
  : GaudiTool(type,name,parent)
{
  declareInterface< IROsFeeTdcTool >(this) ;

  declareProperty("UseMultiHitTdc",m_multiHitTdc=true,
                  "Allow multi-hit TDC functionality?");
  declareProperty("MultiLimit",m_multiLimit=5,
                  "Number of allowed TDC's in one readout");
  declareProperty("TdcResetCycles",m_tdcResetCycles=32,
                  "Reset Time on the TDC, in TDC counts (640 MHz)");
}

ROsFeeTdcTool::~ROsFeeTdcTool(){}

StatusCode ROsFeeTdcTool::initialize()
{
  m_fun = new TF1("f","-4845.33-337.431*x+13.1919*x*x-0.0950649*x*x*x",52,78);
  return StatusCode::SUCCESS;
}

StatusCode ROsFeeTdcTool::finalize()
{
  return StatusCode::SUCCESS;
}

StatusCode ROsFeeTdcTool::readoutTdc(const std::vector<int>& tdc_in,
             unsigned int first_cycle,
             unsigned int last_cycle,
             std::vector<int>& tdc_out)
{
    verbose() << "reading out tdc window ["
              << first_cycle << "," << last_cycle << ")" << endreq;

    if(tdc_out.size()>0) tdc_out.clear();
    
    // Select TDC values for readout
    std::vector<int>::const_iterator tdcIt;
    int prevTdcCycle = -2;
    int prevReadoutTdcCycle = -1000;
    for(tdcIt = tdc_in.begin(); tdcIt != tdc_in.end(); ++tdcIt){
      int tdcCycle= *tdcIt;
      // Skip TDC value if the signal has not dropped below threshold since last crossing
      if(tdcIt != tdc_in.begin())
        if( tdcCycle-*(tdcIt-1) == 1 ) {prevTdcCycle = tdcCycle; continue;} 
      // Skip TDC value if it is too close to previous threshold crossing
      double tdcReset = m_fun->GetRandom()/1.5625;
      if( (tdcCycle-prevReadoutTdcCycle)< (int)tdcReset ) {prevTdcCycle = tdcCycle; continue;}
      //if( (tdcCycle-prevReadoutTdcCycle)<m_tdcResetCycles ) {prevTdcCycle = tdcCycle; continue;}
      prevReadoutTdcCycle = tdcCycle;
      // Check if current TDC is in the readout window
      if( tdcCycle >= (int)first_cycle && tdcCycle < (int)last_cycle ){ 
        tdc_out.push_back(tdcCycle);
        // If multi-hit TDC is disabled, don't consider any more tdc values
        if( !m_multiHitTdc ) break;
        // If it reaches the number of tdc's limit, leave this.
        if( tdc_out.size() >= m_multiLimit ) break;
      }
      prevTdcCycle = tdcCycle;
    }

    verbose() << "      NTDC = " << tdc_out.size() << endreq; 
    
    return StatusCode::SUCCESS;
}