Jpp master_rocky-44-g75b7c4f75
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 "Jeep/JParser.hh"
30#include "Jeep/JMessage.hh"
31
32
33/**
34 * \file
35 *
36 * Auxiliary program to generate noise in Monte Carlo event.\n
37 * This program will be deprecated (JTriggerEfficiency.cc should be used instead).
38 * \author mdejong
39 */
40int main(int argc, char **argv)
41{
42 using namespace std;
43 using namespace JPP;
44
45 JMultipleFileScanner<Evt> inputFile;
46 JFileRecorder<JTYPELIST<JAAnetTypes_t, JMeta, JRootTypes_t>::typelist> outputFile;
47 JLimit_t& numberOfEvents = inputFile.getLimit();
48 string detectorFile;
49 JK40Rates rates_Hz;
50 JTimeRange period;
51 UInt_t seed;
52 int debug;
53
54 try {
55
56 JParser<> zap("Auxiliary program to generate noise in Monte Carlo event.");
57
58 zap['f'] = make_field(inputFile) = JPARSER::initialised();
59 zap['o'] = make_field(outputFile) = "modk40.root";
60 zap['n'] = make_field(numberOfEvents) = JLimit::max();
61 zap['a'] = make_field(detectorFile);
62 zap['T'] = make_field(period) = JTimeRange::DEFAULT_RANGE();
63 zap['B'] = make_field(rates_Hz) = JPARSER::initialised();
64 zap['S'] = make_field(seed) = 0;
65 zap['d'] = make_field(debug) = 1;
66
67 zap(argc, argv);
68 }
69 catch(const exception &error) {
70 FATAL(error.what() << endl);
71 }
72
73
74 gRandom->SetSeed(seed);
75
76 const int NPE = 1;
77
78
79 JDetector detector;
80
81 try {
82 load(detectorFile, detector);
83 }
84 catch(const JException& error) {
85 FATAL(error);
86 }
87
88
89 if (period == JTimeRange::DEFAULT_RANGE()) {
90
91 const double Tmax = getMaximalTime(detector);
92
93 period = JTimeRange(-Tmax, +Tmax);
94 }
95
96
97 const JK40DefaultSimulator modk40(rates_Hz);
98
99 outputFile.open();
100
101 if (!outputFile.is_open()) {
102 FATAL("Error opening file " << outputFile << endl);
103 }
104
105 outputFile.put(JMeta(argc, argv));
106 outputFile.put(*gRandom);
107
108 if (!inputFile.empty()) {
109
110 JMultipleFileScanner<Head> io(inputFile);
111
112 io >> outputFile;
113
114 while (inputFile.hasNext()) {
115
116 STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
117
118 Evt* event = inputFile.next();
119
120 JTimeRange time_range(JTimeRange::DEFAULT_RANGE());
121
122 vector<Hit>::iterator __end = event->mc_hits.end();
123
124 for (vector<Hit>::iterator i = event->mc_hits.begin(); i != __end; ) {
125
126 if (is_noise(*i)) {
127
128 iter_swap(i, --__end);
129
130 } else {
131
132 time_range.include(getTime(*i));
133
134 ++i;
135 }
136 }
137
138 event->mc_hits.erase(__end, event->mc_hits.end());
139
140 if (time_range.is_valid())
141 time_range.add(period);
142 else
143 time_range = period;
144
145
146 JModuleData buffer;
147
148 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
149
150 buffer.reset(module->size());
151
152 modk40.generateHits(*module, time_range, buffer);
153
154 for (unsigned int pmt = 0; pmt != buffer.size(); ++pmt) {
155
156 const JModuleData::value_type& frame = buffer[pmt];
157
158 for (JModuleData::value_type::const_iterator hit = frame.begin(); hit != frame.end(); ++hit) {
159
160 event->mc_hits.push_back(JHit_t(event->mc_hits.size() + 1,
161 module->getPMT(pmt).getID(),
162 HIT_TYPE_NOISE,
163 0,
164 hit->t_ns,
165 NPE));
166 }
167 }
168 }
169
170
171 outputFile.put(*event);
172 }
173 } else {
174
175 for (counter_type counter = 0; counter != numberOfEvents; ++counter) {
176
177 STATUS("event: " << setw(10) << counter << '\r'); DEBUG(endl);
178
179 Evt event;
180
181 JModuleData buffer;
182
183 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
184
185 buffer.reset(module->size());
186
187 modk40.generateHits(*module, period, buffer);
188
189 for (unsigned int pmt = 0; pmt != buffer.size(); ++pmt) {
190
191 const JModuleData::value_type& frame = buffer[pmt];
192
193 for (JModuleData::value_type::const_iterator hit = frame.begin(); hit != frame.end(); ++hit) {
194
195 event.mc_hits.push_back(JHit_t(event.mc_hits.size() + 1,
196 module->getPMT(pmt).getID(),
197 HIT_TYPE_NOISE,
198 0,
199 hit->t_ns,
200 NPE));
201 }
202 }
203 }
204
205
206 outputFile.put(event);
207 }
208 }
209 STATUS(endl);
210
211 outputFile.put(*gRandom);
212 outputFile.close();
213}
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:69
JMeta.hh
ROOT I/O of application specific meta data.
main
int main(int argc, char **argv)
Definition JModK40.cc:40
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:505
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
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