BeamHaloProducer.h

Go to the documentation of this file.

#ifndef BeamHalo_Producer_h
#define BeamHalo_Producer_h

#include <map>
#include <string>

#include "HepMC/GenEvent.h"

#include "CLHEP/Random/RandFlat.h"

#include "FWCore/Framework/interface/one/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

namespace CLHEP {
  class HepRandomEngine;
}

namespace edm
{
  class BeamHaloProducer : public one::EDProducer<EndRunProducer, one::WatchLuminosityBlocks, one::SharedResources> {
  public:

    BeamHaloProducer(const ParameterSet &);
    ~BeamHaloProducer() override;

    void setRandomEngine(CLHEP::HepRandomEngine* v);

  private:
    bool call_ki_bhg_init(long& seed);
        bool call_bh_set_parameters(int* ival, float* fval,const std::string cval_string);
    bool call_ki_bhg_fill(int& iret, float& weight);
    bool call_ki_bhg_stat(int& iret);

  private:

    void produce(Event & e, const EventSetup & es) override;
    void endRunProduce(Run & r, const EventSetup & es) override;
    void beginLuminosityBlock(LuminosityBlock const&, EventSetup const&) override;
    void endLuminosityBlock(LuminosityBlock const&, EventSetup const&) override { }

    void clear();

    HepMC::GenEvent  *evt;

        int GENMOD_;
    int LHC_B1_;
    int LHC_B2_;
    int IW_MUO_;
    int IW_HAD_;
    float EG_MIN_;
    float EG_MAX_;
    std::string G3FNAME_;

    bool isInitialized_;
  };

}

#endif