SiPixelGenErrorDBObjectESProducer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiPixelGenErrorDBObjectESProducer
// Class:      SiPixelGenErrorDBObjectESProducer
//
//
// Original Author:  D.Fehling
//         Created:  Tue Sep 29 14:49:31 CET 2009
//
//
 
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/ESProductTag.h"
#include "FWCore/Utilities/interface/do_nothing_deleter.h"
#include "FWCore/Framework/interface/ModuleFactory.h"
 
#include "CondFormats/SiPixelObjects/interface/SiPixelGenErrorDBObject.h"
#include "CalibTracker/Records/interface/SiPixelGenErrorDBObjectESProducerRcd.h"
#include "MagneticField/Engine/interface/MagneticField.h"
 
#include <memory>
 
using namespace edm;
 
class SiPixelGenErrorDBObjectESProducer : public edm::ESProducer {
public:
  SiPixelGenErrorDBObjectESProducer(const edm::ParameterSet& iConfig);
  std::shared_ptr<const SiPixelGenErrorDBObject> produce(const SiPixelGenErrorDBObjectESProducerRcd&);
 
private:
  edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> magfieldToken_;
  edm::ESGetToken<SiPixelGenErrorDBObject, SiPixelGenErrorDBObjectRcd> genErrorToken_;
};
 
SiPixelGenErrorDBObjectESProducer::SiPixelGenErrorDBObjectESProducer(const edm::ParameterSet& iConfig) {
  auto cc = setWhatProduced(this);
  cc.setMayConsume(
      genErrorToken_,
      [](const auto& get, edm::ESTransientHandle<MagneticField> iMagfield) {
        const GlobalPoint center(0.0, 0.0, 0.0);
        const float theMagField = iMagfield->inTesla(center).mag();
        if (theMagField >= -0.1 && theMagField < 1.0)
          return get("", "0T");
        else if (theMagField >= 1.0 && theMagField < 2.5)
          return get("", "2T");
        else if (theMagField >= 2.5 && theMagField < 3.25)
          return get("", "3T");
        else if (theMagField >= 3.25 && theMagField < 3.65)
          return get("", "35T");
        else if (theMagField >= 3.9 && theMagField < 4.1)
          return get("", "4T");
        else {
          if (theMagField >= 4.1 || theMagField < -0.1)
            edm::LogWarning("UnexpectedMagneticFieldUsingDefaultPixelGenError") << "Magnetic field is " << theMagField;
          //return get("", "3.8T");
          return get("", "");
        }
      },
      edm::ESProductTag<MagneticField, IdealMagneticFieldRecord>("", ""));
  magfieldToken_ = cc.consumes();
}
 
std::shared_ptr<const SiPixelGenErrorDBObject> SiPixelGenErrorDBObjectESProducer::produce(
    const SiPixelGenErrorDBObjectESProducerRcd& iRecord) {
  const GlobalPoint center(0.0, 0.0, 0.0);
  const float theMagField = iRecord.get(magfieldToken_).inTesla(center).mag();
 
  const auto& dbobject = iRecord.get(genErrorToken_);
 
  if (std::fabs(theMagField - dbobject.sVector()[22]) > 0.1)
    edm::LogWarning("UnexpectedMagneticFieldUsingNonIdealPixelGenError")
        << "Magnetic field is " << theMagField << " GenError Magnetic field is " << dbobject.sVector()[22];
 
  return std::shared_ptr<const SiPixelGenErrorDBObject>(&dbobject, edm::do_nothing_deleter());
}
 
DEFINE_FWK_EVENTSETUP_MODULE(SiPixelGenErrorDBObjectESProducer);