genproductions
genproductions copied to clipboard
GammaGamma to tau tau inclusive cards
Dear expert,
I've produced gammagamma->tautau inclusive production madgraph cards at 13TeV, which included BSM weights, and are used for g-2 measurement of the tau lepton.
For inelastic photons. the MSHT20qed_nnlo_inelastic set is used:
lhapdf = pdlabel ! PDF set
26500 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
For elastic photons hadronic survival probabilities with charge forma factors were used based on gamma-UPC model, and implemented using:
2 = lpp1 ! beam 1 type - photon from proton
chff = pdlabel ! charge form factors
Please review the cards, and if there are any issues, please let me know.
Thank you in advance! Michael
Hi @michael-pitt another minor question is that have you had any expectation on the pdf uncertainty, since the current pdf list of UL doesn't include many qed pdf
Hi @menglu21
I've changed the cards. Previously we were using pdf-based photon fluxes, which are incorrect if one uses proton tagging or rapidity gaps to identify photon-induced processes.
A recent MadGraph implementation of the correction factors was introduced in 2207.03012 and will become available in MG3.5.0 for elastic case only (photons from protons using lpp=2
).
I prepared in advance cards in 4 folders: Elastic tautau production with a mass cut of 50 or 300 GeV (two folders) Elastic-Inelastic and Inelastic-Elastic with mtautau>50 GeV cut (two folders)
The Elastic-Inelastic need to be debugged (fails during generate_events pilotrun), but the first two cards are now used in the analysis.
The cards can run with older MadGraph releases if one switches off the corrections by changing pdlabel
from edff
to none
in the run_card.dat.
One can assign systematic uncertainty on the photon flux by changing pdlabel between edff
and chff
, I will cross-check with authors if this is the right way to do it.
Regards, Michael
This is a followup on the error with Inelastic cards. I uploaded a card for fully inelastic gammagamma->tautau production (both photons taken from pdf). Unfortunately, the following command works only with mg265UL
but not with the latest branches (tested with mg27x
and mg33x
), complaining about LHAPDF:
./gridpack_generation.sh GGToTauTau_Inelastic_M_50GeV_TuneCP5_madgraphLO_reweight cards/production/13TeV/GammaGammaToTauTau/GGToTauTau_Inelastic_M_50GeV_TuneCP5_madgraphLO_reweight/
Hi @michael-pitt so currently only fully inelastic cards have issue? what's the error, I can't test since i don't have the model
you can't have tau > all all
in madspin since SMEFT model does not treat mesons as a stable particle. you have to let pythia do the rest (or tauola).
also drll 0.4
is a bit tight as your analysis level will likely give the same cut for the leptons. for better (or realistic) modelling, it would be better to set this to 0
This PR is currently waiting for two issues to be solved:
- addition of MG3.5 version (discussed in https://github.com/cms-sw/genproductions/issues/3326)
- adding kT smearing (in discussion with pythia experts)
Once it is done, then we can proceed with this RP.
Thanks, Michael
We keep this one open, until we manage to get MG3.5 available soon.
Since the analysts are giving a talk right now, I was wondering if the MG3.5 or MG3.X branch is ready to use and have the analysts tried updating to it?