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	(	const TrajectoryMeasurement &	tm,
		const TrackerTopology *	tTopo,
		const TrackerGeometry *	tracker,
		const Propagator &	propagator,
		const unsigned int	mono = 0
	)

Definition at line 19 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(), bookConverter::max, min(), GluedGeomDet::monoDet(), Propagator::propagate(), DetId::rawId(), TrajectoryMeasurement::recHit(), mathSSE::sqrt(), GluedGeomDet::stereoDet(), 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);
    
    GluedGeomDet * gdet=(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 = (GeomDetUnit*)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 165 of file TrajectoryAtInvalidHit.cc.

{
  return globX;
}

double TrajectoryAtInvalidHit::globalY

(

)

const

Definition at line 169 of file TrajectoryAtInvalidHit.cc.

{
  return globY;
}

double TrajectoryAtInvalidHit::globalZ

(

)

const

Definition at line 173 of file TrajectoryAtInvalidHit.cc.

{
  return globZ;
}

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

Definition at line 193 of file TrajectoryAtInvalidHit.cc.

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

double TrajectoryAtInvalidHit::localDxDz

(

)

const

Definition at line 159 of file TrajectoryAtInvalidHit.cc.

                                               {
  return locDxDz;
}

double TrajectoryAtInvalidHit::localDyDz

(

)

const

Definition at line 162 of file TrajectoryAtInvalidHit.cc.

                                               {
  return locDyDz;
}

double TrajectoryAtInvalidHit::localErrorX

(

)

const

Definition at line 151 of file TrajectoryAtInvalidHit.cc.

{
  return locXError;
}

double TrajectoryAtInvalidHit::localErrorY

(

)

const

Definition at line 155 of file TrajectoryAtInvalidHit.cc.

{
  return locYError;
}

double TrajectoryAtInvalidHit::localX

(

)

const

Definition at line 139 of file TrajectoryAtInvalidHit.cc.

{
  return locX;
}

double TrajectoryAtInvalidHit::localY

(

)

const

Definition at line 143 of file TrajectoryAtInvalidHit.cc.

{
  return locY;
}

double TrajectoryAtInvalidHit::localZ

(

)

const

Definition at line 147 of file TrajectoryAtInvalidHit.cc.

{
  return locZ;
}

unsigned int TrajectoryAtInvalidHit::monodet_id

(

)

const

Definition at line 178 of file TrajectoryAtInvalidHit.cc.

{
  return iidd;
}

TrajectoryStateOnSurface TrajectoryAtInvalidHit::tsos

(

)

const

Definition at line 225 of file TrajectoryAtInvalidHit.cc.

                                                            {
  return theCombinedPredictedState;
}

bool TrajectoryAtInvalidHit::validHit

(

)

const

Definition at line 188 of file TrajectoryAtInvalidHit.cc.

{
  return hasValidHit;
}

bool TrajectoryAtInvalidHit::withinAcceptance

(

)

const

Definition at line 183 of file TrajectoryAtInvalidHit.cc.

{
  return acceptance;
}

Member Data Documentation

bool TrajectoryAtInvalidHit::acceptance

private

Definition at line 59 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::globX

private

Definition at line 57 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::globY

private

Definition at line 57 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::globZ

private

Definition at line 57 of file TrajectoryAtInvalidHit.h.

bool TrajectoryAtInvalidHit::hasValidHit

private

Definition at line 60 of file TrajectoryAtInvalidHit.h.

unsigned int TrajectoryAtInvalidHit::iidd

private

Definition at line 58 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::locDxDz

private

Definition at line 56 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::locDyDz

private

Definition at line 56 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::locX

private

Definition at line 54 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::locXError

private

Definition at line 55 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::locY

private

Definition at line 54 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::locYError

private

Definition at line 55 of file TrajectoryAtInvalidHit.h.

float TrajectoryAtInvalidHit::locZ

private

Definition at line 54 of file TrajectoryAtInvalidHit.h.

TrajectoryStateOnSurface TrajectoryAtInvalidHit::theCombinedPredictedState

private

Definition at line 53 of file TrajectoryAtInvalidHit.h.

TrackingRecHit::ConstRecHitPointer TrajectoryAtInvalidHit::theHit

private

Definition at line 62 of file TrajectoryAtInvalidHit.h.