JSummary.cc File Reference

Auxiliary program to blend and write summary data. More...

#include <string>
#include <iostream>
#include <iomanip>
#include "JDAQ/JDAQSummarysliceIO.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JTrigger/JTriggerParameters.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JSupportToolkit.hh"
#include "JSupport/JMeta.hh"
#include "JSummaryslice/JSummaryslice.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

Auxiliary program to blend and write summary data.

This application can be used to convert measured rates of any detector to a consistent set of rates for a simulated complete detector.
The rates as well as the high-rate veto and FIFO (almost) full will then be included in the run-by-run simulations using JTriggerEfficiency.cc.

The option -f <inputFile> corresponds to a list of input files containing summary data (commonly referred to as "measured rates").
These data can originate from any detector.

The option -a <detectorFile> corresponds to the detector to be used in simulations.
This could be a complete detector.

The option -@ <parameters> corresponds to the trigger parameters.
This is needed here because in the run-by-run simulations, the trigger parameters are taken from the input file.

The application JSummaryWriter.cc can be used to write summary data according user defined input rates.

Author

mdejong

Definition in file JSummary.cc.

Function Documentation

main()

int main	(	int	argc,
		char **	argv )

Definition at line 46 of file JSummary.cc.

{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;
 
  JMultipleFileScanner<JDAQSummaryslice> inputFile;
  JFileRecorder       <JTYPELIST<JDAQSummaryslice, JTriggerParameters, JMeta>::typelist> outputFile;
  JDAQUTCTimeRange   UTC;
  string             detectorFile;
  JTriggerParameters parameters;
  int                run_number;
  int                debug;
 
  try { 
 
    JParser<> zap("Auxiliary program to blend and write summary data.");
    
    zap['f'] = make_field(inputFile);
    zap['U'] = make_field(UTC);
    zap['o'] = make_field(outputFile);
    zap['a'] = make_field(detectorFile);
    zap['@'] = make_field(parameters)          = JPARSER::initialised();
    zap['R'] = make_field(run_number)          = 1;
    zap['d'] = make_field(debug)               = 0;
    
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  JDetector detector;
 
  try {
    load(detectorFile, detector);
  }
  catch(const JException& error) {
    FATAL(error);
  }
 
 
  // check availability of sufficient summary data
 
  unsigned int number_of_modules = 0;
 
  while (inputFile.hasNext() && number_of_modules < detector.size()) {
    number_of_modules += inputFile.next()->size();
  }
  
  if (number_of_modules < detector.size()) {
    FATAL("Input summary data insufficient: " << number_of_modules << " < " << detector.size() << endl);
  }
 
  inputFile.rewind();
 
 
  outputFile.open();
 
  outputFile.put(JMeta(argc, argv));  
  outputFile.put(parameters);
 
  int frame_index = 1;
 
  for (JDAQUTCExtended utc = UTC.getLowerLimit(); utc < UTC.getUpperLimit(); ) {
 
    NOTICE("event: " << setw(10) << frame_index << '\r'); DEBUG(endl);
 
    JSummaryslice summary(JDAQChronometer(detector.getID(),
                                          run_number, 
                                          frame_index, 
                                          utc),
                          inputFile, 
                          detector);
    
    if (summary.size() == detector.size()) {
 
      outputFile.put(summary);
 
      utc.addTimeNanoSecond(getFrameTime());
 
      frame_index += 1;
 
    } else {
 
      inputFile.rewind();
    }
  }
  NOTICE(endl);
 
  outputFile.close();
}