|
|
|
#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 |
| ConstReferenceCountingPointer < TransientTrackingRecHit > | theHit |
Definition at line 21 of file TrajectoryAtInvalidHit.h.
| 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 abs, TrajectoryMeasurement::backwardPredictedState(), Surface::bounds(), TrajectoryStateCombiner::combine(), gather_cfg::cout, TrajectoryMeasurement::forwardPredictedState(), GeomDet::geographicalId(), TrackerGeometry::idToDet(), Bounds::inside(), TrajectoryStateOnSurface::isValid(), 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;
}
}
| 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;
}
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.
Definition at line 53 of file TrajectoryAtInvalidHit.h.
Definition at line 62 of file TrajectoryAtInvalidHit.h.
1.7.3