HerwigHepMCFragmentationModule.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <framework/core/Module.h>
#include <mdst/dataobjects/MCParticle.h>
#include <mdst/dataobjects/MCParticleGraph.h>
#include <generators/evtgen/EvtGenInterface.h>
 
#include <string>
#include <map>
#include <memory>
 
namespace HepMC { class GenEvent; }
 
namespace Belle2 {
 
  class HerwigFragHelper;

  class HerwigHepMCFragmentationModule : public Module {
 
  public:
    HerwigHepMCFragmentationModule();
    virtual ~HerwigHepMCFragmentationModule();

    virtual void initialize() override;
    virtual void terminate() override;
    virtual void event() override;
 
  private:
    // Steerable module parameters
    std::string m_workDir{};          
    std::string m_decFile{};          
    std::string m_userDecFile{};      
    bool m_useEvtGen{true};           
    bool m_coherentMixing{true};      
    bool m_keepTempFiles{false};      
 
    // Non-steerable internal state: read-only after initialize()
    std::string m_hepmcPath{};        

    std::map<int, HepMC::GenEvent*> m_hepmcEvents;

    std::map<int, std::string> m_sidecarMap;
 
    std::unique_ptr<HerwigFragHelper> m_helper;    
    MCParticleGraph m_mcParticleGraph;              
    EvtGenInterface m_evtGenInterface;              
 
    int m_eventCounter{0};            
    int m_successCounter{0};          

    bool loadManifest();

    bool loadHepMCEvents();

    bool loadKKMCSidecar(int globalEvtNum, std::vector<int>& quarkIndices);

    bool applyEvtGenDecays();
  };

} // namespace Belle2