Jpp/public/JModK40_8cc_source.html

#include <string>

#include <iostream>

#include <iomanip>

#include <vector>


#include "TROOT.h"

#include "TFile.h"

#include "TRandom3.h"


#include "km3net-dataformat/offline/Head.hh"

#include "km3net-dataformat/offline/MultiHead.hh"

#include "km3net-dataformat/offline/Evt.hh"


#include "JAAnet/JAAnetToolkit.hh"

#include "JSirene/JSirene.hh"


#include "JDetector/JDetector.hh"

#include "JDetector/JDetectorToolkit.hh"

#include "JDetector/JK40DefaultSimulator.hh"

#include "JDetector/JTimeRange.hh"


#include "JSupport/JMultipleFileScanner.hh"

#include "JSupport/JFileRecorder.hh"

#include "JSupport/JMonteCarloFileSupportkit.hh"

#include "JSupport/JSupport.hh"

#include "JSupport/JMeta.hh"


#include "JROOT/JRandom.hh"


#include "Jeep/JParser.hh"

#include "Jeep/JMessage.hh"


/**

 * \file

 *

 * Auxiliary program to generate noise in Monte Carlo event.\n

 * This program will be deprecated (JTriggerEfficiency.cc should be used instead).

 * \author mdejong

 */


int main(int argc, char **argv)

{

  using namespace std;

  using namespace JPP;


  JMultipleFileScanner<Evt>                                              inputFile;

  JFileRecorder<JTYPELIST<JAAnetTypes_t, JMeta, JRootTypes_t>::typelist> outputFile;

  JLimit_t&      numberOfEvents = inputFile.getLimit();

  string         detectorFile;

  JK40Rates      rates_Hz;

  JTimeRange     period;

  JRandom        seed;

  int            debug;


  try {


    JParser<> zap("Auxiliary program to generate noise in Monte Carlo event.");


    zap['f'] = make_field(inputFile)       = JPARSER::initialised();

    zap['o'] = make_field(outputFile)      = "modk40.root";

    zap['n'] = make_field(numberOfEvents)  = JLimit::max();

    zap['a'] = make_field(detectorFile);

    zap['T'] = make_field(period)          = JTimeRange::DEFAULT_RANGE();

    zap['B'] = make_field(rates_Hz)        = JPARSER::initialised();

    zap['S'] = make_field(seed)            = 0;

    zap['d'] = make_field(debug)           = 1;


    zap(argc, argv);

  }

  catch(const exception &error) {

    FATAL(error.what() << endl);

  }


  seed.set(gRandom);


  const int NPE = 1;


  JDetector detector;


  try {

    load(detectorFile, detector);

  }

  catch(const JException& error) {

    FATAL(error);

  }


  if (period == JTimeRange::DEFAULT_RANGE()) {


    const double Tmax = getMaximalTime(detector);


    period = JTimeRange(-Tmax, +Tmax);

  }


  const JK40DefaultSimulator modk40(rates_Hz);


  outputFile.open();


  if (!outputFile.is_open()) {

    FATAL("Error opening file " << outputFile << endl);

  }


  outputFile.put(JMeta(argc, argv));

  outputFile.put(*gRandom);


  if (!inputFile.empty()) {


    JMultipleFileScanner<Head> io(inputFile);


    io >> outputFile;


    while (inputFile.hasNext()) {


      STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);


      Evt* event = inputFile.next();


      JTimeRange time_range(JTimeRange::DEFAULT_RANGE());


      vector<Hit>::iterator __end = event->mc_hits.end();


      for (vector<Hit>::iterator i = event->mc_hits.begin(); i != __end; ) {


        if (is_noise(*i)) {


          iter_swap(i, --__end);


        } else {


          time_range.include(getTime(*i));


          ++i;

        }

      }


      event->mc_hits.erase(__end, event->mc_hits.end());


      if (time_range.is_valid())

        time_range.add(period);

      else

        time_range = period;


      JModuleData buffer;


      for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {


        buffer.reset(module->size());


        modk40.generateHits(*module, time_range, buffer);


        for (unsigned int pmt = 0; pmt != buffer.size(); ++pmt) {


          const JModuleData::value_type& frame = buffer[pmt];


          for (JModuleData::value_type::const_iterator hit = frame.begin(); hit != frame.end(); ++hit) {


            event->mc_hits.push_back(JHit_t(event->mc_hits.size() + 1,

                                            module->getPMT(pmt).getID(),

                                            HIT_TYPE_NOISE,

                                            0,

                                            hit->t_ns,

                                            NPE));

          }

        }

      }


      outputFile.put(*event);

    }

  } else {


    for (counter_type counter = 0; counter != numberOfEvents; ++counter) {


      STATUS("event: " << setw(10) << counter << '\r'); DEBUG(endl);


      Evt event;


      JModuleData buffer;


      for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {


        buffer.reset(module->size());


        modk40.generateHits(*module, period, buffer);


        for (unsigned int pmt = 0; pmt != buffer.size(); ++pmt) {


          const JModuleData::value_type& frame = buffer[pmt];


          for (JModuleData::value_type::const_iterator hit = frame.begin(); hit != frame.end(); ++hit) {


            event.mc_hits.push_back(JHit_t(event.mc_hits.size() + 1,

                                           module->getPMT(pmt).getID(),

                                           HIT_TYPE_NOISE,

                                           0,

                                           hit->t_ns,

                                           NPE));

          }

        }

      }


      outputFile.put(event);

    }

  }

  STATUS(endl);


  outputFile.put(*gRandom);

  outputFile.close();

}


Evt.hh

Head.hh

JAAnetToolkit.hh
Definition of hit and track types and auxiliary methods for handling Monte Carlo data.

outputFile
string outputFile
Definition JDAQTimesliceSelector.cc:37

JDetectorToolkit.hh

JTimeRange.hh

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.

JK40DefaultSimulator.hh

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

JMonteCarloFileSupportkit.hh

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

JRandom.hh

JSirene.hh

JSupport.hh
ROOT TTree parameter settings of various packages.

MultiHead.hh

JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96

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:217

JDETECTOR::JModuleData::reset
void reset(size_t size)
Reset buffers.
Definition JPMTSimulator.hh:232

JLANG::JException
General exception.
Definition JException.hh:24

JOSCPROB::JOscProbInterpolator
Template definition of a multi-dimensional oscillation probability interpolation table.
Definition JOscProbInterpolator.hh:71

JTOOLS::JRange< double >

JTOOLS::JRange< double >::DEFAULT_RANGE
static JRange< double, std::less< double > > DEFAULT_RANGE()
Default range.
Definition JRange.hh:555

std::vector
Definition JSTDTypes.hh:15

JAANET::HIT_TYPE_NOISE
@ HIT_TYPE_NOISE
Random noise.
Definition JAAnetToolkit.hh:89

JAANET::is_noise
bool is_noise(const Hit &hit)
Verify hit origin.
Definition JAAnetToolkit.hh:124

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:520

JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:45

JSYSTEM::getTime
const char * getTime()
Get current local time conform ISO-8601 standard.
Definition JDateAndTime.hh:56

KM3NETDAQ::getMaximalTime
double getMaximalTime(const double R_Hz)
Get maximal time for given rate.
Definition JHighRateVeto.hh:35

std
Definition JSTDTypes.hh:14

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:60

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