CLB simulation based on run-by-run information. More...

#include <JCLBRunByRunSimulator.hh>

Inheritance diagram for JTRIGGER::JCLBRunByRunSimulator:

Public Types
typedef JDAQHit::JPMT_t	JPMT_t

typedef JDAQHit::JTDC_t	JTDC_t

typedef JDAQHit::JTOT_t	JTOT_t

Public Member Functions
	JCLBRunByRunSimulator (const JSummaryRouter &router)
	Constructor.

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

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

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

virtual bool	setHighRateVeto () const
	Set high-rate veto based on number of hits.

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

virtual bool	getFIFOStatus (const JPMTIdentifier &id) const override
	Get FIFO (almost) full of given PMT.

virtual bool	hasCLB (const JModuleIdentifier &id) const
	Check if CLB exist.

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

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

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

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

Protected Attributes
const JSummaryRouter &	summary_router

Private Attributes
std::unique_ptr< JTDC >	TDC

std::unique_ptr< JStateMachine >	state_machine

Detailed Description

CLB simulation based on run-by-run information.

This class overwrites the method JCLBDefaultSimulator::getHighRateVeto. The actual high-rate veto is obtained from summary data.

Definition at line 32 of file JCLBRunByRunSimulator.hh.

Member Typedef Documentation

◆ JPMT_t

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

inherited

Definition at line 62 of file JCLBDefaultSimulatorInterface.hh.

◆ JTDC_t

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

inherited

Definition at line 63 of file JCLBDefaultSimulatorInterface.hh.

◆ JTOT_t

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

inherited

Definition at line 64 of file JCLBDefaultSimulatorInterface.hh.

Constructor & Destructor Documentation

◆ JCLBRunByRunSimulator()

JTRIGGER::JCLBRunByRunSimulator::JCLBRunByRunSimulator ( const JSummaryRouter & router )

inline

Constructor.

Parameters

router summary router

Definition at line 41 of file JCLBRunByRunSimulator.hh.

                                                        :
      JCLBDefaultSimulator(),
      summary_router(router)
    {}

Member Function Documentation

◆ getUDPNumberOfReceivedPackets()

virtual int JTRIGGER::JCLBRunByRunSimulator::getUDPNumberOfReceivedPackets ( const JModuleIdentifier & id ) const

inlineoverridevirtual

Get number of received UDP packets.

Parameters

id	module identifier

Returns: number of received UDP packets

Reimplemented from JDETECTOR::JCLBDefaultSimulatorInterface.

Definition at line 53 of file JCLBRunByRunSimulator.hh.

    {
      using namespace KM3NETDAQ;
 
      if (summary_router.hasSummaryFrame(id.getID())) {
 
        const JDAQSummaryFrame& frame = summary_router.getSummaryFrame(id.getID());
 
        return frame.getUDPNumberOfReceivedPackets();
      }
 
      return JCLBDefaultSimulator::getUDPNumberOfReceivedPackets(id);
    }

◆ getUDPMaximalSequenceNumber()

virtual int JTRIGGER::JCLBRunByRunSimulator::getUDPMaximalSequenceNumber ( const JModuleIdentifier & id ) const

inlineoverridevirtual

Get maximal sequence number of UDP packet.

Parameters

id	module identifier

Returns: maximal sequence number

Reimplemented from JDETECTOR::JCLBDefaultSimulatorInterface.

Definition at line 74 of file JCLBRunByRunSimulator.hh.

    {
      using namespace KM3NETDAQ;
 
      if (summary_router.hasSummaryFrame(id.getID())) {
 
        const JDAQSummaryFrame& frame = summary_router.getSummaryFrame(id.getID());
 
        return frame.getUDPMaximalSequenceNumber();
      }
 
      return JCLBDefaultSimulator::getUDPMaximalSequenceNumber(id);
    }

◆ hasUDPTrailer()

virtual bool JTRIGGER::JCLBRunByRunSimulator::hasUDPTrailer ( const JModuleIdentifier & id ) const

inlineoverridevirtual

Get UDP trailer status.

Parameters

id	module identifier

Returns: true

Reimplemented from JDETECTOR::JCLBDefaultSimulatorInterface.

Definition at line 95 of file JCLBRunByRunSimulator.hh.

    {
      using namespace KM3NETDAQ;
 
      if (summary_router.hasSummaryFrame(id.getID())) {
 
        const JDAQSummaryFrame& frame = summary_router.getSummaryFrame(id.getID());
 
        return frame.hasUDPTrailer();
      }
 
      return JCLBDefaultSimulator::hasUDPTrailer(id);
    }

◆ setHighRateVeto()

virtual bool JTRIGGER::JCLBRunByRunSimulator::setHighRateVeto ( ) const

inlinevirtual

Set high-rate veto based on number of hits.

Returns: false

Reimplemented from JDETECTOR::JCLBDefaultSimulatorInterface.

Definition at line 115 of file JCLBRunByRunSimulator.hh.

    {
      return false;
    }

◆ getHighRateVeto()

virtual bool JTRIGGER::JCLBRunByRunSimulator::getHighRateVeto ( const JPMTIdentifier & id ) const

inlineoverridevirtual

Get high-rate veto of given PMT.

Parameters

id	PMT identifier

Returns: true if high-rate veto; else false

Reimplemented from JDETECTOR::JCLBDefaultSimulatorInterface.

Definition at line 127 of file JCLBRunByRunSimulator.hh.

    {
      using namespace KM3NETDAQ;
 
      if (summary_router.hasSummaryFrame(id.getModuleID())) {
 
        const JDAQSummaryFrame& frame = summary_router.getSummaryFrame(id.getModuleID());
 
        return frame.testHighRateVeto(id.getPMTAddress());
      }
 
      return JCLBDefaultSimulator::getHighRateVeto(id);
    }

◆ getFIFOStatus()

virtual bool JTRIGGER::JCLBRunByRunSimulator::getFIFOStatus ( const JPMTIdentifier & id ) const

inlineoverridevirtual

Get FIFO (almost) full of given PMT.

Parameters

id	PMT identifier

Returns: true if FIFO (almost) full; else false

Reimplemented from JDETECTOR::JCLBDefaultSimulatorInterface.

Definition at line 148 of file JCLBRunByRunSimulator.hh.

    {
      using namespace KM3NETDAQ;
 
      if (summary_router.hasSummaryFrame(id.getModuleID())) {
 
        const JDAQSummaryFrame& frame = summary_router.getSummaryFrame(id.getModuleID());
 
        return frame.testFIFOStatus(id.getPMTAddress());
      }
 
      return JCLBDefaultSimulator::getFIFOStatus(id);
    }

◆ hasCLB()

virtual bool JTRIGGER::JCLBRunByRunSimulator::hasCLB ( const JModuleIdentifier & id ) const

inlinevirtual

Check if CLB exist.

Parameters

id	module identifier

Returns: true if summary data available; else false

Reimplemented from JDETECTOR::JCLBSimulator.

Definition at line 169 of file JCLBRunByRunSimulator.hh.

    {
      return summary_router.hasSummaryFrame(id.getID());
    }

◆ getDAQFrameStatus()

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 144 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)));                     // high-rate veto
        JBit(pmt).set(fifo,   getFIFOStatus  (JPMTIdentifier(id,pmt)));                     // FIFO (almost) full
      }
 
      return JDAQFrameStatus(daq, status, fifo);
    }

◆ processData()

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 169 of file JCLBDefaultSimulatorInterface.hh.

    {      
      using namespace std;
      using namespace JPP;
      using namespace KM3NETDAQ;
 
 
      const double Tmin  =  0.0;                // Minimal TDC value [ns]
      const double Tmax  =  getFrameTime();     // Maximal TDC value [ns]
 
      output.setDAQFrameStatus(getDAQFrameStatus(id));
 
      vector<JDAQHit> buffer;
 
      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()) {
 
          if (setHighRateVeto()) {
            output.setHighRateVeto(pmt, true);
          }
 
          if (output.testHighRateVeto(pmt)) {
            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) {
 
        const int ns = getBayesianMedian(n2, n1, 0);
 
        JTDC_t t0 = (JTDC_t)   Tmin;
        JTDC_t ts = (JTDC_t) ((Tmax - Tmin) / ns);
 
        for (int i = 0; i != ns; ++i, t0 += ts) {
 
          if (gRandom->Rndm() * ns >= (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());
    }

◆ getMinimalToT()

static double JDETECTOR::JCLBDefaultSimulatorInterface::getMinimalToT ( )

inlinestaticinherited

Get minimal pulse length of time-over-threshold measurement.

Returns: TDC value [ns]

Definition at line 351 of file JCLBDefaultSimulatorInterface.hh.

    {
      return 0.5;
    }

Member Data Documentation

◆ summary_router

const JSummaryRouter& JTRIGGER::JCLBRunByRunSimulator::summary_router

protected

Definition at line 176 of file JCLBRunByRunSimulator.hh.

◆ compare

const struct JDETECTOR::JCLBDefaultSimulatorInterface::compare JDETECTOR::JCLBDefaultSimulatorInterface::compare

staticinherited

◆ TDC

std::unique_ptr<JTDC> JDETECTOR::JCLBDefaultSimulatorInterface::TDC

privateinherited

Definition at line 389 of file JCLBDefaultSimulatorInterface.hh.

◆ state_machine

std::unique_ptr<JStateMachine> JDETECTOR::JCLBDefaultSimulatorInterface::state_machine

privateinherited

Definition at line 390 of file JCLBDefaultSimulatorInterface.hh.

The documentation for this class was generated from the following file:

JCLBRunByRunSimulator.hh

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Attributes

Private Attributes

Detailed Description

Member Typedef Documentation

◆ JPMT_t

◆ JTDC_t

◆ JTOT_t

Constructor & Destructor Documentation

◆ JCLBRunByRunSimulator()

Member Function Documentation

◆ getUDPNumberOfReceivedPackets()

◆ getUDPMaximalSequenceNumber()

◆ hasUDPTrailer()

◆ setHighRateVeto()

◆ getHighRateVeto()

◆ getFIFOStatus()

◆ hasCLB()

◆ getDAQFrameStatus()

◆ processData()

◆ getMinimalToT()

Member Data Documentation

◆ summary_router

◆ compare

◆ TDC

◆ state_machine