Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F19274573
D135.1759140957.diff
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Award Token
Flag For Later
Size
6 KB
Referenced Files
None
Subscribers
None
D135.1759140957.diff
View Options
diff --git a/EvtGenExternal/EvtExternalGenFactory.hh b/EvtGenExternal/EvtExternalGenFactory.hh
--- a/EvtGenExternal/EvtExternalGenFactory.hh
+++ b/EvtGenExternal/EvtExternalGenFactory.hh
@@ -46,7 +46,9 @@
void definePythiaGenerator( std::string xmlDir, bool convertPhysCodes,
bool useEvtGenRandom = true );
void defineTauolaGenerator( bool useEvtGenRandom = true,
- bool seedTauolaFortran = true );
+ bool seedTauolaFortran = true,
+ bool useTauolaRadiation = false,
+ double infraredCutOffTauola = 1.0e-7 );
//methods to add configuration commands to the pythia generators
//void addPythiaCommand( std::string generator, std::string module, std::string param, std::string value);
diff --git a/EvtGenExternal/EvtExternalGenList.hh b/EvtGenExternal/EvtExternalGenList.hh
--- a/EvtGenExternal/EvtExternalGenList.hh
+++ b/EvtGenExternal/EvtExternalGenList.hh
@@ -35,7 +35,9 @@
std::string pythiaXmlDir = "",
std::string photonType = "gamma",
bool useEvtGenRandom = true,
- bool seedTauolaFortran = true );
+ bool seedTauolaFortran = true,
+ bool useTauolaRadiation = false,
+ double infraredCutOffTauola = 1.0e-7 );
virtual ~EvtExternalGenList();
diff --git a/EvtGenExternal/EvtTauolaEngine.hh b/EvtGenExternal/EvtTauolaEngine.hh
--- a/EvtGenExternal/EvtTauolaEngine.hh
+++ b/EvtGenExternal/EvtTauolaEngine.hh
@@ -49,7 +49,9 @@
class EvtTauolaEngine : public EvtAbsExternalGen {
public:
- EvtTauolaEngine( bool useEvtGenRandom = true, bool seedTauolaFortran = true );
+ EvtTauolaEngine( bool useEvtGenRandom = true, bool seedTauolaFortran = true,
+ bool useTauolaRadiation = false,
+ double infraredCutOffTauola = 1.0e-7 );
bool doDecay( EvtParticle* theMother ) override;
@@ -68,6 +70,10 @@
bool m_useEvtGenRandom{ true };
bool m_seedTauolaFortran{ true };
+ // Arguments to steer radiation in leptonic tau decays
+ bool m_useTauolaRadiation{ false };
+ double m_infraredCutOffTauola{ 1.0e-7 };
+
// PDG standard code integer ID for tau particle
static constexpr int m_tauPDG{ 15 };
// Number of possible decay modes in Tauola
diff --git a/src/EvtGenExternal/EvtExternalGenFactory.cpp b/src/EvtGenExternal/EvtExternalGenFactory.cpp
--- a/src/EvtGenExternal/EvtExternalGenFactory.cpp
+++ b/src/EvtGenExternal/EvtExternalGenFactory.cpp
@@ -73,12 +73,15 @@
#ifdef EVTGEN_TAUOLA
void EvtExternalGenFactory::defineTauolaGenerator( bool useEvtGenRandom,
- bool seedTauolaFortran )
+ bool seedTauolaFortran,
+ bool useTauolaRadiation,
+ double infraredCutOffTauola )
{
EvtGenReport( EVTGEN_INFO, "EvtGen" ) << "Defining EvtTauolaEngine." << endl;
- m_extGenMap[GenId::TauolaGenId] =
- std::make_unique<EvtTauolaEngine>( useEvtGenRandom, seedTauolaFortran );
+ m_extGenMap[GenId::TauolaGenId] = std::make_unique<EvtTauolaEngine>(
+ useEvtGenRandom, seedTauolaFortran, useTauolaRadiation,
+ infraredCutOffTauola );
}
#else
void EvtExternalGenFactory::defineTauolaGenerator( bool, bool )
diff --git a/src/EvtGenExternal/EvtExternalGenList.cpp b/src/EvtGenExternal/EvtExternalGenList.cpp
--- a/src/EvtGenExternal/EvtExternalGenList.cpp
+++ b/src/EvtGenExternal/EvtExternalGenList.cpp
@@ -28,11 +28,10 @@
#include "EvtGenExternal/EvtSherpaPhotons.hh"
#include "EvtGenExternal/EvtTauola.hh"
-EvtExternalGenList::EvtExternalGenList( bool convertPythiaCodes,
- std::string pythiaXmlDir,
- std::string photonType,
- bool useEvtGenRandom,
- bool seedTauolaFortran ) :
+EvtExternalGenList::EvtExternalGenList(
+ bool convertPythiaCodes, std::string pythiaXmlDir, std::string photonType,
+ bool useEvtGenRandom, bool seedTauolaFortran, bool useTauolaRadiation,
+ double infraredCutOffTauola ) :
m_photonType{ photonType }, m_useEvtGenRandom{ useEvtGenRandom }
{
// Instantiate the external generator factory
@@ -51,7 +50,8 @@
extFactory.definePythiaGenerator( pythiaXmlDir, convertPythiaCodes,
useEvtGenRandom );
- extFactory.defineTauolaGenerator( useEvtGenRandom, seedTauolaFortran );
+ extFactory.defineTauolaGenerator( useEvtGenRandom, seedTauolaFortran,
+ useTauolaRadiation, infraredCutOffTauola );
}
EvtExternalGenList::~EvtExternalGenList()
diff --git a/src/EvtGenExternal/EvtTauolaEngine.cpp b/src/EvtGenExternal/EvtTauolaEngine.cpp
--- a/src/EvtGenExternal/EvtTauolaEngine.cpp
+++ b/src/EvtGenExternal/EvtTauolaEngine.cpp
@@ -49,8 +49,13 @@
bool EvtTauolaEngine::m_initialised = false;
std::mutex EvtTauolaEngine::m_tauola_mutex;
-EvtTauolaEngine::EvtTauolaEngine( bool useEvtGenRandom, bool seedTauolaFortran ) :
- m_useEvtGenRandom{ useEvtGenRandom }, m_seedTauolaFortran{ seedTauolaFortran }
+EvtTauolaEngine::EvtTauolaEngine( bool useEvtGenRandom, bool seedTauolaFortran,
+ bool useTauolaRadiation,
+ double infraredCutOffTauola ) :
+ m_useEvtGenRandom{ useEvtGenRandom },
+ m_seedTauolaFortran{ seedTauolaFortran },
+ m_useTauolaRadiation{ useTauolaRadiation },
+ m_infraredCutOffTauola{ infraredCutOffTauola }
{
}
@@ -92,6 +97,12 @@
Tauolapp::Tauola::setRandomGenerator( EvtRandom::Flat );
}
+ // Switch for radiation in leptonic tau decays
+ Tauolapp::Tauola::setRadiation( m_useTauolaRadiation );
+
+ // Cut-Off parameter of radition. Above that value photon is explicitly generated.
+ Tauolapp::Tauola::setRadiationCutOff( m_infraredCutOffTauola );
+
// Use the BaBar-tuned chiral current calculations by default. Can be changed using the
// TauolaCurrentOption keyword in decay files
Tauolapp::Tauola::setNewCurrents( 1 );
File Metadata
Details
Attached
Mime Type
text/plain
Expires
Mon, Sep 29, 11:15 AM (21 h, 15 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6561980
Default Alt Text
D135.1759140957.diff (6 KB)
Attached To
Mode
D135: Control radiation corrections in TAUOLA
Attached
Detach File
Event Timeline
Log In to Comment