Final sample is not complete after pythia and delphes in ep collision
Dear Experts,
I am trying to generate SM below process with MadGraph 2.7.3, pythia8 and Delphes inside the MadGraph:
define p = p b b~
generate e- p > vl j b b~
In the general SM model (import model sm), it crashes. I tried with both:
import model sm-full
In this case I have used False = use_syst and nb_core = 1 to solve some problems especially pythia parallelization issue
second case:
import model sm-no_b_mass
After all, one problem still remained:
Number of data after delphes in the final root file is less than I generated in madgraph. e.g. about 200 event out of 10000.
I would be appreciated for any help.
Thanks,
Reza
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Revision history for this message
|
#1 |
Hi,
How many event did you have after MG5aMC? and after Pythia?
Did you contact Delphes author?
I would make sense Delphes has detector acceptance, did you change those to march a DIS type of detector?
Cheers,
Olivier
> On 14 Feb 2021, at 08:55, Reza Jafari <email address hidden> wrote:
>
> New question #695526 on MadGraph5_aMC@NLO:
> https:/
>
> Dear Experts,
>
> I am trying to generate SM below process with MadGraph 2.7.3, pythia8 and Delphes inside the MadGraph:
>
> define p = p b b~
> generate e- p > vl j b b~
>
> In the general SM model (import model sm), it crashes. I tried with both:
>
> import model sm-full
> import model sm-no_b_mass
>
> After solving some problems, one problem remained:
> Number of data after delphes in the final root file is less than I generated in madgraph. e.g. about 200 event out of 10000.
>
> I would be apprecieted for any help.
>
> Thanks,
> Reza
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#2 |
Dear Olivier,
Thanks for your prompt answer.
I just checked again now.
In the last run, the number of events after MG5aMC level for sm-full is correct (10k), but after Pythia is incorrect (0 event).
for sm-no_b_mass the number of events after both MG5aMC (90 events out of 10k) and Pythia (68 event) are wrong, as you can see below.
I just turned on Delphes in MG5aMC shell, nothing more.
--
sm-no_b_mass:
*------- PYTHIA Event and Cross Section Statistics -------
| |
| Subprocess Code | Number of events | sigma +- delta |
| | Tried Selected Accepted | (estimated) (mb) |
| | | |
|-----
| | | |
| Les Houches User Process(es) 9999 | 90 90 68 | 2.058e-07 0.000e+00 |
| ... whereof user classification code 1 | 90 90 68 | |
| | | |
| sum | 90 90 68 | 2.058e-07 0.000e+00 |
| |
*------- End PYTHIA Event and Cross Section Statistics -------
sm-full:
*------- PYTHIA Event and Cross Section Statistics -------
| |
| Subprocess Code | Number of events | sigma +- delta |
| | Tried Selected Accepted | (estimated) (mb) |
| | | |
|-----
| | | |
| | | |
| sum | 0 0 0 | 0.000e+00 0.000e+00 |
| |
*------- End PYTHIA Event and Cross Section Statistics -------
Revision history for this message
|
#3 |
What cut/setup did you use for MG5aMC.
Can you give me your full pythia8_card since i do not know how to setup it for dIS type of production.
Cheers,
Olivier
Revision history for this message
|
#4 |
Dear Olivier,
Please see full run_tag_banner file in the following.
Thanks for your time.
Regards,
Reza
For sm-full model:
--
<LesHouchesEvents version="3.0">
<header>
<!--
#******
# *
# MadGraph5_aMC@NLO *
# *
# Going Beyond *
# *
# http://
# http://
# http://
# *
# The MadGraph5_aMC@NLO team *
# *
#......
# *
# This file contains all the information necessary to reproduce *
# the events generated: *
# *
# 1. software version *
# 2. proc_card : code generation info including model *
# 3. param_card : model primary parameters in the LH format *
# 4. run_card : running parameters (collider and cuts) *
# 5. pythia_card : present only if pythia has been run *
# 6. pgs_card : present only if pgs has been run *
# 7. delphes_cards : present only if delphes has been run *
# *
# *
#******
-->
<MGVersion>
2.7.3
</MGVersion>
<MG5ProcCard>
<![CDATA[
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.7.3 2020-06-21 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set group_subprocesses Auto
set ignore_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model sm-full
generate e- p > vl j b b~
output bkg1_test_sm_full
]]>
</MG5ProcCard>
<MGProcCard>
#******
# MadGraph/MadEvent *
# http://
# *
# proc_card.dat *
#******
# *
# This Files is generated by MADGRAPH 5 *
# *
# WARNING: This Files is generated for MADEVENT (compatibility issue)*
# This files is NOT a valid MG4 proc_card.dat *
# Running this in MG4 will NEVER reproduce the result of MG5*
# *
#******
#******
# Process(es) requested : mg2 input *
#******
# Begin PROCESS # This is TAG. Do not modify this line
e- p > vl j b b~ #Process
# Be carefull the coupling are here in MG5 convention
end_coup # End the couplings input
done # this tells MG there are no more procs
# End PROCESS # This is TAG. Do not modify this line
#******
# Model information *
#******
# Begin MODEL # This is TAG. Do not modify this line
sm-full
# End MODEL # This is TAG. Do not modify this line
#******
# Start multiparticle definitions *
#******
# Begin MULTIPARTICLES # This is TAG. Do not modify this line
# End MULTIPARTICLES # This is TAG. Do not modify this line
</MGProcCard>
<MGRunCard>
<![CDATA[
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update full_run_card *
#******
#
#******
# Tag name for the run (one word) *
#******
tag_1 = run_tag ! name of the run
#******
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#******
10000 = nevents ! Number of unweighted events requested
27 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#******
0 = lpp1 ! beam 1 type
1 = lpp2 ! beam 2 type
60.0 = ebeam1 ! beam 1 total energy in GeV
6500.0 = ebeam2 ! beam 2 total energy in GeV
#******
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#******
0.0 = polbeam1 ! beam polarization for beam 1
0.0 = polbeam2 ! beam polarization for beam 2
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
nn23lo1 = pdlabel ! PDF set
230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#******
# Renormalization and factorization scales *
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.188 = scale ! fixed ren scale
91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
-1 = dynamical_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format
#******
False = gridpack !True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"
#******
#
#******
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#******
0 = nhel ! using helicities importance sampling or not.
#******
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.
#******
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#******
# Parton level cuts definition *
#******
#
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# Standard Cuts *
#******
# 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
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 = missetmax ! maximum missing Et (sum of neutrino's momenta)
{} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
{} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#******
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
5.0 = etaj ! max rap for the jets
-1.0 = etab ! max rap for the b
0.0 = etabmin ! min rap for the b
{} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
{} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#******
# Minimum and maximum DeltaR distance *
#******
0.0 = drbb ! min distance between b's
0.0 = drbj ! min distance between b and jet
-1.0 = drbbmax ! max distance between b's
-1.0 = drbjmax ! max distance between b and jet
#******
# Minimum and maximum invariant mass for pairs *
#******
0.0 = mmbb ! min invariant mass of a b pair
-1.0 = mmbbmax ! max invariant mass of a b pair
{} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
{'default': False} = mxx_only_
#******
# Inclusive cuts *
#******
0.0 = xptb ! minimum pt for at least one b
#*****
# Control the Ht(k)=Sum of k leading jets *
#*****
0.0 = ihtmin !inclusive Ht for all partons (including b)
-1.0 = ihtmax !inclusive Ht for all partons (including b)
0.0 = ht2min ! minimum Ht for the two leading jets
0.0 = ht3min ! minimum Ht for the three leading jets
-1.0 = ht2max ! maximum Ht for the two leading jets
-1.0 = ht3max ! maximum Ht for the three leading jets
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
#
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
False = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_
# Syscalc is deprecated but to see the associate options type'update syscalc'
]]>
</MGRunCard>
<slha>
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
#######
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
#######
#######
## INFORMATION FOR MASS
#######
Block mass
4 1.270000e+00 # MC
5 4.700000e+00 # MB
6 1.720000e+02 # MT
11 5.110000e-04 # Me
13 1.056600e-01 # MM
15 1.777000e+00 # MTA
23 9.118800e+01 # MZ
25 1.250000e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
1 0.000000e+00 # d : 0.0
2 0.000000e+00 # u : 0.0
3 0.000000e+00 # s : 0.0
12 0.000000e+00 # ve : 0.0
14 0.000000e+00 # vm : 0.0
16 0.000000e+00 # vt : 0.0
21 0.000000e+00 # g : 0.0
22 0.000000e+00 # a : 0.0
24 8.041900e+01 # w+ : cmath.sqrt(
#######
## INFORMATION FOR SMINPUTS
#######
Block sminputs
1 1.325070e+02 # aEWM1
2 1.166390e-05 # Gf
3 1.180000e-01 # aS
#######
## INFORMATION FOR WOLFENSTEIN
#######
Block wolfenstein
1 2.253000e-01 # lamWS
2 8.080000e-01 # AWS
3 1.320000e-01 # rhoWS
4 3.410000e-01 # etaWS
#######
## INFORMATION FOR YUKAWA
#######
Block yukawa
4 1.270000e+00 # ymc
5 4.200000e+00 # ymb
6 1.645000e+02 # ymt
11 5.110000e-04 # yme
13 1.056600e-01 # ymm
15 1.777000e+00 # ymtau
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.508336e+00 # WT
DECAY 15 2.270000e-12 # WTau
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
</slha>
<MGPythiaCard>
<![CDATA[
!
! It is possible to run this card manually with:
! LD_LIBRARY_
!
!
! Pythia8 cmd card automatically generated by MadGraph5_aMC@NLO
! For more information on the use of the MG5aMC / Pythia8 interface, visit
! https:/
!
! ==================
! General parameters
! ==================
!
Main:numberOfEvents = -1
!
! -------
! Specify the HEPMC output of the Pythia8 shower. You can set it to:
! auto : MG5aMC will automatically place it the run_<i> directory
! /dev/null : to turn off the HEPMC output.
! <path> : to select where the HEPMC file must written. It will
! therefore not be placed in the run_<i> directory. The
! specified path, if not absolute, will be relative to
! the Event/run_<i> directory of the process output.
! fifo : to have MG5aMC setup the piping of the PY8 output to
! analysis tools such as MadAnalysis5.
! fifo@<fifo_path> :
! Same as 'fifo', but selecting a custom path to create the
! fifo pipe. (useful to select a mounted drive that supports
! fifo). Note that the fifo file extension *must* be '.hepmc.fifo'.
! -------
!
HEPMCoutput:file = tag_1_pythia8_
!
! -------
! Parameters relevant only when performing MLM merging, which can be
! turned on by setting ickkw to '1' in the run_card and chosing a
! positive value for the parameter xqcut.
! For details, see section 'Jet Matching' on the left-hand menu of
! http://
! -------
! If equal to -1.0, MadGraph5_aMC@NLO will set it automatically based
! on the parameter 'xqcut' of the run_card.dat
! The following parameter was forced to be commented out by MG5aMC.
! JetMatching:qCut = -1.0
! Use default kt-MLM to match parton level jets to those produced by the
! shower. But the other Shower-kt scheme is available too with this option.
! The following parameter was forced to be commented out by MG5aMC.
! JetMatching:
! A value of -1 means that it is automatically guessed by MadGraph.
! It is however always safer to explicitly set it.
! The following parameter was forced to be commented out by MG5aMC.
! JetMatching:nJetMax = -1
!
! -------
! Parameters relevant only when performing CKKW-L merging, which can
! be turned on by setting the parameter 'ptlund' *or* 'ktdurham' to
! a positive value.
! For details, see section 'CKKW-L Merging' on the left-hand menu of
! http://
! -------
! Central merging scale values you want to be used.
! If equal to -1.0, then MadGraph5_aMC@NLO will set this automatically
! based on the parameter 'ktdurham' of the run_card.dat
! The following parameter was forced to be commented out by MG5aMC.
! Merging:TMS = -1.0
! This must be set manually, according to Pythia8 directives.
! An example of possible value is 'pp>LEPTONS,
! Alternatively, from Pythia v8.223 onwards, the value 'guess' can be
! used to instruct Pythia to guess the hard process. The guess would mean
! that all particles apart from light partons will be considered as a part
! of the hard process. This guess is prone to errors if the desired hard
! process is complicated (i.e. contains light partons). The user should
! then be wary of suspicious error messages in the Pythia log file.
! The following parameter was forced to be commented out by MG5aMC.
! Merging:Process = <set_by_user>
! A value of -1 means that it is automatically guessed by MadGraph.
! It is however always safer to explicitly set it.
! The following parameter was forced to be commented out by MG5aMC.
! Merging:nJetMax = -1
!
! For all merging schemes, decide whehter you want the merging scale
! variation computed for only the central weights or all other
! PDF and scale variation weights as well
SysCalc:
!
! =======
! User customized parameters
! =======
!
! By default, Pythia8 generates multi-parton interaction events. This is
! often irrelevant for phenomenology and very slow. You can turn this
! feature off by uncommenting the line below if so desired.
!partonlevel:mpi = off
!
! Additional technical parameters set by MG5_aMC.
!
! Tell Pythia8 that an LHEF input is used.
Beams:frameType=4
! Specify one must read inputs from the MadGraph banner.
JetMatching:
JetMatching:
! 1.0 corresponds to HEPMC weight given in [mb]. We choose here the [pb] normalization.
HEPMCoutput:
! Be more forgiving with momentum mismatches.
Check:epTolErr=
!
! =======
! Subrun definitions
! =======
!
LHEFInputs:
Main:subrun=0
!
! Definition of subrun 0
!
Beams:LHEF=
]]>
</MGPythiaCard>
<MGDelphesCard>
<![CDATA[
#######
# Order of execution of various modules
#######
set ExecutionPath {
ParticlePropa
ChargedHadron
ElectronTrack
MuonTrackingE
ChargedHadron
ElectronMomen
MuonMomentumS
TrackMerger
ECal
HCal
Calorimeter
EFlowMerger
EFlowFilter
PhotonEfficiency
PhotonIsolation
ElectronFilter
ElectronEffic
ElectronIsolation
ChargedHadron
MuonEfficiency
MuonIsolation
MissingET
NeutrinoFilter
GenJetFinder
GenMissingET
FastJetFinder
FatJetFinder
JetEnergyScale
JetFlavorAsso
BTagging
TauTagging
UniqueObjectF
ScalarHT
TreeWriter
}
#######
# Propagate particles in cylinder
#######
module ParticlePropagator ParticlePropagator {
set InputArray Delphes/
set OutputArray stableParticles
set ChargedHadronOu
set ElectronOutputArray electrons
set MuonOutputArray muons
# radius of the magnetic field coverage, in m
set Radius 1.29
# half-length of the magnetic field coverage, in m
set HalfLength 3.00
# magnetic field
set Bz 3.8
}
#######
# Charged hadron tracking efficiency
#######
module Efficiency ChargedHadronTr
set InputArray ParticlePropaga
set OutputArray chargedHadrons
# add EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for charged hadrons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
}
#######
# Electron tracking efficiency
#######
module Efficiency ElectronTrackin
set InputArray ParticlePropaga
set OutputArray electrons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for electrons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
}
#######
# Muon tracking efficiency
#######
module Efficiency MuonTrackingEff
set InputArray ParticlePropaga
set OutputArray muons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for muons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
}
#######
# Momentum resolution for charged tracks
#######
module MomentumSmearing ChargedHadronMo
set InputArray ChargedHadronTr
set OutputArray chargedHadrons
# set ResolutionFormula {resolution formula as a function of eta and pt}
# resolution formula for charged hadrons
# based on arXiv:1405.6569
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
}
#######
# Momentum resolution for electrons
#######
module MomentumSmearing ElectronMomentu
set InputArray ElectronTrackin
set OutputArray electrons
# set ResolutionFormula {resolution formula as a function of eta and energy}
# resolution formula for electrons
# based on arXiv:1502.02701
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
}
#######
# Momentum resolution for muons
#######
module MomentumSmearing MuonMomentumSme
set InputArray MuonTrackingEff
set OutputArray muons
# set ResolutionFormula {resolution formula as a function of eta and pt}
# resolution formula for muons
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
}
##############
# Track merger
##############
module Merger TrackMerger {
# add InputArray InputArray
add InputArray ChargedHadronMo
add InputArray ElectronMomentu
add InputArray MuonMomentumSme
set OutputArray tracks
}
#############
# ECAL
#############
module SimpleCalorimeter ECal {
set ParticleInputArray ParticlePropaga
set TrackInputArray TrackMerger/tracks
set TowerOutputArray ecalTowers
set EFlowTrackOutpu
set EFlowTowerOutpu
set IsEcal true
set EnergyMin 0.5
set EnergySignifica
set SmearTowerCenter true
set pi [expr {acos(-1)}]
# lists of the edges of each tower in eta and phi
# each list starts with the lower edge of the first tower
# the list ends with the higher edged of the last tower
# assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
set PhiBins {}
for {set i -180} {$i <= 180} {incr i} {
add PhiBins [expr {$i * $pi/180.0}]
}
# 0.02 unit in eta up to eta = 1.5 (barrel)
for {set i -85} {$i <= 86} {incr i} {
set eta [expr {$i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
# assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
set PhiBins {}
for {set i -180} {$i <= 180} {incr i} {
add PhiBins [expr {$i * $pi/180.0}]
}
# 0.02 unit in eta up to eta = 3
for {set i 1} {$i <= 84} {incr i} {
set eta [expr { -2.958 + $i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
for {set i 1} {$i <= 84} {incr i} {
set eta [expr { 1.4964 + $i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
# take present CMS granularity for HF
# 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
set PhiBins {}
for {set i -18} {$i <= 18} {incr i} {
add PhiBins [expr {$i * $pi/18.0}]
}
foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
add EtaPhiBins $eta $PhiBins
}
add EnergyFraction {0} {0.0}
# energy fractions for e, gamma and pi0
add EnergyFraction {11} {1.0}
add EnergyFraction {22} {1.0}
add EnergyFraction {111} {1.0}
# energy fractions for muon, neutrinos and neutralinos
add EnergyFraction {12} {0.0}
add EnergyFraction {13} {0.0}
add EnergyFraction {14} {0.0}
add EnergyFraction {16} {0.0}
add EnergyFraction {1000022} {0.0}
add EnergyFraction {1000023} {0.0}
add EnergyFraction {1000025} {0.0}
add EnergyFraction {1000035} {0.0}
add EnergyFraction {1000045} {0.0}
# energy fractions for K0short and Lambda
add EnergyFraction {310} {0.3}
add EnergyFraction {3122} {0.3}
# set ResolutionFormula {resolution formula as a function of eta and energy}
# for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
# set ECalResolutionF
# Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^
}
#############
# HCAL
#############
module SimpleCalorimeter HCal {
set ParticleInputArray ParticlePropaga
set TrackInputArray ECal/eflowTracks
set TowerOutputArray hcalTowers
set EFlowTrackOutpu
set EFlowTowerOutpu
set IsEcal false
set EnergyMin 1.0
set EnergySignifica
set SmearTowerCenter true
set pi [expr {acos(-1)}]
# lists of the edges of each tower in eta and phi
# each list starts with the lower edge of the first tower
# the list ends with the higher edged of the last tower
# 5 degrees towers
set PhiBins {}
for {set i -36} {$i <= 36} {incr i} {
add PhiBins [expr {$i * $pi/36.0}]
}
foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
add EtaPhiBins $eta $PhiBins
}
# 10 degrees towers
set PhiBins {}
for {set i -18} {$i <= 18} {incr i} {
add PhiBins [expr {$i * $pi/18.0}]
}
foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
add EtaPhiBins $eta $PhiBins
}
# 20 degrees towers
set PhiBins {}
for {set i -9} {$i <= 9} {incr i} {
add PhiBins [expr {$i * $pi/9.0}]
}
foreach eta {-5 -4.7 -4.525 4.7 5} {
add EtaPhiBins $eta $PhiBins
}
# default energy fractions {abs(PDG code)} {Fecal Fhcal}
add EnergyFraction {0} {1.0}
# energy fractions for e, gamma and pi0
add EnergyFraction {11} {0.0}
add EnergyFraction {22} {0.0}
add EnergyFraction {111} {0.0}
# energy fractions for muon, neutrinos and neutralinos
add EnergyFraction {12} {0.0}
add EnergyFraction {13} {0.0}
add EnergyFraction {14} {0.0}
add EnergyFraction {16} {0.0}
add EnergyFraction {1000022} {0.0}
add EnergyFraction {1000023} {0.0}
add EnergyFraction {1000025} {0.0}
add EnergyFraction {1000035} {0.0}
add EnergyFraction {1000045} {0.0}
# energy fractions for K0short and Lambda
add EnergyFraction {310} {0.7}
add EnergyFraction {3122} {0.7}
# set HCalResolutionF
set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^
}
#################
# Electron filter
#################
module PdgCodeFilter ElectronFilter {
set InputArray HCal/eflowTracks
set OutputArray electrons
set Invert true
add PdgCode {11}
add PdgCode {-11}
}
#######
# ChargedHadronFilter
#######
module PdgCodeFilter ChargedHadronFilter {
set InputArray HCal/eflowTracks
set OutputArray chargedHadrons
add PdgCode {11}
add PdgCode {-11}
add PdgCode {13}
add PdgCode {-13}
}
#######
# Tower Merger (in case not using e-flow algorithm)
#######
module Merger Calorimeter {
# add InputArray InputArray
add InputArray ECal/ecalTowers
add InputArray HCal/hcalTowers
set OutputArray towers
}
#######
# Energy flow merger
#######
module Merger EFlowMerger {
# add InputArray InputArray
add InputArray HCal/eflowTracks
add InputArray ECal/eflowPhotons
add InputArray HCal/eflowNeutr
set OutputArray eflow
}
#######
# EFlowFilter
#######
module PdgCodeFilter EFlowFilter {
set InputArray EFlowMerger/eflow
set OutputArray eflow
add PdgCode {11}
add PdgCode {-11}
add PdgCode {13}
add PdgCode {-13}
}
###################
# Photon efficiency
###################
module Efficiency PhotonEfficiency {
set InputArray ECal/eflowPhotons
set OutputArray photons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# efficiency formula for photons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
}
##################
# Photon isolation
##################
module Isolation PhotonIsolation {
set CandidateInputArray PhotonEfficienc
set IsolationInputArray EFlowFilter/eflow
set OutputArray photons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.12
}
#######
# Electron efficiency
#######
module Efficiency ElectronEfficiency {
set InputArray ElectronFilter/
set OutputArray electrons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# efficiency formula for electrons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
}
#######
# Electron isolation
#######
module Isolation ElectronIsolation {
set CandidateInputArray ElectronEfficie
set IsolationInputArray EFlowFilter/eflow
set OutputArray electrons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.12
}
#################
# Muon efficiency
#################
module Efficiency MuonEfficiency {
set InputArray MuonMomentumSme
set OutputArray muons
# set EfficiencyFormula {efficiency as a function of eta and pt}
# efficiency formula for muons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
}
################
# Muon isolation
################
module Isolation MuonIsolation {
set CandidateInputArray MuonEfficiency/
set IsolationInputArray EFlowFilter/eflow
set OutputArray muons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.25
}
###################
# Missing ET merger
###################
module Merger MissingET {
# add InputArray InputArray
add InputArray EFlowMerger/eflow
set MomentumOutputArray momentum
}
##################
# Scalar HT merger
##################
module Merger ScalarHT {
# add InputArray InputArray
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
set EnergyOutputArray energy
}
#######
# Neutrino Filter
#######
module PdgCodeFilter NeutrinoFilter {
set InputArray Delphes/
set OutputArray filteredParticles
set PTMin 0.0
add PdgCode {12}
add PdgCode {14}
add PdgCode {16}
add PdgCode {-12}
add PdgCode {-14}
add PdgCode {-16}
}
#######
# MC truth jet finder
#######
module FastJetFinder GenJetFinder {
set InputArray NeutrinoFilter/
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.5
set JetPTMin 20.0
}
#######
# Gen Missing ET merger
#######
module Merger GenMissingET {
# add InputArray InputArray
add InputArray NeutrinoFilter/
set MomentumOutputArray momentum
}
############
# Jet finder
############
module FastJetFinder FastJetFinder {
# set InputArray Calorimeter/towers
set InputArray EFlowMerger/eflow
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.5
set JetPTMin 20.0
}
##################
# Fat Jet finder
##################
module FastJetFinder FatJetFinder {
set InputArray EFlowMerger/eflow
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.8
set ComputeNsubjett
set Beta 1.0
set AxisMode 4
set ComputeTrimming 1
set RTrim 0.2
set PtFracTrim 0.05
set ComputePruning 1
set ZcutPrun 0.1
set RcutPrun 0.5
set RPrun 0.8
set ComputeSoftDrop 1
set BetaSoftDrop 0.0
set SymmetryCutSoftDrop 0.1
set R0SoftDrop 0.8
set JetPTMin 200.0
}
##################
# Jet Energy Scale
##################
module EnergyScale JetEnergyScale {
set InputArray FastJetFinder/jets
set OutputArray jets
# scale formula for jets
set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
}
#######
# Jet Flavor Association
#######
module JetFlavorAssoci
set PartonInputArray Delphes/partons
set ParticleInputArray Delphes/
set ParticleLHEFInp
set JetInputArray JetEnergyScale/jets
set DeltaR 0.5
set PartonPTMin 1.0
set PartonEtaMax 2.5
}
###########
# b-tagging
###########
module BTagging BTagging {
set JetInputArray JetEnergyScale/jets
set BitNumber 0
# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
# PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
# gluon's PDG code has the lowest priority
# based on arXiv:1211.4462
# default efficiency formula (misidentification rate)
add EfficiencyFormula {0} {0.01+0.000038*pt}
# efficiency formula for c-jets (misidentification rate)
add EfficiencyFormula {4} {0.25*tanh(
# efficiency formula for b-jets
add EfficiencyFormula {5} {0.85*tanh(
}
#############
# tau-tagging
#############
module TauTagging TauTagging {
set ParticleInputArray Delphes/
set PartonInputArray Delphes/partons
set JetInputArray JetEnergyScale/jets
set DeltaR 0.5
set TauPTMin 1.0
set TauEtaMax 2.5
# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
# default efficiency formula (misidentification rate)
add EfficiencyFormula {0} {0.01}
# efficiency formula for tau-jets
add EfficiencyFormula {15} {0.6}
}
#######
# Find uniquely identified photons/
#######
module UniqueObjectFinder UniqueObjectFinder {
# earlier arrays take precedence over later ones
# add InputArray InputArray OutputArray
add InputArray PhotonIsolation
add InputArray ElectronIsolati
add InputArray MuonIsolation/muons muons
add InputArray JetEnergyScale/jets jets
}
##################
# ROOT tree writer
##################
# tracks, towers and eflow objects are not stored by default in the output.
# if needed (for jet constituent or other studies), uncomment the relevant
# "add Branch ..." lines.
module TreeWriter TreeWriter {
# add Branch InputArray BranchName BranchClass
add Branch Delphes/
add Branch TrackMerger/tracks Track Track
add Branch Calorimeter/towers Tower Tower
add Branch HCal/eflowTracks EFlowTrack Track
add Branch ECal/eflowPhotons EFlowPhoton Tower
add Branch HCal/eflowNeutr
add Branch GenJetFinder/jets GenJet Jet
add Branch GenMissingET/
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch FatJetFinder/jets FatJet Jet
add Branch MissingET/momentum MissingET MissingET
add Branch ScalarHT/energy ScalarHT ScalarHT
}
]]>
</MGDelphesCard>
<MGGenerationInfo>
# Number of Events : 10000
# Integrated weight (pb) : 0.497115
</MGGenerationInfo>
</header>
<init>
11 2212 6.000000e+01 6.500000e+03 0 0 247000 247000 -4 1
4.971150e-01 1.497799e-03 4.971150e-01 1
<generator name='MadGraph5
</init>
</LesHouchesEvents>
Revision history for this message
|
#5 |
And here also is for sm-no_b_mass model:
<LesHouchesEvents version="3.0">
<header>
<!--
#******
# *
# MadGraph5_aMC@NLO *
# *
# Going Beyond *
# *
# http://
# http://
# http://
# *
# The MadGraph5_aMC@NLO team *
# *
#......
# *
# This file contains all the information necessary to reproduce *
# the events generated: *
# *
# 1. software version *
# 2. proc_card : code generation info including model *
# 3. param_card : model primary parameters in the LH format *
# 4. run_card : running parameters (collider and cuts) *
# 5. pythia_card : present only if pythia has been run *
# 6. pgs_card : present only if pgs has been run *
# 7. delphes_cards : present only if delphes has been run *
# *
# *
#******
-->
<MGVersion>
2.7.3
</MGVersion>
<MG5ProcCard>
<![CDATA[
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.7.3 2020-06-21 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set group_subprocesses Auto
set ignore_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model sm-no_b_mass
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
define p = p b b~
generate e- p > vl j b b~
output bkg1_ep2vjbb_
]]>
</MG5ProcCard>
<MGProcCard>
#******
# MadGraph/MadEvent *
# http://
# *
# proc_card.dat *
#******
# *
# This Files is generated by MADGRAPH 5 *
# *
# WARNING: This Files is generated for MADEVENT (compatibility issue)*
# This files is NOT a valid MG4 proc_card.dat *
# Running this in MG4 will NEVER reproduce the result of MG5*
# *
#******
#******
# Process(es) requested : mg2 input *
#******
# Begin PROCESS # This is TAG. Do not modify this line
e- p > vl j b b~ #Process
# Be carefull the coupling are here in MG5 convention
end_coup # End the couplings input
done # this tells MG there are no more procs
# End PROCESS # This is TAG. Do not modify this line
#******
# Model information *
#******
# Begin MODEL # This is TAG. Do not modify this line
sm-no_b_mass
# End MODEL # This is TAG. Do not modify this line
#******
# Start multiparticle definitions *
#******
# Begin MULTIPARTICLES # This is TAG. Do not modify this line
# End MULTIPARTICLES # This is TAG. Do not modify this line
</MGProcCard>
<MGRunCard>
<![CDATA[
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update full_run_card *
#******
#
#******
# Tag name for the run (one word) *
#******
tag_1 = run_tag ! name of the run
#******
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#******
10000 = nevents ! Number of unweighted events requested
30 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#******
0 = lpp1 ! beam 1 type
1 = lpp2 ! beam 2 type
60.0 = ebeam1 ! beam 1 total energy in GeV
7000.0 = ebeam2 ! beam 2 total energy in GeV
#******
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#******
0.0 = polbeam1 ! beam polarization for beam 1
0.0 = polbeam2 ! beam polarization for beam 2
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
nn23lo1 = pdlabel ! PDF set
230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#******
# Renormalization and factorization scales *
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.188 = scale ! fixed ren scale
91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
-1 = dynamical_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format
#******
False = gridpack !True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"
#******
#
#******
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#******
0 = nhel ! using helicities importance sampling or not.
#******
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.
#******
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#******
# Parton level cuts definition *
#******
#
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# Standard Cuts *
#******
# 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
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 = missetmax ! maximum missing Et (sum of neutrino's momenta)
{} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
{} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#******
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
5.0 = etaj ! max rap for the jets
-1.0 = etab ! max rap for the b
0.0 = etabmin ! min rap for the b
{} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
{} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#******
# Minimum and maximum DeltaR distance *
#******
0.0 = drbb ! min distance between b's
0.0 = drbj ! min distance between b and jet
-1.0 = drbbmax ! max distance between b's
-1.0 = drbjmax ! max distance between b and jet
#******
# Minimum and maximum invariant mass for pairs *
#******
0.0 = mmbb ! min invariant mass of a b pair
-1.0 = mmbbmax ! max invariant mass of a b pair
{} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
{'default': False} = mxx_only_
#******
# Inclusive cuts *
#******
0.0 = xptb ! minimum pt for at least one b
#*****
# Control the Ht(k)=Sum of k leading jets *
#*****
0.0 = ihtmin !inclusive Ht for all partons (including b)
-1.0 = ihtmax !inclusive Ht for all partons (including b)
0.0 = ht2min ! minimum Ht for the two leading jets
0.0 = ht3min ! minimum Ht for the three leading jets
-1.0 = ht2max ! maximum Ht for the two leading jets
-1.0 = ht3max ! maximum Ht for the three leading jets
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
#
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
True = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_
# Syscalc is deprecated but to see the associate options type'update syscalc'
]]>
</MGRunCard>
<slha>
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
#######
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
#######
#######
## INFORMATION FOR MASS
#######
Block mass
6 1.730000e+02 # MT
15 1.777000e+00 # MTA
23 9.118800e+01 # MZ
25 1.250000e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
1 0.000000e+00 # d : 0.0
2 0.000000e+00 # u : 0.0
3 0.000000e+00 # s : 0.0
4 0.000000e+00 # c : 0.0
5 0.000000e+00 # b : 0.0
11 0.000000e+00 # e- : 0.0
12 0.000000e+00 # ve : 0.0
13 0.000000e+00 # mu- : 0.0
14 0.000000e+00 # vm : 0.0
16 0.000000e+00 # vt : 0.0
21 0.000000e+00 # g : 0.0
22 0.000000e+00 # a : 0.0
24 8.041900e+01 # w+ : cmath.sqrt(
#######
## INFORMATION FOR SMINPUTS
#######
Block sminputs
1 1.325070e+02 # aEWM1
2 1.166390e-05 # Gf
3 1.180000e-01 # aS
#######
## INFORMATION FOR YUKAWA
#######
Block yukawa
6 1.730000e+02 # ymt
15 1.777000e+00 # ymtau
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
</slha>
<MGPythiaCard>
<![CDATA[
!
! It is possible to run this card manually with:
! LD_LIBRARY_
!
!
! Pythia8 cmd card automatically generated by MadGraph5_aMC@NLO
! For more information on the use of the MG5aMC / Pythia8 interface, visit
! https:/
!
! ==================
! General parameters
! ==================
!
Main:numberOfEvents = -1
!
! -------
! Specify the HEPMC output of the Pythia8 shower. You can set it to:
! auto : MG5aMC will automatically place it the run_<i> directory
! /dev/null : to turn off the HEPMC output.
! <path> : to select where the HEPMC file must written. It will
! therefore not be placed in the run_<i> directory. The
! specified path, if not absolute, will be relative to
! the Event/run_<i> directory of the process output.
! fifo : to have MG5aMC setup the piping of the PY8 output to
! analysis tools such as MadAnalysis5.
! fifo@<fifo_path> :
! Same as 'fifo', but selecting a custom path to create the
! fifo pipe. (useful to select a mounted drive that supports
! fifo). Note that the fifo file extension *must* be '.hepmc.fifo'.
! -------
!
HEPMCoutput:file = tag_1_pythia8_
!
! -------
! Parameters relevant only when performing MLM merging, which can be
! turned on by setting ickkw to '1' in the run_card and chosing a
! positive value for the parameter xqcut.
! For details, see section 'Jet Matching' on the left-hand menu of
! http://
! -------
! If equal to -1.0, MadGraph5_aMC@NLO will set it automatically based
! on the parameter 'xqcut' of the run_card.dat
! The following parameter was forced to be commented out by MG5aMC.
! JetMatching:qCut = -1.0
! Use default kt-MLM to match parton level jets to those produced by the
! shower. But the other Shower-kt scheme is available too with this option.
! The following parameter was forced to be commented out by MG5aMC.
! JetMatching:
! A value of -1 means that it is automatically guessed by MadGraph.
! It is however always safer to explicitly set it.
! The following parameter was forced to be commented out by MG5aMC.
! JetMatching:nJetMax = -1
!
! -------
! Parameters relevant only when performing CKKW-L merging, which can
! be turned on by setting the parameter 'ptlund' *or* 'ktdurham' to
! a positive value.
! For details, see section 'CKKW-L Merging' on the left-hand menu of
! http://
! -------
! Central merging scale values you want to be used.
! If equal to -1.0, then MadGraph5_aMC@NLO will set this automatically
! based on the parameter 'ktdurham' of the run_card.dat
! The following parameter was forced to be commented out by MG5aMC.
! Merging:TMS = -1.0
! This must be set manually, according to Pythia8 directives.
! An example of possible value is 'pp>LEPTONS,
! Alternatively, from Pythia v8.223 onwards, the value 'guess' can be
! used to instruct Pythia to guess the hard process. The guess would mean
! that all particles apart from light partons will be considered as a part
! of the hard process. This guess is prone to errors if the desired hard
! process is complicated (i.e. contains light partons). The user should
! then be wary of suspicious error messages in the Pythia log file.
! The following parameter was forced to be commented out by MG5aMC.
! Merging:Process = <set_by_user>
! A value of -1 means that it is automatically guessed by MadGraph.
! It is however always safer to explicitly set it.
! The following parameter was forced to be commented out by MG5aMC.
! Merging:nJetMax = -1
!
! For all merging schemes, decide whehter you want the merging scale
! variation computed for only the central weights or all other
! PDF and scale variation weights as well
SysCalc:
!
! =======
! User customized parameters
! =======
!
! By default, Pythia8 generates multi-parton interaction events. This is
! often irrelevant for phenomenology and very slow. You can turn this
! feature off by uncommenting the line below if so desired.
!partonlevel:mpi = off
!
! Additional technical parameters set by MG5_aMC.
!
! Tell Pythia8 that an LHEF input is used.
Beams:frameType=4
! Specify one must read inputs from the MadGraph banner.
JetMatching:
JetMatching:
! 1.0 corresponds to HEPMC weight given in [mb]. We choose here the [pb] normalization.
HEPMCoutput:
! Be more forgiving with momentum mismatches.
Check:epTolErr=
!
! =======
! Subrun definitions
! =======
!
LHEFInputs:
Main:subrun=0
!
! Definition of subrun 0
!
Beams:LHEF=
]]>
</MGPythiaCard>
<MGDelphesCard>
<![CDATA[
#######
# Order of execution of various modules
#######
set ExecutionPath {
ParticlePropa
ChargedHadron
ElectronTrack
MuonTrackingE
ChargedHadron
ElectronMomen
MuonMomentumS
TrackMerger
ECal
HCal
Calorimeter
EFlowMerger
EFlowFilter
PhotonEfficiency
PhotonIsolation
ElectronFilter
ElectronEffic
ElectronIsolation
ChargedHadron
MuonEfficiency
MuonIsolation
MissingET
NeutrinoFilter
GenJetFinder
GenMissingET
FastJetFinder
FatJetFinder
JetEnergyScale
JetFlavorAsso
BTagging
TauTagging
UniqueObjectF
ScalarHT
TreeWriter
}
#######
# Propagate particles in cylinder
#######
module ParticlePropagator ParticlePropagator {
set InputArray Delphes/
set OutputArray stableParticles
set ChargedHadronOu
set ElectronOutputArray electrons
set MuonOutputArray muons
# radius of the magnetic field coverage, in m
set Radius 1.29
# half-length of the magnetic field coverage, in m
set HalfLength 3.00
# magnetic field
set Bz 3.8
}
#######
# Charged hadron tracking efficiency
#######
module Efficiency ChargedHadronTr
set InputArray ParticlePropaga
set OutputArray chargedHadrons
# add EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for charged hadrons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
}
#######
# Electron tracking efficiency
#######
module Efficiency ElectronTrackin
set InputArray ParticlePropaga
set OutputArray electrons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for electrons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
}
#######
# Muon tracking efficiency
#######
module Efficiency MuonTrackingEff
set InputArray ParticlePropaga
set OutputArray muons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for muons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
}
#######
# Momentum resolution for charged tracks
#######
module MomentumSmearing ChargedHadronMo
set InputArray ChargedHadronTr
set OutputArray chargedHadrons
# set ResolutionFormula {resolution formula as a function of eta and pt}
# resolution formula for charged hadrons
# based on arXiv:1405.6569
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
}
#######
# Momentum resolution for electrons
#######
module MomentumSmearing ElectronMomentu
set InputArray ElectronTrackin
set OutputArray electrons
# set ResolutionFormula {resolution formula as a function of eta and energy}
# resolution formula for electrons
# based on arXiv:1502.02701
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
}
#######
# Momentum resolution for muons
#######
module MomentumSmearing MuonMomentumSme
set InputArray MuonTrackingEff
set OutputArray muons
# set ResolutionFormula {resolution formula as a function of eta and pt}
# resolution formula for muons
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
}
##############
# Track merger
##############
module Merger TrackMerger {
# add InputArray InputArray
add InputArray ChargedHadronMo
add InputArray ElectronMomentu
add InputArray MuonMomentumSme
set OutputArray tracks
}
#############
# ECAL
#############
module SimpleCalorimeter ECal {
set ParticleInputArray ParticlePropaga
set TrackInputArray TrackMerger/tracks
set TowerOutputArray ecalTowers
set EFlowTrackOutpu
set EFlowTowerOutpu
set IsEcal true
set EnergyMin 0.5
set EnergySignifica
set SmearTowerCenter true
set pi [expr {acos(-1)}]
# lists of the edges of each tower in eta and phi
# each list starts with the lower edge of the first tower
# the list ends with the higher edged of the last tower
# assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
set PhiBins {}
for {set i -180} {$i <= 180} {incr i} {
add PhiBins [expr {$i * $pi/180.0}]
}
# 0.02 unit in eta up to eta = 1.5 (barrel)
for {set i -85} {$i <= 86} {incr i} {
set eta [expr {$i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
# assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
set PhiBins {}
for {set i -180} {$i <= 180} {incr i} {
add PhiBins [expr {$i * $pi/180.0}]
}
# 0.02 unit in eta up to eta = 3
for {set i 1} {$i <= 84} {incr i} {
set eta [expr { -2.958 + $i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
for {set i 1} {$i <= 84} {incr i} {
set eta [expr { 1.4964 + $i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
# take present CMS granularity for HF
# 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
set PhiBins {}
for {set i -18} {$i <= 18} {incr i} {
add PhiBins [expr {$i * $pi/18.0}]
}
foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
add EtaPhiBins $eta $PhiBins
}
add EnergyFraction {0} {0.0}
# energy fractions for e, gamma and pi0
add EnergyFraction {11} {1.0}
add EnergyFraction {22} {1.0}
add EnergyFraction {111} {1.0}
# energy fractions for muon, neutrinos and neutralinos
add EnergyFraction {12} {0.0}
add EnergyFraction {13} {0.0}
add EnergyFraction {14} {0.0}
add EnergyFraction {16} {0.0}
add EnergyFraction {1000022} {0.0}
add EnergyFraction {1000023} {0.0}
add EnergyFraction {1000025} {0.0}
add EnergyFraction {1000035} {0.0}
add EnergyFraction {1000045} {0.0}
# energy fractions for K0short and Lambda
add EnergyFraction {310} {0.3}
add EnergyFraction {3122} {0.3}
# set ResolutionFormula {resolution formula as a function of eta and energy}
# for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
# set ECalResolutionF
# Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^
}
#############
# HCAL
#############
module SimpleCalorimeter HCal {
set ParticleInputArray ParticlePropaga
set TrackInputArray ECal/eflowTracks
set TowerOutputArray hcalTowers
set EFlowTrackOutpu
set EFlowTowerOutpu
set IsEcal false
set EnergyMin 1.0
set EnergySignifica
set SmearTowerCenter true
set pi [expr {acos(-1)}]
# lists of the edges of each tower in eta and phi
# each list starts with the lower edge of the first tower
# the list ends with the higher edged of the last tower
# 5 degrees towers
set PhiBins {}
for {set i -36} {$i <= 36} {incr i} {
add PhiBins [expr {$i * $pi/36.0}]
}
foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
add EtaPhiBins $eta $PhiBins
}
# 10 degrees towers
set PhiBins {}
for {set i -18} {$i <= 18} {incr i} {
add PhiBins [expr {$i * $pi/18.0}]
}
foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
add EtaPhiBins $eta $PhiBins
}
# 20 degrees towers
set PhiBins {}
for {set i -9} {$i <= 9} {incr i} {
add PhiBins [expr {$i * $pi/9.0}]
}
foreach eta {-5 -4.7 -4.525 4.7 5} {
add EtaPhiBins $eta $PhiBins
}
# default energy fractions {abs(PDG code)} {Fecal Fhcal}
add EnergyFraction {0} {1.0}
# energy fractions for e, gamma and pi0
add EnergyFraction {11} {0.0}
add EnergyFraction {22} {0.0}
add EnergyFraction {111} {0.0}
# energy fractions for muon, neutrinos and neutralinos
add EnergyFraction {12} {0.0}
add EnergyFraction {13} {0.0}
add EnergyFraction {14} {0.0}
add EnergyFraction {16} {0.0}
add EnergyFraction {1000022} {0.0}
add EnergyFraction {1000023} {0.0}
add EnergyFraction {1000025} {0.0}
add EnergyFraction {1000035} {0.0}
add EnergyFraction {1000045} {0.0}
# energy fractions for K0short and Lambda
add EnergyFraction {310} {0.7}
add EnergyFraction {3122} {0.7}
# set HCalResolutionF
set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^
}
#################
# Electron filter
#################
module PdgCodeFilter ElectronFilter {
set InputArray HCal/eflowTracks
set OutputArray electrons
set Invert true
add PdgCode {11}
add PdgCode {-11}
}
#######
# ChargedHadronFilter
#######
module PdgCodeFilter ChargedHadronFilter {
set InputArray HCal/eflowTracks
set OutputArray chargedHadrons
add PdgCode {11}
add PdgCode {-11}
add PdgCode {13}
add PdgCode {-13}
}
#######
# Tower Merger (in case not using e-flow algorithm)
#######
module Merger Calorimeter {
# add InputArray InputArray
add InputArray ECal/ecalTowers
add InputArray HCal/hcalTowers
set OutputArray towers
}
#######
# Energy flow merger
#######
module Merger EFlowMerger {
# add InputArray InputArray
add InputArray HCal/eflowTracks
add InputArray ECal/eflowPhotons
add InputArray HCal/eflowNeutr
set OutputArray eflow
}
#######
# EFlowFilter
#######
module PdgCodeFilter EFlowFilter {
set InputArray EFlowMerger/eflow
set OutputArray eflow
add PdgCode {11}
add PdgCode {-11}
add PdgCode {13}
add PdgCode {-13}
}
###################
# Photon efficiency
###################
module Efficiency PhotonEfficiency {
set InputArray ECal/eflowPhotons
set OutputArray photons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# efficiency formula for photons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
}
##################
# Photon isolation
##################
module Isolation PhotonIsolation {
set CandidateInputArray PhotonEfficienc
set IsolationInputArray EFlowFilter/eflow
set OutputArray photons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.12
}
#######
# Electron efficiency
#######
module Efficiency ElectronEfficiency {
set InputArray ElectronFilter/
set OutputArray electrons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# efficiency formula for electrons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
}
#######
# Electron isolation
#######
module Isolation ElectronIsolation {
set CandidateInputArray ElectronEfficie
set IsolationInputArray EFlowFilter/eflow
set OutputArray electrons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.12
}
#################
# Muon efficiency
#################
module Efficiency MuonEfficiency {
set InputArray MuonMomentumSme
set OutputArray muons
# set EfficiencyFormula {efficiency as a function of eta and pt}
# efficiency formula for muons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
}
################
# Muon isolation
################
module Isolation MuonIsolation {
set CandidateInputArray MuonEfficiency/
set IsolationInputArray EFlowFilter/eflow
set OutputArray muons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.25
}
###################
# Missing ET merger
###################
module Merger MissingET {
# add InputArray InputArray
add InputArray EFlowMerger/eflow
set MomentumOutputArray momentum
}
##################
# Scalar HT merger
##################
module Merger ScalarHT {
# add InputArray InputArray
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
set EnergyOutputArray energy
}
#######
# Neutrino Filter
#######
module PdgCodeFilter NeutrinoFilter {
set InputArray Delphes/
set OutputArray filteredParticles
set PTMin 0.0
add PdgCode {12}
add PdgCode {14}
add PdgCode {16}
add PdgCode {-12}
add PdgCode {-14}
add PdgCode {-16}
}
#######
# MC truth jet finder
#######
module FastJetFinder GenJetFinder {
set InputArray NeutrinoFilter/
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.5
set JetPTMin 20.0
}
#######
# Gen Missing ET merger
#######
module Merger GenMissingET {
# add InputArray InputArray
add InputArray NeutrinoFilter/
set MomentumOutputArray momentum
}
############
# Jet finder
############
module FastJetFinder FastJetFinder {
# set InputArray Calorimeter/towers
set InputArray EFlowMerger/eflow
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.5
set JetPTMin 20.0
}
##################
# Fat Jet finder
##################
module FastJetFinder FatJetFinder {
set InputArray EFlowMerger/eflow
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.8
set ComputeNsubjett
set Beta 1.0
set AxisMode 4
set ComputeTrimming 1
set RTrim 0.2
set PtFracTrim 0.05
set ComputePruning 1
set ZcutPrun 0.1
set RcutPrun 0.5
set RPrun 0.8
set ComputeSoftDrop 1
set BetaSoftDrop 0.0
set SymmetryCutSoftDrop 0.1
set R0SoftDrop 0.8
set JetPTMin 200.0
}
##################
# Jet Energy Scale
##################
module EnergyScale JetEnergyScale {
set InputArray FastJetFinder/jets
set OutputArray jets
# scale formula for jets
set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
}
#######
# Jet Flavor Association
#######
module JetFlavorAssoci
set PartonInputArray Delphes/partons
set ParticleInputArray Delphes/
set ParticleLHEFInp
set JetInputArray JetEnergyScale/jets
set DeltaR 0.5
set PartonPTMin 1.0
set PartonEtaMax 2.5
}
###########
# b-tagging
###########
module BTagging BTagging {
set JetInputArray JetEnergyScale/jets
set BitNumber 0
# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
# PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
# gluon's PDG code has the lowest priority
# based on arXiv:1211.4462
# default efficiency formula (misidentification rate)
add EfficiencyFormula {0} {0.01+0.000038*pt}
# efficiency formula for c-jets (misidentification rate)
add EfficiencyFormula {4} {0.25*tanh(
# efficiency formula for b-jets
add EfficiencyFormula {5} {0.85*tanh(
}
#############
# tau-tagging
#############
module TauTagging TauTagging {
set ParticleInputArray Delphes/
set PartonInputArray Delphes/partons
set JetInputArray JetEnergyScale/jets
set DeltaR 0.5
set TauPTMin 1.0
set TauEtaMax 2.5
# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
# default efficiency formula (misidentification rate)
add EfficiencyFormula {0} {0.01}
# efficiency formula for tau-jets
add EfficiencyFormula {15} {0.6}
}
#######
# Find uniquely identified photons/
#######
module UniqueObjectFinder UniqueObjectFinder {
# earlier arrays take precedence over later ones
# add InputArray InputArray OutputArray
add InputArray PhotonIsolation
add InputArray ElectronIsolati
add InputArray MuonIsolation/muons muons
add InputArray JetEnergyScale/jets jets
}
##################
# ROOT tree writer
##################
# tracks, towers and eflow objects are not stored by default in the output.
# if needed (for jet constituent or other studies), uncomment the relevant
# "add Branch ..." lines.
module TreeWriter TreeWriter {
# add Branch InputArray BranchName BranchClass
add Branch Delphes/
add Branch TrackMerger/tracks Track Track
add Branch Calorimeter/towers Tower Tower
add Branch HCal/eflowTracks EFlowTrack Track
add Branch ECal/eflowPhotons EFlowPhoton Tower
add Branch HCal/eflowNeutr
add Branch GenJetFinder/jets GenJet Jet
add Branch GenMissingET/
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch FatJetFinder/jets FatJet Jet
add Branch MissingET/momentum MissingET MissingET
add Branch ScalarHT/energy ScalarHT ScalarHT
}
]]>
</MGDelphesCard>
<MGGenerationInfo>
# Number of Events : 10000
# Integrated weight (pb) : 272.419
</MGGenerationInfo>
</header>
<init>
11 2212 6.000000e+01 7.000000e+03 0 0 247000 247000 -4 1
2.724190e+02 8.190691e+01 2.724190e+02 1
<generator name='MadGraph5
</init>
</LesHouchesEvents>
Revision history for this message
|
#6 |
Hi,
The issue with the sm-no_b_mass is that you do not have any cut on the b quark.
Looks like the maxjetflavour parameter is still 4 in this case. and therefore the "b" quark have their own set of cuts different of the ptj cut (and equivalent for other variable. Consequently you have a singular regime, the cross-section is technically not finite and it is not possible to generate events.
I would see what I can do to change the default value of maxjetflavour for such type of process and such model.
Concerning pythia8, I do not see any flag to have assymetric handling of the beam. I have never tried to use pythia8 for dIS type of collision so I have never learn how to configure pythia8 for such type of collision. I would advise you to read the pythia8 manual to learn that.
Cheers,
Olivier
Revision history for this message
|
#7 |
Dear Olivier,
Thanks a lot for your time and checking the banner files.
I am looking forward to hearing from you about the first issue and I will try to solve the pythia8 problem.
Many thanks,
Reza
Revision history for this message
|
#8 |
Dear Olivier,
Before you would see what you can do to change the default value of maxjetflavour for such type of process and such model,
would you please let me I should change maxjetflavour to 5 because in sm-no_b_mass model, b quark is a light quark ?
Revision history for this message
|
#9 |
Hi,
You have to put cut on the b quark.
One easy way is to set maxjetflavour to 5 such that all default cut for light jet are also applied on b quark.
The second option is to keep maxjetflavour to 4 and add one cut preventing soft singularities and another one preventing collinear singularities.
Cheers,
Olivier
Revision history for this message
|
#10 |
Dear Olivier,
Thanks a lot for your help. Both options solved the issue at the MadGraph level.
I am searching about Pythia8 problem now.
Regards,
Reza
Revision history for this message
|
#11 |
Thanks Olivier Mattelaer, that solved my question.
Revision history for this message
|
#12 |
Dear Olivier,
As the problem solved in MadGraph level, now I could not handle Pythia8 and I think you are more familiar with Pythia6, I tried also with Pythia6.
The problem still exists. I would be appreciated if you take a quick look in pythia_card and below part of tag_1_pythia.log file:
Thanks,
Reza
--
some parts of tag_1_pythia.log:
Initiating pythia with external process
1****************** PYINIT: initialization of PYTHIA routines *****************
Reading parameter sminputs 1 132.50700000000001
Reading parameter sminputs 2 1.1663900000000
Reading parameter sminputs 3 0.11799999999999999
Reading model: sm-no_b_
IMSS(21) changed from 0 to 24
Initializing PYR with random seed 102
MSTP(52) changed from 1 to 2
Warning: parameter ickkw not found
setting it to default value 0
MSTP(98) changed from 0 to 1
==== PYTHIA WILL USE LHAPDF ====
WRONG LHAPDF set number = 247000 given! STOP EXE!
WILL USE 10041 INSTEAD
******
* LHAPDF Version 4.0 *
******
>>>>>> PDF description: <<<<<<
CTEQ6L - LO with NLO alpha_s
Reference:
J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky,
STDXWOPEN WARNING: I/O is initialized for stdhep only
W.K. Tung
hep-ph/0201195
>>>>>> <<<<<<
Parametrization: CTEQ6
warning empty CtLhQCDSET call: NAME=IHDN VALUE= 1.0000000000000000
warning empty CtLhQCDSET call: NAME=QMAX VALUE= 100000.00000000000
warning empty CtLhQCDSET call: NAME=XMIN VALUE= 9.9999999999999
warning empty CtLhQCDSET call: NAME=QINI VALUE= 1.3000000000000000
warning empty CtLhQCDSET call: NAME=IKNL VALUE= 1.0000000000000000
======
PDFset name /data/rjafari/
with 1 members
==== initialized. =======
warning empty CtLhQCDSET call: NAME=NX VALUE= 95.000000000000000
warning empty CtLhQCDSET call: NAME=NT VALUE= 18.000000000000000
Strong coupling at Mz for PDF is: 0.11798
...
======
I I I I
I Subprocess I Number of points I Sigma I
I I I I
I-----
I I I I
I N:o Type I Generated Tried I I
I I I I
======
I I I I
I 0 All included subprocesses I 367 10000 I 1.392D-12 I
I 4 User process 1 I 367 10000 I 1.392D-12 I
I I I I
======
********* Total number of errors, excluding junctions = 924 *************
********* Total number of errors, including junctions = 927 *************
********* Total number of warnings = 2 *************
********* Fraction of events that fail fragmentation cuts = 0.96330 *********
--
pythia_card:
!...Parton showering on or off
MSTP(61)=1
MSTP(71)=1
!...Fragmentati
MSTJ(1)=1
!...Multiple interactions on or off
MSTP(81)=20
!...Don't stop execution after 10 errors
MSTU(21)=1
!...PDFset if MG set not supported by pythia-pgs package (set in lhapdf5 or higher)
LHAID= 10041
LHAPATH=
Revision history for this message
|
#13 |
Hi,
I did not run pythia6 since ~10 years so I actually know more about pythia8 than pythia6.
Pythia6 should not be used anymore, it is not supported anymore by the pythia author. So valid reason to use pythia6 today are very limited (reproduction of legacy result or tracking of very old bug).
You should really need to read the pythia8 manual to learn how to setup flag for such type of process.
If something is not clear in their manual, you should then contact the autors of the tool.
Cheers,
Olivier
Revision history for this message
|
#14 |
Dear Olivier,
Thanks a lot for your prompt and clear explanation.
I will read the manual carefully and hope that solves my problem.
Regards,
Reza
Can you help with this problem?
Provide an answer of your own, or ask Reza Jafari for more information if necessary.