Rivet analyses referenceLHCB_2014_I1281685Charged particle multiplicities and densities in $pp$ collisions at $\sqrt{s} = 7$ TeVExperiment: LHCB (LHC) Inspire ID: 1281685 Status: VALIDATED Authors:
Beam energies: (3500.0, 3500.0) GeV Run details:
Measurement of prompt charged particle production in proton-proton collisions at a centre-of-mass energy of 7 TeV at the LHC. The measurement is performed in the kinematic range $p > 2$ GeV/$c$, $p_\perp > 0.2$ GeV/$c$ and $2.0 < \eta < 4.8$. At least one prompt charged particle per event is required. A prompt particle is defined as a particle, that either originates directly from the primary vertex or from a decay chain in which the sum of mean lifetimes does not exceed 10 ps. Charged particle multiplicities are measured in the full kinematic range as well as in regions of pseudorapidity and transverse momentum. In addition, also mean particle densities are presented as functions of pseudorapidity and transverse momentum. Source code: LHCB_2014_I1281685.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/ChargedFinalState.hh"
4
5namespace Rivet {
6
7
8 /// Charged particle multiplicities and densities in $pp$ collisions at $\sqrt{s} = 7$ TeV
9 class LHCB_2014_I1281685 : public Analysis {
10 public:
11
12 /// @name Constructors etc.
13 /// @{
14
15 /// Constructor
16 LHCB_2014_I1281685()
17 : Analysis("LHCB_2014_I1281685"),
18 _p_min(2.0),
19 _pt_min(0.2),
20 _eta_min(2.0),
21 _eta_max(4.8),
22 _maxlft(1.0e-11)
23 { }
24
25 /// @}
26
27
28 /// @name Analysis methods
29 /// @{
30
31 /// Book histograms and initialise projections before the run
32 void init() {
33 fillMap(_partLftMap);
34
35 // Projections
36 declare(ChargedFinalState(Cuts::etaIn(_eta_min, _eta_max) && Cuts::pT >= _pt_min*GeV),
37 "CFS");
38
39 // Book histograms
40 book(_h_mult_total ,"d03-x01-y01", 50, 0.5, 50.5);
41
42 book(_h_mult_eta[0] ,"d04-x01-y01", 21, -0.5, 20.5); //eta=[2.0,2.5]
43 book(_h_mult_eta[1] ,"d04-x01-y02", 21, -0.5, 20.5); //eta=[2.5,3.0]
44 book(_h_mult_eta[2] ,"d04-x01-y03", 21, -0.5, 20.5); //eta=[3.0,3.5]
45 book(_h_mult_eta[3] ,"d04-x01-y04", 21, -0.5, 20.5); //eta=[3.5,4.0]
46 book(_h_mult_eta[4] ,"d04-x01-y05", 21, -0.5, 20.5); //eta=[4.0,4.5]
47
48 book(_h_mult_pt[0] ,"d05-x01-y01", 21, -0.5, 20.5); //pT=[0.2,0.3]GeV
49 book(_h_mult_pt[1] ,"d05-x01-y02", 21, -0.5, 20.5); //pT=[0.3,0.4]GeV
50 book(_h_mult_pt[2] ,"d05-x01-y03", 21, -0.5, 20.5); //pT=[0.4,0.6]GeV
51 book(_h_mult_pt[3] ,"d05-x01-y04", 21, -0.5, 20.5); //pT=[0.6,1.0]GeV
52 book(_h_mult_pt[4] ,"d05-x01-y05", 21, -0.5, 20.5); //pT=[1.0,2.0]GeV
53
54 book(_h_dndeta ,"d01-x01-y01", 14, 2.0, 4.8); //eta=[2,4.8]
55 book(_h_dndpt ,"d02-x01-y01", 18, 0.2, 2.0); //pT =[0,2]GeV
56
57 // Counters
58 book(_sumW, "TMP/sumW");
59 }
60
61
62 /// Perform the per-event analysis
63 void analyze(const Event& event) {
64
65 // Variable to store multiplicities per event
66 int LHCbcountAll = 0; //count particles fulfiling all requirements
67 int LHCbcountEta[8] = {0,0,0,0,0,0,0,0}; //count per eta-bin
68 int LHCbcountPt[7] = {0,0,0,0,0,0,0}; //count per pT-bin
69 vector<double> val_dNdEta;
70 vector<double> val_dNdPt;
71 val_dNdEta.clear();
72 val_dNdPt.clear();
73
74 const ChargedFinalState& cfs = apply<ChargedFinalState>(event, "CFS");
75 for (const Particle& p : cfs.particles()) {
76 int id = p.pid();
77 // continue if particle is not a pion, kaon, proton, muon or electron
78 if ( !( (abs(id) == 211) || (abs(id) == 321) || (abs(id) == 2212) || (abs(id) == 13) || (abs(id) == 11)) ) {
79 continue;
80 }
81
82 const FourMomentum& qmom = p.momentum();
83 const double eta = p.momentum().eta();
84 const double pT = p.momentum().pT();
85 //minimum momentum
86 if (qmom.p3().mod() < _p_min) continue;
87 //minimum tr. momentum
88 if (pT < _pt_min) continue;
89 //eta range
90 if ((eta < _eta_min) || (eta > _eta_max)) continue;
91
92 /* Select only prompt particles via lifetime */
93 //Sum of all mother lifetimes (PDG lifetime) < 10ps
94 double ancestors_sumlft = getAncestorSumLifetime(p);
95 if( (ancestors_sumlft > _maxlft) || (ancestors_sumlft < 0) ) continue;
96
97 //after all cuts;
98 LHCbcountAll++; //count particles in whole kin. range
99
100 //in eta bins
101 if( eta >2.0 && eta <= 2.5) LHCbcountEta[0]++;
102 if( eta >2.5 && eta <= 3.0) LHCbcountEta[1]++;
103 if( eta >3.0 && eta <= 3.5) LHCbcountEta[2]++;
104 if( eta >3.5 && eta <= 4.0) LHCbcountEta[3]++;
105 if( eta >4.0 && eta <= 4.5) LHCbcountEta[4]++;
106 if( eta >2.0 && eta <= 4.8) LHCbcountEta[5]++; //cross-check
107 //in pT bins
108 if( pT > 0.2 && pT <= 0.3) LHCbcountPt[0]++;
109 if( pT > 0.3 && pT <= 0.4) LHCbcountPt[1]++;
110 if( pT > 0.4 && pT <= 0.6) LHCbcountPt[2]++;
111 if( pT > 0.6 && pT <= 1.0) LHCbcountPt[3]++;
112 if( pT > 1.0 && pT <= 2.0) LHCbcountPt[4]++;
113 if( pT > 0.2) LHCbcountPt[5]++; //cross-check
114
115 //particle densities -> need proper normalization (finalize)
116 val_dNdPt.push_back( pT );
117 val_dNdEta.push_back( eta );
118 }//end for
119
120
121 // Fill histograms only, if at least 1 particle pre event was within the
122 // kinematic range of the analysis!
123 if (LHCbcountAll) {
124 _sumW->fill();
125
126 _h_mult_total->fill(LHCbcountAll);
127
128 _h_mult_eta[0]->fill(LHCbcountEta[0]);
129 _h_mult_eta[1]->fill(LHCbcountEta[1]);
130 _h_mult_eta[2]->fill(LHCbcountEta[2]);
131 _h_mult_eta[3]->fill(LHCbcountEta[3]);
132 _h_mult_eta[4]->fill(LHCbcountEta[4]);
133
134 _h_mult_pt[0]->fill(LHCbcountPt[0]);
135 _h_mult_pt[1]->fill(LHCbcountPt[1]);
136 _h_mult_pt[2]->fill(LHCbcountPt[2]);
137 _h_mult_pt[3]->fill(LHCbcountPt[3]);
138 _h_mult_pt[4]->fill(LHCbcountPt[4]);
139
140 for (size_t part = 0; part < val_dNdEta.size(); part++)
141 _h_dndeta->fill(val_dNdEta[part]);
142
143 for (size_t part = 0; part < val_dNdPt.size(); part++)
144 _h_dndpt->fill(val_dNdPt[part]);
145
146 }
147 }
148
149
150 /// Normalise histograms etc., after the run
151 void finalize() {
152 const double scalefactor = 1.0/_sumW->val(); // normalize multiplicity histograms by nEvents
153 const double scale1k = 1000.; // to match '10^3' scale in reference histograms
154
155 scale( _h_dndeta, scalefactor );
156 scale( _h_dndpt, scalefactor*0.1 ); //additional factor 0.1 for [0.1 GeV/c]
157 scale( _h_mult_total, scalefactor*scale1k);
158
159 _h_mult_eta[0]->scaleW( scalefactor*scale1k );
160 _h_mult_eta[1]->scaleW( scalefactor*scale1k );
161 _h_mult_eta[2]->scaleW( scalefactor*scale1k );
162 _h_mult_eta[3]->scaleW( scalefactor*scale1k );
163 _h_mult_eta[4]->scaleW( scalefactor*scale1k );
164
165 _h_mult_pt[0]->scaleW( scalefactor*scale1k );
166 _h_mult_pt[1]->scaleW( scalefactor*scale1k );
167 _h_mult_pt[2]->scaleW( scalefactor*scale1k );
168 _h_mult_pt[3]->scaleW( scalefactor*scale1k );
169 _h_mult_pt[4]->scaleW( scalefactor*scale1k );
170 }
171
172 /// @}
173
174
175 private:
176
177
178 // Get mean PDG lifetime for particle with PID
179 double getLifetime(int pid) {
180 double lft = 0.;
181 map<int, double>::iterator pPartLft = _partLftMap.find(pid);
182 if (pPartLft != _partLftMap.end()) {
183 lft = (*pPartLft).second;
184 } else {
185 // allow identifying missing life times only in debug mode
186 MSG_DEBUG("Could not determine lifetime for particle with PID " << pid << "... Assume non-prompt particle");
187 lft = -1;
188 }
189 return lft;
190 }
191
192
193 // Get sum of all ancestor particles
194 double getAncestorSumLifetime(const Particle& p) {
195 double lftSum = 0.;
196 double plft = 0.;
197 ConstGenParticlePtr part = p.genParticle();
198 if ( 0 == part ) return -1;
199 ConstGenVertexPtr ivtx = part->production_vertex();
200 while(ivtx) {
201 vector<ConstGenParticlePtr> part_in = HepMCUtils::particles(ivtx, Relatives::PARENTS);
202 if (part_in.size() < 1) { lftSum = -1.; break; };
203 part = part_in.at(0);
204 if ( !(part) ) { lftSum = -1.; break; };
205 ivtx = part->production_vertex();
206 if ( (part->pdg_id() == 2212) || !(ivtx) ) break; // reached beam
207 plft = getLifetime(part->pdg_id());
208 if (plft < 0.) { lftSum = -1.; break; };
209 lftSum += plft;
210 }
211 return (lftSum);
212 }
213
214
215 /// Hard-coded map linking PDG ID with PDG lifetime[s] (converted from ParticleTable.txt)
216 bool fillMap(map<int, double>& m) {
217 // PDGID = LIFETIME
218 m[22] = 1.000000e+016;
219 m[-11] = 1.000000e+016;
220 m[11] = 1.000000e+016;
221 m[12] = 1.000000e+016;
222 m[-13] = 2.197036e-006;
223 m[13] = 2.197036e-006;
224 m[111] = 8.438618e-017;
225 m[211] = 2.603276e-008;
226 m[-211] = 2.603276e-008;
227 m[130] = 5.174624e-008;
228 m[321] = 1.238405e-008;
229 m[-321] = 1.238405e-008;
230 m[2112] = 885.646128;
231 m[2212] = 1.000000e+016;
232 m[-2212] = 1.000000e+016;
233 m[310] = 8.934603e-011;
234 m[221] = 5.578070e-019;
235 m[3122] = 2.631796e-010;
236 m[3222] = 8.018178e-011;
237 m[3212] = 7.395643e-020;
238 m[3112] = 1.479129e-010;
239 m[3322] = 2.899613e-010;
240 m[3312] = 1.637344e-010;
241 m[3334] = 8.207135e-011;
242 m[-2112] = 885.646128;
243 m[-3122] = 2.631796e-010;
244 m[-3222] = 8.018178e-011;
245 m[-3212] = 7.395643e-020;
246 m[-3112] = 1.479129e-010;
247 m[-3322] = 2.899613e-010;
248 m[-3312] = 1.637344e-010;
249 m[-3334] = 8.207135e-011;
250 m[113] = 4.411610e-024;
251 m[213] = 4.411610e-024;
252 m[-213] = 4.411610e-024;
253 m[223] = 7.798723e-023;
254 m[333] = 1.545099e-022;
255 m[323] = 1.295693e-023;
256 m[-323] = 1.295693e-023;
257 m[313] = 1.298249e-023;
258 m[-313] = 1.298249e-023;
259 m[20213] = 1.500000e-024;
260 m[-20213] = 1.500000e-024;
261 m[450000000] = 1.000000e+015;
262 m[460000000] = 1.000000e+015;
263 m[470000000] = 1.000000e+015;
264 m[480000000] = 1.000000e+015;
265 m[490000000] = 1.000000e+015;
266 m[20022] = 1.000000e+016;
267 m[-15] = 2.906014e-013;
268 m[15] = 2.906014e-013;
269 m[24] = 3.104775e-025;
270 m[-24] = 3.104775e-025;
271 m[23] = 2.637914e-025;
272 m[411] = 1.051457e-012;
273 m[-411] = 1.051457e-012;
274 m[421] = 4.116399e-013;
275 m[-421] = 4.116399e-013;
276 m[431] = 4.904711e-013;
277 m[-431] = 4.904711e-013;
278 m[4122] = 1.994582e-013;
279 m[-4122] = 1.994582e-013;
280 m[443] = 7.565657e-021;
281 m[413] = 6.856377e-021;
282 m[-413] = 6.856377e-021;
283 m[423] = 1.000003e-019;
284 m[-423] = 1.000003e-019;
285 m[433] = 1.000003e-019;
286 m[-433] = 1.000003e-019;
287 m[521] = 1.671000e-012;
288 m[-521] = 1.671000e-012;
289 m[511] = 1.536000e-012;
290 m[-511] = 1.536000e-012;
291 m[531] = 1.461000e-012;
292 m[-531] = 1.461000e-012;
293 m[541] = 4.600000e-013;
294 m[-541] = 4.600000e-013;
295 m[5122] = 1.229000e-012;
296 m[-5122] = 1.229000e-012;
297 m[4112] = 4.388081e-022;
298 m[-4112] = 4.388081e-022;
299 m[4212] = 3.999999e-022;
300 m[-4212] = 3.999999e-022;
301 m[4222] = 3.291060e-022;
302 m[-4222] = 3.291060e-022;
303 m[25] = 9.400000e-026;
304 m[35] = 9.400000e-026;
305 m[36] = 9.400000e-026;
306 m[37] = 9.400000e-026;
307 m[-37] = 9.400000e-026;
308 m[4312] = 9.800002e-014;
309 m[-4312] = 9.800002e-014;
310 m[4322] = 3.500001e-013;
311 m[-4322] = 3.500001e-013;
312 m[4332] = 6.453061e-014;
313 m[-4332] = 6.453061e-014;
314 m[4132] = 9.824063e-014;
315 m[-4132] = 9.824063e-014;
316 m[4232] = 4.417532e-013;
317 m[-4232] = 4.417532e-013;
318 m[5222] = 1.000000e-019;
319 m[-5222] = 1.000000e-019;
320 m[5212] = 1.000000e-019;
321 m[-5212] = 1.000000e-019;
322 m[5112] = 1.000000e-019;
323 m[-5112] = 1.000000e-019;
324 m[5312] = 1.000000e-019;
325 m[-5312] = 1.000000e-019;
326 m[5322] = 1.000000e-019;
327 m[-5322] = 1.000000e-019;
328 m[5332] = 1.550000e-012;
329 m[-5332] = 1.550000e-012;
330 m[5132] = 1.390000e-012;
331 m[-5132] = 1.390000e-012;
332 m[5232] = 1.390000e-012;
333 m[-5232] = 1.390000e-012;
334 m[100443] = 2.194041e-021;
335 m[331] = 3.258476e-021;
336 m[441] = 4.113826e-023;
337 m[10441] = 4.063038e-023;
338 m[20443] = 7.154480e-022;
339 m[445] = 3.164482e-022;
340 m[9000111] = 1.149997e-023;
341 m[9000211] = 1.149997e-023;
342 m[-9000211] = 1.149997e-023;
343 m[20113] = 1.500000e-024;
344 m[115] = 6.151516e-024;
345 m[215] = 6.151516e-024;
346 m[-215] = 6.151516e-024;
347 m[10323] = 7.313469e-024;
348 m[-10323] = 7.313469e-024;
349 m[10313] = 7.313469e-024;
350 m[-10313] = 7.313469e-024;
351 m[20323] = 3.782829e-024;
352 m[-20323] = 3.782829e-024;
353 m[20313] = 3.782829e-024;
354 m[-20313] = 3.782829e-024;
355 m[10321] = 2.238817e-024;
356 m[-10321] = 2.238817e-024;
357 m[10311] = 2.238817e-024;
358 m[-10311] = 2.238817e-024;
359 m[325] = 6.682357e-024;
360 m[-325] = 6.682357e-024;
361 m[315] = 6.038644e-024;
362 m[-315] = 6.038644e-024;
363 m[10411] = 4.380000e-024;
364 m[20413] = 2.630000e-024;
365 m[10413] = 3.290000e-023;
366 m[-415] = 2.632849e-023;
367 m[-10411] = 4.380000e-024;
368 m[-20413] = 2.630000e-024;
369 m[-10413] = 3.290000e-023;
370 m[415] = 2.632849e-023;
371 m[10421] = 4.380000e-024;
372 m[20423] = 2.630000e-024;
373 m[10423] = 3.482604e-023;
374 m[-425] = 2.861792e-023;
375 m[-10421] = 4.380000e-024;
376 m[-20423] = 2.630000e-024;
377 m[-10423] = 3.482604e-023;
378 m[425] = 2.861792e-023;
379 m[10431] = 6.582100e-022;
380 m[20433] = 6.582100e-022;
381 m[10433] = 6.582100e-022;
382 m[435] = 4.388100e-023;
383 m[-10431] = 6.582100e-022;
384 m[-20433] = 6.582100e-022;
385 m[-10433] = 6.582100e-022;
386 m[-435] = 4.388100e-023;
387 m[2224] = 5.485102e-024;
388 m[2214] = 5.485102e-024;
389 m[2114] = 5.485102e-024;
390 m[1114] = 5.485102e-024;
391 m[-2224] = 5.485102e-024;
392 m[-2214] = 5.485102e-024;
393 m[-2114] = 5.485102e-024;
394 m[-1114] = 5.485102e-024;
395 m[-523] = 1.000019e-019;
396 m[523] = 1.000019e-019;
397 m[513] = 1.000019e-019;
398 m[-513] = 1.000019e-019;
399 m[533] = 1.000000e-019;
400 m[-533] = 1.000000e-019;
401 m[10521] = 4.390000e-024;
402 m[20523] = 2.630000e-024;
403 m[10523] = 1.650000e-023;
404 m[525] = 1.310000e-023;
405 m[-10521] = 4.390000e-024;
406 m[-20523] = 2.630000e-024;
407 m[-10523] = 1.650000e-023;
408 m[-525] = 1.310000e-023;
409 m[10511] = 4.390000e-024;
410 m[20513] = 2.630000e-024;
411 m[10513] = 1.650000e-023;
412 m[515] = 1.310000e-023;
413 m[-10511] = 4.390000e-024;
414 m[-20513] = 2.630000e-024;
415 m[-10513] = 1.650000e-023;
416 m[-515] = 1.310000e-023;
417 m[10531] = 4.390000e-024;
418 m[20533] = 2.630000e-024;
419 m[10533] = 1.650000e-023;
420 m[535] = 1.310000e-023;
421 m[-10531] = 4.390000e-024;
422 m[-20533] = 2.630000e-024;
423 m[-10533] = 1.650000e-023;
424 m[-535] = 1.310000e-023;
425 m[14] = 1.000000e+016;
426 m[-14] = 1.000000e+016;
427 m[-12] = 1.000000e+016;
428 m[1] = 0.000000e+000;
429 m[-1] = 0.000000e+000;
430 m[2] = 0.000000e+000;
431 m[-2] = 0.000000e+000;
432 m[3] = 0.000000e+000;
433 m[-3] = 0.000000e+000;
434 m[4] = 0.000000e+000;
435 m[-4] = 0.000000e+000;
436 m[5] = 0.000000e+000;
437 m[-5] = 0.000000e+000;
438 m[6] = 4.707703e-025;
439 m[-6] = 4.707703e-025;
440 m[7] = 0.000000e+000;
441 m[-7] = 0.000000e+000;
442 m[8] = 0.000000e+000;
443 m[-8] = 0.000000e+000;
444 m[16] = 1.000000e+016;
445 m[-16] = 1.000000e+016;
446 m[17] = 0.000000e+000;
447 m[-17] = 0.000000e+000;
448 m[18] = 0.000000e+000;
449 m[-18] = 0.000000e+000;
450 m[21] = 0.000000e+000;
451 m[32] = 0.000000e+000;
452 m[33] = 0.000000e+000;
453 m[34] = 0.000000e+000;
454 m[-34] = 0.000000e+000;
455 m[39] = 0.000000e+000;
456 m[41] = 0.000000e+000;
457 m[-41] = 0.000000e+000;
458 m[42] = 0.000000e+000;
459 m[-42] = 0.000000e+000;
460 m[43] = 0.000000e+000;
461 m[44] = 0.000000e+000;
462 m[-44] = 0.000000e+000;
463 m[81] = 0.000000e+000;
464 m[82] = 0.000000e+000;
465 m[-82] = 0.000000e+000;
466 m[83] = 0.000000e+000;
467 m[84] = 3.335641e-013;
468 m[-84] = 3.335641e-013;
469 m[85] = 1.290893e-012;
470 m[-85] = 1.290893e-012;
471 m[86] = 0.000000e+000;
472 m[-86] = 0.000000e+000;
473 m[87] = 0.000000e+000;
474 m[-87] = 0.000000e+000;
475 m[88] = 0.000000e+000;
476 m[90] = 0.000000e+000;
477 m[91] = 0.000000e+000;
478 m[92] = 0.000000e+000;
479 m[93] = 0.000000e+000;
480 m[94] = 0.000000e+000;
481 m[95] = 0.000000e+000;
482 m[96] = 0.000000e+000;
483 m[97] = 0.000000e+000;
484 m[98] = 0.000000e+000;
485 m[99] = 0.000000e+000;
486 m[117] = 4.088275e-024;
487 m[119] = 1.828367e-024;
488 m[217] = 4.088275e-024;
489 m[-217] = 4.088275e-024;
490 m[219] = 1.828367e-024;
491 m[-219] = 1.828367e-024;
492 m[225] = 3.555982e-024;
493 m[227] = 3.917930e-024;
494 m[229] = 3.392846e-024;
495 m[311] = 1.000000e+016;
496 m[-311] = 1.000000e+016;
497 m[317] = 4.139699e-024;
498 m[-317] = 4.139699e-024;
499 m[319] = 3.324304e-024;
500 m[-319] = 3.324304e-024;
501 m[327] = 4.139699e-024;
502 m[-327] = 4.139699e-024;
503 m[329] = 3.324304e-024;
504 m[-329] = 3.324304e-024;
505 m[335] = 8.660687e-024;
506 m[337] = 7.565657e-024;
507 m[543] = 0.000000e+000;
508 m[-543] = 0.000000e+000;
509 m[545] = 0.000000e+000;
510 m[-545] = 0.000000e+000;
511 m[551] = 0.000000e+000;
512 m[553] = 1.253738e-020;
513 m[555] = 1.000000e+016;
514 m[557] = 0.000000e+000;
515 m[-450000000] = 0.000000e+000;
516 m[-490000000] = 0.000000e+000;
517 m[-460000000] = 0.000000e+000;
518 m[-470000000] = 0.000000e+000;
519 m[1103] = 0.000000e+000;
520 m[-1103] = 0.000000e+000;
521 m[1112] = 4.388081e-024;
522 m[-1112] = 4.388081e-024;
523 m[1116] = 1.880606e-024;
524 m[-1116] = 1.880606e-024;
525 m[1118] = 2.194041e-024;
526 m[-1118] = 2.194041e-024;
527 m[1212] = 4.388081e-024;
528 m[-1212] = 4.388081e-024;
529 m[1214] = 5.485102e-024;
530 m[-1214] = 5.485102e-024;
531 m[1216] = 1.880606e-024;
532 m[-1216] = 1.880606e-024;
533 m[1218] = 1.462694e-024;
534 m[-1218] = 1.462694e-024;
535 m[2101] = 0.000000e+000;
536 m[-2101] = 0.000000e+000;
537 m[2103] = 0.000000e+000;
538 m[-2103] = 0.000000e+000;
539 m[2116] = 4.388081e-024;
540 m[-2116] = 4.388081e-024;
541 m[2118] = 2.194041e-024;
542 m[-2118] = 2.194041e-024;
543 m[2122] = 4.388081e-024;
544 m[-2122] = 4.388081e-024;
545 m[2124] = 5.485102e-024;
546 m[-2124] = 5.485102e-024;
547 m[2126] = 1.880606e-024;
548 m[-2126] = 1.880606e-024;
549 m[2128] = 1.462694e-024;
550 m[-2128] = 1.462694e-024;
551 m[2203] = 0.000000e+000;
552 m[-2203] = 0.000000e+000;
553 m[2216] = 4.388081e-024;
554 m[-2216] = 4.388081e-024;
555 m[2218] = 2.194041e-024;
556 m[-2218] = 2.194041e-024;
557 m[2222] = 4.388081e-024;
558 m[-2222] = 4.388081e-024;
559 m[2226] = 1.880606e-024;
560 m[-2226] = 1.880606e-024;
561 m[2228] = 2.194041e-024;
562 m[-2228] = 2.194041e-024;
563 m[3101] = 0.000000e+000;
564 m[-3101] = 0.000000e+000;
565 m[3103] = 0.000000e+000;
566 m[-3103] = 0.000000e+000;
567 m[3114] = 1.670589e-023;
568 m[-3114] = 1.670589e-023;
569 m[3116] = 5.485102e-024;
570 m[-3116] = 5.485102e-024;
571 m[3118] = 3.656734e-024;
572 m[-3118] = 3.656734e-024;
573 m[3124] = 4.219309e-023;
574 m[-3124] = 4.219309e-023;
575 m[3126] = 8.227653e-024;
576 m[-3126] = 8.227653e-024;
577 m[3128] = 3.291061e-024;
578 m[-3128] = 3.291061e-024;
579 m[3201] = 0.000000e+000;
580 m[-3201] = 0.000000e+000;
581 m[3203] = 0.000000e+000;
582 m[-3203] = 0.000000e+000;
583 m[3214] = 1.828367e-023;
584 m[-3214] = 1.828367e-023;
585 m[3216] = 5.485102e-024;
586 m[-3216] = 5.485102e-024;
587 m[3218] = 3.656734e-024;
588 m[-3218] = 3.656734e-024;
589 m[3224] = 1.838582e-023;
590 m[-3224] = 1.838582e-023;
591 m[3226] = 5.485102e-024;
592 m[-3226] = 5.485102e-024;
593 m[3228] = 3.656734e-024;
594 m[-3228] = 3.656734e-024;
595 m[3303] = 0.000000e+000;
596 m[-3303] = 0.000000e+000;
597 m[3314] = 6.648608e-023;
598 m[-3314] = 6.648608e-023;
599 m[3324] = 7.233101e-023;
600 m[-3324] = 7.233101e-023;
601 m[4101] = 0.000000e+000;
602 m[-4101] = 0.000000e+000;
603 m[4103] = 0.000000e+000;
604 m[-4103] = 0.000000e+000;
605 m[4114] = 0.000000e+000;
606 m[-4114] = 0.000000e+000;
607 m[4201] = 0.000000e+000;
608 m[-4201] = 0.000000e+000;
609 m[4203] = 0.000000e+000;
610 m[-4203] = 0.000000e+000;
611 m[4214] = 3.291061e-022;
612 m[-4214] = 3.291061e-022;
613 m[4224] = 0.000000e+000;
614 m[-4224] = 0.000000e+000;
615 m[4301] = 0.000000e+000;
616 m[-4301] = 0.000000e+000;
617 m[4303] = 0.000000e+000;
618 m[-4303] = 0.000000e+000;
619 m[4314] = 0.000000e+000;
620 m[-4314] = 0.000000e+000;
621 m[4324] = 0.000000e+000;
622 m[-4324] = 0.000000e+000;
623 m[4334] = 0.000000e+000;
624 m[-4334] = 0.000000e+000;
625 m[4403] = 0.000000e+000;
626 m[-4403] = 0.000000e+000;
627 m[4412] = 3.335641e-013;
628 m[-4412] = 3.335641e-013;
629 m[4414] = 3.335641e-013;
630 m[-4414] = 3.335641e-013;
631 m[4422] = 3.335641e-013;
632 m[-4422] = 3.335641e-013;
633 m[4424] = 3.335641e-013;
634 m[-4424] = 3.335641e-013;
635 m[4432] = 3.335641e-013;
636 m[-4432] = 3.335641e-013;
637 m[4434] = 3.335641e-013;
638 m[-4434] = 3.335641e-013;
639 m[4444] = 3.335641e-013;
640 m[-4444] = 3.335641e-013;
641 m[5101] = 0.000000e+000;
642 m[-5101] = 0.000000e+000;
643 m[5103] = 0.000000e+000;
644 m[-5103] = 0.000000e+000;
645 m[5114] = 0.000000e+000;
646 m[-5114] = 0.000000e+000;
647 m[5142] = 1.290893e-012;
648 m[-5142] = 1.290893e-012;
649 m[5201] = 0.000000e+000;
650 m[-5201] = 0.000000e+000;
651 m[5203] = 0.000000e+000;
652 m[-5203] = 0.000000e+000;
653 m[5214] = 0.000000e+000;
654 m[-5214] = 0.000000e+000;
655 m[5224] = 0.000000e+000;
656 m[-5224] = 0.000000e+000;
657 m[5242] = 1.290893e-012;
658 m[-5242] = 1.290893e-012;
659 m[5301] = 0.000000e+000;
660 m[-5301] = 0.000000e+000;
661 m[5303] = 0.000000e+000;
662 m[-5303] = 0.000000e+000;
663 m[5314] = 0.000000e+000;
664 m[-5314] = 0.000000e+000;
665 m[5324] = 0.000000e+000;
666 m[-5324] = 0.000000e+000;
667 m[5334] = 0.000000e+000;
668 m[-5334] = 0.000000e+000;
669 m[5342] = 1.290893e-012;
670 m[-5342] = 1.290893e-012;
671 m[5401] = 0.000000e+000;
672 m[-5401] = 0.000000e+000;
673 m[5403] = 0.000000e+000;
674 m[-5403] = 0.000000e+000;
675 m[5412] = 1.290893e-012;
676 m[-5412] = 1.290893e-012;
677 m[5414] = 1.290893e-012;
678 m[-5414] = 1.290893e-012;
679 m[5422] = 1.290893e-012;
680 m[-5422] = 1.290893e-012;
681 m[5424] = 1.290893e-012;
682 m[-5424] = 1.290893e-012;
683 m[5432] = 1.290893e-012;
684 m[-5432] = 1.290893e-012;
685 m[5434] = 1.290893e-012;
686 m[-5434] = 1.290893e-012;
687 m[5442] = 1.290893e-012;
688 m[-5442] = 1.290893e-012;
689 m[5444] = 1.290893e-012;
690 m[-5444] = 1.290893e-012;
691 m[5503] = 0.000000e+000;
692 m[-5503] = 0.000000e+000;
693 m[5512] = 1.290893e-012;
694 m[-5512] = 1.290893e-012;
695 m[5514] = 1.290893e-012;
696 m[-5514] = 1.290893e-012;
697 m[5522] = 1.290893e-012;
698 m[-5522] = 1.290893e-012;
699 m[5524] = 1.290893e-012;
700 m[-5524] = 1.290893e-012;
701 m[5532] = 1.290893e-012;
702 m[-5532] = 1.290893e-012;
703 m[5534] = 1.290893e-012;
704 m[-5534] = 1.290893e-012;
705 m[5542] = 1.290893e-012;
706 m[-5542] = 1.290893e-012;
707 m[5544] = 1.290893e-012;
708 m[-5544] = 1.290893e-012;
709 m[5554] = 1.290893e-012;
710 m[-5554] = 1.290893e-012;
711 m[10022] = 0.000000e+000;
712 m[10111] = 2.483820e-024;
713 m[10113] = 4.635297e-024;
714 m[10115] = 2.541360e-024;
715 m[10211] = 2.483820e-024;
716 m[-10211] = 2.483820e-024;
717 m[10213] = 4.635297e-024;
718 m[-10213] = 4.635297e-024;
719 m[10215] = 2.541360e-024;
720 m[-10215] = 2.541360e-024;
721 m[9010221] = 1.316424e-023;
722 m[10223] = 1.828367e-024;
723 m[10225] = 0.000000e+000;
724 m[10315] = 3.538775e-024;
725 m[-10315] = 3.538775e-024;
726 m[10325] = 3.538775e-024;
727 m[-10325] = 3.538775e-024;
728 m[10331] = 5.265698e-024;
729 m[10333] = 0.000000e+000;
730 m[10335] = 0.000000e+000;
731 m[10443] = 0.000000e+000;
732 m[10541] = 0.000000e+000;
733 m[-10541] = 0.000000e+000;
734 m[10543] = 0.000000e+000;
735 m[-10543] = 0.000000e+000;
736 m[10551] = 1.000000e+016;
737 m[10553] = 0.000000e+000;
738 m[10555] = 0.000000e+000;
739 m[11112] = 0.000000e+000;
740 m[-11112] = 0.000000e+000;
741 m[11114] = 2.194041e-024;
742 m[-11114] = 2.194041e-024;
743 m[11116] = 1.880606e-024;
744 m[-11116] = 1.880606e-024;
745 m[11212] = 1.880606e-024;
746 m[-11212] = 1.880606e-024;
747 m[11216] = 0.000000e+000;
748 m[-11216] = 0.000000e+000;
749 m[12112] = 1.880606e-024;
750 m[-12112] = 1.880606e-024;
751 m[12114] = 2.194041e-024;
752 m[-12114] = 2.194041e-024;
753 m[12116] = 5.063171e-024;
754 m[-12116] = 5.063171e-024;
755 m[12118] = 0.000000e+000;
756 m[-12118] = 0.000000e+000;
757 m[12122] = 0.000000e+000;
758 m[-12122] = 0.000000e+000;
759 m[12126] = 1.880606e-024;
760 m[-12126] = 1.880606e-024;
761 m[12212] = 1.880606e-024;
762 m[-12212] = 1.880606e-024;
763 m[12214] = 2.194041e-024;
764 m[-12214] = 2.194041e-024;
765 m[12216] = 5.063171e-024;
766 m[-12216] = 5.063171e-024;
767 m[12218] = 0.000000e+000;
768 m[-12218] = 0.000000e+000;
769 m[12222] = 0.000000e+000;
770 m[-12222] = 0.000000e+000;
771 m[12224] = 2.194041e-024;
772 m[-12224] = 2.194041e-024;
773 m[12226] = 1.880606e-024;
774 m[-12226] = 1.880606e-024;
775 m[13112] = 6.582122e-024;
776 m[-13112] = 6.582122e-024;
777 m[13114] = 1.097020e-023;
778 m[-13114] = 1.097020e-023;
779 m[13116] = 5.485102e-024;
780 m[-13116] = 5.485102e-024;
781 m[13122] = 1.316424e-023;
782 m[-13122] = 1.316424e-023;
783 m[13124] = 1.097020e-023;
784 m[-13124] = 1.097020e-023;
785 m[13126] = 6.928549e-024;
786 m[-13126] = 6.928549e-024;
787 m[13212] = 6.582122e-024;
788 m[-13212] = 6.582122e-024;
789 m[13214] = 1.097020e-023;
790 m[-13214] = 1.097020e-023;
791 m[13216] = 5.485102e-024;
792 m[-13216] = 5.485102e-024;
793 m[13222] = 6.582122e-024;
794 m[-13222] = 6.582122e-024;
795 m[13224] = 1.097020e-023;
796 m[-13224] = 1.097020e-023;
797 m[13226] = 5.485102e-024;
798 m[-13226] = 5.485102e-024;
799 m[13314] = 2.742551e-023;
800 m[-13314] = 2.742551e-023;
801 m[13316] = 0.000000e+000;
802 m[-13316] = 0.000000e+000;
803 m[13324] = 2.742551e-023;
804 m[-13324] = 2.742551e-023;
805 m[13326] = 0.000000e+000;
806 m[-13326] = 0.000000e+000;
807 m[14122] = 1.828367e-022;
808 m[-14122] = 1.828367e-022;
809 m[14124] = 0.000000e+000;
810 m[-14124] = 0.000000e+000;
811 m[10221] = 2.194040e-024;
812 m[20223] = 2.742551e-023;
813 m[20315] = 2.384827e-024;
814 m[-20315] = 2.384827e-024;
815 m[20325] = 2.384827e-024;
816 m[-20325] = 2.384827e-024;
817 m[20333] = 1.185968e-023;
818 m[20543] = 0.000000e+000;
819 m[-20543] = 0.000000e+000;
820 m[20553] = 1.000000e+016;
821 m[20555] = 0.000000e+000;
822 m[21112] = 2.632849e-024;
823 m[-21112] = 2.632849e-024;
824 m[21114] = 3.291061e-024;
825 m[-21114] = 3.291061e-024;
826 m[21212] = 2.632849e-024;
827 m[-21212] = 2.632849e-024;
828 m[21214] = 6.582122e-024;
829 m[-21214] = 6.582122e-024;
830 m[22112] = 4.388081e-024;
831 m[-22112] = 4.388081e-024;
832 m[22114] = 3.291061e-024;
833 m[-22114] = 3.291061e-024;
834 m[22122] = 2.632849e-024;
835 m[-22122] = 2.632849e-024;
836 m[22124] = 6.582122e-024;
837 m[-22124] = 6.582122e-024;
838 m[22212] = 4.388081e-024;
839 m[-22212] = 4.388081e-024;
840 m[22214] = 3.291061e-024;
841 m[-22214] = 3.291061e-024;
842 m[22222] = 2.632849e-024;
843 m[-22222] = 2.632849e-024;
844 m[22224] = 3.291061e-024;
845 m[-22224] = 3.291061e-024;
846 m[23112] = 7.313469e-024;
847 m[-23112] = 7.313469e-024;
848 m[23114] = 2.991874e-024;
849 m[-23114] = 2.991874e-024;
850 m[23122] = 4.388081e-024;
851 m[-23122] = 4.388081e-024;
852 m[23124] = 6.582122e-024;
853 m[-23124] = 6.582122e-024;
854 m[23126] = 3.291061e-024;
855 m[-23126] = 3.291061e-024;
856 m[23212] = 7.313469e-024;
857 m[-23212] = 7.313469e-024;
858 m[23214] = 2.991874e-024;
859 m[-23214] = 2.991874e-024;
860 m[23222] = 7.313469e-024;
861 m[-23222] = 7.313469e-024;
862 m[23224] = 2.991874e-024;
863 m[-23224] = 2.991874e-024;
864 m[23314] = 0.000000e+000;
865 m[-23314] = 0.000000e+000;
866 m[23324] = 0.000000e+000;
867 m[-23324] = 0.000000e+000;
868 m[30113] = 2.742551e-024;
869 m[30213] = 2.742551e-024;
870 m[-30213] = 2.742551e-024;
871 m[30223] = 2.991874e-024;
872 m[30313] = 2.056913e-024;
873 m[-30313] = 2.056913e-024;
874 m[30323] = 2.056913e-024;
875 m[-30323] = 2.056913e-024;
876 m[30343] = 0.000000e+000;
877 m[-30343] = 0.000000e+000;
878 m[30353] = 0.000000e+000;
879 m[-30353] = 0.000000e+000;
880 m[30363] = 0.000000e+000;
881 m[-30363] = 0.000000e+000;
882 m[30411] = 0.000000e+000;
883 m[-30411] = 0.000000e+000;
884 m[30413] = 0.000000e+000;
885 m[-30413] = 0.000000e+000;
886 m[30421] = 0.000000e+000;
887 m[-30421] = 0.000000e+000;
888 m[30423] = 0.000000e+000;
889 m[-30423] = 0.000000e+000;
890 m[30443] = 2.789035e-023;
891 m[30553] = 0.000000e+000;
892 m[31114] = 1.880606e-024;
893 m[-31114] = 1.880606e-024;
894 m[31214] = 4.388081e-024;
895 m[-31214] = 4.388081e-024;
896 m[32112] = 4.388081e-024;
897 m[-32112] = 4.388081e-024;
898 m[32114] = 1.880606e-024;
899 m[-32114] = 1.880606e-024;
900 m[32124] = 4.388081e-024;
901 m[-32124] = 4.388081e-024;
902 m[32212] = 4.388081e-024;
903 m[-32212] = 4.388081e-024;
904 m[32214] = 1.880606e-024;
905 m[-32214] = 1.880606e-024;
906 m[32224] = 1.880606e-024;
907 m[-32224] = 1.880606e-024;
908 m[33122] = 1.880606e-023;
909 m[-33122] = 1.880606e-023;
910 m[33314] = 0.000000e+000;
911 m[-33314] = 0.000000e+000;
912 m[33324] = 0.000000e+000;
913 m[-33324] = 0.000000e+000;
914 m[41214] = 0.000000e+000;
915 m[-41214] = 0.000000e+000;
916 m[42112] = 6.582122e-024;
917 m[-42112] = 6.582122e-024;
918 m[42124] = 0.000000e+000;
919 m[-42124] = 0.000000e+000;
920 m[42212] = 6.582122e-024;
921 m[-42212] = 6.582122e-024;
922 m[43122] = 2.194041e-024;
923 m[-43122] = 2.194041e-024;
924 m[52114] = 0.000000e+000;
925 m[-52114] = 0.000000e+000;
926 m[52214] = 0.000000e+000;
927 m[-52214] = 0.000000e+000;
928 m[53122] = 4.388081e-024;
929 m[-53122] = 4.388081e-024;
930 m[100111] = 1.645531e-024;
931 m[100113] = 2.123265e-024;
932 m[100211] = 1.645531e-024;
933 m[-100211] = 1.645531e-024;
934 m[100213] = 2.123265e-024;
935 m[-100213] = 2.123265e-024;
936 m[100221] = 1.196749e-023;
937 m[100223] = 3.871836e-024;
938 m[100225] = 0.000000e+000;
939 m[100311] = 0.000000e+000;
940 m[-100311] = 0.000000e+000;
941 m[100313] = 2.837122e-024;
942 m[-100313] = 2.837122e-024;
943 m[100315] = 0.000000e+000;
944 m[-100315] = 0.000000e+000;
945 m[100321] = 0.000000e+000;
946 m[-100321] = 0.000000e+000;
947 m[100323] = 2.837122e-024;
948 m[-100323] = 2.837122e-024;
949 m[100325] = 0.000000e+000;
950 m[-100325] = 0.000000e+000;
951 m[100331] = 0.000000e+000;
952 m[100333] = 4.388081e-024;
953 m[100335] = 3.291061e-024;
954 m[100441] = 0.000000e+000;
955 m[100551] = 0.000000e+000;
956 m[100553] = 1.495937e-020;
957 m[100555] = 1.000000e+016;
958 m[100557] = 0.000000e+000;
959 m[110551] = 1.000000e+016;
960 m[110553] = 0.000000e+000;
961 m[110555] = 0.000000e+000;
962 m[120553] = 1.000000e+016;
963 m[120555] = 0.000000e+000;
964 m[130553] = 0.000000e+000;
965 m[200111] = 3.134344e-024;
966 m[200211] = 3.134344e-024;
967 m[-200211] = 3.134344e-024;
968 m[200551] = 0.000000e+000;
969 m[200553] = 2.502708e-020;
970 m[200555] = 0.000000e+000;
971 m[210551] = 0.000000e+000;
972 m[210553] = 0.000000e+000;
973 m[220553] = 0.000000e+000;
974 m[300553] = 4.701516e-023;
975 m[9000221] = 0.000000e+000;
976 m[9000443] = 1.265793e-023;
977 m[9000553] = 5.983747e-024;
978 m[9010443] = 8.438618e-024;
979 m[9010553] = 8.331800e-024;
980 m[9020221] = 6.038644e-024;
981 m[9020443] = 1.530726e-023;
982 m[9060225] = 4.388081e-024;
983 m[9070225] = 2.056913e-024;
984 m[1000001] = 0.000000e+000;
985 m[-1000001] = 0.000000e+000;
986 m[1000002] = 0.000000e+000;
987 m[-1000002] = 0.000000e+000;
988 m[1000003] = 0.000000e+000;
989 m[-1000003] = 0.000000e+000;
990 m[1000004] = 0.000000e+000;
991 m[-1000004] = 0.000000e+000;
992 m[1000005] = 0.000000e+000;
993 m[-1000005] = 0.000000e+000;
994 m[1000006] = 0.000000e+000;
995 m[-1000006] = 0.000000e+000;
996 m[1000011] = 0.000000e+000;
997 m[-1000011] = 0.000000e+000;
998 m[1000012] = 0.000000e+000;
999 m[-1000012] = 0.000000e+000;
1000 m[1000013] = 0.000000e+000;
1001 m[-1000013] = 0.000000e+000;
1002 m[1000014] = 0.000000e+000;
1003 m[-1000014] = 0.000000e+000;
1004 m[1000015] = 0.000000e+000;
1005 m[-1000015] = 0.000000e+000;
1006 m[1000016] = 0.000000e+000;
1007 m[-1000016] = 0.000000e+000;
1008 m[1000021] = 0.000000e+000;
1009 m[1000022] = 0.000000e+000;
1010 m[1000023] = 0.000000e+000;
1011 m[1000024] = 0.000000e+000;
1012 m[-1000024] = 0.000000e+000;
1013 m[1000025] = 0.000000e+000;
1014 m[1000035] = 0.000000e+000;
1015 m[1000037] = 0.000000e+000;
1016 m[-1000037] = 0.000000e+000;
1017 m[1000039] = 0.000000e+000;
1018 m[2000001] = 0.000000e+000;
1019 m[-2000001] = 0.000000e+000;
1020 m[2000002] = 0.000000e+000;
1021 m[-2000002] = 0.000000e+000;
1022 m[2000003] = 0.000000e+000;
1023 m[-2000003] = 0.000000e+000;
1024 m[2000004] = 0.000000e+000;
1025 m[-2000004] = 0.000000e+000;
1026 m[2000005] = 0.000000e+000;
1027 m[-2000005] = 0.000000e+000;
1028 m[2000006] = 0.000000e+000;
1029 m[-2000006] = 0.000000e+000;
1030 m[2000011] = 0.000000e+000;
1031 m[-2000011] = 0.000000e+000;
1032 m[2000012] = 0.000000e+000;
1033 m[-2000012] = 0.000000e+000;
1034 m[2000013] = 0.000000e+000;
1035 m[-2000013] = 0.000000e+000;
1036 m[2000014] = 0.000000e+000;
1037 m[-2000014] = 0.000000e+000;
1038 m[2000015] = 0.000000e+000;
1039 m[-2000015] = 0.000000e+000;
1040 m[2000016] = 0.000000e+000;
1041 m[-2000016] = 0.000000e+000;
1042 m[3000111] = 0.000000e+000;
1043 m[3000113] = 0.000000e+000;
1044 m[3000211] = 0.000000e+000;
1045 m[-3000211] = 0.000000e+000;
1046 m[3000213] = 0.000000e+000;
1047 m[-3000213] = 0.000000e+000;
1048 m[3000221] = 0.000000e+000;
1049 m[3000223] = 0.000000e+000;
1050 m[3000331] = 0.000000e+000;
1051 m[3100021] = 0.000000e+000;
1052 m[3100111] = 0.000000e+000;
1053 m[3100113] = 0.000000e+000;
1054 m[3200111] = 0.000000e+000;
1055 m[3200113] = 0.000000e+000;
1056 m[3300113] = 0.000000e+000;
1057 m[3400113] = 0.000000e+000;
1058 m[4000001] = 0.000000e+000;
1059 m[-4000001] = 0.000000e+000;
1060 m[4000002] = 0.000000e+000;
1061 m[-4000002] = 0.000000e+000;
1062 m[4000011] = 0.000000e+000;
1063 m[-4000011] = 0.000000e+000;
1064 m[4000012] = 0.000000e+000;
1065 m[-4000012] = 0.000000e+000;
1066 m[5000039] = 0.000000e+000;
1067 m[9900012] = 0.000000e+000;
1068 m[9900014] = 0.000000e+000;
1069 m[9900016] = 0.000000e+000;
1070 m[9900023] = 0.000000e+000;
1071 m[9900024] = 0.000000e+000;
1072 m[-9900024] = 0.000000e+000;
1073 m[9900041] = 0.000000e+000;
1074 m[-9900041] = 0.000000e+000;
1075 m[9900042] = 0.000000e+000;
1076 m[-9900042] = 0.000000e+000;
1077 m[1027013000] = 0.000000e+000;
1078 m[1012006000] = 0.000000e+000;
1079 m[1063029000] = 0.000000e+000;
1080 m[1014007000] = 0.000000e+000;
1081 m[1016008000] = 0.000000e+000;
1082 m[1028014000] = 0.000000e+000;
1083 m[1065029000] = 0.000000e+000;
1084 m[1009004000] = 0.000000e+000;
1085 m[1019009000] = 0.000000e+000;
1086 m[1056026000] = 0.000000e+000;
1087 m[1207082000] = 0.000000e+000;
1088 m[1208082000] = 0.000000e+000;
1089 m[1029014000] = 0.000000e+000;
1090 m[1206082000] = 0.000000e+000;
1091 m[1054026000] = 0.000000e+000;
1092 m[1018008000] = 0.000000e+000;
1093 m[1030014000] = 0.000000e+000;
1094 m[1057026000] = 0.000000e+000;
1095 m[1204082000] = 0.000000e+000;
1096 m[-99000000] = 0.000000e+000;
1097 m[1028013000] = 0.000000e+000;
1098 m[1040018000] = 0.000000e+000;
1099 m[1011005000] = 0.000000e+000;
1100 m[1012005000] = 0.000000e+000;
1101 m[1013006000] = 0.000000e+000;
1102 m[1014006000] = 0.000000e+000;
1103 m[1052024000] = 0.000000e+000;
1104 m[1024012000] = 0.000000e+000;
1105 m[1026012000] = 0.000000e+000;
1106 m[1027012000] = 0.000000e+000;
1107 m[1015007000] = 0.000000e+000;
1108 m[1022010000] = 0.000000e+000;
1109 m[1058028000] = 0.000000e+000;
1110 m[1060028000] = 0.000000e+000;
1111 m[1062028000] = 0.000000e+000;
1112 m[1064028000] = 0.000000e+000;
1113 m[1007003000] = 0.000000e+000;
1114 m[1025012000] = 0.000000e+000;
1115 m[1053024000] = 0.000000e+000;
1116 m[1055025000] = 0.000000e+000;
1117 m[1008004000] = 0.000000e+000;
1118 m[1010004000] = 0.000000e+000;
1119 m[1010005000] = 0.000000e+000;
1120 m[1016007000] = 0.000000e+000;
1121 m[1017008000] = 0.000000e+000;
1122 m[1019008000] = 0.000000e+000;
1123 m[1023010000] = 0.000000e+000;
1124 m[1024011000] = 0.000000e+000;
1125 m[1031015000] = 0.000000e+000;
1126 m[1039017000] = 0.000000e+000;
1127 m[1040017000] = 0.000000e+000;
1128 m[1036018000] = 0.000000e+000;
1129 m[1050024000] = 0.000000e+000;
1130 m[1054024000] = 0.000000e+000;
1131 m[1059026000] = 0.000000e+000;
1132 m[1061028000] = 0.000000e+000;
1133 m[1063028000] = 0.000000e+000;
1134 m[1092042000] = 0.000000e+000;
1135 m[1095042000] = 0.000000e+000;
1136 m[1096042000] = 0.000000e+000;
1137 m[1097042000] = 0.000000e+000;
1138 m[1098042000] = 0.000000e+000;
1139 m[1100042000] = 0.000000e+000;
1140 m[1108046000] = 0.000000e+000;
1141
1142 // Added by hand:
1143 m[9902210] = 0.000000e+000; //diffractive p-state -> assume no lifetime
1144 return true;
1145 }
1146
1147
1148 private:
1149
1150 /// @name Histograms
1151 /// @{
1152 Histo1DPtr _h_mult_total; // full kinematic range
1153 Histo1DPtr _h_mult_eta[5]; // in eta bins
1154 Histo1DPtr _h_mult_pt[5]; // in pT bins
1155 Histo1DPtr _h_dndeta; // density dn/deta
1156 Histo1DPtr _h_dndpt; // density dn/dpT
1157 /// @}
1158
1159
1160 /// @name Private variables
1161 double _p_min;
1162 double _pt_min;
1163 double _eta_min;
1164 double _eta_max;
1165 double _maxlft;
1166
1167 /// Count selected events
1168 CounterPtr _sumW;
1169
1170 map<int, double> _partLftMap; // Map <PDGID, PDGLIFETIME>
1171
1172 };
1173
1174
1175 RIVET_DECLARE_PLUGIN(LHCB_2014_I1281685);
1176
1177}
|