HSCPDeDxInfoProducer Class Reference

#include <RecoTracker/HSCPDeDxInfoProducer/src/HSCPDeDxInfoProducer.cc>

Inheritance diagram for HSCPDeDxInfoProducer:

Public Member Functions
	HSCPDeDxInfoProducer (const edm::ParameterSet &)
	~HSCPDeDxInfoProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
void	beginRun (edm::Run const &run, const edm::EventSetup &) override
void	makeCalibrationMap (const TrackerGeometry &tkGeom)
void	processHit (const TrackingRecHit *recHit, float trackMomentum, float &cosine, susybsm::HSCPDeDxInfo &hscpDeDxInfo, LocalPoint HitLocalPos)
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::vector< std::vector< float > >	calibGains
std::string	m_calibrationPath
unsigned int	m_off
edm::EDGetTokenT< reco::TrackCollection >	m_tracksTag
edm::EDGetTokenT< TrajTrackAssociationCollection >	m_trajTrackAssociationTag
unsigned int	MaxNrStrips
float	MaxTrackChiOverNdf
float	MaxTrackEta
float	MaxTrackMomentum
float	meVperADCPixel
float	meVperADCStrip
float	MinTrackEta
unsigned int	MinTrackHits
float	MinTrackMomentum
TH3F *	Prob_ChargePath
std::string	ProbabilityMode
std::string	Reccord
bool	shapetest
bool	useCalibration
bool	usePixel
bool	useStrip
bool	useTrajectory

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

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

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

Definition at line 33 of file HSCPDeDxInfoProducer.h.

Constructor & Destructor Documentation

HSCPDeDxInfoProducer()

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

explicit

Definition at line 29 of file HSCPDeDxInfoProducer.cc.

                                                                         {
  produces<ValueMap<susybsm::HSCPDeDxInfo> >();
 
  MaxNrStrips = iConfig.getUntrackedParameter<unsigned>("maxNrStrips", 255);
  MinTrackHits = iConfig.getUntrackedParameter<unsigned>("MinTrackHits", 4);
  MinTrackMomentum = iConfig.getUntrackedParameter<double>("minTrackMomentum", 0.0);
  MaxTrackMomentum = iConfig.getUntrackedParameter<double>("maxTrackMomentum", 99999.0);
  MinTrackEta = iConfig.getUntrackedParameter<double>("minTrackEta", -5.0);
  MaxTrackEta = iConfig.getUntrackedParameter<double>("maxTrackEta", 5.0);
 
  m_tracksTag = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("tracks"));
  m_trajTrackAssociationTag =
      consumes<TrajTrackAssociationCollection>(iConfig.getParameter<edm::InputTag>("trajectoryTrackAssociation"));
  useTrajectory = iConfig.getUntrackedParameter<bool>("UseTrajectory", true);
 
  usePixel = iConfig.getParameter<bool>("UsePixel");
  useStrip = iConfig.getParameter<bool>("UseStrip");
  meVperADCPixel = iConfig.getParameter<double>("MeVperADCPixel");
  meVperADCStrip = iConfig.getParameter<double>("MeVperADCStrip");
 
  shapetest = iConfig.getParameter<bool>("ShapeTest");
  useCalibration = iConfig.getParameter<bool>("UseCalibration");
  m_calibrationPath = iConfig.getParameter<string>("calibrationPath");
 
  //   Reccord             = iConfig.getUntrackedParameter<std::string>  ("Reccord"            , "SiStripDeDxMip_3D_Rcd");
  //   ProbabilityMode     = iConfig.getUntrackedParameter<std::string>  ("ProbabilityMode"    , "Accumulation");
  //   Prob_ChargePath     = NULL;
 
  if (!usePixel && !useStrip)
    edm::LogWarning("DeDxHitsProducer")
        << "Pixel Hits AND Strip Hits will not be used to estimate dEdx --> BUG, Please Update the config file";
}

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_FULL_cff::MaxTrackEta, HLT_FULL_cff::MinTrackHits, SiStripFineDelayHit_cfi::MinTrackMomentum, dedxEstimators_cff::useCalibration, dedxEstimators_cff::usePixel, dedxEstimators_cff::useStrip, and conversionTrackProducer_cfi::useTrajectory.

~HSCPDeDxInfoProducer()

HSCPDeDxInfoProducer::~HSCPDeDxInfoProducer ( )

override

Definition at line 62 of file HSCPDeDxInfoProducer.cc.

{}

Member Function Documentation

beginRun()

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

overrideprivate

Definition at line 65 of file HSCPDeDxInfoProducer.cc.

                                                                                  {
  if (useCalibration && calibGains.empty()) {
    edm::ESHandle<TrackerGeometry> tkGeom;
    iSetup.get<TrackerDigiGeometryRecord>().get(tkGeom);
    m_off = tkGeom->offsetDU(GeomDetEnumerators::PixelBarrel);  //index start at the first pixel
 
    DeDxTools::makeCalibrationMap(m_calibrationPath, *tkGeom, calibGains, m_off);
  }
 
  //   DeDxTools::buildDiscrimMap(run, iSetup, Reccord,  ProbabilityMode, Prob_ChargePath);
}

References reco::get(), edm::EventSetup::get(), DeDxTools::makeCalibrationMap(), TrackerGeometry::offsetDU(), GeomDetEnumerators::PixelBarrel, and dedxEstimators_cff::useCalibration.

makeCalibrationMap()

void HSCPDeDxInfoProducer::makeCalibrationMap ( const TrackerGeometry &  tkGeom )

private

processHit()

void HSCPDeDxInfoProducer::processHit ( const TrackingRecHit *  recHit, float  trackMomentum, float &  cosine, susybsm::HSCPDeDxInfo &  hscpDeDxInfo, LocalPoint  HitLocalPos  )

private

Definition at line 141 of file HSCPDeDxInfoProducer.cc.

                                                              {
  auto const& thit = static_cast<BaseTrackerRecHit const&>(*recHit);
  if (!thit.isValid())
    return;
 
  auto const& clus = thit.firstClusterRef();
  if (!clus.isValid())
    return;
 
  if (clus.isPixel()) {
    if (!usePixel)
      return;
 
    auto& detUnit = *(recHit->detUnit());
    float pathLen = detUnit.surface().bounds().thickness() / fabs(cosine);
    float chargeAbs = clus.pixelCluster().charge();
    hscpDeDxInfo.charges.push_back(chargeAbs);
    hscpDeDxInfo.pathlengths.push_back(pathLen);
    hscpDeDxInfo.detIds.push_back(thit.geographicalId());
    hscpDeDxInfo.localPosXs.push_back(HitLocalPos.x());
    hscpDeDxInfo.localPosYs.push_back(HitLocalPos.y());
    hscpDeDxInfo.clusterIndices.push_back(clus.key());
  } else if (clus.isStrip() && !thit.isMatched()) {
    if (!useStrip)
      return;
 
    auto& detUnit = *(recHit->detUnit());
    int NSaturating = 0;
    float pathLen = detUnit.surface().bounds().thickness() / fabs(cosine);
    float chargeAbs = DeDxTools::getCharge(&(clus.stripCluster()), NSaturating, detUnit, calibGains, m_off);
    hscpDeDxInfo.charges.push_back(chargeAbs);
    hscpDeDxInfo.pathlengths.push_back(pathLen);
    hscpDeDxInfo.detIds.push_back(thit.geographicalId());
    hscpDeDxInfo.localPosXs.push_back(HitLocalPos.x());
    hscpDeDxInfo.localPosYs.push_back(HitLocalPos.y());
    hscpDeDxInfo.clusterIndices.push_back(clus.key());
  } else if (clus.isStrip() && thit.isMatched()) {
    if (!useStrip)
      return;
    const SiStripMatchedRecHit2D* matchedHit = dynamic_cast<const SiStripMatchedRecHit2D*>(recHit);
    if (!matchedHit)
      return;
 
    auto& detUnitM = *(matchedHit->monoHit().detUnit());
    int NSaturating = 0;
    float pathLen = detUnitM.surface().bounds().thickness() / fabs(cosine);
    float chargeAbs =
        DeDxTools::getCharge(&(matchedHit->monoHit().stripCluster()), NSaturating, detUnitM, calibGains, m_off);
    hscpDeDxInfo.charges.push_back(chargeAbs);
    hscpDeDxInfo.pathlengths.push_back(pathLen);
    hscpDeDxInfo.detIds.push_back(thit.geographicalId());
    hscpDeDxInfo.localPosXs.push_back(HitLocalPos.x());
    hscpDeDxInfo.localPosYs.push_back(HitLocalPos.y());
    const OmniClusterRef monoClusterRef = matchedHit->monoClusterRef();
    hscpDeDxInfo.clusterIndices.push_back(monoClusterRef.key());
 
    auto& detUnitS = *(matchedHit->stereoHit().detUnit());
    NSaturating = 0;
    pathLen = detUnitS.surface().bounds().thickness() / fabs(cosine);
    chargeAbs =
        DeDxTools::getCharge(&(matchedHit->stereoHit().stripCluster()), NSaturating, detUnitS, calibGains, m_off);
    hscpDeDxInfo.charges.push_back(chargeAbs);
    hscpDeDxInfo.pathlengths.push_back(pathLen);
    hscpDeDxInfo.detIds.push_back(thit.geographicalId());
    hscpDeDxInfo.localPosXs.push_back(HitLocalPos.x());
    hscpDeDxInfo.localPosYs.push_back(HitLocalPos.y());
    const OmniClusterRef stereoClusterRef = matchedHit->stereoClusterRef();
    hscpDeDxInfo.clusterIndices.push_back(stereoClusterRef.key());
  }
}

References susybsm::HSCPDeDxInfo::charges, susybsm::HSCPDeDxInfo::clusterIndices, susybsm::HSCPDeDxInfo::detIds, TrackerSingleRecHit::detUnit(), DeDxTools::getCharge(), OmniClusterRef::key(), susybsm::HSCPDeDxInfo::localPosXs, susybsm::HSCPDeDxInfo::localPosYs, SiStripMatchedRecHit2D::monoClusterRef(), SiStripMatchedRecHit2D::monoHit(), susybsm::HSCPDeDxInfo::pathlengths, rpcPointValidation_cfi::recHit, SiStripMatchedRecHit2D::stereoClusterRef(), SiStripMatchedRecHit2D::stereoHit(), TrackerSingleRecHit::stripCluster(), dedxEstimators_cff::usePixel, dedxEstimators_cff::useStrip, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

produce()

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

overrideprivate

Definition at line 77 of file HSCPDeDxInfoProducer.cc.

                                                                                {
  unique_ptr<ValueMap<susybsm::HSCPDeDxInfo> > trackDeDxAssociation(new ValueMap<susybsm::HSCPDeDxInfo>);
  ValueMap<susybsm::HSCPDeDxInfo>::Filler filler(*trackDeDxAssociation);
 
  edm::Handle<reco::TrackCollection> trackCollectionHandle;
  iEvent.getByToken(m_tracksTag, trackCollectionHandle);
 
  Handle<TrajTrackAssociationCollection> trajTrackAssociationHandle;
  if (useTrajectory)
    iEvent.getByToken(m_trajTrackAssociationTag, trajTrackAssociationHandle);
 
  std::vector<susybsm::HSCPDeDxInfo> dEdxInfos(trackCollectionHandle->size());
 
  TrajTrackAssociationCollection::const_iterator cit;
  if (useTrajectory)
    cit = trajTrackAssociationHandle->begin();
  for (unsigned int j = 0; j < trackCollectionHandle->size(); j++) {
    const reco::TrackRef track = reco::TrackRef(trackCollectionHandle.product(), j);
 
    susybsm::HSCPDeDxInfo hscpDeDxInfo;
 
    if (useTrajectory) {  //trajectory allows to take into account the local direction of the particle on the module sensor --> muc much better 'dx' measurement
      const edm::Ref<std::vector<Trajectory> > traj = cit->key;
      cit++;
      const vector<TrajectoryMeasurement>& measurements = traj->measurements();
      for (vector<TrajectoryMeasurement>::const_iterator it = measurements.begin(); it != measurements.end(); it++) {
        TrajectoryStateOnSurface trajState = it->updatedState();
        if (!trajState.isValid())
          continue;
 
        const TrackingRecHit* recHit = (*it->recHit()).hit();
        if (!recHit)
          continue;
        LocalVector trackDirection = trajState.localDirection();
        float cosine = trackDirection.z() / trackDirection.mag();
 
        processHit(recHit, trajState.localMomentum().mag(), cosine, hscpDeDxInfo, trajState.localPosition());
      }
 
    } else {  //assume that the particles trajectory is a straight line originating from the center of the detector  (can be improved)
      for (unsigned int h = 0; h < track->recHitsSize(); h++) {
        const TrackingRecHit* recHit = &(*(track->recHit(h)));
        auto const& thit = static_cast<BaseTrackerRecHit const&>(*recHit);
        if (!thit.isValid())
          continue;  //make sure it's a tracker hit
 
        const GlobalVector& ModuleNormal = recHit->detUnit()->surface().normalVector();
        float cosine =
            (track->px() * ModuleNormal.x() + track->py() * ModuleNormal.y() + track->pz() * ModuleNormal.z()) /
            track->p();
 
        processHit(recHit, track->p(), cosine, hscpDeDxInfo, LocalPoint(0.0, 0.0));
      }
    }
 
    dEdxInfos[j] = hscpDeDxInfo;
  }
 
  filler.insert(trackCollectionHandle, dEdxInfos.begin(), dEdxInfos.end());
  filler.fill();
  iEvent.put(std::move(trackDeDxAssociation));
}

References edm::AssociationMap< Tag >::begin(), edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, trigObjTnPSource_cfi::filler, iEvent, TrajectoryStateOnSurface::isValid(), dqmiolumiharvest::j, edm::Ref< C, T, F >::key(), TrajectoryStateOnSurface::localDirection(), TrajectoryStateOnSurface::localMomentum(), TrajectoryStateOnSurface::localPosition(), PV3DBase< T, PVType, FrameType >::mag(), eostools::move(), edm::Handle< T >::product(), rpcPointValidation_cfi::recHit, HLT_FULL_cff::track, conversionTrackProducer_cfi::useTrajectory, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Member Data Documentation

calibGains

std::vector<std::vector<float> > HSCPDeDxInfoProducer::calibGains

private

Definition at line 72 of file HSCPDeDxInfoProducer.h.

m_calibrationPath

std::string HSCPDeDxInfoProducer::m_calibrationPath

private

Definition at line 68 of file HSCPDeDxInfoProducer.h.

m_off

unsigned int HSCPDeDxInfoProducer::m_off

private

Definition at line 73 of file HSCPDeDxInfoProducer.h.

m_tracksTag

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

private

Definition at line 52 of file HSCPDeDxInfoProducer.h.

m_trajTrackAssociationTag

edm::EDGetTokenT<TrajTrackAssociationCollection> HSCPDeDxInfoProducer::m_trajTrackAssociationTag

private

Definition at line 51 of file HSCPDeDxInfoProducer.h.

MaxNrStrips

unsigned int HSCPDeDxInfoProducer::MaxNrStrips

private

Definition at line 60 of file HSCPDeDxInfoProducer.h.

MaxTrackChiOverNdf

float HSCPDeDxInfoProducer::MaxTrackChiOverNdf

private

Definition at line 66 of file HSCPDeDxInfoProducer.h.

MaxTrackEta

float HSCPDeDxInfoProducer::MaxTrackEta

private

Definition at line 65 of file HSCPDeDxInfoProducer.h.

MaxTrackMomentum

float HSCPDeDxInfoProducer::MaxTrackMomentum

private

Definition at line 63 of file HSCPDeDxInfoProducer.h.

meVperADCPixel

float HSCPDeDxInfoProducer::meVperADCPixel

private

Definition at line 57 of file HSCPDeDxInfoProducer.h.

meVperADCStrip

float HSCPDeDxInfoProducer::meVperADCStrip

private

Definition at line 58 of file HSCPDeDxInfoProducer.h.

MinTrackEta

float HSCPDeDxInfoProducer::MinTrackEta

private

Definition at line 64 of file HSCPDeDxInfoProducer.h.

MinTrackHits

unsigned int HSCPDeDxInfoProducer::MinTrackHits

private

Definition at line 61 of file HSCPDeDxInfoProducer.h.

MinTrackMomentum

float HSCPDeDxInfoProducer::MinTrackMomentum

private

Definition at line 62 of file HSCPDeDxInfoProducer.h.

Prob_ChargePath

TH3F* HSCPDeDxInfoProducer::Prob_ChargePath

private

Definition at line 77 of file HSCPDeDxInfoProducer.h.

ProbabilityMode

std::string HSCPDeDxInfoProducer::ProbabilityMode

private

Definition at line 76 of file HSCPDeDxInfoProducer.h.

Reccord

std::string HSCPDeDxInfoProducer::Reccord

private

Definition at line 75 of file HSCPDeDxInfoProducer.h.

shapetest

bool HSCPDeDxInfoProducer::shapetest

private

Definition at line 70 of file HSCPDeDxInfoProducer.h.

useCalibration

bool HSCPDeDxInfoProducer::useCalibration

private

Definition at line 69 of file HSCPDeDxInfoProducer.h.

usePixel

bool HSCPDeDxInfoProducer::usePixel

private

Definition at line 55 of file HSCPDeDxInfoProducer.h.

useStrip

bool HSCPDeDxInfoProducer::useStrip

private

Definition at line 56 of file HSCPDeDxInfoProducer.h.

useTrajectory

bool HSCPDeDxInfoProducer::useTrajectory

private

Definition at line 54 of file HSCPDeDxInfoProducer.h.