#include <IEvtGen.h>

Public Member Functions
	IEvtGen (const std::string &inputFileName, const std::vector< int > &pattern, TRandom *rand, const double &unitConversion)

virtual	~IEvtGen ()

void	GenerateEvent (std::vector< std::vector< double >> &daughters, double &weight_gen)

std::vector< double >	GetDaughter (const IDalitzEvent *const event, const unsigned int &daughter_id)

Private Attributes
const double	m_unitConversion

SignalGenerator *	m_sg

Detailed Description

Definition at line 14 of file IEvtGen.h.

Constructor & Destructor Documentation

◆ IEvtGen()

IEvtGen::IEvtGen	(	const std::string &	inputFileName,
		const std::vector< int > &	pattern,
		TRandom *	rand,
		const double &	unitConversion
	)

Interface to EvtGen Requires Mint steering file as argument The pattern should be deduced in the EvtGen model The random number generator should point to the EvtGen implementation Requires unit conversion rate between Mint and EvtGen as argument Mint 4-vectors will be divided by this factor and returned to EvtGen

Definition at line 5 of file IEvtGen.cpp.

   : m_unitConversion(unitConversion)
 {
   NamedParameterBase::setDefaultInputFile(inputFileName);
 
   DalitzEventPattern pdg(pattern);
   if( pattern[0] < 0 )
     pdg = pdg.makeCPConjugate();
 
   std::cout << "Mother Particle " << pattern[0] << std::endl
         << " got event pattern: " << pdg << std::endl;
 
   m_sg = new SignalGenerator(pdg, rand);
 }

◆ ~IEvtGen()

virtual MINT::IEvtGen::~IEvtGen ( )

inlinevirtual

Definition at line 29 of file IEvtGen.h.

30 { delete m_sg; }

MINT::IEvtGen::m_sg

SignalGenerator * m_sg

Definition: IEvtGen.h:47

Member Function Documentation

◆ GenerateEvent()

void IEvtGen::GenerateEvent	(	std::vector< std::vector< double >> &	daughters,
		double &	weight_gen
	)

Generate single event in Mint

Definition at line 22 of file IEvtGen.cpp.

 {
   MINT::counted_ptr<IDalitzEvent> newEvt = m_sg->newEvent();
   const IDalitzEvent* const event = newEvt.get();
   //const IDalitzEvent* const event = m_sg->newEvent().get();
 
   //0th entry is the mother particle
   for( unsigned int i=1;
        i<static_cast<unsigned int>(event->eventPattern().numDaughters()+1);
        ++i ){
     const std::vector<double> daughter = GetDaughter(event, i);
     daughters.push_back(daughter);
   }
 
   weight_gen = event->getGeneratorPdfRelativeToPhaseSpace();
 }

◆ GetDaughter()

std::vector< double > IEvtGen::GetDaughter	(	const IDalitzEvent *const	event,
		const unsigned int &	daughter_id
	)

Intermediary to go between Mint and EvtGen 4-vector types

Definition at line 40 of file IEvtGen.cpp.

 {
   std::vector<double> daughter;
   daughter.push_back(event->p(daughter_id).T()/m_unitConversion);
   daughter.push_back(event->p(daughter_id).X()/m_unitConversion);
   daughter.push_back(event->p(daughter_id).Y()/m_unitConversion);
   daughter.push_back(event->p(daughter_id).Z()/m_unitConversion);
 
   return daughter;
 }

Member Data Documentation

◆ m_sg

SignalGenerator* MINT::IEvtGen::m_sg

private

Definition at line 47 of file IEvtGen.h.

◆ m_unitConversion

const double MINT::IEvtGen::m_unitConversion