Simulating muon decay with new particle

Hi Bertrand,

> - what is the best way to implement this decay

The first question is do you have a Lagrangian for your model?
If yes, then the easiest is to use FeynRules to create a UFO model.
If not, then you need to create the model by hand by modifying
either the UFO or the _v4 one.

> - how do I code the cuts (and which ones need to be hardcoded). I should remember enough fortran to implement any receipt.

In fact this is pretty easy:
in the file SubProcesses/setcuts.f
at the line 159, you have the line
if(nincoming.eq.1) do_cuts(i)=.false.
you can pass it to
if(nincoming.eq.1) do_cuts(i)=.true.
And then all the cuts define in the run_card.dat will be applied.
Note that I don't know how the eta cut are done in that case, so you need to put those one to zero or dig into the code to see if they apply for decay or not.

Cheers,

Olivier

On Nov 18, 2013, at 6:36 PM, Bertrand Echenard <email address hidden> wrote:

> New question #239421 on MadGraph5:
> https://answers.launchpad.net/madgraph5/+question/239421
>
> Dear Johan, Olivier,
>
> I am trying to simulate a radiative muon decay involving a new particle:
>
> mu -> e nu_e nu_mu A', A' -> e+e-
>
> where A' is a new massive spin-1 boson (hidden photon).
>
> In a previous post, Olivier mentionned that the muon decays are well simulated within the sm-full model, but the cuts are open and need to be hardcoded. So I have two questions:
> - what is the best way to implement this decay
> - how do I code the cuts (and which ones need to be hardcoded). I should remember enough fortran to implement any receipt.
>
> Thanks for your help.
>
> Best regards,
>
> Bertrand
>
>
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.