#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 std::string	calibrationPath_

const bool	doShapeTest_

const float	lowPtTracksDeDxThreshold_

GenericTruncatedAverageDeDxEstimator	lowPtTracksEstimator_

const unsigned int	lowPtTracksPrescaleFail_

const unsigned int	lowPtTracksPrescalePass_

const float	maxTrackEta_

const unsigned int	minTrackHits_

const float	minTrackPt_

const float	minTrackPtPrescale_

unsigned int	offsetDU_

const float	theMeVperADCPixel_

const float	theMeVperADCStrip_

edm::ESHandle< TrackerGeometry >	tkGeom_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

const edm::EDGetTokenT< reco::TrackCollection >	tracksToken_

const bool	useCalibration_

const bool	usePixel_

const bool	usePixelForPrescales_

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 45 of file DeDxHitInfoProducer.cc.

Constructor & Destructor Documentation

◆ DeDxHitInfoProducer()

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

explicit

Definition at line 98 of file DeDxHitInfoProducer.cc.

References edm::BeginRun, usePixel_, and useStrip_.

     : usePixel_(iConfig.getParameter<bool>("usePixel")),
       useStrip_(iConfig.getParameter<bool>("useStrip")),
       theMeVperADCPixel_(iConfig.getParameter<double>("MeVperADCPixel")),
       theMeVperADCStrip_(iConfig.getParameter<double>("MeVperADCStrip")),
       minTrackHits_(iConfig.getParameter<unsigned>("minTrackHits")),
       minTrackPt_(iConfig.getParameter<double>("minTrackPt")),
       minTrackPtPrescale_(iConfig.getParameter<double>("minTrackPtPrescale")),
       maxTrackEta_(iConfig.getParameter<double>("maxTrackEta")),
       calibrationPath_(iConfig.getParameter<string>("calibrationPath")),
       useCalibration_(iConfig.getParameter<bool>("useCalibration")),
       doShapeTest_(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")),
       usePixelForPrescales_(iConfig.getParameter<bool>("usePixelForPrescales")),
       tracksToken_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("tracks"))),
       tkGeomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()) {
   produces<reco::DeDxHitInfoCollection>();
   produces<reco::DeDxHitInfoAss>();
   produces<edm::ValueMap<int>>("prescale");
 
   if (!usePixel_ && !useStrip_)
     edm::LogError("DeDxHitsProducer") << "No Pixel Hits NOR Strip Hits will be saved.  Running this module is useless";
 }

◆ ~DeDxHitInfoProducer()

DeDxHitInfoProducer::~DeDxHitInfoProducer ( )

overridedefault

Member Function Documentation

◆ beginRun()

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

overrideprivate

Definition at line 128 of file DeDxHitInfoProducer.cc.

References calibGains_, calibrationPath_, edm::EventSetup::getHandle(), deDxTools::makeCalibrationMap(), TrackerGeometry::offsetDU(), offsetDU_, GeomDetEnumerators::PixelBarrel, tkGeom_, tkGeomToken_, and useCalibration_.

                                                                                   {
   tkGeom_ = iSetup.getHandle(tkGeomToken_);
   if (useCalibration_ && calibGains_.empty()) {
     offsetDU_ = tkGeom_->offsetDU(GeomDetEnumerators::PixelBarrel);  //index start at the first pixel
 
     deDxTools::makeCalibrationMap(calibrationPath_, *tkGeom_, calibGains_, offsetDU_);
   }
 }

◆ 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 239 of file DeDxHitInfoProducer.cc.

References funct::abs(), reco::DeDxHitInfo::addHit(), Surface::bounds(), calibGains_, runTheMatrix::const, TrackerSingleRecHit::detUnit(), f, TrackingRecHit::geographicalId(), deDxTools::getCharge(), TrackerGeometry::idToDet(), SiStripPI::max, SiStripMatchedRecHit2D::monoHit(), offsetDU_, rpcPointValidation_cfi::recHit, SiStripMatchedRecHit2D::stereoHit(), TrackerSingleRecHit::stripCluster(), GeomDet::surface(), Bounds::thickness(), tkGeom_, usePixel_, and useStrip_.

Referenced by produce().

                                                                     {
   auto const& thit = static_cast<BaseTrackerRecHit const&>(*recHit);
   if (!thit.isValid())
     return;
 
   //make sure cosine is not 0
   float cosineAbs = std::max(0.00000001f, std::abs(cosine));
 
   auto const& clus = thit.firstClusterRef();
   if (!clus.isValid())
     return;
 
   const auto* detUnit = recHit->detUnit();
   if (detUnit == nullptr) {
     detUnit = tkGeom_->idToDet(thit.geographicalId());
   }
   float pathLen = detUnit->surface().bounds().thickness() / cosineAbs;
 
   if (clus.isPixel()) {
     if (!usePixel_)
       return;
 
     float chargeAbs = clus.pixelCluster().charge();
     hitDeDxInfo.addHit(chargeAbs, pathLen, thit.geographicalId(), hitLocalPos, clus.pixelCluster());
   } else if (clus.isStrip() && !thit.isMatched()) {
     if (!useStrip_)
       return;
 
     int NSaturating = 0;
     float chargeAbs = deDxTools::getCharge(&(clus.stripCluster()), NSaturating, *detUnit, calibGains_, offsetDU_);
     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;
     auto pathLenM = detUnitM->surface().bounds().thickness() / cosineAbs;
     float chargeAbs =
         deDxTools::getCharge(&(matchedHit->monoHit().stripCluster()), NSaturating, *detUnitM, calibGains_, offsetDU_);
     hitDeDxInfo.addHit(chargeAbs, pathLenM, thit.geographicalId(), hitLocalPos, matchedHit->monoHit().stripCluster());
 
     const auto* detUnitS = matchedHit->stereoHit().detUnit();
     if (detUnitS == nullptr)
       detUnitS = tkGeom_->idToDet(matchedHit->stereoHit().geographicalId());
     NSaturating = 0;
     auto pathLenS = detUnitS->surface().bounds().thickness() / cosineAbs;
     chargeAbs =
         deDxTools::getCharge(&(matchedHit->stereoHit().stripCluster()), NSaturating, *detUnitS, calibGains_, offsetDU_);
     hitDeDxInfo.addHit(chargeAbs, pathLenS, thit.geographicalId(), hitLocalPos, matchedHit->stereoHit().stripCluster());
   }
 }

◆ produce()

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

overrideprivate

Definition at line 137 of file DeDxHitInfoProducer.cc.

References funct::abs(), reco::DeDxHitInfo::charge(), GenericTruncatedAverageDeDxEstimator::dedx(), reco::DeDxHitInfo::detId(), doShapeTest_, trigObjTnPSource_cfi::filler, h, hcalSimParameters_cfi::hb, hfClusterShapes_cfi::hits, mps_fire::i, iEvent, dqmdumpme::indices, trackerHitRTTI::isFromDet(), dqmiolumiharvest::j, lowPtTracksDeDxThreshold_, lowPtTracksEstimator_, lowPtTracksPrescaleFail_, lowPtTracksPrescalePass_, maxTrackEta_, 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(), reco::DeDxHitInfo::size(), jetUpdater_cfi::sort, reco::DeDxHitInfo::stripCluster(), DetId::subdetId(), theMeVperADCPixel_, theMeVperADCStrip_, HLT_2022v12_cff::track, duplicaterechits_cfi::trackCollection, tracksToken_, usePixelForPrescales_, and xorshift128p().

                                                                                {
   edm::Handle<reco::TrackCollection> trackCollectionHandle;
   iEvent.getByToken(tracksToken_, 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 && usePixelForPrescales_) {
           hits.push_back(DeDxHit(hitDeDxInfo.charge(i) / hitDeDxInfo.pathlength(i) * theMeVperADCPixel_, 0, 0, 0));
         } else if (hitDeDxInfo.detId(i).subdetId() > 2) {
           if (doShapeTest_ && !deDxTools::shapeSelection(*hitDeDxInfo.stripCluster(i)))
             continue;
           hits.push_back(DeDxHit(hitDeDxInfo.charge(i) / hitDeDxInfo.pathlength(i) * theMeVperADCStrip_, 0, 0, 0));
         }
       }
 
       // In case we have a pixel only track, but usePixelForPrescales_ is false
       if (hits.empty()) {
         indices.push_back(-1);
         continue;
       }
       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");
 }

◆ xorshift128p()

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

inlineprivate

Definition at line 88 of file DeDxHitInfoProducer.cc.

References x.

Referenced by produce().

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

Member Data Documentation

◆ calibGains_

std::vector<std::vector<float> > DeDxHitInfoProducer::calibGains_

private

Definition at line 85 of file DeDxHitInfoProducer.cc.

Referenced by beginRun(), and processHit().

◆ calibrationPath_

const std::string DeDxHitInfoProducer::calibrationPath_

private

Definition at line 72 of file DeDxHitInfoProducer.cc.

Referenced by beginRun().

◆ doShapeTest_

const bool DeDxHitInfoProducer::doShapeTest_

private

Definition at line 74 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ lowPtTracksDeDxThreshold_

const float DeDxHitInfoProducer::lowPtTracksDeDxThreshold_

private

Definition at line 78 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ lowPtTracksEstimator_

GenericTruncatedAverageDeDxEstimator DeDxHitInfoProducer::lowPtTracksEstimator_

private

Definition at line 77 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ lowPtTracksPrescaleFail_

const unsigned int DeDxHitInfoProducer::lowPtTracksPrescaleFail_

private

Definition at line 76 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ lowPtTracksPrescalePass_

const unsigned int DeDxHitInfoProducer::lowPtTracksPrescalePass_

private

Definition at line 76 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ maxTrackEta_

const float DeDxHitInfoProducer::maxTrackEta_

private

Definition at line 70 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ minTrackHits_

const unsigned int DeDxHitInfoProducer::minTrackHits_

private

Definition at line 67 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ minTrackPt_

const float DeDxHitInfoProducer::minTrackPt_

private

Definition at line 68 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ minTrackPtPrescale_

const float DeDxHitInfoProducer::minTrackPtPrescale_

private

Definition at line 69 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ offsetDU_

unsigned int DeDxHitInfoProducer::offsetDU_

private

Definition at line 86 of file DeDxHitInfoProducer.cc.

Referenced by beginRun(), and processHit().

◆ theMeVperADCPixel_

const float DeDxHitInfoProducer::theMeVperADCPixel_

private

Definition at line 64 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ theMeVperADCStrip_

const float DeDxHitInfoProducer::theMeVperADCStrip_

private

Definition at line 65 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ tkGeom_

edm::ESHandle<TrackerGeometry> DeDxHitInfoProducer::tkGeom_

private

Definition at line 82 of file DeDxHitInfoProducer.cc.

Referenced by beginRun(), and processHit().

◆ tkGeomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> DeDxHitInfoProducer::tkGeomToken_

private

Definition at line 83 of file DeDxHitInfoProducer.cc.

Referenced by beginRun().

◆ tracksToken_

const edm::EDGetTokenT<reco::TrackCollection> DeDxHitInfoProducer::tracksToken_

private

Definition at line 81 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ useCalibration_

const bool DeDxHitInfoProducer::useCalibration_

private

Definition at line 73 of file DeDxHitInfoProducer.cc.

Referenced by beginRun().

◆ usePixel_

const bool DeDxHitInfoProducer::usePixel_

private

Definition at line 62 of file DeDxHitInfoProducer.cc.

Referenced by DeDxHitInfoProducer(), and processHit().

◆ usePixelForPrescales_

const bool DeDxHitInfoProducer::usePixelForPrescales_

private

Definition at line 79 of file DeDxHitInfoProducer.cc.

Referenced by produce().

◆ useStrip_

const bool DeDxHitInfoProducer::useStrip_

private

Definition at line 63 of file DeDxHitInfoProducer.cc.

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_

◆ calibrationPath_

◆ doShapeTest_

◆ lowPtTracksDeDxThreshold_

◆ lowPtTracksEstimator_

◆ lowPtTracksPrescaleFail_

◆ lowPtTracksPrescalePass_

◆ maxTrackEta_

◆ minTrackHits_

◆ minTrackPt_

◆ minTrackPtPrescale_

◆ offsetDU_

◆ theMeVperADCPixel_

◆ theMeVperADCStrip_

◆ tkGeom_

◆ tkGeomToken_

◆ tracksToken_

◆ useCalibration_

◆ usePixel_

◆ usePixelForPrescales_

◆ useStrip_