TrackerSeedCleaner Class Reference

#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
bool	useDirection_Cleaner
bool	usePt_Cleaner

Detailed Description

Seeds Cleaner based on direction

Author

A. Grelli - Purdue University, Pavia University

Definition at line 34 of file TrackerSeedCleaner.h.

Member Typedef Documentation

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

Definition at line 38 of file TrackerSeedCleaner.h.

Constructor & Destructor Documentation

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

inline

constructor

Definition at line 40 of file TrackerSeedCleaner.h.

References edm::ConsumesCollector::consumes(), edm::ParameterSet::getParameter(), StGenEvtProducer_cfi::init, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                           : theProxyService(0),theEvent(0) {
                   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);
  }

virtual TrackerSeedCleaner::~TrackerSeedCleaner ( )

inlinevirtual

destructor

Definition at line 53 of file TrackerSeedCleaner.h.

References clean.

{}

Member Function Documentation

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

virtual

the cleaner

Definition at line 61 of file TrackerSeedCleaner.cc.

References MuonCkfTrajectoryBuilder_cfi::deltaPhi, TrackingRegion::direction(), PV3DBase< T, PVType, FrameType >::eta(), MuonErrorMatrixValues_cff::etaRange, RectangularEtaPhiTrackingRegion::etaRange(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::globalMomentum(), LogDebug, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), RectangularEtaPhiTrackingRegion::phiMargin(), ALCARECOTkAlBeamHalo_cff::ptMin, TrackingRegion::ptMin(), rpcPointValidation_cfi::recHit, mps_fire::result, TkTrackingRegionsMargin< T >::right(), SurveyInfoScenario_cff::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().

Referenced by TSGFromL2Muon::produce().

                                                                                                                       {


 // call the shared input cleaner
 if(cleanBySharedHits) theRedundantCleaner->define(seeds);

 theProxyService->eventSetup().get<TransientRecHitRecord>().get(builderName_,theTTRHBuilder);

 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().second-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.size()!=0 && (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 43 of file TrackerSeedCleaner.cc.

References gsfElectronCkfTrackCandidateMaker_cff::RedundantSeedCleaner.

Referenced by TSGFromL2Muon::beginRun().

                                                            {

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

void TrackerSeedCleaner::setEvent ( const edm::Event &  event )

virtual

setEvent

Definition at line 53 of file TrackerSeedCleaner.cc.

Referenced by TSGFromL2Muon::produce().

{
 event.getByToken(beamspotToken_, bsHandle_);
}

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot> TrackerSeedCleaner::beamspotToken_

private

Definition at line 66 of file TrackerSeedCleaner.h.

edm::Handle<reco::BeamSpot> TrackerSeedCleaner::bsHandle_

private

Definition at line 65 of file TrackerSeedCleaner.h.

std::string TrackerSeedCleaner::builderName_

private

Definition at line 69 of file TrackerSeedCleaner.h.

bool TrackerSeedCleaner::cleanBySharedHits

private

Definition at line 71 of file TrackerSeedCleaner.h.

edm::InputTag TrackerSeedCleaner::theBeamSpotTag

private

Definition at line 64 of file TrackerSeedCleaner.h.

const edm::Event* TrackerSeedCleaner::theEvent

private

Definition at line 62 of file TrackerSeedCleaner.h.

const MuonServiceProxy* TrackerSeedCleaner::theProxyService

private

Definition at line 61 of file TrackerSeedCleaner.h.

RedundantSeedCleaner* TrackerSeedCleaner::theRedundantCleaner

private

Definition at line 67 of file TrackerSeedCleaner.h.

edm::ESHandle<TransientTrackingRecHitBuilder> TrackerSeedCleaner::theTTRHBuilder

private

Definition at line 70 of file TrackerSeedCleaner.h.

bool TrackerSeedCleaner::useDirection_Cleaner

private

Definition at line 71 of file TrackerSeedCleaner.h.

bool TrackerSeedCleaner::usePt_Cleaner

private

Definition at line 71 of file TrackerSeedCleaner.h.