LaserAlignmentProducer.cc

Go to the documentation of this file.

// system include files
#include "FWCore/Framework/interface/Event.h"

// user include files
#include "Alignment/LaserAlignmentSimulation/plugins/LaserAlignmentProducer.h"
#include "FWCore/Framework/interface/MakerMacros.h" 

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"

//
// constructors and destructor
//
LaserAlignmentProducer::LaserAlignmentProducer(const edm::ParameterSet&) :
  EDProducer(),
  theEvent(nullptr)
{
  //register your products
  produces<edm::HepMCProduct>("unsmeared");

  //now do what ever other initialization is needed
}


LaserAlignmentProducer::~LaserAlignmentProducer()
{
  // no need to cleanup theEvent since it's done in HepMCProduct
}

// ------------ method called to produce the event  ------------
void LaserAlignmentProducer::produce(edm::Event& iEvent, const edm::EventSetup&)
{
  // create the event
  theEvent = new HepMC::GenEvent();

  // create a primary vertex
  HepMC::GenVertex * theVtx = new HepMC::GenVertex(HepMC::FourVector(0.,0.,0.));

  // add a particle to the vertex; this is needed to avoid crashes in OscarProducer. Use a 
  // electron neutrino, with zero energy and mass
  HepMC::GenParticle * theParticle = new HepMC::GenParticle(HepMC::FourVector(0.,0.,0.,0.),12,1);
  
  theVtx->add_particle_out(theParticle);

  // add the vertex to the event
  theEvent->add_vertex(theVtx);

  // set the event number
  theEvent->set_event_number(iEvent.id().event());
  // set the signal process id
  theEvent->set_signal_process_id(20);

  // create an empty output collection
  auto theOutput = std::make_unique<edm::HepMCProduct>();
  theOutput->addHepMCData(theEvent);
   
  // put the output to the event
  iEvent.put(std::move(theOutput));
}

//define this as a plug-in

DEFINE_FWK_MODULE(LaserAlignmentProducer);