#include <RecoTracker/DeDx/plugins/DeDxHitInfoProducer.cc>

Inheritance diagram for DeDxHitInfoProducer:

Public Member Functions
	DeDxHitInfoProducer (const edm::ParameterSet &)

	~DeDxHitInfoProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	beginRun (edm::Run const &run, const edm::EventSetup &) override

void	makeCalibrationMap (const TrackerGeometry &tkGeom)

void	processHit (const TrackingRecHit *recHit, const float trackMomentum, const float cosine, reco::DeDxHitInfo &hitDeDxInfo, const LocalPoint &hitLocalPos)

void	produce (edm::Event &, const edm::EventSetup &) override

uint64_t	xorshift128p (uint64_t state[2])

Private Attributes
std::vector< std::vector< float > >	calibGains

const float	lowPtTracksDeDxThreshold

GenericTruncatedAverageDeDxEstimator	lowPtTracksEstimator

const unsigned int	lowPtTracksPrescaleFail

const unsigned int	lowPtTracksPrescalePass

const std::string	m_calibrationPath

unsigned int	m_off

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	m_tkGeomToken

edm::EDGetTokenT< reco::TrackCollection >	m_tracksTag

const float	maxTrackEta

const float	MeVperADCPixel

const float	MeVperADCStrip

const unsigned int	minTrackHits

const float	minTrackPt

const float	minTrackPtPrescale

const bool	shapetest

edm::ESHandle< TrackerGeometry >	tkGeom

const bool	useCalibration

const bool	usePixel

const bool	useStrip

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 35 of file DeDxHitInfoProducer.h.

Constructor & Destructor Documentation

◆ DeDxHitInfoProducer()

DeDxHitInfoProducer::DeDxHitInfoProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 26 of file DeDxHitInfoProducer.cc.

     : usePixel(iConfig.getParameter<bool>("usePixel")),
       useStrip(iConfig.getParameter<bool>("useStrip")),
       MeVperADCPixel(iConfig.getParameter<double>("MeVperADCPixel")),
       MeVperADCStrip(iConfig.getParameter<double>("MeVperADCStrip")),
       minTrackHits(iConfig.getParameter<unsigned>("minTrackHits")),
       minTrackPt(iConfig.getParameter<double>("minTrackPt")),
       minTrackPtPrescale(iConfig.getParameter<double>("minTrackPtPrescale")),
       maxTrackEta(iConfig.getParameter<double>("maxTrackEta")),
       m_calibrationPath(iConfig.getParameter<string>("calibrationPath")),
       useCalibration(iConfig.getParameter<bool>("useCalibration")),
       shapetest(iConfig.getParameter<bool>("shapeTest")),
       lowPtTracksPrescalePass(iConfig.getParameter<uint32_t>("lowPtTracksPrescalePass")),
       lowPtTracksPrescaleFail(iConfig.getParameter<uint32_t>("lowPtTracksPrescaleFail")),
       lowPtTracksEstimator(iConfig.getParameter<edm::ParameterSet>("lowPtTracksEstimatorParameters")),
       lowPtTracksDeDxThreshold(iConfig.getParameter<double>("lowPtTracksDeDxThreshold")) {
   produces<reco::DeDxHitInfoCollection>();
   produces<reco::DeDxHitInfoAss>();
   produces<edm::ValueMap<int>>("prescale");
  
   m_tracksTag = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("tracks"));
   m_tkGeomToken = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>();
  
   if (!usePixel && !useStrip)
     edm::LogError("DeDxHitsProducer") << "No Pixel Hits NOR Strip Hits will be saved.  Running this module is useless";
 }

References edm::ParameterSet::getParameter(), m_tkGeomToken, m_tracksTag, usePixel, and useStrip.

◆ ~DeDxHitInfoProducer()

DeDxHitInfoProducer::~DeDxHitInfoProducer ( )

override

Definition at line 53 of file DeDxHitInfoProducer.cc.

53 {}

Member Function Documentation

◆ beginRun()

void DeDxHitInfoProducer::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 56 of file DeDxHitInfoProducer.cc.

                                                                                  {
   tkGeom = iSetup.getHandle(m_tkGeomToken);
   if (useCalibration && calibGains.empty()) {
     m_off = tkGeom->offsetDU(GeomDetEnumerators::PixelBarrel);  //index start at the first pixel
  
     DeDxTools::makeCalibrationMap(m_calibrationPath, *tkGeom, calibGains, m_off);
   }
 }

References calibGains, edm::EventSetup::getHandle(), m_calibrationPath, m_off, m_tkGeomToken, DeDxTools::makeCalibrationMap(), TrackerGeometry::offsetDU(), GeomDetEnumerators::PixelBarrel, tkGeom, and useCalibration.

◆ makeCalibrationMap()

void DeDxHitInfoProducer::makeCalibrationMap ( const TrackerGeometry & tkGeom )

private

◆ processHit()

void DeDxHitInfoProducer::processHit	(	const TrackingRecHit *	recHit,
		const float	trackMomentum,
		const float	cosine,
		reco::DeDxHitInfo &	hitDeDxInfo,
		const LocalPoint &	hitLocalPos
	)

private

Definition at line 161 of file DeDxHitInfoProducer.cc.

                                                                     {
   auto const& thit = static_cast<BaseTrackerRecHit const&>(*recHit);
   if (!thit.isValid())
     return;
  
   float cosineAbs = std::max(0.00000001f, std::abs(cosine));  //make sure cosine is not 0
  
   auto const& clus = thit.firstClusterRef();
   if (!clus.isValid())
     return;
  
   if (clus.isPixel()) {
     if (!usePixel)
       return;
  
     const auto* detUnit = recHit->detUnit();
     if (detUnit == nullptr)
       detUnit = tkGeom->idToDet(thit.geographicalId());
     float pathLen = detUnit->surface().bounds().thickness() / cosineAbs;
     float chargeAbs = clus.pixelCluster().charge();
     hitDeDxInfo.addHit(chargeAbs, pathLen, thit.geographicalId(), hitLocalPos, clus.pixelCluster());
   } else if (clus.isStrip() && !thit.isMatched()) {
     if (!useStrip)
       return;
  
     const auto* detUnit = recHit->detUnit();
     if (detUnit == nullptr)
       detUnit = tkGeom->idToDet(thit.geographicalId());
     int NSaturating = 0;
     float pathLen = detUnit->surface().bounds().thickness() / cosineAbs;
     float chargeAbs = DeDxTools::getCharge(&(clus.stripCluster()), NSaturating, *detUnit, calibGains, m_off);
     hitDeDxInfo.addHit(chargeAbs, pathLen, thit.geographicalId(), hitLocalPos, clus.stripCluster());
   } else if (clus.isStrip() && thit.isMatched()) {
     if (!useStrip)
       return;
     const SiStripMatchedRecHit2D* matchedHit = dynamic_cast<const SiStripMatchedRecHit2D*>(recHit);
     if (!matchedHit)
       return;
  
     const auto* detUnitM = matchedHit->monoHit().detUnit();
     if (detUnitM == nullptr)
       detUnitM = tkGeom->idToDet(matchedHit->monoHit().geographicalId());
     int NSaturating = 0;
     float pathLen = detUnitM->surface().bounds().thickness() / cosineAbs;
     float chargeAbs =
         DeDxTools::getCharge(&(matchedHit->monoHit().stripCluster()), NSaturating, *detUnitM, calibGains, m_off);
     hitDeDxInfo.addHit(chargeAbs, pathLen, thit.geographicalId(), hitLocalPos, matchedHit->monoHit().stripCluster());
  
     const auto* detUnitS = matchedHit->stereoHit().detUnit();
     if (detUnitS == nullptr)
       detUnitS = tkGeom->idToDet(matchedHit->stereoHit().geographicalId());
     NSaturating = 0;
     pathLen = detUnitS->surface().bounds().thickness() / cosineAbs;
     chargeAbs =
         DeDxTools::getCharge(&(matchedHit->stereoHit().stripCluster()), NSaturating, *detUnitS, calibGains, m_off);
     hitDeDxInfo.addHit(chargeAbs, pathLen, thit.geographicalId(), hitLocalPos, matchedHit->stereoHit().stripCluster());
   }
 }

References funct::abs(), reco::DeDxHitInfo::addHit(), Surface::bounds(), calibGains, TrackerSingleRecHit::detUnit(), f, TrackingRecHit::geographicalId(), DeDxTools::getCharge(), TrackerGeometry::idToDet(), m_off, SiStripPI::max, SiStripMatchedRecHit2D::monoHit(), rpcPointValidation_cfi::recHit, SiStripMatchedRecHit2D::stereoHit(), TrackerSingleRecHit::stripCluster(), GeomDet::surface(), Bounds::thickness(), tkGeom, usePixel, and useStrip.

Referenced by produce().

◆ produce()

void DeDxHitInfoProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 65 of file DeDxHitInfoProducer.cc.

                                                                                {
   edm::Handle<reco::TrackCollection> trackCollectionHandle;
   iEvent.getByToken(m_tracksTag, trackCollectionHandle);
   const TrackCollection& trackCollection(*trackCollectionHandle.product());
  
   // creates the output collection
   auto resultdedxHitColl = std::make_unique<reco::DeDxHitInfoCollection>();
  
   std::vector<int> indices;
   std::vector<int> prescales;
   uint64_t state[2] = {iEvent.id().event(), iEvent.id().luminosityBlock()};
   for (unsigned int j = 0; j < trackCollection.size(); j++) {
     const reco::Track& track = trackCollection[j];
  
     //track selection
     bool passPt = (track.pt() >= minTrackPt), passLowDeDx = false, passHighDeDx = false, pass = passPt;
     if (!pass && (track.pt() >= minTrackPtPrescale)) {
       if (lowPtTracksPrescalePass > 0) {
         passHighDeDx = ((xorshift128p(state) % lowPtTracksPrescalePass) == 0);
       }
       if (lowPtTracksPrescaleFail > 0) {
         passLowDeDx = ((xorshift128p(state) % lowPtTracksPrescaleFail) == 0);
       }
       pass = passHighDeDx || passLowDeDx;
     }
     if (!pass || std::abs(track.eta()) > maxTrackEta || track.numberOfValidHits() < minTrackHits) {
       indices.push_back(-1);
       continue;
     }
  
     reco::DeDxHitInfo hitDeDxInfo;
     auto const& trajParams = track.extra()->trajParams();
     auto hb = track.recHitsBegin();
     for (unsigned int h = 0; h < track.recHitsSize(); h++) {
       auto recHit = *(hb + h);
       if (!trackerHitRTTI::isFromDet(*recHit))
         continue;
  
       auto trackDirection = trajParams[h].direction();
       float cosine = trackDirection.z() / trackDirection.mag();
  
       processHit(recHit, track.p(), cosine, hitDeDxInfo, trajParams[h].position());
     }
  
     if (!passPt) {
       std::vector<DeDxHit> hits;
       hits.reserve(hitDeDxInfo.size());
       for (unsigned int i = 0, n = hitDeDxInfo.size(); i < n; ++i) {
         if (hitDeDxInfo.detId(i).subdetId() <= 2) {
           hits.push_back(DeDxHit(hitDeDxInfo.charge(i) / hitDeDxInfo.pathlength(i) * MeVperADCPixel, 0, 0, 0));
         } else {
           if (shapetest && !DeDxTools::shapeSelection(*hitDeDxInfo.stripCluster(i)))
             continue;
           hits.push_back(DeDxHit(hitDeDxInfo.charge(i) / hitDeDxInfo.pathlength(i) * MeVperADCStrip, 0, 0, 0));
         }
       }
       std::sort(hits.begin(), hits.end(), std::less<DeDxHit>());
       if (lowPtTracksEstimator.dedx(hits).first < lowPtTracksDeDxThreshold) {
         if (passLowDeDx) {
           prescales.push_back(lowPtTracksPrescaleFail);
         } else {
           indices.push_back(-1);
           continue;
         }
       } else {
         if (passHighDeDx) {
           prescales.push_back(lowPtTracksPrescalePass);
         } else {
           indices.push_back(-1);
           continue;
         }
       }
     } else {
       prescales.push_back(1);
     }
     indices.push_back(resultdedxHitColl->size());
     resultdedxHitColl->push_back(hitDeDxInfo);
   }
  
   edm::OrphanHandle<reco::DeDxHitInfoCollection> dedxHitCollHandle = iEvent.put(std::move(resultdedxHitColl));
  
   //create map passing the handle to the matched collection
   auto dedxMatch = std::make_unique<reco::DeDxHitInfoAss>(dedxHitCollHandle);
   reco::DeDxHitInfoAss::Filler filler(*dedxMatch);
   filler.insert(trackCollectionHandle, indices.begin(), indices.end());
   filler.fill();
   iEvent.put(std::move(dedxMatch));
  
   auto dedxPrescale = std::make_unique<edm::ValueMap<int>>();
   edm::ValueMap<int>::Filler pfiller(*dedxPrescale);
   pfiller.insert(dedxHitCollHandle, prescales.begin(), prescales.end());
   pfiller.fill();
   iEvent.put(std::move(dedxPrescale), "prescale");
 }

References funct::abs(), reco::DeDxHitInfo::charge(), GenericTruncatedAverageDeDxEstimator::dedx(), reco::DeDxHitInfo::detId(), trigObjTnPSource_cfi::filler, h, hcalSimParameters_cfi::hb, hfClusterShapes_cfi::hits, mps_fire::i, iEvent, dqmdumpme::indices, trackerHitRTTI::isFromDet(), dqmiolumiharvest::j, lowPtTracksDeDxThreshold, lowPtTracksEstimator, lowPtTracksPrescaleFail, lowPtTracksPrescalePass, m_tracksTag, maxTrackEta, MeVperADCPixel, MeVperADCStrip, minTrackHits, minTrackPt, minTrackPtPrescale, eostools::move(), dqmiodumpmetadata::n, reco::DeDxHitInfo::pathlength(), hltrates_dqm_sourceclient-live_cfg::prescales, processHit(), edm::Handle< T >::product(), rpcPointValidation_cfi::recHit, DeDxTools::shapeSelection(), shapetest, reco::DeDxHitInfo::size(), jetUpdater_cfi::sort, reco::DeDxHitInfo::stripCluster(), DetId::subdetId(), HLT_FULL_cff::track, duplicaterechits_cfi::trackCollection, and xorshift128p().

◆ xorshift128p()

uint64_t DeDxHitInfoProducer::xorshift128p ( uint64_t state[2] )

inlineprivate

Definition at line 79 of file DeDxHitInfoProducer.h.

                                            {
     uint64_t x = state[0];
     uint64_t const y = state[1];
     state[0] = y;
     x ^= x << 23;                              // a
     state[1] = x ^ y ^ (x >> 17) ^ (y >> 26);  // b, c
     return state[1] + y;
   }

References x, and y.

Referenced by produce().

Member Data Documentation

◆ calibGains

std::vector<std::vector<float> > DeDxHitInfoProducer::calibGains

private

Definition at line 73 of file DeDxHitInfoProducer.h.

Referenced by beginRun(), and processHit().

◆ lowPtTracksDeDxThreshold

const float DeDxHitInfoProducer::lowPtTracksDeDxThreshold

private

Definition at line 71 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ lowPtTracksEstimator

GenericTruncatedAverageDeDxEstimator DeDxHitInfoProducer::lowPtTracksEstimator

private

Definition at line 70 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ lowPtTracksPrescaleFail

const unsigned int DeDxHitInfoProducer::lowPtTracksPrescaleFail

private

Definition at line 69 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ lowPtTracksPrescalePass

const unsigned int DeDxHitInfoProducer::lowPtTracksPrescalePass

private

Definition at line 69 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ m_calibrationPath

const std::string DeDxHitInfoProducer::m_calibrationPath

private

Definition at line 65 of file DeDxHitInfoProducer.h.

Referenced by beginRun().

◆ m_off

unsigned int DeDxHitInfoProducer::m_off

private

Definition at line 74 of file DeDxHitInfoProducer.h.

Referenced by beginRun(), and processHit().

◆ m_tkGeomToken

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> DeDxHitInfoProducer::m_tkGeomToken

private

Definition at line 77 of file DeDxHitInfoProducer.h.

Referenced by beginRun(), and DeDxHitInfoProducer().

◆ m_tracksTag

edm::EDGetTokenT<reco::TrackCollection> DeDxHitInfoProducer::m_tracksTag

private

Definition at line 53 of file DeDxHitInfoProducer.h.

Referenced by DeDxHitInfoProducer(), and produce().

◆ maxTrackEta

const float DeDxHitInfoProducer::maxTrackEta

private

Definition at line 63 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ MeVperADCPixel

const float DeDxHitInfoProducer::MeVperADCPixel

private

Definition at line 57 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ MeVperADCStrip

const float DeDxHitInfoProducer::MeVperADCStrip

private

Definition at line 58 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ minTrackHits

const unsigned int DeDxHitInfoProducer::minTrackHits

private

Definition at line 60 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ minTrackPt

const float DeDxHitInfoProducer::minTrackPt

private

Definition at line 61 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ minTrackPtPrescale

const float DeDxHitInfoProducer::minTrackPtPrescale

private

Definition at line 62 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ shapetest

const bool DeDxHitInfoProducer::shapetest

private

Definition at line 67 of file DeDxHitInfoProducer.h.

Referenced by produce().

◆ tkGeom

edm::ESHandle<TrackerGeometry> DeDxHitInfoProducer::tkGeom

private

Definition at line 76 of file DeDxHitInfoProducer.h.

Referenced by beginRun(), and processHit().

◆ useCalibration

const bool DeDxHitInfoProducer::useCalibration

private

Definition at line 66 of file DeDxHitInfoProducer.h.

Referenced by beginRun().

◆ usePixel

const bool DeDxHitInfoProducer::usePixel

private

Definition at line 55 of file DeDxHitInfoProducer.h.

Referenced by DeDxHitInfoProducer(), and processHit().

◆ useStrip

const bool DeDxHitInfoProducer::useStrip

private

Definition at line 56 of file DeDxHitInfoProducer.h.

Referenced by DeDxHitInfoProducer(), and processHit().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DeDxHitInfoProducer()

◆ ~DeDxHitInfoProducer()

Member Function Documentation

◆ beginRun()

◆ makeCalibrationMap()

◆ processHit()

◆ produce()

◆ xorshift128p()

Member Data Documentation

◆ calibGains

◆ lowPtTracksDeDxThreshold

◆ lowPtTracksEstimator

◆ lowPtTracksPrescaleFail

◆ lowPtTracksPrescalePass

◆ m_calibrationPath

◆ m_off

◆ m_tkGeomToken

◆ m_tracksTag

◆ maxTrackEta

◆ MeVperADCPixel

◆ MeVperADCStrip

◆ minTrackHits

◆ minTrackPt

◆ minTrackPtPrescale

◆ shapetest

◆ tkGeom

◆ useCalibration

◆ usePixel

◆ useStrip