ATLAS Offline Software
Loading...
Searching...
No Matches
Public Member Functions | Private Attributes | List of all members
Pythia8::PowhegBB4L Class Reference
Inheritance diagram for Pythia8::PowhegBB4L:
Collaboration diagram for Pythia8::PowhegBB4L:

Public Member Functions

 PowhegBB4L ()
 ~PowhegBB4L ()
double findresscale (const int iRes, const Event &event)
bool match_decay (int iparticle, const Event &e, const std::vector< int > &ids, std::vector< int > &positions, std::vector< Vec4 > &momenta, bool exitOnExtraLegs=true)
double qSplittingScale (Vec4 pt, Vec4 p1, Vec4 p2)
double gSplittingScale (Vec4 pt, Vec4 p1, Vec4 p2)
double getdechardness (int topcharge, const Event &e)
bool canVetoPartonLevel ()
bool doVetoPartonLevel (const Event &e)
bool canSetResonanceScale ()
double scaleResonance (int iRes, const Event &e)
bool canVetoFSREmission ()
bool doVetoFSREmission (int, const Event &e, int, bool inResonance)

Private Attributes

double m_topresscale
double m_atopresscale
Pythia8_UserHooks::UserSetting< int > m_onlyDistance1
Pythia8_UserHooks::UserSetting< int > m_useScaleResonanceInstead

Detailed Description

Definition at line 22 of file PowhegBB4L.cxx.

Constructor & Destructor Documentation

◆ PowhegBB4L()

Pythia8::PowhegBB4L::PowhegBB4L ( )
inline

Definition at line 28 of file PowhegBB4L.cxx.

28 : m_topresscale(-1.), m_atopresscale(-1.),
29 m_onlyDistance1("Powheg:bb4l:onlyDistance1", 0),
30 m_useScaleResonanceInstead("Powheg:bb4l:useScaleResonanceInstead", 0){}
Pythia8::PowhegBB4L::m_onlyDistance1
Pythia8_UserHooks::UserSetting< int > m_onlyDistance1
Definition PowhegBB4L.cxx:314
Pythia8::PowhegBB4L::m_useScaleResonanceInstead
Pythia8_UserHooks::UserSetting< int > m_useScaleResonanceInstead
Definition PowhegBB4L.cxx:315

◆ ~PowhegBB4L()

Pythia8::PowhegBB4L::~PowhegBB4L ( )
inline

Definition at line 31 of file PowhegBB4L.cxx.

31{}

Member Function Documentation

◆ canSetResonanceScale()

bool Pythia8::PowhegBB4L::canSetResonanceScale ( )
inline

Definition at line 215 of file PowhegBB4L.cxx.

215{ return true; }

◆ canVetoFSREmission()

bool Pythia8::PowhegBB4L::canVetoFSREmission ( )
inline

Definition at line 233 of file PowhegBB4L.cxx.

233 {
234 if (m_useScaleResonanceInstead(settingsPtr))
235 return false;
236 else
237 return true;
238 }

◆ canVetoPartonLevel()

bool Pythia8::PowhegBB4L::canVetoPartonLevel ( )
inline

Definition at line 196 of file PowhegBB4L.cxx.

196{ return true; }

◆ doVetoFSREmission()

bool Pythia8::PowhegBB4L::doVetoFSREmission ( int ,
const Event & e,
int ,
bool inResonance )
inline

Definition at line 242 of file PowhegBB4L.cxx.

242 {
243
244 if (inResonance) {
245
246 // get the participants of the splitting: the radiator and the emitted
247 int iEmt = e.size() - 2;
248 int iRadAft = e.size() - 3;
249 int iRadBef = e[iEmt].mother1();
250
251 // find the top resonance the radiator originates from
252 int iTop = e[iRadBef].mother1();
253 int distance = 1;
254 while (std::abs(e[iTop].id()) != 6 && iTop > 0) {
255 iTop = e[iTop].mother1();
256 distance ++;
257 }
258 if (iTop == 0) {
259#if PYTHIA_VERSION_INTEGER >= 8310
260 loggerPtr->ERROR_MSG("Warning in PowhegHooksBB4L::doVetoFSREmission: emission in resonance not from top quark, not vetoing");
261#else
262 infoPtr->errorMsg("Warning in PowhegHooksBB4L::doVetoFSREmission: emission in resonance not from top quark, not vetoing");
263#endif
264 return false;
265 }
266 int iTopCharge = (e[iTop].id()>0)?1:-1;
267
268
269 // calculate the scale of the emission
270 Vec4 pr(e[iRadAft].p()), pe(e[iEmt].p()), pt(e[iTop].p());
271 double scale;
272 // gluon splitting into two partons
273 if (e[iRadBef].id() == 21)
274 scale = gSplittingScale(pt, pr, pe);
275 // quark emitting a gluon
276 else if (std::abs(e[iRadBef].id()) <= 5)
277 scale = qSplittingScale(pt, pr, pe);
278 // other stuff (which we should not veto)
279 else {
280 scale = 0;
281 }
282
283 if (iTopCharge > 0) {
284 if (m_onlyDistance1(settingsPtr))
285 return (distance == 1) && scale > m_topresscale;
286 else
287 return scale > m_topresscale;
288 }
289 else {
290 if (m_onlyDistance1(settingsPtr))
291 return (distance == 1) && scale > m_atopresscale;
292 else
293 return scale > m_atopresscale;
294 }
295 } else {
296 return false;
297 }
298
299 return false;
300 }
Pythia8::PowhegBB4L::gSplittingScale
double gSplittingScale(Vec4 pt, Vec4 p1, Vec4 p2)
Definition PowhegBB4L.cxx:119
Pythia8::PowhegBB4L::qSplittingScale
double qSplittingScale(Vec4 pt, Vec4 p1, Vec4 p2)
Definition PowhegBB4L.cxx:113
Amg::distance
float distance(const Amg::Vector3D &p1, const Amg::Vector3D &p2)
calculates the distance between two point in 3D space
Definition GeoPrimitivesHelpers.h:54
test_pyathena.pt
pt
Definition test_pyathena.py:11
yodamerge_tmp.scale
scale
Definition yodamerge_tmp.py:138

◆ doVetoPartonLevel()

bool Pythia8::PowhegBB4L::doVetoPartonLevel ( const Event & e)
inline

Definition at line 197 of file PowhegBB4L.cxx.

197 {
198 double topdechardness = getdechardness(1, e), atopdechardness = getdechardness(-1, e);
199 if ((topdechardness > m_topresscale) or (atopdechardness > m_atopresscale)) {
200
201#if PYTHIA_VERSION_INTEGER >= 8310
202 loggerPtr->ERROR_MSG("Warning in PowhegHooksBB4L::doVetoPartonLevel: passed doVetoFSREmission veto, but wouldn't have passed veto based on the full event listing");
203#else
204 infoPtr->errorMsg("Warning in PowhegHooksBB4L::doVetoPartonLevel: passed doVetoFSREmission veto, but wouldn't have passed veto based on the full event listing");
205#endif
206 }
207// cout << " veto scales: " << fixed << setprecision(17) << setw(30) << topresscale << setw(30) << atopresscale << endl;
208 m_topresscale = -1;
209 m_atopresscale = -1;
210
211 return false;
212 }
Pythia8::PowhegBB4L::getdechardness
double getdechardness(int topcharge, const Event &e)
Definition PowhegBB4L.cxx:126

◆ findresscale()

double Pythia8::PowhegBB4L::findresscale ( const int iRes,
const Event & event )
inline

Definition at line 38 of file PowhegBB4L.cxx.

38 {
39 double scale = 0.;
40
41 int nDau = event[iRes].daughterList().size();
42
43 if (nDau == 0) {
44 // No resonance found, set scale to high value
45 // Pythia will shower any MC generated resonance unrestricted
46 scale = 1e30;
47 }
48 else if (nDau < 3) {
49 // No radiating resonance found
50 scale = 0.8;
51 }
52 else if (std::abs(event[iRes].id()) == 6) {
53 // Find top daughters
54 int idw = -1, idb = -1, idg = -1;
55
56 for (int i = 0; i < nDau; i++) {
57 int iDau = event[iRes].daughterList()[i];
58 if (std::abs(event[iDau].id()) == 24) idw = iDau;
59 if (std::abs(event[iDau].id()) == 5) idb = iDau;
60 if (std::abs(event[iDau].id()) == 21) idg = iDau;
61 }
62
63 // Get daughter 4-vectors in resonance frame
64 Vec4 pw(event[idw].p());
65 pw.bstback(event[iRes].p());
66
67 Vec4 pb(event[idb].p());
68 pb.bstback(event[iRes].p());
69
70 Vec4 pg(event[idg].p());
71 pg.bstback(event[iRes].p());
72
73 // Calculate scale
74 scale = std::sqrt(2*pg*pb*pg.e()/pb.e());
75 }
76 else {
77 scale = 1e30;
78 }
79
80 return scale;
81 }
ParticleGun_EoverP_Config.pg
pg
Definition ParticleGun_EoverP_Config.py:61
keylayer_zslicemap.pb
dict pb
Definition keylayer_zslicemap.py:188

◆ getdechardness()

double Pythia8::PowhegBB4L::getdechardness ( int topcharge,
const Event & e )
inline

Definition at line 126 of file PowhegBB4L.cxx.

126 {
127 int tid = 6*topcharge, wid = 24*topcharge, bid = 5*topcharge, gid = 21, wildcard = 0;
128 // find last top in the record
129 int i_top = -1;
130 Vec4 p_top, p_b, p_g, p_g1, p_g2;
131 for (int i = 0; i < e.size(); i++)
132 if (e[i].id() == tid) {
133 i_top = i;
134 p_top = e[i].p();
135 }
136 if (i_top == -1) return -1.0;
137
138 // summary of cases
139 // 1.) t > W b
140 // a.) b > 3 ... error
141 // b.) b > b g ... h = sqrt(2*p_g*p_b*p_g.e()/p_b.e())
142 // c.) b > other ... h = -1
143 // return h
144 // 2.) t > W b g
145 // a.) b > 3 ... error
146 // b.) b > b g ... h1 = sqrt(2*p_g*p_b*p_g.e()/p_b.e())
147 // c.) b > other ... h1 = -1
148 // i.) g > 3 ... error
149 // ii.) g > 2 ... h2 = sqrt(2*p_g1*p_g2*p_g1.e()*p_g2.e()/(pow(p_g1.e(),2)+pow(p_g2.e(),2))) );
150 // iii.) g > other ... h2 = -1
151 // return max(h1,h2)
152 // 3.) else ... error
153
154 vector<Vec4> momenta;
155 vector<int> positions;
156
157 // 1.) t > b W
158 if ( match_decay(i_top, e, std::vector<int> {wid, bid}, positions, momenta, false) ) {
159 double h;
160 int i_b = positions[1];
161 // a.+b.) b > 3 or b > b g
162 if ( match_decay(i_b, e, std::vector<int> {bid, gid}, positions, momenta) )
163 h = qSplittingScale(e[i_top].p(), momenta[0], momenta[1]);
164 // c.) b > other
165 else
166 h = -1;
167 return h;
168 }
169 // 2.) t > b W g
170 else if ( match_decay(i_top, e, std::vector<int> {wid, bid, gid}, positions, momenta, false) ) {
171 double h1, h2;
172 int i_b = positions[1], i_g = positions[2];
173 // a.+b.) b > 3 or b > b g
174 if ( match_decay(i_b, e, std::vector<int> {bid, gid}, positions, momenta) )
175 h1 = qSplittingScale(e[i_top].p(), momenta[0], momenta[1]);
176 // c.) b > other
177 else
178 h1 = -1;
179 // i.+ii.) g > 3 or g > 2
180 if ( match_decay(i_g, e, std::vector<int> {wildcard, wildcard}, positions, momenta) )
181 h2 = gSplittingScale(e[i_top].p(), momenta[0], momenta[1]);
182 // c.) b > other
183 else
184 h2 = -1;
185 return max(h1, h2);
186 }
187 // 3.) else
188 else {
189 exit(-1);
190 }
191 }
gid
Acts::GeometryIdentifier gid
Definition ActsWriteTrackingGeometryTransforms.cxx:30
max
#define max(a, b)
Definition cfImp.cxx:41
Pythia8::PowhegBB4L::match_decay
bool match_decay(int iparticle, const Event &e, const std::vector< int > &ids, std::vector< int > &positions, std::vector< Vec4 > &momenta, bool exitOnExtraLegs=true)
Definition PowhegBB4L.cxx:89
read_hist_ntuple.h1
h1
Definition read_hist_ntuple.py:21

◆ gSplittingScale()

double Pythia8::PowhegBB4L::gSplittingScale ( Vec4 pt,
Vec4 p1,
Vec4 p2 )
inline

Definition at line 119 of file PowhegBB4L.cxx.

119 {
120 p1.bstback(pt);
121 p2.bstback(pt);
122 return std::sqrt( 2*p1*p2*p1.e()*p2.e()/(std::pow(p1.e(),2)+std::pow(p2.e(),2)) );
123 }
TRTCalib_cfilter.p2
p2
Definition TRTCalib_cfilter.py:131
TRTCalib_cfilter.p1
p1
Definition TRTCalib_cfilter.py:130

◆ match_decay()

bool Pythia8::PowhegBB4L::match_decay ( int iparticle,
const Event & e,
const std::vector< int > & ids,
std::vector< int > & positions,
std::vector< Vec4 > & momenta,
bool exitOnExtraLegs = true )
inline

Definition at line 89 of file PowhegBB4L.cxx.

89 {
90 // compare sizes
91 if (e[iparticle].daughterList().size() != ids.size()) {
92 if (exitOnExtraLegs && e[iparticle].daughterList().size() > ids.size()) exit(-1);
93 return false;
94 }
95 // compare content
96 for (size_t i = 0; i < e[iparticle].daughterList().size(); i++) {
97 int di = e[iparticle].daughterList()[i];
98 if (ids[i] != 0 && e[di].id() != ids[i])
99 return false;
100 }
101 // reset the positions and momenta vectors (because they may be reused)
102 positions.clear();
103 momenta.clear();
104 // construct the array of momenta
105 for (size_t i = 0; i < e[iparticle].daughterList().size(); i++) {
106 int di = e[iparticle].daughterList()[i];
107 positions.push_back(di);
108 momenta.push_back(e[di].p());
109 }
110 return true;
111 }

◆ qSplittingScale()

double Pythia8::PowhegBB4L::qSplittingScale ( Vec4 pt,
Vec4 p1,
Vec4 p2 )
inline

Definition at line 113 of file PowhegBB4L.cxx.

113 {
114 p1.bstback(pt);
115 p2.bstback(pt);
116 return std::sqrt( 2*p1*p2*p2.e()/p1.e() );
117 }

◆ scaleResonance()

double Pythia8::PowhegBB4L::scaleResonance ( int iRes,
const Event & e )
inline

Definition at line 216 of file PowhegBB4L.cxx.

216 {
217 double scale = 0.;
218
219 if (e[iRes].id() == 6){
220 scale = m_topresscale = findresscale(iRes, e);
221 }else if (e[iRes].id() == -6){
222 scale = m_atopresscale = findresscale(iRes, e);
223 }
224
225 if (m_useScaleResonanceInstead(settingsPtr)) return scale;
226
227 return 1e30;
228 }
Pythia8::PowhegBB4L::findresscale
double findresscale(const int iRes, const Event &event)
Definition PowhegBB4L.cxx:38

Member Data Documentation

◆ m_atopresscale

double Pythia8::PowhegBB4L::m_atopresscale
private

Definition at line 313 of file PowhegBB4L.cxx.

◆ m_onlyDistance1

Pythia8_UserHooks::UserSetting<int> Pythia8::PowhegBB4L::m_onlyDistance1
private

Definition at line 314 of file PowhegBB4L.cxx.

◆ m_topresscale

double Pythia8::PowhegBB4L::m_topresscale
private

Definition at line 312 of file PowhegBB4L.cxx.

◆ m_useScaleResonanceInstead

Pythia8_UserHooks::UserSetting<int> Pythia8::PowhegBB4L::m_useScaleResonanceInstead
private

Definition at line 315 of file PowhegBB4L.cxx.

The documentation for this class was generated from the following file:
Generated on for ATLAS Offline Software by doxygen 1.14.0