Page MenuHomeHEPForge

D135.id575.diff
No OneTemporary

D135.id575.diff

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/History.md b/History.md
--- a/History.md
+++ b/History.md
@@ -11,6 +11,10 @@
===
## R03-00-00-betaX
+7 Dec 2025 Fernando Abudinen
+* D135: Control radiation corrections in TAUOLA
+ - Introduced arguments to switch on/off TAUOLA FSR for leptonic decays and control photon cut-off.
+
3 Dec 2024 Fernando Abudinen
* D134: Add support for Sherpa 3
- Modified EvtSherpaPhotons interface to support both Sherpa 2 and 3 releases
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 radiation. 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

Mime Type
text/plain
Expires
Sun, Jan 19, 9:16 PM (4 h, 4 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4211200
Default Alt Text
D135.id575.diff (6 KB)

Event Timeline