StripCPE Class Reference

#include <StripCPE.h>

Inheritance diagram for StripCPE:

Classes
struct	Param

Public Member Functions
LocalVector	driftDirection (const StripGeomDetUnit *det) const
StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, const GeomDetUnit &) const
	StripCPE (edm::ParameterSet &conf, const MagneticField &, const TrackerGeometry &, const SiStripLorentzAngle &, const SiStripBackPlaneCorrection &, const SiStripConfObject &, const SiStripLatency &)
Public Member Functions inherited from StripClusterParameterEstimator
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
typedef std::pair< LocalPoint, LocalError >	LocalValues
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 58 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.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;
  }
}

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

virtual

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.

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

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

inlineprotected

Definition at line 51 of file StripCPE.h.

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

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

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

Member Data Documentation

const SiStripBackPlaneCorrection& StripCPE::BackPlaneCorrectionMap_

protected

Definition at line 37 of file StripCPE.h.

Referenced by fillParams().

const TrackerGeometry& StripCPE::geom_

protected

Definition at line 34 of file StripCPE.h.

Referenced by fillParams().

const SiStripLorentzAngle& StripCPE::LorentzAngleMap_

protected

Definition at line 36 of file StripCPE.h.

Referenced by driftDirection().

unsigned int StripCPE::m_off

private

Definition at line 60 of file StripCPE.h.

Referenced by fillParams(), and param().

Params StripCPE::m_Params

private

Definition at line 59 of file StripCPE.h.

Referenced by fillParams(), and param().

const MagneticField& StripCPE::magfield_

protected

Definition at line 35 of file StripCPE.h.

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

const bool StripCPE::peakMode_

protected

Definition at line 33 of file StripCPE.h.

Referenced by StripCPE().

std::vector<float> StripCPE::xtalk1

protected

Definition at line 38 of file StripCPE.h.

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

std::vector<float> StripCPE::xtalk2

protected

Definition at line 39 of file StripCPE.h.

Referenced by StripCPE().