57 int main(
int argc,
char **argv)
61 using namespace KM3NETDAQ;
71 JCalibration_t calibrationFile;
78 JParser<> zap(
"Auxiliary program to convert fit results to Evt format.\
79 \nThe option -L corresponds to the name of a shared library \
80 \nand function so to rearrange the order of fit results.");
84 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
93 catch(
const exception& error) {
94 FATAL(error.what() << endl);
99 typedef JSingleFileScanner_t::pointer_type pointer_type;
104 if (detectorFile !=
"") {
115 unique_ptr<JDynamics> dynamics;
121 dynamics->load(calibrationFile);
123 catch(
const exception& error) {
124 if (!calibrationFile.empty()) {
140 }
catch(
const exception& error) {}
147 copy(buffer, header);
160 STATUS(
"Processing: " << *i << endl);
165 in.setLimit(inputFile.getLimit());
175 center = get<Vec>(head);
178 }
catch(
const exception& error) {}
181 while (in.hasNext()) {
183 STATUS(
"event: " << setw(10) << in.getCounter() <<
'\r');
DEBUG(endl);
185 pointer_type ps = in.next();
190 dynamics->update(*tev);
195 if (mc.getEntries() != 0) {
215 if (!pmt_router->empty()) {
219 if (pmt_router.hasPMT(i->pmt_id)) {
221 const JPMTAddress address = pmt_router.getAddress(i->pmt_id);
223 const JPMT& pmt = pmt_router.getPMT(address);
225 i->dom_id =
id.
getID();
226 i->channel_id =
id.getPMTAddress();
232 FATAL(
"Missing PMT" << i->pmt_id << endl);
237 if (!mod_router->empty()) {
241 if (mod_router.hasModule(i->dom_id)) {
248 const JPMT& pmt = mod_router.getPMT(
id);
254 i->pmt_id = pmt.
getID();
259 i->a = cpu.
getNPE(i->tot);
263 FATAL(
"Missing module " << i->dom_id << endl);
268 out.
id = ++number_of_events;
Router for direct addressing of PMT data in detector data structure.
double getT() const
Get calibrated time of hit.
Utility class to parse command line options.
int main(int argc, char *argv[])
ROOT TTree parameter settings of various packages.
static void setSlewing(const bool slewing)
Set slewing option.
static JDetectorMechanics getMechanics
Function object to get string mechanics.
const JCalibration & getCalibration() const
Get calibration.
ROOT TTree parameter settings.
bool read(Vec &v, std::istream &is)
Read a Vec(tor) from a stream.
Recording of objects on file according a format that follows from the file name extension.
Router for direct addressing of module data in detector data structure.
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Long64_t counter_type
Type definition for counter.
double livetime_s
Live time [s].
double getTime(const Hit &hit)
Get true time of hit.
Basic data structure for time and time over threshold information of hit.
Data structure for detector geometry and calibration.
ROOT TTree parameter settings.
Auxiliary base class for file name.
Head getHeader(const JMultipleFileScanner_t &file_list)
Get Monte Carlo header.
Scanning of objects from a single file according a format that follows from the extension of each fil...
double getToT() const
Get calibrated time over threshold of hit.
Auxiliary class for defining the range of iterations of objects.
bool slewing
time slewing of analogue signal
The Vec class is a straightforward 3-d vector, which also works in pyroot.
void push(T JHead::*pd)
Push given data member to Head.
JDirection3D getDirection(const Vec &dir)
Get direction.
int mc_run_id
MC run identifier.
Auxiliary class for recursive type list generation.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
int getID() const
Get identifier.
JAANET::start_run start_run
Data structure for PMT geometry, calibration and status.
JPosition3D getPosition(const Vec &pos)
Get position.
Auxiliary class for map of PMT parameters.
Dynamic detector calibration.
void load(const std::string &file_name)
Load mechanical model parameters from file.
Direct access to PMT in detector data structure.
General purpose messaging.
The Head class reflects the header of Monte-Carlo event files, which consists of keys (also referred ...
Direct access to module in detector data structure.
PMT analogue signal processor.
Dynamic detector calibration.
virtual double getNPE(const double tot_ns, const double eps=1.0e-3) const
Get number of photo-electrons.
Address of PMT in detector data structure.
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
double getLivetime(const std::string &file_name)
Get data taking live time.
General purpose class for object reading from a list of file names.
Utility class to parse command line options.
std::vector< Hit > mc_hits
MC: list of MC truth hits.
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
void copy(const Head &from, JHead &to)
Copy header from from to to.
PMT analogue signal processor.
int id
offline event identifier
Object reading from a list of files.
const JLimit & getLimit() const
Get limit.
Data structure for PMT parameters.
std::vector< Hit > hits
list of hits
do set_variable DETECTOR_TXT $WORKDIR detector
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
JTriggerCounter_t getCounter() const
Get trigger counter.
std::vector< Trk > mc_trks
MC: list of MC truth tracks.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.