#include <AlignmentGlobalTrackSelector.h>

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

Public Member Functions
	AlignmentGlobalTrackSelector (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)
	constructor More...

Tracks	select (const Tracks &tracks, const edm::Event &iEvent, const edm::EventSetup &eSetup)
	select tracks More...

bool	useThisFilter ()
	returns if any of the Filters is used. More...

	~AlignmentGlobalTrackSelector ()
	destructor More...

Private Member Functions
Tracks	checkIsolation (const Tracks &cands, const edm::Event &iEvent) const
	returns only isolated tracks in [cands] More...

Tracks	checkJetCount (const Tracks &cands, const edm::Event &iEvent) const
	returns [tracks] if there are less than theMaxCount Jets with theMinJetPt and an empty set if not More...

Tracks	findMuons (const Tracks &tracks, const edm::Event &iEvent) const
	filter for Tracks that match the Track of a global Muon More...

Tracks	matchTracks (const Tracks &src, const Tracks &comp) const
	matches [src] with [comp] returns collection with matching Tracks coming from [src] More...

void	printTracks (const Tracks &col) const
	print Information on Track-Collection More...

Private Attributes
edm::ParameterSet	theConf
	private data members More...

bool	theGMFilterSwitch

bool	theIsoFilterSwitch

bool	theJetCountFilterSwitch

edm::EDGetTokenT< reco::CaloJetCollection >	theJetCountToken

edm::EDGetTokenT< reco::CaloJetCollection >	theJetIsoToken

int	theMaxJetCount

double	theMaxJetPt

double	theMaxTrackDeltaR

int	theMinGlobalMuonCount

int	theMinIsolatedCount

double	theMinJetDeltaR

double	theMinJetPt

edm::EDGetTokenT< reco::MuonCollection >	theMuonToken

Detailed Description

Definition at line 17 of file AlignmentGlobalTrackSelector.h.

Member Typedef Documentation

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

Definition at line 22 of file AlignmentGlobalTrackSelector.h.

Constructor & Destructor Documentation

AlignmentGlobalTrackSelector::AlignmentGlobalTrackSelector	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC
	)

constructor

Definition at line 29 of file AlignmentGlobalTrackSelector.cc.

References edm::ConsumesCollector::consumes(), edm::ParameterSet::getParameter(), LogDebug, theGMFilterSwitch, theIsoFilterSwitch, theJetCountFilterSwitch, theJetCountToken, theJetIsoToken, theMaxJetCount, theMaxJetPt, theMaxTrackDeltaR, theMinGlobalMuonCount, theMinIsolatedCount, theMinJetDeltaR, theMinJetPt, and theMuonToken.

                                                                                                                 :
   theGMFilterSwitch(cfg.getParameter<bool>("applyGlobalMuonFilter")),
   theIsoFilterSwitch(cfg.getParameter<bool>("applyIsolationtest")),
   theJetCountFilterSwitch(cfg.getParameter<bool>("applyJetCountFilter"))
 {
   if (theGMFilterSwitch || theIsoFilterSwitch || theJetCountFilterSwitch)
     LogDebug("Alignment") << "> applying global Trackfilter ...";
 
   if (theGMFilterSwitch) {
     edm::InputTag theMuonSource = cfg.getParameter<InputTag>("muonSource");
     theMuonToken = iC.consumes<reco::MuonCollection>(theMuonSource);
     theMaxTrackDeltaR =cfg.getParameter<double>("maxTrackDeltaR");
     theMinGlobalMuonCount = cfg.getParameter<int>("minGlobalMuonCount");
     LogDebug("Alignment") << ">  GlobalMuonFilter : source, maxTrackDeltaR, min. Count       : "
               << theMuonSource << " , "
               << theMaxTrackDeltaR << " , "
               << theMinIsolatedCount;
   }
   
   if (theIsoFilterSwitch) {
     edm::InputTag theJetIsoSource = cfg.getParameter<InputTag>("jetIsoSource");
     theJetIsoToken = iC.consumes<reco::CaloJetCollection>(theJetIsoSource);
     theMaxJetPt = cfg.getParameter<double>("maxJetPt");
     theMinJetDeltaR = cfg.getParameter<double>("minJetDeltaR");
     theMinIsolatedCount = cfg.getParameter<int>("minIsolatedCount");
     LogDebug("Alignment") << ">  Isolationtest    : source, maxJetPt, minJetDeltaR, min. Count: "
               << theJetIsoSource << " , "
               << theMaxJetPt << " ,"
               << theMinJetDeltaR << " ,"
               << theMinGlobalMuonCount;
   }
   
   if (theJetCountFilterSwitch) {
     edm::InputTag theJetCountSource = cfg.getParameter<InputTag>("jetCountSource");
     theJetCountToken = iC.consumes<reco::CaloJetCollection>(theJetCountSource);
     theMinJetPt = cfg.getParameter<double>("minJetPt");
     theMaxJetCount = cfg.getParameter<int>("maxJetCount");
     LogDebug("Alignment") << ">  JetCountFilter   : source, minJetPt, maxJetCount             : "
               << theJetCountSource << " , "
               << theMinJetPt << " ,"
               << theMaxJetCount;  
   }
 }

AlignmentGlobalTrackSelector::~AlignmentGlobalTrackSelector ( )

destructor

Definition at line 74 of file AlignmentGlobalTrackSelector.cc.

75 {}

Member Function Documentation

AlignmentGlobalTrackSelector::Tracks AlignmentGlobalTrackSelector::checkIsolation	(	const Tracks &	cands,
		const edm::Event &	iEvent
	)		const

private

returns only isolated tracks in [cands]

Definition at line 136 of file AlignmentGlobalTrackSelector.cc.

References deltaR(), edm::Event::getByToken(), GlobalTrackerMuonAlignment_cfi::isolated, edm::HandleBase::isValid(), fwrapper::jets, mps_fire::result, theJetIsoToken, theMaxJetPt, theMinIsolatedCount, and theMinJetDeltaR.

Referenced by select().

 {
   Tracks result; result.clear();
 
   Handle<reco::CaloJetCollection> jets;
   iEvent.getByToken(theJetIsoToken, jets);
 
   if (jets.isValid()) {
     for (Tracks::const_iterator it = cands.begin();
      it!=cands.end();
      ++it) {
       bool isolated = true;
       for (reco::CaloJetCollection::const_iterator itJet = jets->begin();
        itJet!=jets->end();
        ++itJet) 
     isolated &= !((*itJet).pt() > theMaxJetPt && deltaR(*(*it),(*itJet)) < theMinJetDeltaR);
       
       if (isolated)
     result.push_back(*it);
     }
     //    LogDebug("Alignment") << "D  Found "<<result.size()<<" isolated of "<< cands.size()<<" Tracks!";   
     
   } else
     LogError("Alignment")<< "@SUB=AlignmentGlobalTrackSelector::checkIsolation"
              << "> could not optain jetCollection!";
 
   if (static_cast<int>(result.size()) < theMinIsolatedCount) result.clear();
 
   return result;
 }

AlignmentGlobalTrackSelector::Tracks AlignmentGlobalTrackSelector::checkJetCount	(	const Tracks &	cands,
		const edm::Event &	iEvent
	)		const

private

returns [tracks] if there are less than theMaxCount Jets with theMinJetPt and an empty set if not

Definition at line 169 of file AlignmentGlobalTrackSelector.cc.

References edm::Event::getByToken(), edm::HandleBase::isValid(), fwrapper::jets, LogDebug, mps_fire::result, theJetCountToken, theMaxJetCount, theMinJetPt, and l1t::tracks.

Referenced by select().

 {
   Tracks result; result.clear();
 
   Handle<reco::CaloJetCollection> jets;
   iEvent.getByToken(theJetCountToken, jets);
 
   if (jets.isValid()) {
     int jetCount = 0;
     for (reco::CaloJetCollection::const_iterator itJet = jets->begin();
      itJet!=jets->end();
      ++itJet) {
       if ((*itJet).pt() > theMinJetPt)
     jetCount++;
     }
 
     if (jetCount <= theMaxJetCount)
       result = tracks;
 
     LogDebug("Alignment") << "> found " << jetCount << " Jets";
   } else
     LogError("Alignment") << "@SUB=AlignmentGlobalTrackSelector::checkJetCount"
               << ">  could not optain jetCollection!";
 
   return result;
 }

AlignmentGlobalTrackSelector::Tracks AlignmentGlobalTrackSelector::findMuons	(	const Tracks &	tracks,
		const edm::Event &	iEvent
	)		const

private

filter for Tracks that match the Track of a global Muon

Definition at line 101 of file AlignmentGlobalTrackSelector.cc.

References edm::Event::getByToken(), globalMuons_cfi::globalMuons, edm::HandleBase::isValid(), LogDebug, matchTracks(), extraflags_cff::muons, mps_fire::result, theMinGlobalMuonCount, and theMuonToken.

Referenced by select().

 {
   Tracks result;
   Tracks globalMuons;
 
   //fill globalMuons with muons
   Handle<reco::MuonCollection> muons;
   iEvent.getByToken(theMuonToken, muons);
 
   if (muons.isValid()) {
     for (reco::MuonCollection::const_iterator itMuon = muons->begin(); 
      itMuon != muons->end();
      ++itMuon) {
       const reco::Track* muonTrack = (*itMuon).get<reco::TrackRef>().get();
       if (!muonTrack) {
     LogDebug("Alignment") << "@SUB=AlignmentGlobalTrackSelector::findMuons"
                   << "Found muon without track: Standalone Muon!";
       } else {
     globalMuons.push_back(muonTrack);
       }
     }
   } else {
     LogError("Alignment") << "@SUB=AlignmentGlobalTrackSelector::findMuons"
               <<">  could not optain mounCollection!";
   }
 
   result = this->matchTracks(tracks, globalMuons);
   
   if (static_cast<int>(result.size()) < theMinGlobalMuonCount) result.clear();
 
   return result;
 }

AlignmentGlobalTrackSelector::Tracks AlignmentGlobalTrackSelector::matchTracks	(	const Tracks &	src,
		const Tracks &	comp
	)		const

private

matches [src] with [comp] returns collection with matching Tracks coming from [src]

Definition at line 200 of file AlignmentGlobalTrackSelector.cc.

References deltaR(), mps_fire::i, match(), min(), mps_fire::result, and theMaxTrackDeltaR.

Referenced by findMuons().

 {
   Tracks result;
   for (Tracks::const_iterator itComp = comp.begin();
        itComp!=comp.end();
        ++itComp) {
     int match = -1;
     double min = theMaxTrackDeltaR;
     for (unsigned int i=0;i<src.size();i++) {
       // LogDebug("Alignment") << "> Trackmatch dist: "<<deltaR(src.at(i),*itComp);
       if(min > deltaR(*(src.at(i)),*(*itComp))){
     min = deltaR(*(src.at(i)),*(*itComp));
     match = static_cast<int>(i);
       }
     }
     if (match > -1)
       result.push_back(src.at(match)); 
   }
   return result;
 }

void AlignmentGlobalTrackSelector::printTracks ( const Tracks & col ) const

private

print Information on Track-Collection

Definition at line 222 of file AlignmentGlobalTrackSelector.cc.

References KineDebug3::count(), and LogDebug.

 {
   int count = 0;
   LogDebug("Alignment") << ">......................................";
   for (Tracks::const_iterator it = col.begin();
        it < col.end();
        ++it,++count) {
     LogDebug("Alignment") 
       << ">  Track No. " << count << ": p = (" 
       << (*it)->px() << "," 
       << (*it)->py() << ","
       << (*it)->pz() << ")\n"
       << ">                        pT = "
       << (*it)->pt() << " eta = "
       << (*it)->eta() << " charge = "
       << (*it)->charge();    
   }
   LogDebug("Alignment") << ">......................................";
 }

AlignmentGlobalTrackSelector::Tracks AlignmentGlobalTrackSelector::select	(	const Tracks &	tracks,
		const edm::Event &	iEvent,
		const edm::EventSetup &	eSetup
	)

select tracks

Definition at line 85 of file AlignmentGlobalTrackSelector.cc.

References checkIsolation(), checkJetCount(), findMuons(), LogDebug, mps_fire::result, theGMFilterSwitch, theIsoFilterSwitch, theJetCountFilterSwitch, and l1t::tracks.

Referenced by Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Views.LineDecayView.LineDecayContainer::mousePressEvent(), Vispa.Gui.PortConnection.PointToPointConnection::mousePressEvent(), Vispa.Gui.VispaWidget.VispaWidget::mousePressEvent(), Vispa.Views.AbstractView.AbstractView::restoreSelection(), and TrackConfigSelector::select().

 {
   Tracks result = tracks;
 
   if (theGMFilterSwitch) result = findMuons(result, iEvent);
   if (theIsoFilterSwitch) result = checkIsolation(result, iEvent);
   if (theJetCountFilterSwitch) result = checkJetCount(result, iEvent);
   LogDebug("Alignment") << ">  Global: tracks all, kept: " << tracks.size() << ", " << result.size();
 //  LogDebug("Alignment")<<">  o kept:";
 //  printTracks(result);
   
   return result;
 }

bool AlignmentGlobalTrackSelector::useThisFilter ( )

returns if any of the Filters is used.

Definition at line 78 of file AlignmentGlobalTrackSelector.cc.

References theGMFilterSwitch, theIsoFilterSwitch, and theJetCountFilterSwitch.

Referenced by TrackConfigSelector::TrackConfigSelector().

 {
   return theGMFilterSwitch || theIsoFilterSwitch|| theJetCountFilterSwitch;
 }

Member Data Documentation

edm::ParameterSet AlignmentGlobalTrackSelector::theConf

private

private data members

Definition at line 45 of file AlignmentGlobalTrackSelector.h.

bool AlignmentGlobalTrackSelector::theGMFilterSwitch

private

Definition at line 48 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), select(), and useThisFilter().

bool AlignmentGlobalTrackSelector::theIsoFilterSwitch

private

Definition at line 49 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), select(), and useThisFilter().

bool AlignmentGlobalTrackSelector::theJetCountFilterSwitch

private

Definition at line 50 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), select(), and useThisFilter().

edm::EDGetTokenT<reco::CaloJetCollection> AlignmentGlobalTrackSelector::theJetCountToken

private

Definition at line 64 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and checkJetCount().

edm::EDGetTokenT<reco::CaloJetCollection> AlignmentGlobalTrackSelector::theJetIsoToken

private

Definition at line 58 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and checkIsolation().

int AlignmentGlobalTrackSelector::theMaxJetCount

private

Definition at line 66 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and checkJetCount().

double AlignmentGlobalTrackSelector::theMaxJetPt

private

Definition at line 59 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and checkIsolation().

double AlignmentGlobalTrackSelector::theMaxTrackDeltaR

private

Definition at line 54 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and matchTracks().

int AlignmentGlobalTrackSelector::theMinGlobalMuonCount

private

Definition at line 55 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and findMuons().

int AlignmentGlobalTrackSelector::theMinIsolatedCount

private

Definition at line 61 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and checkIsolation().

double AlignmentGlobalTrackSelector::theMinJetDeltaR

private

Definition at line 60 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and checkIsolation().

double AlignmentGlobalTrackSelector::theMinJetPt

private

Definition at line 65 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and checkJetCount().

edm::EDGetTokenT<reco::MuonCollection> AlignmentGlobalTrackSelector::theMuonToken

private

Definition at line 53 of file AlignmentGlobalTrackSelector.h.

Referenced by AlignmentGlobalTrackSelector(), and findMuons().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation