SiPixelVCalDB.cc

Go to the documentation of this file.

// system includes
#include <fstream>
#include <iostream>
#include <map>
#include <memory>
#include <string>

// user includes
#include "CondFormats/SiPixelObjects/interface/SiPixelVCal.h"
#include "CondCore/DBOutputService/interface/PoolDBOutputService.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/TrackerCommon/interface/PixelBarrelName.h"
#include "DataFormats/TrackerCommon/interface/PixelEndcapName.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "Geometry/CommonDetUnit/interface/PixelGeomDetUnit.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"

class SiPixelVCalDB : public edm::one::EDAnalyzer<> {
public:
  explicit SiPixelVCalDB(const edm::ParameterSet& conf);
  explicit SiPixelVCalDB();
  ~SiPixelVCalDB() override;
  void analyze(const edm::Event&, const edm::EventSetup&) override;

private:
  const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
  const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> tkTopoToken_;
  std::string recordName_;
  typedef std::vector<edm::ParameterSet> Parameters;
  Parameters BPixParameters_;
  Parameters FPixParameters_;
};

using namespace std;
using namespace edm;

SiPixelVCalDB::SiPixelVCalDB(edm::ParameterSet const& iConfig)
    : tkGeomToken_(esConsumes()), tkTopoToken_(esConsumes()) {
  recordName_ = iConfig.getUntrackedParameter<std::string>("record", "SiPixelVCalRcd");
  BPixParameters_ = iConfig.getUntrackedParameter<Parameters>("BPixParameters");
  FPixParameters_ = iConfig.getUntrackedParameter<Parameters>("FPixParameters");
}

SiPixelVCalDB::~SiPixelVCalDB() = default;

// Analyzer: Functions that gets called by framework every event
void SiPixelVCalDB::analyze(const edm::Event& e, const edm::EventSetup& iSetup) {
  SiPixelVCal vcal;
  bool phase1 = true;

  // Retrieve tracker topology from geometry
  const TrackerTopology* const tTopo = &iSetup.getData(tkTopoToken_);

  // Retrieve old style tracker geometry from geometry
  const TrackerGeometry* pDD = &iSetup.getData(tkGeomToken_);
  edm::LogPrint("SiPixelVCalDB") << " There are " << pDD->detUnits().size() << " modules" << std::endl;

  for (const auto& it : pDD->detUnits()) {
    if (dynamic_cast<PixelGeomDetUnit const*>(it) != nullptr) {
      const DetId detid = it->geographicalId();
      const unsigned int rawDetId = detid.rawId();
      int subid = detid.subdetId();

      // FILL BPIX
      if (subid == static_cast<int>(PixelSubdetector::PixelBarrel)) {
        int layer = tTopo->pxbLayer(detid);    // 1, 2, 3, 4
        int ladder = tTopo->pxbLadder(detid);  // 1-12/28/44/64
        edm::LogPrint("SiPixelVCalDB") << " pixel barrel:"
                                       << " detId=" << rawDetId << ", layer=" << layer << ", ladder=" << ladder;
        for (Parameters::iterator it = BPixParameters_.begin(); it != BPixParameters_.end(); ++it) {
          if (it->getParameter<int>("layer") == layer && it->getParameter<int>("ladder") == ladder) {
            float slope = (float)it->getParameter<double>("slope");
            float offset = (float)it->getParameter<double>("offset");
            edm::LogPrint("SiPixelVCalDB") << ";  VCal slope " << slope << ", offset " << offset;
            // edm::LogInfo("SiPixelVCalDB")  << "  detId " << rawDetId << " \t
            // VCal slope " << slope << ", offset " << offset;
            vcal.putSlopeAndOffset(detid, slope, offset);
          }
        }
        edm::LogPrint("SiPixelVCalDB") << std::endl;

        // FILL FPIX
      } else if (subid == static_cast<int>(PixelSubdetector::PixelEndcap)) {
        PixelEndcapName fpix(detid, tTopo, phase1);
        int side = tTopo->pxfSide(detid);   // 1 (-z), 2 for (+z)
        int disk = fpix.diskName();         // 1, 2, 3
        int disk2 = tTopo->pxfDisk(detid);  // 1, 2, 3
        int ring = fpix.ringName();         // 1 (lower), 2 (upper)
        if (disk != disk2) {
          edm::LogError("SiPixelVCalDB::analyze")
              << "Found contradicting FPIX disk number: " << disk << " vs." << disk2 << std::endl;
        }
        edm::LogPrint("SiPixelVCalDB") << " pixel endcap:"
                                       << " detId=" << rawDetId << ", side=" << side << ", disk=" << disk
                                       << ", ring=" << ring;
        for (Parameters::iterator it = FPixParameters_.begin(); it != FPixParameters_.end(); ++it) {
          if (it->getParameter<int>("side") == side && it->getParameter<int>("disk") == disk &&
              it->getParameter<int>("ring") == ring) {
            float slope = (float)it->getParameter<double>("slope");
            float offset = (float)it->getParameter<double>("offset");
            edm::LogPrint("SiPixelVCalDB") << ";  VCal slope " << slope << ", offset " << offset;
            // edm::LogInfo("SiPixelVCalDB")  << "  detId " << rawDetId << " \t
            // VCal slope " << slope << ", offset " << offset;
            vcal.putSlopeAndOffset(rawDetId, slope, offset);
          }
        }
        edm::LogPrint("SiPixelVCalDB") << std::endl;

      } else {
        edm::LogError("SiPixelVCalDB::analyze") << "detid is Pixel but neither bpix nor fpix" << std::endl;
      }
    }
  }

  // Save to DB
  edm::Service<cond::service::PoolDBOutputService> mydbservice;
  if (mydbservice.isAvailable()) {
    try {
      if (mydbservice->isNewTagRequest(recordName_)) {
        mydbservice->createOneIOV<SiPixelVCal>(vcal, mydbservice->beginOfTime(), recordName_);
      } else {
        mydbservice->appendOneIOV<SiPixelVCal>(vcal, mydbservice->currentTime(), recordName_);
      }
    } catch (const cond::Exception& er) {
      edm::LogError("SiPixelVCalDB") << er.what() << std::endl;
    } catch (const std::exception& er) {
      edm::LogError("SiPixelVCalDB") << "caught std::exception " << er.what() << std::endl;
    } catch (...) {
      edm::LogError("SiPixelVCalDB") << "Funny error" << std::endl;
    }
  } else {
    edm::LogError("SiPixelVCalDB") << "Service is unavailable" << std::endl;
  }
}
DEFINE_FWK_MODULE(SiPixelVCalDB);