gen::ExternalDecayDriver Class Reference

#include <ExternalDecayDriver.h>

Public Member Functions
HepMC::GenEvent *	decay (HepMC::GenEvent *)
	ExternalDecayDriver (const edm::ParameterSet &)
void	init (const edm::EventSetup &)
const std::vector< int > &	operatesOnParticles ()
const std::vector< std::string > &	specialSettings ()
void	statistics () const
	~ExternalDecayDriver ()

Private Attributes
EvtGenInterfaceBase *	fEvtGenInterface
bool	fIsInitialized
std::vector< int >	fPDGs
PhotosInterfaceBase *	fPhotosInterface
std::vector< std::string >	fSpecialSettings
TauolaInterfaceBase *	fTauolaInterface
bool	hastauolapp

Detailed Description

Definition at line 18 of file ExternalDecayDriver.h.

Constructor & Destructor Documentation

ExternalDecayDriver::ExternalDecayDriver

(

const edm::ParameterSet &

pset

)

Definition at line 21 of file ExternalDecayDriver.cc.

References gen::PhotosInterfaceBase::avoidTauLeptonicDecays(), gen::PhotosInterfaceBase::configureOnlyFor(), decayRandomEngine, edm::hlt::Exception, fEvtGenInterface, fPhotosInterface, fTauolaInterface, reco::get(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hastauolapp, gen::TauolaInterfaceBase::SetDecayRandomEngine(), gen::PhotosInterfaceBase::SetDecayRandomEngine(), and gen::EvtGenInterfaceBase::SetPhotosDecayRandomEngine().

  : hastauolapp(false),
     fIsInitialized(false),
     fTauolaInterface(0),
     fEvtGenInterface(0),
     fPhotosInterface(0)
{
  
  std::vector<std::string> extGenNames = pset.getParameter< std::vector<std::string> >("parameterSets");
  
  Service<RandomNumberGenerator> rng;
  if(!rng.isAvailable()) {
    throw cms::Exception("Configuration")
      << "The RandomNumberProducer module requires the RandomNumberGeneratorService\n"
      "which appears to be absent.  Please add that service to your configuration\n"
      "or remove the modules that require it." << std::endl;
  } 
  decayRandomEngine = &rng->getEngine();   
  
  for (unsigned int ip=0; ip<extGenNames.size(); ++ip ){
    std::string curSet = extGenNames[ip];
    if ( curSet == "EvtGen" || curSet == "EvtGenLHC91"){
      fEvtGenInterface = (EvtGenInterfaceBase*)(EvtGenFactory::get()->create("EvtGenLHC91", pset.getUntrackedParameter< ParameterSet >(curSet)));
      fEvtGenInterface->SetPhotosDecayRandomEngine(decayRandomEngine);
    }
    if ( curSet == "Tauola" || curSet == "Tauola271215" ){
      fTauolaInterface = (TauolaInterfaceBase*)(TauolaFactory::get()->create("Tauola271215", pset.getUntrackedParameter< ParameterSet >(curSet)));
      fTauolaInterface->SetDecayRandomEngine(decayRandomEngine);
    }
    if(curSet =="Tauolapp111a"){
      fTauolaInterface = (TauolaInterfaceBase*)(TauolaFactory::get()->create("Tauolapp111a", pset.getUntrackedParameter< ParameterSet >(curSet)));
      fPhotosInterface = (PhotosInterfaceBase*)(PhotosFactory::get()->create("Photos2155", pset.getUntrackedParameter< ParameterSet >(curSet)));
      fPhotosInterface->SetDecayRandomEngine(decayRandomEngine);
      fPhotosInterface->configureOnlyFor(15); 
      fPhotosInterface->avoidTauLeptonicDecays();
      hastauolapp=true;
    }
    if ( curSet == "Photos" || curSet == "Photos2155Legacy"){
      if ( !fPhotosInterface ){
    fPhotosInterface = (PhotosInterfaceBase*)(PhotosFactory::get()->create("Photos2155Legacy", pset.getUntrackedParameter< ParameterSet >(curSet)));
    fPhotosInterface->SetDecayRandomEngine(decayRandomEngine);
      }
    }
    if (curSet == "Photos2155"){
      if ( !fPhotosInterface ){
    fPhotosInterface = (PhotosInterfaceBase*)(PhotosFactory::get()->create("Photos2155", pset.getUntrackedParameter< ParameterSet>(curSet)));
    fPhotosInterface->SetDecayRandomEngine(decayRandomEngine);
      }
    }
  }
}

ExternalDecayDriver::~ExternalDecayDriver

(

)

Definition at line 73 of file ExternalDecayDriver.cc.

References fEvtGenInterface, fPhotosInterface, and fTauolaInterface.

                                         {
  if ( fEvtGenInterface ) delete fEvtGenInterface;
  if ( fTauolaInterface ) delete fTauolaInterface;
  if ( fPhotosInterface ) delete fPhotosInterface;
}

Member Function Documentation

HepMC::GenEvent * ExternalDecayDriver::decay

(

HepMC::GenEvent *

evt

)

Definition at line 79 of file ExternalDecayDriver.cc.

References gen::PhotosInterfaceBase::apply(), gen::TauolaInterfaceBase::decay(), gen::EvtGenInterfaceBase::decay(), fEvtGenInterface, fIsInitialized, fPhotosInterface, and fTauolaInterface.

                                                             {
  
  if ( !fIsInitialized ) return evt;
  
  if ( fEvtGenInterface ){  
    evt = fEvtGenInterface->decay( evt ); 
    if ( !evt ) return 0;
  }
  
  if ( fTauolaInterface ){
    evt = fTauolaInterface->decay( evt ); 
    if ( !evt ) return 0;
  }
  
  if ( fPhotosInterface ){
    evt = fPhotosInterface->apply( evt );
    if ( !evt ) return 0;
  }

  return evt;
}

void ExternalDecayDriver::init

(

const edm::EventSetup &

es

)

Definition at line 101 of file ExternalDecayDriver.cc.

References fEvtGenInterface, fIsInitialized, fPDGs, fPhotosInterface, fSpecialSettings, fTauolaInterface, hastauolapp, i, gen::EvtGenInterfaceBase::init(), gen::TauolaInterfaceBase::init(), gen::PhotosInterfaceBase::init(), gen::TauolaInterfaceBase::operatesOnParticles(), gen::EvtGenInterfaceBase::operatesOnParticles(), and gen::PhotosInterfaceBase::specialSettings().

                                                       {

   if ( fIsInitialized ) return;
   
   if ( fTauolaInterface ) {
     fTauolaInterface->init( es );
     for ( std::vector<int>::const_iterator i=fTauolaInterface->operatesOnParticles().begin();
       i!=fTauolaInterface->operatesOnParticles().end(); i++ ) 
       fPDGs.push_back( *i );
   }
   if ( fEvtGenInterface ){
     fEvtGenInterface->init();
     for ( std::vector<int>::const_iterator i=fEvtGenInterface->operatesOnParticles().begin();
       i!=fEvtGenInterface->operatesOnParticles().end(); i++ )
       fPDGs.push_back( *i );
   }
   
   
   if( fPhotosInterface){
     fPhotosInterface->init();
     if(hastauolapp){
       if ( fPhotosInterface ){
     for ( unsigned int iss=0; iss<fPhotosInterface->specialSettings().size(); iss++ ){
       fSpecialSettings.push_back( fPhotosInterface->specialSettings()[iss]);
     }
       }
     }
   }
   // now put in hack for TauolaFortran settings
   if(!hastauolapp){
     if ( fPhotosInterface ){
       fSpecialSettings.push_back( "QED-brem-off:all" );
     }
     if ( fTauolaInterface ){
       // override !
       fSpecialSettings.clear();
       fSpecialSettings.push_back( "QED-brem-off:15" );
     }
   }
   fIsInitialized = true;
   
   return;
}

const std::vector<int>& gen::ExternalDecayDriver::operatesOnParticles ( )

inline

Definition at line 26 of file ExternalDecayDriver.h.

References fPDGs.

{ return fPDGs; }

const std::vector<std::string>& gen::ExternalDecayDriver::specialSettings ( )

inline

Definition at line 27 of file ExternalDecayDriver.h.

References fSpecialSettings.

{ return fSpecialSettings; }

void ExternalDecayDriver::statistics

(

)

const

Definition at line 145 of file ExternalDecayDriver.cc.

References fTauolaInterface, and gen::TauolaInterfaceBase::statistics().

{
   if ( fTauolaInterface ) fTauolaInterface->statistics();
   // similar for EvtGen and/or Photos, if needs be
   return;
}

Member Data Documentation

EvtGenInterfaceBase* gen::ExternalDecayDriver::fEvtGenInterface

private

Definition at line 35 of file ExternalDecayDriver.h.

Referenced by decay(), ExternalDecayDriver(), init(), and ~ExternalDecayDriver().

bool gen::ExternalDecayDriver::fIsInitialized

private

Definition at line 33 of file ExternalDecayDriver.h.

Referenced by decay(), and init().

std::vector<int> gen::ExternalDecayDriver::fPDGs

private

Definition at line 37 of file ExternalDecayDriver.h.

Referenced by init(), and operatesOnParticles().

PhotosInterfaceBase* gen::ExternalDecayDriver::fPhotosInterface

private

Definition at line 36 of file ExternalDecayDriver.h.

Referenced by decay(), ExternalDecayDriver(), init(), and ~ExternalDecayDriver().

std::vector<std::string> gen::ExternalDecayDriver::fSpecialSettings

private

Definition at line 38 of file ExternalDecayDriver.h.

Referenced by init(), and specialSettings().

TauolaInterfaceBase* gen::ExternalDecayDriver::fTauolaInterface

private

Definition at line 34 of file ExternalDecayDriver.h.

Referenced by decay(), ExternalDecayDriver(), init(), statistics(), and ~ExternalDecayDriver().

bool gen::ExternalDecayDriver::hastauolapp

private

Definition at line 32 of file ExternalDecayDriver.h.

Referenced by ExternalDecayDriver(), and init().