62 int main(
int argc,
char **argv)
66 using namespace KM3NETDAQ;
76 JCalibration_t calibrationFile;
85 JParser<> zap(
"Auxiliary program to convert fit results to Evt format.\
86 \nThe option -L corresponds to the name of a shared library \
87 \nand function so to rearrange the order of fit results.");
91 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
102 catch(
const exception& error) {
103 FATAL(error.what() << endl);
108 typedef JParallelFileScanner_t::multi_pointer_type multi_pointer_type;
113 if (detectorFile !=
"") {
124 unique_ptr<JDynamics> dynamics;
130 dynamics->load(calibrationFile);
132 catch(
const exception& error) {
133 if (!calibrationFile.empty()) {
149 }
catch(
const exception& error) {}
156 copy(buffer, header);
169 STATUS(
"Processing: " << *i << endl);
171 JParallelFileScanner_t
in(*i);
174 in.setLimit(inputFile.getLimit());
184 center = get<Vec>(head);
187 }
catch(
const exception& error) {}
190 while (in.hasNext()) {
192 STATUS(
"event: " << setw(10) << in.getCounter() <<
'\r');
DEBUG(endl);
194 multi_pointer_type ps = in.next();
200 dynamics->update(*tev);
203 JFIT::JEvt::iterator __end = evt->end();
205 if (numberOfFits > 0) {
206 advance(__end = evt->begin(), min(numberOfFits, evt->size()));
210 partial_sort(evt->begin(), __end, evt->end(),
qualitySorter);
215 if (mc.getEntries() != 0) {
235 if (!pmt_router->empty()) {
239 if (pmt_router.hasPMT(i->pmt_id)) {
241 const JPMTAddress address = pmt_router.getAddress(i->pmt_id);
243 const JPMT& pmt = pmt_router.getPMT(address);
245 i->dom_id =
id.
getID();
246 i->channel_id =
id.getPMTAddress();
252 FATAL(
"Missing PMT" << i->pmt_id << endl);
257 if (!mod_router->empty()) {
261 if (mod_router.hasModule(i->dom_id)) {
268 const JPMT& pmt = mod_router.getPMT(
id);
274 i->pmt_id = pmt.
getID();
279 i->a = cpu.
getNPE(i->tot);
283 FATAL(
"Missing module " << i->dom_id << endl);
288 for (JFIT::JEvt::const_iterator fit = evt->begin(); fit != __end; ++fit) {
292 trk.
id = out.
trks.size() + 1;
293 trk.
pos =
Vec(fit->getX(), fit->getY(), fit->getZ());
294 trk.
dir =
Vec(fit->getDX(), fit->getDY(), fit->getDZ());
297 trk.
lik = fit->getQ();
300 for (JHistory::const_iterator i = fit->getHistory().begin(); i != fit->getHistory().end(); ++i) {
304 for (
int i = 0; i != fit->getN(); ++i) {
305 trk.
fitinf.push_back(fit->getW(i));
310 out.
trks.push_back(trk);
313 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.
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
static void setSlewing(const bool slewing)
Set slewing option.
double t
track time [ns] (when the particle is at pos )
static JDetectorMechanics getMechanics
Function object to get string mechanics.
const JCalibration & getCalibration() const
Get calibration.
ROOT TTree parameter settings.
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.
General purpose sorter of fit results.
*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
General purpose class for parallel reading of objects from a single file or multiple files...
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.
double E
Energy [GeV] (either MC truth or reconstructed)
Data structure for detector geometry and calibration.
bool qualitySorter(const JFit &first, const JFit &second)
Comparison of fit results.
ROOT TTree parameter settings.
Head getHeader(const JMultipleFileScanner_t &file_list)
Get Monte Carlo header.
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
Parallel scanning of objects from a single file or multiple files according a format that follows fro...
Data structure for PMT geometry, calibration and status.
JPosition3D getPosition(const Vec &pos)
Get position.
std::vector< double > fitinf
place to store additional fit info, for jgandalf, see JFitParameters.hh
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.
counter_type advance(counter_type &counter, const counter_type value, const counter_type limit=std::numeric_limits< counter_type >::max())
Advance counter.
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.
Vec pos
postion of the track at time t [m]
Address of PMT in detector data structure.
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Data structure for set of track fit results.
std::vector< Trk > trks
list of reconstructed tracks (can be several because of prefits,showers, etc).
double getLivetime(const std::string &file_name)
Get data taking live time.
std::vector< int > rec_stages
list of identifyers of succesfull fitting stages resulting in this track
General purpose class for object reading from a list of file names.
std::vector< double > error_matrix
(NxN) error covariance matrix for fit parameters (stored as linear vector)
double lik
likelihood or lambda value (for aafit, lambda)
Utility class to parse command line options.
std::vector< Hit > mc_hits
MC: list of MC truth hits.
static const int JPP_RECONSTRUCTION_TYPE
KM3NeT Data Definitions v2.1.0-7-g97845ea https://git.km3net.de/common/km3net-dataformat.
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
const JLimit & getLimit() const
Get limit.
Compass event data types.
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
The Trk class represents a Monte Carlo (MC) particle as well as a reconstructed track/shower.
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.
then usage $script< detector file >< tripodfile >< stage > input file nInput files correspond to the output of JAcousticsEventBuilder[.sh] nFirst stage eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY eval JPrintDetector a $DETECTOR O CAN source JAcoustics sh $DETECTOR_ID typeset A CONFIGURATION for key in Tmax_s
int rec_type
identifyer for the overall fitting algorithm/chain/strategy