Fork me on GitHub

Changeset 44 in svn for trunk/Delphes.cpp


Ignore:
Timestamp:
Nov 19, 2008, 8:42:24 PM (16 years ago)
Author:
severine ovyn
Message:

add log file

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/Delphes.cpp

    r43 r44  
    8888        exit(1);
    8989  }
     90  //create output log-file name
     91  string LogName = outputfilename.erase(outputfilename.find(".root"));
     92  LogName = LogName+"_run.log";
    9093
    9194  TFile *outputFile = TFile::Open(outputfilename.c_str(), "RECREATE"); // Creates the file, but should be closed just after
    … …  
    9598  ifstream infile(inputFileList.c_str());
    9699  infile >> line; // the first line determines the type of input files
     100
     101  //read the datacard input file
     102  string DetDatacard("");
     103  if(argc==4)  DetDatacard =argv[3];
     104  RESOLution *DET = new RESOLution();
     105  DET->ReadDataCard(DetDatacard);
     106  DET->Logfile(LogName);
     107
     108  todo(LogName.c_str());
    97109 
    98110  DataConverter *converter=0;
    … …  
    100112  if(strstr(line.c_str(),".hep"))
    101113    {                           
    102       cout<<"*************************************************************************"<<endl;
    103       cout<<"************         StdHEP file format detected           **************"<<endl;
    104       cout<<"************     Starting convertion to TRoot format       **************"<<endl;
    105       cout<<"*************************************************************************"<<endl;
     114      cout<<"#**********************************************************************"<<endl;
     115      cout<<"#**********         StdHEP file format detected           *************"<<endl;
     116      cout<<"#***********     Starting convertion to TRoot format     **************"<<endl;
     117      cout<<"#**********************************************************************"<<endl;
    106118      converter = new STDHEPConverter(inputFileList,outputfilename);//case ntpl file in input list
    107119    }
    108120  else if(strstr(line.c_str(),".lhe"))
    109121    {
    110       cout<<"*************************************************************************"<<endl;
    111       cout<<"************          LHEF file format detected            **************"<<endl;
    112       cout<<"************     Starting convertion to TRoot format       **************"<<endl;
    113       cout<<"*************************************************************************"<<endl;
     122      cout<<"#**********************************************************************"<<endl;
     123      cout<<"#***********          LHEF file format detected            ************"<<endl;
     124      cout<<"#***********     Starting convertion to TRoot format       ************"<<endl;
     125      cout<<"#**********************************************************************"<<endl;
    114126      converter = new LHEFConverter(inputFileList,outputfilename);//case ntpl file in input list
    115127    }
    116128  else if(strstr(line.c_str(),".root"))
    117129    {
    118       cout<<"*************************************************************************"<<endl;
    119       cout<<"************         h2root file format detected           **************"<<endl;
    120       cout<<"************     Starting convertion to TRoot format       **************"<<endl;
    121       cout<<"*************************************************************************"<<endl;
     130      cout<<"#**********************************************************************"<<endl;
     131      cout<<"#**********         h2root file format detected           *************"<<endl;
     132      cout<<"#**********     Starting convertion to TRoot format       *************"<<endl;
     133      cout<<"#**********************************************************************"<<endl;
    122134      converter = new HEPEVTConverter(inputFileList,outputfilename);//case ntpl file in input list
    123135    }
    … …  
    157169  TRootRomanPotHits *elementRP220, *elementFP420;
    158170 
    159   //read the datacard input file
    160   string DetDatacard("");
    161   if(argc==4)  DetDatacard =argv[3];
    162   RESOLution *DET = new RESOLution();   
    163   DET->ReadDataCard(DetDatacard);
    164  
    165171  TLorentzVector genMomentum(0,0,0,0);
    166172  TLorentzVector genMomentumCalo(0,0,0,0);
    … …  
    235241    // set up a siscone jet definition
    236242#ifdef ENABLE_PLUGIN_SISCONE
    237     int npass = 0;               // do infinite number of passes
    238     double protojet_ptmin = 0.0; // use all protojets
    239     plugins = new fastjet::SISConePlugin (DET->CONERADIUS,DET->OVERLAPTHRESHOLD,npass, protojet_ptmin);
     243    plugins = new fastjet::SISConePlugin (DET->CONERADIUS,DET->OVERLAPTHRESHOLD,DET->NPASS, DET->PROTOJET_PTMIN);
    240244    jet_def = fastjet::JetDefinition(plugins);
    241245#else
    … …  
    485489     
    486490      // run the jet clustering with the above jet definition
    487       cout<<input_particles.size()<<endl;
    488491      if(input_particles.size()!=0)
    489492        {
    … …  
    492495          double ptmin = 5.0;
    493496          inclusive_jets = clust_seq.inclusive_jets(ptmin);
    494 cout<<"inclusive_jets "<<inclusive_jets.size()<<endl;     
    495497          // sort jets into increasing pt
    496498          sorted_jets = sorted_by_pt(inclusive_jets);
    497499        }
    498 cout<<"sorted_jets.size() "<<sorted_jets.size()<<endl;
    499500      for (unsigned int i = 0; i < sorted_jets.size(); i++) {
    500501        TLorentzVector JET;
Note: See TracChangeset for help on using the changeset viewer.