TrajectoryAtInvalidHit Class Reference

#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::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.

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

References funct::abs(), TrajectoryMeasurement::backwardPredictedState(), Surface::bounds(), TrajectoryStateCombiner::combine(), gather_cfg::cout, TrajectoryMeasurement::forwardPredictedState(), GeomDet::geographicalId(), Bounds::inside(), TrajectoryStateOnSurface::isValid(), SiStripPI::max, min(), GluedGeomDet::monoDet(), TrackCandidateProducer_cfi::propagator, DetId::rawId(), TrajectoryMeasurement::recHit(), fftjetvertexadder_cfi::sigmaX, fftjetvertexadder_cfi::sigmaY, mathSSE::sqrt(), GluedGeomDet::stereoDet(), digitizers_cfi::strip, GeomDet::surface(), StripSubdetector::TOB, PbPb_ZMuSkimMuonDPG_cff::tracker, geometryCSVtoXML::xx, geometryCSVtoXML::yy, and geometryCSVtoXML::zz.

Member Function Documentation

globalX()

double TrajectoryAtInvalidHit::globalX

(

)

const

Definition at line 150 of file TrajectoryAtInvalidHit.cc.

{ return globX; }

globalY()

double TrajectoryAtInvalidHit::globalY

(

)

const

Definition at line 151 of file TrajectoryAtInvalidHit.cc.

{ return globY; }

globalZ()

double TrajectoryAtInvalidHit::globalZ

(

)

const

Definition at line 152 of file TrajectoryAtInvalidHit.cc.

{ return globZ; }

isDoubleSided()

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

Definition at line 160 of file TrajectoryAtInvalidHit.cc.

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

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

localDxDz()

double TrajectoryAtInvalidHit::localDxDz

(

)

const

Definition at line 148 of file TrajectoryAtInvalidHit.cc.

{ return locDxDz; }

localDyDz()

double TrajectoryAtInvalidHit::localDyDz

(

)

const

Definition at line 149 of file TrajectoryAtInvalidHit.cc.

{ return locDyDz; }

localErrorX()

double TrajectoryAtInvalidHit::localErrorX

(

)

const

Definition at line 146 of file TrajectoryAtInvalidHit.cc.

{ return locXError; }

localErrorY()

double TrajectoryAtInvalidHit::localErrorY

(

)

const

Definition at line 147 of file TrajectoryAtInvalidHit.cc.

{ return locYError; }

localX()

double TrajectoryAtInvalidHit::localX

(

)

const

Definition at line 143 of file TrajectoryAtInvalidHit.cc.

{ return locX; }

localY()

double TrajectoryAtInvalidHit::localY

(

)

const

Definition at line 144 of file TrajectoryAtInvalidHit.cc.

{ return locY; }

localZ()

double TrajectoryAtInvalidHit::localZ

(

)

const

Definition at line 145 of file TrajectoryAtInvalidHit.cc.

{ return locZ; }

monodet_id()

unsigned int TrajectoryAtInvalidHit::monodet_id

(

)

const

Definition at line 154 of file TrajectoryAtInvalidHit.cc.

{ return iidd; }

tsos()

TrajectoryStateOnSurface TrajectoryAtInvalidHit::tsos

(

)

const

Definition at line 192 of file TrajectoryAtInvalidHit.cc.

{ return theCombinedPredictedState; }

validHit()

bool TrajectoryAtInvalidHit::validHit

(

)

const

Definition at line 158 of file TrajectoryAtInvalidHit.cc.

{ return hasValidHit; }

withinAcceptance()

bool TrajectoryAtInvalidHit::withinAcceptance

(

)

const

Definition at line 156 of file TrajectoryAtInvalidHit.cc.

{ return acceptance; }

Member Data Documentation

acceptance

bool TrajectoryAtInvalidHit::acceptance

private

Definition at line 57 of file TrajectoryAtInvalidHit.h.

globX

float TrajectoryAtInvalidHit::globX

private

Definition at line 55 of file TrajectoryAtInvalidHit.h.

globY

float TrajectoryAtInvalidHit::globY

private

Definition at line 55 of file TrajectoryAtInvalidHit.h.

globZ

float TrajectoryAtInvalidHit::globZ

private

Definition at line 55 of file TrajectoryAtInvalidHit.h.

hasValidHit

bool TrajectoryAtInvalidHit::hasValidHit

private

Definition at line 58 of file TrajectoryAtInvalidHit.h.

iidd

unsigned int TrajectoryAtInvalidHit::iidd

private

Definition at line 56 of file TrajectoryAtInvalidHit.h.

locDxDz

float TrajectoryAtInvalidHit::locDxDz

private

Definition at line 54 of file TrajectoryAtInvalidHit.h.

locDyDz

float TrajectoryAtInvalidHit::locDyDz

private

Definition at line 54 of file TrajectoryAtInvalidHit.h.

locX

float TrajectoryAtInvalidHit::locX

private

Definition at line 52 of file TrajectoryAtInvalidHit.h.

locXError

float TrajectoryAtInvalidHit::locXError

private

Definition at line 53 of file TrajectoryAtInvalidHit.h.

locY

float TrajectoryAtInvalidHit::locY

private

Definition at line 52 of file TrajectoryAtInvalidHit.h.

locYError

float TrajectoryAtInvalidHit::locYError

private

Definition at line 53 of file TrajectoryAtInvalidHit.h.

locZ

float TrajectoryAtInvalidHit::locZ

private

Definition at line 52 of file TrajectoryAtInvalidHit.h.

theCombinedPredictedState

TrajectoryStateOnSurface TrajectoryAtInvalidHit::theCombinedPredictedState

private

Definition at line 51 of file TrajectoryAtInvalidHit.h.

theHit

TrackingRecHit::ConstRecHitPointer TrajectoryAtInvalidHit::theHit

private

Definition at line 60 of file TrajectoryAtInvalidHit.h.