edm::GeneratorFilter< HAD, DEC > Class Template Reference

#include <GeneratorFilter.h>

Inheritance diagram for edm::GeneratorFilter< HAD, DEC >:

Public Types
typedef DEC	Decayer
typedef HAD	Hadronizer
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
typedef WorkerT< EDFilter >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Member Functions
virtual bool	beginLuminosityBlock (LuminosityBlock &, EventSetup const &)
virtual bool	beginRun (Run &, EventSetup const &)
virtual void	endJob ()
virtual bool	endLuminosityBlock (LuminosityBlock &, EventSetup const &)
virtual bool	endRun (Run &, EventSetup const &)
virtual bool	filter (Event &e, EventSetup const &es)
	GeneratorFilter (ParameterSet const &ps)
virtual void	respondToCloseInputFile (FileBlock const &fb)
virtual void	respondToCloseOutputFiles (FileBlock const &fb)
virtual void	respondToOpenInputFile (FileBlock const &fb)
virtual void	respondToOpenOutputFiles (FileBlock const &fb)
virtual	~GeneratorFilter ()
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
virtual	~EDFilter ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Attributes
Decayer *	decayer_
Hadronizer	hadronizer_
unsigned int	nEventsInLumiBlock_

Additional Inherited Members
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

template<class HAD, class DEC>
class edm::GeneratorFilter< HAD, DEC >

Definition at line 35 of file GeneratorFilter.h.

Member Typedef Documentation

template<class HAD , class DEC >

typedef DEC edm::GeneratorFilter< HAD, DEC >::Decayer

Definition at line 39 of file GeneratorFilter.h.

template<class HAD , class DEC >

typedef HAD edm::GeneratorFilter< HAD, DEC >::Hadronizer

Definition at line 38 of file GeneratorFilter.h.

Constructor & Destructor Documentation

template<class HAD , class DEC >

edm::GeneratorFilter< HAD, DEC >::GeneratorFilter ( ParameterSet const &  ps )

explicit

Definition at line 70 of file GeneratorFilter.h.

References edm::GeneratorFilter< HAD, DEC >::decayer_, edm::ParameterSet::exists(), and edm::ParameterSet::getParameter().

                                                                  :
    EDFilter(),
    hadronizer_(ps),
    decayer_(0),
    nEventsInLumiBlock_(0)
  {
    // TODO:
    // Put the list of types produced by the filters here.
    // The current design calls for:
    //   * GenRunInfoProduct
    //   * HepMCProduct
    //
    // other maybe added as needs be
    //
    
    if ( ps.exists("ExternalDecays") )
    {
       //decayer_ = new gen::ExternalDecayDriver(ps.getParameter<ParameterSet>("ExternalDecays"));
       ParameterSet ps1 = ps.getParameter<ParameterSet>("ExternalDecays");
       decayer_ = new Decayer(ps1);
    }
    
    produces<edm::HepMCProduct>();
    produces<GenEventInfoProduct>();
    produces<GenLumiInfoProduct, edm::InLumi>();
    produces<GenRunInfoProduct, edm::InRun>();
 
  }

template<class HAD , class DEC >

edm::GeneratorFilter< HAD, DEC >::~GeneratorFilter ( )

virtual

Definition at line 100 of file GeneratorFilter.h.

  { if ( decayer_ ) delete decayer_;}

Member Function Documentation

template<class HAD , class DEC >

bool edm::GeneratorFilter< HAD, DEC >::beginLuminosityBlock ( LuminosityBlock &  , EventSetup const &  es  )

virtual

Reimplemented from edm::EDFilter.

Definition at line 214 of file GeneratorFilter.h.

References edm::errors::Configuration, and edm::hlt::Exception.

  {
    nEventsInLumiBlock_ = 0;
    
    if ( !hadronizer_.readSettings(0) )
       throw edm::Exception(errors::Configuration) 
     << "Failed to read settings for the hadronizer "
     << hadronizer_.classname() << " \n";
        
    if ( decayer_ )
    {
       decayer_->init(es);
       if ( !hadronizer_.declareStableParticles( decayer_->operatesOnParticles() ) )
          throw edm::Exception(errors::Configuration)
            << "Failed to declare stable particles in hadronizer "
            << hadronizer_.classname()
            << "\n";
       if ( !hadronizer_.declareSpecialSettings( decayer_->specialSettings() ) )
          throw edm::Exception(errors::Configuration)
            << "Failed to declare special settings in hadronizer "
            << hadronizer_.classname()
            << "\n";
    }

    if ( !hadronizer_.initializeForInternalPartons() )
       throw edm::Exception(errors::Configuration) 
     << "Failed to initialize hadronizer "
     << hadronizer_.classname()
     << " for internal parton generation\n";

    return true;

  }

template<class HAD , class DEC >

bool edm::GeneratorFilter< HAD, DEC >::beginRun ( Run &  , EventSetup const &  es  )

virtual

Reimplemented from edm::EDFilter.

Definition at line 186 of file GeneratorFilter.h.

  {
    
    return true;

  }

template<class HAD , class DEC >

void edm::GeneratorFilter< HAD, DEC >::endJob ( void  )

virtual

Reimplemented from edm::EDFilter.

Definition at line 181 of file GeneratorFilter.h.

  { }

template<class HAD , class DEC >

bool edm::GeneratorFilter< HAD, DEC >::endLuminosityBlock ( LuminosityBlock &  lumi, EventSetup const &    )

virtual

Reimplemented from edm::EDFilter.

Definition at line 250 of file GeneratorFilter.h.

References GenRunInfoProduct::XSec::error(), GenRunInfoProduct::internalXSec(), edm::LuminosityBlock::put(), GenLumiInfoProduct::ProcessInfo::setAccepted(), GenLumiInfoProduct::ProcessInfo::setAcceptedBr(), GenLumiInfoProduct::ProcessInfo::setKilled(), GenLumiInfoProduct::ProcessInfo::setLheXSec(), GenLumiInfoProduct::ProcessInfo::setNPassNeg(), GenLumiInfoProduct::ProcessInfo::setNPassPos(), GenLumiInfoProduct::ProcessInfo::setNTotalNeg(), GenLumiInfoProduct::ProcessInfo::setNTotalPos(), GenLumiInfoProduct::ProcessInfo::setProcess(), GenLumiInfoProduct::ProcessInfo::setSelected(), GenLumiInfoProduct::ProcessInfo::setTried(), groupFilesInBlocks::temp, and GenRunInfoProduct::XSec::value().

  {
    
    hadronizer_.statistics();    
    if ( decayer_ ) decayer_->statistics();

    GenRunInfoProduct genRunInfo = GenRunInfoProduct(hadronizer_.getGenRunInfo());
    std::vector<GenLumiInfoProduct::ProcessInfo> GenLumiProcess;
    GenRunInfoProduct::XSec xsec = genRunInfo.internalXSec();
    GenLumiInfoProduct::ProcessInfo temp;      
    temp.setProcess(0);
    temp.setLheXSec(xsec.value(), xsec.error()); // Pythia gives error of -1
    temp.setNPassPos(nEventsInLumiBlock_);
    temp.setNPassNeg(0);
    temp.setNTotalPos(nEventsInLumiBlock_);
    temp.setNTotalNeg(0);
    temp.setTried(nEventsInLumiBlock_, nEventsInLumiBlock_, nEventsInLumiBlock_);
    temp.setSelected(nEventsInLumiBlock_, nEventsInLumiBlock_, nEventsInLumiBlock_);
    temp.setKilled(nEventsInLumiBlock_, nEventsInLumiBlock_, nEventsInLumiBlock_);
    temp.setAccepted(0,-1,-1);
    temp.setAcceptedBr(0,-1,-1);
    GenLumiProcess.push_back(temp);

    std::auto_ptr<GenLumiInfoProduct> genLumiInfo(new GenLumiInfoProduct());
    genLumiInfo->setHEPIDWTUP(-1);
    genLumiInfo->setProcessInfo( GenLumiProcess );
    lumi.put(genLumiInfo);

    nEventsInLumiBlock_ = 0;
    
    return true;
  }

template<class HAD , class DEC >

bool edm::GeneratorFilter< HAD, DEC >::endRun ( Run &  r, EventSetup const &    )

virtual

Reimplemented from edm::EDFilter.

Definition at line 195 of file GeneratorFilter.h.

References edm::Run::put().

  {
    // If relevant, record the integrated luminosity for this run
    // here.  To do so, we would need a standard function to invoke on
    // the contained hadronizer that would report the integrated
    // luminosity.

    hadronizer_.statistics();
    
    if ( decayer_ ) decayer_->statistics();
    
    std::auto_ptr<GenRunInfoProduct> griproduct(new GenRunInfoProduct(hadronizer_.getGenRunInfo()));
    r.put(griproduct);

    return true;
  }

template<class HAD , class DEC >

bool edm::GeneratorFilter< HAD, DEC >::filter ( Event &  e, EventSetup const &  es  )

virtual

Implements edm::EDFilter.

Definition at line 105 of file GeneratorFilter.h.

References edm::EventID::event(), event(), edm::EventBase::id(), and edm::Event::put().

  {
    //added for selecting/filtering gen events, in the case of hadronizer+externalDecayer
      
    bool passEvtGenSelector = false;
    std::auto_ptr<HepMC::GenEvent> event(0);
   
    while(!passEvtGenSelector)
      {
    event.reset();
    hadronizer_.setEDMEvent(ev);
    
    if ( !hadronizer_.generatePartonsAndHadronize() ) return false;
    
    //  this is "fake" stuff
    // in principle, decays are done as part of full event generation,
    // except for particles that are marked as to be kept stable
    // but we currently keep in it the design, because we might want
    // to use such feature for other applications
    //
    if ( !hadronizer_.decay() ) return false;
    
    event = std::auto_ptr<HepMC::GenEvent>(hadronizer_.getGenEvent());
    if ( !event.get() ) return false; 
    
    // The external decay driver is being added to the system,
    // it should be called here
    //
    if ( decayer_ ) 
      {
        event.reset( decayer_->decay( event.get() ) );
      }
    if ( !event.get() ) return false;
    
    passEvtGenSelector = hadronizer_.select( event.get() );
    
      }
    // check and perform if there're any unstable particles after 
    // running external decay packages
    //
    // fisrt of all, put back modified event tree (after external decay)
    //
    hadronizer_.resetEvent( event.release() );
    
    //
    // now run residual decays
    //
    if ( !hadronizer_.residualDecay() ) return false;
        
    hadronizer_.finalizeEvent();
    
    event.reset( hadronizer_.getGenEvent() );
    if ( !event.get() ) return false;
    
    event->set_event_number( ev.id().event() );
    
    //
    // tutto bene - finally, form up EDM products !
    //
    std::auto_ptr<GenEventInfoProduct> genEventInfo(hadronizer_.getGenEventInfo());
    if (!genEventInfo.get())
      { 
    // create GenEventInfoProduct from HepMC event in case hadronizer didn't provide one
    genEventInfo.reset(new GenEventInfoProduct(event.get()));
      }
    ev.put(genEventInfo);
   
    std::auto_ptr<HepMCProduct> bare_product(new HepMCProduct());
    bare_product->addHepMCData( event.release() );
    ev.put(bare_product);
    nEventsInLumiBlock_ ++;
    return true;
  }

template<class HAD , class DEC >

void edm::GeneratorFilter< HAD, DEC >::respondToCloseInputFile ( FileBlock const &  fb )

virtual

Reimplemented from edm::EDFilter.

Definition at line 290 of file GeneratorFilter.h.

  { }

template<class HAD , class DEC >

void edm::GeneratorFilter< HAD, DEC >::respondToCloseOutputFiles ( FileBlock const &  fb )

virtual

Reimplemented from edm::EDFilter.

Definition at line 300 of file GeneratorFilter.h.

  { }

template<class HAD , class DEC >

void edm::GeneratorFilter< HAD, DEC >::respondToOpenInputFile ( FileBlock const &  fb )

virtual

Reimplemented from edm::EDFilter.

Definition at line 285 of file GeneratorFilter.h.

  { }

template<class HAD , class DEC >

void edm::GeneratorFilter< HAD, DEC >::respondToOpenOutputFiles ( FileBlock const &  fb )

virtual

Reimplemented from edm::EDFilter.

Definition at line 295 of file GeneratorFilter.h.

  { }

Member Data Documentation

template<class HAD , class DEC >

Decayer* edm::GeneratorFilter< HAD, DEC >::decayer_

private

Definition at line 61 of file GeneratorFilter.h.

Referenced by edm::GeneratorFilter< HAD, DEC >::GeneratorFilter().

template<class HAD , class DEC >

Hadronizer edm::GeneratorFilter< HAD, DEC >::hadronizer_

private

Definition at line 59 of file GeneratorFilter.h.

template<class HAD , class DEC >

unsigned int edm::GeneratorFilter< HAD, DEC >::nEventsInLumiBlock_

private

Definition at line 62 of file GeneratorFilter.h.