How to generate QCD multijet event in PP collision with HTbin like HT[100-200,200-300,300-400,--------1500-2000,2000-inf] and merge them with proper weight

Asked by suman

Here I try to generate p p >j j and add process p p > j j j for MLM matching.
 but for proper pt distribution I have to use the ptj_bias as
{'ptj_bias_target_ptj': 7000.0, 'ptj_bias_enhancement_power': 4.0} = bias_parameters

now if I want to generate without bias then how to control the phase space bined sample generation

like HT=100-200,200-300,300-500,500-700, --- --,1500-2000,2000-inf

further how to take care of weight for each HT bin ?

Thanks in advance.

suman kundu

Question information

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

hi,

> now if I want to generate without bias then how to control the phase space bined sample generation
>
> like HT=100-200,200-300,300-500,500-700, --- --,1500-2000,2000-inf

You can use the cut in the run_card for that.

> further how to take care of weight for each HT bin ?

Nothing specific here.each event will have his weight (the cross-section of the associated sample)
that you have to keep when you merge the sample.

Cheers,

Olivier

> On 17 Aug 2018, at 15:52, suman <email address hidden> wrote:
>
> New question #672421 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/672421
>
> Here I try to generate p p >j j and add process p p > j j j for MLM matching.
> but for proper pt distribution I have to use the ptj_bias as
> {'ptj_bias_target_ptj': 7000.0, 'ptj_bias_enhancement_power': 4.0} = bias_parameters
>
> now if I want to generate without bias then how to control the phase space bined sample generation
>
> like HT=100-200,200-300,300-500,500-700, --- --,1500-2000,2000-inf
>
> further how to take care of weight for each HT bin ?
>
> Thanks in advance.
>
> suman kundu
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
suman (sumanmg5) said :
#2

Hello olivier ,

You mean I can control the phase space with the following parameters

# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
-1.0 = htjmax ! maximum jet HT=Sum(jet pt)

or I have to update it in any portion of the run_card.dat

Thanks for early reply
best regards,

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

> You mean I can control the phase space with the following parameters
>
> # Control the Ht(k)=Sum of k leading jets *
> #*********************************************************************
> 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1.0 = htjmax ! maximum jet HT=Sum(jet pt)

Yes those are the parameter that you can edit.
Obviously those are cut BEFORE the parton-shower so after you run pythia those cuts
will not be full-filled anymore.

Cheers,

Olivier

> On 17 Aug 2018, at 17:02, suman <email address hidden> wrote:
>
> Question #672421 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/672421
>
> Status: Answered => Open
>
> suman is still having a problem:
> Hello olivier ,
>
> You mean I can control the phase space with the following parameters
>
> # Control the Ht(k)=Sum of k leading jets *
> #*********************************************************************
> 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
>
>
> or I have to update it in any portion of the run_card.dat
>
> Thanks for early reply
> best regards,
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
suman (sumanmg5) said :
#4

Ok thats fine.

Thanks 0livier for your kind help and suggestion.

Best Regards