genproductions icon indicating copy to clipboard operation
genproductions copied to clipboard

GammaGamma to tau tau inclusive cards

Open michael-pitt opened this issue 2 years ago • 8 comments

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

michael-pitt avatar Nov 29 '22 15:11 michael-pitt

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

menglu21 avatar Dec 07 '22 05:12 menglu21

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

michael-pitt avatar Jan 15 '23 04:01 michael-pitt

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/

michael-pitt avatar Jan 15 '23 12:01 michael-pitt

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

menglu21 avatar Feb 01 '23 08:02 menglu21

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

sihyunjeon avatar Jun 30 '23 11:06 sihyunjeon

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

michael-pitt avatar Aug 15 '23 13:08 michael-pitt

We keep this one open, until we manage to get MG3.5 available soon.

bbilin avatar Mar 06 '24 13:03 bbilin

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?

Saptaparna avatar Jun 25 '24 12:06 Saptaparna