SiPixelTemplateDBObjectESProducer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiPixelTemplateDBObjectESProducer
// Class:      SiPixelTemplateDBObjectESProducer
//
//
// Original Author:  D.Fehling
//         Created:  Tue Sep 29 14:49:31 CET 2009
//
//
 
#include <memory>
 
#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 "MagneticField/Engine/interface/MagneticField.h"
 
#include "CondFormats/SiPixelObjects/interface/SiPixelTemplateDBObject.h"
#include "CalibTracker/Records/interface/SiPixelTemplateDBObjectESProducerRcd.h"
 
using namespace edm;
 
class SiPixelTemplateDBObjectESProducer : public edm::ESProducer {
public:
  SiPixelTemplateDBObjectESProducer(const edm::ParameterSet& iConfig);
  std::shared_ptr<const SiPixelTemplateDBObject> produce(const SiPixelTemplateDBObjectESProducerRcd&);
 
private:
  edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> magfieldToken_;
  edm::ESGetToken<SiPixelTemplateDBObject, SiPixelTemplateDBObjectRcd> templateToken_;
};
 
SiPixelTemplateDBObjectESProducer::SiPixelTemplateDBObjectESProducer(const edm::ParameterSet& iConfig) {
  auto cc = setWhatProduced(this);
  cc.setMayConsume(
      templateToken_,
      [](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("UnexpectedMagneticFieldUsingDefaultPixelTemplate") << "Magnetic field is " << theMagField;
          //return get("", "3.8T");
          return get("", "");
        }
      },
      edm::ESProductTag<MagneticField, IdealMagneticFieldRecord>("", ""));
  magfieldToken_ = cc.consumes();
}
 
std::shared_ptr<const SiPixelTemplateDBObject> SiPixelTemplateDBObjectESProducer::produce(
    const SiPixelTemplateDBObjectESProducerRcd& iRecord) {
  const GlobalPoint center(0.0, 0.0, 0.0);
  const float theMagField = iRecord.get(magfieldToken_).inTesla(center).mag();
 
  const auto& dbobject = iRecord.get(templateToken_);
 
  if (std::fabs(theMagField - dbobject.sVector()[22]) > 0.1)
    edm::LogWarning("UnexpectedMagneticFieldUsingNonIdealPixelTemplate")
        << "Magnetic field is " << theMagField << " Template Magnetic field is " << dbobject.sVector()[22];
 
  return std::shared_ptr<const SiPixelTemplateDBObject>(&dbobject, edm::do_nothing_deleter());
}
 
DEFINE_FWK_EVENTSETUP_MODULE(SiPixelTemplateDBObjectESProducer);