#include <StripCPE.h>

Inheritance diagram for StripCPE:

Classes
struct	AlgoParam

struct	Param

Public Member Functions
LocalVector	driftDirection (const StripGeomDetUnit *det) const

AlgoParam	getAlgoParam (const GeomDetUnit &det, const LocalTrajectoryParameters &ltp) const

StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, const GeomDetUnit &) const override

virtual StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, AlgoParam const &ap) const

	StripCPE (edm::ParameterSet &conf, const MagneticField &, const TrackerGeometry &, const SiStripLorentzAngle &, const SiStripBackPlaneCorrection &, const SiStripConfObject &, const SiStripLatency &)

Public Member Functions inherited from StripClusterParameterEstimator
virtual void	localParameters (AClusters const &clusters, ALocalValues &retValues, const GeomDetUnit &gd, const LocalTrajectoryParameters &ltp) const

virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const LocalTrajectoryParameters &) const

virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const

virtual VLocalValues	localParametersV (const SiStripCluster &cluster, const GeomDetUnit &gd) const

virtual VLocalValues	localParametersV (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const

virtual	~StripClusterParameterEstimator ()

Protected Member Functions
Param const &	param (const GeomDetUnit &det) const

Protected Attributes
const SiStripBackPlaneCorrection &	BackPlaneCorrectionMap_

const TrackerGeometry &	geom_

const SiStripLorentzAngle &	LorentzAngleMap_

const MagneticField &	magfield_

const bool	peakMode_

std::vector< float >	xtalk1

std::vector< float >	xtalk2

Private Types
typedef std::vector< Param >	Params

Private Member Functions
void	fillParams ()

Private Attributes
unsigned int	m_off

Params	m_Params

Additional Inherited Members
Public Types inherited from StripClusterParameterEstimator
using	AClusters = DynArray< SiStripCluster const * >

using	ALocalValues = DynArray< LocalValues >

using	LocalValues = std::pair< LocalPoint, LocalError >

typedef std::vector< LocalValues >	VLocalValues

Detailed Description

Definition at line 16 of file StripCPE.h.

Member Typedef Documentation

typedef std::vector<Param> StripCPE::Params

private

Definition at line 88 of file StripCPE.h.

Constructor & Destructor Documentation

StripCPE::StripCPE	(	edm::ParameterSet &	conf,
		const MagneticField &	mag,
		const TrackerGeometry &	geom,
		const SiStripLorentzAngle &	LorentzAngle,
		const SiStripBackPlaneCorrection &	BackPlaneCorrection,
		const SiStripConfObject &	confObj,
		const SiStripLatency &	latency
	)

Definition at line 11 of file StripCPE.cc.

References fillParams(), SiStripConfObject::get(), SiStripDetId::IB1, SiStripDetId::IB2, SiStripConfObject::isParameter(), edmConvertToStreamModule::modules, SiStripDetId::OB1, SiStripDetId::OB2, peakMode_, edm::second(), findQualityFiles::size, AlCaHLTBitMon_QueryRunRegistry::string, SiStripDetId::W1A, SiStripDetId::W1B, SiStripDetId::W2A, SiStripDetId::W2B, SiStripDetId::W3A, SiStripDetId::W3B, SiStripDetId::W4, SiStripDetId::W5, SiStripDetId::W6, SiStripDetId::W7, xtalk1, and xtalk2.

   : peakMode_(latency.singleReadOutMode() == 1),
     geom_(geom),
     magfield_(mag),
     LorentzAngleMap_(LorentzAngle),
     BackPlaneCorrectionMap_(BackPlaneCorrection)
 {
   typedef std::map<std::string,SiStripDetId::ModuleGeometry> map_t;
   map_t modules;
   modules["IB1"]=SiStripDetId::IB1;  
   modules["IB2"]=SiStripDetId::IB2;
   modules["OB1"]=SiStripDetId::OB1;  
   modules["OB2"]=SiStripDetId::OB2;
   modules["W1A"]=SiStripDetId::W1A;    
   modules["W2A"]=SiStripDetId::W2A;  
   modules["W3A"]=SiStripDetId::W3A;  
   modules["W1B"]=SiStripDetId::W1B;
   modules["W2B"]=SiStripDetId::W2B;
   modules["W3B"]=SiStripDetId::W3B;
   modules["W4"] =SiStripDetId::W4;
   modules["W5"] =SiStripDetId::W5;
   modules["W6"] =SiStripDetId::W6;
   modules["W7"] =SiStripDetId::W7;
   
   const unsigned size = max_element( modules.begin(),modules.end(),
                      boost::bind(&map_t::value_type::second,boost::lambda::_1) < 
                      boost::bind(&map_t::value_type::second,boost::lambda::_2) )->second + 1;
   xtalk1.resize(size);
   xtalk2.resize(size);
 
   for(map_t::const_iterator it=modules.begin(); it!=modules.end(); it++) {
     const std::string 
       modeS(peakMode_?"Peak":"Deco"),
       xtalk1S("xtalk1_" + it->first + modeS ),
       xtalk2S("xtalk2_" + it->first + modeS );
 
     if(!confObj.isParameter(xtalk1S)) throw cms::Exception("SiStripConfObject does not contain: ") << xtalk1S;
     if(!confObj.isParameter(xtalk2S)) throw cms::Exception("SiStripConfObject does not contain: ") << xtalk2S;
 
     xtalk1[it->second] = confObj.get<double>(xtalk1S);
     xtalk2[it->second] = confObj.get<double>(xtalk2S);
   }
   
   fillParams();
 
 }

Member Function Documentation

LocalVector StripCPE::driftDirection ( const StripGeomDetUnit * det ) const

virtual

Implements StripClusterParameterEstimator.

Definition at line 81 of file StripCPE.cc.

References GeomDet::geographicalId(), SiStripLorentzAngle::getLorentzAngle(), MagneticField::inTesla(), LorentzAngleMap_, magfield_, GloballyPositioned< T >::position(), DetId::rawId(), GeomDet::surface(), toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by fillParams().

                                                   { 
   LocalVector lbfield = (det->surface()).toLocal(magfield_.inTesla(det->surface().position()));  
   
   float tanLorentzAnglePerTesla = LorentzAngleMap_.getLorentzAngle(det->geographicalId().rawId());
   
   float dir_x = -tanLorentzAnglePerTesla * lbfield.y();
   float dir_y =  tanLorentzAnglePerTesla * lbfield.x();
   float dir_z =  1.f; // E field always in z direction
   
   return LocalVector(dir_x,dir_y,dir_z);
 }

void StripCPE::fillParams ( )

private

Definition at line 95 of file StripCPE.cc.

References TkRadialStripTopology::angularWidth(), assert(), BackPlaneCorrectionMap_, Surface::bounds(), TrackerGeometry::detUnits(), driftDirection(), f, GeomDet::geographicalId(), geom_, SiStripBackPlaneCorrection::getBackPlaneCorrection(), i, GeomDet::index(), GeomDetType::isTrackerStrip(), TkRadialStripTopology::localPitch(), LogDebug, m_off, m_Params, TrackerGeometry::offsetDU(), AlCaHLTBitMon_ParallelJobs::p, funct::pow(), DetId::rawId(), GeomDet::specificSurface(), StripGeomDetType::specificTopology(), StripGeomDetUnit::specificType(), mathSSE::sqrt(), TkRadialStripTopology::stripLength(), GeomDetEnumerators::tkDetEnum, StripGeomDetUnit::topology(), and StripGeomDetUnit::type().

Referenced by StripCPE().

                      {  
   auto const & dus = geom_.detUnits();
   m_off = dus.size();
   for(unsigned int i=1;i<7;++i) {
     LogDebug("LookingForFirstStrip") << "Subdetector " << i 
                      << " GeomDetEnumerator " << GeomDetEnumerators::tkDetEnum[i] 
                      << " offset " << geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) 
                      << " is it strip? " << (geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) != dus.size() ? dus[geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i])]->type().isTrackerStrip() : false);
     if(geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) != dus.size() && 
        dus[geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i])]->type().isTrackerStrip()) {
       if(geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) < m_off) m_off = geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]);
     }
   } 
   LogDebug("LookingForFirstStrip") << " Chosen offset: " << m_off;
   m_Params.resize(dus.size()-m_off);
   for (auto i=m_off; i!=dus.size();++i) {
     auto & p= m_Params[i-m_off];
     const StripGeomDetUnit * stripdet=(const StripGeomDetUnit*)(dus[i]);
     assert(stripdet->index()==int(i));
     //    assert(stripdet->geographicalId().subdetId()>1); // not pixel..
     assert(stripdet->type().isTrackerStrip()); // not pixel
 
     const Bounds& bounds = stripdet->specificSurface().bounds();
     p.maxLength = std::sqrt( std::pow(bounds.length(),2.f)+std::pow(bounds.width(),2.f) );
     p.thickness = bounds.thickness();
     p.invThickness = 1.f/p.thickness;
     p.drift = driftDirection(stripdet) * p.thickness;
     p.topology=(StripTopology*)(&stripdet->topology());    
     p.nstrips = p.topology->nstrips(); 
     p.moduleGeom = SiStripDetId(stripdet->geographicalId()).moduleGeometry();
     p.backplanecorrection = BackPlaneCorrectionMap_.getBackPlaneCorrection(stripdet->geographicalId().rawId());
     
     const TkRadialStripTopology* rtop = dynamic_cast<const TkRadialStripTopology*>(&stripdet->specificType().specificTopology());
     p.pitch_rel_err2 = (rtop) 
       ? pow( 0.5f * rtop->angularWidth() * rtop->stripLength()/rtop->localPitch(LocalPoint(0,0,0)), 2.f) / 12.f
       : 0.f;
   }
 }

AlgoParam StripCPE::getAlgoParam	(	const GeomDetUnit &	det,
		const LocalTrajectoryParameters &	ltp
	)		const

inline

Definition at line 55 of file StripCPE.h.

References funct::abs(), StripCPE::Param::backplanecorrection, corr, StripCPE::Param::coveredStrips(), StripCPE::Param::drift, f, GeomDet::geographicalId(), LocalTrajectoryParameters::momentum(), AlCaHLTBitMon_ParallelJobs::p, param(), LocalTrajectoryParameters::position(), StripCPE::Param::thickness, and PV3DBase< T, PVType, FrameType >::z().

Referenced by StripCPEfromTrackAngle::localParameters(), TkStripMeasurementDet::recHits(), and TkStripMeasurementDet::simpleRecHits().

                                                                                               {
 
     StripCPE::Param const & p = param(det);
     SiStripDetId::SubDetector loc = SiStripDetId( det.geographicalId() ).subDetector();  
  
     LocalVector track = ltp.momentum();
     track *= -p.thickness/track.z();
 
     const float fullProjection = p.coveredStrips( track+p.drift, ltp.position());
 
     auto const corr = -  0.5f*(1.f-p.backplanecorrection) * fullProjection
       + 0.5f*p.coveredStrips(track, ltp.position());
 
     return  AlgoParam{p,ltp,loc,std::abs(fullProjection),corr};
   }

StripClusterParameterEstimator::LocalValues StripCPE::localParameters	(	const SiStripCluster &	cl,
		const GeomDetUnit &	det
	)		const

overridevirtual

Reimplemented from StripClusterParameterEstimator.

Definition at line 65 of file StripCPE.cc.

References StripCPE::Param::backplanecorrection, SiStripCluster::barycenter(), StripCPE::Param::coveredStrips(), StripCPE::Param::drift, f, StripTopology::localError(), StripTopology::localPosition(), AlCaHLTBitMon_ParallelJobs::p, param(), StripCPE::Param::thickness, and StripCPE::Param::topology.

Referenced by TkStripMeasurementDet::hit().

                                                                               {
   StripCPE::Param const & p = param(det);
   const float barycenter = cluster.barycenter();
   const float fullProjection = p.coveredStrips( p.drift + LocalVector(0,0,-p.thickness), p.topology->localPosition(barycenter));
   const float strip = barycenter - 0.5f * (1.f-p.backplanecorrection) * fullProjection;
 
   return std::make_pair( p.topology->localPosition(strip),
              p.topology->localError(strip, 1.f/12.f) );
 }

virtual StripClusterParameterEstimator::LocalValues StripCPE::localParameters	(	const SiStripCluster &	cl,
		AlgoParam const &	ap
	)		const

inlinevirtual

Reimplemented in StripCPEfromTrackAngle.

Definition at line 51 of file StripCPE.h.

                                                                          {
     return std::make_pair(LocalPoint(), LocalError());
   }

Param const& StripCPE::param ( const GeomDetUnit & det ) const

inlineprotected

Definition at line 81 of file StripCPE.h.

References GeomDet::index(), m_off, and m_Params.

Referenced by getAlgoParam(), StripCPEgeometric::localParameters(), StripCPEfromTemplate::localParameters(), and localParameters().

                                                     {
     return m_Params[det.index()-m_off];
   }

Member Data Documentation

const SiStripBackPlaneCorrection& StripCPE::BackPlaneCorrectionMap_

protected

Definition at line 77 of file StripCPE.h.

Referenced by fillParams().

const TrackerGeometry& StripCPE::geom_

protected

Definition at line 74 of file StripCPE.h.

Referenced by fillParams().

const SiStripLorentzAngle& StripCPE::LorentzAngleMap_

protected

Definition at line 76 of file StripCPE.h.

Referenced by driftDirection().

unsigned int StripCPE::m_off

private

Definition at line 90 of file StripCPE.h.

Referenced by fillParams(), and param().

Params StripCPE::m_Params

private

Definition at line 89 of file StripCPE.h.

Referenced by fillParams(), and param().

const MagneticField& StripCPE::magfield_

protected

Definition at line 75 of file StripCPE.h.

Referenced by driftDirection(), and StripCPEfromTemplate::localParameters().

const bool StripCPE::peakMode_

protected

Definition at line 73 of file StripCPE.h.

Referenced by StripCPE().

std::vector<float> StripCPE::xtalk1

protected

Definition at line 78 of file StripCPE.h.

Referenced by StripCPEgeometric::localParameters(), and StripCPE().

std::vector<float> StripCPE::xtalk2

protected

Definition at line 79 of file StripCPE.h.

Referenced by StripCPE().

Classes

Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation