HLTMuonPointingFilter Class Reference

#include <HLTMuonPointingFilter.h>

Inheritance diagram for HLTMuonPointingFilter:

Public Member Functions
virtual bool	filter (edm::Event &, edm::EventSetup const &) override
	HLTMuonPointingFilter (const edm::ParameterSet &)
	Constructor. More...
	~HLTMuonPointingFilter ()
	Destructor. More...
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
ModuleDescription const &	moduleDescription () const
virtual	~EDFilter ()
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 (const std::string &iProcessName, std::vector< const char * > &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 ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)

Private Attributes
unsigned long long	m_cacheRecordId
Cylinder::CylinderPointer	theCyl
const double	theMaxZ
Plane::PlanePointer	theNegPlane
Plane::PlanePointer	thePosPlane
Propagator *	thePropagator
const std::string	thePropagatorName
const double	theRadius
const edm::EDGetTokenT < reco::TrackCollection >	theSTAMuonToken

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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

EDFilter to select muons that points to a cylinder of configurable radius and lenght.

Author

Stefano Lacaprara - INFN Legnaro stefa.nosp@m.no.l.nosp@m.acapr.nosp@m.ara@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t

Definition at line 32 of file HLTMuonPointingFilter.h.

Constructor & Destructor Documentation

HLTMuonPointingFilter::HLTMuonPointingFilter

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 37 of file HLTMuonPointingFilter.cc.

References Cylinder::build(), newFWLiteAna::build, edm::ParameterSet::getParameter(), LogDebug, theCyl, theMaxZ, theNegPlane, thePosPlane, and theRadius.

                                                                        :
  theSTAMuonToken(  consumes<reco::TrackCollection>( pset.getParameter<edm::InputTag>("SALabel") ) ),      // token to read the muons
  thePropagatorName(pset.getParameter<std::string>("PropagatorName") ),
  theRadius(        pset.getParameter<double>("radius") ),              // cyl's radius (cm)
  theMaxZ(          pset.getParameter<double>("maxZ") ),                // cyl's half lenght (cm)
  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);

  LogDebug("HLTMuonPointing") << " SALabel : " << pset.getParameter<edm::InputTag>("SALabel")
    << " Radius : " << theRadius
    << " Half lenght : " << theMaxZ;
}

HLTMuonPointingFilter::~HLTMuonPointingFilter

(

)

Destructor.

Definition at line 62 of file HLTMuonPointingFilter.cc.

                                              {
}

Member Function Documentation

void HLTMuonPointingFilter::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 137 of file HLTMuonPointingFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLT_25ns14e33_v1_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                        {
  edm::ParameterSetDescription desc;

  desc.add<edm::InputTag>("SALabel", edm::InputTag("hltCosmicMuonBarrelOnly"));
  desc.add<std::string>("PropagatorName", "SteppingHelixPropagatorAny");
  desc.add<double>("radius", 90.0);
  desc.add<double>("maxZ", 280.0);

  descriptions.add("hltMuonPointingFilter", desc);
}

bool HLTMuonPointingFilter::filter ( edm::Event &  event, edm::EventSetup const &  eventSetup  )

overridevirtual

Implements edm::EDFilter.

Definition at line 66 of file HLTMuonPointingFilter.cc.

References accept(), anyDirection, edm::eventsetup::EventSetupRecord::cacheIdentifier(), TrajectoryStateOnSurface::freeState(), edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), TrajectoryStateOnSurface::globalPosition(), reco::TransientTrack::innermostMeasurementState(), TrajectoryStateOnSurface::isValid(), LogDebug, m_cacheRecordId, or, PV3DBase< T, PVType, FrameType >::perp(), edm::ESHandle< class >::product(), Propagator::propagate(), Propagator::propagationDirection(), muonTagProbeFilters_cff::staTracks, theCyl, theMaxZ, theNegPlane, thePosPlane, thePropagator, thePropagatorName, theRadius, theSTAMuonToken, and PV3DBase< T, PVType, FrameType >::z().

                                                                                   {
  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();
  }

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

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

  // Get the RecTrack collection from the event
  Handle<reco::TrackCollection> staTracks;
  event.getByToken(theSTAMuonToken, staTracks);

  reco::TrackCollection::const_iterator staTrack;

  for (staTrack = staTracks->begin(); staTrack != staTracks->end(); ++staTrack){
    reco::TransientTrack track(*staTrack,&*theMGField,theTrackingGeometry);

    TrajectoryStateOnSurface innerTSOS = track.innermostMeasurementState();

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

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

    if ( tsosAtCyl.isValid() ) {
      LogDebug("HLTMuonPointing") << " extrap TSOS " << tsosAtCyl;
      if (fabs(tsosAtCyl.globalPosition().z())<theMaxZ ) {
        accept=true;
        return accept;
      }
      else {
        LogDebug("HLTMuonPointing") << " extrap TSOS z too big " << tsosAtCyl.globalPosition().z();
    TrajectoryStateOnSurface tsosAtPlane;
    if (tsosAtCyl.globalPosition().z()>0)
      tsosAtPlane=thePropagator->propagate(*innerTSOS.freeState(), *thePosPlane);
    else
      tsosAtPlane=thePropagator->propagate(*innerTSOS.freeState(), *theNegPlane);

    if (tsosAtPlane.isValid()){
      if (tsosAtPlane.globalPosition().perp()< theRadius){
        accept=true;
        return accept;
      }
    }
    else
      LogDebug("HLTMuonPointing") << " extrap to plane failed ";
      }
    } else {
      LogDebug("HLTMuonPointing") << " extrap to cyl failed ";
    }

  }

  return accept;
}

Member Data Documentation

unsigned long long HLTMuonPointingFilter::m_cacheRecordId

private

Definition at line 57 of file HLTMuonPointingFilter.h.

Referenced by filter().

Cylinder::CylinderPointer HLTMuonPointingFilter::theCyl

private

Definition at line 53 of file HLTMuonPointingFilter.h.

Referenced by filter(), and HLTMuonPointingFilter().

const double HLTMuonPointingFilter::theMaxZ

private

Definition at line 51 of file HLTMuonPointingFilter.h.

Referenced by filter(), and HLTMuonPointingFilter().

Plane::PlanePointer HLTMuonPointingFilter::theNegPlane

private

Definition at line 54 of file HLTMuonPointingFilter.h.

Referenced by filter(), and HLTMuonPointingFilter().

Plane::PlanePointer HLTMuonPointingFilter::thePosPlane

private

Definition at line 54 of file HLTMuonPointingFilter.h.

Referenced by filter(), and HLTMuonPointingFilter().

Propagator* HLTMuonPointingFilter::thePropagator

mutableprivate

Definition at line 56 of file HLTMuonPointingFilter.h.

Referenced by filter().

const std::string HLTMuonPointingFilter::thePropagatorName

private

Definition at line 49 of file HLTMuonPointingFilter.h.

Referenced by filter().

const double HLTMuonPointingFilter::theRadius

private

Definition at line 50 of file HLTMuonPointingFilter.h.

Referenced by filter(), and HLTMuonPointingFilter().

const edm::EDGetTokenT<reco::TrackCollection> HLTMuonPointingFilter::theSTAMuonToken

private

Definition at line 48 of file HLTMuonPointingFilter.h.

Referenced by filter().