BeamSpotOnlineProducer Class Reference

#include <BeamSpotOnlineProducer.h>

Inheritance diagram for BeamSpotOnlineProducer:

Public Types
typedef std::vector< edm::ParameterSet >	Parameters
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Member Functions
	BeamSpotOnlineProducer (const edm::ParameterSet &iConf)
	constructor More...
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	produce a beam spot class More...
	~BeamSpotOnlineProducer () override
	destructor More...
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
const bool	changeFrame_
const edm::EDGetTokenT< L1GlobalTriggerEvmReadoutRecord >	l1GtEvmReadoutRecordToken_
const edm::EDGetTokenT< BeamSpotOnlineCollection >	scalerToken_
const unsigned int	theBeamShoutMode
double	theMaxR2
const double	theMaxZ
const double	theSetSigmaZ

Detailed Description

---

class: BeamSpotOnlineProducer.h package: RecoVertex/BeamSpotProducer

author: Francisco Yumiceva, Fermilab (yumic.nosp@m.eva@.nosp@m.fnal..nosp@m.gov)

---

Definition at line 23 of file BeamSpotOnlineProducer.h.

Member Typedef Documentation

typedef std::vector<edm::ParameterSet> BeamSpotOnlineProducer::Parameters

Definition at line 25 of file BeamSpotOnlineProducer.h.

Constructor & Destructor Documentation

BeamSpotOnlineProducer::BeamSpotOnlineProducer ( const edm::ParameterSet &  iConf )

explicit

constructor

Definition at line 15 of file BeamSpotOnlineProducer.cc.

References edm::ParameterSet::getParameter(), and theMaxR2.

    : changeFrame_(iconf.getParameter<bool>("changeToCMSCoordinates")),
      theMaxZ(iconf.getParameter<double>("maxZ")),
      theSetSigmaZ(iconf.getParameter<double>("setSigmaZ")),
      scalerToken_(consumes<BeamSpotOnlineCollection>(iconf.getParameter<InputTag>("src"))),
      l1GtEvmReadoutRecordToken_(consumes<L1GlobalTriggerEvmReadoutRecord>(iconf.getParameter<InputTag>("gtEvmLabel"))),
      theBeamShoutMode(iconf.getUntrackedParameter<unsigned int>("beamMode", 11)) {
  theMaxR2 = iconf.getParameter<double>("maxRadius");
  theMaxR2 *= theMaxR2;

  produces<reco::BeamSpot>();
}

BeamSpotOnlineProducer::~BeamSpotOnlineProducer ( )

override

destructor

Definition at line 28 of file BeamSpotOnlineProducer.cc.

{}

Member Function Documentation

void BeamSpotOnlineProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

produce a beam spot class

Definition at line 30 of file BeamSpotOnlineProducer.cc.

References L1GtfeExtWord::beamMode(), align::BeamSpot, changeFrame_, DEFINE_FWK_MODULE, BeamSpotOnline::dxdz(), BeamSpotOnline::dydz(), BeamSpotOnline::err_sigma_z(), BeamSpotOnline::err_x(), BeamSpotOnline::err_y(), BeamSpotOnline::err_z(), f, edm::EventSetup::get(), BeamSpotObjects::GetBeamWidthX(), BeamSpotObjects::GetBeamWidthY(), BeamSpotObjects::GetBetaStar(), edm::Event::getByToken(), BeamSpotObjects::GetCovariance(), BeamSpotObjects::Getdxdz(), BeamSpotObjects::Getdydz(), BeamSpotObjects::GetEmittanceX(), BeamSpotObjects::GetEmittanceY(), BeamSpotObjects::GetSigmaZ(), BeamSpotObjects::GetX(), BeamSpotObjects::GetY(), BeamSpotObjects::GetZ(), L1GlobalTriggerEvmReadoutRecord::gtfeWord(), mps_fire::i, dqmiolumiharvest::j, l1GtEvmReadoutRecordToken_, reco::BeamSpot::LHC, makeMuonMisalignmentScenario::matrix, eostools::move(), edm::ESHandle< T >::product(), edm::Event::put(), diffTwoXMLs::r2, mps_fire::result, scalerToken_, BeamSpotOnline::sigma_z(), fftjetvertexadder_cfi::sigmaZ, theBeamShoutMode, theMaxR2, theMaxZ, theSetSigmaZ, reco::BeamSpot::Tracker, BeamSpotOnline::width_x(), BeamSpotOnline::width_y(), BeamSpotOnline::x(), BeamSpotOnline::y(), and BeamSpotOnline::z().

                                                                            {
  //shout MODE only in stable beam
  bool shoutMODE = false;
  edm::Handle<L1GlobalTriggerEvmReadoutRecord> gtEvmReadoutRecord;
  if (iEvent.getByToken(l1GtEvmReadoutRecordToken_, gtEvmReadoutRecord)) {
    if (gtEvmReadoutRecord->gtfeWord().beamMode() == theBeamShoutMode)
      shoutMODE = true;
  } else {
    shoutMODE = true;
  }

  // get scalar collection
  Handle<BeamSpotOnlineCollection> handleScaler;
  iEvent.getByToken(scalerToken_, handleScaler);

  // beam spot scalar object
  BeamSpotOnline spotOnline;

  // product is a reco::BeamSpot object
  auto result = std::make_unique<reco::BeamSpot>();

  reco::BeamSpot aSpot;

  bool fallBackToDB = false;
  if (!handleScaler->empty()) {
    // get one element
    spotOnline = *(handleScaler->begin());

    // in case we need to switch to LHC reference frame
    // ignore for the moment rotations, and translations
    double f = 1.;
    if (changeFrame_)
      f = -1.;

    reco::BeamSpot::Point apoint(f * spotOnline.x(), spotOnline.y(), f * spotOnline.z());

    reco::BeamSpot::CovarianceMatrix matrix;
    matrix(0, 0) = spotOnline.err_x() * spotOnline.err_x();
    matrix(1, 1) = spotOnline.err_y() * spotOnline.err_y();
    matrix(2, 2) = spotOnline.err_z() * spotOnline.err_z();
    matrix(3, 3) = spotOnline.err_sigma_z() * spotOnline.err_sigma_z();

    double sigmaZ = spotOnline.sigma_z();
    if (theSetSigmaZ > 0)
      sigmaZ = theSetSigmaZ;

    aSpot = reco::BeamSpot(apoint, sigmaZ, spotOnline.dxdz(), f * spotOnline.dydz(), spotOnline.width_x(), matrix);

    aSpot.setBeamWidthY(spotOnline.width_y());
    aSpot.setEmittanceX(0.);
    aSpot.setEmittanceY(0.);
    aSpot.setbetaStar(0.);
    aSpot.setType(reco::BeamSpot::LHC);  // flag value from scalars

    // check if we have a valid beam spot fit result from online DQM
    if (spotOnline.x() == 0 && spotOnline.y() == 0 && spotOnline.z() == 0 && spotOnline.width_x() == 0 &&
        spotOnline.width_y() == 0) {
      if (shoutMODE) {
        edm::LogWarning("BeamSpotFromDB")
            << "Online Beam Spot producer falls back to DB value because the scaler values are zero ";
      }
      fallBackToDB = true;
    }
    double r2 = spotOnline.x() * spotOnline.x() + spotOnline.y() * spotOnline.y();
    if (fabs(spotOnline.z()) >= theMaxZ || r2 >= theMaxR2) {
      if (shoutMODE) {
        edm::LogError("BeamSpotFromDB")
            << "Online Beam Spot producer falls back to DB value because the scaler values are too big to be true :"
            << spotOnline.x() << " " << spotOnline.y() << " " << spotOnline.z();
      }
      fallBackToDB = true;
    }
  } else {
    //empty online beamspot collection: FED data was empty
    //the error should probably have been send at unpacker level
    fallBackToDB = true;
  }

  if (fallBackToDB) {
    edm::ESHandle<BeamSpotObjects> beamhandle;
    iSetup.get<BeamSpotObjectsRcd>().get(beamhandle);
    const BeamSpotObjects* spotDB = beamhandle.product();

    // translate from BeamSpotObjects to reco::BeamSpot
    reco::BeamSpot::Point apoint(spotDB->GetX(), spotDB->GetY(), spotDB->GetZ());

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

    // this assume beam width same in x and y
    aSpot = reco::BeamSpot(
        apoint, spotDB->GetSigmaZ(), spotDB->Getdxdz(), spotDB->Getdydz(), spotDB->GetBeamWidthX(), matrix);
    aSpot.setBeamWidthY(spotDB->GetBeamWidthY());
    aSpot.setEmittanceX(spotDB->GetEmittanceX());
    aSpot.setEmittanceY(spotDB->GetEmittanceY());
    aSpot.setbetaStar(spotDB->GetBetaStar());
    aSpot.setType(reco::BeamSpot::Tracker);
  }

  *result = aSpot;

  iEvent.put(std::move(result));
}

Member Data Documentation

const bool BeamSpotOnlineProducer::changeFrame_

private

Definition at line 36 of file BeamSpotOnlineProducer.h.

Referenced by produce().

const edm::EDGetTokenT<L1GlobalTriggerEvmReadoutRecord> BeamSpotOnlineProducer::l1GtEvmReadoutRecordToken_

private

Definition at line 40 of file BeamSpotOnlineProducer.h.

Referenced by produce().

const edm::EDGetTokenT<BeamSpotOnlineCollection> BeamSpotOnlineProducer::scalerToken_

private

Definition at line 39 of file BeamSpotOnlineProducer.h.

Referenced by produce().

const unsigned int BeamSpotOnlineProducer::theBeamShoutMode

private

Definition at line 42 of file BeamSpotOnlineProducer.h.

Referenced by produce().

double BeamSpotOnlineProducer::theMaxR2

private

Definition at line 38 of file BeamSpotOnlineProducer.h.

Referenced by BeamSpotOnlineProducer(), and produce().

const double BeamSpotOnlineProducer::theMaxZ

private

Definition at line 37 of file BeamSpotOnlineProducer.h.

Referenced by produce().

const double BeamSpotOnlineProducer::theSetSigmaZ

private

Definition at line 37 of file BeamSpotOnlineProducer.h.

Referenced by produce().