MCFileSource.cc

Go to the documentation of this file.

#include <iostream>
#include <string>


#include "FWCore/Framework/interface/Event.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "IOMC/Input/interface/HepMCFileReader.h" 
#include "IOMC/Input/interface/MCFileSource.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/GenEventInfoProduct.h"

namespace edm {

//-------------------------------------------------------------------------
MCFileSource::MCFileSource(const ParameterSet & pset, InputSourceDescription const& desc) :
  ProducerSourceFromFiles(pset, desc, false),
  reader_(HepMCFileReader::instance()), evt_(0)
{
  LogInfo("MCFileSource") << "Reading HepMC file:" << fileNames()[0];
  std::string fileName = fileNames()[0];
  // strip the file: 
  if (fileName.find("file:") == 0){
    fileName.erase(0,5);
  }  
  
  reader_->initialize(fileName);  
  produces<HepMCProduct>("generator");
  produces<GenEventInfoProduct>("generator");
}


//-------------------------------------------------------------------------
MCFileSource::~MCFileSource(){
}

//-------------------------------------------------------------------------
bool MCFileSource::setRunAndEventInfo(EventID&, TimeValue_t&) {
  // Read one HepMC event
  LogInfo("MCFileSource") << "Start Reading";
  evt_ = reader_->fillCurrentEventData(); 
  return(evt_ != nullptr);
}

//-------------------------------------------------------------------------
void MCFileSource::produce(Event &e) {
  // Store one HepMC event in the Event.

  std::auto_ptr<HepMCProduct> bare_product(new HepMCProduct());  
  bare_product->addHepMCData(evt_);
  e.put(bare_product,"generator");
  std::auto_ptr<GenEventInfoProduct> info ( new GenEventInfoProduct( evt_ ) );
  e.put(std::move(info),"generator" );
}

}