Alignment/TrackerAlignment/test/ApeAdder.cpp - master - githubmirror/cmssw

Alignment/TrackerAlignment/test/ApeAdder.cpp
108 строк4 KB
mmusich
get rid of CMS deprecation warnings in Alignment/TrackerAlignment
20 ноя 2021, 20:30
20 ноя 2021, 20:30e16b9af
// -*- C++ -*-
//
//
// Description: Module to add a custom APE to the alignment object
//
//
// Original Author:  Frederic Ronga
//         Created:  December 5, 2006
//

// system include files

// user include files
//#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

// Conditions database
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CondCore/DBOutputService/interface/PoolDBOutputService.h"

#include "Alignment/TrackerAlignment/interface/AlignableTracker.h"

#include "CondFormats/AlignmentRecord/interface/TrackerAlignmentRcd.h"
#include "CondFormats/Alignment/interface/AlignmentErrorsExtended.h"
#include "CondFormats/AlignmentRecord/interface/TrackerAlignmentErrorExtendedRcd.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/AlignmentPositionError.h"

#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"

// Class declaration
class ApeAdder : public edm::one::EDAnalyzer<> {
public:
  explicit ApeAdder(const edm::ParameterSet&);
  ~ApeAdder() override = default;

  void analyze(const edm::Event&, const edm::EventSetup&) override;

private:
  // methods
  void addApe(const align::Alignables& alignables);

private:
  // members
  const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> tTopoToken_;
  const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
  std::string theErrorRecordName;
  std::vector<double> theApe;  // Amount of APE to add (from config.)
};

ApeAdder::ApeAdder(const edm::ParameterSet& iConfig)
    : tTopoToken_(esConsumes()), tkGeomToken_(esConsumes()), theErrorRecordName("TrackerAlignmentErrorExtendedRcd") {
  // The APE to set to all GeomDets
  theApe = iConfig.getUntrackedParameter<std::vector<double> >("apeVector");
}

//__________________________________________________________________________________________________
void ApeAdder::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)

{
  //Retrieve tracker topology from geometry
  const TrackerTopology* const tTopo = &iSetup.getData(tTopoToken_);
  // Get geometry from ES
  const TrackerGeometry* trackerGeometry = &iSetup.getData(tkGeomToken_);

  // Create the alignable hierarchy
  AlignableTracker* theAlignableTracker = new AlignableTracker(&(*trackerGeometry), tTopo);

  // Now loop on alignable dets and add alignment error
  if (theAlignableTracker->barrelGeomDets().size())
    this->addApe(theAlignableTracker->barrelGeomDets());
  if (theAlignableTracker->pixelHalfBarrelGeomDets().size())
    this->addApe(theAlignableTracker->pixelHalfBarrelGeomDets());
  if (theAlignableTracker->endcapGeomDets().size())
    this->addApe(theAlignableTracker->endcapGeomDets());
  if (theAlignableTracker->TIDGeomDets().size())
    this->addApe(theAlignableTracker->TIDGeomDets());
  if (theAlignableTracker->pixelEndcapGeomDets().size())
    this->addApe(theAlignableTracker->pixelEndcapGeomDets());

  // Store to DB
  AlignmentErrorsExtended* alignmentErrors = theAlignableTracker->alignmentErrors();

  // Call service
  edm::Service<cond::service::PoolDBOutputService> poolDbService;
  if (!poolDbService.isAvailable())  // Die if not available
    throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";

  // Save to DB (assuming that the delete below will also delete the concerned payload )
  poolDbService->writeOneIOV<AlignmentErrorsExtended>(
      *alignmentErrors, poolDbService->beginOfTime(), theErrorRecordName);

  delete theAlignableTracker;
}

void ApeAdder::addApe(const align::Alignables& alignables) {
  AlignmentPositionError ape(theApe[0], theApe[1], theApe[2]);
  for (const auto& iDet : alignables) {
    iDet->setAlignmentPositionError(ape, true);  // true: propagate to components
  }
}

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