Default CLB simulation. More...

#include <JCLBDefaultSimulatorInterface.hh>

Inheritance diagram for JDETECTOR::JCLBDefaultSimulatorInterface:

Classes
struct	compare
	Auxiliary data structure for sorting of hits. More...

class	JStateMachine
	Interface to mimic hit ordering effects due to state machine inside CLB. More...

class	JTDC
	Interface for TDC. More...

Public Types
typedef JDAQHit::JPMT_t	JPMT_t

typedef JDAQHit::JTDC_t	JTDC_t

typedef JDAQHit::JTOT_t	JTOT_t

Public Member Functions
	JCLBDefaultSimulatorInterface (JTDC TDC, JStateMachine state_machine)
	Constructor. More...

JDAQFrameStatus	getDAQFrameStatus (const JModuleIdentifier &id) const
	Get DAQ frame status of given module. More...

virtual void	processData (const JModuleIdentifier &id, const JCLBInput &input, JDAQSuperFrame &output) const override
	Process data. More...

virtual int	getUDPNumberOfReceivedPackets (const JModuleIdentifier &id) const
	Get number of received UDP packets. More...

virtual int	getUDPMaximalSequenceNumber (const JModuleIdentifier &id) const
	Get maximal sequence number of UDP packet. More...

virtual bool	hasUDPTrailer (const JModuleIdentifier &id) const
	Get UDP trailer status. More...

virtual bool	getHighRateVeto (const JPMTIdentifier &id) const
	Get high-rate veto of given PMT. More...

Static Public Member Functions
static double	getMinimalToT ()
	Get minimal pulse length of time-over-threshold measurement. More...

Static Public Attributes
static struct JDETECTOR::JCLBDefaultSimulatorInterface::compare	compare

Private Attributes
JLANG::JSinglePointer< JTDC >	TDC

JLANG::JSinglePointer < JStateMachine >	state_machine

std::vector< JDAQHit >	buffer

Detailed Description

Default CLB simulation.

This class provides for a default implementation of the JCLBSimulator interface which is based on a simulation of the TDC and the state machine inside the CLB. The actual number of hits may change due to the high-rate veto and UDP packet loss (loss of hits) and the dynamic range of the time-over-threshold (gain of hits).

The nested class JStateMachine constitutes a user interface for the simulation of the state machine through method JStateMachine::maybeSwapped(). With the default implementation, the overall time ordering if hits is maintained. For a realistic simulation of the CLB, a pointer to a designated implementation of this interface should be provided.

The implementation of the virtual method JCLBDefaultSimulatorInterface::processData provides for the settings of the status of the data frame.
In this, the high-rate veto is set when:

virtual method JCLBDefaultSimulatorInterface::getHighRateVeto returns true for the given PMT; or
number of hits from the PMT exceeds KM3NETDAQ::getMaximalNumberOfHits.

Definition at line 55 of file JCLBDefaultSimulatorInterface.hh.

Member Typedef Documentation

typedef JDAQHit::JPMT_t JDETECTOR::JCLBDefaultSimulatorInterface::JPMT_t

Definition at line 60 of file JCLBDefaultSimulatorInterface.hh.

typedef JDAQHit::JTDC_t JDETECTOR::JCLBDefaultSimulatorInterface::JTDC_t

Definition at line 61 of file JCLBDefaultSimulatorInterface.hh.

typedef JDAQHit::JTOT_t JDETECTOR::JCLBDefaultSimulatorInterface::JTOT_t

Definition at line 62 of file JCLBDefaultSimulatorInterface.hh.

Constructor & Destructor Documentation

JDETECTOR::JCLBDefaultSimulatorInterface::JCLBDefaultSimulatorInterface	(	JTDC *	TDC,
		JStateMachine *	state_machine
	)

inline

Constructor.

This class owns the objects pointed to.

Parameters

TDC	pointer to TDC simulator
state_machine	pointer to state machine

Definition at line 128 of file JCLBDefaultSimulatorInterface.hh.

                                                                            :
       JCLBSimulator(),
       TDC          (TDC),
       state_machine(state_machine)
     {}

Member Function Documentation

JDAQFrameStatus JDETECTOR::JCLBDefaultSimulatorInterface::getDAQFrameStatus ( const JModuleIdentifier & id ) const

inline

Get DAQ frame status of given module.

Parameters

id	module identifier

Returns: DAQ frame status

Definition at line 142 of file JCLBDefaultSimulatorInterface.hh.

     {
       using namespace KM3NETDAQ;
 
       int  daq    = (DAQ_UDP_RECEIVED_PACKETS.write(getUDPNumberOfReceivedPackets(id)) |    //
                      DAQ_UDP_SEQUENCE_NUMBER .write(getUDPMaximalSequenceNumber  (id)));    // DAQ  status
       int  status =  DAQ_WHITE_RABBIT.write(1);                                             // TDC  status
       int  fifo   =  DAQ_UDP_TRAILER .write(hasUDPTrailer(id) ? 1 : 0);                     // FIFO status
 
       for (size_t pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
         JBit(pmt).set(status, getHighRateVeto(JPMTIdentifier(id,pmt)) ? 1 : 0) ;            // high-rate veto
       }
 
       return JDAQFrameStatus(daq, status, fifo);
     }

virtual void JDETECTOR::JCLBDefaultSimulatorInterface::processData	(	const JModuleIdentifier &	id,
		const JCLBInput &	input,
		JDAQSuperFrame &	output
	)		const

inlineoverridevirtual

Process data.

Parameters

id	module identifier
input	PMT data
output	CLB data

Implements JDETECTOR::JCLBSimulator.

Definition at line 166 of file JCLBDefaultSimulatorInterface.hh.

     {      
       using namespace std;
       using namespace KM3NETDAQ;
 
 
       const double Tmin  =  0.0;                // Minimal TDC value [ns]
       const double Tmax  =  getFrameTime();     // Maximal TDC value [ns]
 
       output.setDAQFrameStatus(getDAQFrameStatus(id));
 
       buffer.clear();
 
       size_t ns = 0;                            // current number of hits in data frame
 
       for (size_t pmt = 0; pmt != input.size(); ++pmt) {
 
         const JPMTData<JPMTPulse>& in = input[pmt];
 
         // TDC
 
         for (JPMTData<JPMTPulse>::const_iterator hit = in.begin(); hit != in.end(); ++hit) {
 
           if (hit->t_ns >= Tmin && hit->t_ns <= Tmax) {
 
             double t1  = hit->t_ns;
             double tot = hit->tot_ns;
 
             // generate multiple hits if time-over-threshold exceeds dynamic range
           
             while (tot > JDAQHit::getMaximalToT()) {
 
               buffer.push_back(TDC->makeHit((JPMT_t) pmt, t1, (JTOT_t) JDAQHit::getMaximalToT()));
 
               t1  += JDAQHit::getMaximalToT();
               tot -= JDAQHit::getMaximalToT();
             }
 
             if (tot > getMinimalToT()) {
               buffer.push_back(TDC->makeHit((JPMT_t) pmt, t1, (JTOT_t) (tot + 0.5)));
             }
           }
         }
 
         if (buffer.size() > ns + getMaximalNumberOfHits()) {
 
           output.setHighRateVeto(pmt, true);
           
           buffer.resize(ns + getMaximalNumberOfHits());
         }
 
         inplace_merge(buffer.begin(), buffer.begin() + ns, buffer.end());
 
         ns = buffer.size();
       }
 
 
       // simulate UDP packet loss
 
       const int n1  =  getUDPNumberOfReceivedPackets(id);
       const int n2  =  getUDPMaximalSequenceNumber  (id);
 
       if (n1 < n2 + 1) {
 
         JTDC_t t0 = (JTDC_t)   Tmin;
         JTDC_t ts = (JTDC_t) ((Tmax - Tmin) / (n2 + 1));
 
         for (int i = 0; i != n2 + 1; ++i, t0 += ts) {
           
           if (gRandom->Rndm() * (n2 + 1) >= (double) n1) {
 
             std::vector<JDAQHit>::iterator p = lower_bound(buffer.begin(), buffer.end(), t0,      compare);
             std::vector<JDAQHit>::iterator q = lower_bound(buffer.begin(), buffer.end(), t0 + ts, compare);
           
             buffer.erase(p,q);
           }
         }
       }
 
 
       // process data through state machine
 
       if (buffer.size() > 1) {
 
         for (std::vector<JDAQHit>::iterator q = buffer.begin(), p = q++; q != buffer.end(); ++p, ++q) {
 
           if (state_machine->maybeSwapped(*p,*q)) {
             iter_swap(p,q);
           }
         }
       }
 
 
       // store data
 
       output.add(buffer.size(), buffer.data());
     }