BeamSpotProducer.cc

Go to the documentation of this file.

#include "RecoVertex/BeamSpotProducer/interface/BeamSpotProducer.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"

#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/IOVSyncValue.h"
#include "CondFormats/DataRecord/interface/BeamSpotObjectsRcd.h"
#include "CondFormats/BeamSpotObjects/interface/BeamSpotObjects.h"

#include "DataFormats/Math/interface/Error.h"
#include "DataFormats/Math/interface/Point3D.h"

//
// constructors and destructor
//
BeamSpotProducer::BeamSpotProducer(const edm::ParameterSet& iConf) {
    
    edm::LogInfo("RecoVertex/BeamSpotProducer") 
        << "Initializing Beam Spot producer " << "\n";
  
    //fVerbose=conf.getUntrackedParameter<bool>("verbose", false);
    
    produces<reco::BeamSpot>();

}


BeamSpotProducer::~BeamSpotProducer() {}

//
// member functions
//

// ------------ method called to produce the data  ------------
void
BeamSpotProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
    
    using namespace edm;

    auto result = std::make_unique<reco::BeamSpot>();

    reco::BeamSpot aSpot;

    
    //typedef math::XYZPoint Point;
    //enum { dimension = 7 };
    //typedef math::Error<dimension>::type CovarianceMatrix;

    
    //try {
    edm::LogInfo("RecoVertex/BeamSpotProducer") 
      << "Reconstructing event number: " << iEvent.id() << "\n";

    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() );
        
    //}
    //
    //catch (std::exception & err) {
    //  edm::LogInfo("RecoVertex/BeamSpotProducer") 
    //      << "Exception during event number: " << iEvent.id() 
    //      << "\n" << err.what() << "\n";
    //}

    *result = aSpot;
    
    iEvent.put(std::move(result));
    
}

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