APDShape.cc

Go to the documentation of this file.

#include <cmath>

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

#include<cassert>

void
APDShape::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->apd_shape; 
    time_interval = esps->time_interval;
    m_thresh = esps->apd_thresh;
    }
    else
    {
    m_thresh = 0.0;
        time_interval = 1.0; 
        aVec.reserve(500);
        const double m_tStart=74.5;
        const double m_tau=40.5;

    for( unsigned int i ( 0 ) ; i != 500 ; ++i )
    {
        const double ctime ( ( 1.*i + 0.5 - m_tStart )/m_tau ) ;
            double val =  0 > ctime ? 0 : ctime * exp( 1. - ctime );
        aVec.push_back(val);
    }
    }
}