9#include <analysis/modules/TauDecayMode/TauDecayModeModule.h>
10#include <framework/logging/Logger.h>
11#include <framework/particledb/EvtGenDatabasePDG.h>
14#include <TParticlePDG.h>
17#include <Math/LorentzRotation.h>
18#include <Math/Boost.h>
22std::map<std::string, int> make_map(
const std::string& file,
int chg)
27 B2INFO(
"Missing input mapping file: use mp_file_minus=basf2.find_file('data/analysis/modules/TauDecayMode/map_tauminus.txt') TauDecayMode.param('file_minus', mp_file_minus) to classify with default mapping.");
29 B2INFO(
"Missing input mapping file: use mp_file_plus=basf2.find_file('data/analysis/modules/TauDecayMode/map_tauplus.txt') TauDecayMode.param('file_plus', mp_file_plus) to classify with default mapping.");
40 std::map <std::string, int> map_tau;
43 istringstream ss(line);
66 setDescription(
"Module to identify generated tau pair decays, using MCParticle information.\n"
67 "By default, each tau decay is numbered as TauolaBelle2DecayMode [Ref: BELLE2-NOTE-PH-2020-055]");
69 addParam(
"printmode",
m_printmode,
"Printout more information from each event", std::string(
"default"));
70 addParam(
"file_minus",
m_file_minus,
"Path for an alternative mapping for tau- decays", std::string(
""));
71 addParam(
"file_plus",
m_file_plus,
"Path for an alternative mapping for tau+ decays", std::string(
""));
198 map<int, std::vector<int>> map_vec;
241 map_vec[9010221] =
vec_f0;
247 bool TauMinusIsLeptonicDecay =
false;
248 bool TauPlusIsLeptonicDecay =
false;
250 bool elecFirst =
true;
251 bool muonFirst =
true;
253 bool isPiPizGamTauMinusFirst =
true;
254 bool isPiPizGamTauPlusFirst =
true;
256 bool isLFVTauMinus2BodyDecayFirst =
true;
257 bool isLFVTauPlus2BodyDecayFirst =
true;
259 bool isChargedRhoFromTauMinusFirst =
true;
260 bool isChargedRhoFromTauPlusFirst =
true;
262 bool isChargedA1FromTauMinusFirst =
true;
263 bool isChargedA1FromTauPlusFirst =
true;
265 bool isEtaPPGFromTauMinusFirst =
true;
266 bool isEtaPPGFromTauPlusFirst =
true;
268 bool isOmegaPizGamFromTauMinusFirst =
true;
269 bool isOmegaPizGamFromTauPlusFirst =
true;
271 bool isEtaPizPizPizFromTauMinusFirst =
true;
272 bool isEtaPizPizPizFromTauPlusFirst =
true;
277 bool isOmegaPimPipFromTauMinusFirst =
true;
278 bool isOmegaPimPipFromTauPlusFirst =
true;
284 for (
int i = 0; i <
MCParticles.getEntries(); i++) {
288 int pdgid = p.getPDG();
290 if (pdgid == -
Const::Klong.getPDGCode()) pdgid = -pdgid;
293 bool accept_photon =
false;
297 if (p.isInitial())
continue;
300 if (abs(pdgid) == 12 || abs(pdgid) == 14) {
301 int leptonicdecay = 0;
303 const vector<MCParticle*> daughters = mother->
getDaughters();
305 int dauid = abs(d->getPDG());
306 if (dauid == 11 || dauid == 12 || dauid == 13 || dauid == 14 || dauid == 16) {
310 if (leptonicdecay == 3) {
311 if (mother->
getPDG() == 15) TauMinusIsLeptonicDecay =
true;
312 if (mother->
getPDG() == -15) TauPlusIsLeptonicDecay =
true;
319 if (not mother)
continue;
320 const vector<MCParticle*> daughters = mother->
getDaughters();
323 stringstream elec_ss;
326 elec_ss << d->getPDG() <<
" ";
330 if (nElMinus == 1 && nElPlus == 1) {
331 B2DEBUG(19,
"Mother of elec pair is = " << mother->
getPDG() <<
" which has daughters : " << elec_ss.str());
335 else if (pdgid ==
Const::muon.getPDGCode() && muonFirst) {
338 if (not mother)
continue;
339 const vector<MCParticle*> daughters = mother->
getDaughters();
342 stringstream muon_ss;
345 muon_ss << d->getPDG() <<
" ";
346 if (d->getPDG() ==
Const::muon.getPDGCode()) nMuMinus++;
347 if (d->getPDG() == -
Const::muon.getPDGCode()) nMuPlus++;
349 if (nMuMinus == 1 && nMuPlus == 1) {
350 B2DEBUG(19,
"Mother of muon pair is = " << mother->
getPDG() <<
" which has daughters : " << muon_ss.str());
364 if (not mother)
continue;
365 int mothid = abs(mother->
getPDG());
368 bool isRadiationfromFinalStateChargedParticle =
false;
375 isRadiationfromFinalStateChargedParticle =
true;
379 bool isRadiationFromChargedRho =
false;
382 if (chg < 0 && isChargedRhoFromTauMinusFirst) {
383 isChargedRhoFromTauMinusFirst =
false;
384 isRadiationFromChargedRho =
true;
386 if (chg > 0 && isChargedRhoFromTauPlusFirst) {
387 isChargedRhoFromTauPlusFirst =
false;
388 isRadiationFromChargedRho =
true;
393 bool isRadiationFromChargedA1 =
false;
394 if (mothid == 20213) {
396 if (chg < 0 && isChargedA1FromTauMinusFirst) {
397 isChargedA1FromTauMinusFirst =
false;
398 isRadiationFromChargedA1 =
true;
400 if (chg > 0 && isChargedA1FromTauPlusFirst) {
401 isChargedA1FromTauPlusFirst =
false;
402 isRadiationFromChargedA1 =
true;
407 bool isRadiationfromIntermediateWBoson =
false;
408 if (mothid == 24) isRadiationfromIntermediateWBoson =
true;
413 bool isPiPizGam =
false;
414 if (isRadiationfromIntermediateWBoson) {
415 const vector<MCParticle*> daughters = mother->
getDaughters();
419 int nOtherSisters = 0;
423 if (abs(d->getPDG()) ==
Const::pion.getPDGCode()) {
425 }
else if (abs(d->getPDG()) ==
Const::pi0.getPDGCode()) {
432 if (nTotSisters >= 3 && nPiSisters == 1 && nPizSisters == 1 && nOtherSisters == 0) {
434 if (chg < 0 && isPiPizGamTauMinusFirst) {
435 isPiPizGamTauMinusFirst =
false;
438 if (chg > 0 && isPiPizGamTauPlusFirst) {
439 isPiPizGamTauPlusFirst =
false;
446 bool isRadiationfromTau =
false;
447 if (mothid == 15) isRadiationfromTau =
true;
450 if (isRadiationfromTau) {
457 bool isLFVTau2BodyDecay =
false;
458 if (isRadiationfromTau) {
459 bool hasNeutrinoAsSister =
false;
460 int numChargedSister = 0;
461 int numNeutralNonNeutrinoNonPhotonSister = 0;
462 const vector<MCParticle*> daughters = mother->
getDaughters();
466 if (hasNeutrinoAsSister)
break;
468 if (!hasNeutrinoAsSister) {
472 if (isChargedFinalState) {
475 numNeutralNonNeutrinoNonPhotonSister++;
478 if (numChargedSister == 1 && numNeutralNonNeutrinoNonPhotonSister == 0) {
479 if (mother->
getPDG() == 15 && isLFVTauMinus2BodyDecayFirst) {
480 isLFVTauMinus2BodyDecayFirst =
false;
481 isLFVTau2BodyDecay =
true;
483 if (mother->
getPDG() == -15 && isLFVTauPlus2BodyDecayFirst) {
484 isLFVTauPlus2BodyDecayFirst =
false;
485 isLFVTau2BodyDecay =
true;
489 B2DEBUG(19,
"hasNeutrinoAsSister = " << hasNeutrinoAsSister
490 <<
" numChargedSister = " << numChargedSister
491 <<
" numNeutralNonNeutrinoNonPhotonSister = " << numNeutralNonNeutrinoNonPhotonSister
492 <<
" isLFVTau2BodyDecay = " << isLFVTau2BodyDecay);
495 bool isPi0GG =
false;
496 bool isEtaGG =
false;
497 bool isEtpGG =
false;
498 bool isPi0GEE =
false;
499 bool isEtaPPG =
false;
500 bool isOmPizG =
false;
501 if (mothid == 111 || mothid == 221 || mothid == 331 || mothid == 223) {
502 const vector<MCParticle*> daughters = mother->
getDaughters();
503 int numberofTotalDaughters = 0;
504 int numberOfPhotonDaughters = 0;
505 int numberOfElectronDaughters = 0;
506 int numberOfPionDaughters = 0;
507 int numberOfPizDaughters = 0;
510 numberofTotalDaughters ++;
511 if (abs(d->getPDG()) ==
Const::photon.getPDGCode()) numberOfPhotonDaughters++;
512 if (abs(d->getPDG()) ==
Const::electron.getPDGCode()) numberOfElectronDaughters++;
513 if (abs(d->getPDG()) ==
Const::pion.getPDGCode()) numberOfPionDaughters++;
514 if (abs(d->getPDG()) ==
Const::pi0.getPDGCode()) numberOfPizDaughters++;
516 if (numberofTotalDaughters == 2 && numberOfPhotonDaughters == 2) {
517 if (mothid ==
Const::pi0.getPDGCode()) isPi0GG =
true;
518 if (mothid == 221) isEtaGG =
true;
519 if (mothid == 331) isEtpGG =
true;
521 if (numberofTotalDaughters >= 3 && numberOfPhotonDaughters >= 1 && numberOfElectronDaughters == 2 && numberOfPionDaughters == 0
522 && numberOfPizDaughters == 0) {
523 if (mothid ==
Const::pi0.getPDGCode()) isPi0GEE =
true;
525 if (numberofTotalDaughters >= 3 && numberOfPhotonDaughters >= 1 && numberOfPizDaughters == 0 && numberOfPionDaughters == 2) {
528 if (chg < 0 && isEtaPPGFromTauMinusFirst) {
529 isEtaPPGFromTauMinusFirst =
false;
532 if (chg > 0 && isEtaPPGFromTauPlusFirst) {
533 isEtaPPGFromTauPlusFirst =
false;
538 if (numberofTotalDaughters >= 2 && numberOfPhotonDaughters >= 1 && numberOfPizDaughters == 1 && numberOfPionDaughters == 0) {
541 if (chg < 0 && isOmegaPizGamFromTauMinusFirst) {
542 isOmegaPizGamFromTauMinusFirst =
false;
545 if (chg > 0 && isOmegaPizGamFromTauPlusFirst) {
546 isOmegaPizGamFromTauPlusFirst =
false;
553 B2DEBUG(19,
"isRadiationfromFinalStateChargedParticle = " << isRadiationfromFinalStateChargedParticle);
554 B2DEBUG(19,
"isRadiationFromChargedRho = " << isRadiationFromChargedRho);
555 B2DEBUG(19,
"isRadiationFromChargedA1 = " << isRadiationFromChargedA1);
556 B2DEBUG(19,
"isRadiationfromIntermediateWBoson = " << isRadiationfromIntermediateWBoson <<
" isPiPizGam = " << isPiPizGam);
557 B2DEBUG(19,
"isRadiationfromTau = " << isRadiationfromTau <<
" isLFVTau2BodyDecay = " << isLFVTau2BodyDecay);
558 B2DEBUG(19,
"isPi0GG = " << isPi0GG <<
" isEtaGG = " << isEtaGG <<
" isEtpGG = " << isEtpGG <<
" isPi0GEE = " << isPi0GEE);
559 B2DEBUG(19,
"isEtaPPG = " << isEtaPPG <<
" isOmPizG = " << isOmPizG);
565 if (isRadiationfromFinalStateChargedParticle) {
566 }
else if (isRadiationFromChargedRho) {
567 accept_photon =
true;
568 }
else if (isRadiationFromChargedA1) {
569 accept_photon =
true;
570 }
else if (isRadiationfromIntermediateWBoson) {
572 accept_photon =
true;
574 }
else if (isRadiationfromTau) {
575 if (isLFVTau2BodyDecay) {
576 accept_photon =
true;
578 }
else if (isPi0GG) {
579 }
else if (isEtaGG) {
580 }
else if (isEtpGG) {
581 }
else if (isPi0GEE) {
582 }
else if (isEtaPPG) {
583 accept_photon =
true;
584 }
else if (isOmPizG) {
585 accept_photon =
true;
593 else if (pdgid ==
Const::pi0.getPDGCode() && (isEtaPizPizPizFromTauMinusFirst || isEtaPizPizPizFromTauPlusFirst)) {
596 if (mother and (mother->
getPDG() == 221)) {
597 const vector<MCParticle*> daughters = mother->
getDaughters();
601 if (d->getPDG() ==
Const::pi0.getPDGCode()) nPizSisters++;
603 B2DEBUG(19,
"nPizSisters = " << nPizSisters);
604 if (nPizSisters == 3) {
606 if (chg < 0 && isEtaPizPizPizFromTauMinusFirst) {
607 isEtaPizPizPizFromTauMinusFirst =
false;
610 if (chg > 0 && isEtaPizPizPizFromTauPlusFirst) {
611 isEtaPizPizPizFromTauPlusFirst =
false;
631 else if (pdgid == -
Const::pion.getPDGCode() && (isOmegaPimPipFromTauMinusFirst || isOmegaPimPipFromTauPlusFirst)) {
634 if (mother and (mother->
getPDG() == 223)) {
635 const vector<MCParticle*> daughters = mother->
getDaughters();
636 int nOmegaDaughters = 0;
643 if (d->getPDG() == -
Const::pion.getPDGCode()) nPimSisters++;
644 if (d->getPDG() ==
Const::pion.getPDGCode()) nPipSisters++;
645 if (d->getPDG() ==
Const::pi0.getPDGCode()) nPizSisters++;
648 "nOmegaDaughters = " << nOmegaDaughters <<
" "
649 "nPimSisters = " << nPimSisters <<
" "
650 "nPipSisters = " << nPipSisters <<
" "
651 "nPizSisters = " << nPizSisters);
652 if (nOmegaDaughters >= 2 && nPimSisters == 1 && nPipSisters == 1 && nPizSisters == 0) {
654 if (chg < 0 && isOmegaPimPipFromTauMinusFirst) {
655 isOmegaPimPipFromTauMinusFirst =
false;
658 if (chg > 0 && isOmegaPimPipFromTauPlusFirst) {
659 isOmegaPimPipFromTauPlusFirst =
false;
666 "isEtaPizPizPizFromTauMinusFirst = " << isEtaPizPizPizFromTauMinusFirst <<
" "
668 "isEtaPizPizPizFromTauPlusFirst = " << isEtaPizPizPizFromTauPlusFirst <<
" "
672 "isOmegaPimPipFromTauMinusFirst = " << isOmegaPimPipFromTauMinusFirst <<
" "
674 "isOmegaPimPipFromTauPlusFirst = " << isOmegaPimPipFromTauPlusFirst <<
" "
679 map<int, std::vector<int>>::iterator ite ;
680 for (ite = map_vec.begin(); ite != map_vec.end(); ++ite) {
681 if (pdgid == ite->first) {
684 B2DEBUG(19,
"Photon accepted");
685 ite-> second.push_back(i);
688 ite-> second.push_back(i);
699 map<int, std::vector<int>>::iterator itr ;
700 for (itr = map_vec.begin(); itr != map_vec.end(); ++itr) {
701 for (
unsigned int i = 0; i < itr-> second.size(); i++) {
702 int ii = itr-> second[i];
714 for (
unsigned iorder = 0; iorder <
OrderedList.size(); ++iorder) {
719 int pdg = p->getPDG();
721 if (pdg != ii)
continue;
726 pdgname = databasePDG->GetParticle(pdg)->GetName();
736 for (
unsigned iorder = 0; iorder <
OrderedList.size(); ++iorder) {
741 int pdg = p->getPDG();
743 if (pdg != ii)
continue;
748 pdgname = databasePDG->GetParticle(pdg)->GetName();
763 for (itr = mode_decay.begin(); itr != mode_decay.end(); ++itr) {
764 if (state == itr-> first) {
765 int mode = itr-> second;
792 if (mother ==
nullptr) {
804 }
else if (abs(p->getPDG()) ==
Const::muon.getPDGCode() && mother->
getPDG() == 221) {
806 }
else if (abs(p->getPDG()) ==
Const::muon.getPDGCode() && mother->
getPDG() == 113) {
808 }
else if (abs(p->getPDG()) ==
Const::muon.getPDGCode() && mother->
getPDG() == 333) {
810 }
else if (mother->
getPDG() == 15) {
812 }
else if (mother->
getPDG() == -15) {
826 for (
int i = 0; i <
MCParticles.getEntries(); i++) {
829 if (p.getStatus() == 1 && p.getPDG() == 15) {
833 if (p.getStatus() == 1 && p.getPDG() == -15) {
847 const vector<MCParticle*> daughters = p.getDaughters();
848 if (daughters.empty())
return ret;
855 if (isChargedFinalState) ret++;
865 ROOT::Math::PxPyPzEVector p4_tau = p.get4Vector();
866 B2DEBUG(19,
"p4_tau: " << p4_tau <<
" " << p4_tau.P());
867 ROOT::Math::Boost boost_to_mother_rest_frame(p4_tau.BoostToCM());
868 for (
auto i : vec_radgam) {
869 ROOT::Math::PxPyPzEVector p4_gamma =
MCParticles[i]->get4Vector();
870 B2DEBUG(19,
"p4_gamma: " << p4_gamma <<
" " << p4_gamma.P());
871 ROOT::Math::PxPyPzEVector p4_gamma_rest = boost_to_mother_rest_frame * p4_gamma;
872 B2DEBUG(19,
"p4_gamma_rest: " << p4_gamma_rest <<
" " << p4_gamma_rest.P());
873 double energy_rest = p4_gamma_rest.E();
874 if (energy_rest > egstar) egstar = energy_rest;
876 B2DEBUG(19,
"egstar: " << egstar);
static const ParticleType pi0
neutral pion particle
static const ChargedStable muon
muon particle
static const ChargedStable pion
charged pion particle
static const ParticleType Klong
K^0_L particle.
static const ParticleType photon
photon particle
static const ChargedStable electron
electron particle
Replacement for TDatabasePDG that is filled from EvtGen's evt.pdl.
static EvtGenDatabasePDG * Instance()
Instance method that loads the EvtGen table.
A Class to store the Monte Carlo particle information.
@ c_IsPHOTOSPhoton
bit 8: Particle is an radiative photon from PHOTOS
@ c_PrimaryParticle
bit 0: Particle is primary particle.
std::vector< Belle2::MCParticle * > getDaughters() const
Get vector of all daughter particles, empty vector if none.
int getPDG() const
Return PDG code of particle.
void setDescription(const std::string &description)
Sets the description of the module.
std::vector< int > vec_numu
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_em
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_omega
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_k0s
Variable name of the vector where particles identified in the event are stored.
Int_t m_mprong
Prong of the decay channel of tau-.
std::vector< int > vec_b1p
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_f1
Variable name of the vector where particles identified in the event are stored.
int numOfTauPlus
Number of tau+ in the event.
Int_t m_pprong
Prong of the decay channel of tau+.
StoreObjPtr< EventMetaData > m_event_metadata
event number
void IdentifyTauPair()
Identifies if the event is a generated tau pair.
std::vector< int > vec_apro
Variable name of the vector where particles identified in the event are stored.
double getEgstar(const std::vector< int > &vec_radgam, const MCParticle &mc)
Energy of the radiative photon in tau rest frame.
std::vector< int > vec_pi0
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_alpha
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_radgam_taum
Variable name of the vector where particles identified in the event are stored.
StoreArray< MCParticle > MCParticles
StoreArray of MCParticles.
std::vector< int > vec_anue
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_nue
Variable name of the vector where particles identified in the event are stored.
virtual void initialize() override
Initializes the module.
std::vector< int > vec_a0p_980
Variable name of the vector where particles identified in the event are stored.
std::string m_file_plus
Alternative mapping for tau+.
static constexpr std::array< int, 46 > OrderedList
PDG codes of ORDERED particles.
std::vector< int > vec_a0m_1450
Variable name of the vector where particles identified in the event are stored.
virtual void event() override
Method is called for each event.
bool m_isEtaPizPizPizFromTauMinus
Flag for eta->pi0pi0pi0 decays from tau-.
static constexpr std::array< int, 3 > Neutrinos
PDG codes of neutrinos in final state particles in generation: {nu_e, nu_mu, mu_tau}...
double m_pegstar
Energy of radiative photon from tau+.
std::vector< int > vec_lmb_br
Variable name of the vector where particles identified in the event are stored.
Int_t m_mmode
ID of the decay channel of tau-.
void AnalyzeTauPairEvent()
Analyze a generated tau pair event.
bool m_isOmegaPimPipFromTauMinus
Flag for omega->pi-pi+ decays from tau-.
int getRecursiveMotherCharge(const MCParticle *mc)
Identifies particles coming from tau decays.
std::vector< int > vec_phi
Variable name of the vector where particles identified in the event are stored.
int m_taum_no
number of tau- unclassified events
std::string m_tauminusdecaymode
Variable name for the decay mode of the tau-.
std::vector< int > vec_a1p
Variable name of the vector where particles identified in the event are stored.
bool m_isOmegaPimPipFromTauPlus
Flag for omega->pi-pi+ decays from tau+.
bool tauPair
Boolean variable used to identify tau event.
int TauolaBelle2DecayMode(const std::string &s, int chg)
Classifies the decays of the event and assigns a decay mode.
int idOfTauMinus
Index of the generated tau-.
std::vector< int > vec_lambda
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_anut
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_f0
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_mum
Variable name of the vector where particles identified in the event are stored.
Int_t m_pmode
ID of the decay channel of tau+.
int numOfTauMinus
Number of tau- in the event.
std::vector< int > vec_nut
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_kstar0_br
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_radgam_taup
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_a00_1450
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_kstar0
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_a0m_980
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_a00_980
Variable name of the vector where particles identified in the event are stored.
int getProngOfDecay(const MCParticle &mc)
Identifies the number of charged final state particles in the decay.
std::vector< int > vec_a1m
Variable name of the vector where particles identified in the event are stored.
int idOfTauPlus
Index of the generated tau+.
StoreObjPtr< TauPairDecay > m_tauDecay
pointer to tau pair decay objects
std::vector< int > vec_rhop
Variable name of the vector where particles identified in the event are stored.
std::string m_printmode
Parameter passed by the user to indicate the information to be printed.
std::vector< int > vec_pro
Variable name of the vector where particles identified in the event are stored.
bool m_isEtaPizPizPizFromTauPlus
Flag for eta->pi0pi0pi0 decays from tau+.
std::string m_tauplusdecaymode
Variable name for the decay mode of the tau+.
std::vector< int > vec_pim
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_mup
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_pip
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_kstarp
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_kstarm
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_etapr
Variable name of the vector where particles identified in the event are stored.
double m_megstar
Energy of radiative photon from tau-.
std::vector< int > vec_ep
Variable name of the vector where particles identified in the event are stored.
TauDecayModeModule()
Constructor: Sets the description, the properties and the parameters of the module.
std::vector< int > vec_k0l
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_eta
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_km
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_anumu
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_dau_tauplus
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_rhom
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_gam
Variable name of the vector where particles identified in the event are stored.
std::string m_file_minus
Alternative mapping for tau-.
int m_taup_no
number of tau+ unclassified events
std::vector< int > vec_dau_tauminus
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_kp
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_b1m
Variable name of the vector where particles identified in the event are stored.
std::vector< int > vec_rho0
Variable name of the vector where particles identified in the event are stored.
static constexpr std::array< int, 5 > finalStatePDGs
PDG codes accepted as charged final state particles in generation: {e, mu, pi, K, p}...
std::map< std::string, int > mode_decay_minus
Mapping for tau- decays.
std::map< std::string, int > mode_decay_plus
Mapping for tau+ decays.
std::vector< int > vec_a0p_1450
Variable name of the vector where particles identified in the event are stored.
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.