#include <TrajectoryAtInvalidHit.h>

Public Member Functions
double	globalX () const

double	globalY () const

double	globalZ () const

bool	isDoubleSided (unsigned int iidd, const TrackerTopology *tTopo) const

double	localDxDz () const

double	localDyDz () const

double	localErrorX () const

double	localErrorY () const

double	localX () const

double	localY () const

double	localZ () const

unsigned int	monodet_id () const

	TrajectoryAtInvalidHit (const TrajectoryMeasurement &, const TrackerTopology tTopo, const TrackerGeometry tracker, const Propagator &propagator, const unsigned int mono=0)

TrajectoryStateOnSurface	tsos () const

bool	validHit () const

bool	withinAcceptance () const

Private Attributes
bool	acceptance

float	globX

float	globY

float	globZ

bool	hasValidHit

unsigned int	iidd

float	locDxDz

float	locDyDz

float	locX

float	locXError

float	locY

float	locYError

float	locZ

TrajectoryStateOnSurface	theCombinedPredictedState

TrackingRecHit::ConstRecHitPointer	theHit

Detailed Description

Definition at line 21 of file TrajectoryAtInvalidHit.h.

Constructor & Destructor Documentation

TrajectoryAtInvalidHit::TrajectoryAtInvalidHit	(	const TrajectoryMeasurement &	tm,
		const TrackerTopology *	tTopo,
		const TrackerGeometry *	tracker,
		const Propagator &	propagator,
		const unsigned int	mono = `0`
	)

Definition at line 18 of file TrajectoryAtInvalidHit.cc.

References funct::abs(), TrajectoryMeasurement::backwardPredictedState(), Surface::bounds(), TrajectoryStateCombiner::combine(), gather_cfg::cout, TrajectoryMeasurement::forwardPredictedState(), GeomDet::geographicalId(), TrackerGeometry::idToDet(), Bounds::inside(), TrajectoryStateOnSurface::isValid(), SiStripPI::max, min(), GluedGeomDet::monoDet(), Propagator::propagate(), DetId::rawId(), TrajectoryMeasurement::recHit(), BeamSpotPI::sigmaX, BeamSpotPI::sigmaY, mathSSE::sqrt(), GluedGeomDet::stereoDet(), digitizers_cfi::strip, DetId::subdetId(), GeomDet::surface(), and StripSubdetector::TOB.

                                                                         {
   if (tm.backwardPredictedState().isValid())
     theCombinedPredictedState =
         TrajectoryStateCombiner().combine(tm.forwardPredictedState(), tm.backwardPredictedState());
   else
     theCombinedPredictedState = tm.forwardPredictedState();
 
   if (!theCombinedPredictedState.isValid()) {
     return;
   }
   theHit = tm.recHit();
   iidd = theHit->geographicalId().rawId();
   StripSubdetector strip = StripSubdetector(iidd);
   unsigned int subid = strip.subdetId();
   // xB and yB are for absolute borders on the trajectories included in the study, sigmaX sigmaY are
   // significance cuts on the distance from the detector surface
   float xB = 0.;
   float sigmaX = 5.0;
   float yB = 0.;
   float sigmaY = 5.0;
   float sigmaYBond = 0.;
   //set bounds for point to be within to be counted in the study
   if (subid == StripSubdetector::TOB) {
     sigmaYBond = 5.0;
   }
   const GeomDetUnit* monodet;
 
   // if module is from a double sided layer, write out info for either the
   // rphi surface (mono = 1) or the stereo surface (mono = 2)--not the matched hit surface
   if ((mono > 0) && isDoubleSided(iidd, tTopo)) {
     // find matched det id, that is the matched hit surface between the two sensors
     unsigned int matched_iidd = iidd - (iidd & 0x3);
     DetId matched_id(matched_iidd);
 
     const GluedGeomDet* gdet = static_cast<const GluedGeomDet*>(tracker->idToDet(matched_id));
 
     // get the sensor det indicated by mono
     if (mono == 1)
       monodet = gdet->stereoDet();
     else
       monodet = gdet->monoDet();  // this should only be mono == 2
 
     // set theCombinedPredictedState to be on the sensor surface, not the matched surface
     DetId mono_id = monodet->geographicalId();
     const Surface& surface = tracker->idToDet(mono_id)->surface();
     theCombinedPredictedState = propagator.propagate(theCombinedPredictedState, surface);
 
     if (!theCombinedPredictedState.isValid()) {
       cout << "found invalid combinedpredictedstate after propagation" << endl;
       return;
     }
 
     //check if hit on trajectory is valid for this sensor
     if ((iidd & 0x3) == 0) {
       // passed id is on matched surface (ie neither rphi, nor stereo surface)
       // so the rechits on both sensors should be invalid
       hasValidHit = theHit->isValid();
       //cout << "should always be invalid = " << theHit->isValid() << endl;
     } else {
       if (iidd == mono_id) {
         //mono sensor for this TrajAtInvalidHit is same as originally passed in from trajectory
         //normally should only happen for valid hits, but sometimes is invalid
         hasValidHit = theHit->isValid();
         //cout << "check always true = " << theHit->isValid() << endl;
       } else {
         // mono sensor for this TrajAtInvalidHit is not as originally passed in from trajectory
         // so this is a partner missing from the trajectory and is thus invalid
         hasValidHit = false;
         //cout << "setting invalid by fiat" << endl;
       }
     }
 
     //set module id to be mono det
     iidd = monodet->geographicalId().rawId();
   } else {
     monodet = theHit->det();
     hasValidHit = theHit->isValid();
   }
 
   locX = theCombinedPredictedState.localPosition().x();
   locY = theCombinedPredictedState.localPosition().y();
   locZ = theCombinedPredictedState.localPosition().z();
   locXError = sqrt(theCombinedPredictedState.localError().positionError().xx());
   locYError = sqrt(theCombinedPredictedState.localError().positionError().yy());
   locDxDz = theCombinedPredictedState.localParameters().vector()[1];
   locDyDz = theCombinedPredictedState.localParameters().vector()[2];
   globX = theCombinedPredictedState.globalPosition().x();
   globY = theCombinedPredictedState.globalPosition().y();
   globZ = theCombinedPredictedState.globalPosition().z();
 
   // this should never be a glued det, only rphi or stero
   //cout << "From TrajAtValidHit module " << iidd << "   matched/stereo/rphi = " << ((iidd & 0x3)==0) << "/" << ((iidd & 0x3)==1) << "/" << ((iidd & 0x3)==2) << endl;
 
   // Restrict the bound regions for better understanding of the modul assignment.
 
   LocalPoint BoundedPoint;
   float xx, yy, zz;
 
   // Insert the bounded values
   if (locX < 0.)
     xx = min(locX - xB, locX - sigmaX * locXError);
   else
     xx = max(locX + xB, locX + sigmaX * locXError);
 
   if (locY < 0.)
     yy = min(locY - yB, locY - sigmaY * locYError);
   else
     yy = max(locY + yB, locY + sigmaY * locYError);
 
   zz = theCombinedPredictedState.localPosition().z();
 
   BoundedPoint = LocalPoint(xx, yy, zz);
 
   if (monodet->surface().bounds().inside(BoundedPoint) && abs(locY) > sigmaYBond * locYError) {
     acceptance = true;
   } else {
     // hit is within xB, yB from the edge of the detector, so throw it out
     acceptance = false;
   }
 }

Member Function Documentation

double TrajectoryAtInvalidHit::globalX ( ) const

Definition at line 150 of file TrajectoryAtInvalidHit.cc.

150 { return globX; }

TrajectoryAtInvalidHit::globX

float globX

Definition: TrajectoryAtInvalidHit.h:55

double TrajectoryAtInvalidHit::globalY ( ) const

Definition at line 151 of file TrajectoryAtInvalidHit.cc.

151 { return globY; }

TrajectoryAtInvalidHit::globY

float globY

Definition: TrajectoryAtInvalidHit.h:55

double TrajectoryAtInvalidHit::globalZ ( ) const

Definition at line 152 of file TrajectoryAtInvalidHit.cc.

152 { return globZ; }

TrajectoryAtInvalidHit::globZ

float globZ

Definition: TrajectoryAtInvalidHit.h:55

bool TrajectoryAtInvalidHit::isDoubleSided	(	unsigned int	iidd,
		const TrackerTopology *	tTopo
	)		const

Definition at line 160 of file TrajectoryAtInvalidHit.cc.

References phase1PixelTopology::layer, digitizers_cfi::strip, DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecRing(), StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, and TrackerTopology::tobLayer().

                                                                                                 {
   StripSubdetector strip = StripSubdetector(iidd);
   unsigned int subid = strip.subdetId();
   unsigned int layer = 0;
   if (subid == StripSubdetector::TIB) {
     layer = tTopo->tibLayer(iidd);
     if (layer == 1 || layer == 2)
       return true;
     else
       return false;
   } else if (subid == StripSubdetector::TOB) {
     layer = tTopo->tobLayer(iidd) + 4;
     if (layer == 5 || layer == 6)
       return true;
     else
       return false;
   } else if (subid == StripSubdetector::TID) {
     layer = tTopo->tidRing(iidd) + 10;
     if (layer == 11 || layer == 12)
       return true;
     else
       return false;
   } else if (subid == StripSubdetector::TEC) {
     layer = tTopo->tecRing(iidd) + 13;
     if (layer == 14 || layer == 15 || layer == 18)
       return true;
     else
       return false;
   } else
     return false;
 }