CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
AlcaBeamSpotFromDB.cc
Go to the documentation of this file.
1 
12 // C++ standard
13 #include <string>
14 // CMS
17 
19 
25 
27  : beamSpotToken_(esConsumes<BeamSpotObjects, BeamSpotObjectsRcd, edm::Transition::EndLuminosityBlock>()) {
28  produces<reco::BeamSpot, edm::Transition::EndLuminosityBlock>("alcaBeamSpot");
29 }
30 
32 
33 //--------------------------------------------------------------------------------------------------
35 
36 //--------------------------------------------------------------------------------------------------
38  // read DB object
39  const BeamSpotObjects *spotDB = &iSetup.getData(beamSpotToken_);
40 
41  // translate from BeamSpotObjects to reco::BeamSpot
42  reco::BeamSpot::Point apoint(spotDB->GetX(), spotDB->GetY(), spotDB->GetZ());
43 
45  for (int i = 0; i < 7; ++i) {
46  for (int j = 0; j < 7; ++j) {
47  matrix(i, j) = spotDB->GetCovariance(i, j);
48  }
49  }
50 
51  reco::BeamSpot aSpot;
52  // this assume beam width same in x and y
53  aSpot = reco::BeamSpot(
54  apoint, spotDB->GetSigmaZ(), spotDB->Getdxdz(), spotDB->Getdydz(), spotDB->GetBeamWidthX(), matrix);
55  aSpot.setBeamWidthY(spotDB->GetBeamWidthY());
56  aSpot.setEmittanceX(spotDB->GetEmittanceX());
57  aSpot.setEmittanceY(spotDB->GetEmittanceY());
58  aSpot.setbetaStar(spotDB->GetBetaStar());
59 
60  if (spotDB->GetBeamType() == 2) {
62  } else {
64  }
65 
66  auto result = std::make_unique<reco::BeamSpot>();
67  *result = aSpot;
68  lumiSeg.put(std::move(result), std::string("alcaBeamSpot"));
69 
70  // std::cout << " for runs: " << iEvent.id().run() << " - " <<
71  // iEvent.id().run() << std::endl;
72  std::cout << aSpot << std::endl;
73 }
74 
76 
78 
79 // define this as a plug-in
double Getdydz() const
get dydz slope, crossing angle in YZ
math::Error< dimension >::type CovarianceMatrix
Definition: BeamSpot.h:29
const edm::ESGetToken< BeamSpotObjects, BeamSpotObjectsRcd > beamSpotToken_
double GetY() const
get Y beam position
void setbetaStar(double v)
Definition: BeamSpot.h:138
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
double GetSigmaZ() const
get sigma Z, RMS bunch length
double GetBeamWidthX() const
get average transverse beam width
math::XYZPoint Point
point in the space
Definition: BeamSpot.h:27
double GetBeamWidthY() const
get average transverse beam width
double GetEmittanceX() const
get emittance
tuple result
Definition: mps_fire.py:311
void setType(BeamType type)
set beam type
Definition: BeamSpot.h:124
bool getData(T &iHolder) const
Definition: EventSetup.h:128
int iEvent
Definition: GenABIO.cc:224
void setBeamWidthY(double v)
Definition: BeamSpot.h:105
void setEmittanceY(double v)
Definition: BeamSpot.h:137
int GetBeamType() const
get beam type
void put(std::unique_ptr< PROD > product)
Put a new product.
def move
Definition: eostools.py:511
void endLuminosityBlockProduce(edm::LuminosityBlock &lumiSeg, const edm::EventSetup &iSetup) final
Transition
Definition: Transition.h:12
double GetZ() const
get Z beam position
AlcaBeamSpotFromDB(const edm::ParameterSet &)
~AlcaBeamSpotFromDB() override
double Getdxdz() const
get dxdz slope, crossing angle in XZ
void produce(edm::Event &iEvent, const edm::EventSetup &iSetup) final
double GetX() const
get X beam position
double GetBetaStar() const
get beta star
double GetCovariance(int i, int j) const
get i,j element of the full covariance matrix 7x7
tuple cout
Definition: gather_cfg.py:144
double GetEmittanceY() const
get emittance
void setEmittanceX(double v)
Definition: BeamSpot.h:136
ESGetTokenH3DDVariant esConsumes(std::string const &Reccord, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283