Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-09-02 03:46:08

0001 #include <cmath>
0002 
0003 #include "SimCalorimetry/EcalSimAlgos/interface/EEShape.h"
0004 
0005 void EEShape::fillShape(float& time_interval,
0006                         double& m_thresh,
0007                         EcalShapeBase::DVec& aVec,
0008                         const edm::EventSetup* es) const {
0009   if (m_useDBShape) {
0010     if (es == nullptr) {
0011       throw cms::Exception("EcalShapeBase:: DB conditions are not available, const edm::EventSetup* es == nullptr ");
0012     }
0013     auto const& esps = es->getData(espsToken_);
0014 
0015     aVec = esps.endcap_shape;
0016     time_interval = esps.time_interval;
0017     m_thresh = esps.endcap_thresh;
0018   } else {
0019     m_thresh = 0.00025;
0020     time_interval = 1.0;
0021     aVec.reserve(500);
0022     aVec = {
0023         9.09091e-05,  9.96016e-05,  8.08219e-05,  0.000119685,  9.56522e-05,  0.000143969,  0.000101639,  0.000115789,
0024         8.31461e-05,  0.000117091,  0.000141224,  5.38745e-05,  0.000158079,  0.000104906,  7.91946e-05,  8e-05,
0025         0.000119703,  9.96516e-05,  8.82591e-05,  0.000146988,  0.000120588,  9.89474e-05,  0.000116788,  0.000151292,
0026         0.000108434,  0.000108651,  8.5259e-05,   0.000146575,  8.66142e-05,  0.000103679,  0.000153307,  7.80328e-05,
0027         8.42105e-05,  0.000114607,  6.90909e-05,  9.22449e-05,  8.92989e-05,  0.000105677,  9.58491e-05,  0.000107383,
0028         8.85714e-05,  0.000119703,  0.000141463,  0.000122267,  6.0241e-05,   0.000120588,  0.000115789,  0.000121168,
0029         8.04428e-05,  9.39759e-05,  3.391e-05,    -7.72908e-05, 5.61644e-05,  0.000181102,  -3.67893e-05, 6.45914e-05,
0030         2.68852e-05,  -1.57895e-05, -2.47191e-05, -5.30909e-05, -1.06122e-05, -3.69004e-06, 5.85153e-05,  -3.54717e-05,
0031         -9.39597e-06, 5e-05,        -3.19703e-05, 3.69338e-05,  -0.000169231, -4.09639e-05, -1.61765e-05, -8.21053e-05,
0032         6.86131e-05,  -6.12546e-05, -6.50602e-05, 2.56055e-05,  -3.90438e-05, -5.89041e-05, 0.000195276,  -7.29097e-05,
0033         1.32296e-05,  3.27869e-06,  -6.31579e-05, -3.37079e-05, 5.16364e-05,  -4.97959e-05, 8.92989e-05,  0.000561572,
0034         0.00171698,   0.0035906,    0.00742143,   0.01249,      0.020922,     0.030336,     0.043588,     0.0593529,
0035         0.0764947,    0.0958978,    0.11567,      0.137087,     0.160643,     0.188973,     0.218062,     0.240682,
0036         0.265773,     0.293327,     0.320832,     0.348689,     0.374411,     0.403575,     0.431147,     0.458644,
0037         0.488217,     0.515254,     0.53902,      0.566819,     0.589752,     0.615587,     0.63808,      0.664137,
0038         0.687959,     0.707379,     0.726826,     0.747043,     0.762952,     0.781342,     0.799515,     0.818559,
0039         0.832301,     0.846658,     0.860791,     0.874009,     0.886242,     0.896982,     0.908713,     0.9189,
0040         0.9283,       0.938453,     0.945995,     0.95186,      0.958876,     0.963283,     0.968989,     0.972846,
0041         0.977824,     0.981341,     0.982743,     0.984082,     0.986609,     0.985925,     0.984707,     0.981384,
0042         0.980941,     0.980132,     0.978176,     0.975987,     0.973396,     0.970395,     0.967422,     0.96346,
0043         0.959622,     0.955396,     0.950846,     0.945995,     0.940835,     0.935373,     0.929544,     0.923836,
0044         0.91787,      0.910812,     0.903924,     0.897669,     0.89127,      0.884613,     0.878051,     0.867734,
0045         0.856756,     0.848477,     0.842189,     0.833787,     0.825903,     0.817931,     0.8102,       0.802771,
0046         0.79388,      0.785348,     0.777295,     0.767854,     0.759619,     0.752358,     0.743394,     0.736674,
0047         0.726982,     0.719321,     0.708937,     0.700003,     0.693202,     0.685961,     0.676889,     0.670513,
0048         0.658788,     0.648343,     0.639381,     0.632761,     0.623614,     0.616164,     0.607884,     0.601437,
0049         0.593944,     0.585817,     0.577248,     0.570328,     0.56147,      0.553509,     0.548765,     0.539964,
0050         0.53343,      0.525787,     0.519868,     0.510817,     0.503274,     0.498478,     0.493058,     0.484176,
0051         0.479896,     0.469998,     0.462893,     0.45576,      0.450375,     0.442683,     0.437887,     0.430736,
0052         0.427216,     0.420775,     0.41485,      0.407769,     0.403546,     0.397028,     0.389534,     0.388569,
0053         0.380296,     0.375312,     0.369457,     0.365884,     0.358889,     0.353344,     0.350364,     0.347136,
0054         0.33872,      0.336706};
0055 
0056     for (unsigned int i(250); i != 500; ++i)
0057       aVec.push_back(exp(2.63163e-01 - 1.37795e-02 * (i - 151.0)));
0058   }
0059 }