12 #include "TProfile2D.h"
35 int main(
int argc,
char **argv)
39 using namespace KM3NETDAQ;
54 JParser<> zap(
"Auxiliary program to determine time slewing from K40 data.");
59 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
69 catch(
const exception &error) {
70 FATAL(error.what() << endl);
74 using namespace KM3NETDAQ;
76 if (!multiplicity.is_valid()) {
FATAL(
"Invalid multiplicity " << multiplicity << endl); }
77 if ( multiplicity.getLowerLimit() < 2) {
FATAL(
"Invalid multiplicity " << multiplicity << endl); }
78 if (!totVeto_ns .
is_valid()) {
FATAL(
"Invalid time over threshold " << totVeto_ns << endl); }
92 FATAL(
"Empty detector." << endl);
100 TH1D h0(
"h0", NULL, 41, -TMax_ns, +TMax_ns);
101 TH1D h1(
"h1", NULL, 50, -0.5, 49.5);
102 TProfile2D hx(
"hx", NULL, 50, -0.5, 49.5, 50, -0.5, 49.5, -2*TMax_ns, +2*TMax_ns);
103 TProfile h2(
"h2", NULL, 100, 0.5, 100.5, -2*TMax_ns, +2*TMax_ns);
111 while (inputFile.hasNext()) {
113 STATUS(
"event: " << setw(10) << inputFile.getCounter() <<
'\r');
DEBUG(endl);
117 for (JDAQTimeslice::const_iterator super_frame = timeslice->begin(); super_frame != timeslice->end(); ++super_frame) {
119 if (router.
hasModule(super_frame->getModuleID()) && !super_frame->empty()) {
125 buffer.resize(super_frame->size());
134 sort(buffer.begin(), buffer.end());
140 while (++q != buffer.end() && q->getT() - p->getT() <= TMax_ns ) {}
150 module.
getPMT(__q->getPMT()));
154 h0.Fill(__q->getT() - __p->getT());
155 h1.Fill(__q->getToT());
156 hx.Fill(__p->getToT(), __q->getToT(), __q->getT() - __p->getT());
158 if (totVeto_ns(__p->getToT())) {
159 h2.Fill(__q->getToT(), __q->getT() - __p->getT());
178 TF1
f1(
"f1",
"[0]*exp([1]*sqrt(x) + [2]*x) + [3]");
180 f1.SetParameter(0, h2.GetMaximum());
181 f1.SetParameter(1, -0.01);
182 f1.SetParameter(2, -0.05);
183 f1.SetParameter(3, h2.GetMinimum());
185 h2.ProjectionX()->Fit(&f1,
"",
"same");
187 for (
int i = 0; i != f1.GetNpar(); ++i) {
188 cout <<
"\tstatic double p" << i <<
"() { return " << setw(9) << setprecision(5) << f1.GetParameter(i) <<
"; }" << endl;
Utility class to parse command line options.
int main(int argc, char *argv[])
ROOT TTree parameter settings of various packages.
Basic data structure for L0 hit.
const JModule & getModule(const JObjectID &id) const
Get module parameters.
Data structure for a composite optical module.
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
bool is_valid(const json &js)
Check validity of JSon data.
Router for direct addressing of module data in detector data structure.
Data structure for detector geometry and calibration.
const JPolynome f1(1.0, 2.0, 3.0)
Function.
Auxiliary class for defining the range of iterations of objects.
Reduced data structure for L0 hit.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
double getDot(const JFirst_t &first, const JSecond_t &second)
Get dot product of objects.
then $JPP_DIR examples JDetector JSlewing P slewing
const JPMT & getPMT(const int index) const
Get PMT.
General purpose messaging.
Scanning of objects from multiple files according a format that follows from the extension of each fi...
Direct access to module in detector data structure.
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Auxiliary class to define a range between two values.
General purpose class for object reading from a list of file names.
Utility class to parse command line options.
bool hasModule(const JObjectID &id) const
Has module.
const JLimit & getLimit() const
Get limit.
do set_variable DETECTOR_TXT $WORKDIR detector
KM3NeT DAQ constants, bit handling, etc.
#define DEBUG(A)
Message macros.