LHCB_2014_I1281685.cc
Go to the documentation of this file.
00001 // -*- C++ -*- 00002 #include "Rivet/Analysis.hh" 00003 #include "Rivet/Projections/ChargedFinalState.hh" 00004 00005 namespace Rivet { 00006 00007 00008 /// Charged particle multiplicities and densities in $pp$ collisions at $\sqrt{s} = 7$ TeV 00009 class LHCB_2014_I1281685 : public Analysis { 00010 public: 00011 00012 /// @name Constructors etc. 00013 //@{ 00014 00015 /// Constructor 00016 LHCB_2014_I1281685() 00017 : Analysis("LHCB_2014_I1281685"), 00018 _p_min(2.0), 00019 _pt_min(0.2), 00020 _eta_min(2.0), 00021 _eta_max(4.8), 00022 _maxlft(1.0e-11) 00023 { } 00024 00025 //@} 00026 00027 00028 /// @name Analysis methods 00029 //@{ 00030 00031 /// Book histograms and initialise projections before the run 00032 void init() { 00033 fillMap(_partLftMap); 00034 00035 // Projections 00036 addProjection(ChargedFinalState(_eta_min, _eta_max, _pt_min*GeV), "CFS"); 00037 00038 // Book histograms 00039 _h_mult_total = bookHisto1D("d03-x01-y01", 50, 0.5, 50.5); 00040 00041 _h_mult_eta[0] = bookHisto1D("d04-x01-y01", 21, -0.5, 20.5); //eta=[2.0,2.5] 00042 _h_mult_eta[1] = bookHisto1D("d04-x01-y02", 21, -0.5, 20.5); //eta=[2.5,3.0] 00043 _h_mult_eta[2] = bookHisto1D("d04-x01-y03", 21, -0.5, 20.5); //eta=[3.0,3.5] 00044 _h_mult_eta[3] = bookHisto1D("d04-x01-y04", 21, -0.5, 20.5); //eta=[3.5,4.0] 00045 _h_mult_eta[4] = bookHisto1D("d04-x01-y05", 21, -0.5, 20.5); //eta=[4.0,4.5] 00046 00047 _h_mult_pt[0] = bookHisto1D("d05-x01-y01", 21, -0.5, 20.5); //pT=[0.2,0.3]GeV 00048 _h_mult_pt[1] = bookHisto1D("d05-x01-y02", 21, -0.5, 20.5); //pT=[0.3,0.4]GeV 00049 _h_mult_pt[2] = bookHisto1D("d05-x01-y03", 21, -0.5, 20.5); //pT=[0.4,0.6]GeV 00050 _h_mult_pt[3] = bookHisto1D("d05-x01-y04", 21, -0.5, 20.5); //pT=[0.6,1.0]GeV 00051 _h_mult_pt[4] = bookHisto1D("d05-x01-y05", 21, -0.5, 20.5); //pT=[1.0,2.0]GeV 00052 00053 _h_dndeta = bookHisto1D("d01-x01-y01", 14, 2.0, 4.8); //eta=[2,4.8] 00054 _h_dndpt = bookHisto1D("d02-x01-y01", 18, 0.2, 2.0); //pT =[0,2]GeV 00055 00056 // Counters 00057 _sumW = 0; 00058 } 00059 00060 00061 /// Perform the per-event analysis 00062 void analyze(const Event& event) { 00063 00064 // Variable to store multiplicities per event 00065 int LHCbcountAll = 0; //count particles fulfiling all requirements 00066 int LHCbcountEta[8] = {0,0,0,0,0,0,0,0}; //count per eta-bin 00067 int LHCbcountPt[7] = {0,0,0,0,0,0,0}; //count per pT-bin 00068 vector<double> val_dNdEta; 00069 vector<double> val_dNdPt; 00070 val_dNdEta.clear(); 00071 val_dNdPt.clear(); 00072 00073 const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(event, "CFS"); 00074 foreach (const Particle& p, cfs.particles()) { 00075 int id = p.pdgId(); 00076 // continue if particle is not a pion, kaon, proton, muon or electron 00077 if ( !( (abs(id) == 211) || (abs(id) == 321) || (abs(id) == 2212) || (abs(id) == 13) || (abs(id) == 11)) ) { 00078 continue; 00079 } 00080 00081 const FourMomentum& qmom = p.momentum(); 00082 const double eta = p.momentum().eta(); 00083 const double pT = p.momentum().pT(); 00084 //minimum momentum 00085 if (qmom.p3().mod() < _p_min) continue; 00086 //minimum tr. momentum 00087 if (pT < _pt_min) continue; 00088 //eta range 00089 if ((eta < _eta_min) || (eta > _eta_max)) continue; 00090 00091 /* Select only prompt particles via lifetime */ 00092 //Sum of all mother lifetimes (PDG lifetime) < 10ps 00093 double ancestors_sumlft = getAncestorSumLifetime(p); 00094 if( (ancestors_sumlft > _maxlft) || (ancestors_sumlft < 0) ) continue; 00095 00096 //after all cuts; 00097 LHCbcountAll++; //count particles in whole kin. range 00098 00099 //in eta bins 00100 if( eta >2.0 && eta <= 2.5) LHCbcountEta[0]++; 00101 if( eta >2.5 && eta <= 3.0) LHCbcountEta[1]++; 00102 if( eta >3.0 && eta <= 3.5) LHCbcountEta[2]++; 00103 if( eta >3.5 && eta <= 4.0) LHCbcountEta[3]++; 00104 if( eta >4.0 && eta <= 4.5) LHCbcountEta[4]++; 00105 if( eta >2.0 && eta <= 4.8) LHCbcountEta[5]++; //cross-check 00106 //in pT bins 00107 if( pT > 0.2 && pT <= 0.3) LHCbcountPt[0]++; 00108 if( pT > 0.3 && pT <= 0.4) LHCbcountPt[1]++; 00109 if( pT > 0.4 && pT <= 0.6) LHCbcountPt[2]++; 00110 if( pT > 0.6 && pT <= 1.0) LHCbcountPt[3]++; 00111 if( pT > 1.0 && pT <= 2.0) LHCbcountPt[4]++; 00112 if( pT > 0.2) LHCbcountPt[5]++; //cross-check 00113 00114 //particle densities -> need proper normalization (finalize) 00115 val_dNdPt.push_back( pT ); 00116 val_dNdEta.push_back( eta ); 00117 }//end foreach 00118 00119 00120 // Fill histograms only, if at least 1 particle pre event was within the 00121 // kinematic range of the analysis! 00122 if (LHCbcountAll) { 00123 const double weight = event.weight(); 00124 _sumW += weight; 00125 00126 _h_mult_total->fill(LHCbcountAll); 00127 00128 _h_mult_eta[0]->fill(LHCbcountEta[0], weight); 00129 _h_mult_eta[1]->fill(LHCbcountEta[1], weight); 00130 _h_mult_eta[2]->fill(LHCbcountEta[2], weight); 00131 _h_mult_eta[3]->fill(LHCbcountEta[3], weight); 00132 _h_mult_eta[4]->fill(LHCbcountEta[4], weight); 00133 00134 _h_mult_pt[0]->fill(LHCbcountPt[0], weight); 00135 _h_mult_pt[1]->fill(LHCbcountPt[1], weight); 00136 _h_mult_pt[2]->fill(LHCbcountPt[2], weight); 00137 _h_mult_pt[3]->fill(LHCbcountPt[3], weight); 00138 _h_mult_pt[4]->fill(LHCbcountPt[4], weight); 00139 00140 for (size_t part = 0; part < val_dNdEta.size(); part++) 00141 _h_dndeta->fill(val_dNdEta[part], weight); 00142 00143 for (size_t part = 0; part < val_dNdPt.size(); part++) 00144 _h_dndpt->fill(val_dNdPt[part], weight); 00145 00146 } 00147 } 00148 00149 00150 /// Normalise histograms etc., after the run 00151 void finalize() { 00152 const double scalefactor = 1.0/_sumW; // normalize multiplicity histograms by nEvents 00153 const double scale1k = 1000.; // to match '10^3' scale in reference histograms 00154 00155 scale( _h_dndeta, scalefactor ); 00156 scale( _h_dndpt, scalefactor*0.1 ); //additional factor 0.1 for [0.1 GeV/c] 00157 scale( _h_mult_total, scalefactor*scale1k); 00158 00159 _h_mult_eta[0]->scaleW( scalefactor*scale1k ); 00160 _h_mult_eta[1]->scaleW( scalefactor*scale1k ); 00161 _h_mult_eta[2]->scaleW( scalefactor*scale1k ); 00162 _h_mult_eta[3]->scaleW( scalefactor*scale1k ); 00163 _h_mult_eta[4]->scaleW( scalefactor*scale1k ); 00164 00165 _h_mult_pt[0]->scaleW( scalefactor*scale1k ); 00166 _h_mult_pt[1]->scaleW( scalefactor*scale1k ); 00167 _h_mult_pt[2]->scaleW( scalefactor*scale1k ); 00168 _h_mult_pt[3]->scaleW( scalefactor*scale1k ); 00169 _h_mult_pt[4]->scaleW( scalefactor*scale1k ); 00170 } 00171 00172 //@} 00173 00174 00175 private: 00176 00177 00178 // Get mean PDG lifetime for particle with PID 00179 double getLifetime(int pid) { 00180 double lft = 0.; 00181 map<int, double>::iterator pPartLft = _partLftMap.find(pid); 00182 if (pPartLft != _partLftMap.end()) { 00183 lft = (*pPartLft).second; 00184 } else { 00185 // allow identifying missing life times only in debug mode 00186 MSG_DEBUG("Could not determine lifetime for particle with PID " << pid << "... Assume non-prompt particle"); 00187 lft = -1; 00188 } 00189 return lft; 00190 } 00191 00192 00193 // Get sum of all ancestor particles 00194 const double getAncestorSumLifetime(const Particle& p) { 00195 double lftSum = 0.; 00196 double plft = 0.; 00197 const GenParticle* part = p.genParticle(); 00198 if ( 0 == part ) return -1; 00199 const GenVertex* ivtx = part->production_vertex(); 00200 while(ivtx) { 00201 if (ivtx->particles_in_size() < 1) { lftSum = -1.; break; }; 00202 const GenVertex::particles_in_const_iterator iPart_invtx = ivtx->particles_in_const_begin(); 00203 part = (*iPart_invtx); 00204 if ( !(part) ) { lftSum = -1.; break; }; 00205 ivtx = part->production_vertex(); 00206 if ( (part->pdg_id() == 2212) || !(ivtx) ) break; // reached beam 00207 plft = getLifetime(part->pdg_id()); 00208 if (plft < 0.) { lftSum = -1.; break; }; 00209 lftSum += plft; 00210 } 00211 return (lftSum); 00212 } 00213 00214 00215 /// Hard-coded map linking PDG ID with PDG lifetime[s] (converted from ParticleTable.txt) 00216 bool fillMap(map<int, double>& m) { 00217 // PDGID = LIFETIME 00218 m[22] = 1.000000e+016; 00219 m[-11] = 1.000000e+016; 00220 m[11] = 1.000000e+016; 00221 m[12] = 1.000000e+016; 00222 m[-13] = 2.197036e-006; 00223 m[13] = 2.197036e-006; 00224 m[111] = 8.438618e-017; 00225 m[211] = 2.603276e-008; 00226 m[-211] = 2.603276e-008; 00227 m[130] = 5.174624e-008; 00228 m[321] = 1.238405e-008; 00229 m[-321] = 1.238405e-008; 00230 m[2112] = 885.646128; 00231 m[2212] = 1.000000e+016; 00232 m[-2212] = 1.000000e+016; 00233 m[310] = 8.934603e-011; 00234 m[221] = 5.578070e-019; 00235 m[3122] = 2.631796e-010; 00236 m[3222] = 8.018178e-011; 00237 m[3212] = 7.395643e-020; 00238 m[3112] = 1.479129e-010; 00239 m[3322] = 2.899613e-010; 00240 m[3312] = 1.637344e-010; 00241 m[3334] = 8.207135e-011; 00242 m[-2112] = 885.646128; 00243 m[-3122] = 2.631796e-010; 00244 m[-3222] = 8.018178e-011; 00245 m[-3212] = 7.395643e-020; 00246 m[-3112] = 1.479129e-010; 00247 m[-3322] = 2.899613e-010; 00248 m[-3312] = 1.637344e-010; 00249 m[-3334] = 8.207135e-011; 00250 m[113] = 4.411610e-024; 00251 m[213] = 4.411610e-024; 00252 m[-213] = 4.411610e-024; 00253 m[223] = 7.798723e-023; 00254 m[333] = 1.545099e-022; 00255 m[323] = 1.295693e-023; 00256 m[-323] = 1.295693e-023; 00257 m[313] = 1.298249e-023; 00258 m[-313] = 1.298249e-023; 00259 m[20213] = 1.500000e-024; 00260 m[-20213] = 1.500000e-024; 00261 m[450000000] = 1.000000e+015; 00262 m[460000000] = 1.000000e+015; 00263 m[470000000] = 1.000000e+015; 00264 m[480000000] = 1.000000e+015; 00265 m[490000000] = 1.000000e+015; 00266 m[20022] = 1.000000e+016; 00267 m[-15] = 2.906014e-013; 00268 m[15] = 2.906014e-013; 00269 m[24] = 3.104775e-025; 00270 m[-24] = 3.104775e-025; 00271 m[23] = 2.637914e-025; 00272 m[411] = 1.051457e-012; 00273 m[-411] = 1.051457e-012; 00274 m[421] = 4.116399e-013; 00275 m[-421] = 4.116399e-013; 00276 m[431] = 4.904711e-013; 00277 m[-431] = 4.904711e-013; 00278 m[4122] = 1.994582e-013; 00279 m[-4122] = 1.994582e-013; 00280 m[443] = 7.565657e-021; 00281 m[413] = 6.856377e-021; 00282 m[-413] = 6.856377e-021; 00283 m[423] = 1.000003e-019; 00284 m[-423] = 1.000003e-019; 00285 m[433] = 1.000003e-019; 00286 m[-433] = 1.000003e-019; 00287 m[521] = 1.671000e-012; 00288 m[-521] = 1.671000e-012; 00289 m[511] = 1.536000e-012; 00290 m[-511] = 1.536000e-012; 00291 m[531] = 1.461000e-012; 00292 m[-531] = 1.461000e-012; 00293 m[541] = 4.600000e-013; 00294 m[-541] = 4.600000e-013; 00295 m[5122] = 1.229000e-012; 00296 m[-5122] = 1.229000e-012; 00297 m[4112] = 4.388081e-022; 00298 m[-4112] = 4.388081e-022; 00299 m[4212] = 3.999999e-022; 00300 m[-4212] = 3.999999e-022; 00301 m[4222] = 3.291060e-022; 00302 m[-4222] = 3.291060e-022; 00303 m[25] = 9.400000e-026; 00304 m[35] = 9.400000e-026; 00305 m[36] = 9.400000e-026; 00306 m[37] = 9.400000e-026; 00307 m[-37] = 9.400000e-026; 00308 m[4312] = 9.800002e-014; 00309 m[-4312] = 9.800002e-014; 00310 m[4322] = 3.500001e-013; 00311 m[-4322] = 3.500001e-013; 00312 m[4332] = 6.453061e-014; 00313 m[-4332] = 6.453061e-014; 00314 m[4132] = 9.824063e-014; 00315 m[-4132] = 9.824063e-014; 00316 m[4232] = 4.417532e-013; 00317 m[-4232] = 4.417532e-013; 00318 m[5222] = 1.000000e-019; 00319 m[-5222] = 1.000000e-019; 00320 m[5212] = 1.000000e-019; 00321 m[-5212] = 1.000000e-019; 00322 m[5112] = 1.000000e-019; 00323 m[-5112] = 1.000000e-019; 00324 m[5312] = 1.000000e-019; 00325 m[-5312] = 1.000000e-019; 00326 m[5322] = 1.000000e-019; 00327 m[-5322] = 1.000000e-019; 00328 m[5332] = 1.550000e-012; 00329 m[-5332] = 1.550000e-012; 00330 m[5132] = 1.390000e-012; 00331 m[-5132] = 1.390000e-012; 00332 m[5232] = 1.390000e-012; 00333 m[-5232] = 1.390000e-012; 00334 m[100443] = 2.194041e-021; 00335 m[331] = 3.258476e-021; 00336 m[441] = 4.113826e-023; 00337 m[10441] = 4.063038e-023; 00338 m[20443] = 7.154480e-022; 00339 m[445] = 3.164482e-022; 00340 m[9000111] = 1.149997e-023; 00341 m[9000211] = 1.149997e-023; 00342 m[-9000211] = 1.149997e-023; 00343 m[20113] = 1.500000e-024; 00344 m[115] = 6.151516e-024; 00345 m[215] = 6.151516e-024; 00346 m[-215] = 6.151516e-024; 00347 m[10323] = 7.313469e-024; 00348 m[-10323] = 7.313469e-024; 00349 m[10313] = 7.313469e-024; 00350 m[-10313] = 7.313469e-024; 00351 m[20323] = 3.782829e-024; 00352 m[-20323] = 3.782829e-024; 00353 m[20313] = 3.782829e-024; 00354 m[-20313] = 3.782829e-024; 00355 m[10321] = 2.238817e-024; 00356 m[-10321] = 2.238817e-024; 00357 m[10311] = 2.238817e-024; 00358 m[-10311] = 2.238817e-024; 00359 m[325] = 6.682357e-024; 00360 m[-325] = 6.682357e-024; 00361 m[315] = 6.038644e-024; 00362 m[-315] = 6.038644e-024; 00363 m[10411] = 4.380000e-024; 00364 m[20413] = 2.630000e-024; 00365 m[10413] = 3.290000e-023; 00366 m[-415] = 2.632849e-023; 00367 m[-10411] = 4.380000e-024; 00368 m[-20413] = 2.630000e-024; 00369 m[-10413] = 3.290000e-023; 00370 m[415] = 2.632849e-023; 00371 m[10421] = 4.380000e-024; 00372 m[20423] = 2.630000e-024; 00373 m[10423] = 3.482604e-023; 00374 m[-425] = 2.861792e-023; 00375 m[-10421] = 4.380000e-024; 00376 m[-20423] = 2.630000e-024; 00377 m[-10423] = 3.482604e-023; 00378 m[425] = 2.861792e-023; 00379 m[10431] = 6.582100e-022; 00380 m[20433] = 6.582100e-022; 00381 m[10433] = 6.582100e-022; 00382 m[435] = 4.388100e-023; 00383 m[-10431] = 6.582100e-022; 00384 m[-20433] = 6.582100e-022; 00385 m[-10433] = 6.582100e-022; 00386 m[-435] = 4.388100e-023; 00387 m[2224] = 5.485102e-024; 00388 m[2214] = 5.485102e-024; 00389 m[2114] = 5.485102e-024; 00390 m[1114] = 5.485102e-024; 00391 m[-2224] = 5.485102e-024; 00392 m[-2214] = 5.485102e-024; 00393 m[-2114] = 5.485102e-024; 00394 m[-1114] = 5.485102e-024; 00395 m[-523] = 1.000019e-019; 00396 m[523] = 1.000019e-019; 00397 m[513] = 1.000019e-019; 00398 m[-513] = 1.000019e-019; 00399 m[533] = 1.000000e-019; 00400 m[-533] = 1.000000e-019; 00401 m[10521] = 4.390000e-024; 00402 m[20523] = 2.630000e-024; 00403 m[10523] = 1.650000e-023; 00404 m[525] = 1.310000e-023; 00405 m[-10521] = 4.390000e-024; 00406 m[-20523] = 2.630000e-024; 00407 m[-10523] = 1.650000e-023; 00408 m[-525] = 1.310000e-023; 00409 m[10511] = 4.390000e-024; 00410 m[20513] = 2.630000e-024; 00411 m[10513] = 1.650000e-023; 00412 m[515] = 1.310000e-023; 00413 m[-10511] = 4.390000e-024; 00414 m[-20513] = 2.630000e-024; 00415 m[-10513] = 1.650000e-023; 00416 m[-515] = 1.310000e-023; 00417 m[10531] = 4.390000e-024; 00418 m[20533] = 2.630000e-024; 00419 m[10533] = 1.650000e-023; 00420 m[535] = 1.310000e-023; 00421 m[-10531] = 4.390000e-024; 00422 m[-20533] = 2.630000e-024; 00423 m[-10533] = 1.650000e-023; 00424 m[-535] = 1.310000e-023; 00425 m[14] = 1.000000e+016; 00426 m[-14] = 1.000000e+016; 00427 m[-12] = 1.000000e+016; 00428 m[1] = 0.000000e+000; 00429 m[-1] = 0.000000e+000; 00430 m[2] = 0.000000e+000; 00431 m[-2] = 0.000000e+000; 00432 m[3] = 0.000000e+000; 00433 m[-3] = 0.000000e+000; 00434 m[4] = 0.000000e+000; 00435 m[-4] = 0.000000e+000; 00436 m[5] = 0.000000e+000; 00437 m[-5] = 0.000000e+000; 00438 m[6] = 4.707703e-025; 00439 m[-6] = 4.707703e-025; 00440 m[7] = 0.000000e+000; 00441 m[-7] = 0.000000e+000; 00442 m[8] = 0.000000e+000; 00443 m[-8] = 0.000000e+000; 00444 m[16] = 1.000000e+016; 00445 m[-16] = 1.000000e+016; 00446 m[17] = 0.000000e+000; 00447 m[-17] = 0.000000e+000; 00448 m[18] = 0.000000e+000; 00449 m[-18] = 0.000000e+000; 00450 m[21] = 0.000000e+000; 00451 m[32] = 0.000000e+000; 00452 m[33] = 0.000000e+000; 00453 m[34] = 0.000000e+000; 00454 m[-34] = 0.000000e+000; 00455 m[39] = 0.000000e+000; 00456 m[41] = 0.000000e+000; 00457 m[-41] = 0.000000e+000; 00458 m[42] = 0.000000e+000; 00459 m[-42] = 0.000000e+000; 00460 m[43] = 0.000000e+000; 00461 m[44] = 0.000000e+000; 00462 m[-44] = 0.000000e+000; 00463 m[81] = 0.000000e+000; 00464 m[82] = 0.000000e+000; 00465 m[-82] = 0.000000e+000; 00466 m[83] = 0.000000e+000; 00467 m[84] = 3.335641e-013; 00468 m[-84] = 3.335641e-013; 00469 m[85] = 1.290893e-012; 00470 m[-85] = 1.290893e-012; 00471 m[86] = 0.000000e+000; 00472 m[-86] = 0.000000e+000; 00473 m[87] = 0.000000e+000; 00474 m[-87] = 0.000000e+000; 00475 m[88] = 0.000000e+000; 00476 m[90] = 0.000000e+000; 00477 m[91] = 0.000000e+000; 00478 m[92] = 0.000000e+000; 00479 m[93] = 0.000000e+000; 00480 m[94] = 0.000000e+000; 00481 m[95] = 0.000000e+000; 00482 m[96] = 0.000000e+000; 00483 m[97] = 0.000000e+000; 00484 m[98] = 0.000000e+000; 00485 m[99] = 0.000000e+000; 00486 m[117] = 4.088275e-024; 00487 m[119] = 1.828367e-024; 00488 m[217] = 4.088275e-024; 00489 m[-217] = 4.088275e-024; 00490 m[219] = 1.828367e-024; 00491 m[-219] = 1.828367e-024; 00492 m[225] = 3.555982e-024; 00493 m[227] = 3.917930e-024; 00494 m[229] = 3.392846e-024; 00495 m[311] = 1.000000e+016; 00496 m[-311] = 1.000000e+016; 00497 m[317] = 4.139699e-024; 00498 m[-317] = 4.139699e-024; 00499 m[319] = 3.324304e-024; 00500 m[-319] = 3.324304e-024; 00501 m[327] = 4.139699e-024; 00502 m[-327] = 4.139699e-024; 00503 m[329] = 3.324304e-024; 00504 m[-329] = 3.324304e-024; 00505 m[335] = 8.660687e-024; 00506 m[337] = 7.565657e-024; 00507 m[543] = 0.000000e+000; 00508 m[-543] = 0.000000e+000; 00509 m[545] = 0.000000e+000; 00510 m[-545] = 0.000000e+000; 00511 m[551] = 0.000000e+000; 00512 m[553] = 1.253738e-020; 00513 m[555] = 1.000000e+016; 00514 m[557] = 0.000000e+000; 00515 m[-450000000] = 0.000000e+000; 00516 m[-490000000] = 0.000000e+000; 00517 m[-460000000] = 0.000000e+000; 00518 m[-470000000] = 0.000000e+000; 00519 m[1103] = 0.000000e+000; 00520 m[-1103] = 0.000000e+000; 00521 m[1112] = 4.388081e-024; 00522 m[-1112] = 4.388081e-024; 00523 m[1116] = 1.880606e-024; 00524 m[-1116] = 1.880606e-024; 00525 m[1118] = 2.194041e-024; 00526 m[-1118] = 2.194041e-024; 00527 m[1212] = 4.388081e-024; 00528 m[-1212] = 4.388081e-024; 00529 m[1214] = 5.485102e-024; 00530 m[-1214] = 5.485102e-024; 00531 m[1216] = 1.880606e-024; 00532 m[-1216] = 1.880606e-024; 00533 m[1218] = 1.462694e-024; 00534 m[-1218] = 1.462694e-024; 00535 m[2101] = 0.000000e+000; 00536 m[-2101] = 0.000000e+000; 00537 m[2103] = 0.000000e+000; 00538 m[-2103] = 0.000000e+000; 00539 m[2116] = 4.388081e-024; 00540 m[-2116] = 4.388081e-024; 00541 m[2118] = 2.194041e-024; 00542 m[-2118] = 2.194041e-024; 00543 m[2122] = 4.388081e-024; 00544 m[-2122] = 4.388081e-024; 00545 m[2124] = 5.485102e-024; 00546 m[-2124] = 5.485102e-024; 00547 m[2126] = 1.880606e-024; 00548 m[-2126] = 1.880606e-024; 00549 m[2128] = 1.462694e-024; 00550 m[-2128] = 1.462694e-024; 00551 m[2203] = 0.000000e+000; 00552 m[-2203] = 0.000000e+000; 00553 m[2216] = 4.388081e-024; 00554 m[-2216] = 4.388081e-024; 00555 m[2218] = 2.194041e-024; 00556 m[-2218] = 2.194041e-024; 00557 m[2222] = 4.388081e-024; 00558 m[-2222] = 4.388081e-024; 00559 m[2226] = 1.880606e-024; 00560 m[-2226] = 1.880606e-024; 00561 m[2228] = 2.194041e-024; 00562 m[-2228] = 2.194041e-024; 00563 m[3101] = 0.000000e+000; 00564 m[-3101] = 0.000000e+000; 00565 m[3103] = 0.000000e+000; 00566 m[-3103] = 0.000000e+000; 00567 m[3114] = 1.670589e-023; 00568 m[-3114] = 1.670589e-023; 00569 m[3116] = 5.485102e-024; 00570 m[-3116] = 5.485102e-024; 00571 m[3118] = 3.656734e-024; 00572 m[-3118] = 3.656734e-024; 00573 m[3124] = 4.219309e-023; 00574 m[-3124] = 4.219309e-023; 00575 m[3126] = 8.227653e-024; 00576 m[-3126] = 8.227653e-024; 00577 m[3128] = 3.291061e-024; 00578 m[-3128] = 3.291061e-024; 00579 m[3201] = 0.000000e+000; 00580 m[-3201] = 0.000000e+000; 00581 m[3203] = 0.000000e+000; 00582 m[-3203] = 0.000000e+000; 00583 m[3214] = 1.828367e-023; 00584 m[-3214] = 1.828367e-023; 00585 m[3216] = 5.485102e-024; 00586 m[-3216] = 5.485102e-024; 00587 m[3218] = 3.656734e-024; 00588 m[-3218] = 3.656734e-024; 00589 m[3224] = 1.838582e-023; 00590 m[-3224] = 1.838582e-023; 00591 m[3226] = 5.485102e-024; 00592 m[-3226] = 5.485102e-024; 00593 m[3228] = 3.656734e-024; 00594 m[-3228] = 3.656734e-024; 00595 m[3303] = 0.000000e+000; 00596 m[-3303] = 0.000000e+000; 00597 m[3314] = 6.648608e-023; 00598 m[-3314] = 6.648608e-023; 00599 m[3324] = 7.233101e-023; 00600 m[-3324] = 7.233101e-023; 00601 m[4101] = 0.000000e+000; 00602 m[-4101] = 0.000000e+000; 00603 m[4103] = 0.000000e+000; 00604 m[-4103] = 0.000000e+000; 00605 m[4114] = 0.000000e+000; 00606 m[-4114] = 0.000000e+000; 00607 m[4201] = 0.000000e+000; 00608 m[-4201] = 0.000000e+000; 00609 m[4203] = 0.000000e+000; 00610 m[-4203] = 0.000000e+000; 00611 m[4214] = 3.291061e-022; 00612 m[-4214] = 3.291061e-022; 00613 m[4224] = 0.000000e+000; 00614 m[-4224] = 0.000000e+000; 00615 m[4301] = 0.000000e+000; 00616 m[-4301] = 0.000000e+000; 00617 m[4303] = 0.000000e+000; 00618 m[-4303] = 0.000000e+000; 00619 m[4314] = 0.000000e+000; 00620 m[-4314] = 0.000000e+000; 00621 m[4324] = 0.000000e+000; 00622 m[-4324] = 0.000000e+000; 00623 m[4334] = 0.000000e+000; 00624 m[-4334] = 0.000000e+000; 00625 m[4403] = 0.000000e+000; 00626 m[-4403] = 0.000000e+000; 00627 m[4412] = 3.335641e-013; 00628 m[-4412] = 3.335641e-013; 00629 m[4414] = 3.335641e-013; 00630 m[-4414] = 3.335641e-013; 00631 m[4422] = 3.335641e-013; 00632 m[-4422] = 3.335641e-013; 00633 m[4424] = 3.335641e-013; 00634 m[-4424] = 3.335641e-013; 00635 m[4432] = 3.335641e-013; 00636 m[-4432] = 3.335641e-013; 00637 m[4434] = 3.335641e-013; 00638 m[-4434] = 3.335641e-013; 00639 m[4444] = 3.335641e-013; 00640 m[-4444] = 3.335641e-013; 00641 m[5101] = 0.000000e+000; 00642 m[-5101] = 0.000000e+000; 00643 m[5103] = 0.000000e+000; 00644 m[-5103] = 0.000000e+000; 00645 m[5114] = 0.000000e+000; 00646 m[-5114] = 0.000000e+000; 00647 m[5142] = 1.290893e-012; 00648 m[-5142] = 1.290893e-012; 00649 m[5201] = 0.000000e+000; 00650 m[-5201] = 0.000000e+000; 00651 m[5203] = 0.000000e+000; 00652 m[-5203] = 0.000000e+000; 00653 m[5214] = 0.000000e+000; 00654 m[-5214] = 0.000000e+000; 00655 m[5224] = 0.000000e+000; 00656 m[-5224] = 0.000000e+000; 00657 m[5242] = 1.290893e-012; 00658 m[-5242] = 1.290893e-012; 00659 m[5301] = 0.000000e+000; 00660 m[-5301] = 0.000000e+000; 00661 m[5303] = 0.000000e+000; 00662 m[-5303] = 0.000000e+000; 00663 m[5314] = 0.000000e+000; 00664 m[-5314] = 0.000000e+000; 00665 m[5324] = 0.000000e+000; 00666 m[-5324] = 0.000000e+000; 00667 m[5334] = 0.000000e+000; 00668 m[-5334] = 0.000000e+000; 00669 m[5342] = 1.290893e-012; 00670 m[-5342] = 1.290893e-012; 00671 m[5401] = 0.000000e+000; 00672 m[-5401] = 0.000000e+000; 00673 m[5403] = 0.000000e+000; 00674 m[-5403] = 0.000000e+000; 00675 m[5412] = 1.290893e-012; 00676 m[-5412] = 1.290893e-012; 00677 m[5414] = 1.290893e-012; 00678 m[-5414] = 1.290893e-012; 00679 m[5422] = 1.290893e-012; 00680 m[-5422] = 1.290893e-012; 00681 m[5424] = 1.290893e-012; 00682 m[-5424] = 1.290893e-012; 00683 m[5432] = 1.290893e-012; 00684 m[-5432] = 1.290893e-012; 00685 m[5434] = 1.290893e-012; 00686 m[-5434] = 1.290893e-012; 00687 m[5442] = 1.290893e-012; 00688 m[-5442] = 1.290893e-012; 00689 m[5444] = 1.290893e-012; 00690 m[-5444] = 1.290893e-012; 00691 m[5503] = 0.000000e+000; 00692 m[-5503] = 0.000000e+000; 00693 m[5512] = 1.290893e-012; 00694 m[-5512] = 1.290893e-012; 00695 m[5514] = 1.290893e-012; 00696 m[-5514] = 1.290893e-012; 00697 m[5522] = 1.290893e-012; 00698 m[-5522] = 1.290893e-012; 00699 m[5524] = 1.290893e-012; 00700 m[-5524] = 1.290893e-012; 00701 m[5532] = 1.290893e-012; 00702 m[-5532] = 1.290893e-012; 00703 m[5534] = 1.290893e-012; 00704 m[-5534] = 1.290893e-012; 00705 m[5542] = 1.290893e-012; 00706 m[-5542] = 1.290893e-012; 00707 m[5544] = 1.290893e-012; 00708 m[-5544] = 1.290893e-012; 00709 m[5554] = 1.290893e-012; 00710 m[-5554] = 1.290893e-012; 00711 m[10022] = 0.000000e+000; 00712 m[10111] = 2.483820e-024; 00713 m[10113] = 4.635297e-024; 00714 m[10115] = 2.541360e-024; 00715 m[10211] = 2.483820e-024; 00716 m[-10211] = 2.483820e-024; 00717 m[10213] = 4.635297e-024; 00718 m[-10213] = 4.635297e-024; 00719 m[10215] = 2.541360e-024; 00720 m[-10215] = 2.541360e-024; 00721 m[9010221] = 1.316424e-023; 00722 m[10223] = 1.828367e-024; 00723 m[10225] = 0.000000e+000; 00724 m[10315] = 3.538775e-024; 00725 m[-10315] = 3.538775e-024; 00726 m[10325] = 3.538775e-024; 00727 m[-10325] = 3.538775e-024; 00728 m[10331] = 5.265698e-024; 00729 m[10333] = 0.000000e+000; 00730 m[10335] = 0.000000e+000; 00731 m[10443] = 0.000000e+000; 00732 m[10541] = 0.000000e+000; 00733 m[-10541] = 0.000000e+000; 00734 m[10543] = 0.000000e+000; 00735 m[-10543] = 0.000000e+000; 00736 m[10551] = 1.000000e+016; 00737 m[10553] = 0.000000e+000; 00738 m[10555] = 0.000000e+000; 00739 m[11112] = 0.000000e+000; 00740 m[-11112] = 0.000000e+000; 00741 m[11114] = 2.194041e-024; 00742 m[-11114] = 2.194041e-024; 00743 m[11116] = 1.880606e-024; 00744 m[-11116] = 1.880606e-024; 00745 m[11212] = 1.880606e-024; 00746 m[-11212] = 1.880606e-024; 00747 m[11216] = 0.000000e+000; 00748 m[-11216] = 0.000000e+000; 00749 m[12112] = 1.880606e-024; 00750 m[-12112] = 1.880606e-024; 00751 m[12114] = 2.194041e-024; 00752 m[-12114] = 2.194041e-024; 00753 m[12116] = 5.063171e-024; 00754 m[-12116] = 5.063171e-024; 00755 m[12118] = 0.000000e+000; 00756 m[-12118] = 0.000000e+000; 00757 m[12122] = 0.000000e+000; 00758 m[-12122] = 0.000000e+000; 00759 m[12126] = 1.880606e-024; 00760 m[-12126] = 1.880606e-024; 00761 m[12212] = 1.880606e-024; 00762 m[-12212] = 1.880606e-024; 00763 m[12214] = 2.194041e-024; 00764 m[-12214] = 2.194041e-024; 00765 m[12216] = 5.063171e-024; 00766 m[-12216] = 5.063171e-024; 00767 m[12218] = 0.000000e+000; 00768 m[-12218] = 0.000000e+000; 00769 m[12222] = 0.000000e+000; 00770 m[-12222] = 0.000000e+000; 00771 m[12224] = 2.194041e-024; 00772 m[-12224] = 2.194041e-024; 00773 m[12226] = 1.880606e-024; 00774 m[-12226] = 1.880606e-024; 00775 m[13112] = 6.582122e-024; 00776 m[-13112] = 6.582122e-024; 00777 m[13114] = 1.097020e-023; 00778 m[-13114] = 1.097020e-023; 00779 m[13116] = 5.485102e-024; 00780 m[-13116] = 5.485102e-024; 00781 m[13122] = 1.316424e-023; 00782 m[-13122] = 1.316424e-023; 00783 m[13124] = 1.097020e-023; 00784 m[-13124] = 1.097020e-023; 00785 m[13126] = 6.928549e-024; 00786 m[-13126] = 6.928549e-024; 00787 m[13212] = 6.582122e-024; 00788 m[-13212] = 6.582122e-024; 00789 m[13214] = 1.097020e-023; 00790 m[-13214] = 1.097020e-023; 00791 m[13216] = 5.485102e-024; 00792 m[-13216] = 5.485102e-024; 00793 m[13222] = 6.582122e-024; 00794 m[-13222] = 6.582122e-024; 00795 m[13224] = 1.097020e-023; 00796 m[-13224] = 1.097020e-023; 00797 m[13226] = 5.485102e-024; 00798 m[-13226] = 5.485102e-024; 00799 m[13314] = 2.742551e-023; 00800 m[-13314] = 2.742551e-023; 00801 m[13316] = 0.000000e+000; 00802 m[-13316] = 0.000000e+000; 00803 m[13324] = 2.742551e-023; 00804 m[-13324] = 2.742551e-023; 00805 m[13326] = 0.000000e+000; 00806 m[-13326] = 0.000000e+000; 00807 m[14122] = 1.828367e-022; 00808 m[-14122] = 1.828367e-022; 00809 m[14124] = 0.000000e+000; 00810 m[-14124] = 0.000000e+000; 00811 m[10221] = 2.194040e-024; 00812 m[20223] = 2.742551e-023; 00813 m[20315] = 2.384827e-024; 00814 m[-20315] = 2.384827e-024; 00815 m[20325] = 2.384827e-024; 00816 m[-20325] = 2.384827e-024; 00817 m[20333] = 1.185968e-023; 00818 m[20543] = 0.000000e+000; 00819 m[-20543] = 0.000000e+000; 00820 m[20553] = 1.000000e+016; 00821 m[20555] = 0.000000e+000; 00822 m[21112] = 2.632849e-024; 00823 m[-21112] = 2.632849e-024; 00824 m[21114] = 3.291061e-024; 00825 m[-21114] = 3.291061e-024; 00826 m[21212] = 2.632849e-024; 00827 m[-21212] = 2.632849e-024; 00828 m[21214] = 6.582122e-024; 00829 m[-21214] = 6.582122e-024; 00830 m[22112] = 4.388081e-024; 00831 m[-22112] = 4.388081e-024; 00832 m[22114] = 3.291061e-024; 00833 m[-22114] = 3.291061e-024; 00834 m[22122] = 2.632849e-024; 00835 m[-22122] = 2.632849e-024; 00836 m[22124] = 6.582122e-024; 00837 m[-22124] = 6.582122e-024; 00838 m[22212] = 4.388081e-024; 00839 m[-22212] = 4.388081e-024; 00840 m[22214] = 3.291061e-024; 00841 m[-22214] = 3.291061e-024; 00842 m[22222] = 2.632849e-024; 00843 m[-22222] = 2.632849e-024; 00844 m[22224] = 3.291061e-024; 00845 m[-22224] = 3.291061e-024; 00846 m[23112] = 7.313469e-024; 00847 m[-23112] = 7.313469e-024; 00848 m[23114] = 2.991874e-024; 00849 m[-23114] = 2.991874e-024; 00850 m[23122] = 4.388081e-024; 00851 m[-23122] = 4.388081e-024; 00852 m[23124] = 6.582122e-024; 00853 m[-23124] = 6.582122e-024; 00854 m[23126] = 3.291061e-024; 00855 m[-23126] = 3.291061e-024; 00856 m[23212] = 7.313469e-024; 00857 m[-23212] = 7.313469e-024; 00858 m[23214] = 2.991874e-024; 00859 m[-23214] = 2.991874e-024; 00860 m[23222] = 7.313469e-024; 00861 m[-23222] = 7.313469e-024; 00862 m[23224] = 2.991874e-024; 00863 m[-23224] = 2.991874e-024; 00864 m[23314] = 0.000000e+000; 00865 m[-23314] = 0.000000e+000; 00866 m[23324] = 0.000000e+000; 00867 m[-23324] = 0.000000e+000; 00868 m[30113] = 2.742551e-024; 00869 m[30213] = 2.742551e-024; 00870 m[-30213] = 2.742551e-024; 00871 m[30223] = 2.991874e-024; 00872 m[30313] = 2.056913e-024; 00873 m[-30313] = 2.056913e-024; 00874 m[30323] = 2.056913e-024; 00875 m[-30323] = 2.056913e-024; 00876 m[30343] = 0.000000e+000; 00877 m[-30343] = 0.000000e+000; 00878 m[30353] = 0.000000e+000; 00879 m[-30353] = 0.000000e+000; 00880 m[30363] = 0.000000e+000; 00881 m[-30363] = 0.000000e+000; 00882 m[30411] = 0.000000e+000; 00883 m[-30411] = 0.000000e+000; 00884 m[30413] = 0.000000e+000; 00885 m[-30413] = 0.000000e+000; 00886 m[30421] = 0.000000e+000; 00887 m[-30421] = 0.000000e+000; 00888 m[30423] = 0.000000e+000; 00889 m[-30423] = 0.000000e+000; 00890 m[30443] = 2.789035e-023; 00891 m[30553] = 0.000000e+000; 00892 m[31114] = 1.880606e-024; 00893 m[-31114] = 1.880606e-024; 00894 m[31214] = 4.388081e-024; 00895 m[-31214] = 4.388081e-024; 00896 m[32112] = 4.388081e-024; 00897 m[-32112] = 4.388081e-024; 00898 m[32114] = 1.880606e-024; 00899 m[-32114] = 1.880606e-024; 00900 m[32124] = 4.388081e-024; 00901 m[-32124] = 4.388081e-024; 00902 m[32212] = 4.388081e-024; 00903 m[-32212] = 4.388081e-024; 00904 m[32214] = 1.880606e-024; 00905 m[-32214] = 1.880606e-024; 00906 m[32224] = 1.880606e-024; 00907 m[-32224] = 1.880606e-024; 00908 m[33122] = 1.880606e-023; 00909 m[-33122] = 1.880606e-023; 00910 m[33314] = 0.000000e+000; 00911 m[-33314] = 0.000000e+000; 00912 m[33324] = 0.000000e+000; 00913 m[-33324] = 0.000000e+000; 00914 m[41214] = 0.000000e+000; 00915 m[-41214] = 0.000000e+000; 00916 m[42112] = 6.582122e-024; 00917 m[-42112] = 6.582122e-024; 00918 m[42124] = 0.000000e+000; 00919 m[-42124] = 0.000000e+000; 00920 m[42212] = 6.582122e-024; 00921 m[-42212] = 6.582122e-024; 00922 m[43122] = 2.194041e-024; 00923 m[-43122] = 2.194041e-024; 00924 m[52114] = 0.000000e+000; 00925 m[-52114] = 0.000000e+000; 00926 m[52214] = 0.000000e+000; 00927 m[-52214] = 0.000000e+000; 00928 m[53122] = 4.388081e-024; 00929 m[-53122] = 4.388081e-024; 00930 m[100111] = 1.645531e-024; 00931 m[100113] = 2.123265e-024; 00932 m[100211] = 1.645531e-024; 00933 m[-100211] = 1.645531e-024; 00934 m[100213] = 2.123265e-024; 00935 m[-100213] = 2.123265e-024; 00936 m[100221] = 1.196749e-023; 00937 m[100223] = 3.871836e-024; 00938 m[100225] = 0.000000e+000; 00939 m[100311] = 0.000000e+000; 00940 m[-100311] = 0.000000e+000; 00941 m[100313] = 2.837122e-024; 00942 m[-100313] = 2.837122e-024; 00943 m[100315] = 0.000000e+000; 00944 m[-100315] = 0.000000e+000; 00945 m[100321] = 0.000000e+000; 00946 m[-100321] = 0.000000e+000; 00947 m[100323] = 2.837122e-024; 00948 m[-100323] = 2.837122e-024; 00949 m[100325] = 0.000000e+000; 00950 m[-100325] = 0.000000e+000; 00951 m[100331] = 0.000000e+000; 00952 m[100333] = 4.388081e-024; 00953 m[100335] = 3.291061e-024; 00954 m[100441] = 0.000000e+000; 00955 m[100551] = 0.000000e+000; 00956 m[100553] = 1.495937e-020; 00957 m[100555] = 1.000000e+016; 00958 m[100557] = 0.000000e+000; 00959 m[110551] = 1.000000e+016; 00960 m[110553] = 0.000000e+000; 00961 m[110555] = 0.000000e+000; 00962 m[120553] = 1.000000e+016; 00963 m[120555] = 0.000000e+000; 00964 m[130553] = 0.000000e+000; 00965 m[200111] = 3.134344e-024; 00966 m[200211] = 3.134344e-024; 00967 m[-200211] = 3.134344e-024; 00968 m[200551] = 0.000000e+000; 00969 m[200553] = 2.502708e-020; 00970 m[200555] = 0.000000e+000; 00971 m[210551] = 0.000000e+000; 00972 m[210553] = 0.000000e+000; 00973 m[220553] = 0.000000e+000; 00974 m[300553] = 4.701516e-023; 00975 m[9000221] = 0.000000e+000; 00976 m[9000443] = 1.265793e-023; 00977 m[9000553] = 5.983747e-024; 00978 m[9010443] = 8.438618e-024; 00979 m[9010553] = 8.331800e-024; 00980 m[9020221] = 6.038644e-024; 00981 m[9020443] = 1.530726e-023; 00982 m[9060225] = 4.388081e-024; 00983 m[9070225] = 2.056913e-024; 00984 m[1000001] = 0.000000e+000; 00985 m[-1000001] = 0.000000e+000; 00986 m[1000002] = 0.000000e+000; 00987 m[-1000002] = 0.000000e+000; 00988 m[1000003] = 0.000000e+000; 00989 m[-1000003] = 0.000000e+000; 00990 m[1000004] = 0.000000e+000; 00991 m[-1000004] = 0.000000e+000; 00992 m[1000005] = 0.000000e+000; 00993 m[-1000005] = 0.000000e+000; 00994 m[1000006] = 0.000000e+000; 00995 m[-1000006] = 0.000000e+000; 00996 m[1000011] = 0.000000e+000; 00997 m[-1000011] = 0.000000e+000; 00998 m[1000012] = 0.000000e+000; 00999 m[-1000012] = 0.000000e+000; 01000 m[1000013] = 0.000000e+000; 01001 m[-1000013] = 0.000000e+000; 01002 m[1000014] = 0.000000e+000; 01003 m[-1000014] = 0.000000e+000; 01004 m[1000015] = 0.000000e+000; 01005 m[-1000015] = 0.000000e+000; 01006 m[1000016] = 0.000000e+000; 01007 m[-1000016] = 0.000000e+000; 01008 m[1000021] = 0.000000e+000; 01009 m[1000022] = 0.000000e+000; 01010 m[1000023] = 0.000000e+000; 01011 m[1000024] = 0.000000e+000; 01012 m[-1000024] = 0.000000e+000; 01013 m[1000025] = 0.000000e+000; 01014 m[1000035] = 0.000000e+000; 01015 m[1000037] = 0.000000e+000; 01016 m[-1000037] = 0.000000e+000; 01017 m[1000039] = 0.000000e+000; 01018 m[2000001] = 0.000000e+000; 01019 m[-2000001] = 0.000000e+000; 01020 m[2000002] = 0.000000e+000; 01021 m[-2000002] = 0.000000e+000; 01022 m[2000003] = 0.000000e+000; 01023 m[-2000003] = 0.000000e+000; 01024 m[2000004] = 0.000000e+000; 01025 m[-2000004] = 0.000000e+000; 01026 m[2000005] = 0.000000e+000; 01027 m[-2000005] = 0.000000e+000; 01028 m[2000006] = 0.000000e+000; 01029 m[-2000006] = 0.000000e+000; 01030 m[2000011] = 0.000000e+000; 01031 m[-2000011] = 0.000000e+000; 01032 m[2000012] = 0.000000e+000; 01033 m[-2000012] = 0.000000e+000; 01034 m[2000013] = 0.000000e+000; 01035 m[-2000013] = 0.000000e+000; 01036 m[2000014] = 0.000000e+000; 01037 m[-2000014] = 0.000000e+000; 01038 m[2000015] = 0.000000e+000; 01039 m[-2000015] = 0.000000e+000; 01040 m[2000016] = 0.000000e+000; 01041 m[-2000016] = 0.000000e+000; 01042 m[3000111] = 0.000000e+000; 01043 m[3000113] = 0.000000e+000; 01044 m[3000211] = 0.000000e+000; 01045 m[-3000211] = 0.000000e+000; 01046 m[3000213] = 0.000000e+000; 01047 m[-3000213] = 0.000000e+000; 01048 m[3000221] = 0.000000e+000; 01049 m[3000223] = 0.000000e+000; 01050 m[3000331] = 0.000000e+000; 01051 m[3100021] = 0.000000e+000; 01052 m[3100111] = 0.000000e+000; 01053 m[3100113] = 0.000000e+000; 01054 m[3200111] = 0.000000e+000; 01055 m[3200113] = 0.000000e+000; 01056 m[3300113] = 0.000000e+000; 01057 m[3400113] = 0.000000e+000; 01058 m[4000001] = 0.000000e+000; 01059 m[-4000001] = 0.000000e+000; 01060 m[4000002] = 0.000000e+000; 01061 m[-4000002] = 0.000000e+000; 01062 m[4000011] = 0.000000e+000; 01063 m[-4000011] = 0.000000e+000; 01064 m[4000012] = 0.000000e+000; 01065 m[-4000012] = 0.000000e+000; 01066 m[5000039] = 0.000000e+000; 01067 m[9900012] = 0.000000e+000; 01068 m[9900014] = 0.000000e+000; 01069 m[9900016] = 0.000000e+000; 01070 m[9900023] = 0.000000e+000; 01071 m[9900024] = 0.000000e+000; 01072 m[-9900024] = 0.000000e+000; 01073 m[9900041] = 0.000000e+000; 01074 m[-9900041] = 0.000000e+000; 01075 m[9900042] = 0.000000e+000; 01076 m[-9900042] = 0.000000e+000; 01077 m[1027013000] = 0.000000e+000; 01078 m[1012006000] = 0.000000e+000; 01079 m[1063029000] = 0.000000e+000; 01080 m[1014007000] = 0.000000e+000; 01081 m[1016008000] = 0.000000e+000; 01082 m[1028014000] = 0.000000e+000; 01083 m[1065029000] = 0.000000e+000; 01084 m[1009004000] = 0.000000e+000; 01085 m[1019009000] = 0.000000e+000; 01086 m[1056026000] = 0.000000e+000; 01087 m[1207082000] = 0.000000e+000; 01088 m[1208082000] = 0.000000e+000; 01089 m[1029014000] = 0.000000e+000; 01090 m[1206082000] = 0.000000e+000; 01091 m[1054026000] = 0.000000e+000; 01092 m[1018008000] = 0.000000e+000; 01093 m[1030014000] = 0.000000e+000; 01094 m[1057026000] = 0.000000e+000; 01095 m[1204082000] = 0.000000e+000; 01096 m[-99000000] = 0.000000e+000; 01097 m[1028013000] = 0.000000e+000; 01098 m[1040018000] = 0.000000e+000; 01099 m[1011005000] = 0.000000e+000; 01100 m[1012005000] = 0.000000e+000; 01101 m[1013006000] = 0.000000e+000; 01102 m[1014006000] = 0.000000e+000; 01103 m[1052024000] = 0.000000e+000; 01104 m[1024012000] = 0.000000e+000; 01105 m[1026012000] = 0.000000e+000; 01106 m[1027012000] = 0.000000e+000; 01107 m[1015007000] = 0.000000e+000; 01108 m[1022010000] = 0.000000e+000; 01109 m[1058028000] = 0.000000e+000; 01110 m[1060028000] = 0.000000e+000; 01111 m[1062028000] = 0.000000e+000; 01112 m[1064028000] = 0.000000e+000; 01113 m[1007003000] = 0.000000e+000; 01114 m[1025012000] = 0.000000e+000; 01115 m[1053024000] = 0.000000e+000; 01116 m[1055025000] = 0.000000e+000; 01117 m[1008004000] = 0.000000e+000; 01118 m[1010004000] = 0.000000e+000; 01119 m[1010005000] = 0.000000e+000; 01120 m[1016007000] = 0.000000e+000; 01121 m[1017008000] = 0.000000e+000; 01122 m[1019008000] = 0.000000e+000; 01123 m[1023010000] = 0.000000e+000; 01124 m[1024011000] = 0.000000e+000; 01125 m[1031015000] = 0.000000e+000; 01126 m[1039017000] = 0.000000e+000; 01127 m[1040017000] = 0.000000e+000; 01128 m[1036018000] = 0.000000e+000; 01129 m[1050024000] = 0.000000e+000; 01130 m[1054024000] = 0.000000e+000; 01131 m[1059026000] = 0.000000e+000; 01132 m[1061028000] = 0.000000e+000; 01133 m[1063028000] = 0.000000e+000; 01134 m[1092042000] = 0.000000e+000; 01135 m[1095042000] = 0.000000e+000; 01136 m[1096042000] = 0.000000e+000; 01137 m[1097042000] = 0.000000e+000; 01138 m[1098042000] = 0.000000e+000; 01139 m[1100042000] = 0.000000e+000; 01140 m[1108046000] = 0.000000e+000; 01141 01142 // Added by hand: 01143 m[9902210] = 0.000000e+000; //diffractive p-state -> assume no lifetime 01144 return true; 01145 } 01146 01147 01148 private: 01149 01150 /// @name Histograms 01151 //@{ 01152 Histo1DPtr _h_mult_total; // full kinematic range 01153 Histo1DPtr _h_mult_eta[5]; // in eta bins 01154 Histo1DPtr _h_mult_pt[5]; // in pT bins 01155 Histo1DPtr _h_dndeta; // density dn/deta 01156 Histo1DPtr _h_dndpt; // density dn/dpT 01157 //@} 01158 01159 01160 /// @name Private variables 01161 double _p_min; 01162 double _pt_min; 01163 double _eta_min; 01164 double _eta_max; 01165 double _maxlft; 01166 01167 /// Count selected events 01168 double _sumW; 01169 01170 map<int, double> _partLftMap; // Map <PDGID, PDGLIFETIME> 01171 01172 }; 01173 01174 01175 // The hook for the plugin system 01176 DECLARE_RIVET_PLUGIN(LHCB_2014_I1281685); 01177 01178 } Generated on Thu Mar 10 2016 08:29:51 for The Rivet MC analysis system by 
1.7.6.1
|
