Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Run3 bbH cards for producing bbHToGG samples #3774

Open
wants to merge 4 commits into
base: master
Choose a base branch
from

Conversation

nitishkumarkv
Copy link
Contributor

@nitishkumarkv nitishkumarkv commented Oct 8, 2024

Dear experts,

I have added the bbH cards for yb2, ybyt and yt2 contributions. The cards are a copy from Run2 but only the beam energy is changed to 6800.0 GeV. I assume decay products should be defined in the fragment file.

For the yb2 contribution, I unpacked one of the gridpacks of bbHToTauTau_yb2 samples and I found the proc_card (in folder process/Cards/proc_card_mg5.dat) matching with SUSYGluGluToBBH process card. But there are also separate cards for 2017 bbH_4FS_yb2. Please let me know which one should be used.

Please review the cards and mention any changes if needed

Comment on lines +12 to +15
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

These can be removed as they are not used

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Dear @lviliani,

Thanks a lot for the detailed comments. As I mentioned before, the cards are the exact copy of what was used for producing Run2 bbHToTauTau samples, just that I have changed the beam energy.

I will remove these lines if they do not affect anything.

Comment on lines +64 to +67
lhapdf = pdlabel ! PDF set
320900 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
! numbers for central PDF sets are allowed. Can be a list;
! PDF sets beyond the first are included via reweighting.
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is there a reason for not using the default PDF set? You can get the default automatically adding these lines:

     'lhapdf'    = pdlabel     ! PDF set                                     
     $DEFAULT_PDF_SETS    = lhaid     ! if pdlabel=lhapdf, this is the lhapdf number
     $DEFAULT_PDF_MEMBERS = reweight_PDF

Copy link
Contributor Author

@nitishkumarkv nitishkumarkv Nov 6, 2024

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Sorry for my limited knowledge on this part. Could you please explain me what would change if we keep it to default?

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

With the default settings you get the recommended PDF sets automatically.

set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do you need this line?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would keep it like that as it was used also for Run2.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This looks the same as the previous card, so you can drop it.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You can drop this card as well

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do you need this param_card? What do you need change wrt default values?

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Same comments as previous run_card

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I have kept the param_card and run_card for each contributions(yb2, yt2 and ybyt), even though they might be same. I did this so as to be consistent with what was done for Run2. Please let me know if the repeated param_card and run_card still needs to be removed.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, please remove all the duplicate cards.

@@ -0,0 +1,6 @@
import model HC_NLO_X0_UFO-heft_bmass_noyb

generate p p > x0 b b~ / t [QCD]
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

with the / t syntax you are removing all the diagrams containing top quarks. Is it intentional?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think it is intentional, anyway it does not restrict top quarks loops

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants