mergeProphecy4f.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
#include "Prophecy4fMerger.h"
#include <iostream>
#include <AsgMessaging/MessageCheck.h>
 
using namespace asg::msgUserCode;
 
int main(int argc, char*argv[]){
 
    if(argc<5){
        ANA_MSG_INFO(" Prophecy4fMerger: merge Powheg and Prophecy LHE files! ");
        ANA_MSG_INFO(" Usage:  ./Prophecy4fMerger [Options]");
        ANA_MSG_INFO("         Nedeed options:");
        ANA_MSG_INFO("               --inPowheg [Powheg lhe file]");
        ANA_MSG_INFO("               --inProphecy4e    [Prophecy4f lhe file with 4e decays]");
        ANA_MSG_INFO("               --inProphecy4mu   [Prophecy4f lhe file with 4mu decays]");
        ANA_MSG_INFO("               --inProphecy2e2mu [Prophecy4f lhe file with 2e2mu decays]");
        ANA_MSG_INFO("               --outLHE [Outupt lhe file]");
        ANA_MSG_INFO("               --randomSeed [the seed]");
        ANA_MSG_INFO("         Additional options:");
        ANA_MSG_INFO("               --debug [Enable debug mode]");
        return 0;
    }
    std::string powheg, prophecy4e, prophecy4mu, prophecy2e2mu, out;
    unsigned long long seed = 0;
    bool debug=false;
    for(int a=1; a<argc; a++){
        if(!strcmp(argv[a],"--inPowheg")){
            powheg = argv[a+1];
        }
        else if(!strcmp(argv[a],"--inProphecy4e")){
            prophecy4e = argv[a+1];
        }
        else if(!strcmp(argv[a],"--inProphecy4mu")){
            prophecy4mu = argv[a+1];
        }
        else if(!strcmp(argv[a],"--inProphecy2e2mu")){
            prophecy2e2mu = argv[a+1];
        }
        else if(!strcmp(argv[a],"--outLHE")){
            out = argv[a+1];
        }
        else if(!strcmp(argv[a],"--randomSeed")){
            seed = std::stoull(argv[a+1]);
        }
        else if(!strcmp(argv[a],"--debug")){
            debug = true;
        }
    }
 
    Prophecy4fMerger merger;
    merger.setIO(powheg, prophecy4e, prophecy4mu, prophecy2e2mu, out, debug);
    merger.setRandomSeed(seed);
  try{
    merger.merge(); 
  } catch (const std::runtime_error & e){
    ANA_MSG_ERROR("Prophecy4fMerger: Exception thrown in merger.merge() step: " << e.what());
    return 1;
  }
 
    return 0;
  
}