210 JParser<> zap(
"Auxiliary program to edit PMT parameters map.");
212 zap[
'a'] =
make_field(detectorFile,
"detector file.") =
"";
213 zap[
'D'] =
make_field(detectorID,
"detector identifier (in absence of detector file).") = 0;
217 zap[
'E'] =
make_field(mu,
"expectation value for npe given two-fold coincidence") = 0.0;
218 zap[
'T'] =
make_field(T_ns,
"time-over-threshold range.") = JRange_t();
224 catch(
const exception &error) {
225 FATAL(error.what() << endl);
229 if (detectorFile !=
"") {
236 load(detectorFile, detector);
242 if (detectorID == 0) {
244 detectorID = detector.
getID();
246 }
else if (detectorID != detector.
getID()) {
248 FATAL(
"Inconsistent detector identifier " << detectorID <<
" != " << detector.
getID() << endl);
252 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
254 for (
unsigned int pmt = 0; pmt != module->size(); ++pmt) {
260 DEBUG(
"Setting default parameters for PMT " <<
id << endl);
271 FATAL(
"Invalid detector identifier " << detectorID << endl);
280 for (
vector< JModifier<JPMTPhysicalAddress> >::const_iterator i = mod.begin(); i != mod.end(); ++i) {
282 if (compare(i->address, address)) {
284 DEBUG(
"Modifying parameters for PMT " << ps->first <<
' ' << i->action <<
' ' << i->key <<
' ' << i->value << endl);
286 if (!i->apply(ps->second)) {
287 ERROR(
"No valid action: " << *i << endl);
298 for (
vector< JModifier<JPMTIdentifier> >::const_iterator i = daq.begin(); i != daq.end(); ++i) {
300 if (compare(ps->first, i->address)) {
302 DEBUG(
"Modifying parameters for PMT " << ps->first <<
' ' << i->action <<
' ' << i->key <<
' ' << i->value << endl);
304 if (!i->apply(ps->second)) {
305 ERROR(
"No valid action: " << *i << endl);
315 DEBUG(
"Correct measured QE for two-hit probability " << mu << endl);
324 }
else if (mu < 0.0) {
326 FATAL(
"Invalid expection value for two-hit probability " << mu << endl);
330 if (T_ns != JRange_t()) {
332 DEBUG(
"Correct measured QE for time-over-threshold range " << T_ns << endl);
340 i->second.QE *= (cpu.getIntegralOfChargeProbability(i->second.threshold,
341 cpu.getNPE(T_ns.getUpperLimit(), NPE),
344 cpu.getIntegralOfChargeProbability(cpu.getNPE(T_ns.getLowerLimit(), NPE),
345 cpu.getNPE(T_ns.getUpperLimit(), NPE),
Utility class to parse command line options.
bool hasDetectorAddressMap(const int id)
Check if detector address map is available.
*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
Lookup table for PMT addresses in detector.
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
const JModuleAddressMap & get(const int id) const
Get module address map.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
int getID() const
Get identifier.
Auxiliary class for map of PMT parameters.
JDetectorAddressMap & getDetectorAddressMap()
Get detector address map.
PMT analogue signal processor.
virtual const char * what() const override
Get error message.
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Data structure for PMT physical address.
do set_variable DETECTOR_TXT $WORKDIR detector
#define DEBUG(A)
Message macros.