MCParticleReplacer Class Reference

#include <MCParticleReplacer.h>

Inheritance diagram for MCParticleReplacer:

Public Member Functions
virtual void	beginJob ()
virtual void	beginRun (edm::Run &, const edm::EventSetup &)
template<typename T >
void	call_produces (const std::string &instanceName)
template<typename T >
void	call_put (T &product, const std::string &instanceName)
virtual void	endJob ()
virtual void	endRun ()
edm::StreamID	getStreamID () const
	MCParticleReplacer (const edm::ParameterSet &)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~MCParticleReplacer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Types
enum	HepMcMode { kInvalid = 0, kNew, kReplace }

Static Private Member Functions
static HepMcMode	stringToHepMcMode (const std::string &name)

Private Attributes
edm::Event *	evt_
HepMcMode	hepMcMode_
ParticleReplacerBase *	replacer_
edm::InputTag	src_
edm::InputTag	srcHepMC_
int	verbosity_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Replace muons reconstructed in selected Z –> mu+ mu- events by generator level particles, which will be passed to detector simulation & reconstruction modules to create "hybrid" events ("embedded" leptons from Monte Carlo simulation, rest of the event taken from data)

Per default, the reconstructed muons are replaced by generator level tau leptons, which are passed to TAUOLA in order to produce generator level tau decay products.

For systematic/background studies, it is possible also to:

• replace generator level muons
• "embed" electrons or muons

Author

Manuel Zeise

Version

Revision:

1.3

Id:

MCParticleReplacer.h,v 1.3 2013/01/31 09:07:18 veelken Exp

Definition at line 35 of file MCParticleReplacer.h.

Member Enumeration Documentation

enum MCParticleReplacer::HepMcMode

private

Enumerator
kInvalid
kNew
kReplace

Definition at line 62 of file MCParticleReplacer.h.

{ kInvalid = 0, kNew, kReplace };

Constructor & Destructor Documentation

MCParticleReplacer::MCParticleReplacer ( const edm::ParameterSet &  cfg )

explicit

Definition at line 19 of file MCParticleReplacer.cc.

References HLT_25ns14e33_v1_cff::algorithm, edm::ParameterSet::exists(), reco::get(), edm::ParameterSet::getParameter(), replacer_, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity_.

  : src_(cfg.getParameter<edm::InputTag>("src")),
    srcHepMC_(cfg.getParameter<edm::InputTag>("hepMcSrc")),
    hepMcMode_(stringToHepMcMode(cfg.getParameter<std::string>("hepMcMode"))),
    replacer_(0),
    evt_(0)
{
  std::string algorithm = cfg.getParameter<std::string>("algorithm");
  edm::ParameterSet cfgAlgorithm = cfg.getParameter<edm::ParameterSet>(algorithm);
  std::string pluginType = cfg.getParameter<std::string>("pluginType");
  replacer_ = ParticleReplacerPluginFactory::get()->create(pluginType, cfgAlgorithm);

  verbosity_ = ( cfg.exists("verbosity") ) ?
    cfg.getParameter<int>("verbosity") : 0;

  produces<edm::HepMCProduct>();
  produces<GenFilterInfo>("minVisPtFilter");
  replacer_->declareExtraProducts(this);
}

MCParticleReplacer::~MCParticleReplacer

(

)

Definition at line 39 of file MCParticleReplacer.cc.

References replacer_.

{
  delete replacer_;
}

Member Function Documentation

void MCParticleReplacer::beginJob ( void  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 104 of file MCParticleReplacer.cc.

References ParticleReplacerBase::beginJob(), and replacer_.

{
  
  replacer_->beginJob();
}

void MCParticleReplacer::beginRun ( edm::Run &  run, const edm::EventSetup &  es  )

virtual

Definition at line 93 of file MCParticleReplacer.cc.

References ParticleReplacerBase::beginRun(), and replacer_.

{
  replacer_->beginRun(run, es);
}

template<typename T >

void MCParticleReplacer::call_produces ( const std::string &  instanceName )

inline

Definition at line 48 of file MCParticleReplacer.h.

Referenced by ParticleReplacerZtautau::declareExtraProducts().

  {
    produces<T>(instanceName);
  }

template<typename T >

void MCParticleReplacer::call_put ( T &  product, const std::string &  instanceName  )

inline

Definition at line 54 of file MCParticleReplacer.h.

References evt_, and edm::Event::put().

Referenced by ParticleReplacerZtautau::produce().

  {
    evt_->put(product, instanceName);
  }

void MCParticleReplacer::endJob ( void  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 111 of file MCParticleReplacer.cc.

References ParticleReplacerBase::endJob(), and replacer_.

{
  replacer_->endJob();
}

void MCParticleReplacer::endRun ( void  )

virtual

Definition at line 98 of file MCParticleReplacer.cc.

References ParticleReplacerBase::endRun(), and replacer_.

{
  replacer_->endRun();
}

edm::StreamID MCParticleReplacer::getStreamID ( ) const

inline

Definition at line 59 of file MCParticleReplacer.h.

References assert(), evt_, and edm::Event::streamID().

Referenced by ParticleReplacerZtautau::produce().

{ assert(evt_ != nullptr); return evt_->streamID(); }

void MCParticleReplacer::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

virtual

Implements edm::EDProducer.

Definition at line 53 of file MCParticleReplacer.cc.

References reco::Candidate::charge(), evt_, edm::hlt::Exception, edm::Event::getByLabel(), getSelMuons(), getTheMuMinus(), getTheMuPlus(), hepMcMode_, info(), edm::RefToBase< T >::isNonnull(), kNew, kReplace, patZpeak::muons, reco::Candidate::p4(), ParticleReplacerBase::passed_, ParticleReplacerBase::produce(), edm::Event::put(), replacer_, src_, srcHepMC_, AlCaHLTBitMon_QueryRunRegistry::string, ParticleReplacerBase::tried_, and reco::Candidate::vertex().

{
  std::vector<reco::CandidateBaseRef> selMuons = getSelMuons(evt, src_);
  const reco::CandidateBaseRef muPlus  = getTheMuPlus(selMuons);
  const reco::CandidateBaseRef muMinus = getTheMuMinus(selMuons);

  std::vector<reco::Particle> muons;
  if ( muPlus.isNonnull()  ) muons.push_back(reco::Particle(muPlus->charge(), muPlus->p4(), muPlus->vertex(), -15, 0, true));
  if ( muMinus.isNonnull() ) muons.push_back(reco::Particle(muMinus->charge(), muMinus->p4(), muMinus->vertex(), +15, 0, true));

  if ( muons.size() == 0 ) {
    edm::LogError("MCParticleReplacer") 
      << "No Z->mumu candidates or muons found !!" << std::endl;
    return;
  }
    
  evt_ = &evt;

  std::auto_ptr<HepMC::GenEvent> hepMC;
  if ( hepMcMode_ == kReplace ) {
    edm::Handle<edm::HepMCProduct> HepMCHandle;  
    evt.getByLabel(srcHepMC_, HepMCHandle);
    hepMC = replacer_->produce(muons, 0, HepMCHandle->GetEvent(), this);
  } else if( hepMcMode_ == kNew ) {
    hepMC = replacer_->produce(muons, 0, 0, this);
  } else
    throw cms::Exception("LogicError") 
      << "Invalid hepMcMode " << hepMcMode_ << " !!" << std::endl;

  if ( hepMC.get() != 0 ) {
    std::auto_ptr<edm::HepMCProduct> bare_product(new edm::HepMCProduct());  
    bare_product->addHepMCData(hepMC.release()); // transfer ownership of the HepMC:GenEvent to bare_product
    
    evt.put(bare_product);
    
    std::auto_ptr<GenFilterInfo> info(new GenFilterInfo(replacer_->tried_, replacer_->passed_));
    evt.put(info, std::string("minVisPtFilter"));
  }
}

MCParticleReplacer::HepMcMode MCParticleReplacer::stringToHepMcMode ( const std::string &  name )

staticprivate

Definition at line 44 of file MCParticleReplacer.cc.

References edm::hlt::Exception, kNew, and kReplace.

{
  if ( name == "new" ) return kNew;
  else if ( name == "replace" ) return kReplace;
  else throw cms::Exception("Configuration") 
    << "Unsupported hepMcMode " << name << ": should be either 'new' or 'replace' !!\n" << std::endl;
}

Member Data Documentation

edm::Event* MCParticleReplacer::evt_

private

Definition at line 69 of file MCParticleReplacer.h.

Referenced by call_put(), getStreamID(), and produce().

HepMcMode MCParticleReplacer::hepMcMode_

private

Definition at line 67 of file MCParticleReplacer.h.

Referenced by produce().

ParticleReplacerBase* MCParticleReplacer::replacer_

private

Definition at line 68 of file MCParticleReplacer.h.

Referenced by beginJob(), beginRun(), endJob(), endRun(), MCParticleReplacer(), produce(), and ~MCParticleReplacer().

edm::InputTag MCParticleReplacer::src_

private

Definition at line 65 of file MCParticleReplacer.h.

Referenced by produce().

edm::InputTag MCParticleReplacer::srcHepMC_

private

Definition at line 66 of file MCParticleReplacer.h.

Referenced by produce().

int MCParticleReplacer::verbosity_

private

Definition at line 71 of file MCParticleReplacer.h.

Referenced by MCParticleReplacer().