Default CLB simulation. More...

#include <JCLBDefaultSimulator.hh>

Inheritance diagram for JDETECTOR::JCLBDefaultSimulator:

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

class	JTDC
	Implementation of non-linearity of 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
	JCLBDefaultSimulator (const JDAQHit::JTDC_t Tmax=0)
	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

Detailed Description

Default CLB simulation.

This class provides for the implementation of the conversion of a time-over-threshold pulse to a so-called hit.
A hit consists of:

PMT readout channel;
time stamp of the leading edge of the time-over-threshold pulse;
length of the time-over-threshold pulse;

The time stamping and pulse length determination is done by a TDC which is integrated inside the FPGA. The observed non-linearity of the TDC is implemented in this class.

Although the hits from the same PMT are striclty time sorted, hits from different PMTs could be reordered due to the occupancy of the buffers inside the FPGA. This reordering is also implemented in the this class.

Definition at line 135 of file JCLBDefaultSimulator.hh.

Member Typedef Documentation

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

inherited

Definition at line 60 of file JCLBDefaultSimulatorInterface.hh.

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

inherited

Definition at line 61 of file JCLBDefaultSimulatorInterface.hh.

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

inherited

Definition at line 62 of file JCLBDefaultSimulatorInterface.hh.

Constructor & Destructor Documentation

JDETECTOR::JCLBDefaultSimulator::JCLBDefaultSimulator ( const JDAQHit::JTDC_t Tmax = 0 )

inline

Constructor.

Parameters

Tmax	maximal time difference for swap [ns]

Definition at line 224 of file JCLBDefaultSimulator.hh.

                                                        : 
       JCLBDefaultSimulatorInterface(new JTDC(), 
                                     new JStateMachine(Tmax))
     {}

Member Function Documentation

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

inlineinherited

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

inlineoverridevirtualinherited

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());
     }