CMS 3D CMS Logo

/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/SimGeneral/GFlash/interface/GflashNameSpace.h

Go to the documentation of this file.
00001 #ifndef GflashNameSpace_H
00002 #define GflashNameSpace_H 
00003 
00004 #include "SimGeneral/GFlash/interface/Gflash3Vector.h"
00005 
00006 namespace Gflash {
00007 
00008   //unit convention: energy in GeV, length in cm, time in nanosecond
00009 
00010   enum CalorimeterNumber {
00011     kNULL = -1,         // ECAL Barrel - ESPM
00012     kESPM,              // ECAL Barrel - ESPM
00013     kHB,                // HCAL Barrel - HB
00014     kENCA,              // ECAL Endcap - ENCA
00015     kHE,                // HCAL Endcap - HE
00016     kHO,                // HCAL Outer  - HO
00017     kHF,                // HCAL Forward - HF
00018     kNumberCalorimeter
00019   };
00020   
00021   CalorimeterNumber getCalorimeterNumber(const Gflash3Vector position);
00022   double rhoBackEB(const Gflash3Vector position);
00023   double zBackEE(const Gflash3Vector position);
00024   int findShowerType(const Gflash3Vector position);
00025 
00026   //                                            EB     HB     EE     HE      HO,     HF
00027   const double Zmin[kNumberCalorimeter]   = {0.000, 0.000, 317.2, 391.95, 0.000, 1110.0}; // in cm
00028   const double Zmax[kNumberCalorimeter]   = {317.0, 433.2, 399.1, 554.10, 661.0, 1275.0};
00029   const double Rmin[kNumberCalorimeter]   = {123.8, 177.5,  31.6,  31.6,  382.0,  12.5};
00030   const double Rmax[kNumberCalorimeter]   = {177.5, 287.7, 171.1, 263.9,  407.0, 140.0};
00031 
00032   const double EtaMin[kNumberCalorimeter] = {0.000, 0.000, 1.479,  1.479, 0.000, 3.000};
00033   const double EtaMax[kNumberCalorimeter] = {1.479, 1.479, 3.000,  3.000, 1.262, 5.000};
00034 
00035   const double RFrontCrystalEB = 129.0; // rho in cm
00036   const double ZFrontCrystalEE = 320.9; // z in cm
00037   const double ROffCrystalEB = 5.2; // rho in cm
00038   const double ZOffCrystalEE = 3.7; // z in cm
00039   const double LengthCrystalEB  = 23.0;
00040   const double LengthCrystalEE  = 22.0;
00041 
00042   //constants needed for GflashHadronShowerModel
00043   const double energyCutOff                 = 1.0;
00044 
00045   //cut value for quasi-elastic like interactions on the secondary energy, 
00046   //ratio = (Leading Energy of Secondaries)/(Total Energy)
00047   const double QuasiElasticLike = 0.95;
00048 
00049   //minimum distance to the back of parameterized envelopes in [cm]
00050   const double MinDistanceToOut = 10.0;
00051     
00052   //constants needed for GflashHadronShowerProfile
00053 
00054   const double rMoliere[kNumberCalorimeter]  = { 2.19, 2.19, 2.19, 2.19, 2.19, 1.72}; // in cm
00055   const double radLength[kNumberCalorimeter] = { 0.89, 1.49, 0.89, 1.49, 1.49, 1.76}; // in cm
00056   const double intLength[kNumberCalorimeter] = { 22.4,16.42, 22.4, 16.42, 16.42, 16.77}; // in cm
00057   const double Z[kNumberCalorimeter]         = { 68.36, 68.36, 68.36, 68.36, 68.36, 55.845}; // mass of molicule
00058   const double criticalEnergy                = 0.0086155;
00059 
00060   // The step size of showino along the helix trajectory in cm unit
00061   const double divisionStep = 1.0;
00062 
00063   //maximum shower depth for the hadronic lateral in interaction length unit
00064   const double maxShowerDepthforR50 = 6.0;
00065 
00066   //maximum size of the lateral arm for hadronic showers = 12 * interaction length
00067   const double maxLateralArmforR50 = 197.0;
00068 
00069   //minimum particle energy to trigger HO parameterization in [GeV]
00070   const double MinEnergyCutOffForHO = 2.5;
00071 
00072   //additional energy scale for the Hcal sensitive detector
00073   const double scaleSensitiveHB = 0.200;
00074   const double scaleSensitiveHE = 0.200;
00075 
00076   // properties for different sub-detectors (calorimeters)
00077 
00078   const double SAMHAD[3][kNumberCalorimeter] = {{0.0,0.89,0.0,0.89,0.89,0.0},
00079                                                   {0.0,0.00,0.0,0.00,0.00,0.0},
00080                                                   {0.0,0.00,0.0,0.00,0.00,0.0}};
00081   const double RLTHAD[kNumberCalorimeter] = {32.7,23.7,32.7,23.7,23.7,23.7};
00082 
00083 
00084   //parameters for logitudinal parameterizations and energy responses
00085 
00086   const int NPar  = 5;
00087 
00088   const double fdep[2][4] = {
00089     {  8.5073e-01,  2.9650e-02,  8.1327e-01,  2.7547e+00},
00090     {  9.7806e-02,  4.6989e-02,  0.0000e+00,  0.0000e+00}
00091   };
00092 
00093   const double emscale[4][4] = {
00094     {  5.4463e-01, -4.1210e-02,  1.8231e+00,  4.1472e+00 },
00095     { -3.4608e-01, -1.7809e-01,  1.1329e+00,  2.0367e+00 }, //not used 
00096     {  1.7327e-01,  3.1683e-02,  1.0484e+00,  2.3877e+00 },
00097     {  3.5000e-02,  4.5000e-02,  1.5000e+00,  2.5000e+00 }
00098   };
00099 
00100   const double hadscale[7][4] = {
00101     //based on LogNormal of Had Energy + (0.05*Einc)% swift 
00102     { -3.1358e+01,  5.3531e+01,  3.8817e-02, -1.5548e+01 },
00103     {  8.2970e-01,  2.6359e-01, -1.2500e+00,  3.7566e+00 },
00104     {  6.1175e-01,  1.3195e-01, -2.0953e+00,  2.9990e+00 },
00105     { -3.1528e-01,  8.0146e-02, -6.1015e+00,  1.8951e+00 },
00106     //LogNormal of Had Energy for mip-like response based on 2006 H2 test beam data
00107     {  4.4527e+00,  5.7913e+00,  1.6584e-01,  4.6872e+00 },
00108     {  2.4081e-01,  1.6057e-01, -6.1439e-01,  1.9368e+00 },
00109     //correction of response if ssp is between Ecal and Hcal
00110     {  1.8736e-01,  1.7391e-01,  1.0991e+00,  3.3557e+00 }
00111   };
00112 
00113   const double correl_hadem[4] = { -7.8255e-01,  1.7976e-01, -8.8001e-01,  2.3474e+00 };
00114 
00115   //Proton paramters
00116   const double pro_emscale[4][4] = {
00117     {  5.0463e-01, -8.1210e-02,  1.8231e+00,  2.7472e+00 },
00118     { -3.4608e-01, -1.7809e-01,  1.1329e+00,  2.0367e+00 }, //not used 
00119     {  1.7327e-01,  3.1683e-02,  1.0484e+00,  2.3877e+00 },
00120     {  3.5000e-02,  4.5000e-02,  1.5000e+00,  2.5000e+00 }
00121   };
00122 
00123   const double pro_hadscale[7][4] = {
00124     //based on LogNormal of Had Energy 
00125     { -5.0239e+00,  1.2308e+01,  1.9628e-01,  0.0000e+00 },
00126     {  8.2970e-01,  2.3590e-01, -8.0000e+00,  4.0000e+00 },
00127     {  9.9975e-01,  6.0775e-01, -9.2005e-01,  2.8056e+00 },
00128     { -6.2579e-01,  2.9609e-01,  2.0549e+00,  3.7047e+00 },
00129 
00130     {  3.7385e+00,  3.3977e+00,  3.0798e-01,  3.8725e+00 },
00131     {  2.4081e-01,  1.6057e-01, -6.1439e-01,  1.9368e+00 },
00132     //correction of response if ssp is between Ecal and Hcal
00133     {  1.8736e-01,  1.7391e-01,  1.0991e+00,  3.3557e+00 }
00134   };
00135 
00136   const double emcorr_pid[10][4] = {
00137     {  9.3665e-02,  1.0012e-01, -1.1824e+00,  1.1446e+00 }, //pi eb
00138     {  5.6892e-02,  5.6048e-02, -9.0275e-01,  1.5974e+00 }, //pi ee
00139     {  1.0012e-01,  1.1394e-01, -2.2312e+00,  1.3881e+00 }, //K- eb
00140     {  6.1321e-02,  7.0704e-02, -1.3792e+00,  1.6116e+00 }, //K- ee
00141     {  3.9237e-02,  9.8766e-02, -1.6805e+00,  1.5693e+00 }, //K+ eb
00142     {  4.3418e-03,  5.4551e-02, -2.0302e+00,  2.0768e+00 }, //K+ ee
00143     {  2.4292e-01,  2.8905e-01, -2.3882e+00,  6.9453e-01 }, //pbar eb
00144     {  2.3073e-01,  2.9771e-01, -5.8904e-01, -4.1660e-01 }, //pbar ee
00145     {  5.1063e-02,  7.7481e-02, -3.4550e+00,  1.2754e+00 }, //p eb
00146     {  1.3621e-02,  3.9391e-02, -2.9119e+00,  1.6276e+00 }  //p ee
00147   };
00148 
00149   const double hadcorr_pid[10][4] = {
00150     {  2.0727e-02,  8.4981e-02, -4.4552e+00,  1.2077e+00 }, //pi eb
00151     {  3.9438e-03,  5.3269e-02, -1.3471e+01,  1.9124e+00 }, //pi ee
00152     {  5.5982e-02,  1.0265e-01, -3.4828e+00,  1.3052e+00 }, //K- eb
00153     {  4.9414e-02,  9.4409e-02, -1.7006e+00,  1.6530e+00 }, //K- ee
00154     {  9.7977e-02,  1.1940e-01, -4.8080e+00,  1.1847e+00 }, //K+ eb
00155     {  1.1574e-01,  1.1238e-01, -2.6366e+00,  1.5198e+00 }, //K+ ee
00156     {  2.2534e-02,  1.1664e-01, -4.6129e+00,  9.6495e-01 }, //pbar eb
00157     { -8.6348e-02,  4.5002e-02, -3.4570e+01,  1.1149e+00 }, //pbar ee
00158     { -5.3266e-02,  1.3402e-01, -8.8320e-01,  3.5483e+00 }, //p eb
00159     { -8.0026e-03,  6.9235e-02, -2.8974e+00,  3.4128e+00 }  //p ee
00160   };
00161 
00162   const double mipcorr_pid[10][4] = {
00163     {  9.5677e-02,  1.3784e-01, -4.2091e+00,  1.2392e+00 }, //pi eb
00164     {  6.2573e-02,  1.5033e-01, -1.5482e+00,  1.5297e+00 }, //pi eb for pi ee
00165     {  1.4989e-01,  1.9360e-01, -3.1175e+00,  1.4048e+00 }, //K- eb
00166     {  1.0899e-01,  1.8504e-01, -2.4958e+00,  1.6699e+00 }, //K- eb for K- ee
00167     {  1.4712e-01,  1.9871e-01, -3.4212e+00,  1.3119e+00 }, //K+ eb
00168     {  1.0670e-01,  1.9217e-01, -3.0280e+00,  1.4374e+00 }, //K+ eb for K+ ee
00169     {  1.2268e-01,  2.9728e-01, -2.7384e+01,  6.9221e-01 }, //pbar eb
00170     {  1.6795e-02,  2.0576e-01, -2.5716e+01,  7.1464e-01 }, //pbar ee
00171     {  7.0818e-02,  2.0474e-01, -7.2400e-01,  1.6688e+00 }, //p eb
00172     { -5.3189e-03,  2.0254e-01, -7.1329e-01,  2.2570e+00 }  //p eb for p ee
00173   };
00174 
00175   //longitudinal parameters
00176 
00177   const double par[8*NPar][4] = {
00178     {  6.7438e-01,  1.7413e-01,  8.1317e-01,  2.6713e+00 }, //Mean Ecal ShowerType1
00179     {  9.5916e-01,  1.0649e+00,  1.8075e-01,  3.0428e+00 },
00180     { -1.2427e+00,  2.9128e-01,  6.2858e-01,  2.6237e+00 },
00181     {  7.9194e-01, -4.9592e-01, -1.1214e+00,  2.6263e+00 },
00182     {  1.8503e+00,  1.7336e-01,  1.6952e+00,  2.9239e+00 },
00183 
00184     {  6.8778e-01,  1.3956e-01, -1.7557e+00,  1.5919e+00 }, //Mean Hcal ShowerType1
00185     {  7.0000e-01, -5.0000e-02,  4.2028e+00,  1.4877e+00 },
00186     { -1.5000e+00, -9.5593e-01,  2.3739e-01,  5.8041e-01 },
00187     {  3.0602e+00, -2.3879e+00,  4.6035e-01, -4.5823e-01 },
00188     {  2.0000e+00,  2.0000e+00, -4.0535e-01,  8.4006e-01 },
00189 
00190     {  6.0622e-01,  3.4491e-02,  1.9601e+00,  4.0597e+00 }, //Mean Hcal ShowerType2
00191     {  6.1856e-01, -5.3972e-02,  4.2128e+00,  3.0421e+00 },
00192     { -1.9157e+00,  5.9675e-01, -3.1110e-01,  2.3064e+00 },
00193     {  1.1704e+00,  4.1968e-01, -7.8091e-01,  1.6688e+00 },
00194     {  2.0000e+00, -2.0000e+00,  4.3481e-01,  7.7922e-01 },
00195 
00196     {  6.3159e-01,  1.4557e-01,  4.7778e-01,  2.5534e+00 }, //Mean Hcal ShowerType3
00197     {  6.4682e-01, -7.3109e-02,  3.4524e+01,  2.9954e+00 },
00198     { -1.9768e+00,  3.1293e-01, -5.6383e-01,  3.0850e+00 },
00199     {  9.6181e-01,  1.1237e-01, -1.7285e+00,  2.3366e+00 },
00200     {  1.5000e+00, -1.5000e+00,  4.2321e-01,  8.9595e-01 },
00201 
00202     { -4.4895e+00,  5.2126e+00, -2.3116e-02,  6.4173e+01 }, //Sigma Ecal ShowerType1
00203     {  3.3856e-01,  7.4435e-02, -1.4309e+00,  2.6740e+00 },
00204     {  2.6759e-01,  2.5205e-02, -2.3270e+00,  3.3109e+00 },
00205     {  6.7557e-01,  3.2703e-01, -7.6511e-01,  1.7957e+00 },
00206     {  4.3271e-01,  2.0000e-01, -4.7512e-01,  1.0227e+00 },
00207 
00208     {  3.0677e-01,  2.2520e-01, -6.2077e-01,  1.1081e+00 }, //Sigma Hcal ShowerType1
00209     {  4.0648e-01,  7.7278e-03,  2.4117e+02,  2.2477e+00 },
00210     {  4.4076e-01,  2.4788e-01, -2.5523e-01,  3.3124e+00 },
00211     {  1.1130e+00,  9.1704e-01, -3.6457e-01,  1.1608e+00 },
00212     {  1.5000e+00, -1.3000e+00,  3.2784e-01,  3.5213e-01 },
00213 
00214     {  2.5000e-01,  1.6444e-01, -3.9152e-01,  8.4963e-01 }, //Sigma Hcal ShowerType2
00215     {  4.9724e-01, -3.6006e-02,  1.1210e+00,  3.9498e+00 },
00216     {  4.2496e-01, -1.6087e-01,  4.0393e-01,  3.9751e+00 },
00217     {  8.5041e-01,  5.1464e-01, -6.6302e-01,  1.9777e+00 },
00218     {  1.3960e+00, -1.2054e+00,  3.3610e-01,  7.0376e-01 },
00219 
00220     {  2.5000e-01,  2.1781e-01, -2.2974e-01,  6.2906e-01 }, //Sigma Hcal ShowerType3
00221     {  5.2411e-01, -1.0607e-01,  7.4782e-01,  2.7602e+00 },
00222     {  1.0462e+00, -1.9208e+01,  2.7239e-03, -6.7807e+00 },
00223     {  6.6619e-01,  3.1562e-01, -7.8957e-01,  2.2168e+00 },
00224     {  8.2528e-01, -6.2222e-01,  3.9534e-01,  1.8290e+00 }
00225   };
00226 
00227   const double rho[8*NPar][4] = {
00228     {  1.4190e-01,  2.6308e-01,  1.2979e+00,  1.9024e+00 }, //Ecal ShowerType1
00229     {  4.7805e-01,  1.3546e-01,  1.5815e+00,  1.6697e+00 },
00230     {  8.2088e-01,  2.1833e-02,  2.2896e+00,  1.5505e+00 },
00231     { -5.2167e-02,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00232     {  2.7072e-01, -6.9839e-02,  1.0111e+02,  2.2077e+00 },
00233     {  1.3216e-01, -1.0180e-01,  9.8382e+00,  2.2553e+00 },
00234     {  5.3132e-02, -6.8272e-02,  1.3674e+01,  2.3076e+00 },
00235     {  2.7417e-01,  0.0000e+00,  4.0000e+00,  1.0000e+00 },
00236     {  2.1588e-01, -1.0111e-01,  1.1345e+01,  2.3139e+00 },
00237     {  8.0280e-01,  1.4780e-01,  8.5541e-01,  2.2995e+00 },
00238 
00239     {  5.9673e-02,  2.6255e-01,  7.3541e-01,  2.7660e+00 }, //Hcal ShowerType1
00240     {  2.7126e-01,  1.5810e-01,  1.3596e+00,  2.4690e+00 },
00241     {  7.5047e-01,  6.8428e-02,  8.7936e-01,  5.0565e+00 },
00242     {  8.9725e-02,  8.0812e-02,  2.2949e+00,  4.0534e+00 },
00243     {  5.0000e-01,  6.6307e+01,  1.6838e-03,  3.2369e+00 },
00244     {  2.7641e-01,  1.5000e-01,  1.0194e+00,  1.2387e+00 },
00245     {  1.2386e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00246     {  3.8881e-01,  1.8000e-01,  8.3829e-01,  5.1864e+00 },
00247     {  4.6535e-01,  1.5000e-01,  9.3439e-01,  9.1905e-01 },
00248     {  8.0000e-01,  2.0000e-01, -5.7854e-01,  2.0167e+00 },
00249 
00250     { -2.2610e-02,  8.7079e-02,  1.3255e+01,  2.2260e+00 }, //Hcal ShowerType2
00251     {  2.7147e-01,  1.2332e-01,  3.2099e+00,  1.9717e+00 },
00252     {  6.6040e-01, -1.1842e-02,  6.8650e+00,  2.0085e+00 },
00253     {  1.0205e-01,  1.0000e-01,  5.4438e-01,  3.6551e+00 },
00254     {  3.0000e-01,  2.0000e-01,  8.2230e-01,  1.3549e+00 },
00255     {  3.4682e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00256     {  1.7160e-01,  3.9781e-02,  6.7383e+00,  1.9256e+00 },
00257     {  2.0000e-01,  1.3237e-01,  2.5413e-01,  1.2529e+00 },
00258     {  6.1090e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00259     {  7.2769e-01, -1.2600e-01,  8.1386e-01,  2.4818e+00 },
00260 
00261     {  1.0345e-01,  4.8906e-02,  6.6775e+01,  2.0782e+00 }, //Hcal ShowerType3
00262     {  4.2870e-01,  1.0980e-01,  4.7327e+00,  1.7190e+00 },
00263     {  7.0416e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00264     { -1.0516e-02, -1.0000e-01,  1.0603e+00,  4.5948e+00 },
00265     {  3.9288e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00266     {  4.0569e-01, -3.4859e+01,  2.7762e-03,  9.9758e-01 },
00267     {  2.1099e-01, -1.0207e-01,  1.6655e+00,  4.2542e+00 },
00268     {  2.0000e-01,  2.7118e-01,  6.1110e-01,  1.0454e+00 },
00269     {  5.5325e-01, -7.0746e-02,  3.5839e+00,  4.3084e+00 },
00270     {  7.3839e-01, -1.1350e-01,  6.3423e-01,  2.1142e+00 }
00271   };
00272 
00273   //lateral parameters
00274   const int Nrpar = 4;
00275   const double rpar[4*Nrpar][2] = {
00276     {  1.6065e+00, -1.9118e-01}, //Ecal for showerType = 1
00277     {  8.3070e+00, -1.2512e+00},
00278     {  6.2098e-01, -2.9482e-02},
00279     { -1.6002e-01,  1.9410e-02},
00280     {  2.1779e+01, -2.6719e+00}, //Hcal for showerType = 1
00281     { -1.1468e-01, -1.2217e-01},
00282     {  4.4093e-01, -1.5766e-02},
00283     { -2.3089e-02,  9.3034e-03},
00284     {  1.3713e+01, -1.9910e+00}, //Hcal for showerType = 2
00285     { -1.4097e+00,  7.2315e-02},
00286     {  4.4531e-01, -2.8746e-03},
00287     {  3.1108e-02, -4.1437e-04},
00288     {  1.9392e+00,  3.4218e-02}, //Hcal for showerType = 3
00289     {  7.8122e+00, -1.5978e+00},
00290     {  6.6561e-01, -3.9986e-02},
00291     { -2.5633e-01,  5.7796e-02}
00292   };
00293 
00294   //Anti-proton paramters
00295 
00296   const double pbar_emscale[2][4] = {
00297     {  5.6058e-01, -2.5201e-01,  7.3427e-01,  2.7537e+00 },
00298     {  2.4500e-01,  7.2200e-02,  5.2100e-01,  4.4500e+00 }
00299   };
00300 
00301   const double pbar_hadscale[7][4] = {
00302     //based on LogNormal of Had Energy 
00303     { -5.0239e+00,  1.2308e+01,  1.9628e-01,  0.0000e+00 },
00304     {  1.8635e+00,  1.1818e+00, -6.6838e-01,  1.8306e+00 },
00305     {  9.9975e-01,  6.0775e-01, -9.2005e-01,  2.8056e+00 },
00306     { -6.2579e-01,  2.9609e-01,  2.0549e+00,  3.7047e+00 },
00307     {  3.4519e+00,  2.9611e+00,  3.4592e-01,  3.5619e+00 },
00308     {  4.0909e-01,  3.6049e-01, -8.1641e-01,  1.0249e+00 },
00309     //correction of response if ssp is between Ecal and Hcal
00310     {  0.9800e-01,  0.8400e-01,  1.6000e+00,  4.2000e+00 }
00311   };
00312 
00313   const double pbar_par[8*NPar][4] = {
00314     {  1.4815e+00,  3.5510e+00,  1.5773e-02,  1.8235e+01 }, //Mean Ecal ShowerType1
00315     {  1.0727e+00,  4.7796e-01,  5.0331e-01,  3.2954e+00 },
00316     { -1.0941e+00,  2.3135e-01,  4.8154e-01,  3.8732e+00 },
00317     { -1.2430e-01,  1.5000e+00,  3.9037e-01,  7.8237e-01 },
00318     { -1.3253e+00,  3.6674e+00,  1.0286e-01, -9.7204e+00 },
00319 
00320     {  5.6542e-01, -5.0471e-02, -1.1991e+00,  4.8874e+00 }, //Mean Hcal ShowerType1
00321     {  9.3080e-01, -8.0336e-02,  6.0412e+00,  1.5598e+00 },
00322     { -1.5000e+00, -1.0766e+00,  2.3959e-01,  1.9593e+00 },
00323     {  1.4521e+00, -6.2515e-01,  1.2788e+00,  1.6684e+00 },
00324     {  2.0452e+00, -2.1115e+00,  4.4611e-01,  1.3040e+00 },
00325 
00326     {  6.0078e-01,  6.8636e-02,  8.6106e-01,  4.8556e+00 }, //Mean Hcal ShowerType2
00327     {  7.8057e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00328     { -1.4680e+00,  1.3644e+00, -1.8247e-01,  1.5390e+00 },
00329     {  1.3517e+00,  4.7462e-01, -1.9144e+00,  1.5090e+00 },
00330     {  2.0000e+00, -2.0000e+00,  4.8238e-01,  1.1081e+00 },
00331 
00332     {  6.6993e-01,  6.9898e-02,  9.8165e-01,  5.0369e+00 }, //Mean Hcal ShowerType3
00333     {  8.0659e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00334     { -1.7897e+00,  5.7671e-01, -5.0421e-01,  2.8906e+00 },
00335     {  9.9329e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00336     {  1.5000e+00, -1.5000e+00,  4.2544e-01,  1.1097e+00 },
00337 
00338     {  4.7782e-02,  1.1050e+00, -1.2951e-02,  1.0204e+01 }, //Sigma Ecal ShowerType1
00339     {  4.5173e-01,  4.0025e-02, -6.2352e+00,  1.7565e+00 },
00340     {  2.3372e-01,  7.2348e-01, -1.3646e-02,  1.2775e+01 },
00341     {  7.0000e-01,  3.5000e-01, -2.8551e-01,  3.2593e-01 },
00342     {  4.8372e-01,  2.0000e-01, -7.1584e-01,  1.0953e+00 },
00343 
00344     {  1.9193e-01,  1.2269e-01, -5.2125e-01,  2.4453e+00 }, //Sigma Hcal ShowerType1
00345     {  5.5009e-01, -6.2468e-02,  1.8419e+00,  4.7151e+00 },
00346     {  5.0000e-01,  3.0000e-01, -3.4911e-01,  3.8876e+00 },
00347     {  1.7835e+00,  1.4497e+00, -4.7752e-01,  8.1104e-01 },
00348     {  1.8000e+00, -1.6000e+00,  3.9842e-01,  7.3509e-01 },
00349 
00350     {  2.5000e-01,  1.7633e-01, -4.5662e-01,  1.4590e+00 }, //Sigma Hcal ShowerType2
00351     {  5.4673e-01,  7.1990e-02, -1.5031e+00,  4.4666e+00 },
00352     {  2.3564e+00, -2.3653e+00,  9.5422e-02, -7.6442e+00 },
00353     {  1.5848e+01,  1.5517e+01, -3.9489e-01, -3.0760e+00 },
00354     {  2.7839e+00,  2.5688e+00, -4.0084e-01, -2.1336e-01 },
00355 
00356     {  2.5000e-01,  1.9134e-01, -3.4221e-01,  1.0576e+00 }, //Sigma Hcal ShowerType3
00357     {  1.3747e-01,  1.5399e+00, -2.2446e-02,  1.4971e+01 },
00358     {  7.8008e-01, -2.4582e+00,  2.8835e-02, -7.0443e-01 },
00359     {  2.5066e+01,  2.4727e+01, -3.6163e-01, -4.6946e+00 },
00360     {  1.8343e+00,  1.6107e+00, -3.9030e-01,  1.9760e-02 }
00361   };
00362 
00363   const double pbar_rho[8*NPar][4] = {
00364     {  2.3079e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 }, //Ecal ShowerType1
00365     {  4.3290e-01,  6.7327e-02, -4.5756e+01,  3.2191e+00 },
00366     {  8.4086e-01,  4.8027e-02,  4.0865e+01,  3.1625e+00 },
00367     { -5.6374e-02,  1.0346e-01,  6.3263e+00,  3.1695e+00 },
00368     { -4.4616e+00, -5.0152e+00, -2.1702e-01, -5.4528e+00 },
00369     {  1.4026e-01,  1.3193e-01,  8.6221e+01,  3.2452e+00 },
00370     {  1.1522e-01, -7.5725e-02,  2.1818e+00,  1.8879e+00 },
00371     { -8.2358e+00,  8.7201e+00,  2.7548e-01, -5.6808e+00 },
00372     {  1.0653e+00, -1.8383e+00, -2.0251e-02,  2.5977e+01 },
00373     {  8.0521e-02,  8.8452e-01,  4.5910e-01,  1.1158e-02 },
00374 
00375     {  1.5216e-01,  1.4560e-01,  1.9375e+00,  2.7285e+00 }, //Hcal ShowerType1
00376     {  2.6865e-01,  1.0966e-01,  3.0722e+00,  2.2432e+00 },
00377     {  7.4470e-01,  0.0000e+00,  0.0000e+00,  0.0000e+00 },
00378     {  1.7389e-01,  1.6479e-01,  8.6705e-01,  4.2026e+00 },
00379     {  5.0000e-01,  2.7037e-01,  7.2902e-01,  2.7075e+00 },
00380     {  3.0480e-01,  1.5000e-01,  1.3249e+00,  1.4761e+00 },
00381     {  7.2012e-02,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00382     {  4.6644e-01,  1.8000e-01,  4.1184e-01,  5.2569e+00 },
00383     {  4.7765e-01,  1.5000e-01,  1.1606e+00,  1.0662e+00 },
00384     {  7.3257e-01,  1.1639e-01, -1.7675e+00,  2.8941e+00 },
00385 
00386     {  9.0554e-02,  8.2545e-02,  4.6059e+00,  2.2499e+00 }, //Hcal ShowerType2
00387     {  2.6407e-01,  6.4073e-02,  2.7583e+00,  1.6532e+00 },
00388     {  7.3779e-01, -8.2733e-03,  2.7742e+01,  2.6352e+00 },
00389     {  1.3900e-01,  1.0000e-01,  1.1953e+00,  4.2759e+00 },
00390     {  4.8511e-01,  2.1766e-01,  9.5147e-01,  2.7201e+00 },
00391     {  3.6706e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00392     {  2.0711e+01, -2.0659e+01,  7.2447e-02, -3.9754e+01 },
00393     {  2.0000e-01,  5.0335e+00,  1.0157e-02, -7.1268e-02 },
00394     {  5.9565e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00395     {  7.4264e-01, -9.7084e-02,  1.9748e+00,  2.5358e+00 },
00396 
00397     { -2.2297e+00,  4.4590e+00, -2.1916e-02,  3.0011e+01 }, //Hcal ShowerType3
00398     {  2.8783e-01,  1.7653e-01, -3.5284e+00,  3.3965e+00 },
00399     {  7.4446e-01,  0.0000e+00,  1.0000e+00,  1.0000e+00 },
00400     {  1.4560e-01, -1.0000e-01,  1.2909e+00,  1.2426e+00 },
00401     {  4.9473e-01,  1.1084e-01,  1.7734e+01,  3.3464e+00 },
00402     {  2.0462e+00, -1.8330e+00,  1.1470e+00,  1.0094e-01 },
00403     {  1.6943e-01, -1.1284e-01,  1.8951e+01,  3.3884e+00 },
00404     {  2.0000e-01,  4.2824e+00,  7.5278e-03, -4.3992e+00 },
00405     {  5.2574e-01, -7.5330e-02,  6.4920e-01,  5.0040e+00 },
00406     {  7.7100e-01, -7.6616e-02,  1.6244e+00,  1.9033e+00 }
00407   };
00408 
00409   const double pbar_rpar[4*Nrpar][2] = {
00410     {  1.3206e+00, -2.0591e-02}, //Ecal for showerType = 1
00411     {  8.7262e+00, -1.2565e+00},
00412     {  4.8914e-01,  3.1070e-02},
00413     {  2.8581e-03, -3.7858e-02},
00414     {  2.1375e+01, -2.9364e+00}, //Hcal for showerType = 1
00415     { -1.2541e+00,  1.2558e-01},
00416     {  4.9430e-01, -1.9658e-02},
00417     { -7.1026e-03,  6.7101e-03},
00418     {  1.2999e+01, -1.8209e+00}, //Hcal for showerType = 2
00419     { -1.5361e+00,  1.8317e-01},
00420     {  5.1010e-01, -2.6852e-02},
00421     { -2.1509e-02,  8.9109e-03},
00422     {  1.9082e+00, -1.1602e-02}, //Hcal for showerType = 3
00423     {  6.3817e+00, -1.1000e+00},
00424     {  6.0838e-01, -2.8689e-02},
00425     { -8.5457e-02, -1.7099e-03}
00426   };
00427 
00428   //HcalOuter parameters
00429   const double  ho_nonzero[4] = {4.79943e-01,4.61158e-01,7.09011e-01,4.86440e+00};   
00430 }
00431 
00432 #endif