9#include <analysis/modules/TauDecayMarker/TauDecayMarkerModule.h>
11#include <framework/logging/Logger.h>
12#include <framework/gearbox/Const.h>
27 idOfTauMinus(-1), m_pmode(-2), m_mmode(-2), m_pprong(0), m_mprong(0)
30 setDescription(
"Module to identify generated tau pair decays, using MCParticle information. Each tau lepton decay channel "
31 "is numbered following the order in the default KKMC decay table. Using this module, "
32 "the channel number will be stored in the variables `tauPlusMCMode`, and `tauMinusMCMode`. "
33 "Further details and usage can be found at `TauDecayMCModes`. ");
58 B2INFO(
"Decay ID: " <<
m_pmode <<
" (tau+), " <<
m_mmode <<
" (tau-)." <<
84 if (p.getStatus() == 1 && p.getPDG() == 15) {
88 if (p.getStatus() == 1 && p.getPDG() == -15) {
100 if (s == 0 || !
tauPair)
return -1;
107 for (
int i = p.getFirstDaughter(); i <= p.getLastDaughter(); ++i) {
109 if (abs(d.getPDG()) == 24)
110 ret += d.getLastDaughter() - d.getFirstDaughter() + 1;
114 for (
int i = p.getFirstDaughter(); i <= p.getLastDaughter(); ++i) {
116 int pdg = d.getPDG();
117 if (pdg ==
id || (sign == 0 && abs(pdg) == abs(
id))) ret++;
118 if (abs(pdg) == 24) {
119 for (
int j = d.getFirstDaughter(); j <= d.getLastDaughter(); ++j) {
121 int pdg2 = e.getPDG();
133 if (s == 0 || !
tauPair)
return -1;
140 for (
int i = p.getFirstDaughter(); i <= p.getLastDaughter(); ++i) {
142 if (abs(d.getPDG()) == 24) {
143 for (
int j = d.getFirstDaughter(); j <= d.getLastDaughter(); ++j) {
150 for (
int i = p.getFirstDaughter(); i <= p.getLastDaughter(); ++i) {
152 int pdg = d.getPDG();
153 if (abs(pdg) == 24) {
154 for (
int j = d.getFirstDaughter(); j <= d.getLastDaughter(); ++j) {
156 int pdg2 = e.getPDG();
158 (sign == 0 && abs(pdg2) == abs(
id))) ret++;
160 }
else if (pdg ==
id || (sign == 0 && abs(pdg) == abs(
id))) ret++;
564 const vector<MCParticle*> daughters = p.getDaughters();
565 if (daughters.empty())
return ret;
572 if (isChargedFinalState) ret++;
int getPDGCode() const
PDG code.
static const ParticleType photon
photon particle
A Class to store the Monte Carlo particle information.
@ c_PrimaryParticle
bit 0: Particle is primary particle.
void setDescription(const std::string &description)
Sets the description of the module.
void setPropertyFlags(unsigned int propertyFlags)
Sets the flags for the module properties.
@ c_ParallelProcessingCertified
This module can be run in parallel processing mode safely (All I/O must be done through the data stor...
int getEntries() const
Get the number of objects in the array.
Int_t m_mprong
Prong of the decay channel of negative tau.
const int finalStatePDGs[5]
PDG codes accepted as charged final state particles in generation: {e, mu, pi, K, p}.
int numOfTauPlus
Number of positive tau leptons in the event.
Int_t m_pprong
Prong of the decay channel of positive tau.
void IdentifyTauPair()
Identifies if the event is a generated tau pair.
virtual void initialize() override
Initializes the module.
virtual void event() override
Method is called for each event.
int getNumDaughterOfTau(int s=0, int id=0, int sign=0)
Count the number of daughers of the generated tau.
Int_t m_mmode
ID of the decay channel of negative tau.
int getNumDaughterOfTauExceptGamma(int s=0, int id=0, int sign=0)
Count the number of daughers of the generated tau except gammas.
bool tauPair
True if the generated event is a tau pair event.
int idOfTauMinus
Index of the generated negative tau.
Int_t m_pmode
ID of the decay channel of positive tau.
int numOfTauMinus
Number of negative tau leptons in the event.
int getProngOfDecay(const MCParticle &mc)
Get the prong of the generated tau decay.
int idOfTauPlus
Index of the generated positive tau.
StoreObjPtr< TauPairDecay > m_tauDecay
pointer to tau pair decay objects
int getDecayChannelOfTau(int s=0)
Gets the id of the decay channel.
bool m_printDecayInfo
If true, prints ID and prong of each tau lepton in the event.
StoreArray< MCParticle > m_MCParticles
StoreArray of MCParticles.
TauDecayMarkerModule()
Constructor: Sets the description, the properties and the parameters of the module.
void addParam(const std::string &name, T ¶mVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module.
#define REG_MODULE(moduleName)
Register the given module (without 'Module' suffix) with the framework.
Abstract base class for different kinds of events.
1.9.6