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/APDShape.h"
0004 
0005 #include <cassert>
0006 
0007 void APDShape::fillShape(float& time_interval,
0008                          double& m_thresh,
0009                          EcalShapeBase::DVec& aVec,
0010                          const edm::EventSetup* es) const {
0011   if (m_useDBShape) {
0012     if (es == nullptr) {
0013       throw cms::Exception("EcalShapeBase:: DB conditions are not available, const edm::EventSetup* es == nullptr ");
0014     }
0015     auto const& esps = es->getData(espsToken_);
0016 
0017     aVec = esps.apd_shape;
0018     time_interval = esps.time_interval;
0019     m_thresh = esps.apd_thresh;
0020   } else {
0021     m_thresh = 0.0;
0022     time_interval = 1.0;
0023     aVec.reserve(500);
0024     const double m_tStart = 74.5;
0025     const double m_tau = 40.5;
0026 
0027     for (unsigned int i(0); i != 500; ++i) {
0028       const double ctime((1. * i + 0.5 - m_tStart) / m_tau);
0029       double val = 0 > ctime ? 0 : ctime * exp(1. - ctime);
0030       aVec.push_back(val);
0031     }
0032   }
0033 }