MadGraph5_aMC@NLO

Using standard cuts for non-standard particles in MadGraph

Asked by J. Heisig

Hi,

I'd like to cut on the invariant mass of two DM particles (a majorana fermion with some lengthy particle ID, say a SUSY like '1000022'). As I don't intend to cut on any photons (I also tried b-quarks) I though that the easiest way to achieve this is to declare that 1000022 is a photon (or b-quark) by writing

 if (idup(i,1,iproc).eq.1000022) is_a_a(i)=.true. ! photon

around line 217 in setcuts.f and using the standard mmaa cut from the run_card. But this does not work. Why?

At least the cut on the photon does not work anymore, meaning that my modification does something... I didn't found any other place where the code could actually recognize that 1000022 is not a photon.

Thanks a lot
Jan

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Solved by:
marco zaro
Solved:
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said : 		#1

Hi,
indeed, that should be a good way to proceed... Are you sure everything is recompiled afterwads (time stamp of the .o files > time stamp of setcuts.f)?
Another thing you can do is to edit SubProcesses/cuts.f and add your own cuts in there.
Let us know!

cheers,

Marco

Revision history for this message
J. Heisig (heisig) said : 		#2

Hi Marco,

I used a "fresh" version of MG5_aMC_v2.2.3 and changed setcuts.f in the very beginning. (Before running MG the first time, i.e. generating the process and launching it.) So I don't see how this could be a compiler issue. Somehow it seems to recognize that 1000022 is not a photon, but how?

Cheers
Jan

Revision history for this message
marco zaro (marco-zaro) said : 		#3

Hi,
are you changing setcuts inside Template/LO/SubProcesses/ ?
How are you changing it?
can you attach it or send to my private email address?

Thanks,

Marco

On 14 Jul 2015, at 18:21, J. Heisig <email address hidden> wrote:

> Question #269098 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/269098
>
> Status: Answered => Open
>
> J. Heisig is still having a problem:
> Hi Marco,
>
> I used a "fresh" version of MG5_aMC_v2.2.3 and changed setcuts.f in the
> very beginning. (Before running MG the first time, i.e. generating the
> process and launching it.) So I don't see how this could be a compiler
> issue. Somehow it seems to recognize that 1000022 is not a photon, but
> how?
>
> Cheers
> Jan
>
> --
> You received this question notification because you are a direct
> subscriber of the question.

Revision history for this message
J. Heisig (heisig) said : 		#4

Hi Marco,

> are you changing setcuts inside Template/LO/SubProcesses/ ?

Yes, exactly.

> How are you changing it?

I added

         if (idup(i,1,iproc).eq.1000022) is_a_a(i)=.true. ! neutralino

below

         if (idup(i,1,iproc).eq.22) is_a_a(i)=.true. ! photon

> can you attach it or send to my private email address?

Here you go:

Thanks a lot!
Cheers,
Jan

--
Jan Heisig
Institute for Theoretical Particle Physics and Cosmology
RWTH Aachen University
52056 Aachen
Germany

Phone +49 (0)241 80 27060
Email <email address hidden>

Revision history for this message
Best marco zaro (marco-zaro) said : 		#5

Hi Jan,
so I add here the trick that solved your problem

inside setcuts.f there is this line (187)

        if (pmass(i).gt.20d0) do_cuts(i)=.false. ! no cuts on top,W,Z,H

where you should put an .and.idup(i,1,iproc).ne.1000022 because otherwise cuts on heavy particles won’t be applied…

Can you please confirm that the question is solved?

Thanks again!

Marco

Revision history for this message
J. Heisig (heisig) said : 		#6

Thanks marco zaro, that solved my question.