Skip to content

Conversation

@SabrinaBrickner
Copy link

Description

Adds 4 files to the MeVPrtl generator. 2 files are the new tools required to create and decay LLPs, and the other 2 files are the relevant CMakeLists. The Kaon2LLPFlux_tool.cc is in the process of validation and is being pushed so that other collaborators can work on it. The LLPMakeDecay_tool.cc is not complete, but is pushed to allow the rest of the workflow to operate properly.

  • Have you added a label? (bug/enhancement/physics etc.)
  • Have you assigned at least 1 reviewer?
  • [] Is this PR related to an open issue / project?
  • [] Does this PR affect CAF data format? If so, please assign a CAF maintainer as additional reviewer.
  • Does this PR require merging another PR in a different repository (such as sbnanobj/sbnobj etc.)? If so, please link it in the description. -> in sbndcode it requires "MeVPrtl LLP (fcls and files)"
  • [] Are you submitting this PR on behalf of someone else who made the code changes? If so, please mention them in the description.

@SabrinaBrickner SabrinaBrickner added the enhancement New feature or request label Oct 7, 2025
* @brief Kaon2LLPFlux class definiton
*
* Implementation of Kaon->LLP branching ratio taken from:
* arXiv:1912.07622
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

What equation do you use from arXiv:1912.07622?


// Branch the parent Kaon Decay
double llp_mass = fM;
double br = 1;
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This should be configurable with a fcl parameter in case we want to add multiple decay modes: K --> mu + LLP, K --> e + LLP...

// Branch the parent Kaon Decay
double llp_mass = fM;
double br = 1;
bool is_muon = true;
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This should be configurable with a fcl parameter, and add the support in case it is a pion instead.

llp.mass = fM;

llp.meson_pdg = kaon.meson_pdg;
llp.secondary_pdg = (is_muon ? 13 : 11) * (kaon.meson_pdg > 0 ? 1 : -1);
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Add support for pion case


llp.meson_pdg = kaon.meson_pdg;
llp.secondary_pdg = (is_muon ? 13 : 11) * (kaon.meson_pdg > 0 ? 1 : -1);
llp.generator = 1; // kLLP
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This code corresponds to the HNL. Add a new one to sbnobj/Common/EventGen/MeVPrtl/MeVPrtlTruth.h

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I did not review this module.

@jicrespo
Copy link
Member

For the record, this is an ongoing review and it is not required to converge for the SBND Fall Production.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request

Projects

Status: Open pull requests

Development

Successfully merging this pull request may close these issues.

4 participants