AlcaBeamSpotFromDB.cc

Go to the documentation of this file.

// C++ standard
#include <string>
// CMS
#include "Calibration/TkAlCaRecoProducers/interface/AlcaBeamSpotFromDB.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/IOVSyncValue.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"

AlcaBeamSpotFromDB::AlcaBeamSpotFromDB(const edm::ParameterSet &iConfig)
    : beamSpotToken_(esConsumes<BeamSpotObjects, BeamSpotObjectsRcd, edm::Transition::EndLuminosityBlock>()) {
  produces<reco::BeamSpot, edm::Transition::EndLuminosityBlock>("alcaBeamSpot");
}

AlcaBeamSpotFromDB::~AlcaBeamSpotFromDB() {}

//--------------------------------------------------------------------------------------------------
void AlcaBeamSpotFromDB::produce(edm::Event &iEvent, const edm::EventSetup &iSetup) {}

//--------------------------------------------------------------------------------------------------
void AlcaBeamSpotFromDB::endLuminosityBlockProduce(edm::LuminosityBlock &lumiSeg, const edm::EventSetup &iSetup) {
  // read DB object
  const BeamSpotObjects *spotDB = &iSetup.getData(beamSpotToken_);

  // translate from BeamSpotObjects to reco::BeamSpot
  reco::BeamSpot::Point apoint(spotDB->x(), spotDB->y(), spotDB->z());

  reco::BeamSpot::CovarianceMatrix matrix;
  for (int i = 0; i < 7; ++i) {
    for (int j = 0; j < 7; ++j) {
      matrix(i, j) = spotDB->covariance(i, j);
    }
  }

  reco::BeamSpot aSpot;
  // this assume beam width same in x and y
  aSpot = reco::BeamSpot(apoint, spotDB->sigmaZ(), spotDB->dxdz(), spotDB->dydz(), spotDB->beamWidthX(), matrix);
  aSpot.setBeamWidthY(spotDB->beamWidthY());
  aSpot.setEmittanceX(spotDB->emittanceX());
  aSpot.setEmittanceY(spotDB->emittanceY());
  aSpot.setbetaStar(spotDB->betaStar());

  if (spotDB->beamType() == 2) {
    aSpot.setType(reco::BeamSpot::Tracker);
  } else {
    aSpot.setType(reco::BeamSpot::Fake);
  }

  auto result = std::make_unique<reco::BeamSpot>();
  *result = aSpot;
  lumiSeg.put(std::move(result), std::string("alcaBeamSpot"));

  // std::cout << " for runs: " << iEvent.id().run() << " - " <<
  // iEvent.id().run() << std::endl;
  std::cout << aSpot << std::endl;
}

void AlcaBeamSpotFromDB::beginJob() {}

void AlcaBeamSpotFromDB::endJob() {}

// define this as a plug-in
DEFINE_FWK_MODULE(AlcaBeamSpotFromDB);