#include <CosmicTrackingParticleSelector.h>

Public Types
typedef TrackingParticleCollection	collection

typedef container::const_iterator	const_iterator

typedef std::vector< const TrackingParticle * >	container

Public Member Functions
const_iterator	begin () const

	CosmicTrackingParticleSelector ()

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

	CosmicTrackingParticleSelector (double ptMin, double minRapidity, double maxRapidity, double tip, double lip, int minHit, bool chargedOnly, const std::vector< int > &pdgId=std::vector< int >())

const_iterator	end () const

void	initEvent (edm::Handle< SimHitTPAssociationProducer::SimHitTPAssociationList > simHitsTPAssocToSet) const

bool	operator() (const TrackingParticleRef tpr, const reco::BeamSpot *bs, const edm::Event &iEvent, const edm::EventSetup &iSetup) const

void	select (const edm::Handle< collection > &c, const edm::Event &event, const edm::EventSetup &setup)

size_t	size () const

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

bool	chargedOnly_

edm::ESGetToken< GlobalTrackingGeometry, GlobalTrackingGeometryRecord >	globalTrackingGeomToken_

double	lip_

double	maxRapidity_

int	minHit_

double	minRapidity_

std::vector< int >	pdgId_

double	ptMin_

container	selected_

edm::Handle< SimHitTPAssociationProducer::SimHitTPAssociationList >	simHitsTPAssoc

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theMFToken_

double	tip_

Detailed Description

Definition at line 39 of file CosmicTrackingParticleSelector.h.

Member Typedef Documentation

◆ collection

typedef TrackingParticleCollection CosmicTrackingParticleSelector::collection

Definition at line 41 of file CosmicTrackingParticleSelector.h.

◆ const_iterator

typedef container::const_iterator CosmicTrackingParticleSelector::const_iterator

Definition at line 43 of file CosmicTrackingParticleSelector.h.

◆ container

typedef std::vector<const TrackingParticle*> CosmicTrackingParticleSelector::container

Definition at line 42 of file CosmicTrackingParticleSelector.h.

Constructor & Destructor Documentation

◆ CosmicTrackingParticleSelector() [1/3]

CosmicTrackingParticleSelector::CosmicTrackingParticleSelector ( )

inline

Definition at line 45 of file CosmicTrackingParticleSelector.h.

45 {}

◆ CosmicTrackingParticleSelector() [2/3]

CosmicTrackingParticleSelector::CosmicTrackingParticleSelector	(	double	ptMin,
		double	minRapidity,
		double	maxRapidity,
		double	tip,
		double	lip,
		int	minHit,
		bool	chargedOnly,
		const std::vector< int > &	pdgId = `std::vector<int>()`
	)

inline

Definition at line 47 of file CosmicTrackingParticleSelector.h.

       : ptMin_(ptMin),
         minRapidity_(minRapidity),
         maxRapidity_(maxRapidity),
         tip_(tip),
         lip_(lip),
         minHit_(minHit),
         chargedOnly_(chargedOnly),
         pdgId_(pdgId) {}

◆ CosmicTrackingParticleSelector() [3/3]

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

inline

Definition at line 64 of file CosmicTrackingParticleSelector.h.

       : ptMin_(cfg.getParameter<double>("ptMin")),
         minRapidity_(cfg.getParameter<double>("minRapidity")),
         maxRapidity_(cfg.getParameter<double>("maxRapidity")),
         tip_(cfg.getParameter<double>("tip")),
         lip_(cfg.getParameter<double>("lip")),
         minHit_(cfg.getParameter<int>("minHit")),
         chargedOnly_(cfg.getParameter<bool>("chargedOnly")),
         pdgId_(cfg.getParameter<std::vector<int> >("pdgId")),
         beamSpotToken_(iC.consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"))),
         globalTrackingGeomToken_(iC.esConsumes()),
         theMFToken_(iC.esConsumes()) {}

Member Function Documentation

◆ begin()

const_iterator CosmicTrackingParticleSelector::begin ( void ) const

inline

Definition at line 87 of file CosmicTrackingParticleSelector.h.

87 { return selected_.begin(); }

References selected_.

◆ end()

const_iterator CosmicTrackingParticleSelector::end ( void ) const

inline

Definition at line 88 of file CosmicTrackingParticleSelector.h.

88 { return selected_.end(); }

References selected_.

Referenced by Types.LuminosityBlockRange::cppID(), and Types.EventRange::cppID().

◆ initEvent()

void CosmicTrackingParticleSelector::initEvent ( edm::Handle< SimHitTPAssociationProducer::SimHitTPAssociationList > simHitsTPAssocToSet ) const

inline

Definition at line 90 of file CosmicTrackingParticleSelector.h.

                                                                                                           {
     simHitsTPAssoc = simHitsTPAssocToSet;
   }

References simHitsTPAssoc.

◆ operator()()

bool CosmicTrackingParticleSelector::operator()	(	const TrackingParticleRef	tpr,
		const reco::BeamSpot *	bs,
		const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

inline

Definition at line 95 of file CosmicTrackingParticleSelector.h.

                                                      {
     if (chargedOnly_ && tpr->charge() == 0)
       return false;  //select only if charge!=0
     //bool testId = false;
     //unsigned int idSize = pdgId_.size();
     //if (idSize==0) testId = true;
     //else for (unsigned int it=0;it!=idSize;++it){
     //if (tpr->pdgId()==pdgId_[it]) testId = true;
     //}
  
     edm::ESHandle<GlobalTrackingGeometry> theGeometry = iSetup.getHandle(globalTrackingGeomToken_);
  
     GlobalVector finalGV(0, 0, 0);
     GlobalPoint finalGP(0, 0, 0);
     GlobalVector momentum(0, 0, 0);  //At the PCA
     GlobalPoint vertex(0, 0, 0);     //At the PCA
     double radius(9999);
     bool found(false);
  
     int ii = 0;
     DetId::Detector det;
     int subdet;
  
     edm::LogVerbatim("CosmicTrackingParticleSelector")
         << "TOT Number of PSimHits = " << tpr->numberOfHits()
         << ", Number of Tracker PSimHits = " << tpr->numberOfTrackerHits() << "\n";
  
     if (simHitsTPAssoc.isValid() == 0) {
       edm::LogError("CosmicTrackingParticleSelector") << "Invalid handle!";
       return false;
     }
     std::pair<TrackingParticleRef, TrackPSimHitRef> clusterTPpairWithDummyTP(
         tpr, TrackPSimHitRef());  //SimHit is dummy: for simHitTPAssociationListGreater
                                   // sorting only the cluster is needed
     auto range = std::equal_range(simHitsTPAssoc->begin(),
                                   simHitsTPAssoc->end(),
                                   clusterTPpairWithDummyTP,
                                   SimHitTPAssociationProducer::simHitTPAssociationListGreater);
     for (auto ip = range.first; ip != range.second; ++ip) {
       TrackPSimHitRef it = ip->second;
       ++ii;
       const GeomDet* tmpDet = theGeometry->idToDet(DetId(it->detUnitId()));
       if (!tmpDet) {
         edm::LogVerbatim("CosmicTrackingParticleSelector")
             << "***WARNING:  PSimHit " << ii << ", no GeomDet for: " << it->detUnitId() << ". Skipping it.";
         continue;
       } else {
         det = DetId(it->detUnitId()).det();
         subdet = DetId(it->detUnitId()).subdetId();
       }
  
       LocalVector lv = it->momentumAtEntry();
       Local3DPoint lp = it->localPosition();
       GlobalVector gv = tmpDet->surface().toGlobal(lv);
       GlobalPoint gp = tmpDet->surface().toGlobal(lp);
       edm::LogVerbatim("CosmicTrackingParticleSelector")
           << "PSimHit " << ii << ", Detector = " << det << ", subdet = " << subdet << "\t Radius = " << gp.perp()
           << ", z = " << gp.z() << "\t     pt = " << gv.perp() << ", pz = " << gv.z();
       edm::LogVerbatim("CosmicTrackingParticleSelector")
           << "\t trackId = " << it->trackId() << ", particleType = " << it->particleType()
           << ", processType = " << it->processType();
  
       // discard hits related to low energy debris from the primary particle
       if (it->processType() != 0)
         continue;
  
       if (gp.perp() < radius) {
         found = true;
         radius = gp.perp();
         finalGV = gv;
         finalGP = gp;
       }
     }
     edm::LogVerbatim("CosmicTrackingParticleSelector")
         << "\n"
         << "FINAL State at InnerMost Hit:        Radius = " << finalGP.perp() << ", z = " << finalGP.z()
         << ", pt = " << finalGV.perp() << ", pz = " << finalGV.z();
  
     if (!found)
       return false;
     else {
       FreeTrajectoryState ftsAtProduction(finalGP, finalGV, TrackCharge(tpr->charge()), &iSetup.getData(theMFToken_));
       TSCBLBuilderNoMaterial tscblBuilder;
       //as in TrackProducerAlgorithm
       TrajectoryStateClosestToBeamLine tsAtClosestApproach = tscblBuilder(ftsAtProduction, *bs);
       if (!tsAtClosestApproach.isValid()) {
         edm::LogVerbatim("CosmicTrackingParticleSelector")
             << "*** WARNING in CosmicTrackingParticleSelector: tsAtClosestApproach is not valid."
             << "\n";
         return false;
       } else {
         momentum = tsAtClosestApproach.trackStateAtPCA().momentum();
         vertex = tsAtClosestApproach.trackStateAtPCA().position();
  
         edm::LogVerbatim("CosmicTrackingParticleSelector")
             << "FINAL State extrapolated at PCA: Radius = " << vertex.perp() << ", z = " << vertex.z()
             << ", pt = " << momentum.perp() << ", pz = " << momentum.z() << "\n";
  
         return (tpr->numberOfTrackerLayers() >= minHit_ && sqrt(momentum.perp2()) >= ptMin_ &&
                 momentum.eta() >= minRapidity_ && momentum.eta() <= maxRapidity_ && sqrt(vertex.perp2()) <= tip_ &&
                 fabs(vertex.z()) <= lip_);
       }
     }
   }

References cms::cuda::bs, chargedOnly_, PV3DBase< T, PVType, FrameType >::eta(), newFWLiteAna::found, edm::EventSetup::getData(), edm::EventSetup::getHandle(), globalTrackingGeomToken_, runTauDisplay::gp, GlobalTrackingGeometry::idToDet(), cuy::ii, TrajectoryStateClosestToBeamLine::isValid(), edm::HandleBase::isValid(), lip_, maxRapidity_, minHit_, minRapidity_, FreeTrajectoryState::momentum(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::perp2(), FreeTrajectoryState::position(), ptMin_, CosmicsPD_Skims::radius, FastTimerService_cff::range, simHitsTPAssoc, SimHitTPAssociationProducer::simHitTPAssociationListGreater(), mathSSE::sqrt(), GeomDet::surface(), theMFToken_, tip_, Surface::toGlobal(), TrajectoryStateClosestToBeamLine::trackStateAtPCA(), bphysicsOniaDQM_cfi::vertex, and PV3DBase< T, PVType, FrameType >::z().

◆ select()

void CosmicTrackingParticleSelector::select	(	const edm::Handle< collection > &	c,
		const edm::Event &	event,
		const edm::EventSetup &	setup
	)

inline

Definition at line 77 of file CosmicTrackingParticleSelector.h.

                                                                                                  {
     selected_.clear();
     edm::Handle<reco::BeamSpot> beamSpot;
     event.getByToken(beamSpotToken_, beamSpot);
     for (TrackingParticleCollection::const_iterator itp = c->begin(); itp != c->end(); ++itp)
       if (operator()(TrackingParticleRef(c, itp - c->begin()), beamSpot.product(), event, setup)) {
         selected_.push_back(&*itp);
       }
   }