37 {
38
39
40 m_kCoup = settingsPtr->parm(
"LeptoQuark:kCoup");
41
42
43#if PYTHIA_VERSION_INTEGER >= 8307
45#else
46 ParticleDataEntry* particleSPtr = particlePtr;
47#endif
48
49
50 int id1Now = particleSPtr->channel(0).product(0);
51 int id2Now = particleSPtr->channel(0).product(1);
52
53
54
55
56
57 if (id1Now < 1 || id1Now > 6) {
58 std::cout << "ERROR in ResonanceLQ::init: unallowed input quark flavour reset to u" << std::endl;
59 id1Now = 2;
60 particleSPtr->channel(0).product(0, id1Now);
61 }
62 if (std::abs(id2Now) < 11 || std::abs(id2Now) > 16) {
63 std::cout << "ERROR in ResonanceLQ::init:unallowed input lepton flavour reset to e-" << std::endl;
64 id2Now = 11;
65 particleSPtr->channel(0).product(1, id2Now);
66 }
67
68
69 bool changed = particleSPtr->hasChanged();
70
71 int chargeLQ = particleDataPtr->chargeType(id1Now)
72 + particleDataPtr->chargeType(id2Now);
73
74 particleSPtr->setChargeType(chargeLQ);
75
76 std::string nameLQ = "LQ_" + particleDataPtr->name(id1Now) + ","
77 + particleDataPtr->name(id2Now);
78
79 particleSPtr->setNames(nameLQ, nameLQ + "bar");
80 if (!changed) particleSPtr->setHasChanged(false);
81
82 return;
83 }
Pythia8::ParticleDataEntry * ParticleDataEntryPtr