Closed mafrahm closed 3 months ago
Really nice, and a lot of work!
I'll approve but I would wait for @nprouvost to agree as well.
Just a side note: some of the blocks in which you define higgs subprocesses look fairly similar and just a detail like the name of the decay map changes. Could there be a way to also loop over that? It might not be particularly clean, but to add the subprocess with their correct names to the module, one could do
locals()[the_name] = add_decay_process(...)
.
Yes, it's always the same block adding all the Higgs decay channels and their subsequent decay channels. I propose we do this in a separate PR after this.
This process adds all combinations of single Higgs production and decay channel processes (around 450 processes in total). To accomplish this, I implemented helper functions (e.g.
add_decay_process
) that allows creating sub-processes with a single line of code. This makes our process highly generalized and allows us to easily accomplish changes in the naming, IDs, and labels of lots of processes by changing a single value in a dictionary.Open questions:
H_{ggf}, H \rightarrow bb
for Process h_ggf_hbb. We could skip the, H
part of the label, but I'm not sure how to generalize this)