filter high pt events

Asked by Sergei Chekanov

 Hello,

 For our studies, we would like to generate ttbar events and keep (filter) only events with large pT of top quarks, before filling lhe files. Is it possible to do this? Or even better, generate events with a min cut on inv. mass of 2->2 system, as this can be done in Pythia or Herwig. thanks , Sergei

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

You can use ptheavy in the run_card for that.

Cheers,

Olivier
On Apr 8, 2014, at 9:51 PM, Sergei Chekanov <email address hidden> wrote:

> New question #246709 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/246709
>
>
> Hello,
>
> For our studies, we would like to generate ttbar events and keep (filter) only events with large pT of top quarks, before filling lhe files. Is it possible to do this? Or even better, generate events with a min cut on inv. mass of 2->2 system, as this can be done in Pythia or Herwig. thanks , Sergei
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Sergei Chekanov (chekanov) said :
#2

 Hi, Oliver

 I've tried this using this for MG5_aMC_v2_0_0 and "generate p p > t t~ [QCD]" (NLO QCD)
 and setting the cut pT(top)>2500 GeV as:

  0 = pta ! minimum pt for the photons
  0 = ptl ! minimum pt for the charged leptons
  0 = misset ! minimum missing Et (sum of neutrino's momenta)
 2500.0 = ptheavy ! minimum pt for one heavy final state

this did not change the output files that still have all top pT's

best, Sergei

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

Hi Sergei,

Indeed this cut is only define for LO process.
You need to modify by hand the cuts.f file if you want to implement this cut.
This file is now documented to try to simplify this.

Cheers,

Olivier

On Apr 9, 2014, at 4:36 PM, Sergei Chekanov <email address hidden> wrote:

> Question #246709 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246709
>
> Status: Answered => Open
>
> Sergei Chekanov is still having a problem:
>
> Hi, Oliver
>
> I've tried this using this for MG5_aMC_v2_0_0 and "generate p p > t t~ [QCD]" (NLO QCD)
> and setting the cut pT(top)>2500 GeV as:
>
> 0 = pta ! minimum pt for the photons
> 0 = ptl ! minimum pt for the charged leptons
> 0 = misset ! minimum missing Et (sum of neutrino's momenta)
> 2500.0 = ptheavy ! minimum pt for one heavy final state
>
> this did not change the output files that still have all top pT's
>
> best, Sergei
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Can you help with this problem?

Provide an answer of your own, or ask Sergei Chekanov for more information if necessary.

To post a message you must log in.