CSCStripConditions.cc

Go to the documentation of this file.

#include "SimMuon/CSCDigitizer/src/CSCStripConditions.h"

#include "SimGeneral/NoiseGenerators/interface/CorrelatedNoisifier.icc"

#include "CLHEP/Random/RandGaussQ.h"

template class CorrelatedNoisifier< CSCStripConditions::CSCCorrelatedNoiseMatrix > ;

template
void
CorrelatedNoisifier< CSCStripConditions::CSCCorrelatedNoiseMatrix
             >::noisify(CSCAnalogSignal&,
                                CLHEP::HepRandomEngine*,
                const std::vector<double>* rangau) const ;


CSCStripConditions::CSCStripConditions()
  : theNoisifier(nullptr)
 {}


CSCStripConditions::~CSCStripConditions() 
{
}


float CSCStripConditions::smearedGain(const CSCDetId & detId, int channel, CLHEP::HepRandomEngine* engine) const
{
  return CLHEP::RandGaussQ::shoot(engine, gain(detId, channel), gainSigma(detId, channel) );
}


void CSCStripConditions::noisify(const CSCDetId & detId, CSCAnalogSignal & signal, CLHEP::HepRandomEngine* engine)
{
  const int nScaBins = 8;
  const float scaBinSize = 50.;
  int channel = signal.getElement();
  std::vector<float> binValues(nScaBins, 0.);
  // use a temporary signal, in case we have to rebin
  CSCAnalogSignal tmpSignal(channel, scaBinSize, binValues,
                            0., signal.getTimeOffset());

  fetchNoisifier(detId, channel );
  theNoisifier->noisify(tmpSignal, engine);
  // noise matrix is in ADC counts. onvert to fC
  tmpSignal *= 1./gain(detId, channel);
  signal.superimpose(tmpSignal);
}

 
float CSCStripConditions::analogNoise(const CSCDetId & detId, int channel) const
{
  return sqrt(2) * pedestalSigma(detId, channel) / gain(detId, channel);
}