#include <CosmicParametersDefinerForTP.h>

Inheritance diagram for CosmicParametersDefinerForTP:

Public Member Functions
std::unique_ptr < ParametersDefinerForTP >	clone () const override

	CosmicParametersDefinerForTP (edm::ConsumesCollector iC)

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

TrackingParticle::Vector	momentum (const edm::Event &iEvent, const edm::EventSetup &iSetup, const TrackingParticleRef &tpr) const override

TrackingParticle::Vector	momentum (const edm::Event &iEvent, const edm::EventSetup &iSetup, const Charge ch, const Point &vertex, const LorentzVector &lv) const override

std::tuple < TrackingParticle::Vector, TrackingParticle::Point >	momentumAndVertex (const edm::Event &iEvent, const edm::EventSetup &iSetup, const TrackingParticleRef &tpr) const override

TrackingParticle::Point	vertex (const edm::Event &iEvent, const edm::EventSetup &iSetup, const TrackingParticleRef &tpr) const override

TrackingParticle::Point	vertex (const edm::Event &iEvent, const edm::EventSetup &iSetup, const Charge ch, const Point &vertex, const LorentzVector &lv) const override

	~CosmicParametersDefinerForTP () override

Public Member Functions inherited from ParametersDefinerForTP
virtual TrackingParticle::Vector	momentum (const edm::Event &iEvent, const edm::EventSetup &iSetup, const reco::Candidate &tp) const

std::tuple < TrackingParticle::Vector, TrackingParticle::Point >	momentumAndVertex (const edm::Event &iEvent, const edm::EventSetup &iSetup, const Charge ch, const Point &vtx, const LorentzVector &lv) const

	ParametersDefinerForTP (const edm::InputTag &beamspot, edm::ConsumesCollector iC)

virtual TrackingParticle::Point	vertex (const edm::Event &iEvent, const edm::EventSetup &iSetup, const reco::Candidate &tp) const

virtual	~ParametersDefinerForTP ()

Private Attributes
const edm::ESGetToken < GlobalTrackingGeometry, GlobalTrackingGeometryRecord >	geometryToken_

edm::Handle < SimHitTPAssociationProducer::SimHitTPAssociationList >	simHitsTPAssoc

Additional Inherited Members
Public Types inherited from ParametersDefinerForTP
typedef int	Charge
	electric charge type More...

typedef math::XYZTLorentzVectorD	LorentzVector
	Lorentz vector. More...

typedef math::XYZPointD	Point
	point in the space More...

Protected Attributes inherited from ParametersDefinerForTP
const edm::EDGetTokenT < reco::BeamSpot >	bsToken_

const edm::ESGetToken < MagneticField, IdealMagneticFieldRecord >	mfToken_

Detailed Description

Author: Boris Mangano (UCSD) 5/7/2009

Definition at line 18 of file CosmicParametersDefinerForTP.h.

Constructor & Destructor Documentation

CosmicParametersDefinerForTP::CosmicParametersDefinerForTP ( edm::ConsumesCollector iC )

Definition at line 15 of file CosmicParametersDefinerForTP.cc.

16 : ParametersDefinerForTP(edm::InputTag("offlineBeamSpot"), iC), geometryToken_(iC.esConsumes()) {}

edm::ConsumesCollector::esConsumes

auto esConsumes()

Definition: ConsumesCollector.h:97

ParametersDefinerForTP::ParametersDefinerForTP

ParametersDefinerForTP(const edm::InputTag &beamspot, edm::ConsumesCollector iC)

Definition: ParametersDefinerForTP.cc:12

CosmicParametersDefinerForTP::geometryToken_

const edm::ESGetToken< GlobalTrackingGeometry, GlobalTrackingGeometryRecord > geometryToken_

Definition: CosmicParametersDefinerForTP.h:60

edm::InputTag

Definition: InputTag.h:15

CosmicParametersDefinerForTP::~CosmicParametersDefinerForTP ( )

overridedefault

Member Function Documentation

std::unique_ptr<ParametersDefinerForTP> CosmicParametersDefinerForTP::clone ( void ) const

inlineoverridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 55 of file CosmicParametersDefinerForTP.h.

                                                                {
     return std::make_unique<CosmicParametersDefinerForTP>(*this);
   }

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

inlineoverridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 51 of file CosmicParametersDefinerForTP.h.

References simHitsTPAssoc.

                                                                                                              {
     simHitsTPAssoc = simHitsTPAssocToSet;
   }

TrackingParticle::Vector CosmicParametersDefinerForTP::momentum	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const TrackingParticleRef &	tpr
	)		const

overridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 19 of file CosmicParametersDefinerForTP.cc.

References cms::cuda::bs, ParametersDefinerForTP::bsToken_, newFWLiteAna::found, relativeConstraints::geometry, geometryToken_, edm::Event::get(), edm::EventSetup::getData(), runTauDisplay::gp, TrajectoryStateClosestToBeamLine::isValid(), edm::HandleBase::isValid(), ParametersDefinerForTP::mfToken_, FreeTrajectoryState::momentum(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), CosmicsPD_Skims::radius, sistrip::SpyUtilities::range(), dt_dqm_sourceclient_common_cff::reco, simHitsTPAssoc, SimHitTPAssociationProducer::simHitTPAssociationListGreater(), mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), TrajectoryStateClosestToBeamLine::trackStateAtPCA(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by momentumAndVertex().

                                                                                                       {
   // to add a new implementation for cosmic. For the moment, it is just as for
   // the base class:
   using namespace edm;
   using namespace std;
   using namespace reco;
 
   auto const &bs = iEvent.get(bsToken_);
   auto const &geometry = iSetup.getData(geometryToken_);
   auto const &mf = iSetup.getData(mfToken_);
 
   GlobalVector finalGV(0, 0, 0);
   GlobalPoint finalGP(0, 0, 0);
   double radius(9999);
   bool found(false);
   TrackingParticle::Vector momentum(0, 0, 0);
 
   edm::LogVerbatim("CosmicParametersDefinerForTP") << "\t in CosmicParametersDefinerForTP::momentum";
   edm::LogVerbatim("CosmicParametersDefinerForTP")
       << "\t \t Original TP state:          pt = " << tpr->pt() << ", pz = " << tpr->pz();
 
   if (simHitsTPAssoc.isValid() == 0) {
     LogError("TrackAssociation") << "Invalid handle!";
     return momentum;
   }
   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;
     const GeomDet *tmpDet = geometry.idToDet(DetId(it->detUnitId()));
     if (!tmpDet) {
       edm::LogVerbatim("CosmicParametersDefinerForTP")
           << "***WARNING in CosmicParametersDefinerForTP::momentum: no GeomDet "
              "for: "
           << it->detUnitId() << ". Skipping it."
           << "\n";
       continue;
     }
 
     LocalVector lv = it->momentumAtEntry();
     Local3DPoint lp = it->localPosition();
     GlobalVector gv = tmpDet->surface().toGlobal(lv);
     GlobalPoint gp = tmpDet->surface().toGlobal(lp);
 
     // 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("CosmicParametersDefinerForTP")
       //   <<"\t FINAL State at InnerMost Hit: Radius = "<< finalGP.perp() << ",
       //   z = "<< finalGP.z()
       //  <<", pt = "<< finalGV.perp() << ", pz = "<< finalGV.z();
       << "\t \t FINAL State at InnerMost Hit:   pt = " << finalGV.perp() << ", pz = " << finalGV.z();
 
   if (found) {
     FreeTrajectoryState ftsAtProduction(finalGP, finalGV, TrackCharge(tpr->charge()), &mf);
     TSCBLBuilderNoMaterial tscblBuilder;
     TrajectoryStateClosestToBeamLine tsAtClosestApproach =
         tscblBuilder(ftsAtProduction, bs);  // as in TrackProducerAlgorithm
 
     if (tsAtClosestApproach.isValid()) {
       GlobalVector p = tsAtClosestApproach.trackStateAtPCA().momentum();
       momentum = TrackingParticle::Vector(p.x(), p.y(), p.z());
     } else {
       edm::LogVerbatim("CosmicParametersDefinerForTP") << "*** WARNING in CosmicParametersDefinerForTP::momentum: "
                                                           "tsAtClosestApproach is not valid."
                                                        << "\n";
     }
 
     edm::LogVerbatim("CosmicParametersDefinerForTP")
         << "\t \t FINAL State extrap. at PCA: pt = " << sqrt(momentum.x() * momentum.x() + momentum.y() * momentum.y())
         << ", pz = " << momentum.z() << "\n";
 
     return momentum;
   }
 
   edm::LogVerbatim("CosmicParametersDefinerForTP")
       << "*** WARNING in CosmicParametersDefinerForTP::momentum: NOT found the "
          "innermost TP point"
       << "\n";
   edm::LogVerbatim("CosmicParametersDefinerForTP")
       << "*** FINAL Reference MOMENTUM TP (px,py,pz) = " << momentum.x() << momentum.y() << momentum.z() << "\n";
   return momentum;
 }

TrackingParticle::Vector CosmicParametersDefinerForTP::momentum	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const Charge	ch,
		const Point &	vertex,
		const LorentzVector &	lv
	)		const

inlineoverridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 35 of file CosmicParametersDefinerForTP.h.

                                                                             {
     return TrackingParticle::Vector();
   }

std::tuple<TrackingParticle::Vector, TrackingParticle::Point> CosmicParametersDefinerForTP::momentumAndVertex	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const TrackingParticleRef &	tpr
	)		const

inlineoverridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 30 of file CosmicParametersDefinerForTP.h.

References momentum(), and vertex().

                                                                                                           {
     return std::make_tuple(momentum(iEvent, iSetup, tpr), vertex(iEvent, iSetup, tpr));
   }

TrackingParticle::Point CosmicParametersDefinerForTP::vertex	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const TrackingParticleRef &	tpr
	)		const

overridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 120 of file CosmicParametersDefinerForTP.cc.

References cms::cuda::bs, ParametersDefinerForTP::bsToken_, newFWLiteAna::found, relativeConstraints::geometry, geometryToken_, edm::Event::get(), edm::EventSetup::getData(), runTauDisplay::gp, TrajectoryStateClosestToBeamLine::isValid(), edm::HandleBase::isValid(), ParametersDefinerForTP::mfToken_, PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), CosmicsPD_Skims::radius, sistrip::SpyUtilities::range(), dt_dqm_sourceclient_common_cff::reco, simHitsTPAssoc, SimHitTPAssociationProducer::simHitTPAssociationListGreater(), mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), TrajectoryStateClosestToBeamLine::trackStateAtPCA(), findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by Tau.Tau::dxy(), and momentumAndVertex().

                                                                                                    {
   using namespace edm;
   using namespace std;
   using namespace reco;
 
   auto const &bs = iEvent.get(bsToken_);
   auto const &geometry = iSetup.getData(geometryToken_);
   auto const &mf = iSetup.getData(mfToken_);
 
   GlobalVector finalGV(0, 0, 0);
   GlobalPoint finalGP(0, 0, 0);
   double radius(9999);
   bool found(false);
   TrackingParticle::Point vertex(0, 0, 0);
 
   edm::LogVerbatim("CosmicParametersDefinerForTP") << "\t in CosmicParametersDefinerForTP::vertex";
   edm::LogVerbatim("CosmicParametersDefinerForTP")
       << "\t \t Original TP state:          radius = "
       << sqrt(tpr->vertex().x() * tpr->vertex().x() + tpr->vertex().y() * tpr->vertex().y())
       << ", z = " << tpr->vertex().z();
 
   if (simHitsTPAssoc.isValid() == 0) {
     LogError("TrackAssociation") << "Invalid handle!";
     return vertex;
   }
   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;
     const GeomDet *tmpDet = geometry.idToDet(DetId(it->detUnitId()));
     if (!tmpDet) {
       edm::LogVerbatim("CosmicParametersDefinerForTP")
           << "***WARNING in CosmicParametersDefinerForTP::vertex: no GeomDet "
              "for: "
           << it->detUnitId() << ". Skipping it."
           << "\n";
       continue;
     }
 
     LocalVector lv = it->momentumAtEntry();
     Local3DPoint lp = it->localPosition();
     GlobalVector gv = tmpDet->surface().toGlobal(lv);
     GlobalPoint gp = tmpDet->surface().toGlobal(lp);
 
     // 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("CosmicParametersDefinerForTP")
       << "\t \t FINAL State at InnerMost Hit:   radius = " << finalGP.perp() << ", z = " << finalGP.z();
 
   if (found) {
     FreeTrajectoryState ftsAtProduction(finalGP, finalGV, TrackCharge(tpr->charge()), &mf);
     TSCBLBuilderNoMaterial tscblBuilder;
     TrajectoryStateClosestToBeamLine tsAtClosestApproach =
         tscblBuilder(ftsAtProduction, bs);  // as in TrackProducerAlgorithm
 
     if (tsAtClosestApproach.isValid()) {
       GlobalPoint v = tsAtClosestApproach.trackStateAtPCA().position();
       vertex = TrackingParticle::Point(v.x(), v.y(), v.z());
     } else {
       // to preserve old behaviour
       // would be better to flag this somehow to allow ignoring in downstream
       vertex = TrackingParticle::Point(bs.x0(), bs.y0(), bs.z0());
       edm::LogVerbatim("CosmicParametersDefinerForTP") << "*** WARNING in CosmicParametersDefinerForTP::vertex: "
                                                           "tsAtClosestApproach is not valid."
                                                        << "\n";
     }
     edm::LogVerbatim("CosmicParametersDefinerForTP")
         << "\t \t FINAL State extrap. at PCA: radius = " << sqrt(vertex.x() * vertex.x() + vertex.y() * vertex.y())
         << ", z = " << vertex.z() << "\n";
 
     return vertex;
   }
 
   edm::LogVerbatim("CosmicParametersDefinerForTP")
       << "*** WARNING in CosmicParametersDefinerForTP::vertex: NOT found the "
          "innermost TP point"
       << "\n";
   edm::LogVerbatim("CosmicParametersDefinerForTP")
       << "*** FINAL Reference VERTEX TP   V(x,y,z) = " << vertex.x() << vertex.y() << vertex.z() << "\n";
 
   return vertex;
 }

TrackingParticle::Point CosmicParametersDefinerForTP::vertex	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const Charge	ch,
		const Point &	vertex,
		const LorentzVector &	lv
	)		const

inlineoverridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 43 of file CosmicParametersDefinerForTP.h.

Referenced by Tau.Tau::dxy().

                                                                          {
     return TrackingParticle::Point();
   }

Member Data Documentation

const edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> CosmicParametersDefinerForTP::geometryToken_

private

Definition at line 60 of file CosmicParametersDefinerForTP.h.

Referenced by momentum(), and vertex().

edm::Handle<SimHitTPAssociationProducer::SimHitTPAssociationList> CosmicParametersDefinerForTP::simHitsTPAssoc

private

Definition at line 61 of file CosmicParametersDefinerForTP.h.

Referenced by initEvent(), momentum(), and vertex().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation