Configuration.cc

Go to the documentation of this file.

#include <iostream>
#include <sstream>
#include <stdio.h>
#include "HepMC/GenEvent.h"
#include "HepMC/IO_GenEvent.h"
#include "HepMC/GenParticle.h"
#include "HepMC/GenVertex.h"
#include "HepMC/IO_AsciiParticles.h"
#include "HepMC/SimpleVector.h"
#include "CLHEP/Vector/LorentzVector.h"

using namespace std;

// ROOT headers
#include "TH1.h"
#include <TH2.h>
#include "TFile.h"
#include "TMath.h"
#include "TLorentzVector.h"

//tau analysis header
#include "../include/Configuration.h"

//**********************

//starting functions
Configuration::Configuration()
{
  // the default randomseed will be generated from the time
  srand(time(0)%900000000);
  rseed=rand();
  
  // default: 1000 events
  nevents=1000;
  
  // default Outputfilename
  OutputFileName="output.root";
  
  // default: no Algorithm specified
  dijet = false;
  z = false;
  wplusjet = false;
  tau = false;
  top = false;
  ue = false;
  
  // default Jet parameter
  coneRadius = 1.0;
  overlapThreshold = 0.1;
  jet_ptmin = 5.0;
  max_eta = 100;
  min_pt = 0.1;
  
  FSR=true;
  ISR=true;
  MI=true;
  
  IsOK=true;
}

Configuration::Configuration( std::string configfile )
{
  Configuration();
  IsOK=read(configfile);
}

Configuration::~Configuration()
{
  while(!ConfigFileNames.empty())
    ConfigFileNames.pop_back();
}


int  Configuration::read(std::string configfile){
  string line;
  ifstream myfile (configfile.c_str());
  std::cout << "readConfigFile: Read in Configuration from file " << configfile << std::endl;
  if (myfile.is_open())
    {
      // read line by line of the file
      while (! myfile.eof() )
        {
          std::string parameter;
          std::string str1;
          int mylong=-9999;
          double myfloat=-9999.9;
          bool mybool=false;
          
          // write read lines into string line
          getline (myfile,line);
          
          // now split the line into the subwords
          std::istringstream instring(line);
          std::istringstream instring2(line);
          std::istringstream instring3(line);
          
          // and fill string or float
          instring >> parameter >> str1 ;
          instring2 >> parameter >> myfloat ;
          instring3 >> parameter >> mybool ;
          mylong=long(myfloat);
          
          // check for Comments. If the line starts with # than skip the line
          std::string::size_type pos = parameter.find("#");
          if(pos != std::string::npos) continue;
          // check for Comments. If the line starts with * than skip the line
          std::string::size_type pos2 = parameter.find("*");
          if(pos2 != std::string::npos) continue;
          // check for Comments. If the line starts with '/' than skip the line
          std::string::size_type pos3 = parameter.find("/");
          if(pos3 != std::string::npos) continue;
          
          // skip if the string is empty
          if(parameter.empty()) continue;
          
          
          // read pysub entries
          if ((parameter == "rseed") && (mylong!=-9999)) {
            rseed=mylong;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << rseed << std::endl;
          }
          else if ((parameter == "nevents") && (mylong!=-9999)){
            nevents=(unsigned int)mylong;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << nevents << std::endl;
          }
          else if ((parameter == "OutputFileName") && (!str1.empty())){
            OutputFileName=str1;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << OutputFileName << std::endl;
          }
          else if ((parameter == "ConfigFileNames") && (!str1.empty())){
            ConfigFileNames.push_back(str1);
            std::cout << "readConfigFile: Parameter " << parameter << "=" << ConfigFileNames.back() << std::endl;
          }
          else if ((parameter == "dijet_analysis") && (!str1.empty())){
            if( (str1 == "on") ||  (str1 == "true") ||  (str1 == "True")) dijet=true;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << dijet << std::endl;
          }
          else if ((parameter == "z_analysis")  && (!str1.empty())){
            if( (str1 == "on") ||  (str1 == "true") ||  (str1 == "True")) z=true;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << z << std::endl;
          }
          else if ((parameter == "wplusjet_analysis")  && (!str1.empty())){
            if( (str1 == "on")  ||  (str1 == "true") ||  (str1 == "True") ) wplusjet=true;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << wplusjet << std::endl;
          }
          else if ((parameter == "tau_analysis")  && (!str1.empty())){
            if( (str1 == "on")  ||  (str1 == "true") ||  (str1 == "True") ) tau=true;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << tau << std::endl;
          }
          else if ((parameter == "top_analysis")  && (!str1.empty())){
            if( (str1 == "on")  ||  (str1 == "true") ||  (str1 == "True") ) top=true;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << top << std::endl;
          }
          else if ((parameter == "ue_analysis")  && (!str1.empty())){
            if( (str1 == "on")  ||  (str1 == "true") ||  (str1 == "True") ) ue=true;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << ue << std::endl;
          }
          else if ((parameter == "jet_coneRadius") && (myfloat!=-9999.9)){
            coneRadius=myfloat;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << coneRadius << std::endl;
          }
          else if ((parameter == "jet_overlapThreshold") && (myfloat!=-9999.9)){
            overlapThreshold=myfloat;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << overlapThreshold << std::endl;
          }
          else if ((parameter == "jet_ptmin") && (myfloat!=-9999.9)){
            jet_ptmin=myfloat;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << jet_ptmin << std::endl;
          }
          else if ((parameter == "max_eta") && (myfloat!=-9999.9)){
            max_eta=myfloat;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << max_eta << std::endl;
          }
          else if ((parameter == "min_pt") && (myfloat!=-9999.9)){
            min_pt=myfloat;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << min_pt << std::endl;
          }
          else if ((parameter == "FSR")  && (!str1.empty())){
            if( (str1 == "on")  ||  (str1 == "true") ||  (str1 == "True") ) FSR=true;
            else FSR=false;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << FSR << std::endl;
          }
          else if ((parameter == "ISR")  && (!str1.empty())){
            if( (str1 == "on")  ||  (str1 == "true") ||  (str1 == "True") ) ISR=true;
            else ISR=false;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << ISR << std::endl;
          }
          else if ((parameter == "MI")  && (!str1.empty())){
            if( (str1 == "on")  ||  (str1 == "true") ||  (str1 == "True") ) MI=true;
            else MI=false;
            std::cout << "readConfigFile: Parameter " << parameter << "=" << MI << std::endl;
          }
          else{
            std::cout << "readConfigFile: Parameter " << parameter << " not recognized!" << std::endl;
            return 1;
          }
        }
      myfile.close();
    }
  else 
    {
      cout << "readConfigFile: Unable to open file"; 
      return 1;
    }
  return 0;
}

---