#include <VectorHitBuilderAlgorithmBase.h>

Inheritance diagram for VectorHitBuilderAlgorithmBase:

Public Types
typedef Detset::const_iterator	const_iterator

typedef edmNew::DetSet < Phase2TrackerCluster1D >	Detset

typedef edmNew::DetSetVector < VectorHit >	output_t

typedef edm::Ref < edmNew::DetSetVector < Phase2TrackerCluster1D > , Phase2TrackerCluster1D >	Phase2TrackerCluster1DRef

typedef std::pair < StackGeomDet, std::vector < Phase2TrackerCluster1D > >	StackClusters

Public Member Functions
virtual VectorHit	buildVectorHit (const StackGeomDet *stack, Phase2TrackerCluster1DRef lower, Phase2TrackerCluster1DRef upper) const =0

virtual void	buildVectorHits (VectorHitCollection &vhAcc, VectorHitCollection &vhRej, DetId detIdStack, const StackGeomDet *stack, edm::Handle< edmNew::DetSetVector< Phase2TrackerCluster1D >> clusters, const Detset &DSVinner, const Detset &DSVouter, const std::vector< bool > &phase2OTClustersToSkip=std::vector< bool >()) const =0

double	computeParallaxCorrection (const PixelGeomDetUnit , const Point3DBase< float, LocalTag > &, const PixelGeomDetUnit , const Point3DBase< float, LocalTag > &) const

void	printCluster (const GeomDet geomDetUnit, const Phase2TrackerCluster1D cluster) const

void	printClusters (const edmNew::DetSetVector< Phase2TrackerCluster1D > &clusters) const

virtual void	run (edm::Handle< edmNew::DetSetVector< Phase2TrackerCluster1D >> clusters, VectorHitCollection &vhAcc, VectorHitCollection &vhRej, edmNew::DetSetVector< Phase2TrackerCluster1D > &clustersAcc, edmNew::DetSetVector< Phase2TrackerCluster1D > &clustersRej) const =0

	VectorHitBuilderAlgorithmBase (const edm::ParameterSet &, const TrackerGeometry , const TrackerTopology , const ClusterParameterEstimator< Phase2TrackerCluster1D > *)

virtual	~VectorHitBuilderAlgorithmBase ()

Public Attributes
std::vector< double >	barrelCut_

const ClusterParameterEstimator < Phase2TrackerCluster1D > *	cpe_

std::vector< double >	endcapCut_

unsigned int	nMaxVHforeachStack_

const TrackerGeometry *	tkGeom_

const TrackerTopology *	tkTopo_

Private Attributes
edm::ESInputTag	cpeTag_

Detailed Description

Definition at line 18 of file VectorHitBuilderAlgorithmBase.h.

Member Typedef Documentation

typedef Detset::const_iterator VectorHitBuilderAlgorithmBase::const_iterator

Definition at line 22 of file VectorHitBuilderAlgorithmBase.h.

typedef edmNew::DetSet<Phase2TrackerCluster1D> VectorHitBuilderAlgorithmBase::Detset

Definition at line 21 of file VectorHitBuilderAlgorithmBase.h.

typedef edmNew::DetSetVector<VectorHit> VectorHitBuilderAlgorithmBase::output_t

Definition at line 23 of file VectorHitBuilderAlgorithmBase.h.

typedef edm::Ref<edmNew::DetSetVector<Phase2TrackerCluster1D>, Phase2TrackerCluster1D> VectorHitBuilderAlgorithmBase::Phase2TrackerCluster1DRef

Definition at line 20 of file VectorHitBuilderAlgorithmBase.h.

typedef std::pair<StackGeomDet, std::vector<Phase2TrackerCluster1D> > VectorHitBuilderAlgorithmBase::StackClusters

Definition at line 24 of file VectorHitBuilderAlgorithmBase.h.

Constructor & Destructor Documentation

VectorHitBuilderAlgorithmBase::VectorHitBuilderAlgorithmBase	(	const edm::ParameterSet &	conf,
		const TrackerGeometry *	tkGeomProd,
		const TrackerTopology *	tkTopoProd,
		const ClusterParameterEstimator< Phase2TrackerCluster1D > *	cpeProd
	)

Definition at line 12 of file VectorHitBuilderAlgorithmBase.cc.

     : tkGeom_(tkGeomProd),
       tkTopo_(tkTopoProd),
       cpe_(cpeProd),
       nMaxVHforeachStack_(conf.getParameter<int>("maxVectorHitsInAStack")),
       barrelCut_(conf.getParameter<std::vector<double> >("BarrelCut")),
       endcapCut_(conf.getParameter<std::vector<double> >("EndcapCut")),
       cpeTag_(conf.getParameter<edm::ESInputTag>("CPE")) {}

virtual VectorHitBuilderAlgorithmBase::~VectorHitBuilderAlgorithmBase ( )

inlinevirtual

Definition at line 30 of file VectorHitBuilderAlgorithmBase.h.

30 {}

Member Function Documentation

virtual VectorHit VectorHitBuilderAlgorithmBase::buildVectorHit	(	const StackGeomDet *	stack,
		Phase2TrackerCluster1DRef	lower,
		Phase2TrackerCluster1DRef	upper
	)		const

pure virtual

Implemented in VectorHitBuilderAlgorithm.

virtual void VectorHitBuilderAlgorithmBase::buildVectorHits	(	VectorHitCollection &	vhAcc,
		VectorHitCollection &	vhRej,
		DetId	detIdStack,
		const StackGeomDet *	stack,
		edm::Handle< edmNew::DetSetVector< Phase2TrackerCluster1D >>	clusters,
		const Detset &	DSVinner,
		const Detset &	DSVouter,
		const std::vector< bool > &	phase2OTClustersToSkip = `std::vector< bool >()`
	)		const

pure virtual

Implemented in VectorHitBuilderAlgorithm.

double VectorHitBuilderAlgorithmBase::computeParallaxCorrection	(	const PixelGeomDetUnit *	geomDetUnit_low,
		const Point3DBase< float, LocalTag > &	lPosClu_low,
		const PixelGeomDetUnit *	geomDetUnit_upp,
		const Point3DBase< float, LocalTag > &	lPosClu_upp
	)		const

Definition at line 25 of file VectorHitBuilderAlgorithmBase.cc.

References LogTrace, GeomDet::surface(), Surface::toGlobal(), GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by VectorHitBuilderAlgorithm::buildVectorHits().

                                                                                                                        {
   double parallCorr = 0.0;
   Global3DPoint origin(0, 0, 0);
   Global3DPoint gPosClu_low = geomDetUnit_low->surface().toGlobal(lPosClu_low);
   GlobalVector gV = gPosClu_low - origin;
   LogTrace("VectorHitsBuilderValidation") << " global vector passing to the origin:" << gV;
 
   LocalVector lV = geomDetUnit_low->surface().toLocal(gV);
   LogTrace("VectorHitsBuilderValidation")
       << " local vector passing to the origin (in the lower detector system of reference):" << lV;
   LocalVector lV_norm = lV / lV.z();
   LogTrace("VectorHitsBuilderValidation")
       << " normalized local vector passing to the origin (in low the lower detector system of reference):" << lV_norm;
 
   Global3DPoint gPosClu_upp = geomDetUnit_upp->surface().toGlobal(lPosClu_upp);
   Local3DPoint lPosClu_uppInLow = geomDetUnit_low->surface().toLocal(gPosClu_upp);
   parallCorr = lV_norm.x() * lPosClu_uppInLow.z();
 
   return parallCorr;
 }

void VectorHitBuilderAlgorithmBase::printCluster	(	const GeomDet *	geomDetUnit,
		const Phase2TrackerCluster1D *	cluster
	)		const

Definition at line 67 of file VectorHitBuilderAlgorithmBase.cc.

References Surface::bounds(), cpe_, GeomDet::geographicalId(), TrackerGeometry::getDetectorType(), phase1PixelTopology::layer, TrackerTopology::layer(), Bounds::length(), ClusterParameterEstimator< T >::localParameters(), LogTrace, callgraph::module, TrackerTopology::module(), TrackerGeometry::Ph2PSP, TrackerGeometry::Ph2PSS, TrackerGeometry::Ph2SS, PixelTopology::pitch(), DetId::rawId(), PixelGeomDetUnit::specificTopology(), GeomDet::surface(), tkGeom_, tkTopo_, Surface::toGlobal(), and Bounds::width().

Referenced by VectorHitBuilderAlgorithm::buildVectorHit(), VectorHitBuilderAlgorithm::buildVectorHits(), and printClusters().

                                                                                               {
   if (!geomDetUnit)
     return;
   const PixelGeomDetUnit* pixelGeomDetUnit = dynamic_cast<const PixelGeomDetUnit*>(geomDetUnit);
   const PixelTopology& topol = pixelGeomDetUnit->specificTopology();
   if (!pixelGeomDetUnit)
     return;
 
   unsigned int layer = tkTopo_->layer(geomDetUnit->geographicalId());
   unsigned int module = tkTopo_->module(geomDetUnit->geographicalId());
   LogTrace("VectorHitBuilder") << "Layer:" << layer << " and DetId: " << geomDetUnit->geographicalId().rawId()
                                << std::endl;
   TrackerGeometry::ModuleType mType = tkGeom_->getDetectorType(geomDetUnit->geographicalId());
   if (mType == TrackerGeometry::ModuleType::Ph2PSP)
     LogTrace("VectorHitBuilder") << "Pixel cluster (module:" << module << ") " << std::endl;
   else if (mType == TrackerGeometry::ModuleType::Ph2SS || mType == TrackerGeometry::ModuleType::Ph2PSS)
     LogTrace("VectorHitBuilder") << "Strip cluster (module:" << module << ") " << std::endl;
   else
     LogTrace("VectorHitBuilder") << "no module?!" << std::endl;
   LogTrace("VectorHitBuilder") << "with pitch:" << topol.pitch().first << " , " << topol.pitch().second << std::endl;
   LogTrace("VectorHitBuilder") << " and width:" << pixelGeomDetUnit->surface().bounds().width()
                                << " , lenght:" << pixelGeomDetUnit->surface().bounds().length() << std::endl;
 
   auto&& lparams = cpe_->localParameters(*clustIt, *pixelGeomDetUnit);
   Global3DPoint gparams = pixelGeomDetUnit->surface().toGlobal(lparams.first);
 
   LogTrace("VectorHitBuilder") << "\t global pos " << gparams << std::endl;
   LogTrace("VectorHitBuilder") << "\t local  pos " << lparams.first << "with err " << lparams.second << std::endl;
   LogTrace("VectorHitBuilder") << std::endl;
 
   return;
 }

void VectorHitBuilderAlgorithmBase::printClusters ( const edmNew::DetSetVector< Phase2TrackerCluster1D > & clusters ) const

Definition at line 49 of file VectorHitBuilderAlgorithmBase.cc.

References TrackerGeometry::idToDetUnit(), LogDebug, printCluster(), and tkGeom_.

Referenced by VectorHitBuilderEDProducer::produce().

                                                                                                                   {
   int nCluster = 0;
   int numberOfDSV = 0;
   for (const auto& DSViter : clusters) {
     ++numberOfDSV;
     // Loop over the clusters in the detector unit
     for (const auto& clustIt : DSViter) {
       nCluster++;
       // get the detector unit's id
       const GeomDetUnit* geomDetUnit(tkGeom_->idToDetUnit(DSViter.detId()));
       if (!geomDetUnit)
         return;
       printCluster(geomDetUnit, &clustIt);
     }
   }
   LogDebug("VectorHitBuilder") << " Number of input clusters: " << nCluster << std::endl;
 }

virtual void VectorHitBuilderAlgorithmBase::run	(	edm::Handle< edmNew::DetSetVector< Phase2TrackerCluster1D >>	clusters,
		VectorHitCollection &	vhAcc,
		VectorHitCollection &	vhRej,
		edmNew::DetSetVector< Phase2TrackerCluster1D > &	clustersAcc,
		edmNew::DetSetVector< Phase2TrackerCluster1D > &	clustersRej
	)		const