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 Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	BeamSpotOnlineProducer (const edm::ParameterSet &iConf)
	constructor More...
virtual void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	produce a beam spot class More...
	~BeamSpotOnlineProducer ()
	destructor More...
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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

Additional Inherited Members
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 24 of file BeamSpotOnlineProducer.h.

Member Typedef Documentation

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

Definition at line 27 of file BeamSpotOnlineProducer.h.

Constructor & Destructor Documentation

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

explicit

constructor

Definition at line 17 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

(

)

destructor

Definition at line 33 of file BeamSpotOnlineProducer.cc.

{}

Member Function Documentation

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

overridevirtual

produce a beam spot class

Implements edm::stream::EDProducerBase.

Definition at line 36 of file BeamSpotOnlineProducer.cc.

References align::BeamSpot, changeFrame_, 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(), i, j, l1GtEvmReadoutRecordToken_, reco::BeamSpot::LHC, makeMuonMisalignmentScenario::matrix, edm::ESHandle< class >::product(), edm::Event::put(), diffTwoXMLs::r2, query::result, scalerToken_, reco::BeamSpot::setBeamWidthY(), reco::BeamSpot::setbetaStar(), reco::BeamSpot::setEmittanceX(), reco::BeamSpot::setEmittanceY(), reco::BeamSpot::setType(), BeamSpotOnline::sigma_z(), 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
  std::auto_ptr<reco::BeamSpot> result(new reco::BeamSpot);
  
  reco::BeamSpot aSpot;

  bool fallBackToDB=false;
  if (handleScaler->size()!=0){
    // 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(result);

}

Member Data Documentation

const bool BeamSpotOnlineProducer::changeFrame_

private

Definition at line 39 of file BeamSpotOnlineProducer.h.

Referenced by produce().

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

private

Definition at line 43 of file BeamSpotOnlineProducer.h.

Referenced by produce().

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

private

Definition at line 42 of file BeamSpotOnlineProducer.h.

Referenced by produce().

const unsigned int BeamSpotOnlineProducer::theBeamShoutMode

private

Definition at line 45 of file BeamSpotOnlineProducer.h.

Referenced by produce().

double BeamSpotOnlineProducer::theMaxR2

private

Definition at line 41 of file BeamSpotOnlineProducer.h.

Referenced by BeamSpotOnlineProducer(), and produce().

const double BeamSpotOnlineProducer::theMaxZ

private

Definition at line 40 of file BeamSpotOnlineProducer.h.

Referenced by produce().

const double BeamSpotOnlineProducer::theSetSigmaZ

private

Definition at line 40 of file BeamSpotOnlineProducer.h.

Referenced by produce().