Run Card Cuts are not working for charged leptons

Hi,

What is the value that you choose for the following three variables:
ickkw
auto_ptj_mjj
xqcut

The default for ickkw should be 1 (MLM merging) for your process definition
but 0 is also possible if you do CKKW-L merging.
If it is 1 and the auto_ptj_mjj is on t=True (default as well)
then the value of ptj is not used during the computation and you only havve to specify xqcut.

So this can explains the issue for jet.

For the lepton issue this is more weird.

> Similarly for jet and photons also. It is taking the default numbers.

I do not understand how you have photon plot since you do not have photon.
Did you make the plot after parton-shower?
The cut in the run_card are performed before parton-shower.

Which version of MG5aMC are you using? in principle the line
> 20.0 = pta ! minimum pt for the photons
should not even be present in the run_card for your process.

Cheers,

Olivier

> On 26 Aug 2021, at 10:50, Samadrita Mukherjee <email address hidden> wrote:
>
> New question #698516 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/698516
>
> Dear Madgraph Team,
>
> I am generating a process like the following:
>
> generate p p > mu+ vm
> add process p p > mu+ vm j
> add process p p > mu+ vm j j
>
> and similarly for mu- separately.
>
> Now I want to impose cuts on my final state muons. For this in the run_card.dat, I did the following,
>
> (I changed the jet and photon pT cuts also to check whether it is working)
>
>
> ----------------------------------------------------------------------------------------------------------------------------------------------
> True = cut_decays ! Cut decay products
> #*********************************************************************
> # Standard Cuts *
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 30.0 = ptj ! minimum pt for the jets
> 0.0 = ptb ! minimum pt for the b
> 20.0 = pta ! minimum pt for the photons
> 30.0 = ptl ! minimum pt for the charged leptons
> 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
> -1.0 = ptjmax ! maximum pt for the jets
> -1.0 = ptbmax ! maximum pt for the b
> -1.0 = ptamax ! maximum pt for the photons
>
> ---------------------------------------------------------------------------------------------------------------------------------------------
>
> But when I am making the pT distribution plots it shows the default cuts had been applied only i.e. for jets it is 20 GeV and for charged leptons, it is 10 GeV.
>
> I mean, I have events with pT of charged leptons < 30 GeV, even after putting
>
> 30.0 = ptl ! minimum pt for the charged leptons
>
> Similarly for jet and photons also. It is taking the default numbers.
>
> Can you please tell me if we need to do any other change somewhere to implement the cuts?
>
> Sincerely,
> Samadrita
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Additionally,

I have run the following code (in the Long Term Stable version: currently 2.9.5):
import model sm
generate p p > mu+ vm
add process p p > mu+ vm j
add process p p > mu+ vm j j
output
launch
set ptl 30

and the cut on the muon works as expected.

Cheers,

Olivier

Dear Olivier,

Thanks for your reply. I am using v2_6_7 (a bit older version)

I am using the matching blocks as the default one:

#*********************************************************************
# Matching parameter (MLM only)
#*********************************************************************
 1 = ickkw ! 0 no matching, 1 MLM
 1.0 = alpsfact ! scale factor for QCD emission vx
 False = chcluster ! cluster only according to channel diag
 4 = asrwgtflavor ! highest quark flavor for a_s reweight
 True = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
                                   ! (turn off for VBF and single top processes)
 30.0 = xqcut ! minimum kt jet measure between partons

______________________________________________________________________________

Thanks,
Samadrita

So ptj is not used and the important parameter is xqcut for the jet.

For the rest, can you test with 2.9.5?
I do not have that issue there so that should solve your issue.

Cheers,

Olivier

> On 26 Aug 2021, at 12:00, Samadrita Mukherjee <email address hidden> wrote:
>
> Question #698516 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/698516
>
> Status: Answered => Open
>
> Samadrita Mukherjee is still having a problem:
> Dear Olivier,
>
> Thanks for your reply. I am using v2_6_7 (a bit older version)
>
> I am using the matching blocks as the default one:
>
> #*********************************************************************
> # Matching parameter (MLM only)
> #*********************************************************************
> 1 = ickkw ! 0 no matching, 1 MLM
> 1.0 = alpsfact ! scale factor for QCD emission vx
> False = chcluster ! cluster only according to channel diag
> 4 = asrwgtflavor ! highest quark flavor for a_s reweight
> True = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
> ! (turn off for VBF and single top processes)
> 30.0 = xqcut ! minimum kt jet measure between partons
>
> ______________________________________________________________________________
>
>
> Thanks,
> Samadrita
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Thanks. I am going to try it soon.

Did it work?

Olivier

Yeah, it worked properly, sorry I forgot to update you. Thanks a lot.