BeamSpotOnlineProducer Class Reference

Inheritance diagram for BeamSpotOnlineProducer:

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...
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
const edm::ESGetToken< BeamSpotObjects, BeamSpotObjectsRcd >	beamToken_
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

Additional Inherited Members
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

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 26 of file BeamSpotOnlineProducer.cc.

Constructor & Destructor Documentation

BeamSpotOnlineProducer()

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

explicit

constructor

Definition at line 47 of file BeamSpotOnlineProducer.cc.

    : 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"))),
      beamToken_(esConsumes<BeamSpotObjects, BeamSpotObjectsRcd>()),
      theBeamShoutMode(iconf.getUntrackedParameter<unsigned int>("beamMode", 11)) {
  theMaxR2 = iconf.getParameter<double>("maxRadius");
  theMaxR2 *= theMaxR2;
 
  produces<reco::BeamSpot>();
}

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

Member Function Documentation

produce()

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

override

produce a beam spot class

Definition at line 61 of file BeamSpotOnlineProducer.cc.

                                                                            {
  //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.getHandle(beamToken_);
    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));
}

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

Member Data Documentation

beamToken_

const edm::ESGetToken<BeamSpotObjects, BeamSpotObjectsRcd> BeamSpotOnlineProducer::beamToken_

private

Definition at line 40 of file BeamSpotOnlineProducer.cc.

Referenced by produce().

changeFrame_

const bool BeamSpotOnlineProducer::changeFrame_

private

Definition at line 35 of file BeamSpotOnlineProducer.cc.

Referenced by produce().

l1GtEvmReadoutRecordToken_

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

private

Definition at line 39 of file BeamSpotOnlineProducer.cc.

Referenced by produce().

scalerToken_

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

private

Definition at line 38 of file BeamSpotOnlineProducer.cc.

Referenced by produce().

theBeamShoutMode

const unsigned int BeamSpotOnlineProducer::theBeamShoutMode

private

Definition at line 42 of file BeamSpotOnlineProducer.cc.

Referenced by produce().

theMaxR2

double BeamSpotOnlineProducer::theMaxR2

private

Definition at line 37 of file BeamSpotOnlineProducer.cc.

Referenced by BeamSpotOnlineProducer(), and produce().

theMaxZ

const double BeamSpotOnlineProducer::theMaxZ

private

Definition at line 36 of file BeamSpotOnlineProducer.cc.

Referenced by produce().

theSetSigmaZ

const double BeamSpotOnlineProducer::theSetSigmaZ

private

Definition at line 36 of file BeamSpotOnlineProducer.cc.

Referenced by produce().