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 Member Functions
	JCLBDefaultSimulator (const JDAQHit::JTDC_t Tmax=0)
	Constructor. More...

virtual void	processData (const JCLBInput &input, JDAQFrame &output) const
	Process data. More...

Static Public Member Functions
static JTOT_t	getMaximalToT ()
	Get maximal pulse length of time-over-threshold measurement. More...

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

Protected Types
typedef JDAQHit::JPMT_t	JPMT_t

typedef JDAQHit::JTDC_t	JTDC_t

typedef JDAQHit::JTOT_t	JTOT_t

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 also implemented in this class.

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

Definition at line 139 of file JCLBDefaultSimulator.hh.

Member Typedef Documentation

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

protectedinherited

Definition at line 45 of file JCLBDefaultSimulatorInterface.hh.

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

protectedinherited

Definition at line 46 of file JCLBDefaultSimulatorInterface.hh.

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

protectedinherited

Definition at line 47 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 228 of file JCLBDefaultSimulator.hh.

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

Member Function Documentation

virtual void JDETECTOR::JCLBDefaultSimulatorInterface::processData	(	const JCLBInput &	input,
		JDAQFrame &	output
	)		const

inlinevirtualinherited

Process data.

Parameters

input	PMT data
output	CLB data

Implements JDETECTOR::JCLBSimulator.

Definition at line 126 of file JCLBDefaultSimulatorInterface.hh.

     {      
       using namespace std;
 
       const JTDC_t TMIN = 0;                   // [ns]
       const JTDC_t TMAX = (JTDC_t) getRTS();   // [ns]
 
       buffer.clear();
 
       for (JPMT_t pmt = 0; pmt != (JPMT_t) input.size(); ++pmt) {
 
         const JPMTData<JPMTPulse>& in = input[pmt];
 
         // TDC
 
         JPMTData<JPMTPulse>::const_iterator hit = in.begin(); 
 
         for ( ; hit != in.end() && hit->t_ns < TMIN; ++hit) {}
 
         for ( ; hit != in.end() && hit->t_ns < TMAX; ++hit) {
 
           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(pmt, t1, JDAQHit::getMaximalToT()));
             
             t1  += JDAQHit::getMaximalToT();
             tot -= JDAQHit::getMaximalToT();
           }
           
           if (tot > getMinimalToT()) {
             buffer.push_back(TDC->makeHit(pmt, t1, (JTOT_t) (tot + 0.5)));
           }
         }
       }
 
       if (buffer.size() > 1) {
 
         // process data through state machine
 
         sort(buffer.begin(), buffer.end());
 
         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);
           }
         }
       }
 
       output.add(buffer.size(), buffer.data());
     }