MadGraph5_aMC@NLO

deltaeta cut not applied?

Asked by Will Flanagan

Hi Madgraph Experts,

It seems that the deltaeta cut is not being applied for me. I am doing some simple SUSY production:
mg5>import model mssm; generate p p > p p t1 t1~; output testdir; launch

I am running with the 'default' run card except that I require an invariant dijet mass cut (mmjj) and a delta eta cut between jets (deltaeta). Since I am stopping at the parton level, the dijet definition is unambiguous; there are only two jets and two non-jet SUSY particles in my final state.

Nonetheless, when I plot dETA(j[1] j[2]) in MadAnalysis, I am seeing events in between -deltaeta and +deltaeta... So it seems the deltaeta cut is not being fully applied.

Does anyone have suggestions for some simple mistakes I might be making? I have checked my LHE files and I have checked that MadAnalysis agrees that each event has two and only two jets; the SUSY particles are not contributing to this distribution.

Any suggestions are much appreciated!
Will

P.S. My MadAnalysis script looks something like this:

>bin/ma5 script.txt

######
import testdir/Events/run_01/events.lhe.gz as dummy
plot dETA(j[1] j[2])
set selection[1].logY = true
set selection[1].stacking_method = normalize2one
submit dummy
open dummy/DVI
#######

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:

This question was reopened

Revision history for this message
Will Flanagan (will-flanagan) said : 		#1

Hi All,

It was a silly mistake that I was able to solve myself. I had not set 'auto_ptj_mjj' to 'F' as the comments say to do.

By the way, is there any documentation on what this switch does?

Thanks,
Will

Revision history for this message
Will Flanagan (will-flanagan) said : 		#2

Hi All,

Actually, I am still getting events with dETA(j[1] j[2]) in between -deltaeta and +deltaeta even with setting auto_ptj_mjj to 'F'. This didn't fix the problem. Are there any suggestions on how to actually apply this deltaeta cut?

Thanks,
Will

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

Hi Will,

The cut implemented in MG is in fact the following:
            if (abs(rap(p(0,hardj1))-rap(p(0,hardj2))) .lt. deltaeta) then
             passcuts=.false.
             return

I don't know how dETA(j[1] j[2]) is define in MadAnalysis.

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Will Flanagan for more information if necessary.

To post a message you must log in.