GflashNameSpace.h

Go to the documentation of this file.

#ifndef GflashNameSpace_H
#define GflashNameSpace_H

#include "SimGeneral/GFlash/interface/Gflash3Vector.h"

namespace Gflash {

  // unit convention: energy in GeV, length in cm, time in nanosecond

  enum CalorimeterNumber {
    kNULL = -1,  // ECAL Barrel - ESPM
    kESPM,       // ECAL Barrel - ESPM
    kHB,         // HCAL Barrel - HB
    kENCA,       // ECAL Endcap - ENCA
    kHE,         // HCAL Endcap - HE
    kHO,         // HCAL Outer  - HO
    kHF,         // HCAL Forward - HF
    kNumberCalorimeter
  };

  CalorimeterNumber getCalorimeterNumber(const Gflash3Vector &position);
  double rhoBackEB(const Gflash3Vector &position);
  double zBackEE(const Gflash3Vector &position);
  int findShowerType(const Gflash3Vector &position);

  //                                            EB     HB     EE     HE      HO,
  //                                            HF
  const double Zmin[kNumberCalorimeter] = {0.000, 0.000, 317.2, 391.95, 0.000, 1110.0};  // in cm
  const double Zmax[kNumberCalorimeter] = {317.0, 433.2, 399.1, 554.10, 661.0, 1275.0};
  const double Rmin[kNumberCalorimeter] = {123.8, 177.5, 31.6, 31.6, 382.0, 12.5};
  const double Rmax[kNumberCalorimeter] = {177.5, 287.7, 171.1, 263.9, 407.0, 140.0};

  const double EtaMin[kNumberCalorimeter] = {0.000, 0.000, 1.479, 1.479, 0.000, 3.000};
  const double EtaMax[kNumberCalorimeter] = {1.479, 1.479, 3.000, 3.000, 1.262, 5.000};

  const double RFrontCrystalEB = 129.0;  // rho in cm
  const double ZFrontCrystalEE = 320.9;  // z in cm
  const double ROffCrystalEB = 5.2;      // rho in cm
  const double ZOffCrystalEE = 3.7;      // z in cm
  const double LengthCrystalEB = 23.0;
  const double LengthCrystalEE = 22.0;

  // constants needed for GflashHadronShowerModel
  const double energyCutOff = 1.0;

  // cut value for quasi-elastic like interactions on the secondary energy,
  // ratio = (Leading Energy of Secondaries)/(Total Energy)
  const double QuasiElasticLike = 0.95;

  // minimum distance to the back of parameterized envelopes in [cm]
  const double MinDistanceToOut = 10.0;

  // constants needed for GflashHadronShowerProfile

  const double criticalEnergy = 0.0086155;
  const double rMoliere[kNumberCalorimeter] = {2.19, 2.19, 2.19, 2.19, 2.19, 1.72};       // in cm
  const double radLength[kNumberCalorimeter] = {0.89, 1.49, 0.89, 1.49, 1.49, 1.76};      // in cm
  const double intLength[kNumberCalorimeter] = {22.4, 16.42, 22.4, 16.42, 16.42, 16.77};  // in cm
  const double Z[kNumberCalorimeter] = {68.36, 68.36, 68.36, 68.36, 68.36, 55.845};       // mass of molicule

  // The step size of showino along the helix trajectory in cm unit
  const double divisionStep = 1.0;

  // maximum shower depth for the hadronic lateral in interaction length unit
  const double maxShowerDepthforR50 = 6.0;

  // maximum size of the lateral arm for hadronic showers = 12 * interaction
  // length
  const double maxLateralArmforR50 = 197.0;

  // minimum particle energy to trigger HO parameterization in [GeV]
  const double MinEnergyCutOffForHO = 2.5;

  // additional energy scale for the Hcal sensitive detector
  const double scaleSensitiveHB = 0.180;
  //  const double scaleSensitiveHE = 0.220;
  const double scaleSensitiveHE = 0.180;

  // properties for different sub-detectors (calorimeters)

  const double SAMHAD[3][kNumberCalorimeter] = {
      {0.0, 0.89, 0.0, 0.89, 0.89, 0.0}, {0.0, 0.00, 0.0, 0.00, 0.00, 0.0}, {0.0, 0.00, 0.0, 0.00, 0.00, 0.0}};
  const double RLTHAD[kNumberCalorimeter] = {32.7, 23.7, 32.7, 23.7, 23.7, 23.7};

  // parameters for logitudinal/lateral parameterizations and energy responses

  const int NPar = 5;
  const int Nrpar = 4;

  // Pion(-/+) parameters

  const double pion_emscale[2][5] = {{8.4338e-01, 1.8663e-01, -1.1906e+00, 6.8340e-01, -3.1237e-02},
                                     {1.8918e-01, 2.2620e-02, 4.2555e+00, 1.8317e+00, -6.1358e-03}};

  const double pion_hadscale[8][5] = {
      // based on LogNormal of Had Energy - RegionI
      {-4.4752e+00, 3.8163e+00, 2.7962e-01, -7.6887e-01, 9.6721e-01},
      {1.1228e+00, 6.5039e-01, 4.5585e+00, -8.0282e+00, -1.7630e-01},
      {2.8223e+00, 5.9029e+00, -1.5873e-01, -8.7078e-01, 3.8598e-01},
      {-8.2726e-01, 1.0359e-01, 2.2052e+01, 2.9719e+00, 6.4103e-02},
      // LogNormal of Had Energy for mip-like response
      {-8.1787e-02, 8.3243e+00, 3.7488e-03, 3.2653e+00, 9.5662e-01},
      {1.0567e-01, -2.0485e-01, 2.3345e-01, 4.5563e-01, 2.3694e-02},
      // correction of response if ssp is between Ecal and Hcal
      {1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00},
      {5.3155e-02, 9.3601e-02, 8.8143e-01, 2.3251e+00, -2.1999e-02}};

  const double pion_correl_hadem[5] = {-7.3137e-01, 1.2367e-01, -8.5607e-01, 1.0989e+00, -1.2858e-02};

  // pion longitudinal parameters

  const double pion_par[8 * NPar][5] = {
      {4.4600e-01, 9.8040e-02, 1.4319e+00, 2.3920e+00, 5.8971e-02},  // Mean Ecal ShowerType1
      {3.2392e+00, 5.4986e+00, 1.5979e-01, 3.2165e+00, -6.7426e-01},
      {-2.7502e-01, -2.2450e+00, -2.0636e-01, 2.7225e+00, -3.2656e-01},
      {5.3763e-01, 3.1128e-01, 1.5000e+00, 2.7038e+00, 9.7550e-02},
      {5.8269e+00, 5.9630e+00, 2.3568e-01, 3.2485e+00, -1.2201e+00},

      {1.6324e+00, 1.1370e+00, 3.0111e-01, 4.4501e+00, -2.5984e-01},  // Mean Hcal ShowerType1
      {1.2732e+00, -3.8571e-01, 1.2127e+00, 9.7707e-01, -1.4888e-02},
      {-4.0078e-01, -1.0428e+00, 6.6820e-01, 2.3706e-01, -1.3898e-01},
      {1.9612e+00, -1.4451e+00, 5.2179e-01, 7.7221e-01, 4.1455e-02},
      {2.8888e+00, -3.2074e+00, 3.4082e-01, 3.3507e-01, 2.7346e-02},

      {3.5029e+00, 5.5470e+00, 1.5174e-01, 3.9215e+00, -7.5815e-01},  // Mean Hcal ShowerType2
      {1.2035e+00, -4.3918e-01, 9.8431e-01, 5.3485e-01, -3.4481e-02},
      {2.4186e+01, 2.5653e+01, -6.8909e-01, -2.7461e+00, -1.5702e-01},
      {1.0827e+01, 1.0659e+01, -3.3366e-01, -3.1979e+00, 1.1158e-01},
      {1.5010e+01, -1.5988e+01, 2.2580e-01, -4.3233e+00, 1.1448e-01},

      {4.0761e-01, 5.0999e-02, 1.7222e+00, 2.1085e+00, 5.4258e-02},  // Mean Hcal ShowerType3
      {1.0394e+01, -9.5700e+00, 1.1167e+00, -1.6884e+00, -3.3264e-02},
      {-1.1882e+00, 1.6460e-01, -2.5897e+00, 5.4484e-01, -1.5845e-01},
      {2.7227e+01, 2.6921e+01, -3.5664e-01, -4.8584e+00, 1.1698e-01},
      {2.1544e+01, -2.2048e+01, 2.5227e-01, -5.2186e+00, 6.7578e-02},

      {3.6152e-02, 1.5236e-01, -5.4226e-01, 3.4917e+00, 3.1429e-02},  // Sigma Ecal ShowerType1
      {4.9213e-01, 1.0093e-01, -7.8460e-01, 2.2417e+00, -1.9762e-02},
      {3.6555e-01, 2.2955e-02, -3.3972e+01, 2.7112e+00, -2.0412e-02},
      {7.1526e-01, 3.1537e-01, -4.8117e-01, 1.3809e+00, -1.1832e-02},
      {5.3792e-01, 1.7556e-01, -5.9724e-01, 3.9307e-01, -2.1343e-02},

      {2.7662e-01, 5.6892e-02, -1.0211e+00, 1.5592e+00, -2.0545e-02},  // Sigma Hcal ShowerType1
      {2.1221e+00, 1.7139e+00, 2.9685e-01, 3.0293e+00, -4.9449e-01},
      {9.8715e-01, 1.7377e-01, -1.5715e+00, 1.1948e+00, -7.4157e-02},
      {1.2424e+00, 5.8446e-01, -7.3994e-01, 1.6481e+00, -6.4211e-02},
      {1.2797e+00, 9.5722e-01, -4.1961e-01, 1.5458e+00, -1.9276e-02},

      {5.6503e-01, 3.2484e-01, -6.8870e-01, -9.6941e-01, -2.4937e-02},  // Sigma Hcal ShowerType2
      {7.2561e-01, 5.3369e-02, 9.1553e+00, 1.9652e+00, -5.4515e-02},
      {3.3058e+01, 8.7333e+01, -2.4388e-02, -1.5815e+01, 1.6538e+00},
      {1.1171e+00, 6.8039e-01, -6.0860e-01, 1.2764e+00, -2.1545e-02},
      {2.2385e+00, 2.4309e+00, -2.4251e-01, -8.3314e-01, 4.4941e-02},

      {2.5595e-01, 2.0105e-02, -1.7804e+00, 2.6746e+00, -2.6966e-02},  // Sigma Hcal ShowerType3
      {3.0882e-01, 3.1336e-01, -5.1932e-01, 3.3669e+00, 6.9642e-02},
      {6.1781e-01, 1.1997e-01, 9.2596e-01, -1.3085e+01, -5.1192e-02},
      {1.4680e+00, -1.1644e+00, 4.1119e-01, -2.3604e-01, 3.1277e-03},
      {3.9081e+00, 4.2913e+00, -1.7321e-01, -3.9169e+00, 6.0504e-02}};

  const double pion_rho[8 * NPar][5] = {
      {1.2318e-01, 3.0000e-01, 8.8725e-01, 1.6928e+00, 0.0000e+00},  // Ecal ShowerType1
      {4.1546e-01, 1.6375e-01, 2.0719e+00, 1.2931e+00, -2.3695e-03},
      {-9.4588e-01, 1.7680e+00, 1.1180e+00, -1.7623e+00, 6.6812e-03},
      {-1.8668e-01, 1.5000e-01, -5.0000e+00, 2.4000e+00, 6.0000e-02},
      {1.1945e-01, 1.5393e-01, -1.6983e+00, 2.2457e+00, 5.1813e-02},
      {3.4322e-02, 2.0564e-01, -1.4297e+00, 2.1601e+00, 4.0797e-02},
      {-1.0894e-01, 1.4305e-01, -1.5543e+00, 2.5491e+00, 5.6509e-02},
      {8.1763e-02, 8.0722e-02, -5.8023e+00, 2.4338e+00, 5.4228e-02},
      {9.4796e-02, 1.5509e-01, -1.8610e+00, 2.3388e+00, 3.9771e-02},
      {7.8380e-01, 1.4165e-01, 7.7385e-01, 1.9248e+00, 4.2539e-03},

      {4.2393e-01, 2.7909e-01, 9.2577e-01, 2.7916e+00, -6.8012e-02},  // Hcal ShowerType1
      {5.6826e-01, 2.9822e-01, 8.1516e-01, 2.3029e+00, -9.2674e-02},
      {8.0676e-01, 1.0045e+00, -4.2949e-03, -1.1090e+01, -9.6997e-04},
      {-1.8671e-01, 7.4935e-01, 3.6716e-04, 4.3540e+00, 8.4585e-02},
      {3.7378e-01, 1.8905e+00, 1.3206e-01, 8.3958e-01, -1.4278e-01},
      {2.4541e-01, 1.4579e-01, 1.3386e+00, 1.1209e+00, 1.5902e-03},
      {4.4394e-02, 1.1718e-01, 2.0436e+00, 1.0125e+00, -2.0101e-02},
      {1.5076e-01, 8.8556e-02, -1.8958e+00, 2.2495e+00, 4.6895e-02},
      {5.0817e-01, 1.2550e-01, 2.3684e+00, 1.0669e+00, -1.3888e-02},
      {8.3803e-01, 5.6501e-02, -1.2531e+00, 2.7326e+00, -4.1174e-02},

      {3.0277e-01, 2.1400e-01, 1.2867e+00, 2.3086e+00, -8.7281e-02},  // Hcal ShowerType2
      {5.2508e-01, 2.6391e-01, 1.3357e+00, 1.9172e+00, -9.2390e-02},
      {7.5387e-01, -3.2142e-02, 4.7437e+00, 8.5346e-01, -1.1566e-02},
      {-7.5633e-02, 8.8550e-02, 2.2340e+00, 9.7441e-01, 3.2514e-02},
      {2.0314e-01, 1.2932e-01, 9.9352e-01, 6.6276e-01, 3.2726e-02},
      {-2.6996e+00, 3.9944e+00, -2.3836e-01, 3.4896e+00, 8.5793e-01},
      {1.2888e-01, 2.4009e-01, 9.8654e-01, 9.8778e-01, -4.3265e-02},
      {2.2281e-01, 8.4319e-02, 4.8402e-04, -2.0111e-01, 1.9263e-02},
      {-9.3521e+00, 1.0283e+01, 3.5236e-01, -4.8949e+00, -7.7675e-02},
      {7.6844e-01, -9.0198e-02, 1.0364e+00, 2.3443e+00, -1.4107e-02},

      {3.0154e-01, 3.6790e-01, 7.5730e-01, 1.7676e+00, -9.7870e-02},  // Hcal ShowerType3
      {4.2110e-01, 3.6328e-01, 9.7741e-01, 1.1211e+00, -6.4038e-02},
      {7.6092e-01, 3.5780e-02, 9.7086e-01, 1.9511e+00, -1.2363e-02},
      {-1.4874e+00, 2.4273e+00, -1.0552e-01, 7.3591e+00, 1.7302e-01},
      {-1.6966e+00, 2.9215e+00, -2.2195e-01, 3.3990e+00, 6.0529e-01},
      {1.6144e-01, 2.1615e-01, -7.4824e-01, 3.1003e+00, 0.0000e+00},
      {-2.2588e+00, 3.9894e+01, 5.5166e-02, -1.1180e+00, -2.1166e+00},
      {6.8954e-01, 1.1032e+00, 2.4371e-01, 1.8565e+00, -1.9127e-01},
      {4.7802e-02, 3.0961e+00, 1.4566e-01, -1.0462e+00, -3.6505e-01},
      {6.3815e-01, -5.1494e-01, 3.3634e-01, 1.4517e+00, 8.8357e-02}};

  // pion lateral parameters

  const double pion_rpar[4 * Nrpar][5] = {
      {1.7239e+00, 9.3646e-01, -7.6784e-01, 1.6319e+00, 0.0000e+00},  // Ecal for showerType = 1
      {8.0152e+00, 1.3573e+00, -2.2490e+00, 2.0799e+00, -9.1649e-01},
      {7.3494e-01, 8.7912e-02, -1.9433e+00, 1.2044e+00, -3.4421e-02},
      {-2.2025e-01, 7.3525e-02, 3.1836e+00, 1.4077e+00, 2.3818e-02},

      {2.0617e+01, 2.1354e+00, -1.3979e+00, 2.4424e+00, -2.0928e+00},  // Hcal for showerType = 1
      {1.4732e+02, 1.4804e+02, -1.0534e+00, -2.1180e+00, 1.9513e-02},
      {5.2939e-01, 2.3488e-01, 4.2874e-01, 1.4435e+00, -6.7974e-02},
      {-2.2405e-02, 1.7781e-01, 5.8629e-03, -2.6258e+01, 2.0595e-03},

      {1.3184e+01, 5.3308e+00, -6.5926e-01, 1.4287e+00, -9.1948e-01},  // Hcal for showerType = 2
      {6.0658e+00, 7.6593e+00, -6.9664e-01, -1.7824e+00, 1.5619e-01},
      {5.5211e-01, 0.0000e+00, 0.0000e+00, 0.0000e+00, -2.9011e-02},
      {2.1551e-02, -5.4231e-02, 1.6615e+00, 2.7025e-01, 1.1807e-02},

      {9.0230e-01, 8.7071e-01, -5.9984e-01, 2.7973e+00, 2.4553e-01},  // Hcal for showerType = 3
      {7.5508e+00, 3.6556e+00, -8.3028e-01, 1.1326e+00, -6.4124e-01},
      {2.5737e+01, 2.5165e+01, -5.6146e-01, -4.0345e+00, -2.9592e-02},
      {-3.1200e-01, 2.0000e-01, 2.8915e+00, 7.5519e-01, 1.5230e-02}};

  // Proton paramters

  const double pro_emscale[2][5] = {{6.2614e-01, 1.7633e-01, -1.2677e+00, 1.3193e+00, 2.4973e-03},
                                    {1.6413e-01, 5.2570e-03, 1.8876e+02, 1.2887e+00, -4.0266e-03}};

  const double pro_hadscale[8][5] = {
      // based on LogNormal of Had Energy
      {-1.2767e+02, 1.2691e+02, 2.6892e-01, -8.4621e+00, 9.7956e-01},
      {1.2283e+00, 1.6307e-01, -3.6863e+00, 1.8677e+00, -4.5990e-02},
      {1.8917e+00, 2.0947e+00, -3.6371e-01, 3.7861e-01, 8.8760e-02},
      {-7.6830e-01, 3.1319e-01, 8.2844e-01, 1.5865e+00, 2.4735e-02},
      // LogNormal of Had Energy for mip-like response
      {-1.0927e+01, 1.0727e+01, 2.9922e+00, -8.2073e-01, 9.8722e-01},
      {7.2797e-02, 2.3601e-02, -1.4900e+00, 1.8989e+00, 3.4294e-04},
      // correction of response if ssp is between Ecal and Hcal
      {1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00},
      {4.7858e-01, 5.8552e-01, 3.9112e-01, 2.6945e+00, -1.7791e-01}};

  const double pro_correl_hadem[5] = {-6.7711e-01, 2.5760e-01, -4.9826e-01, 3.4791e-01, 0.0000e+00};

  // proton longitudinal parameters

  const double pro_par[8 * NPar][5] = {
      {3.5805e-01, 8.2073e-02, 2.4040e+00, 2.8935e+00, 7.1200e-02},  // Mean Ecal ShowerType1
      {8.1351e-01, 9.3288e-02, 4.5959e+00, 4.1930e+00, 8.2825e-02},
      {-1.2342e+00, 1.3826e-01, 8.1966e-01, 3.7141e+00, 2.7493e-02},
      {5.6140e-01, 2.1720e-01, 2.0061e+00, 2.8141e+00, 1.1165e-01},
      {3.9409e+00, 2.6083e+00, 3.0382e-01, 3.2679e+00, -6.3072e-01},

      {1.5450e+00, 1.1047e+00, 2.9732e-01, 4.5325e+00, -2.4096e-01},  // Mean Hcal ShowerType1
      {1.5098e+00, -6.8331e-01, 8.7004e-01, 5.0915e-01, -9.3536e-03},
      {6.9904e+00, -8.9057e+00, 3.8745e-01, -2.9089e+00, -6.0353e-02},
      {1.8465e+00, -1.1548e+00, 7.1662e-01, 9.2467e-01, 7.0536e-03},
      {2.1807e+00, -1.4464e+00, 7.2033e-01, 1.0696e+00, -1.3923e-01},

      {3.9185e+00, 6.1651e+00, 1.4738e-01, 4.2079e+00, -8.1423e-01},  // Mean Hcal ShowerType2
      {8.3119e+00, 7.5628e+00, -4.9561e-01, -2.9352e+00, -3.0061e-02},
      {2.7150e+00, 4.5234e+00, -4.6494e-01, -1.9846e+00, -9.9347e-02},
      {2.2762e+00, 1.7668e+00, -5.3421e-01, -6.4613e-03, 5.3533e-02},
      {2.6187e+00, -2.3479e+00, 4.9602e-01, 1.8344e-01, -5.3695e-02},

      {1.9666e+00, 3.0810e+00, 1.7874e-01, 3.2441e+00, -4.4566e-01},  // Mean Hcal ShowerType3
      {1.6223e+01, 1.7175e+01, -1.3214e-01, -1.0103e+01, 1.9030e-01},
      {8.1636e+01, 8.3368e+01, -5.0505e-01, -5.2201e+00, -1.0340e-01},
      {5.5674e+01, 5.5241e+01, -3.8856e-01, -5.5563e+00, 9.2886e-02},
      {6.1932e+01, 6.2078e+01, -3.2123e-01, -5.9991e+00, 2.2640e-02},

      {7.9984e-02, 1.0177e-01, -4.9195e-01, 3.6549e+00, 1.3986e-02},  // Sigma Ecal ShowerType1
      {4.6155e-01, 1.1266e-01, -1.1120e+00, 2.8973e+00, -7.7394e-03},
      {3.6587e-01, 3.0463e-02, -1.8552e+00, 2.9503e+00, -1.9858e-02},
      {7.0777e-01, 1.5344e-01, -9.7573e-01, 2.2658e+00, -4.3013e-02},
      {3.7488e-01, 3.3388e-01, -3.2414e-01, 1.7099e+00, 2.7456e-02},

      {2.3330e-01, 8.3848e-02, -1.4841e+00, 1.6539e+00, -8.9738e-03},  // Sigma Hcal ShowerType1
      {9.0862e-01, 1.9284e-01, -2.9385e+00, 5.5159e-01, -4.4707e-02},
      {1.0615e+00, 4.0817e-01, -9.2562e-01, 7.7349e-01, -5.2196e-02},
      {1.3044e+00, 6.9523e-01, -8.0370e-01, 1.3130e+00, -6.0052e-02},
      {1.3648e+00, 6.7204e-01, -7.7137e-01, 1.3261e+00, -8.1459e-02},

      {2.3378e-01, 8.4699e-02, -1.2727e+00, 1.3170e+00, -9.0939e-03},  // Sigma Hcal ShowerType2
      {6.7311e-01, 2.6483e-02, -2.3889e+01, 3.3071e+00, -3.8701e-02},
      {1.9157e+01, 1.8483e+01, -6.2208e-01, -3.5409e+00, -6.3584e-02},
      {1.2720e+00, 7.7241e-01, -6.6997e-01, 8.5122e-01, -3.7883e-02},
      {3.2209e+00, 2.8976e+00, -3.4428e-01, -1.5930e+00, -2.2328e-02},

      {2.3064e+00, 2.1994e+00, -2.4570e-01, -6.2275e+00, -6.5378e-03},  // Sigma Hcal ShowerType3
      {6.2645e-01, 1.0763e-01, -2.4655e+00, 3.0694e+00, -2.2801e-02},
      {3.9143e+00, -3.1821e+00, 2.2346e+00, -8.5515e-01, -7.0756e-02},
      {8.8440e+00, 8.6952e+00, -2.8769e-01, -4.4932e+00, 2.0374e-02},
      {7.1745e+00, 7.0713e+00, -2.5061e-01, -4.6686e+00, 1.0573e-02}};

  const double pro_rho[8 * NPar][5] = {
      {1.9717e+00, 6.6780e+00, 1.4425e-01, 2.2917e+00, -8.1369e-01},  // Ecal ShowerType1
      {-6.3565e+00, 7.0824e+00, 1.9663e-01, -7.8033e+00, -2.3854e-02},
      {8.0629e-01, 6.5672e-02, 1.4501e+00, 3.3136e-01, -5.6348e-03},
      {-1.6000e-01, 2.3904e-01, 1.1818e+00, 6.1474e-01, -1.0088e-02},
      {1.0966e-01, 1.8834e-01, -1.4410e+00, 2.6445e+00, 9.0217e-02},
      {1.0565e-02, 2.4379e-01, -1.3677e+00, 2.7338e+00, 7.6137e-02},
      {-3.1898e-02, 1.1508e-01, 1.5916e+00, 6.2431e-01, -9.8921e-03},
      {1.5641e-01, 1.8429e-01, 1.8207e+00, 4.9841e-01, 1.3907e-03},
      {-1.5532e-01, 3.4349e-01, -8.5195e-01, 2.8581e+00, 1.4698e-01},
      {6.9661e-01, 2.0247e-01, 7.0948e-01, 6.1266e-01, 1.0035e-02},

      {2.8900e-01, 1.9972e-01, 1.3997e+00, 2.8253e+00, -3.0914e-02},  // Hcal ShowerType1
      {3.2869e-01, 8.6843e-02, 6.0682e+01, 2.1983e+00, -1.0146e-02},
      {7.8200e-01, 6.8963e-02, -1.9113e+00, 7.5866e-01, 5.9595e-03},
      {9.5356e-01, 1.1469e+00, 3.0959e-01, 1.0238e+01, 7.6437e-02},
      {3.4378e-01, 1.6417e-01, 1.2061e+00, 1.6405e+00, 2.3441e-02},
      {2.8005e-01, 2.0804e-01, 1.3118e+00, 1.1119e+00, -1.9884e-02},
      {-2.7624e-02, 1.0661e-01, 2.7942e+00, 5.8531e-01, -8.9599e-03},
      {2.3749e-01, 2.6046e-02, 1.8467e+01, 6.2294e-01, 9.8564e-03},
      {4.6006e-01, 2.0859e-01, 1.2562e+00, 7.6745e-01, -1.9148e-02},
      {8.1573e-01, 6.8433e-02, -2.1358e+00, 2.2568e+00, -3.8336e-02},

      {1.1195e-01, 1.1879e-01, 2.9551e+00, 1.9034e+00, -2.0542e-02},  // Hcal ShowerType2
      {4.1285e-01, 4.6288e-01, 5.0506e-01, 1.4100e+00, -1.2918e-01},
      {7.7938e-01, -4.7534e-02, 1.8256e+01, 6.8910e-02, -8.8481e-03},
      {2.1819e-01, 2.1512e-01, 7.0238e-01, 3.6517e+00, -2.4946e-02},
      {3.5413e-01, 1.8842e-01, 1.0927e+00, 1.2173e+00, 3.4229e-04},
      {-1.9019e+00, 2.8987e+00, -2.5328e-01, 3.4091e+00, 6.5499e-01},
      {-7.1433e-02, 7.6591e-02, -6.8698e+00, 2.1165e+00, 5.6828e-02},
      {2.2651e-01, 7.3121e-01, 6.9897e-04, 3.5656e+00, 2.7391e-02},
      {-9.1867e+00, 1.0087e+01, 3.4287e-01, -5.2458e+00, -7.1044e-02},
      {7.7621e-01, -9.5476e-02, 1.0892e+00, 1.8212e+00, -1.6587e-02},

      {1.4812e-01, 3.7207e-01, 6.7394e-01, 1.6932e+00, -8.9577e-02},  // Hcal ShowerType3
      {-1.0980e+01, 1.1814e+01, 2.9231e-01, -5.5518e+00, -9.7033e-02},
      {7.8865e-01, -4.6499e-02, 1.0071e+01, -1.7465e-01, -7.8028e-03},
      {-4.8628e-03, 6.8062e-02, 4.0379e-04, -2.9541e+00, 1.4925e-02},
      {3.2583e-01, 2.2048e-01, 1.3796e+00, 7.7019e-01, -1.9870e-02},
      {-4.1067e+00, 8.5741e+00, -1.8161e-01, 3.0661e+00, 1.4121e+00},
      {1.3314e-01, 1.3727e-01, -2.5181e-04, -4.8352e-01, -7.1482e-03},
      {1.8448e-01, -1.0892e-01, 2.0336e-04, -3.5803e+01, 6.2453e-02},
      {-4.3394e+00, 6.7933e+00, 1.2905e-01, -6.7132e+00, -2.7930e-01},
      {7.6150e-01, -1.3613e-01, 9.9347e-01, 1.1661e+00, 5.4822e-03}};

  // proton lateral parameters

  const double pro_rpar[4 * Nrpar][5] = {
      {2.0473e+01, 1.8293e+01, -2.7258e-03, -1.4304e+03, -2.7740e-01},  // Ecal for showerType = 1
      {4.3030e+00, 3.1796e+01, -1.4979e-01, 1.2390e+00, 2.8661e+00},
      {-7.1506e-01, 3.4815e+00, -1.8282e-01, 2.4101e+00, 5.4402e-01},
      {2.5011e+00, 5.3549e+00, 1.7772e-01, 2.9511e+00, -8.8951e-01},

      {2.0252e+01, 1.6754e+00, -3.3596e+00, 2.1412e+00, -2.1051e+00},  // Hcal for showerType = 1
      {-1.6270e+00, 5.0493e+00, -1.6916e-01, 1.5961e+00, 7.2367e-01},
      {5.0396e-01, 7.1958e-02, 1.4915e+00, 1.4498e+00, -4.0554e-02},
      {-1.3402e-01, 1.3657e-01, 2.3093e+00, -3.1556e+00, 5.2742e-03},

      {1.2771e+01, 2.8591e+00, -7.5931e-01, 1.5374e+00, -1.2623e+00},  // Hcal for showerType = 2
      {-4.1051e+00, 3.1735e+00, -3.5820e-01, 3.0133e+00, 1.0538e+00},
      {4.9758e-01, 5.9756e-02, -1.9593e+00, 2.9959e+00, -1.3451e-02},
      {2.0201e-01, 1.0967e-01, 1.2310e+00, 3.3239e+00, -6.2154e-02},

      {-4.8143e+00, 1.0479e+01, -2.1974e-01, 3.0933e+00, 2.0855e+00},  // Hcal for showerType = 3
      {2.5276e+01, 2.1056e+01, -4.4710e-01, -1.7116e+00, -6.7854e-01},
      {-5.8886e-01, 2.6766e+00, -1.9974e-01, 2.6381e+00, 4.3146e-01},
      {-1.7710e-02, 9.4732e-02, 4.2625e+00, 3.3545e+00, -2.9552e-02}};

  // Anti-Proton paramters

  const double pbar_emscale[2][5] = {{1.9246e+01, 1.8870e+01, -3.3127e-01, -4.8509e+00, -2.3785e-02},
                                     {3.5744e+00, 3.3661e+00, -4.9667e-01, -4.1388e+00, -8.8824e-03}};

  const double pbar_hadscale[8][5] = {
      // based on LogNormal of Had Energy
      {-1.8857e+00, 9.4549e-01, 1.0173e+00, 2.3654e+00, 1.0381e+00},
      {1.4276e+00, 8.3835e-01, -9.1466e-01, 2.3945e+00, 0.0000e+00},
      {1.0005e+00, 5.0454e-01, -9.6777e-01, 2.4216e+00, -3.0214e-02},
      {-7.1494e-01, 3.7825e-01, 1.1627e+00, 2.7411e+00, 5.5459e-03},
      // LogNormal of Had Energy for mip-like response
      {1.5770e+00, 1.1215e+00, 4.2471e-01, 3.8321e+00, 5.2852e-01},
      {-6.5624e-02, 2.7895e-01, -3.2639e-01, 2.3605e+00, 6.3871e-02},
      // correction of response if ssp is between Ecal and Hcal
      {1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00},
      {9.1334e-02, 8.4321e-02, 9.1042e-01, 2.7686e+00, -3.0344e-02}};

  const double pbar_correl_hadem[5] = {-5.8253e-01, 3.6344e-01, -4.6371e-01, 2.8284e-01, 0.0000e+00};

  // anti-proton longitudinal parameters

  const double pbar_par[8 * NPar][5] = {
      {2.2962e+00, 1.8237e+00, 1.5856e-01, 1.1724e+01, -2.5137e-02},  // Mean Ecal ShowerType1
      {5.1688e-01, 2.9636e-02, 3.9787e+01, 3.1512e+00, 1.7428e-01},
      {-1.3058e+00, 4.9530e-02, 1.5425e+00, 4.1539e+00, 5.7473e-02},
      {2.9921e-01, 2.0357e-01, 2.9853e+00, 1.8306e+00, 1.5198e-01},
      {3.8358e+00, 5.0343e+00, 1.3900e-01, 3.2366e+00, -6.1138e-01},

      {9.2918e-01, 3.4905e-01, 4.4705e-01, 4.3199e+00, -1.0055e-01},  // Mean Hcal ShowerType1
      {1.1081e+00, -3.4156e-01, 1.2140e+00, 1.0438e+00, 2.1996e-02},
      {2.5596e+01, -4.7613e+01, 3.8236e-02, -1.6258e+01, 8.3889e-01},
      {1.4731e+00, -8.7282e-01, 7.6302e-01, 1.4135e+00, 4.0020e-02},
      {1.5854e+00, -1.5158e+00, 5.3554e-01, 1.7953e+00, -2.1115e-02},

      {1.6397e+00, 1.4348e+00, 1.8597e-01, 5.0391e+00, -2.1112e-01},  // Mean Hcal ShowerType2
      {1.9561e+00, 9.1382e-01, 5.2437e-01, 3.2069e+00, -3.6565e-01},
      {8.4435e-01, 2.0674e+00, -4.2496e-01, -3.2228e+00, -2.0347e-01},
      {1.3139e+00, 6.8898e-01, -9.0400e-01, 1.1612e+00, 5.2261e-02},
      {1.3894e+00, -2.3487e+00, 3.5770e-01, 1.2953e+00, 1.4271e-01},

      {6.4737e-01, 7.3992e-02, 9.2454e-01, 4.8070e+00, -3.3204e-03},  // Mean Hcal ShowerType3
      {9.6156e-01, 1.5021e-01, 2.1087e+00, 2.7135e+00, -5.6608e-02},
      {-2.9797e+00, 3.7631e+00, -1.8261e-01, 2.8053e+00, 4.3358e-01},
      {1.0332e+00, 4.4820e-01, -1.2826e+00, 9.7747e-01, 9.4973e-02},
      {1.6588e+00, 3.6401e+00, -2.4077e-01, 3.4642e-01, 2.7807e-01},

      {7.9348e-01, 6.0791e-01, -5.2480e-01, -4.1979e+00, -9.7594e-03},  // Sigma Ecal ShowerType1
      {6.7623e-01, 1.1835e-01, 2.5074e+01, 2.9602e+00, -7.6399e-02},
      {4.3953e-01, 5.0370e-02, 4.3043e+01, 2.9560e+00, -3.5965e-02},
      {6.7091e-01, 1.1626e-01, -5.7986e-01, 2.9029e-01, -2.8586e-02},
      {4.2821e-01, 8.0538e-02, -9.3285e-01, 1.2197e+00, -1.3245e-02},

      {2.3715e-01, 6.1849e-02, -7.9327e-01, 2.4567e+00, -1.5569e-02},  // Sigma Hcal ShowerType1
      {7.7529e-01, 8.4253e-02, -2.8453e+00, 6.6960e-01, -2.3133e-02},
      {9.7261e-01, 1.1777e-01, -1.2671e+00, 1.3352e+00, -9.0527e-02},
      {1.1170e+00, 5.7366e-01, -7.9215e-01, 1.8696e+00, -3.2848e-02},
      {1.0912e+00, 5.3449e-01, -8.3799e-01, 2.2648e+00, -5.3589e-02},

      {2.0364e-01, 1.1109e-01, -5.4911e-01, 2.2371e+00, -2.7723e-03},  // Sigma Hcal ShowerType2
      {5.0058e-01, 8.3324e-02, -1.5426e+00, 4.3092e+00, 1.7710e-02},
      {4.8258e+00, -5.3363e+00, -6.9399e-02, 1.3639e+01, -3.1232e-01},
      {1.0654e+00, 4.3313e-01, -8.1864e-01, 1.5794e+00, -5.0301e-02},
      {9.7937e-01, 4.9011e-01, -7.5635e-01, 2.0913e+00, -4.3251e-02},

      {2.4410e-01, 6.9464e-02, -5.7121e-01, 1.5400e+00, -1.8889e-02},  // Sigma Hcal ShowerType3
      {6.8356e-01, 4.4101e-02, 2.0167e+01, 2.9570e+00, -4.3937e-02},
      {2.9627e-01, -4.3021e-01, 1.3628e-01, 1.3540e+01, -5.0582e-02},
      {9.2673e-01, 3.5636e-01, -6.8759e-01, 1.1528e+00, -3.6487e-02},
      {8.1179e-01, 5.0238e-01, -5.7945e-01, 1.8664e+00, -9.9021e-03}};

  const double pbar_rho[8 * NPar][5] = {
      {1.6115e-01, 7.7749e-02, 1.3576e+01, 6.4878e-01, -3.2481e-03},  // Ecal ShowerType1
      {4.1746e-01, 6.4980e-02, -1.8986e+01, 2.9812e+00, 8.2744e-03},
      {8.3730e-01, 5.2519e-02, 1.9080e+01, 2.9734e+00, -7.2551e-04},
      {4.7186e-01, 2.3447e+00, 1.5313e-01, 2.1195e+00, -2.7598e-01},
      {2.6387e-01, 1.2494e-01, 2.0479e+01, 2.9541e+00, 3.9418e-03},
      {1.4487e-01, 1.7347e-01, 2.0829e+01, 2.9708e+00, -1.7575e-02},
      {-3.2116e-02, 1.0593e-01, -3.6256e+00, 2.0405e+00, 3.9701e-02},
      {3.2368e-01, 3.4897e-01, 4.0152e-01, 1.9683e+00, -4.3362e-02},
      {3.3095e-01, 8.8705e-02, 3.0619e+01, 3.2302e+00, -2.6705e-02},
      {6.8646e-01, 1.3220e-01, 1.8172e+00, 1.8252e+00, 2.5609e-02},

      {2.5579e-01, 1.3158e-01, 2.0117e+00, 2.6834e+00, -1.8803e-02},  // Hcal ShowerType1
      {5.4665e-01, 3.5282e-01, 7.0197e-01, 1.8240e+00, -1.2513e-01},
      {7.8407e-01, 4.6267e-02, -1.9313e+00, 1.2055e+00, 3.1770e-03},
      {1.9774e-01, 1.4336e+00, 1.0833e-01, 2.5249e+00, -6.1467e-02},
      {4.8669e-01, 2.5554e-01, 7.8888e-01, 2.7417e+00, 0.0000e+00},
      {2.5270e-01, 1.3412e-01, 7.7256e-01, 1.3281e+00, 1.2302e-02},
      {4.3116e-02, 1.3977e-01, 1.2806e+00, 8.8424e-01, -2.2153e-02},
      {-1.1513e-01, 1.1712e+00, -1.9271e-01, 1.8406e+00, 2.3827e-01},
      {4.8509e-01, 1.9358e-01, 6.8277e-01, 7.8174e-01, -1.2847e-02},
      {6.8735e-01, 1.7148e-01, -8.8434e-01, 2.6390e+00, 1.8336e-02},

      {2.6738e-01, 1.6698e-01, 1.6242e+00, 2.3194e+00, -5.4439e-02},  // Hcal ShowerType2
      {5.1565e-01, 3.1755e-01, 7.6866e-01, 1.4314e+00, -1.2355e-01},
      {1.4488e+00, -1.0738e+00, 1.5373e-01, -4.4945e+00, 5.0825e-02},
      {-5.0531e-02, 5.4899e-01, 6.4605e-04, 6.2060e+00, 5.0480e-02},
      {5.0045e-01, 2.8609e-01, 7.7954e-01, 2.4851e+00, -1.7178e-02},
      {2.4183e-01, 1.1418e-01, 1.1548e+00, 6.8739e-01, 1.4071e-03},
      {-2.7267e-01, 2.9060e-01, -9.2639e-01, 3.3275e+00, 1.3254e-01},
      {3.6510e-01, 6.8755e-02, 1.7529e+00, 3.3932e+00, 3.5968e-03},
      {-2.2256e+00, 4.0953e+00, 1.2237e-01, -6.3486e+00, -1.6487e-01},
      {6.7014e-01, -1.8220e-01, 8.3413e-01, 2.3705e+00, 3.1919e-02},

      {-9.2169e-02, 2.7435e-01, -1.4954e+00, 3.4936e+00, 5.4085e-02},  // Hcal ShowerType3
      {-5.7426e-02, 4.2696e-01, -9.1241e-01, 3.5629e+00, 1.0904e-01},
      {7.4201e-01, 4.7628e-02, 1.5747e+00, -4.6694e+00, -1.4296e-02},
      {1.8949e-02, 1.0836e-01, -1.2242e+00, 2.7149e+00, 1.3031e-02},
      {-1.5217e+00, 2.4872e+00, 8.1875e-02, -1.0666e+01, -1.4661e-03},
      {2.2412e-01, 9.4590e-02, -2.4942e+00, 1.8584e+00, 1.6461e-02},
      {-5.4793e-01, 8.2034e-01, -4.2542e-01, 3.6711e+00, 1.9032e-01},
      {-2.5079e-01, -8.8412e-01, 1.5082e-01, 5.3762e+00, 1.3426e-01},
      {-7.3299e+00, 8.1099e+00, 4.0558e-01, -5.0694e+00, -5.9030e-02},
      {5.6633e-01, -4.7278e-01, 4.1026e-01, 1.7565e+00, 1.1369e-01}};

  // anti-proton lateral parameters
  const double pbar_rpar[4 * Nrpar][5] = {
      {1.8776e+00, 2.8426e-01, 2.3530e+01, 3.0047e+00, -1.7332e-01},  // Ecal for showerType = 1
      {9.1705e+00, 4.0028e+00, -4.4536e-01, 6.4847e-01, -5.9678e-01},
      {5.9148e-01, 7.5552e-02, -3.0172e+01, 1.3417e+00, 2.0513e-02},
      {-1.3812e-01, 8.0542e-02, 5.5994e+01, 1.2251e+00, -1.6587e-02},

      {1.7618e+01, 4.3991e+00, -9.7543e-01, 2.7789e+00, -1.2115e+00},  // Hcal for showerType = 1
      {-7.6558e-02, 1.3021e+00, -1.4519e+00, 1.3966e+00, 1.6778e-01},
      {5.9899e-01, 1.6407e-01, 8.4000e-01, 2.1795e+00, -7.0504e-02},
      {-1.4042e-01, 1.3416e-01, 2.3339e+00, -2.7260e+00, 5.8295e-03},

      {1.2340e+01, 5.3893e+00, -4.9219e-01, 1.5722e+00, -6.8433e-01},  // Hcal for showerType = 2
      {1.0814e-01, 1.6610e+00, -1.1174e+00, -5.0663e+00, 1.9605e-01},
      {4.9994e-01, 0.0000e+00, 0.0000e+00, 0.0000e+00, -2.0081e-02},
      {4.6168e-01, 3.5308e-01, 5.4512e-01, 3.2036e+00, -1.4086e-01},

      {-1.6306e+00, 3.5609e+00, -1.8903e-01, 6.6580e+00, 4.6936e-01},  // Hcal for showerType = 3
      {6.5560e+00, 2.3432e+00, -9.8700e-01, 1.3115e+00, -6.5919e-01},
      {7.4461e-01, 6.6885e-02, -8.9140e-01, -1.2050e+01, -4.5248e-02},
      {-2.4275e-01, 1.6289e-01, -2.9514e-04, -1.1208e+00, 4.3240e-02}};

  // Kaon-Minus paramters

  const double kminus_emscale[2][5] = {{8.3935e-01, 2.4180e-01, -1.4484e+00, 1.0397e+00, -2.2546e-02},
                                       {2.3334e-01, 3.6810e-02, 1.0247e+01, 1.8556e+00, -1.9288e-02}};

  const double kminus_hadscale[8][5] = {
      // based on LogNormal of Had Energy
      {-5.9393e+01, 6.4013e+01, 1.0843e-01, -1.2441e+01, 4.2589e-01},
      {1.3649e+01, 1.2266e+01, -4.9381e-01, -3.6921e+00, -1.1441e-01},
      {1.2102e+00, 4.7458e-01, -9.4128e-01, 2.1482e+00, -5.1031e-02},
      {-6.2018e-01, 2.0288e-01, 1.1271e+00, 2.8535e+00, 0.0000e+00},
      // LogNormal of Had Energy for mip-like response
      {1.7106e-01, 1.3766e-01, 7.3972e-01, 4.1351e+00, 9.0284e-01},
      {5.7520e-01, -1.2499e+00, 8.0351e-02, -4.7481e+00, 6.0937e-02},
      // correction of response if ssp is between Ecal and Hcal
      {1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00},
      {4.2896e-02, 6.7448e-02, 7.4893e-01, 2.8837e+00, -1.1509e-02}};

  const double kminus_correl_hadem[5] = {-7.3480e-01, 1.8541e-01, -7.6211e-01, 9.4878e-01, 2.5519e-03};

  // kaon-minus longitudinal parameters

  const double kminus_par[8 * NPar][5] = {
      {4.5986e-01, 8.0360e-02, 1.9103e+00, 2.7576e+00, 5.8024e-02},  // Mean Ecal ShowerType1
      {5.6306e-01, 8.6175e-02, 7.5699e+00, 1.9142e+00, 1.4928e-01},
      {-1.4029e+00, 8.6442e-02, 4.7447e+00, 2.0026e+00, 7.8132e-02},
      {5.1022e-01, 3.1763e-01, 1.4969e+00, 2.4374e+00, 1.0470e-01},
      {4.4840e+00, 4.0858e+00, 2.4837e-01, 3.0941e+00, -8.4496e-01},

      {2.0880e+00, 1.9024e+00, 2.2650e-01, 4.5397e+00, -3.5337e-01},  // Mean Hcal ShowerType1
      {1.2654e+00, -3.2827e-01, 1.4255e+00, 9.2278e-01, -2.8213e-02},
      {-8.5073e-01, 4.8949e-01, -1.1802e+00, 8.3284e-01, -1.6258e-01},
      {2.4486e+00, -2.0627e+00, 4.1771e-01, 2.8534e-02, 5.9279e-02},
      {1.0307e+01, -1.2125e+01, 1.8949e-01, -3.5864e+00, 1.9567e-01},

      {3.0833e+00, 4.5830e+00, 1.6380e-01, 3.8496e+00, -6.6617e-01},  // Mean Hcal ShowerType2
      {1.5334e+01, -1.4598e+01, 4.5394e-01, -4.3694e+00, -3.0904e-02},
      {1.1765e+01, 1.3252e+01, -6.1629e-01, -2.8795e+00, -1.5726e-01},
      {1.1343e+01, 1.1615e+01, -2.4147e-01, -4.5126e+00, 1.6850e-01},
      {1.5856e+01, -1.6981e+01, 2.0512e-01, -5.0463e+00, 1.2586e-01},

      {4.3696e-01, 5.9373e-02, 1.5520e+00, 2.8488e+00, 4.4038e-02},  // Mean Hcal ShowerType3
      {5.4176e-01, 4.6653e+00, -1.1337e-01, 6.3126e-01, 4.5058e-01},
      {-1.3576e+00, 1.3274e-01, -2.0240e+00, 7.4848e-01, -1.3586e-01},
      {8.4917e+00, 1.0264e+01, -1.5633e-01, -5.1903e+00, 3.7948e-01},
      {2.4254e+01, -2.5119e+01, 2.1971e-01, -6.2339e+00, 1.1753e-01},

      {6.6849e-02, 1.1461e-01, -6.0693e-01, 3.5849e+00, 2.1510e-02},  // Sigma Ecal ShowerType1
      {4.4751e-01, 6.7414e-02, -1.7351e+00, 3.1800e+00, -1.4204e-02},
      {2.3277e-01, 9.8312e-02, -8.1359e-01, 3.3423e+00, 2.0069e-02},
      {6.0942e-01, 2.5435e-01, -5.7408e-01, 2.2588e+00, -2.7305e-03},
      {4.7227e-01, 6.9822e-01, -1.8614e-01, 5.0224e-01, 5.3256e-02},

      {2.7710e-01, 5.4557e-02, -9.1426e-01, 1.3652e+00, -2.1545e-02},  // Sigma Hcal ShowerType1
      {8.5475e-01, 9.0583e-02, 3.6811e+00, 2.1694e+00, -8.1316e-02},
      {9.9153e-01, 1.7054e-01, -2.3509e+00, 1.0596e+00, -7.7110e-02},
      {1.2402e+00, 6.9760e-01, -5.8527e-01, 1.3975e+00, -4.4764e-02},
      {1.1041e+00, 9.4012e-01, -4.0793e-01, 1.7098e+00, 6.9772e-03},

      {4.5757e-01, 2.2138e-01, -8.3016e-01, -6.5474e-01, -2.5193e-02},  // Sigma Hcal ShowerType2
      {4.6612e-01, 1.1899e-01, -8.7496e-01, 4.0150e+00, 2.3778e-02},
      {3.2115e+01, 3.1746e+01, -2.1227e-01, -1.1088e+01, 0.0000e+00},
      {1.4220e+00, 1.3834e+00, -3.3056e-01, 2.2772e-01, 3.5460e-02},
      {1.0770e+00, 2.3134e+00, -2.0993e-01, 7.6896e-01, 1.7177e-01},

      {1.0333e+00, 8.4267e-01, -2.2255e-01, -6.2499e+00, -2.0389e-02},  // Sigma Hcal ShowerType3
      {4.5615e-01, 1.4719e-01, -1.0516e+00, 3.2503e+00, 2.9697e-02},
      {7.1504e-01, 4.0454e-01, 6.8555e-04, -1.7412e+01, -5.0839e-02},
      {3.4149e+00, 3.8890e+00, -1.8545e-01, -3.4381e+00, 9.7396e-02},
      {3.7946e+00, 4.0784e+00, -1.7896e-01, -3.9979e+00, 5.1072e-02}};

  const double kminus_rho[8 * NPar][5] = {
      {-4.5951e-02, 3.6373e-01, 4.9181e-01, 6.9207e-01, 1.4429e-02},  // Ecal ShowerType1
      {1.6996e-01, 4.2946e-01, 6.2046e-01, -3.7396e-02, -1.4572e-02},
      {8.3137e-01, 2.9433e-02, 3.2221e+00, 1.6200e+00, -8.8843e-04},
      {-1.3483e-01, 9.3124e-02, -4.7555e+01, 1.9746e+00, 4.2961e-02},
      {1.4051e-01, 9.3445e-02, -3.3320e+00, 2.3752e+00, 3.7645e-02},
      {7.8256e-02, 1.2162e-01, -2.9812e+00, 2.1452e+00, 1.7049e-02},
      {-8.5479e-03, 9.6577e-02, -2.5794e+00, 2.1738e+00, 2.0528e-02},
      {1.9550e-01, 9.8667e-02, 2.3510e+00, 9.3278e-01, -1.0451e-02},
      {1.3706e-01, -1.1500e-01, 2.4339e+00, 2.3990e+00, 2.4792e-02},
      {7.8427e-01, 1.6946e-01, 6.7162e-01, 1.6939e+00, -6.8641e-04},

      {3.0252e-01, 2.0969e-01, 1.2233e+00, 2.5744e+00, -3.5561e-02},  // Hcal ShowerType1
      {5.0366e-01, 2.6734e-01, 9.1620e-01, 2.1240e+00, -7.5807e-02},
      {7.6809e-01, 5.2611e-02, -3.3231e+00, 1.0126e+00, 6.1410e-03},
      {-1.7444e-01, 4.5874e-01, 1.0344e-03, 2.5811e+00, 7.8505e-02},
      {2.5192e-01, 6.2636e-02, 1.3799e+00, 1.4204e+00, 6.1673e-02},
      {2.9799e-01, 1.5329e-01, 1.2920e+00, 1.3530e+00, -1.0006e-02},
      {-1.8200e-01, 4.7254e-01, 5.0000e-01, -3.6152e-01, -4.6662e-02},
      {1.6888e-01, -6.7814e-02, 6.5507e+00, 2.1306e+00, 4.1608e-02},
      {4.8725e-01, 1.3573e-01, 1.7021e+00, 8.4936e-01, -1.0376e-02},
      {8.3064e-01, 4.7621e-02, -1.4663e+00, 2.5843e+00, -4.1124e-02},

      {2.5820e-01, 2.1436e-01, 1.3064e+00, 2.1257e+00, -7.9308e-02},  // Hcal ShowerType2
      {4.9118e-01, 2.3914e-01, 2.0019e+00, 1.8543e+00, -8.3488e-02},
      {9.7119e-01, -2.4913e-01, 2.1162e+00, -2.2685e-01, -1.1841e-02},
      {-3.1042e-02, 8.5169e-02, 1.6265e+00, 1.3944e+00, 1.6725e-02},
      {2.8430e-01, 1.0624e-01, 1.2654e+00, 8.8490e-01, 2.1335e-02},
      {-1.7883e+00, 2.5865e+00, -2.6200e-01, 3.5317e+00, 5.8945e-01},
      {-3.4648e-01, 3.0452e-01, -8.5428e-01, 3.7070e+00, 1.4478e-01},
      {9.6655e-01, 1.1316e+00, -1.9048e-01, -3.5009e+00, 8.1655e-02},
      {5.3247e-01, 1.1047e-01, 2.1538e+00, 7.0882e-01, -2.0301e-02},
      {7.7207e-01, -8.5138e-02, 8.8082e-01, 2.0751e+00, -1.7783e-02},

      {4.4074e-02, 1.1797e-01, 9.8745e-01, 1.2754e+00, 0.0000e+00},  // Hcal ShowerType3
      {3.0829e-01, 3.7222e-01, 6.4042e-01, 7.1567e-01, -4.8572e-02},
      {7.2438e-01, 1.3993e-02, 2.8178e+01, 1.1987e+00, -1.8113e-03},
      {3.3772e-02, 5.2144e-02, -2.1020e+01, 2.2369e+00, 0.0000e+00},
      {-8.7740e+01, 8.8231e+01, 5.0000e-01, -6.4849e+00, -2.5865e-02},
      {1.6919e-01, 1.5093e-01, -1.6970e+00, 3.1185e+00, 0.0000e+00},
      {-2.7414e-01, 4.2637e-01, -4.5968e-01, 4.0680e+00, 1.1120e-01},
      {3.5779e-01, 1.3748e+00, 1.5201e-01, 4.5045e-01, -1.4935e-01},
      {-4.4323e+01, 4.5858e+01, 1.4803e-01, -1.5183e+01, -1.7599e-01},
      {4.4308e-01, -1.9361e+00, 1.7101e-01, 1.2267e+00, 2.6132e-01}};

  // kaon-minus lateral parameters

  const double kminus_rpar[4 * Nrpar][5] = {
      {1.5255e+00, 7.1460e-01, -9.6109e-01, 2.1485e+00, 0.0000e+00},  // Ecal for showerType = 1
      {8.1702e+00, 1.1618e+00, -3.2615e+00, 1.9965e+00, -9.8151e-01},
      {8.0394e-01, 1.1379e-01, -5.7818e+01, 1.8315e+00, -3.1074e-02},
      {-3.4073e-01, 1.5998e-01, 5.6939e+01, 1.8305e+00, 1.5185e-02},

      {1.8546e+01, 3.9927e+00, -6.4386e-01, 2.5076e+00, -1.3842e+00},  // Hcal for showerType = 1
      {-2.9410e-01, 6.9922e-01, -7.1734e+01, 1.3073e+00, 8.2441e-02},
      {4.8144e-01, 6.6668e-02, 1.4955e+00, 1.5032e+00, -2.8181e-02},
      {-5.4174e-03, 1.7362e-02, 1.2282e+01, 2.4035e-01, 6.1898e-04},

      {2.0741e+01, 3.9654e+01, -1.8029e-01, -4.2558e-01, 2.5375e+00},  // Hcal for showerType = 2
      {-1.1271e+00, 5.6222e-01, -5.1185e+01, 9.9351e-01, 1.9139e-01},
      {4.8514e-01, 2.6816e-02, -3.7361e+01, 2.7567e+00, -8.2626e-03},
      {-9.0682e-02, 7.2902e-02, 2.3265e+00, -6.5459e+00, 8.8080e-03},

      {-2.9248e+00, 5.6213e+00, -3.2437e-01, 3.0222e+00, 1.5066e+00},  // Hcal for showerType = 3
      {7.4285e+00, 3.1786e+00, -9.1437e-01, 1.1013e+00, -7.0831e-01},
      {9.6057e-01, 3.6162e-01, -3.6673e-01, -8.9691e-01, -3.2462e-02},
      {-2.9505e+01, 2.9362e+01, 1.0665e+00, -2.3946e+00, 1.6595e-02}};

  // Kaon-Plus paramters

  const double kplus_emscale[2][5] = {{7.8775e-01, 2.7022e-01, -1.4291e+00, 1.3980e+00, -1.2795e-02},
                                      {2.2224e-01, 1.7480e-02, 1.1984e+01, 1.8265e+00, -1.6718e-02}};

  const double kplus_hadscale[8][5] = {
      // based on LogNormal of Had Energy
      {-1.1629e+02, 1.1506e+02, 4.8753e-01, -4.4057e+00, 1.0940e+00},
      {2.2018e+00, 1.1327e+00, -6.4208e-01, 1.2240e-01, -7.2109e-02},
      {1.5217e+00, 1.9949e+00, -4.0047e-01, 1.0897e+00, 1.4482e-01},
      {-1.2902e+01, 1.5504e+01, 1.1257e-01, -8.4732e+00, -3.0451e-01},
      // LogNormal of Had Energy for mip-like response
      {3.4887e-01, 1.9201e-01, 8.5205e-01, 3.8457e+00, 8.5097e-01},
      {5.1252e-01, 7.7619e-01, -1.3001e-01, -4.1296e+00, 3.8339e-02},
      // correction of response if ssp is between Ecal and Hcal
      {1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00},
      {8.2313e-01, 1.3027e+00, 2.5263e-01, 2.8203e+00, -2.8548e-01}};

  const double kplus_correl_hadem[5] = {-6.9227e-01, 1.8762e-01, -7.1395e-01, 7.1711e-01, -7.3016e-03};

  // kaon-plus longitudinal parameters

  const double kplus_par[8 * NPar][5] = {
      {3.8772e-01, 7.5923e-02, 2.3857e+00, 2.8135e+00, 7.0675e-02},  // Mean Ecal ShowerType1
      {1.4815e+00, 9.8518e-01, 3.4748e-01, 3.3609e+00, -1.2142e-01},
      {-7.4688e-01, 7.4632e-01, 3.7387e-01, 3.1767e+00, -1.2898e-01},
      {6.4780e-01, 3.3263e-01, 1.6380e+00, 2.8163e+00, 7.5372e-02},
      {4.8648e+00, 4.4851e+00, 2.5096e-01, 3.1485e+00, -9.4875e-01},

      {1.6983e+00, 1.3064e+00, 2.8039e-01, 4.3610e+00, -2.8089e-01},  // Mean Hcal ShowerType1
      {1.4443e+00, -5.0550e-01, 1.0680e+00, 5.6710e-01, -1.5739e-02},
      {-5.3822e-01, 1.0584e+00, -6.4612e-01, 3.5104e-01, -1.1463e-01},
      {2.3248e+00, -2.2851e+00, 3.6709e-01, 1.3197e-01, 1.2563e-01},
      {2.1837e+00, -2.2043e+00, 4.6015e-01, 9.2583e-01, -7.4904e-03},

      {3.3203e+00, 5.0310e+00, 1.6001e-01, 3.9441e+00, -7.1274e-01},  // Mean Hcal ShowerType2
      {1.1513e+00, -1.8085e-01, 3.3413e+00, 7.0993e-01, -6.3645e-02},
      {1.3898e+01, 1.5638e+01, -4.1890e-01, -3.9674e+00, -1.1401e-01},
      {2.0012e+01, 1.9866e+01, -2.9212e-01, -5.0006e+00, 1.1974e-01},
      {2.4295e+00, -2.2834e+00, 4.8288e-01, 2.9977e-01, -3.0573e-02},

      {3.8408e-01, 6.1344e-02, 1.0667e+00, 3.0242e+00, 5.0278e-02},  // Mean Hcal ShowerType3
      {1.8964e+00, 5.1610e+00, -8.8169e-02, -1.9272e+00, 3.1513e-01},
      {1.4200e+01, 1.6127e+01, -2.8959e-01, -6.1679e+00, -7.4159e-02},
      {2.7497e+01, 2.7062e+01, -3.3540e-01, -5.5799e+00, 1.0200e-01},
      {2.3118e+01, -2.3752e+01, 2.6134e-01, -5.2177e+00, 9.7844e-02},

      {-3.5273e-02, 2.6492e-01, -2.7154e-01, 4.1599e+00, 3.7807e-02},  // Sigma Ecal ShowerType1
      {5.4462e-01, 4.9507e-02, -1.8513e+01, 2.9250e+00, -3.6791e-02},
      {3.7860e-01, 1.8549e-02, -2.9438e+00, 2.4898e+00, -2.2790e-02},
      {6.6597e-01, 6.5499e-02, -2.1924e+01, 2.9525e+00, -4.7980e-02},
      {1.0753e+00, 9.2583e-01, -1.4431e-01, -5.0882e+00, 2.9234e-03},

      {2.6008e-01, 6.1267e-02, -1.7484e+00, 1.6088e+00, -1.7668e-02},  // Sigma Hcal ShowerType1
      {8.6750e-01, 8.7283e-02, 2.5881e+00, 2.5618e+00, -8.5318e-02},
      {1.0368e+00, 2.5737e-01, -1.3413e+00, 8.7206e-01, -7.0973e-02},
      {1.2517e+00, 6.4973e-01, -7.3033e-01, 1.4095e+00, -5.3997e-02},
      {1.9193e+00, 1.9482e+00, -3.0639e-01, 6.6016e-02, 2.9074e-02},

      {2.7361e-01, 5.8267e-02, -1.5294e+00, 1.0185e+00, -2.2741e-02},  // Sigma Hcal ShowerType2
      {5.3653e-01, 1.0764e-01, -6.6807e-01, 4.7074e+00, -3.4054e-03},
      {5.8520e+01, 5.8162e+01, -2.3280e-01, -1.1136e+01, 0.0000e+00},
      {1.2192e+00, 6.9207e-01, -6.8090e-01, 8.9059e-01, -3.8146e-02},
      {3.2972e+00, 2.9315e+00, -3.4433e-01, -1.6434e+00, -2.8242e-02},

      {1.4349e+00, 1.4443e+00, -1.2105e-01, -9.0252e+00, 2.3150e-03},  // Sigma Hcal ShowerType3
      {4.9465e-01, 1.8568e-01, -9.3747e-01, 3.2665e+00, 2.3940e-02},
      {2.3407e+00, -3.8601e+00, 2.9374e-02, -1.4549e+01, 1.4344e-02},
      {4.8451e+00, 4.6687e+00, -2.7170e-01, -3.7684e+00, 1.8616e-02},
      {5.0127e+00, 5.0025e+00, -2.3093e-01, -4.0513e+00, 2.3821e-02}};

  const double kplus_rho[8 * NPar][5] = {
      {3.8772e-01, 7.5923e-02, 2.3857e+00, 2.8135e+00, 7.0675e-02},  // Ecal ShowerType1
      {1.4815e+00, 9.8518e-01, 3.4748e-01, 3.3609e+00, -1.2142e-01},
      {-7.4688e-01, 7.4632e-01, 3.7387e-01, 3.1767e+00, -1.2898e-01},
      {6.4780e-01, 3.3263e-01, 1.6380e+00, 2.8163e+00, 7.5372e-02},
      {4.8648e+00, 4.4851e+00, 2.5096e-01, 3.1485e+00, -9.4875e-01},
      {-3.5273e-02, 2.6492e-01, -2.7154e-01, 4.1599e+00, 3.7807e-02},
      {5.4462e-01, 4.9507e-02, -1.8513e+01, 2.9250e+00, -3.6791e-02},
      {3.7860e-01, 1.8549e-02, -2.9438e+00, 2.4898e+00, -2.2790e-02},
      {6.6597e-01, 6.5499e-02, -2.1924e+01, 2.9525e+00, -4.7980e-02},
      {1.0753e+00, 9.2583e-01, -1.4431e-01, -5.0882e+00, 2.9234e-03},

      {1.6983e+00, 1.3064e+00, 2.8039e-01, 4.3610e+00, -2.8089e-01},  // Hcal ShowerType1
      {1.4443e+00, -5.0550e-01, 1.0680e+00, 5.6710e-01, -1.5739e-02},
      {-5.3822e-01, 1.0584e+00, -6.4612e-01, 3.5104e-01, -1.1463e-01},
      {2.3248e+00, -2.2851e+00, 3.6709e-01, 1.3197e-01, 1.2563e-01},
      {2.1837e+00, -2.2043e+00, 4.6015e-01, 9.2583e-01, -7.4904e-03},
      {2.6008e-01, 6.1267e-02, -1.7484e+00, 1.6088e+00, -1.7668e-02},
      {8.6750e-01, 8.7283e-02, 2.5881e+00, 2.5618e+00, -8.5318e-02},
      {1.0368e+00, 2.5737e-01, -1.3413e+00, 8.7206e-01, -7.0973e-02},
      {1.2517e+00, 6.4973e-01, -7.3033e-01, 1.4095e+00, -5.3997e-02},
      {1.9193e+00, 1.9482e+00, -3.0639e-01, 6.6016e-02, 2.9074e-02},

      {2.4423e-01, 2.2650e-01, 1.3877e+00, 1.8939e+00, -7.5389e-02},  // Hcal ShowerType2
      {4.6760e-01, 2.6361e-01, 1.4374e+00, 1.6434e+00, -9.9844e-02},
      {7.4796e-01, -4.4661e-02, 1.7292e+01, 7.0879e-01, 2.4763e-03},
      {-6.4806e+00, 6.4231e+00, 1.0400e+00, -2.2650e+00, 4.7778e-02},
      {3.4440e-01, 1.1933e-01, 2.1685e+00, 1.3074e+00, 1.5225e-02},
      {-1.1990e+00, 1.8273e+00, -2.2560e-01, 4.3776e+00, 3.3393e-01},
      {-1.3851e+00, 1.7655e+00, -2.6978e-01, 3.6071e+00, 4.2187e-01},
      {2.8248e-01, 4.8579e-02, 1.0095e+00, 3.5771e+00, 1.1249e-02},
      {-4.8075e+00, 5.5604e+00, 5.3060e-01, -3.1000e+00, -4.5298e-02},
      {7.2139e-01, -1.2988e-01, 1.3460e+00, 1.7526e+00, 8.9762e-03},

      {4.7434e-02, 3.1872e+00, 1.6477e-01, 3.0236e-01, -4.2486e-01},  // Hcal ShowerType3
      {9.4033e-02, 5.7792e-01, 6.8058e-01, 2.1194e-01, -8.2046e-02},
      {7.4700e-01, 2.0561e-02, 4.2268e+01, 1.6399e+00, -1.1949e-02},
      {5.2255e-02, 4.6034e-02, -1.7110e+00, 1.8247e+00, 0.0000e+00},
      {-1.0684e+02, 1.0739e+02, 5.0000e-01, -6.3340e+00, -2.6291e-02},
      {2.0912e-01, 1.3043e-01, -4.9417e+00, 2.9813e+00, 0.0000e+00},
      {-1.8831e+00, 3.4253e+00, -1.8757e-01, 3.5963e+00, 5.6274e-01},
      {1.4591e-03, 1.6858e+00, 1.2556e-01, -7.6405e-01, -1.1943e-01},
      {-6.2423e+00, 7.1383e+00, 3.8814e-01, -4.3872e+00, -9.4735e-02},
      {7.2200e-01, -1.1568e-01, 1.3604e+00, 1.3338e+00, 1.6450e-02}};

  // kaon-plus lateral parameters

  const double kplus_rpar[4 * Nrpar][5] = {
      {2.0917e+00, 1.2528e+00, -6.9494e-01, 1.7600e+00, 0.0000e+00},  // Ecal for showerType = 1
      {7.8866e+00, 6.3777e-01, -4.3657e+00, 2.0492e+00, -1.0208e+00},
      {7.3150e-01, 1.6104e-01, -8.4279e+00, 2.0670e+00, -3.3696e-03},
      {-2.4312e-01, 1.8371e-01, 1.5265e+01, 2.0458e+00, -1.6358e-02},

      {2.0008e+01, 2.1622e+00, -3.4949e+00, 1.9954e+00, -2.0569e+00},  // Hcal for showerType = 1
      {-4.1668e-02, 5.4366e-01, -1.3780e+01, 7.0957e-01, -1.1259e-02},
      {5.0855e-01, 1.0979e-01, 1.2233e+00, 1.4132e+00, -4.5390e-02},
      {-5.9436e-02, -1.0396e-01, 3.7419e-01, 1.9661e+00, 3.6075e-02},

      {1.2970e+01, 2.5925e+00, -1.5404e+00, 1.6391e+00, -1.3978e+00},  // Hcal for showerType = 2
      {1.9183e-01, 1.3373e+00, -1.1090e+00, -7.6909e+00, 6.3587e-02},
      {4.3892e-01, 7.0247e-02, -1.8419e+01, 2.9160e+00, 1.2592e-02},
      {6.2448e-02, 6.3318e-02, 6.1143e-01, 2.4479e+00, -2.4459e-02},

      {2.0954e+00, 4.0916e-01, -6.9750e+01, 1.9870e+00, -5.3104e-02},  // Hcal for showerType = 3
      {6.4162e+00, 1.8264e+00, -9.6305e-01, 1.5012e+00, -7.5696e-01},
      {4.3430e+00, 1.3126e+01, -2.5342e-02, -1.0939e+01, 2.1736e-01},
      {-5.0784e+01, 5.0870e+01, 2.0000e-01, -1.3669e+01, -1.2081e-02}};

  // HcalOuter parameters
  const double ho_nonzero[5] = {4.79943e-01, 4.61158e-01, 7.09011e-01, 4.86440e+00, 0.0000e+00};
}  // namespace Gflash

#endif