Timeslice with random data. More...

#include <JRandomTimeslice.hh>

Inheritance diagram for KM3NETDAQ::JRandomTimeslice:

Public Member Functions
	JRandomTimeslice ()
	Default constructor. More...

	JRandomTimeslice (const JDAQChronometer &chronometer, const JDetectorSimulator &simbad)
	Constructor. More...

void	recycle (const double T_ns)
	Recycle time slice by randomly shuffling time intervals of data. More...

	ClassDef (JDAQTimesliceL0, 1)

	ClassDef (JDAQTimeslice, 4)

	ClassDef (JDAQPreamble, 1)

	ClassDef (JDAQTimesliceHeader, 2)

	ClassDef (JDAQHeader, 2)

	ClassDef (JDAQChronometer, 3)

void	clear ()
	Clear data. More...

JDAQTimeslice &	add (const JDAQTimeslice &timeslice)
	Add another timeslice. More...

std::ostream &	print (std::ostream &out, const bool lpr=false) const
	Print DAQ Timeslice. More...

int	getLength () const
	Get length. More...

int	getDataType () const
	Get data type. More...

	ClassDefNV (JDAQAbstractPreamble, 1)

const JDAQTimesliceHeader &	getDAQTimesliceHeader () const
	Get DAQ time slice header. More...

const JDAQHeader &	getDAQHeader () const
	Get DAQ header. More...

void	setDAQHeader (const JDAQHeader &header)
	Set DAQ header. More...

const JDAQChronometer &	getDAQChronometer () const
	Get DAQ chronometer. More...

void	setDAQChronometer (const JDAQChronometer &chronometer)
	Set DAQ chronometer. More...

int	getDetectorID () const
	Get detector identifier. More...

int	getRunNumber () const
	Get run number. More...

int	getFrameIndex () const
	Get frame index. More...

JDAQUTCExtended	getTimesliceStart () const
	Get start of timeslice. More...

void	setRunNumber (const int run)
	Set run number. More...

void	setFrameIndex (const int frame_index)
	Set frame index. More...

void	setTimesliceStart (const JDAQUTCExtended &timeslice_start)
	Set timeslice start time. More...

Static Public Member Functions
template<class T >
static JDAQPreamble	getDAQPreamble (const T &object)
	Get DAQ preamble. More...

Protected Attributes
int	length

int	type

int	detector_id

int	run

int	frame_index

JDAQUTCExtended	timeslice_start

Detailed Description

Timeslice with random data.

Definition at line 29 of file JRandomTimeslice.hh.

Constructor & Destructor Documentation

KM3NETDAQ::JRandomTimeslice::JRandomTimeslice ( )

inline

Default constructor.

Definition at line 35 of file JRandomTimeslice.hh.

36 {}

KM3NETDAQ::JRandomTimeslice::JRandomTimeslice	(	const JDAQChronometer &	chronometer,
		const JDetectorSimulator &	simbad
	)

inline

Constructor.

Parameters

chronometer	chronometer
simbad	detector simulator

Definition at line 45 of file JRandomTimeslice.hh.

     {
       using namespace JPP;
 
       setDAQChronometer(chronometer);
 
       if (simbad.hasK40Simulator() &&
           simbad.hasPMTSimulator() &&
           simbad.hasCLBSimulator()) {
 
         const double     Tmin = getTimeSinceRTS(chronometer.getFrameIndex());   // [ns]
 
         const JTimeRange period(Tmin, Tmin + getFrameTime());                   // [ns]
         
         JModuleData buffer;
         
         for (JDetector::const_iterator module = simbad->begin(); module != simbad->end(); ++module) {
 
           if (!module->empty()) {
 
             buffer.reset(module->size());
 
             simbad.generateHits(*module, period, buffer);
 
             this->push_back(JDAQSuperFrame(JDAQSuperFrameHeader(chronometer, module->getID())));
 
             simbad(*module, buffer, *(this->rbegin()));
           }
         }
       }
     }

Member Function Documentation

void KM3NETDAQ::JRandomTimeslice::recycle ( const double T_ns )

inline

Recycle time slice by randomly shuffling time intervals of data.

Hits within one time interval are swapped with hits within another -randomly selected- time interval.

Note that the time interval should be:

(much) larger than time differences of hits in L1 coincidence; and
(much) smaller than time covered by data (as per KM3NETDAQ::getFrameTime()).

Parameters

T_ns	time interval

Definition at line 90 of file JRandomTimeslice.hh.

     {
       using namespace std;
       using namespace JPP;
 
       typedef JDAQHit::JTDC_t  JTDC_t;
 
 
       size_t N;                                                  // number of time intervals
 
       if      (T_ns < 1.0)
         N = (size_t)  getFrameTime();
       else if (T_ns < getFrameTime())
         N = (size_t) (getFrameTime() / T_ns + 0.5);
       else
         N = 1;
 
       const JTDC_t T_step = (JTDC_t) (getFrameTime() / N);       // TDC interval
 
 
       typedef vector<JDAQHit>  buffer_type;
 
       vector<size_t>      index (N);                             // indices of time intervals for swapping
       vector<buffer_type> buffer(N);                             // data per time interval
 
       for (iterator frame = this->begin(); frame != this->end(); ++frame) {
 
         if (!frame->empty()) {
 
 
           for (size_t i = 0; i != N; ++i) {
             index [i] = i;
             buffer[i].clear();
           }
 
 
           vector<size_t> N_max(NUMBER_OF_PMTS, N);               // number of time intervals per PMT
           vector<JTDC_t> T_max(NUMBER_OF_PMTS, 0);               // maximal time per PMT, e.g.\ due to high-rate veto
 
           for (JDAQSuperFrame::const_iterator hit = frame->begin(); hit != frame->end(); ++hit) {
             T_max[hit->getPMT()] = hit->getT();
           }
 
           for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
 
            const size_t N_i = (T_max[i] + T_step) / T_step;
 
             if (N_i < N_max[i]) {
               N_max[i] = N_i;
             }
           }
 
           sort(N_max.begin(), N_max.end());                      // swaps by allowed time range
 
           for (size_t i = 0, L = 0; i != N_max.size(); ++i) {
 
             size_t M = N_max[i];
 
             if (L != M) {
               for (size_t i = M - L - 1; i != 0; --i) {
                 std::swap(index[L + i], index[L + gRandom->Integer(i+1)]);
               }
             }
 
             L = M;
           }
 
           const size_t M = N_max[NUMBER_OF_PMTS - 1];
 
 
           // Wall street shuflle
 
           JDAQSuperFrame::const_iterator hit = frame->begin();
 
           for (size_t in = 0; in != M; ++in) {
 
             const size_t  out   = index[in];
             buffer_type&  zbuf  = buffer[out];
 
             const JTDC_t  T_in  = in  * T_step;
             const JTDC_t  T_out = out * T_step;
 
             for ( ; hit != frame->end() && hit->getT() < T_in + T_step; ++hit) {
               zbuf.push_back(JDAQHit(hit->getPMT(), (hit->getT() - T_in) + T_out, hit->getToT()));
             }
           }
 
 
           // copy back data
 
           for (size_t i = 0, number_of_hits = 0; i != M; ++i) {
 
             memcpy(frame->data() + number_of_hits, buffer[i].data(), buffer[i].size() * sizeof(JDAQHit)); 
 
             number_of_hits += buffer[i].size();
           }
         }
       }
     }

KM3NETDAQ::JDAQTimesliceL0::ClassDef	(	JDAQTimesliceL0	,
		1
	)

inherited

KM3NETDAQ::JDAQTimeslice::ClassDef	(	JDAQTimeslice	,
		4
	)

inherited

KM3NETDAQ::JDAQPreamble::ClassDef	(	JDAQPreamble	,
		1
	)

inherited

KM3NETDAQ::JDAQTimesliceHeader::ClassDef	(	JDAQTimesliceHeader	,
		2
	)

inherited

KM3NETDAQ::JDAQHeader::ClassDef	(	JDAQHeader	,
		2
	)

inherited

KM3NETDAQ::JDAQChronometer::ClassDef	(	JDAQChronometer	,
		3
	)

inherited

void KM3NETDAQ::JDAQTimeslice::clear ( )

inlineinherited

Clear data.

Definition at line 97 of file JDAQTimeslice.hh.

     {
       for (iterator i = this->begin(); i != this->end(); ++i) {
         i->clear();
       }
 
       std::vector<JDAQSuperFrame>::clear();
     }

JDAQTimeslice& KM3NETDAQ::JDAQTimeslice::add ( const JDAQTimeslice & timeslice )

inlineinherited

Add another timeslice.

Parameters

timeslice timeslice

Returns: this timeslice

Definition at line 133 of file JDAQTimeslice.hh.

     {
       using namespace std;
 
       map<JDAQModuleIdentifier, int> buffer;
 
       for (const_iterator i = this->begin(); i != this->end(); ++i) {
         buffer[i->getModuleIdentifier()] = distance(static_cast<const JDAQTimeslice&>(*this).begin(),i);
       }
 
       for (JDAQTimeslice::const_iterator i = timeslice.begin(); i != timeslice.end(); ++i) {
 
         map<JDAQModuleIdentifier, int>::const_iterator p = buffer.find(i->getModuleIdentifier());
         
         if (p != buffer.end()) {
 
           JDAQSuperFrame& frame = this->at(p->second);
 
           frame.add(*i);
 
           sort(frame.begin(), frame.end());
 
         } else {
 
           this->push_back(*i);
         }
       }
       
       return *this;
     }

std::ostream& KM3NETDAQ::JDAQTimeslice::print	(	std::ostream &	out,
		const bool	lpr = `false`
	)		const

inlineinherited

Print DAQ Timeslice.

Parameters

out	output stream
lpr	long print

Returns: output stream

Definition at line 172 of file JDAQTimeslice.hh.

     {
       using namespace std;
       
       out << this->ClassName() << endl;
       out << dynamic_cast<const JDAQPreamble&>   (*this) << endl;
       out << dynamic_cast<const JDAQChronometer&>(*this) << endl;
       
       for (JDAQTimeslice::const_iterator frame = this->begin(); frame != this->end(); ++frame) {
         
         out << ' ' << setw(10) << frame->getModuleID();
         out << ' ' << setw(6)  << frame->getLength();
         out << ' ' << setw(6)  << frame->getDataType();
         out << ' ' << setw(6)  << frame->getTimesliceStart();
         out << ' ' << setw(8)  << setfill('0') << hex << frame->getStatus()     << dec << setfill(' ');
         out << '|' << setw(8)  << setfill('0') << hex << frame->getFIFOStatus() << dec << setfill(' ');
         out << ' ' << setw(2)  << frame->getUDPNumberOfReceivedPackets();
         out << '/' << setw(2)  << frame->getUDPMaximalSequenceNumber();
         out << ' ' << setw(6)  << frame->size();
         
         if (!lpr) {
           
           if (!frame->empty()) {
             
             out << ' ' << setw(10) << frame-> begin()->getT();
             out << " ... ";
             out << ' ' << setw(10) << frame->rbegin()->getT();
           }
           
           out << endl;
           
         } else {
           
           out << endl;
           
           int n = 1;
           
           for (JDAQFrame::const_iterator hit = frame->begin(); hit != frame->end(); ++hit, ++n) {
             out << setw(2) << (int) hit->getPMT() << ' ' 
                 << setw(8) << (int) hit->getT()   << ' ' 
                 << setw(3) << (int) hit->getToT() << (n%10 == 0 ? '\n' : ' ');
           }
 
           out << endl;
         }
       }
 
       return out;
     }

template<class T >

static JDAQPreamble KM3NETDAQ::JDAQPreamble::getDAQPreamble ( const T & object )

inlinestaticinherited

Get DAQ preamble.

This method should be used for binary I/O to get the actual data for the given object.
To this end, the following method should be overloaded for the corresponding data type.

   size_t  getSizeof(const T&);

Parameters

object object

Returns: preamble

Definition at line 76 of file JDAQPreamble.hh.

     {
       static JDAQPreamble preamble;
 
       preamble.length = getSizeof(object);
       preamble.type   = KM3NETDAQ::getDataType<T>();
 
       return preamble;
     }

int KM3NETDAQ::JDAQAbstractPreamble::getLength ( ) const

inlineinherited

Get length.

Returns: number of bytes

Definition at line 49 of file JDAQAbstractPreamble.hh.

     { 
       return length;
     }

int KM3NETDAQ::JDAQAbstractPreamble::getDataType ( ) const

inlineinherited

Get data type.

Returns: data type

Definition at line 60 of file JDAQAbstractPreamble.hh.

     { 
       return type; 
     }

KM3NETDAQ::JDAQAbstractPreamble::ClassDefNV	(	JDAQAbstractPreamble	,
		1
	)

inherited

const JDAQTimesliceHeader& KM3NETDAQ::JDAQTimesliceHeader::getDAQTimesliceHeader ( ) const

inlineinherited

Get DAQ time slice header.

Returns: DAQ time slice header

Definition at line 43 of file JDAQTimesliceHeader.hh.

     {
       return static_cast<const JDAQTimesliceHeader&>(*this);
     }

const JDAQHeader& KM3NETDAQ::JDAQHeader::getDAQHeader ( ) const

inlineinherited

Get DAQ header.

Returns: DAQ header

Definition at line 49 of file JDAQHeader.hh.

     { 
       return static_cast<const JDAQHeader&>(*this); 
     }

void KM3NETDAQ::JDAQHeader::setDAQHeader ( const JDAQHeader & header )

inlineinherited

Set DAQ header.

Parameters

header DAQ header

Definition at line 60 of file JDAQHeader.hh.

     { 
       static_cast<JDAQHeader&>(*this) = header;
     }

const JDAQChronometer& KM3NETDAQ::JDAQChronometer::getDAQChronometer ( ) const

inlineinherited

Get DAQ chronometer.

Returns: DAQ chronometer

Definition at line 88 of file JDAQChronometer.hh.

     { 
       return static_cast<const JDAQChronometer&>(*this); 
     }

void KM3NETDAQ::JDAQChronometer::setDAQChronometer ( const JDAQChronometer & chronometer )

inlineinherited

Set DAQ chronometer.

Parameters

chronometer DAQ chronometer

Definition at line 99 of file JDAQChronometer.hh.

     { 
       static_cast<JDAQChronometer&>(*this) = chronometer;
     }

int KM3NETDAQ::JDAQChronometer::getDetectorID ( ) const

inlineinherited

Get detector identifier.

Returns: detector identifier

Definition at line 110 of file JDAQChronometer.hh.

     { 
       return detector_id; 
     }

int KM3NETDAQ::JDAQChronometer::getRunNumber ( ) const

inlineinherited

Get run number.

Returns: run number

Definition at line 121 of file JDAQChronometer.hh.

     { 
       return run; 
     }

int KM3NETDAQ::JDAQChronometer::getFrameIndex ( ) const

inlineinherited

Get frame index.

Returns: frame index

Definition at line 132 of file JDAQChronometer.hh.

     { 
       return frame_index; 
     }

JDAQUTCExtended KM3NETDAQ::JDAQChronometer::getTimesliceStart ( ) const

inlineinherited

Get start of timeslice.

Returns: timeslice start

Definition at line 144 of file JDAQChronometer.hh.

     {
       return timeslice_start;
     }

void KM3NETDAQ::JDAQChronometer::setRunNumber ( const int run )

inlineinherited

Set run number.

Parameters

run	run number

Definition at line 155 of file JDAQChronometer.hh.

     { 
       this->run = run;
     }

void KM3NETDAQ::JDAQChronometer::setFrameIndex ( const int frame_index )

inlineinherited

Set frame index.

Parameters

frame_index frame index

Definition at line 166 of file JDAQChronometer.hh.

     { 
       this->frame_index = frame_index;
     }

void KM3NETDAQ::JDAQChronometer::setTimesliceStart ( const JDAQUTCExtended & timeslice_start )

inlineinherited

Set timeslice start time.

Parameters

timeslice_start timeslice start time

Definition at line 177 of file JDAQChronometer.hh.

     {
       this->timeslice_start = timeslice_start;
     }

Member Data Documentation

int KM3NETDAQ::JDAQAbstractPreamble::length

protectedinherited

Definition at line 69 of file JDAQAbstractPreamble.hh.

int KM3NETDAQ::JDAQAbstractPreamble::type

protectedinherited

Definition at line 70 of file JDAQAbstractPreamble.hh.

int KM3NETDAQ::JDAQChronometer::detector_id

protectedinherited

Definition at line 187 of file JDAQChronometer.hh.

int KM3NETDAQ::JDAQChronometer::run

protectedinherited

Definition at line 188 of file JDAQChronometer.hh.

int KM3NETDAQ::JDAQChronometer::frame_index

protectedinherited

Definition at line 189 of file JDAQChronometer.hh.

JDAQUTCExtended KM3NETDAQ::JDAQChronometer::timeslice_start

protectedinherited

Definition at line 190 of file JDAQChronometer.hh.

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

JRandomTimeslice.hh

Public Member Functions

Static Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation