CosmicParametersDefinerForTP Class Reference

#include <CosmicParametersDefinerForTP.h>

Inheritance diagram for CosmicParametersDefinerForTP:

Public Member Functions
std::unique_ptr< ParametersDefinerForTP >	clone () const override
	CosmicParametersDefinerForTP ()
void	initEvent (edm::Handle< SimHitTPAssociationProducer::SimHitTPAssociationList > simHitsTPAssocToSet) override
TrackingParticle::Vector	momentum (const edm::Event &iEvent, const edm::EventSetup &iSetup, const Charge ch, const Point &vertex, const LorentzVector &lv) const override
TrackingParticle::Vector	momentum (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
TrackingParticle::Point	vertex (const edm::Event &iEvent, const edm::EventSetup &iSetup, const TrackingParticleRef &tpr) 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
	ParametersDefinerForTP ()
	ParametersDefinerForTP (const edm::ParameterSet &iConfig)
virtual TrackingParticle::Point	vertex (const edm::Event &iEvent, const edm::EventSetup &iSetup, const reco::Candidate &tp) const
virtual	~ParametersDefinerForTP ()

Private Attributes
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...
Public Attributes inherited from ParametersDefinerForTP
edm::InputTag	beamSpotInputTag_

Detailed Description

Author

Boris Mangano (UCSD) 5/7/2009

Definition at line 13 of file CosmicParametersDefinerForTP.h.

Constructor & Destructor Documentation

CosmicParametersDefinerForTP()

CosmicParametersDefinerForTP::CosmicParametersDefinerForTP ( )

inline

Definition at line 15 of file CosmicParametersDefinerForTP.h.

{};

Referenced by clone().

~CosmicParametersDefinerForTP()

CosmicParametersDefinerForTP::~CosmicParametersDefinerForTP ( )

inlineoverride

Definition at line 16 of file CosmicParametersDefinerForTP.h.

{};

Member Function Documentation

clone()

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

inlineoverridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 45 of file CosmicParametersDefinerForTP.h.

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

References CosmicParametersDefinerForTP().

Referenced by MuonTrackValidator::analyze().

initEvent()

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

inlineoverridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 41 of file CosmicParametersDefinerForTP.h.

                                                                                                             {
    simHitsTPAssoc = simHitsTPAssocToSet;
  }

References simHitsTPAssoc.

momentum() [1/2]

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 25 of file CosmicParametersDefinerForTP.h.

                                                                            {
    return TrackingParticle::Vector();
  }

momentum() [2/2]

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

overridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 23 of file CosmicParametersDefinerForTP.cc.

                                                                                                      {
  // 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;
 
  ESHandle<TrackerGeometry> tracker;
  iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
  edm::ESHandle<GlobalTrackingGeometry> theGeometry;
  iSetup.get<GlobalTrackingGeometryRecord>().get(theGeometry);
 
  edm::ESHandle<MagneticField> theMF;
  iSetup.get<IdealMagneticFieldRecord>().get(theMF);
 
  edm::Handle<reco::BeamSpot> bs;
  iEvent.getByLabel(InputTag("offlineBeamSpot"), bs);
 
  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 = theGeometry->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()), theMF.product());
    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;
}

References cms::cuda::bs, newFWLiteAna::found, edm::EventSetup::get(), get, runTauDisplay::gp, GlobalTrackingGeometry::idToDet(), iEvent, HLT_FULL_cff::InputTag, TrajectoryStateClosestToBeamLine::isValid(), edm::HandleBase::isValid(), FreeTrajectoryState::momentum(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), edm::ESHandle< T >::product(), CosmicsPD_Skims::radius, FastTimerService_cff::range, simHitsTPAssoc, SimHitTPAssociationProducer::simHitTPAssociationListGreater(), mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), PbPb_ZMuSkimMuonDPG_cff::tracker, TrajectoryStateClosestToBeamLine::trackStateAtPCA(), and PV3DBase< T, PVType, FrameType >::z().

vertex() [1/2]

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 33 of file CosmicParametersDefinerForTP.h.

                                                                         {
    return TrackingParticle::Point();
  }

Referenced by Tau.Tau::dxy().

vertex() [2/2]

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

overridevirtual

Reimplemented from ParametersDefinerForTP.

Definition at line 131 of file CosmicParametersDefinerForTP.cc.

                                                                                                   {
  using namespace edm;
  using namespace std;
  using namespace reco;
 
  ESHandle<TrackerGeometry> tracker;
  iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
  edm::ESHandle<GlobalTrackingGeometry> theGeometry;
  iSetup.get<GlobalTrackingGeometryRecord>().get(theGeometry);
 
  edm::ESHandle<MagneticField> theMF;
  iSetup.get<IdealMagneticFieldRecord>().get(theMF);
 
  edm::Handle<reco::BeamSpot> bs;
  iEvent.getByLabel(InputTag("offlineBeamSpot"), bs);
 
  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 = theGeometry->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()), theMF.product());
    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;
}

References cms::cuda::bs, newFWLiteAna::found, edm::EventSetup::get(), get, runTauDisplay::gp, GlobalTrackingGeometry::idToDet(), iEvent, HLT_FULL_cff::InputTag, TrajectoryStateClosestToBeamLine::isValid(), edm::HandleBase::isValid(), PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), edm::ESHandle< T >::product(), CosmicsPD_Skims::radius, FastTimerService_cff::range, simHitsTPAssoc, SimHitTPAssociationProducer::simHitTPAssociationListGreater(), mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), PbPb_ZMuSkimMuonDPG_cff::tracker, TrajectoryStateClosestToBeamLine::trackStateAtPCA(), findQualityFiles::v, and PV3DBase< T, PVType, FrameType >::z().

Referenced by Tau.Tau::dxy().

Member Data Documentation

simHitsTPAssoc

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

private

Definition at line 50 of file CosmicParametersDefinerForTP.h.

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