makeFluxHistogramsL.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include "TROOT.h"
#include "eventRates.h"
 
int main(int argc, char* argv[]) {

  std::string user = std::string("beam"); 
  std::string version = std::string("v3r0p1");
  std::string macro = std::string("Nominal"); 
  std::string current = std::string("200kA"); 
  std::string location = std::string("LBNEFD"); 
  std::string physics_list = std::string("QGSP_BERT"); 
  int n_files = 250; 
  int start_index = 1; 
  double pot_per_file = 100000.;
  bool on_grid = false;
  bool dk2nu=false;
  
   char* myFile, myPath, myOutPath;

   for (int i = 1; i < argc; i++) { 
      std::string tmpStr = std::string(argv[i]);
         if (tmpStr.find("--user=") == 0) {
                        //01234567
             user = tmpStr.substr(7);
         } else if (tmpStr.find("--version=") == 0) {
                               //012345678910123456789
             version = tmpStr.substr(10);
         } else if (tmpStr.find("--macro=") == 0) {
                               //012345678910123456789
             macro = tmpStr.substr(8);
         } else if (tmpStr.find("--current=") == 0) {
                               //012345678910123456789
             current = tmpStr.substr(10);
         } else if (tmpStr.find("--location=") == 0) {
                               //012345678910123456789
             location = tmpStr.substr(11);
             std::cout<<"location is "<<location<<std::endl;
         } else if (tmpStr.find("--physics_list=") == 0) {
                               //012345678910123456789
             physics_list = tmpStr.substr(15);
             std::cout<<"Physics list is "<<physics_list<<std::endl;
         } else if (tmpStr.find("--n_files=") == 0) {
                               //012345678910123456789
             n_files = atoi((tmpStr.substr(10)).c_str());
         } else if (tmpStr.find("--start_index=") == 0) {
                               //012345678910123456789
             start_index = atoi((tmpStr.substr(14)).c_str());
         } else if (tmpStr.find("--pot_per_file=") == 0) {
                               //012345678910123456789
             pot_per_file = atof((tmpStr.substr(15)).c_str());
             std::cout<<"pot_per_file "<<pot_per_file<<std::endl;
         } else if (tmpStr.find("--on_grid=") == 0) {
                               //012345678910123456789
             std::string log_on_grid = tmpStr.substr(10);
             if ((log_on_grid == std::string("false")) || (log_on_grid == std::string("False"))) on_grid=false;
             if ((log_on_grid == std::string("true")) || (log_on_grid == std::string("True"))) on_grid=true;
         } else if (tmpStr.find("--dk2nu=") == 0) {
                               //012345678910123456789
             std::string log_dk2nu = tmpStr.substr(8);
             if ((log_dk2nu == std::string("false")) || (log_dk2nu == std::string("False"))) dk2nu=false;
             if ((log_dk2nu == std::string("true")) || (log_dk2nu == std::string("True"))) dk2nu=true;
         } else {
             std::cout << " Unrecognized argument .\n" << tmpStr << " Fatal, exit here! \n " << std::endl;
             exit(2); 
        }
   }
   std::string G4LBNEWORKDIR = getenv("G4LBNEWORKDIR");
   std::string DK2NU = getenv("DK2NU");
   
  eventRates t(user, version, macro, current, location, physics_list, n_files, start_index, pot_per_file, on_grid, dk2nu);
                        
  std::cout << " Event Rate class loaded.. " << std::endl;
  std::cout << " Arguments are " << std::endl;
  std::cout << " user = " << user << std::endl; 
  std::cout << "  version = " << version << std::endl;
  std::cout << "  macro = " << macro << std::endl;
  std::cout << "  current = " << current << std::endl;
  std::cout << "  location = " << location << std::endl;
  std::cout << "  physics_list = " << physics_list << std::endl;
  std::cout << "  n_files " << n_files << std::endl; 
  std::cout << "  start_index " <<  start_index << std::endl; 
  std::cout << "  pot_per_file " <<  pot_per_file << std::endl;
  if (on_grid) std::cout << " We are on the grid .. " << std::endl; 
  else  std::cout << " We are NOT on the grid .. " << std::endl;
  if (dk2nu) std::cout << " We are using Dk2nu .. " << std::endl; 
  else  std::cout << " We are using the old NuMI style Ntuple .. " << std::endl;
 
  t.Loop();  

}