rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

LHCB_2014_I1281685

Charged particle multiplicities and densities in $pp$ collisions at $\sqrt{s} = 7$ TeV
Experiment: LHCB (LHC)
Inspire ID: 1281685
Status: VALIDATED
Authors:
  • Marco Meissner
References: Beams: p+ p+
Beam energies: (3500.0, 3500.0) GeV
Run details:
  • Proton-proton collisions at $\sqrt(s)=7$ TeV, minimum bias data with a fixed interaction rate of 1, no pile-up events. Kinematic cuts are, minimum momentum of 2 GeV/$c$, minimum transverse momentum of 0.2 GeV/$c$ and a pseudorapidity range of $2.0 < \eta < 4.8$.

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  // The hook for the plugin system
1176  RIVET_DECLARE_PLUGIN(LHCB_2014_I1281685);
1177
1178}