Param and Run Card Edits
This is so much easier to see specifically what options you're running with now. So this ticket is for discussing some of the things I see at a quick glance.
-
https://gitlab.cern.ch/jkishine/chargino-proj/-/blob/652fc0b81eaf0edbaa257151ce2bf30b8986f62c/cards/run_card.dat#L33-36
- Good.
lpp1
andlpp2
are2
, so the incoming particle from both beams will be considered a photon emitted from a proton. - However, I see that both
ebeam1
andebeam2
are set to 14 TeV. This represents a 28 TeV proton collider. For the HL-LHC colliding Pb, the collision energy should be 14 TeV * (Np/(Np+Nn)) = 14 TeV * (82/208) = 5.519 TeV. But that's then split between the two beam energies, soebeam1
andebeam2
should be2760
.
- Good.
-
https://gitlab.cern.ch/jkishine/chargino-proj/-/blob/652fc0b81eaf0edbaa257151ce2bf30b8986f62c/cards/param_card.dat#L315-318
- Here, don't forget to fix the decay block to what it should be. The relevant code snippets from the notes I posted last week are:
masses['1000022'] = n1_mass
# C1 relevant
masses['1000024'] = c1_mass
c1_dict = chargino_width_GeV[(c1_mass,n1_mass)]
decays["1000024"]="""DECAY 1000024 %E # chargino1+ decays
# BR NDA ID1 ID2 ID3
0.00000000E+00 2 1000022 24 # BR(~chi_1+ -> ~chi_10 W+)
%.8E 2 1000022 211 # BR(~chi_1+ -> ~chi_10 pi+)
%.8E 3 1000022 211 111 # BR(~chi_1+ -> ~chi_10 pi+)
%.8E 3 1000022 -11 12 # BR(~chi_1+ -> ~chi_10 e+ nu_e)
%.8E 3 1000022 -13 14 # BR(~chi_1+ -> ~chi_10 mu+ nu_mu)
0.00000000E+00 3 1000022 -15 16 # BR(~chi_1+ -> ~chi_10 tau+ nu_tau)
#
"""%(c1_dict["width"], c1_dict["BR_1pi"], c1_dict["BR_2pi"], c1_dict["BR_e"], c1_dict["BR_mu"])
and
chargino_width_GeV = {
(150.35, 150): {"width": 1.5220922294648483e-14, "BR_1pi": 0.9171530337320101, "BR_2pi": 0.00001655377967743894, "BR_e": 0.050291041479660833, "BR_mu": 0.03253937100865154},
(150.5, 150): {"width": 5.089812868085344e-14, "BR_1pi": 0.8368779656992351, "BR_2pi": 0.0016986105266600206, "BR_e": 0.08934930351612039, "BR_mu": 0.07207412025798424},
(151, 150): {"width": 7.881725041405977e-13, "BR_1pi": 0.44510054134817434, "BR_2pi": 0.1972427980805653, "BR_e": 0.18372486722969916, "BR_mu": 0.17393179334156117}
}