Cuts: add maximum ptheavy (run_card.dat, cuts.f)

Asked by Jory Sonneveld

I would like to be able to set a maximum for ptheavy in the run card (run_card.dat). In the run card there is only the option for a ptheavy minimum. I guess I could change the cuts.f file for this in the SubProcesses directory, but perhaps there is another way. What is the best way to add a maximum ptheavy for a run in which I generate events?

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Johan Alwall
Solved:
Last query:
Last reply:
Revision history for this message
Best Johan Alwall (johan-alwall) said :
#1

Hello Jory,

Unfortunately we don't have a maxcut for ptheavy implemented. You can fairly easily implement one yourself however. If you want to be able to set the cut in the run_card.dat, you need to change the files Source/setrun.f and SubProcesses/cuts.f. You can model the new cut on the existing ptheavy cut.

All the best,
Johan

Revision history for this message
Jory Sonneveld (jory) said :
#2

Thanks, I found the files and edited them.

Revision history for this message
Jory Sonneveld (jory) said :
#3

Thanks Johan Alwall, that solved my question.

Revision history for this message
Roberto Franceschini (franceschini-roberto) said :
#4

Hi, I just wanted to add my 2 cents. At the card currently stands there is only ptheavy missing to have max cut - and it would be very useful if it were present in the template process out of the box.
Cheers,
Roberto

# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 20.0 = ptj ! minimum pt for the jets
 0.0 = ptb ! minimum pt for the b
 10.0 = pta ! minimum pt for the photons
 10.0 = ptl ! minimum pt for the charged leptons
 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
 0.0001 = ptheavy ! minimum pt for one heavy final state
 -1.0 = ptjmax ! maximum pt for the jets
 -1.0 = ptbmax ! maximum pt for the b
 -1.0 = ptamax ! maximum pt for the photons
 -1.0 = ptlmax ! maximum pt for the charged leptons
 -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#5

Hi,

Ptheavy is actually already different of the other cuts in this category and should not be in this category of cut.
In the next version of the code, you will find it under the "Inclusive cuts" category.

Now, for a while, I would like to implement [simple] cuts depending of the pdg value. Something like:
   0. = ptmin_25
   -1. = ptmax_25

But I can not tell when I will implement this idea.

Cheers,

Olivier

> On 28 Jul 2017, at 10:46, Roberto Franceschini <email address hidden> wrote:
>
> Question #208439 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/208439
>
> Roberto Franceschini posted a new comment:
> Hi, I just wanted to add my 2 cents. At the card currently stands there is only ptheavy missing to have max cut - and it would be very useful if it were present in the template process out of the box.
> Cheers,
> Roberto
>
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 20.0 = ptj ! minimum pt for the jets
> 0.0 = ptb ! minimum pt for the b
> 10.0 = pta ! minimum pt for the photons
> 10.0 = ptl ! minimum pt for the charged leptons
> 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
> 0.0001 = ptheavy ! minimum pt for one heavy final state
> -1.0 = ptjmax ! maximum pt for the jets
> -1.0 = ptbmax ! maximum pt for the b
> -1.0 = ptamax ! maximum pt for the photons
> -1.0 = ptlmax ! maximum pt for the charged leptons
> -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.