#include <TrackerSeedCleaner.h>

Public Types
typedef std::vector < TrajectorySeed >	tkSeeds

Public Member Functions
virtual void	clean (const reco::TrackRef &, const RectangularEtaPhiTrackingRegion &region, tkSeeds &)
	the cleaner More...

virtual void	init (const MuonServiceProxy *service)
	intizialization More...

virtual void	setEvent (const edm::Event &)
	setEvent More...

	TrackerSeedCleaner (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)
	constructor More...

virtual	~TrackerSeedCleaner ()
	destructor More...

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamspotToken_

edm::Handle< reco::BeamSpot >	bsHandle_

std::string	builderName_

bool	cleanBySharedHits

edm::InputTag	theBeamSpotTag

const edm::Event *	theEvent

const MuonServiceProxy *	theProxyService

RedundantSeedCleaner *	theRedundantCleaner

edm::ESHandle < TransientTrackingRecHitBuilder >	theTTRHBuilder

edm::ESGetToken < TransientTrackingRecHitBuilder, TransientRecHitRecord >	theTTRHBuilderToken

bool	useDirection_Cleaner

bool	usePt_Cleaner

Detailed Description

Seeds Cleaner based on direction

Author: A. Grelli - Purdue University, Pavia University

Definition at line 39 of file TrackerSeedCleaner.h.

Member Typedef Documentation

typedef std::vector<TrajectorySeed> TrackerSeedCleaner::tkSeeds

Definition at line 41 of file TrackerSeedCleaner.h.

Constructor & Destructor Documentation

TrackerSeedCleaner::TrackerSeedCleaner	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &	iC
	)

inline

constructor

Definition at line 43 of file TrackerSeedCleaner.h.

References beamspotToken_, builderName_, cleanBySharedHits, edm::ConsumesCollector::consumes(), edm::ConsumesCollector::esConsumes(), edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theBeamSpotTag, theTTRHBuilderToken, useDirection_Cleaner, and usePt_Cleaner.

       : theProxyService(nullptr), theEvent(nullptr) {
     builderName_ = pset.getParameter<std::string>("TTRHBuilder");
     theBeamSpotTag = pset.getParameter<edm::InputTag>("beamSpot");
     useDirection_Cleaner = pset.getParameter<bool>("directionCleaner");
     usePt_Cleaner = pset.getParameter<bool>("ptCleaner");
     cleanBySharedHits = pset.getParameter<bool>("cleanerFromSharedHits");
     beamspotToken_ = iC.consumes<reco::BeamSpot>(theBeamSpotTag);
     theTTRHBuilderToken = iC.esConsumes(edm::ESInputTag("", builderName_));
   }

virtual TrackerSeedCleaner::~TrackerSeedCleaner ( )

inlinevirtual

destructor

Definition at line 58 of file TrackerSeedCleaner.h.

58 {}

Member Function Documentation

void TrackerSeedCleaner::clean	(	const reco::TrackRef &	muR,
		const RectangularEtaPhiTrackingRegion &	region,
		tkSeeds &	seeds
	)

virtual

the cleaner

Definition at line 54 of file TrackerSeedCleaner.cc.

References cms::cuda::bs, srCondWrite_cfg::deltaPhi, TrackingRegion::direction(), PV3DBase< T, PVType, FrameType >::eta(), generateTowerEtThresholdLUT::etaRange, RectangularEtaPhiTrackingRegion::etaRange(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::globalMomentum(), LogDebug, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), RectangularEtaPhiTrackingRegion::phiMargin(), TrackingRegion::ptMin(), ptMin, mps_fire::result, TkTrackingRegionsMargin< T >::right(), fileCollector::seed, trajectoryStateTransform::transientState(), PV3DBase< T, PVType, FrameType >::x(), reco::BeamSpot::x0(), PV3DBase< T, PVType, FrameType >::y(), reco::BeamSpot::y0(), PV3DBase< T, PVType, FrameType >::z(), and reco::BeamSpot::z0().

                                                {
   // call the shared input cleaner
   if (cleanBySharedHits)
     theRedundantCleaner->define(seeds);
 
   theTTRHBuilder = theProxyService->eventSetup().getHandle(theTTRHBuilderToken);
 
   LogDebug("TrackerSeedCleaner") << seeds.size() << " trajectory seeds to the events before cleaning" << endl;
 
   //check the validity otherwise vertexing
   const reco::BeamSpot& bs = *bsHandle_;
   /*reco track and seeds as arguments. Seeds eta and phi are checked and 
    based on deviation from L2 eta and phi seed is accepted or not*/
 
   std::vector<TrajectorySeed> result;
 
   TSCBLBuilderNoMaterial tscblBuilder;
   // PerigeeConversions tspConverter;
   for (TrajectorySeedCollection::iterator seed = seeds.begin(); seed < seeds.end(); ++seed) {
     if (seed->nHits() < 2)
       continue;
     //get parameters and errors from the seed state
     TransientTrackingRecHit::RecHitPointer recHit = theTTRHBuilder->build(&*(seed->recHits().end() - 1));
     TrajectoryStateOnSurface state = trajectoryStateTransform::transientState(
         seed->startingState(), recHit->surface(), theProxyService->magneticField().product());
 
     TrajectoryStateClosestToBeamLine tsAtClosestApproachSeed =
         tscblBuilder(*state.freeState(), bs);  //as in TrackProducerAlgorithms
     if (!tsAtClosestApproachSeed.isValid())
       continue;
     GlobalPoint vSeed1 = tsAtClosestApproachSeed.trackStateAtPCA().position();
     GlobalVector pSeed = tsAtClosestApproachSeed.trackStateAtPCA().momentum();
     GlobalPoint vSeed(vSeed1.x() - bs.x0(), vSeed1.y() - bs.y0(), vSeed1.z() - bs.z0());
 
     //eta,phi info from seeds
     double etaSeed = state.globalMomentum().eta();
     double phiSeed = pSeed.phi();
 
     //if the limits are too stringent rescale limits
     typedef PixelRecoRange<float> Range;
     typedef TkTrackingRegionsMargin<float> Margin;
 
     Range etaRange = region.etaRange();
     double etaLimit = (fabs(fabs(etaRange.max()) - fabs(etaRange.mean())) < 0.1)
                           ? 0.1
                           : fabs(fabs(etaRange.max()) - fabs(etaRange.mean()));
 
     Margin phiMargin = region.phiMargin();
     double phiLimit = (phiMargin.right() < 0.1) ? 0.1 : phiMargin.right();
 
     double ptSeed = pSeed.perp();
     double ptMin = (region.ptMin() > 3.5) ? 3.5 : region.ptMin();
     // Clean
     bool inEtaRange = etaSeed >= (etaRange.mean() - etaLimit) && etaSeed <= (etaRange.mean() + etaLimit);
     bool inPhiRange = (fabs(deltaPhi(phiSeed, double(region.direction().phi()))) < phiLimit);
     // pt cleaner
     bool inPtRange = ptSeed >= ptMin && ptSeed <= 2 * (muR->pt());
 
     // save efficiency don't clean triplets with pt cleaner
     if (seed->nHits() == 3)
       inPtRange = true;
 
     // use pt and angle cleaners
     if (inPtRange && usePt_Cleaner && !useDirection_Cleaner) {
       result.push_back(*seed);
       LogDebug("TrackerSeedCleaner") << " Keeping the seed : this seed passed pt selection";
     }
 
     // use only angle default option
     if (inEtaRange && inPhiRange && !usePt_Cleaner && useDirection_Cleaner) {
       result.push_back(*seed);
       LogDebug("TrackerSeedCleaner") << " Keeping the seed : this seed passed direction selection";
     }
 
     // use all the cleaners
     if (inEtaRange && inPhiRange && inPtRange && usePt_Cleaner && useDirection_Cleaner) {
       result.push_back(*seed);
       LogDebug("TrackerSeedCleaner") << " Keeping the seed : this seed passed direction and pt selection";
     }
 
     LogDebug("TrackerSeedCleaner") << " eta for current seed " << etaSeed << "\n"
                                    << " phi for current seed " << phiSeed << "\n"
                                    << " eta for L2 track  " << muR->eta() << "\n"
                                    << " phi for L2 track  " << muR->phi() << "\n";
   }
 
   //the new seeds collection
   if (!result.empty() && (useDirection_Cleaner || usePt_Cleaner))
     seeds.swap(result);
 
   LogDebug("TrackerSeedCleaner") << seeds.size() << " trajectory seeds to the events after cleaning" << endl;
 
   return;
 }

void TrackerSeedCleaner::init ( const MuonServiceProxy * service )

virtual

intizialization

Definition at line 40 of file TrackerSeedCleaner.cc.

References HLT_FULL_cff::RedundantSeedCleaner.

                                                              {
   theProxyService = service;
 
   theRedundantCleaner = new RedundantSeedCleaner();
 }