Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
KM3NETDAQ::JEventTimeslice Class Reference

Timeslice with Monte Carlo event. More...

#include <JEventTimeslice.hh>

Inheritance diagram for KM3NETDAQ::JEventTimeslice:
KM3NETDAQ::JDAQTimeslice KM3NETDAQ::JDAQPreamble KM3NETDAQ::JDAQTimesliceHeader std::vector< JDAQSuperFrame > KM3NETDAQ::JDAQAbstractPreamble TObject KM3NETDAQ::JDAQHeader KM3NETDAQ::JDAQChronometer

Public Member Functions

 JEventTimeslice (const JDAQChronometer &chronometer, const JDetectorSimulator &simbad, const Evt &event, const JTimeRange &period=JTimeRange::DEFAULT_RANGE)
 Constructor. More...
 
void clear ()
 Clear data. More...
 
JDAQTimesliceadd (const JDAQTimeslice &timeslice)
 Add another timeslice. More...
 
std::ostream & print (std::ostream &out, const bool lpr=false) const
 Print DAQ Timeslice. More...
 
virtual int getSize () const
 Get size of object. More...
 
 ClassDef (JDAQTimeslice, 4)
 
 ClassDef (JDAQPreamble, 1)
 
 ClassDef (JDAQTimesliceHeader, 2)
 
 ClassDef (JDAQHeader, 2)
 
 ClassDef (JDAQChronometer, 3)
 
int getLength () const
 Get length. More...
 
int getDataType () const
 Get data type. More...
 
 ClassDefNV (JDAQAbstractPreamble, 1)
 
const JDAQTimesliceHeadergetDAQTimesliceHeader () const
 Get DAQ time slice header. More...
 
const JDAQHeadergetDAQHeader () const
 Get DAQ header. More...
 
void setDAQHeader (const JDAQHeader &header)
 Set DAQ header. More...
 
const JDAQChronometergetDAQChronometer () 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...
 
const JDAQUTCExtendedgetTimesliceStart () 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

static int sizeOf ()
 Get size of object. More...
 
static int sizeOf ()
 Get size of object. More...
 

Protected Member Functions

void setLength ()
 Set length. More...
 

Protected Attributes

int length
 
int type
 
int detector_id
 
int run
 
int frame_index
 
JDAQUTCExtended timeslice_start
 

Detailed Description

Timeslice with Monte Carlo event.

Definition at line 34 of file JEventTimeslice.hh.

Constructor & Destructor Documentation

KM3NETDAQ::JEventTimeslice::JEventTimeslice ( const JDAQChronometer chronometer,
const JDetectorSimulator simbad,
const Evt &  event,
const JTimeRange period = JTimeRange::DEFAULT_RANGE 
)
inline

Constructor.

This constructor converts the Monte Carlo event to a DAQ timeslice. In this, both the PMT and CLB simulations are applied. If the simulator is available and the given time window is valid, the K40 background is generated and added to the timeslice data. N.B. The given time window extends the time range of hits in the event, if any.

Parameters
chronometerchronometer
simbaddetector simulator
eventMonte Carlo event
periodtime window [ns]

Definition at line 52 of file JEventTimeslice.hh.

55  :
57  {
58  using namespace std;
59  using namespace JDETECTOR;
60  using namespace JAANET;
62 
63  setDAQChronometer(chronometer);
64 
65 
66  if (simbad.hasPMTSimulator() &&
67  simbad.hasCLBSimulator()) {
68 
69  const JTimeConverter converter(event, chronometer);
70 
71  const JTimeRange timeRange = (period.is_valid() ? getTimeRange(event, period) : getTimeRange(event));
72 
73  typedef map<int, JModuleData> JMap_t; // map module index to data
74 
75  JMap_t data;
76 
77  for (vector<Hit>::const_iterator hit = event.mc_hits.begin(); hit != event.mc_hits.end(); ++hit) {
78 
79  if (simbad.hasPMT(hit->pmt_id)) {
80 
81  const JPMTAddress& address = simbad.getAddress(hit->pmt_id);
82 
83  if (!period.is_valid() || timeRange(getTime(*hit))) {
84 
85  data[address.first].resize(NUMBER_OF_PMTS);
86 
87  data[address.first][address.second].push_back(JPMTSignal(converter.putTime(getTime(*hit)), (int) getNPE(*hit)));
88  }
89  }
90  }
91 
92 
93  if (simbad.hasK40Simulator() && period.is_valid()) {
94 
95 
97 
98  for (JMap_t::const_iterator i = data.begin(); i != data.end(); ++i) {
99  time_range.combine(getTimeRange(i->second));
100  }
101 
102  time_range.add(period);
103 
104 
105  JModuleData buffer;
106 
107  for (JDetector::const_iterator module = simbad->begin(); module != simbad->end(); ++module) {
108 
109 
110  // signal
111 
112  JMap_t::iterator i = data.find(distance(simbad->begin(),module));
113 
114  if (i != data.end())
115  buffer.swap(i->second);
116  else
117  buffer.reset(module->size());
118 
119 
120  // background
121 
122  simbad.generateHits(*module, time_range, buffer);
123 
124  this->push_back(JDAQSuperFrame(JDAQSuperFrameHeader(chronometer, module->getID())));
125 
126  simbad(*module, buffer, *(this->rbegin()));
127  }
128 
129  } else {
130 
131  for (JMap_t::iterator i = data.begin(); i != data.end(); ++i) {
132 
133  const JModule& module = simbad.getModule(JModuleAddress(i->first));
134 
135  this->push_back(JDAQSuperFrame(JDAQSuperFrameHeader(chronometer, module.getID())));
136 
137  simbad(module, i->second, *(this->rbegin()));
138  }
139  }
140  }
141  }
static const JRange< double, std::less< double > > DEFAULT_RANGE
Default range.
Definition: JRange.hh:506
void reset(size_t size)
Reset buffers.
Data structure for PMT analogue signal.
Data structure for a composite optical module.
Definition: JModule.hh:47
Data structure for PMT data corresponding to a detector module.
int second
index of PMT in module data structure.
Definition: JPMTAddress.hh:100
JDAQTimeslice()
Default constructor.
Structure to store the ToT mean and standard deviation of the hits produced by a nanobeacon in a sour...
bool hasPMT(const JObjectID &id) const
Has PMT.
Definition: JPMTRouter.hh:114
double getTime(const Hit &hit)
Get true time of hit.
JRange< double > JTimeRange
Type definition for time range.
const JModule & getModule(const JModuleAddress &address) const
Get module.
Definition: JPMTRouter.hh:150
JTimeRange getTimeRange(const Evt &event)
Get time range (i.e.
virtual void generateHits(const JModule &module, const JTimeRange &period, JModuleData &output) const
Generate hits.
int first
index of module in detector data structure
bool hasPMTSimulator() const
Check availability of PMT simulator.
int getID() const
Get identifier.
Definition: JObjectID.hh:54
void setDAQChronometer(const JDAQChronometer &chronometer)
Set DAQ chronometer.
Address of module in detector data structure.
bool hasK40Simulator() const
Check availability of K40 simulator.
bool hasCLBSimulator() const
Check availability of CLB simulator.
Address of PMT in detector data structure.
Definition: JPMTAddress.hh:32
Auxiliary class to convert DAQ/trigger hit time to/from Monte Carlo hit time.
double getNPE(const Hit &hit)
Get true charge of hit.
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition: JDAQ.hh:26
const JPMTAddress & getAddress(const JObjectID &id) const
Get address of PMT.
Definition: JPMTRouter.hh:78

Member Function Documentation

void KM3NETDAQ::JDAQTimeslice::clear ( )
inlineinherited

Clear data.

Definition at line 98 of file JDAQTimeslice.hh.

99  {
100  for (iterator i = this->begin(); i != this->end(); ++i) {
101  i->clear();
102  }
103 
105  }
JDAQTimeslice& KM3NETDAQ::JDAQTimeslice::add ( const JDAQTimeslice timeslice)
inlineinherited

Add another timeslice.

Parameters
timeslicetimeslice
Returns
this timeslice

Definition at line 134 of file JDAQTimeslice.hh.

135  {
136  using namespace std;
137 
139 
140  for (const_iterator i = this->begin(); i != this->end(); ++i) {
141  buffer[i->getModuleIdentifier()] = distance(static_cast<const JDAQTimeslice&>(*this).begin(),i);
142  }
143 
144  for (JDAQTimeslice::const_iterator i = timeslice.begin(); i != timeslice.end(); ++i) {
145 
146  map<JDAQModuleIdentifier, int>::const_iterator p = buffer.find(i->getModuleIdentifier());
147 
148  if (p != buffer.end()) {
149 
150  JDAQSuperFrame& frame = this->at(p->second);
151 
152  frame.add(*i);
153 
154  sort(frame.begin(), frame.end());
155 
156  } else {
157 
158  this->push_back(*i);
159  }
160  }
161 
162  return *this;
163  }
std::ostream& KM3NETDAQ::JDAQTimeslice::print ( std::ostream &  out,
const bool  lpr = false 
) const
inlineinherited

Print DAQ Timeslice.

Parameters
outoutput stream
lprlong print
Returns
output stream

Definition at line 209 of file JDAQTimeslice.hh.

210  {
211  using namespace std;
212 
213  out << this->ClassName() << endl;
214  out << dynamic_cast<const JDAQPreamble&> (*this) << endl;
215  out << dynamic_cast<const JDAQChronometer&>(*this) << endl;
216 
217  for (JDAQTimeslice::const_iterator frame = this->begin(); frame != this->end(); ++frame) {
218 
219  out << ' ' << setw(5) << frame->getModuleID();
220  out << ' ' << setw(6) << frame->getLength();
221  out << ' ' << setw(6) << frame->getDataType();
222  out << ' ' << setw(6) << frame->getTimesliceStart();
223  out << ' ' << setw(8) << hex << frame->getStatus() << dec;
224  out << '|' << setw(8) << hex << frame->getFIFOStatus() << dec;
225  out << ' ' << setw(6) << frame->size();
226 
227  if (!lpr) {
228 
229  if (!frame->empty()) {
230 
231  out << ' ' << setw(10) << frame-> begin()->getT();
232  out << " ... ";
233  out << ' ' << setw(10) << frame->rbegin()->getT();
234  }
235 
236  out << endl;
237 
238  } else {
239 
240  out << endl;
241 
242  int n = 1;
243 
244  for (JDAQFrame::const_iterator hit = frame->begin(); hit != frame->end(); ++hit, ++n) {
245  out << setw(2) << (int) hit->getPMT() << ' ' << setw(8) << hit->getT() << (n%10 == 0 ? '\n' : ' ');
246  }
247 
248  out << endl;
249  }
250  }
251 
252  return out;
253  }
const JDAQHit * const_iterator
Definition: JDAQFrame.hh:74
virtual int KM3NETDAQ::JDAQTimeslice::getSize ( ) const
inlinevirtualinherited

Get size of object.

Returns
number of bytes

Reimplemented from KM3NETDAQ::JDAQPreamble.

Definition at line 260 of file JDAQTimeslice.hh.

261  {
262  int len = 0;
263 
264  len += JDAQPreamble ::sizeOf();
266 
267  len += sizeof(int);
268 
269  for (const_iterator frame = begin(); frame != end(); ++frame) {
270  len += frame->getSize();
271  }
272 
273  return len;
274  }
static int sizeOf()
Get size of object.
Definition: JDAQHeader.hh:103
static int sizeOf()
Get size of object.
KM3NETDAQ::JDAQTimeslice::ClassDef ( JDAQTimeslice  ,
 
)
inherited
KM3NETDAQ::JDAQPreamble::ClassDef ( JDAQPreamble  ,
 
)
inherited
KM3NETDAQ::JDAQTimesliceHeader::ClassDef ( JDAQTimesliceHeader  ,
 
)
inherited
KM3NETDAQ::JDAQHeader::ClassDef ( JDAQHeader  ,
 
)
inherited
KM3NETDAQ::JDAQChronometer::ClassDef ( JDAQChronometer  ,
 
)
inherited
static int KM3NETDAQ::JDAQPreamble::sizeOf ( )
inlinestaticinherited

Get size of object.

Returns
number of bytes

Definition at line 110 of file JDAQPreamble.hh.

111  {
113  }
static int sizeOf()
Get size of object.
static int KM3NETDAQ::JDAQHeader::sizeOf ( )
inlinestaticinherited

Get size of object.

Returns
number of bytes

Definition at line 103 of file JDAQHeader.hh.

104  {
105  return (JDAQChronometer::sizeOf());
106  }
static int sizeOf()
Get size of object.
void KM3NETDAQ::JDAQPreamble::setLength ( )
inlineprotectedinherited

Set length.

Definition at line 134 of file JDAQPreamble.hh.

135  {
136  length = getSize();
137  }
virtual int getSize() const
Get size of object.
int KM3NETDAQ::JDAQAbstractPreamble::getLength ( ) const
inlineinherited

Get length.

Returns
number of bytes

Definition at line 48 of file JDAQAbstractPreamble.hh.

49  {
50  return length;
51  }
int KM3NETDAQ::JDAQAbstractPreamble::getDataType ( ) const
inlineinherited

Get data type.

Returns
data type

Definition at line 59 of file JDAQAbstractPreamble.hh.

60  {
61  return type;
62  }
KM3NETDAQ::JDAQAbstractPreamble::ClassDefNV ( JDAQAbstractPreamble  ,
 
)
inherited
const JDAQTimesliceHeader& KM3NETDAQ::JDAQTimesliceHeader::getDAQTimesliceHeader ( ) const
inlineinherited

Get DAQ time slice header.

Returns
DAQ time slice header

Definition at line 49 of file JDAQTimesliceHeader.hh.

50  {
51  return static_cast<const JDAQTimesliceHeader&>(*this);
52  }
JDAQTimesliceHeader()
Default constructor.
const JDAQHeader& KM3NETDAQ::JDAQHeader::getDAQHeader ( ) const
inlineinherited

Get DAQ header.

Returns
DAQ header

Definition at line 51 of file JDAQHeader.hh.

52  {
53  return static_cast<const JDAQHeader&>(*this);
54  }
JDAQHeader()
Default constructor.
Definition: JDAQHeader.hh:31
void KM3NETDAQ::JDAQHeader::setDAQHeader ( const JDAQHeader header)
inlineinherited

Set DAQ header.

Parameters
headerDAQ header

Definition at line 62 of file JDAQHeader.hh.

63  {
64  static_cast<JDAQHeader&>(*this) = header;
65  }
JDAQHeader()
Default constructor.
Definition: JDAQHeader.hh:31
const JDAQChronometer& KM3NETDAQ::JDAQChronometer::getDAQChronometer ( ) const
inlineinherited

Get DAQ chronometer.

Returns
DAQ chronometer

Definition at line 88 of file JDAQChronometer.hh.

89  {
90  return static_cast<const JDAQChronometer&>(*this);
91  }
JDAQChronometer()
Default constructor.
void KM3NETDAQ::JDAQChronometer::setDAQChronometer ( const JDAQChronometer chronometer)
inlineinherited

Set DAQ chronometer.

Parameters
chronometerDAQ chronometer

Definition at line 99 of file JDAQChronometer.hh.

100  {
101  static_cast<JDAQChronometer&>(*this) = chronometer;
102  }
JDAQChronometer()
Default constructor.
int KM3NETDAQ::JDAQChronometer::getDetectorID ( ) const
inlineinherited

Get detector identifier.

Returns
detector identifier

Definition at line 110 of file JDAQChronometer.hh.

111  {
112  return detector_id;
113  }
int KM3NETDAQ::JDAQChronometer::getRunNumber ( ) const
inlineinherited

Get run number.

Returns
run number

Definition at line 121 of file JDAQChronometer.hh.

122  {
123  return run;
124  }
int KM3NETDAQ::JDAQChronometer::getFrameIndex ( ) const
inlineinherited

Get frame index.

Returns
frame index

Definition at line 132 of file JDAQChronometer.hh.

133  {
134  return frame_index;
135  }
const JDAQUTCExtended& KM3NETDAQ::JDAQChronometer::getTimesliceStart ( ) const
inlineinherited

Get start of timeslice.

Returns
timeslice start

Definition at line 144 of file JDAQChronometer.hh.

145  {
146  return timeslice_start;
147  }
void KM3NETDAQ::JDAQChronometer::setRunNumber ( const int  run)
inlineinherited

Set run number.

Parameters
runrun number

Definition at line 155 of file JDAQChronometer.hh.

156  {
157  this->run = run;
158  }
void KM3NETDAQ::JDAQChronometer::setFrameIndex ( const int  frame_index)
inlineinherited

Set frame index.

Parameters
frame_indexframe index

Definition at line 166 of file JDAQChronometer.hh.

167  {
168  this->frame_index = frame_index;
169  }
void KM3NETDAQ::JDAQChronometer::setTimesliceStart ( const JDAQUTCExtended timeslice_start)
inlineinherited

Set timeslice start time.

Parameters
timeslice_starttimeslice start time

Definition at line 177 of file JDAQChronometer.hh.

178  {
180  }

Member Data Documentation

int KM3NETDAQ::JDAQAbstractPreamble::length
protectedinherited

Definition at line 79 of file JDAQAbstractPreamble.hh.

int KM3NETDAQ::JDAQAbstractPreamble::type
protectedinherited

Definition at line 80 of file JDAQAbstractPreamble.hh.

int KM3NETDAQ::JDAQChronometer::detector_id
protectedinherited

Definition at line 237 of file JDAQChronometer.hh.

int KM3NETDAQ::JDAQChronometer::run
protectedinherited

Definition at line 238 of file JDAQChronometer.hh.

int KM3NETDAQ::JDAQChronometer::frame_index
protectedinherited

Definition at line 239 of file JDAQChronometer.hh.

JDAQUTCExtended KM3NETDAQ::JDAQChronometer::timeslice_start
protectedinherited

Definition at line 240 of file JDAQChronometer.hh.


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