SiStripNoisesFakeESSource.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    CalibTracker/SiStripESProducers
// Class:      SiStripNoisesFakeESSource
//
// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"

#include "CondFormats/SiStripObjects/interface/SiStripNoises.h"
#include "CondFormats/DataRecord/interface/SiStripCondDataRecords.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "SiStripFakeAPVParameters.h"

class SiStripNoisesFakeESSource : public edm::ESProducer, public edm::EventSetupRecordIntervalFinder {
public:
  SiStripNoisesFakeESSource(const edm::ParameterSet&);
  ~SiStripNoisesFakeESSource() override;

  void setIntervalFor( const edm::eventsetup::EventSetupRecordKey&, const edm::IOVSyncValue& iov, edm::ValidityInterval& iValidity ) override;

  typedef std::shared_ptr<SiStripNoises> ReturnType;
  ReturnType produce(const SiStripNoisesRcd&);

private:
  bool m_stripLengthMode;
  double m_noisePar0;
  SiStripFakeAPVParameters m_noisePar1;
  SiStripFakeAPVParameters m_noisePar2;
  edm::FileInPath m_file;
  uint32_t m_printDebug;
};

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"

#include "CLHEP/Random/RandGauss.h"

namespace { // helper methods
  inline void printLog(const uint32_t detId, const unsigned short strip, const double & noise)
  {
    edm::LogInfo("SiStripNoisesDummyCalculator") << "detid: " << detId << " strip: " << strip <<  " noise: " << noise;
  }
}

SiStripNoisesFakeESSource::SiStripNoisesFakeESSource(const edm::ParameterSet& iConfig)
{
  setWhatProduced(this);
  findingRecord<SiStripNoisesRcd>();

  m_stripLengthMode = iConfig.getParameter<bool>("StripLengthMode");

  if ( ! m_stripLengthMode ) {
    //parameters for random noise generation. not used if Strip length mode is chosen
    m_noisePar0 = iConfig.getParameter<double>("MinPositiveNoise");
    m_noisePar1 = SiStripFakeAPVParameters(iConfig, "MeanNoise");
    m_noisePar2 = SiStripFakeAPVParameters(iConfig, "SigmaNoise");
  } else {
    //parameters for strip length proportional noise generation. not used if random mode is chosen
    m_noisePar0 = iConfig.getParameter<double>("electronPerAdc");
    m_noisePar1 = SiStripFakeAPVParameters(iConfig, "NoiseStripLengthSlope");
    m_noisePar2 = SiStripFakeAPVParameters(iConfig, "NoiseStripLengthQuote");
  }

  m_file = iConfig.getParameter<edm::FileInPath>("file");
  m_printDebug = iConfig.getUntrackedParameter<uint32_t>("printDebug", 5);
}

SiStripNoisesFakeESSource::~SiStripNoisesFakeESSource() {}

void SiStripNoisesFakeESSource::setIntervalFor( const edm::eventsetup::EventSetupRecordKey&, const edm::IOVSyncValue& iov, edm::ValidityInterval& iValidity )
{
  iValidity = edm::ValidityInterval{iov.beginOfTime(), iov.endOfTime()};
}

// ------------ method called to produce the data  ------------
SiStripNoisesFakeESSource::ReturnType
SiStripNoisesFakeESSource::produce(const SiStripNoisesRcd& iRecord)
{
  using namespace edm::es;

  edm::ESHandle<TrackerTopology> tTopo;
  iRecord.getRecord<TrackerTopologyRcd>().get(tTopo);

  std::shared_ptr<SiStripNoises> noises{new SiStripNoises};

  SiStripDetInfoFileReader reader{m_file.fullPath()};

  uint32_t count{0};
  for ( const auto& elm : reader.getAllData() ) {
    //Generate Noises for det detid
    SiStripNoises::InputVector theSiStripVector;
    SiStripFakeAPVParameters::index sl = SiStripFakeAPVParameters::getIndex(tTopo.product(), elm.first);

    if ( m_stripLengthMode ) {
      // Use strip length
      const double linearSlope{m_noisePar1.get(sl)};
      const double linearQuote{m_noisePar2.get(sl)};
      const double stripLength{elm.second.stripLength};
      for ( unsigned short j{0}; j < 128*elm.second.nApvs; ++j ) {
        const float noise = (linearSlope*stripLength + linearQuote) / m_noisePar0;
        if ( count < m_printDebug ) printLog(elm.first, j, noise);
        noises->setData(noise, theSiStripVector);
      }
    } else {
      // Use random generator
      const double meanN {m_noisePar1.get(sl)};
      const double sigmaN{m_noisePar2.get(sl)};
      for ( unsigned short j{0}; j < 128*elm.second.nApvs; ++j ) {
        const float noise = std::max(CLHEP::RandGauss::shoot(meanN, sigmaN), m_noisePar0);
        if ( count < m_printDebug ) printLog(elm.first, j, noise);
        noises->setData(noise, theSiStripVector);
      }
    }
    ++count;

    if ( ! noises->put(elm.first, theSiStripVector) ) {
      edm::LogError("SiStripNoisesFakeESSource::produce ") << " detid already exists";
    }
  }

  return noises;
}

//define this as a plug-in
#include "FWCore/Framework/interface/SourceFactory.h"
DEFINE_FWK_EVENTSETUP_SOURCE(SiStripNoisesFakeESSource);