Jpp 19.3.0-rc.2
the software that should make you happy
Loading...
Searching...
No Matches
JModK40.cc
Go to the documentation of this file.
1
2#include <string>
3#include <iostream>
4#include <iomanip>
5#include <vector>
6
7#include "TROOT.h"
8#include "TFile.h"
9#include "TRandom3.h"
10
11#include "km3net-dataformat/offline/Head.hh"
12#include "km3net-dataformat/offline/MultiHead.hh"
13#include "km3net-dataformat/offline/Evt.hh"
14
15#include "JAAnet/JAAnetToolkit.hh"
16#include "JSirene/JSirene.hh"
17
18#include "JDetector/JDetector.hh"
19#include "JDetector/JDetectorToolkit.hh"
20#include "JDetector/JK40DefaultSimulator.hh"
21#include "JDetector/JTimeRange.hh"
22
23#include "JSupport/JMultipleFileScanner.hh"
24#include "JSupport/JFileRecorder.hh"
25#include "JSupport/JMonteCarloFileSupportkit.hh"
26#include "JSupport/JSupport.hh"
27#include "JSupport/JMeta.hh"
28
29#include "JROOT/JRandom.hh"
30
31#include "Jeep/JParser.hh"
32#include "Jeep/JMessage.hh"
33
34
35/**
36 * \file
37 *
38 * Auxiliary program to generate noise in Monte Carlo event.\n
39 * This program will be deprecated (JTriggerEfficiency.cc should be used instead).
40 * \author mdejong
41 */
42int main(int argc, char **argv)
43{
44 using namespace std;
45 using namespace JPP;
46
47 JMultipleFileScanner<Evt> inputFile;
48 JFileRecorder<JTYPELIST<JAAnetTypes_t, JMeta, JRootTypes_t>::typelist> outputFile;
49 JLimit_t& numberOfEvents = inputFile.getLimit();
50 string detectorFile;
51 JK40Rates rates_Hz;
52 JTimeRange period;
53 JRandom seed;
54 int debug;
55
56 try {
57
58 JParser<> zap("Auxiliary program to generate noise in Monte Carlo event.");
59
60 zap['f'] = make_field(inputFile) = JPARSER::initialised();
61 zap['o'] = make_field(outputFile) = "modk40.root";
62 zap['n'] = make_field(numberOfEvents) = JLimit::max();
63 zap['a'] = make_field(detectorFile);
64 zap['T'] = make_field(period) = JTimeRange::DEFAULT_RANGE();
65 zap['B'] = make_field(rates_Hz) = JPARSER::initialised();
66 zap['S'] = make_field(seed) = 0;
67 zap['d'] = make_field(debug) = 1;
68
69 zap(argc, argv);
70 }
71 catch(const exception &error) {
72 FATAL(error.what() << endl);
73 }
74
75
76 seed.set(gRandom);
77
78 const int NPE = 1;
79
80
81 JDetector detector;
82
83 try {
84 load(detectorFile, detector);
85 }
86 catch(const JException& error) {
87 FATAL(error);
88 }
89
90
91 if (period == JTimeRange::DEFAULT_RANGE()) {
92
93 const double Tmax = getMaximalTime(detector);
94
95 period = JTimeRange(-Tmax, +Tmax);
96 }
97
98
99 const JK40DefaultSimulator modk40(rates_Hz);
100
101 outputFile.open();
102
103 if (!outputFile.is_open()) {
104 FATAL("Error opening file " << outputFile << endl);
105 }
106
107 outputFile.put(JMeta(argc, argv));
108 outputFile.put(*gRandom);
109
110 if (!inputFile.empty()) {
111
112 JMultipleFileScanner<Head> io(inputFile);
113
114 io >> outputFile;
115
116 while (inputFile.hasNext()) {
117
118 STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
119
120 Evt* event = inputFile.next();
121
122 JTimeRange time_range(JTimeRange::DEFAULT_RANGE());
123
124 vector<Hit>::iterator __end = event->mc_hits.end();
125
126 for (vector<Hit>::iterator i = event->mc_hits.begin(); i != __end; ) {
127
128 if (is_noise(*i)) {
129
130 iter_swap(i, --__end);
131
132 } else {
133
134 time_range.include(getTime(*i));
135
136 ++i;
137 }
138 }
139
140 event->mc_hits.erase(__end, event->mc_hits.end());
141
142 if (time_range.is_valid())
143 time_range.add(period);
144 else
145 time_range = period;
146
147
148 JModuleData buffer;
149
150 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
151
152 buffer.reset(module->size());
153
154 modk40.generateHits(*module, time_range, buffer);
155
156 for (unsigned int pmt = 0; pmt != buffer.size(); ++pmt) {
157
158 const JModuleData::value_type& frame = buffer[pmt];
159
160 for (JModuleData::value_type::const_iterator hit = frame.begin(); hit != frame.end(); ++hit) {
161
162 event->mc_hits.push_back(JHit_t(event->mc_hits.size() + 1,
163 module->getPMT(pmt).getID(),
164 HIT_TYPE_NOISE,
165 0,
166 hit->t_ns,
167 NPE));
168 }
169 }
170 }
171
172
173 outputFile.put(*event);
174 }
175 } else {
176
177 for (counter_type counter = 0; counter != numberOfEvents; ++counter) {
178
179 STATUS("event: " << setw(10) << counter << '\r'); DEBUG(endl);
180
181 Evt event;
182
183 JModuleData buffer;
184
185 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
186
187 buffer.reset(module->size());
188
189 modk40.generateHits(*module, period, buffer);
190
191 for (unsigned int pmt = 0; pmt != buffer.size(); ++pmt) {
192
193 const JModuleData::value_type& frame = buffer[pmt];
194
195 for (JModuleData::value_type::const_iterator hit = frame.begin(); hit != frame.end(); ++hit) {
196
197 event.mc_hits.push_back(JHit_t(event.mc_hits.size() + 1,
198 module->getPMT(pmt).getID(),
199 HIT_TYPE_NOISE,
200 0,
201 hit->t_ns,
202 NPE));
203 }
204 }
205 }
206
207
208 outputFile.put(event);
209 }
210 }
211 STATUS(endl);
212
213 outputFile.put(*gRandom);
214 outputFile.close();
215}
JAAnetToolkit.hh
Definition of hit and track types and auxiliary methods for handling Monte Carlo data.
outputFile
string outputFile
Definition JDAQTimesliceSelector.cc:37
JDetector.hh
Data structure for detector geometry and calibration.
JFileRecorder.hh
Recording of objects on file according a format that follows from the file name extension.
JMessage.hh
General purpose messaging.
DEBUG
#define DEBUG(A)
Message macros.
Definition JMessage.hh:62
STATUS
#define STATUS(A)
Definition JMessage.hh:63
FATAL
#define FATAL(A)
Definition JMessage.hh:67
debug
int debug
debug level
Definition JSirene.cc:72
JMeta.hh
ROOT I/O of application specific meta data.
main
int main(int argc, char **argv)
Definition JModK40.cc:42
JMultipleFileScanner.hh
Scanning of objects from multiple files according a format that follows from the extension of each fi...
JParser.hh
Utility class to parse command line options.
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
JSupport.hh
ROOT TTree parameter settings of various packages.
JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96
JDETECTOR::JK40DefaultSimulatorInterface::generateHits
virtual void generateHits(const JModule &module, const JTimeRange &period, JModuleData &output) const
Generate hits.
Definition JK40DefaultSimulatorInterface.hh:134
JDETECTOR::JK40DefaultSimulator
Default implementation of the simulation of K40 background.
Definition JK40DefaultSimulator.hh:33
JDETECTOR::JModuleData
Data structure for PMT data corresponding to a detector module.
Definition JPMTSimulator.hh:216
JDETECTOR::JModuleData::reset
void reset(size_t size)
Reset buffers.
Definition JPMTSimulator.hh:231
JLANG::JException
General exception.
Definition JException.hh:24
JPARSER::JParser
Utility class to parse command line options.
Definition JParser.hh:1698
JSUPPORT::JFileRecorder
Object writing to file.
Definition JFileRecorder.hh:44
JSUPPORT::JMultipleFileScanner
General purpose class for object reading from a list of file names.
Definition JMultipleFileScanner.hh:306
JSUPPORT::JMultipleFileScanner::hasNext
virtual bool hasNext() override
Check availability of next element.
Definition JMultipleFileScanner.hh:409
JSUPPORT::JMultipleFileScanner::getCounter
counter_type getCounter() const
Get counter.
Definition JMultipleFileScanner.hh:382
JSUPPORT::JMultipleFileScanner::next
virtual const pointer_type & next() override
Get next element.
Definition JMultipleFileScanner.hh:457
JTOOLS::JRange::is_valid
bool is_valid() const
Check validity of range.
Definition JRange.hh:311
JTOOLS::JRange< double >::DEFAULT_RANGE
static JRange< double, std::less< double > > DEFAULT_RANGE()
Definition JRange.hh:555
JTOOLS::JRange::include
range_type & include(argument_type x)
Include given value to range.
Definition JRange.hh:397
JTOOLS::JRange::add
range_type & add(argument_type x)
Add offset.
Definition JRange.hh:446
JAANET::HIT_TYPE_NOISE
@ HIT_TYPE_NOISE
Random noise.
Definition JAAnetToolkit.hh:85
JAANET::is_noise
bool is_noise(const Hit &hit)
Verify hit origin.
Definition JAAnetToolkit.hh:120
JACOUSTICS::JTimeRange
JTOOLS::JRange< double > JTimeRange
Type definition for time range (unit [s]).
Definition JAcoustics/JTimeRange.hh:19
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition JDetectorToolkit.hh:508
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:43
JROOT::counter_type
Long64_t counter_type
Type definition for counter.
Definition JROOT/JCounter.hh:24
JSYSTEM::getTime
const char * getTime()
Get current local time conform ISO-8601 standard.
Definition JDateAndTime.hh:55
KM3NETDAQ::getMaximalTime
double getMaximalTime(const double R_Hz)
Get maximal time for given rate.
Definition JHighRateVeto.hh:35
Evt
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Definition Evt.hh:21
Evt::mc_hits
std::vector< Hit > mc_hits
MC: list of MC truth hits.
Definition Evt.hh:48
JAANET::detector
Detector file.
Definition JHead.hh:227
JMATH::JRandom
Template definition of random value generator.
Definition JMath/JRandom.hh:26
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
JPHYSICS::JK40Rates
Auxiliary class for K40 rates.
Definition JK40Rates.hh:41
JSIRENE::JHit_t
Auxiliary class to set-up Hit.
Definition JSirene.hh:58
JSUPPORT::JLimit
Auxiliary class for defining the range of iterations of objects.
Definition JLimit.hh:45
JSUPPORT::JLimit::max
static counter_type max()
Get maximum counter value.
Definition JLimit.hh:128
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition JMeta.hh:72
Generated by doxygen 1.12.0