Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Modernized HLTMuonPointingFilter #30987

Merged
merged 1 commit into from Aug 3, 2020
Merged

Modernized HLTMuonPointingFilter #30987

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
59 changes: 19 additions & 40 deletions HLTrigger/special/plugins/HLTMuonPointingFilter.cc
View file
Open in desktop
Expand Up @@ -19,13 +19,8 @@

#include "DataFormats/TrackReco/interface/Track.h"

#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
#include "MagneticField/Engine/interface/MagneticField.h"
#include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"
#include "TrackingTools/GeomPropagators/interface/Propagator.h"
#include "TrackingTools/Records/interface/TrackingComponentsRecord.h"

/* C++ Headers */
using namespace std;
Expand All @@ -38,24 +33,18 @@ HLTMuonPointingFilter::HLTMuonPointingFilter(const edm::ParameterSet& pset)
: theSTAMuonToken(
consumes<reco::TrackCollection>(pset.getParameter<edm::InputTag>("SALabel"))), // token to read the muons
thePropagatorName(pset.getParameter<std::string>("PropagatorName")),
thePropagatorToken(esConsumes<Propagator, TrackingComponentsRecord>(edm::ESInputTag("", thePropagatorName))),
theMGFieldToken(esConsumes<MagneticField, IdealMagneticFieldRecord>()),
theTrackingGeometryToken(esConsumes<GlobalTrackingGeometry, GlobalTrackingGeometryRecord>()),
theRadius(pset.getParameter<double>("radius")), // cyl's radius (cm)
theMaxZ(pset.getParameter<double>("maxZ")), // cyl's half lenght (cm)
thePixHits(pset.getParameter<unsigned int>("PixHits")), // pixel hits
theTkLayers(pset.getParameter<unsigned int>("TkLayers")), // tracker layers with measurements
theMuonHits(pset.getParameter<unsigned int>("MuonHits")), // muon hits
thePropagator(nullptr),
m_cacheRecordId(0) {
// Get a surface (here a cylinder of radius 1290mm) ECAL
Cylinder::PositionType pos0;
Cylinder::RotationType rot0;
theCyl = Cylinder::build(theRadius, pos0, rot0);

Plane::PositionType posPos(0, 0, theMaxZ);
Plane::PositionType posNeg(0, 0, -theMaxZ);

thePosPlane = Plane::build(posPos, rot0);
theNegPlane = Plane::build(posNeg, rot0);

// Get a surface (here a cylinder of radius 1290mm) ECAL
theCyl(Cylinder::build(theRadius, Cylinder::PositionType{}, Cylinder::RotationType{})),
thePosPlane(Plane::build(Plane::PositionType{0, 0, static_cast<float>(theMaxZ)}, Plane::RotationType{})),
theNegPlane(Plane::build(Plane::PositionType{0, 0, static_cast<float>(-theMaxZ)}, Plane::RotationType{})) {
LogDebug("HLTMuonPointing") << " SALabel : " << pset.getParameter<edm::InputTag>("SALabel")
<< " Radius : " << theRadius << " Half lenght : " << theMaxZ
<< " Min pixel hits : " << thePixHits << " Min tk layers measurements : " << theTkLayers
Expand All @@ -66,30 +55,20 @@ HLTMuonPointingFilter::HLTMuonPointingFilter(const edm::ParameterSet& pset)
HLTMuonPointingFilter::~HLTMuonPointingFilter() = default;

/* Operations */
bool HLTMuonPointingFilter::filter(edm::Event& event, const edm::EventSetup& eventSetup) {
bool HLTMuonPointingFilter::filter(edm::StreamID, edm::Event& event, const edm::EventSetup& eventSetup) const {
bool accept = false;

const TrackingComponentsRecord& tkRec = eventSetup.get<TrackingComponentsRecord>();
if (not thePropagator or tkRec.cacheIdentifier() != m_cacheRecordId) {
// delete the old propagator
delete thePropagator;

// get the new propagator from the EventSetup and clone it (for thread safety)
ESHandle<Propagator> propagatorHandle;
tkRec.get(thePropagatorName, propagatorHandle);
thePropagator = propagatorHandle.product()->clone();
if (thePropagator->propagationDirection() != anyDirection)
throw cms::Exception("Configuration")
<< "the propagator " << thePropagatorName
<< " should be configured with PropagationDirection = \"anyDirection\"" << std::endl;
m_cacheRecordId = tkRec.cacheIdentifier();
auto const& propagator = eventSetup.getData(thePropagatorToken);

if (propagator.propagationDirection() != anyDirection) {
throw cms::Exception("Configuration")
<< "the propagator " << thePropagatorName
<< " should be configured with PropagationDirection = \"anyDirection\"" << std::endl;
}

ESHandle<MagneticField> theMGField;
eventSetup.get<IdealMagneticFieldRecord>().get(theMGField);
ESHandle<MagneticField> theMGField = eventSetup.getHandle(theMGFieldToken);

ESHandle<GlobalTrackingGeometry> theTrackingGeometry;
eventSetup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);
ESHandle<GlobalTrackingGeometry> theTrackingGeometry = eventSetup.getHandle(theTrackingGeometryToken);

// Get the RecTrack collection from the event
Handle<reco::TrackCollection> staTracks;
Expand All @@ -110,7 +89,7 @@ bool HLTMuonPointingFilter::filter(edm::Event& event, const edm::EventSetup& eve

LogDebug("HLTMuonPointing") << " InnerTSOS " << innerTSOS;

TrajectoryStateOnSurface tsosAtCyl = thePropagator->propagate(*innerTSOS.freeState(), *theCyl);
TrajectoryStateOnSurface tsosAtCyl = propagator.propagate(*innerTSOS.freeState(), *theCyl);

if (tsosAtCyl.isValid()) {
LogDebug("HLTMuonPointing") << " extrap TSOS " << tsosAtCyl << " number of pixel hits " << pixelHits
Expand All @@ -132,9 +111,9 @@ bool HLTMuonPointingFilter::filter(edm::Event& event, const edm::EventSetup& eve
<< " number of muon hits " << nMuonHits;
TrajectoryStateOnSurface tsosAtPlane;
if (tsosAtCyl.globalPosition().z() > 0)
tsosAtPlane = thePropagator->propagate(*innerTSOS.freeState(), *thePosPlane);
tsosAtPlane = propagator.propagate(*innerTSOS.freeState(), *thePosPlane);
else
tsosAtPlane = thePropagator->propagate(*innerTSOS.freeState(), *theNegPlane);
tsosAtPlane = propagator.propagate(*innerTSOS.freeState(), *theNegPlane);

if (tsosAtPlane.isValid()) {
if (tsosAtPlane.globalPosition().perp() < theRadius) {
Expand Down
41 changes: 26 additions & 15 deletions HLTrigger/special/plugins/HLTMuonPointingFilter.h
View file
Open in desktop
Expand Up @@ -11,7 +11,7 @@
*/

/* Base Class Headers */
#include "FWCore/Framework/interface/EDFilter.h"
#include "FWCore/Framework/interface/global/EDFilter.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
Expand All @@ -22,14 +22,21 @@ class Propagator;
#include "DataFormats/GeometrySurface/interface/Plane.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"

#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
#include "MagneticField/Engine/interface/MagneticField.h"
#include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
#include "TrackingTools/Records/interface/TrackingComponentsRecord.h"

/* C++ Headers */
#include <string>
#include <memory>

/* ====================================================================== */

/* Class HLTMuonPointingFilter Interface */

class HLTMuonPointingFilter : public edm::EDFilter {
class HLTMuonPointingFilter : public edm::global::EDFilter<> {
public:
/// Constructor
HLTMuonPointingFilter(const edm::ParameterSet &);
Expand All @@ -38,23 +45,27 @@ class HLTMuonPointingFilter : public edm::EDFilter {
~HLTMuonPointingFilter() override;

/* Operations */
bool filter(edm::Event &, edm::EventSetup const &) override;
bool filter(edm::StreamID, edm::Event &, edm::EventSetup const &) const override;

static void fillDescriptions(edm::ConfigurationDescriptions &descriptions);

private:
const edm::EDGetTokenT<reco::TrackCollection> theSTAMuonToken;
const std::string thePropagatorName; // name of propagator to be used
const double theRadius; // radius of cylinder
const double theMaxZ; // half length of cylinder
const unsigned int thePixHits; // number of pixel hits
const unsigned int theTkLayers; // number of tracker layers with measurements
const unsigned int theMuonHits; // number of valid muon hits

Cylinder::CylinderPointer theCyl;
Plane::PlanePointer thePosPlane, theNegPlane;

mutable Propagator *thePropagator;
unsigned long long m_cacheRecordId;

const std::string
thePropagatorName; // name of propagator to be used const edm::ESGetToken<Propagator, TrackingComponentsRecord> thePropagatorToken;
const edm::ESGetToken<Propagator, TrackingComponentsRecord> thePropagatorToken;
const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> theMGFieldToken;
const edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> theTrackingGeometryToken;

const double theRadius; // radius of cylinder
const double theMaxZ; // half length of cylinder
const unsigned int thePixHits; // number of pixel hits
const unsigned int theTkLayers; // number of tracker layers with measurements
const unsigned int theMuonHits; // number of valid muon hits

const Cylinder::CylinderPointer theCyl;
const Plane::PlanePointer thePosPlane;
const Plane::PlanePointer theNegPlane;
};
#endif // Muon_HLTMuonPointingFilter_h