CosmicTrackingParticleSelector Class Reference

#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 (double ptMin, double minRapidity, double maxRapidity, double tip, double lip, int minHit, bool chargedOnly, const std::vector< int > &pdgId=std::vector< int >())
	CosmicTrackingParticleSelector (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
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

typedef TrackingParticleCollection CosmicTrackingParticleSelector::collection

Definition at line 41 of file CosmicTrackingParticleSelector.h.

typedef container::const_iterator CosmicTrackingParticleSelector::const_iterator

Definition at line 43 of file CosmicTrackingParticleSelector.h.

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

Definition at line 42 of file CosmicTrackingParticleSelector.h.

Constructor & Destructor Documentation

CosmicTrackingParticleSelector::CosmicTrackingParticleSelector ( )

inline

Definition at line 45 of file CosmicTrackingParticleSelector.h.

{}

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::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

const_iterator CosmicTrackingParticleSelector::begin ( void  ) const

inline

Definition at line 87 of file CosmicTrackingParticleSelector.h.

References selected_.

{ return selected_.begin(); }

const_iterator CosmicTrackingParticleSelector::end ( void  ) const

inline

Definition at line 88 of file CosmicTrackingParticleSelector.h.

References selected_.

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

{ return selected_.end(); }

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

inline

Definition at line 90 of file CosmicTrackingParticleSelector.h.

References simHitsTPAssoc.

Referenced by MultiTrackValidator::dqmAnalyze().

                                                                                                          {
    simHitsTPAssoc = simHitsTPAssocToSet;
  }

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.

References chargedOnly_, PV3DBase< T, PVType, FrameType >::eta(), newFWLiteAna::found, edm::EventSetup::getData(), edm::EventSetup::getHandle(), globalTrackingGeomToken_, runTauDisplay::gp, 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, sistrip::SpyUtilities::range(), simHitsTPAssoc, SimHitTPAssociationProducer::simHitTPAssociationListGreater(), mathSSE::sqrt(), GeomDet::surface(), theMFToken_, tip_, Surface::toGlobal(), TrajectoryStateClosestToBeamLine::trackStateAtPCA(), and PV3DBase< T, PVType, FrameType >::z().

                                                     {
    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_);
      }
    }
  }

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.

References beam_dqm_sourceclient-live_cfg::beamSpot, beamSpotToken_, edmPickEvents::event, edm::Handle< T >::product(), selected_, and GeneralSetup::setup().

                                                                                                 {
    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);
      }
  }

size_t CosmicTrackingParticleSelector::size ( void  ) const

inline

Definition at line 203 of file CosmicTrackingParticleSelector.h.

References selected_.

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

{ return selected_.size(); }

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot> CosmicTrackingParticleSelector::beamSpotToken_

private

Definition at line 215 of file CosmicTrackingParticleSelector.h.

Referenced by select().

bool CosmicTrackingParticleSelector::chargedOnly_

private

Definition at line 212 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().

edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> CosmicTrackingParticleSelector::globalTrackingGeomToken_

private

Definition at line 216 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().

double CosmicTrackingParticleSelector::lip_

private

Definition at line 210 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().

double CosmicTrackingParticleSelector::maxRapidity_

private

Definition at line 208 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().

int CosmicTrackingParticleSelector::minHit_

private

Definition at line 211 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().

double CosmicTrackingParticleSelector::minRapidity_

private

Definition at line 207 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().

std::vector<int> CosmicTrackingParticleSelector::pdgId_

private

Definition at line 213 of file CosmicTrackingParticleSelector.h.

double CosmicTrackingParticleSelector::ptMin_

private

Definition at line 206 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().

container CosmicTrackingParticleSelector::selected_

private

Definition at line 214 of file CosmicTrackingParticleSelector.h.

Referenced by begin(), end(), select(), and size().

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

mutableprivate

Definition at line 219 of file CosmicTrackingParticleSelector.h.

Referenced by initEvent(), and operator()().

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> CosmicTrackingParticleSelector::theMFToken_

private

Definition at line 217 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().

double CosmicTrackingParticleSelector::tip_

private

Definition at line 209 of file CosmicTrackingParticleSelector.h.

Referenced by operator()().