JNuisance.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>
#include <memory>
#include <map>
 
#include "TROOT.h"
#include "TFile.h"
#include "TH1D.h"
#include "TRandom3.h"
 
#include "JAstronomy/JNuisance.hh"
 
#include "JLang/JVectorize.hh"
 
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
/**
 * Toy simulation to study nuisance.
 */
int main(int argc, char**argv)
{
  using namespace std;
  using namespace JPP;
  
  size_t           numberOfEvents;
  string           outputFile;
  nuisance_type    nuisance;
  UInt_t           seed;
  int              debug;
 
  try {
    
    JParser<> zap;
 
    zap['o'] = make_field(outputFile)                                  = "nuisance.root";
    zap['n'] = make_field(numberOfEvents);
    zap['u'] = make_field(nuisance,     "nuisance: \"<type> (parameters)+\"\n\twhere <type> can be:" << get_keys(nuisance_helper));
    zap['S'] = make_field(seed)                = 0;
    zap['d'] = make_field(debug)               = 1;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  gRandom->SetSeed(seed);
 
  TFile out(outputFile.c_str(), "recreate");
 
  TH1D h0("h0", NULL,  100, -5.0, +5.0);
  
  for (size_t counter = 0; counter != numberOfEvents; ++counter) {
 
    STATUS("counter: "<< setw(8) << counter << '\r'); DEBUG(endl);
 
    h0.Fill(nuisance->get());
  }
  
  out.Write();
  out.Close();
}