#include <AlignmentCSCTrackSelector.h>

Public Types
typedef std::vector< const reco::Track * >	Tracks

Public Member Functions
	AlignmentCSCTrackSelector (const edm::ParameterSet &cfg)
	constructor More...

Tracks	select (const Tracks &tracks, const edm::Event &evt) const
	select tracks More...

	~AlignmentCSCTrackSelector ()
	destructor More...

Private Attributes
int	m_maxHitsPerStation

int	m_minHitsDT

int	m_minHitsPerStation

edm::InputTag	m_src

int	m_stationA

int	m_stationB

Detailed Description

Definition at line 15 of file AlignmentCSCTrackSelector.h.

Member Typedef Documentation

◆ Tracks

typedef std::vector<const reco::Track*> AlignmentCSCTrackSelector::Tracks

Definition at line 17 of file AlignmentCSCTrackSelector.h.

Constructor & Destructor Documentation

◆ AlignmentCSCTrackSelector()

AlignmentCSCTrackSelector::AlignmentCSCTrackSelector ( const edm::ParameterSet & cfg )

constructor

Definition at line 15 of file AlignmentCSCTrackSelector.cc.

     : m_src(cfg.getParameter<edm::InputTag>("src")),
       m_stationA(cfg.getParameter<int>("stationA")),
       m_stationB(cfg.getParameter<int>("stationB")),
       m_minHitsDT(cfg.getParameter<int>("minHitsDT")),
       m_minHitsPerStation(cfg.getParameter<int>("minHitsPerStation")),
       m_maxHitsPerStation(cfg.getParameter<int>("maxHitsPerStation")) {}

◆ ~AlignmentCSCTrackSelector()

AlignmentCSCTrackSelector::~AlignmentCSCTrackSelector ( )

destructor

Definition at line 25 of file AlignmentCSCTrackSelector.cc.

25 {}

Member Function Documentation

◆ select()

AlignmentCSCTrackSelector::Tracks AlignmentCSCTrackSelector::select	(	const Tracks &	tracks,
		const edm::Event &	evt
	)		const

select tracks

Definition at line 29 of file AlignmentCSCTrackSelector.cc.

                                                                                                                  {
   Tracks result;
  
   for (auto const& track : tracks) {
     int hitsOnStationA = 0;
     int hitsOnStationB = 0;
  
     for (auto const& hit : track->recHits()) {
       DetId id = hit->geographicalId();
  
       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::DT) {
         if (m_stationA == 0)
           hitsOnStationA++;
         if (m_stationB == 0)
           hitsOnStationB++;
       } else if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::CSC) {
         CSCDetId cscid(id.rawId());
         int station = (cscid.endcap() == 1 ? 1 : -1) * cscid.station();
  
         if (station == m_stationA)
           hitsOnStationA++;
         if (station == m_stationB)
           hitsOnStationB++;
  
       }  // end if CSC
     }    // end loop over hits
  
     bool stationAokay;
     if (m_stationA == 0)
       stationAokay = (m_minHitsDT <= hitsOnStationA);
     else
       stationAokay = (m_minHitsPerStation <= hitsOnStationA && hitsOnStationA <= m_maxHitsPerStation);
  
     bool stationBokay;
     if (m_stationB == 0)
       stationBokay = (m_minHitsDT <= hitsOnStationB);
     else
       stationBokay = (m_minHitsPerStation <= hitsOnStationB && hitsOnStationB <= m_maxHitsPerStation);
  
     if (stationAokay && stationBokay) {
       result.push_back(track);
     }
   }  // end loop over tracks
  
   return result;
 }