From 5c864b9051af0366f5aca07fe8feef9da4194c72 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Rafa=C5=82=20Lalik?= Date: Thu, 8 Jun 2023 13:51:50 +0200 Subject: [PATCH 1/4] Strangeness plugin with new resonances and decay branches --- .../strangeness/include/PStrangenessPlugin.h | 11 + plugins/strangeness/src/PStrangenessPlugin.cc | 266 ++++++++++++++---- 2 files changed, 218 insertions(+), 59 deletions(-) diff --git a/plugins/strangeness/include/PStrangenessPlugin.h b/plugins/strangeness/include/PStrangenessPlugin.h index d20b8bd..4a2963e 100644 --- a/plugins/strangeness/include/PStrangenessPlugin.h +++ b/plugins/strangeness/include/PStrangenessPlugin.h @@ -1,6 +1,8 @@ // Author: I. Froehlich // Written: 12.10.2008 // Revised: +// Updates: Rafał Lalik, 05.06.2023 - New Y* decays channels +// - Software switches to modify decays // @@ -27,10 +29,19 @@ class PStrangenessPlugin : public PDistributionCollection { Bool_t ExecCommand(const char *command, Double_t value); Bool_t Activate(void); + static void EnableHadronDecays(Bool_t enabled) { hyperon_hadron_decays = enabled; } + static void EnablePhotonDecays(Bool_t enabled) { hyperon_photon_decays = enabled; } + static void EnableDalitzDecays(Bool_t enabled) { hyperon_dalitz_decays = enabled; } + static void EnableExperimentalDecays(Bool_t enabled) { hyperon_experimental_decays = enabled; } private: Int_t is_initialized; + static Bool_t hyperon_hadron_decays; + static Bool_t hyperon_photon_decays; + static Bool_t hyperon_dalitz_decays; + static Bool_t hyperon_experimental_decays; + ClassDef(PStrangenessPlugin, 0) //Plugin to add new particles with strangeness and their decays }; diff --git a/plugins/strangeness/src/PStrangenessPlugin.cc b/plugins/strangeness/src/PStrangenessPlugin.cc index 96ba5de..57ff9fb 100644 --- a/plugins/strangeness/src/PStrangenessPlugin.cc +++ b/plugins/strangeness/src/PStrangenessPlugin.cc @@ -14,9 +14,14 @@ #include "PDistributionManager.h" #include "PResonanceDalitz.h" -PStrangenessPlugin::PStrangenessPlugin(const Char_t *id, const Char_t *de): - PDistributionCollection(id, de) { - is_initialized = 0; + +Bool_t PStrangenessPlugin::hyperon_hadron_decays{kTRUE}; +Bool_t PStrangenessPlugin::hyperon_photon_decays{kTRUE}; +Bool_t PStrangenessPlugin::hyperon_dalitz_decays{kTRUE}; +Bool_t PStrangenessPlugin::hyperon_experimental_decays{kFALSE}; + +PStrangenessPlugin::PStrangenessPlugin(const Char_t *id, const Char_t *de): PDistributionCollection(id, de), + is_initialized(0) { // RequiresPlugin("elementary"); } @@ -31,80 +36,223 @@ PStrangenessPlugin::~PStrangenessPlugin() { Bool_t PStrangenessPlugin::ExecCommand(const char *command, Double_t) { if (strcmp (command,"init") == 0) { - if (!is_initialized) { - is_initialized = 1; + if (!is_initialized) { + is_initialized = 1; + + fprintf(stderr, "Strangeness plugin - flags status:\n" + "- Hadron: %d\n" + "- Photon: %d\n" + "- Dalitz: %d\n" + "- Experimental: %d\n", + hyperon_hadron_decays, hyperon_photon_decays, hyperon_dalitz_decays, hyperon_experimental_decays); + + Double_t *system_alloc_verbosity = makeStaticData()->GetBatchValue("_system_alloc_verbosity"); + Double_t old = (*system_alloc_verbosity); + (*system_alloc_verbosity) = 0; + + if (makeStaticData()->AddParticle(70, "Lambda1405", 1.406) > 0) { + makeStaticData()->AddAlias("Lambda1405", "Lambda(1405)"); + + makeStaticData()->SetParticleTotalWidth("Lambda1405", 0.05); + makeStaticData()->SetParticleBaryon("Lambda1405", 1); + makeStaticData()->SetParticleSpin("Lambda1405", 1); + makeStaticData()->SetParticleParity("Lambda1405", -1); + + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Lambda(1405) --> Sigma+ + pi-", "Lambda1405", "Sigma+,pi-", 1. / 3); + makeStaticData()->AddDecay("Lambda(1405) --> Sigma- + pi+", "Lambda1405", "Sigma-,pi+", 1. / 3); + makeStaticData()->AddDecay("Lambda(1405) --> Sigma0 + pi0", "Lambda1405", "Sigma0,pi0", 1. / 3); + } + + if (hyperon_photon_decays) { + makeStaticData()->AddDecay("Lambda(1405) --> Lambda + gamma", "Lambda1405", "Lambda,g", 0.0085); + } + + if (hyperon_dalitz_decays) { + makeStaticData()->AddDecay("Lambda(1405) --> Lambda + dilepton", "Lambda1405", "Lambda,dilepton", 0.0085 / 137.); + + PResonanceDalitz * dalitz = new PResonanceDalitz("Lambda1405_dalitz@Lambda1405_to_Lambda_dilepton", "dgdm from Zetenyi/Wolf", -1); + dalitz->setGm(0.719); + makeDistributionManager()->Add(dalitz); + } + } else { + Error("ExecCommand", "PIDs blocked, plugin disabled"); + } + + + if (makeStaticData()->AddParticle(-1, "Sigma1385+", 1.3828) > 0) { + makeStaticData()->AddAlias("Sigma1385+", "Sigma(1385)+"); + + makeStaticData()->SetParticleTotalWidth("Sigma1385+", 0.0358); + makeStaticData()->SetParticleBaryon("Sigma1385+", 1); + makeStaticData()->SetParticleLMass("Sigma1385+", 1.2); + + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Sigma(1385)+ --> Lambda + pi+", "Sigma1385+", "Lambda,pi+", 0.8815); + makeStaticData()->AddDecay("Sigma(1385)+ --> Sigma+ + pi0", "Sigma1385+", "Sigma+,pi0", 0.05925); + makeStaticData()->AddDecay("Sigma(1385)+ --> Sigma0 + pi+", "Sigma1385+", "Sigma0,pi+", 0.05925); + } + } else { + Error("ExecCommand", "PIDs blocked, plugin disabled"); + } + + + if (makeStaticData()->AddParticle(-1, "Sigma1385-", 1.3872) > 0) { + makeStaticData()->AddAlias("Sigma1385-", "Sigma(1385)-"); + + makeStaticData()->SetParticleTotalWidth("Sigma1385-", 0.0394); + makeStaticData()->SetParticleBaryon("Sigma1385-", 1); + makeStaticData()->SetParticleLMass("Sigma1385-", 1.2); + + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Sigma(1385)- --> Lambda + pi-", "Sigma1385-", "Lambda,pi-", 0.8815 ); + makeStaticData()->AddDecay("Sigma(1385)- --> Sigma- + pi0", "Sigma1385-", "Sigma-,pi0", 0.05925); + makeStaticData()->AddDecay("Sigma(1385)- --> Sigma0 + pi-", "Sigma1385-", "Sigma0,pi-", 0.05925); + } + } else { + Error("ExecCommand", "PIDs blocked, plugin disabled"); + } + + + if (makeStaticData()->AddParticle(-1, "Sigma13850", 1.3837) > 0) { + makeStaticData()->AddAlias("Sigma13850", "Sigma(1385)0"); + + makeStaticData()->SetParticleTotalWidth("Sigma13850", 0.0365); + makeStaticData()->SetParticleBaryon("Sigma13850", 1); + makeStaticData()->SetParticleLMass("Sigma13850", 1.2); + makeStaticData()->SetParticleSpin("Sigma13850", 3); + makeStaticData()->SetParticleParity("Sigma13850", 1); + + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + pi0", "Sigma13850", "Lambda,pi0", 0.87000); + makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + g", "Sigma13850", "Lambda,g", 0.01300 * (1. - 1./137.)); + makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma+ + pi-", "Sigma13850", "Sigma+,pi-", 0.117 / 3); + makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma- + pi+", "Sigma13850", "Sigma-,pi+", 0.117 / 3); + makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma0 + pi0", "Sigma13850", "Sigma0,pi0", 0.117 / 3); + } + + if (hyperon_dalitz_decays) { + makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + dilepton", "Sigma13850", "Lambda,dilepton", 0.013 / 137.); + PResonanceDalitz * dalitz = new PResonanceDalitz("Sigma13850_dalitz@Sigma13850_to_Lambda_dilepton", "dgdm from Zetenyi/Wolf", -1); + dalitz->setGm(0.719); + makeDistributionManager()->Add(dalitz); + } + } else { + Error("ExecCommand", "PIDs blocked, plugin disabled"); + } + + + if (makeStaticData()->AddParticle(-1,"Lambda1520", 1.5195) > 0) { + makeStaticData()->AddAlias("Lambda1520","Lambda(1520)"); + + makeStaticData()->SetParticleTotalWidth("Lambda1520", 0.0156); + makeStaticData()->SetParticleBaryon("Lambda1520", 1); + makeStaticData()->SetParticleSpin("Lambda1520", 3); + makeStaticData()->SetParticleParity("Lambda1520", 1); + + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi0 + pi0", "Lambda1520", "Lambda,pi0,pi0", 0.2 / 3); // - + makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi+ + pi-", "Lambda1520", "Lambda,pi+,pi-", 0.1 / 3); // old BR .014638 + } + + if (hyperon_photon_decays) { + makeStaticData()->AddDecay("Lambda(1520) --> Lambda + g", "Lambda1520", "Lambda,g", 0.0085); + } + + if (hyperon_experimental_decays) { + makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi0", "Lambda1520", "Sigma0,pi0", 0.139096); // - + makeStaticData()->AddDecay("Lambda(1520) --> Sigma+ + pi-", "Lambda1520", "Sigma+,pi-", 0.139096); // - + makeStaticData()->AddDecay("Lambda(1520) --> Sigma- + pi+", "Lambda1520", "Sigma-,pi+", 0.139096); // - + makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi0 + pi0", "Lambda1520", "Sigma0,pi0,pi0", 0.009/5.0); // - + makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi+ + pi-", "Lambda1520", "Sigma0,pi+,pi-", 0.009/5.0); // - + } + + if (hyperon_dalitz_decays) { + makeStaticData()->AddDecay("Lambda(1520) --> Lambda + dilepton", "Lambda1520", "Lambda,dilepton", .007948/137.); + + PResonanceDalitz * dalitz = new PResonanceDalitz("Lambda1520_dalitz@Lambda1520_to_Lambda_dilepton", "dgdm from Zetenyi/Wolf", -1); + dalitz->setGm(0.719); + makeDistributionManager()->Add(dalitz); + } + } else { + Error("ExecCommand", "PIDs blocked, plugin disabled"); + } + - Double_t *system_alloc_verbosity = makeStaticData()->GetBatchValue("_system_alloc_verbosity"); - Double_t old = (*system_alloc_verbosity); - (*system_alloc_verbosity) = 0; + if (makeStaticData()->AddParticle(75,"Xi1530-", 1.5352) > 0) { + makeStaticData()->AddAlias("Xi1530-","Xi(1530)-"); - if (makeStaticData()->AddParticle(70, "Lambda1405", 1.406) > 0) { - makeStaticData()->AddAlias("Lambda1405", "Lambda(1405)"); + makeStaticData()->SetParticleTotalWidth("Xi1530-", 0.091); + makeStaticData()->SetParticleBaryon("Xi1530-", 1); + makeStaticData()->SetParticleSpin("Xi1530-", 2); + makeStaticData()->SetParticleParity("Xi1530-", 1); - makeStaticData()->SetParticleTotalWidth("Lambda1405", 0.05); - makeStaticData()->SetParticleBaryon("Lambda1405", 1); + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Xi(1530)- --> Xi0 + pi-", "Xi1530-", "Xi0,pi-", 1. / 2.0); + makeStaticData()->AddDecay("Xi(1530)- --> Xi- + pi0", "Xi1530-", "Xi-,pi0", 1. / 2.0); + } + if (hyperon_photon_decays) { + makeStaticData()->AddDecay("Xi(1530)- --> Xi- + g", "Xi1530-", "Xi-,g", .037); + } - makeStaticData()->AddDecay("Lambda(1405)0 --> Sigma+ + pi-", "Lambda1405", - "Sigma+,pi-", .3333 ); - makeStaticData()->AddDecay("Lambda(1405)0 --> Sigma- + pi+", "Lambda1405", - "Sigma-,pi+", .3333 ); - makeStaticData()->AddDecay("Lambda(1405)0 --> Sigma0 + pi0", "Lambda1405", - "Sigma0,pi0", .3333 ); + if (hyperon_dalitz_decays) { + makeStaticData()->AddDecay("Xi(1530)- --> Xi- + dilepton", "Xi1530-", "Xi-,dilepton", .037/137.); + PResonanceDalitz * dalitz = new PResonanceDalitz("Xi1530-_dalitz@Xi1530-_to_Xi-_dilepton", "dgdm from Zetenyi/Wolf", -1); + dalitz->setGm(0.719); + makeDistributionManager()->Add(dalitz); + } + } else { + Error("ExecCommand", "PIDs blocked, plugin disabled"); + } - makeStaticData()->AddParticle(71, "Sigma1385+", 1.3828); - makeStaticData()->AddAlias("Sigma1385+", "Sigma(1385)+"); - makeStaticData()->SetParticleTotalWidth("Sigma1385+", 0.0358); - makeStaticData()->SetParticleBaryon("Sigma1385+", 1); - makeStaticData()->SetParticleLMass("Sigma1385+", 1.2); - makeStaticData()->AddDecay("Sigma(1385)+ --> Lambda + pi+", "Sigma1385+", - "Lambda,pi+", .8815 ); - makeStaticData()->AddDecay("Sigma(1385)+ --> Sigma+ + pi0", "Sigma1385+", - "Sigma+, pi0", .05925); - makeStaticData()->AddDecay("Sigma(1385)+ --> Sigma0 + pi+", "Sigma1385+", - "Sigma0, pi+", .05925); + if (makeStaticData()->AddParticle(76,"Xi15300", 1.53178) > 0) { + makeStaticData()->AddAlias("Xi15300","Xi(1530)0"); + makeStaticData()->SetParticleTotalWidth("Xi15300", 0.099); + makeStaticData()->SetParticleBaryon("Xi15300", 1); + makeStaticData()->SetParticleSpin("Xi15300", 2); + makeStaticData()->SetParticleParity("Xi15300", 1); + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Xi(1530)0 --> Xi0 + pi0", "Xi15300", "Xi0,pi0", 1. / 2.0); + makeStaticData()->AddDecay("Xi(1530)0 --> Xi- + pi+", "Xi15300", "Xi-,pi+", 1. / 2.0); + } - makeStaticData()->AddParticle(72, "Sigma1385-", 1.3872); - makeStaticData()->AddAlias("Sigma1385-", "Sigma(1385)-"); - makeStaticData()->SetParticleTotalWidth("Sigma1385-", 0.0394); - makeStaticData()->SetParticleBaryon("Sigma1385-", 1); - makeStaticData()->SetParticleLMass("Sigma1385-", 1.2); + if (hyperon_photon_decays) { + makeStaticData()->AddDecay("Xi(1530)0 --> Xi0 + g", "Xi15300", "Xi0,g", .037); + } - makeStaticData()->AddDecay("Sigma(1385)- --> Lambda + pi-", "Sigma1385-", - "Lambda,pi-", .8815 ); - makeStaticData()->AddDecay("Sigma(1385)- --> Sigma- + pi0", "Sigma1385-", - "Sigma-, pi0", .05925); - makeStaticData()->AddDecay("Sigma(1385)- --> Sigma0 + pi-", "Sigma1385-", - "Sigma0, pi-", .05925); + if (hyperon_dalitz_decays) { + makeStaticData()->AddDecay("Xi(1530)0 --> Xi0 + dilepton", "Xi15300", "Xi0,dilepton", .037/137.); + PResonanceDalitz * dalitz = new PResonanceDalitz("Xi15300_dalitz@Xi15300_to_Xi0_dilepton", "dgdm from Zetenyi/Wolf", -1); + dalitz->setGm(0.719); + makeDistributionManager()->Add(dalitz); + } + } else { + Error("ExecCommand", "PIDs blocked, plugin disabled"); + } - makeStaticData()->AddParticle(73, "Sigma13850", 1.3837); - makeStaticData()->AddAlias("Sigma13850", "Sigma(1385)0"); - makeStaticData()->SetParticleTotalWidth("Sigma13850", 0.0365); - makeStaticData()->SetParticleBaryon("Sigma13850", 1); - makeStaticData()->SetParticleLMass("Sigma13850", 1.2); + if (makeStaticData()->AddParticle(79, "Delta2050++", 2.05) > 0) { + makeStaticData()->AddAlias("Delta2050++", "Delta(2050)++"); - makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + pi0", "Sigma13850", - "Lambda,pi0", .87000 ); - makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + g", "Sigma13850", - "Lambda,g", .01300 ); - makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma+ + pi-", "Sigma13850", - "Sigma+, pi-", .05850); - makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma- + pi+", "Sigma13850", - "Sigma-, pi+", .05850); + makeStaticData()->SetParticleTotalWidth("Delta2050++", 0.25); + makeStaticData()->SetParticleBaryon("Delta2050++", 1); - } else { - Error("ExecCommand", "PIDs blocked, plugin disabled"); - } + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Delta(2050)++ --> Sigma(1385)+ + K+", "Delta2050++", "Sigma1385+,K+", 1.0); + } + } else { + Error("ExecCommand", "PIDs blocked, plugin disabled"); + } - (*system_alloc_verbosity) = old; - } - return kTRUE; + (*system_alloc_verbosity) = old; + } + return kTRUE; } return kFALSE; From e372ad0f1d4d15627170c0eb69a5026baab2716b Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Rafa=C5=82=20Lalik?= Date: Thu, 8 Jun 2023 18:24:19 +0200 Subject: [PATCH 2/4] Add NK decay channel for L1520 with BR=0.45 as of PDG --- plugins/strangeness/src/PStrangenessPlugin.cc | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/plugins/strangeness/src/PStrangenessPlugin.cc b/plugins/strangeness/src/PStrangenessPlugin.cc index 57ff9fb..c48f35f 100644 --- a/plugins/strangeness/src/PStrangenessPlugin.cc +++ b/plugins/strangeness/src/PStrangenessPlugin.cc @@ -151,8 +151,10 @@ Bool_t PStrangenessPlugin::ExecCommand(const char *command, Double_t) { makeStaticData()->SetParticleParity("Lambda1520", 1); if (hyperon_hadron_decays) { - makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi0 + pi0", "Lambda1520", "Lambda,pi0,pi0", 0.2 / 3); // - - makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi+ + pi-", "Lambda1520", "Lambda,pi+,pi-", 0.1 / 3); // old BR .014638 + makeStaticData()->AddDecay("Lambda(1520) --> n + K0S", "Lambda1520", "n,K0S", 0.45 / 2); // N + /K from PDG + makeStaticData()->AddDecay("Lambda(1520) --> n + K0L", "Lambda1520", "n,K0L", 0.45 / 2); // N + /K from PDG + makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi0 + pi0", "Lambda1520", "Lambda,pi0,pi0", 0.2 / 3); // - + makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi+ + pi-", "Lambda1520", "Lambda,pi+,pi-", 0.1 / 3); // old BR .014638 } if (hyperon_photon_decays) { From 3799c3050b1b2d4118ef8b158c665f08cb6659cb Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Rafa=C5=82=20Lalik?= Date: Fri, 9 Jun 2023 16:30:06 +0200 Subject: [PATCH 3/4] Make strangeness flags printing as optional --- .../strangeness/include/PStrangenessPlugin.h | 2 ++ plugins/strangeness/src/PStrangenessPlugin.cc | 17 ++++++++++------- 2 files changed, 12 insertions(+), 7 deletions(-) diff --git a/plugins/strangeness/include/PStrangenessPlugin.h b/plugins/strangeness/include/PStrangenessPlugin.h index 4a2963e..268206c 100644 --- a/plugins/strangeness/include/PStrangenessPlugin.h +++ b/plugins/strangeness/include/PStrangenessPlugin.h @@ -33,6 +33,8 @@ class PStrangenessPlugin : public PDistributionCollection { static void EnablePhotonDecays(Bool_t enabled) { hyperon_photon_decays = enabled; } static void EnableDalitzDecays(Bool_t enabled) { hyperon_dalitz_decays = enabled; } static void EnableExperimentalDecays(Bool_t enabled) { hyperon_experimental_decays = enabled; } + static void ShowFlags(); + private: Int_t is_initialized; diff --git a/plugins/strangeness/src/PStrangenessPlugin.cc b/plugins/strangeness/src/PStrangenessPlugin.cc index c48f35f..bece206 100644 --- a/plugins/strangeness/src/PStrangenessPlugin.cc +++ b/plugins/strangeness/src/PStrangenessPlugin.cc @@ -20,6 +20,16 @@ Bool_t PStrangenessPlugin::hyperon_photon_decays{kTRUE}; Bool_t PStrangenessPlugin::hyperon_dalitz_decays{kTRUE}; Bool_t PStrangenessPlugin::hyperon_experimental_decays{kFALSE}; +void PStrangenessPlugin::ShowFlags() +{ + fprintf(stderr, "Strangeness plugin - flags status:\n" + "- Hadron: %d\n" + "- Photon: %d\n" + "- Dalitz: %d\n" + "- Experimental: %d\n", + hyperon_hadron_decays, hyperon_photon_decays, hyperon_dalitz_decays, hyperon_experimental_decays); +} + PStrangenessPlugin::PStrangenessPlugin(const Char_t *id, const Char_t *de): PDistributionCollection(id, de), is_initialized(0) { // RequiresPlugin("elementary"); @@ -39,13 +49,6 @@ Bool_t PStrangenessPlugin::ExecCommand(const char *command, Double_t) { if (!is_initialized) { is_initialized = 1; - fprintf(stderr, "Strangeness plugin - flags status:\n" - "- Hadron: %d\n" - "- Photon: %d\n" - "- Dalitz: %d\n" - "- Experimental: %d\n", - hyperon_hadron_decays, hyperon_photon_decays, hyperon_dalitz_decays, hyperon_experimental_decays); - Double_t *system_alloc_verbosity = makeStaticData()->GetBatchValue("_system_alloc_verbosity"); Double_t old = (*system_alloc_verbosity); (*system_alloc_verbosity) = 0; From fa8b4b041c0ccbcd082e9e1a1da33c49ad4526c3 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Rafa=C5=82=20Lalik?= Date: Mon, 12 Jun 2023 15:42:48 +0200 Subject: [PATCH 4/4] Fix hyperons order and move experimental to main production group --- plugins/strangeness/src/PStrangenessPlugin.cc | 84 +++++++++---------- 1 file changed, 42 insertions(+), 42 deletions(-) diff --git a/plugins/strangeness/src/PStrangenessPlugin.cc b/plugins/strangeness/src/PStrangenessPlugin.cc index bece206..5a31fd9 100644 --- a/plugins/strangeness/src/PStrangenessPlugin.cc +++ b/plugins/strangeness/src/PStrangenessPlugin.cc @@ -53,28 +53,27 @@ Bool_t PStrangenessPlugin::ExecCommand(const char *command, Double_t) { Double_t old = (*system_alloc_verbosity); (*system_alloc_verbosity) = 0; - if (makeStaticData()->AddParticle(70, "Lambda1405", 1.406) > 0) { - makeStaticData()->AddAlias("Lambda1405", "Lambda(1405)"); - makeStaticData()->SetParticleTotalWidth("Lambda1405", 0.05); - makeStaticData()->SetParticleBaryon("Lambda1405", 1); - makeStaticData()->SetParticleSpin("Lambda1405", 1); - makeStaticData()->SetParticleParity("Lambda1405", -1); + if (makeStaticData()->AddParticle(70, "Sigma13850", 1.3837) > 0) { + makeStaticData()->AddAlias("Sigma13850", "Sigma(1385)0"); - if (hyperon_hadron_decays) { - makeStaticData()->AddDecay("Lambda(1405) --> Sigma+ + pi-", "Lambda1405", "Sigma+,pi-", 1. / 3); - makeStaticData()->AddDecay("Lambda(1405) --> Sigma- + pi+", "Lambda1405", "Sigma-,pi+", 1. / 3); - makeStaticData()->AddDecay("Lambda(1405) --> Sigma0 + pi0", "Lambda1405", "Sigma0,pi0", 1. / 3); - } + makeStaticData()->SetParticleTotalWidth("Sigma13850", 0.0365); + makeStaticData()->SetParticleBaryon("Sigma13850", 1); + makeStaticData()->SetParticleLMass("Sigma13850", 1.2); + makeStaticData()->SetParticleSpin("Sigma13850", 3); + makeStaticData()->SetParticleParity("Sigma13850", 1); - if (hyperon_photon_decays) { - makeStaticData()->AddDecay("Lambda(1405) --> Lambda + gamma", "Lambda1405", "Lambda,g", 0.0085); + if (hyperon_hadron_decays) { + makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + pi0", "Sigma13850", "Lambda,pi0", 0.87000); + makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + g", "Sigma13850", "Lambda,g", 0.01300 * (1. - 1./137.)); + makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma+ + pi-", "Sigma13850", "Sigma+,pi-", 0.117 / 3); + makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma- + pi+", "Sigma13850", "Sigma-,pi+", 0.117 / 3); + makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma0 + pi0", "Sigma13850", "Sigma0,pi0", 0.117 / 3); } if (hyperon_dalitz_decays) { - makeStaticData()->AddDecay("Lambda(1405) --> Lambda + dilepton", "Lambda1405", "Lambda,dilepton", 0.0085 / 137.); - - PResonanceDalitz * dalitz = new PResonanceDalitz("Lambda1405_dalitz@Lambda1405_to_Lambda_dilepton", "dgdm from Zetenyi/Wolf", -1); + makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + dilepton", "Sigma13850", "Lambda,dilepton", 0.013 / 137.); + PResonanceDalitz * dalitz = new PResonanceDalitz("Sigma13850_dalitz@Sigma13850_to_Lambda_dilepton", "dgdm from Zetenyi/Wolf", -1); dalitz->setGm(0.719); makeDistributionManager()->Add(dalitz); } @@ -117,26 +116,29 @@ Bool_t PStrangenessPlugin::ExecCommand(const char *command, Double_t) { } - if (makeStaticData()->AddParticle(-1, "Sigma13850", 1.3837) > 0) { - makeStaticData()->AddAlias("Sigma13850", "Sigma(1385)0"); + if (makeStaticData()->AddParticle(-1, "Lambda1405", 1.406) > 0) { + makeStaticData()->AddAlias("Lambda1405", "Lambda(1405)"); - makeStaticData()->SetParticleTotalWidth("Sigma13850", 0.0365); - makeStaticData()->SetParticleBaryon("Sigma13850", 1); - makeStaticData()->SetParticleLMass("Sigma13850", 1.2); - makeStaticData()->SetParticleSpin("Sigma13850", 3); - makeStaticData()->SetParticleParity("Sigma13850", 1); + makeStaticData()->SetParticleTotalWidth("Lambda1405", 0.05); + makeStaticData()->SetParticleBaryon("Lambda1405", 1); + makeStaticData()->SetParticleSpin("Lambda1405", 1); + makeStaticData()->SetParticleParity("Lambda1405", -1); if (hyperon_hadron_decays) { - makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + pi0", "Sigma13850", "Lambda,pi0", 0.87000); - makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + g", "Sigma13850", "Lambda,g", 0.01300 * (1. - 1./137.)); - makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma+ + pi-", "Sigma13850", "Sigma+,pi-", 0.117 / 3); - makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma- + pi+", "Sigma13850", "Sigma-,pi+", 0.117 / 3); - makeStaticData()->AddDecay("Sigma(1385)0 --> Sigma0 + pi0", "Sigma13850", "Sigma0,pi0", 0.117 / 3); + makeStaticData()->AddDecay("Lambda(1405) --> Sigma+ + pi-", "Lambda1405", "Sigma+,pi-", 1. / 3); + makeStaticData()->AddDecay("Lambda(1405) --> Sigma- + pi+", "Lambda1405", "Sigma-,pi+", 1. / 3); + makeStaticData()->AddDecay("Lambda(1405) --> Sigma0 + pi0", "Lambda1405", "Sigma0,pi0", 1. / 3); + } + + if (hyperon_photon_decays) { + makeStaticData()->AddDecay("Lambda(1405) --> Lambda + gamma", "Lambda1405", "Lambda,g", 0.0085); + makeStaticData()->AddDecay("Lambda(1405) --> Sigma(1385)0 + gamma", "Lambda1405", "Sigma13850,g", 0.0085); } if (hyperon_dalitz_decays) { - makeStaticData()->AddDecay("Sigma(1385)0 --> Lambda + dilepton", "Sigma13850", "Lambda,dilepton", 0.013 / 137.); - PResonanceDalitz * dalitz = new PResonanceDalitz("Sigma13850_dalitz@Sigma13850_to_Lambda_dilepton", "dgdm from Zetenyi/Wolf", -1); + makeStaticData()->AddDecay("Lambda(1405) --> Lambda + dilepton", "Lambda1405", "Lambda,dilepton", 0.0085 / 137.); + + PResonanceDalitz * dalitz = new PResonanceDalitz("Lambda1405_dalitz@Lambda1405_to_Lambda_dilepton", "dgdm from Zetenyi/Wolf", -1); dalitz->setGm(0.719); makeDistributionManager()->Add(dalitz); } @@ -154,24 +156,22 @@ Bool_t PStrangenessPlugin::ExecCommand(const char *command, Double_t) { makeStaticData()->SetParticleParity("Lambda1520", 1); if (hyperon_hadron_decays) { - makeStaticData()->AddDecay("Lambda(1520) --> n + K0S", "Lambda1520", "n,K0S", 0.45 / 2); // N + /K from PDG - makeStaticData()->AddDecay("Lambda(1520) --> n + K0L", "Lambda1520", "n,K0L", 0.45 / 2); // N + /K from PDG - makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi0 + pi0", "Lambda1520", "Lambda,pi0,pi0", 0.2 / 3); // - - makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi+ + pi-", "Lambda1520", "Lambda,pi+,pi-", 0.1 / 3); // old BR .014638 + makeStaticData()->AddDecay("Lambda(1520) --> n + K0S", "Lambda1520", "n,K0S", 0.45 / 3); // N + /K from PDG + makeStaticData()->AddDecay("Lambda(1520) --> n + K0L", "Lambda1520", "n,K0L", 0.45 / 3); // N + /K from PDG + makeStaticData()->AddDecay("Lambda(1520) --> p + K-", "Lambda1520", "p,K-", 0.45 / 3); // N + /K from PDG + makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi0 + pi0", "Lambda1520", "Lambda,pi0,pi0", 0.2 / 3); // 10% * 2/3 + makeStaticData()->AddDecay("Lambda(1520) --> Lambda + pi+ + pi-", "Lambda1520", "Lambda,pi+,pi-", 0.1 / 3); // 10% * 1/3 , old BR .014638 + makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi0", "Lambda1520", "Sigma0,pi0", 0.139096); // 42% / 3 + makeStaticData()->AddDecay("Lambda(1520) --> Sigma+ + pi-", "Lambda1520", "Sigma+,pi-", 0.139096); // 42% / 3 + makeStaticData()->AddDecay("Lambda(1520) --> Sigma- + pi+", "Lambda1520", "Sigma-,pi+", 0.139096); // 42% / 3 + makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi0 + pi0", "Lambda1520", "Sigma0,pi0,pi0", 0.009/5.0); // + makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi+ + pi-", "Lambda1520", "Sigma0,pi+,pi-", 0.009/5.0); // } if (hyperon_photon_decays) { makeStaticData()->AddDecay("Lambda(1520) --> Lambda + g", "Lambda1520", "Lambda,g", 0.0085); } - if (hyperon_experimental_decays) { - makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi0", "Lambda1520", "Sigma0,pi0", 0.139096); // - - makeStaticData()->AddDecay("Lambda(1520) --> Sigma+ + pi-", "Lambda1520", "Sigma+,pi-", 0.139096); // - - makeStaticData()->AddDecay("Lambda(1520) --> Sigma- + pi+", "Lambda1520", "Sigma-,pi+", 0.139096); // - - makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi0 + pi0", "Lambda1520", "Sigma0,pi0,pi0", 0.009/5.0); // - - makeStaticData()->AddDecay("Lambda(1520) --> Sigma0 + pi+ + pi-", "Lambda1520", "Sigma0,pi+,pi-", 0.009/5.0); // - - } - if (hyperon_dalitz_decays) { makeStaticData()->AddDecay("Lambda(1520) --> Lambda + dilepton", "Lambda1520", "Lambda,dilepton", .007948/137.);