FortranCallback.h

Go to the documentation of this file.

#ifndef gen_FortranCallback_h
#define gen_FortranCallback_h
 
#include "SimDataFormats/GeneratorProducts/interface/LHECommonBlocks.h"
 
#include "GeneratorInterface/LHEInterface/interface/LHERunInfo.h"
#include "GeneratorInterface/LHEInterface/interface/LHEEvent.h"
 
namespace HepMC {
  class GenEvent;
}
 
namespace gen {
 
  class FortranCallback {
  public:
    static FortranCallback* getInstance();
 
    //void setLHEEvent(lhef::LHEEvent* lhev) { fPartonLevel = lhev; }
    void setLHERunInfo(lhef::LHERunInfo* lheri) { fRunBlock = lheri; }
    void setLHEEvent(lhef::LHEEvent* lhee) { fEventBlock = lhee; }
 
    void resetIterationsPerEvent() { fIterationsPerEvent = 0; }
 
    void fillHeader();
    void fillEvent();
 
    int getIterationsPerEvent() const { return fIterationsPerEvent; }
 
  private:
    // ctor
 
    FortranCallback();
 
    // data member(s)
 
    lhef::LHERunInfo* fRunBlock;
    lhef::LHEEvent* fEventBlock;
    int fIterationsPerEvent;
 
    static FortranCallback* fInstance;
  };
 
  // --** Implementation **---
 
  inline FortranCallback::FortranCallback()
      //   : fPartonLevel(0)
      : fRunBlock(nullptr), fEventBlock(nullptr), fIterationsPerEvent(0) {}
 
  inline void FortranCallback::fillHeader() {
    if (fRunBlock == nullptr)
      return;
 
    //const lhef::HEPRUP* heprup = &(fRunBlock->heprup());
    const lhef::HEPRUP* heprup = fRunBlock->getHEPRUP();
 
    lhef::CommonBlocks::fillHEPRUP(heprup);
 
    return;
  }
 
  inline void FortranCallback::fillEvent() {
    //if ( fPartonLevel == 0 ) return;
    //const lhef::HEPEUP* hepeup = fPartonLevel->getHEPEUP();
 
    if (fEventBlock == nullptr)
      return;
 
    const lhef::HEPEUP* hepeup = fEventBlock->getHEPEUP();
 
    if (fIterationsPerEvent++) {
      hepeup_.nup = 0;
      return;
    }
 
    lhef::CommonBlocks::fillHEPEUP(hepeup);
 
    return;
  }
 
}  // namespace gen
 
#endif