EBShape.cc

Go to the documentation of this file.

#include <cmath>

#include "SimCalorimetry/EcalSimAlgos/interface/EBShape.h"



void
EBShape::fillShape( float &time_interval, double &m_thresh, EcalShapeBase::DVec& aVec, const edm::EventSetup* es) const
{
    if(m_useDBShape) 
    {
        if(es == nullptr)
        {
        throw cms::Exception("EcalShapeBase:: DB conditions are not available, const edm::EventSetup* es == nullptr ");
    }
    edm::ESHandle<EcalSimPulseShape> esps;
    es->get<EcalSimPulseShapeRcd>().get(esps);

    aVec = esps->barrel_shape; 
    time_interval = esps->time_interval;
    m_thresh = esps->barrel_thresh;

    } else   // use old hardcoded arrays
    {
        m_thresh = 0.00013;
        time_interval = 1.0; 
        aVec.reserve(500);
        aVec = {6.94068e-05,
                -5.03304e-05, 
        -2.13404e-05, 
        6.017e-05, 
            2.01697e-05, 
            0.000114845, 
            2.13998e-05, 
            2.74476e-05, 
            5.2824e-05, 
            8.754e-05, 
            2.95346e-06, 
            -7.58699e-05, 
            -2.72224e-05, 
            3.10997e-06, 
            -3.97771e-05, 
            -1.06916e-05, 
            -0.000113865, 
            6.05044e-05, 
            -5.81202e-05, 
            -6.58974e-06, 
            5.37494e-05, 
            -0.000123729, 
            7.50938e-06, 
            -1.35628e-05, 
            8.33725e-05, 
            3.19299e-05, 
            -3.09232e-05, 
            -7.0086e-05, 
            1.78937e-06, 
            -2.20365e-05, 
            7.68054e-05, 
            -2.5368e-05, 
            5.67291e-06, 
            5.87096e-05, 
            -2.62771e-06, 
            4.31832e-05, 
            8.33616e-06, 
            7.27813e-05, 
            7.6159e-05, 
            -1.60446e-05, 
            -4.12127e-06, 
            -5.93381e-05, 
            1.61444e-05, 
            -5.49559e-05, 
            5.55254e-05, 
            3.32251e-05, 
            -3.15897e-05, 
            7.86588e-05, 
            -2.9704e-05, 
            5.66838e-05, 
            2.85281e-05, 
            -3.02436e-05, 
            -4.16265e-05, 
            -1.63191e-05, 
            6.61193e-05, 
            9.23766e-05, 
            6.68903e-05, 
            -3.20994e-05, 
            0.00011082, 
            -4.07997e-05, 
            -8.29046e-06, 
            -7.42197e-05, 
            -1.64386e-05, 
            1.02508e-05, 
            7.10995e-06, 
            -5.87486e-05, 
            -0.000101201, 
            1.62003e-05, 
            -2.53093e-05, 
            2.65239e-05, 
            -2.68722e-05, 
            -4.02001e-05, 
            5.0674e-05, 
            -1.75884e-05, 
            4.7902e-05, 
            -1.01079e-05, 
            1.08427e-05, 
            -0.000112906, 
            3.33076e-05, 
            0.000181201, 
            0.000426875, 
            0.00114222, 
            0.00237804, 
            0.00541858, 
            0.0089021, 
            0.0149157, 
            0.0231397, 
            0.0344671, 
            0.0471013, 
            0.0625517, 
            0.0857351, 
            0.108561, 
            0.133481, 
            0.163557, 
            0.200243, 
            0.225919, 
            0.269213, 
            0.302929, 
            0.342722, 
            0.378522, 
            0.436563, 
            0.467581, 
            0.510133, 
            0.550063, 
            0.583509, 
            0.619187, 
            0.653245, 
            0.686101, 
            0.721178, 
            0.745129, 
            0.774163, 
            0.799011, 
            0.822177, 
            0.838315, 
            0.858847, 
            0.875559, 
            0.891294, 
            0.90537, 
            0.919617, 
            0.930632, 
            0.936216, 
            0.947739, 
            0.955306, 
            0.961876, 
            0.968124, 
            0.97327, 
            0.977513, 
            0.984885, 
            0.986497, 
            0.990039, 
            0.994798, 
            0.994884, 
            0.99795, 
            0.99834, 
            0.999607, 
            1, 
            0.999047, 
            0.998745, 
            0.999219, 
            0.99814, 
            0.995082, 
            0.992449, 
            0.990418, 
            0.985032, 
            0.982308, 
            0.978696, 
            0.975656, 
            0.971027, 
            0.964811, 
            0.959428, 
            0.95096, 
            0.947428, 
            0.9419, 
            0.933223, 
            0.926482, 
            0.922172, 
            0.912777, 
            0.907388, 
            0.897289, 
            0.891889, 
            0.882056, 
            0.873382, 
            0.865442, 
            0.860032, 
            0.85202, 
            0.841013, 
            0.833802, 
            0.825259, 
            0.815013, 
            0.807465, 
            0.799428, 
            0.792165, 
            0.783088, 
            0.773392, 
            0.764982, 
            0.752174, 
            0.746487, 
            0.737678, 
            0.727396, 
            0.718692, 
            0.712737, 
            0.702738, 
            0.69559, 
            0.684389, 
            0.677989, 
            0.667643, 
            0.659009, 
            0.650217, 
            0.644479, 
            0.636017, 
            0.625257, 
            0.618507, 
            0.609798, 
            0.600097, 
            0.592788, 
            0.584895, 
            0.578228, 
            0.569299, 
            0.560576, 
            0.552404, 
            0.541405, 
            0.536271, 
            0.528734, 
            0.519813, 
            0.512264, 
            0.507001, 
            0.49828, 
            0.492416, 
            0.483181, 
            0.477907, 
            0.469623, 
            0.462528, 
            0.455099, 
            0.45055, 
            0.443576, 
            0.435364, 
            0.429789, 
            0.422724, 
            0.415621, 
            0.409469, 
            0.40401, 
            0.398121, 
            0.391079, 
            0.384414, 
            0.378214, 
            0.369851, 
            0.365966, 
            0.359865, 
            0.353505, 
            0.347899, 
            0.343829, 
            0.337585, 
            0.333089, 
            0.326289, 
            0.322249, 
            0.316079, 
            0.31061, 
            0.305426, 
            0.301885, 
            0.296753, 
            0.290931, 
            0.286877, 
            0.281831, 
            0.276633, 
            0.272283, 
            0.268069, 
            0.26399, 
            0.258457, 
            0.253549, 
            0.249493}; 

    for( unsigned int i ( 250 ) ; i != 500 ; ++i ) aVec.push_back( exp(2.39735 - 0.0151053* ((double)i+1.0)) );
    }
}