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

Inheritance diagram for HSCPDeDxInfoProducer:

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

	~HSCPDeDxInfoProducer ()

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
virtual 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)

virtual 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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

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

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

Definition at line 33 of file HSCPDeDxInfoProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 30 of file HSCPDeDxInfoProducer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), QuarkoniaTrackSelector_cfi::MaxTrackEta, dedxEstimators_cff::useCalibration, dedxEstimators_cff::usePixel, dedxEstimators_cff::useStrip, and conversionTrackProducer_cfi::useTrajectory.

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

HSCPDeDxInfoProducer::~HSCPDeDxInfoProducer ( )

Definition at line 65 of file HSCPDeDxInfoProducer.cc.

65 {}

Member Function Documentation

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

overrideprivatevirtual

Definition at line 68 of file HSCPDeDxInfoProducer.cc.

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

 {
    if(useCalibration && calibGains.size()==0){
       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);
 }

void HSCPDeDxInfoProducer::makeCalibrationMap ( const TrackerGeometry & tkGeom )

private

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

private

Definition at line 140 of file HSCPDeDxInfoProducer.cc.

References Surface::bounds(), susybsm::HSCPDeDxInfo::charges, susybsm::HSCPDeDxInfo::clusterIndices, DEFINE_FWK_MODULE, susybsm::HSCPDeDxInfo::detIds, TrackerSingleRecHit::detUnit(), TrackingRecHit::detUnit(), BaseTrackerRecHit::firstClusterRef(), 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(), BaseTrackerRecHit::surface(), Bounds::thickness(), dedxEstimators_cff::usePixel, dedxEstimators_cff::useStrip, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

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

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

overrideprivatevirtual

Definition at line 83 of file HSCPDeDxInfoProducer.cc.

References edm::AssociationMap< Tag >::begin(), edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, TrackingRecHit::detUnit(), edm::helper::Filler< Map >::fill(), objects.autophobj::filler, edm::Event::getByToken(), edm::helper::Filler< Map >::insert(), TrajectoryStateOnSurface::isValid(), edm::Ref< C, T, F >::key(), TrajectoryStateOnSurface::localDirection(), TrajectoryStateOnSurface::localMomentum(), TrajectoryStateOnSurface::localPosition(), PV3DBase< T, PVType, FrameType >::mag(), eostools::move(), Plane::normalVector(), edm::Handle< T >::product(), edm::Event::put(), rpcPointValidation_cfi::recHit, GeomDet::surface(), HiIsolationCommonParameters_cff::track, conversionTrackProducer_cfi::useTrajectory, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

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

Member Data Documentation

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

private

Definition at line 68 of file HSCPDeDxInfoProducer.h.

std::string HSCPDeDxInfoProducer::m_calibrationPath

private

Definition at line 64 of file HSCPDeDxInfoProducer.h.

unsigned int HSCPDeDxInfoProducer::m_off

private

Definition at line 69 of file HSCPDeDxInfoProducer.h.

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

private

Definition at line 48 of file HSCPDeDxInfoProducer.h.

edm::EDGetTokenT<TrajTrackAssociationCollection> HSCPDeDxInfoProducer::m_trajTrackAssociationTag

private

Definition at line 47 of file HSCPDeDxInfoProducer.h.

unsigned int HSCPDeDxInfoProducer::MaxNrStrips

private

Definition at line 56 of file HSCPDeDxInfoProducer.h.

float HSCPDeDxInfoProducer::MaxTrackChiOverNdf

private

Definition at line 62 of file HSCPDeDxInfoProducer.h.

float HSCPDeDxInfoProducer::MaxTrackEta

private

Definition at line 61 of file HSCPDeDxInfoProducer.h.

float HSCPDeDxInfoProducer::MaxTrackMomentum

private

Definition at line 59 of file HSCPDeDxInfoProducer.h.

float HSCPDeDxInfoProducer::meVperADCPixel

private

Definition at line 53 of file HSCPDeDxInfoProducer.h.

float HSCPDeDxInfoProducer::meVperADCStrip

private

Definition at line 54 of file HSCPDeDxInfoProducer.h.

float HSCPDeDxInfoProducer::MinTrackEta

private

Definition at line 60 of file HSCPDeDxInfoProducer.h.

unsigned int HSCPDeDxInfoProducer::MinTrackHits

private

Definition at line 57 of file HSCPDeDxInfoProducer.h.

float HSCPDeDxInfoProducer::MinTrackMomentum

private

Definition at line 58 of file HSCPDeDxInfoProducer.h.

TH3F* HSCPDeDxInfoProducer::Prob_ChargePath

private

Definition at line 73 of file HSCPDeDxInfoProducer.h.

std::string HSCPDeDxInfoProducer::ProbabilityMode

private

Definition at line 72 of file HSCPDeDxInfoProducer.h.

std::string HSCPDeDxInfoProducer::Reccord

private

Definition at line 71 of file HSCPDeDxInfoProducer.h.

bool HSCPDeDxInfoProducer::shapetest

private

Definition at line 66 of file HSCPDeDxInfoProducer.h.

bool HSCPDeDxInfoProducer::useCalibration

private

Definition at line 65 of file HSCPDeDxInfoProducer.h.

bool HSCPDeDxInfoProducer::usePixel

private

Definition at line 51 of file HSCPDeDxInfoProducer.h.

bool HSCPDeDxInfoProducer::useStrip

private

Definition at line 52 of file HSCPDeDxInfoProducer.h.

bool HSCPDeDxInfoProducer::useTrajectory

private

Definition at line 50 of file HSCPDeDxInfoProducer.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation