#include <MultipleScatteringGeometry.h>

Public Member Functions
std::vector< MSLayer >	detLayers (float eta, float z, const MagneticField &bfield) const

std::vector< MSLayer >	detLayers () const

	MultipleScatteringGeometry (const GeometricSearchTracker &tracker)

std::vector< MSLayer >	otherLayers (float eta) const

Protected Attributes
std::vector< const DetLayer * >	theLayers

Static Protected Attributes
static const float	beamPipeR = 2.94

static const float	endflangesZ = 30

static const float	supportR = 19.

Detailed Description

Definition at line 12 of file MultipleScatteringGeometry.h.

Constructor & Destructor Documentation

◆ MultipleScatteringGeometry()

MultipleScatteringGeometry::MultipleScatteringGeometry ( const GeometricSearchTracker & tracker )

Definition at line 17 of file MultipleScatteringGeometry.cc.

References cuy::ib, and PbPb_ZMuSkimMuonDPG_cff::tracker.

                                                                                             {
   vector<BarrelDetLayer const *> barrelLayers = tracker.barrelLayers();
   vector<BarrelDetLayer const *>::const_iterator ib;
   vector<ForwardDetLayer const *> forwardPosLayers = tracker.posForwardLayers();
   vector<ForwardDetLayer const *> forwardNegLayers = tracker.negForwardLayers();
   vector<ForwardDetLayer const *>::const_iterator ie;
   // barrelLayers = accessor.barrelLayers();
   for (ib = barrelLayers.begin(); ib != barrelLayers.end(); ib++)
     theLayers.push_back(*ib);
 
   //   forwardLayers = accessor.negativeEndcapLayers();
   for (ie = forwardPosLayers.begin(); ie != forwardPosLayers.end(); ie++)
     theLayers.push_back(*ie);
   for (ie = forwardNegLayers.begin(); ie != forwardNegLayers.end(); ie++)
     theLayers.push_back(*ie);
 }

Member Function Documentation

◆ detLayers() [1/2]

vector< MSLayer > MultipleScatteringGeometry::detLayers	(	float	eta,
		float	z,
		const MagneticField &	bfield
	)		const

Definition at line 44 of file MultipleScatteringGeometry.cc.

References alongMomentum, GeomDetEnumerators::barrel, edm::contains(), DeadROC_duringRun::dir, flavorHistoryFilter_cfi::dr, PVValHelper::dz, GeomDetEnumerators::endcap, PVValHelper::eta, TrajectoryStateOnSurface::globalPosition(), PixelRecoRange< T >::inside(), TrajectoryStateOnSurface::isValid(), PV3DBase< T, PVType, FrameType >::perp(), GeometricSearchDet::position(), TrackCandidateProducer_cfi::propagator, alignCSCRings::r, mps_fire::result, BarrelDetLayer::specificSurface(), ForwardDetLayer::specificSurface(), PV3DBase< T, PVType, FrameType >::z(), SiPixelPI::zero, photonAnalyzer_cfi::zMax, photonAnalyzer_cfi::zMin, and beampixel_dqm_sourceclient-live_cfg::zRange.

                                                                                                            {
   vector<MSLayer> result;
   GlobalPoint zero(0, 0, z);
   float r = 1;
   float dirZ = r * sinh(eta);
   GlobalVector dir(r, 0., dirZ);
   FreeTrajectoryState fts(GlobalTrajectoryParameters(zero, dir, 1, &bfield));
   StraightLinePropagator propagator(&bfield, alongMomentum);
   vector<const DetLayer *>::const_iterator il;
   TrajectoryStateOnSurface tsos;
   for (il = theLayers.begin(); il != theLayers.end(); il++) {
     bool contains = false;
     //  if ((*il)->subDetector() != PixelBarrel && (*il)->subDetector()!= PixelEndcap) continue;
 
     if ((*il)->location() == barrel) {
       const BarrelDetLayer *bl = dynamic_cast<const BarrelDetLayer *>(*il);
       if (!bl)
         continue;
       tsos = propagator.propagate(fts, bl->specificSurface());
       if (!tsos.isValid())
         continue;
       float r = bl->specificSurface().radius();
       float dr = bl->specificSurface().bounds().thickness();
       float z = bl->position().z();
       float dz = bl->specificSurface().bounds().length();
       PixelRecoRange<float> rRange(r - dr / 2., r + dr / 2.);
       PixelRecoRange<float> zRange(z - dz / 2., z + dz / 2.);
       contains = rRange.inside(tsos.globalPosition().perp()) && zRange.inside(tsos.globalPosition().z());
       //     contains = bl->contains((*il)->toLocal(tsos.globalPosition()));
     } else if ((*il)->location() == endcap) {
       const ForwardDetLayer *fl = dynamic_cast<const ForwardDetLayer *>(*il);
       if (!fl)
         continue;
       if (fl->position().z() * eta < 0)
         continue;
       const BoundDisk &disk = fl->specificSurface();
       tsos = propagator.propagate(fts, disk);
       if (!tsos.isValid())
         continue;
       float zMin = disk.position().z() - disk.bounds().thickness() / 2;
       float zMax = disk.position().z() + disk.bounds().thickness() / 2;
       PixelRecoRange<float> rRange(disk.innerRadius(), disk.outerRadius());
       PixelRecoRange<float> zRange(zMin, zMax);
       contains = rRange.inside(tsos.globalPosition().perp()) && zRange.inside(tsos.globalPosition().z());
       //     contains = fl->contains(fl->toLocal(tsos.globalPosition()));
     }
     if (contains)
       result.push_back(MSLayer(*il));
   }
   return result;
 }

◆ detLayers() [2/2]

vector< MSLayer > MultipleScatteringGeometry::detLayers ( ) const

Definition at line 35 of file MultipleScatteringGeometry.cc.

References mps_fire::result.

                                                             {
   vector<MSLayer> result;
   vector<const DetLayer *>::const_iterator il;
   for (il = theLayers.begin(); il != theLayers.end(); il++)
     result.push_back(MSLayer(*il));
   return result;
 }

◆ otherLayers()

vector< MSLayer > MultipleScatteringGeometry::otherLayers ( float eta ) const

Definition at line 96 of file MultipleScatteringGeometry.cc.

References GeomDetEnumerators::barrel, GeomDetEnumerators::endcap, PVValHelper::eta, PixelRecoPointRZ::r(), mps_fire::result, and PixelRecoPointRZ::z().

                                                                        {
   vector<MSLayer> result;
   // zero
   //  MSLayer zero(barrel, 0., MSLayer::Range(-15,15));
   //  result.push_back(zero);
 
   // beampipe
   MSLayer beampipe(barrel, beamPipeR, MSLayer::Range(-100, 100));
 
   result.push_back(beampipe);
 
   // endflanges
   PixelRecoPointRZ endfPoint = (eta > 0) ? PixelRecoPointRZ(endflangesZ / sinh(eta), endflangesZ)
                                          : PixelRecoPointRZ(-endflangesZ / sinh(eta), -endflangesZ);
   if (0 < endfPoint.r() && endfPoint.r() < supportR) {
     MSLayer endflanges(endcap, endfPoint.z(), MSLayer::Range(0.1, supportR - 0.1));
     result.push_back(endflanges);
   }
 
   // support
   MSLayer support(barrel, supportR, MSLayer::Range(-280, 280));
   result.push_back(support);
 
   return result;
 }