|
|
|
#include <GSRecHitMatcher.h>
Public Types | |
| typedef std::pair< LocalPoint, LocalPoint > | StripPosition |
Public Member Functions | |
| GSRecHitMatcher () | |
| SiTrackerGSMatchedRecHit2D * | match (const SiTrackerGSRecHit2D *monoRH, const SiTrackerGSRecHit2D *stereoRH, const GluedGeomDet *gluedDet, LocalVector &trackdirection) const |
| StripPosition | project (const GeomDetUnit *det, const GluedGeomDet *glueddet, const StripPosition &strip, const LocalVector &trackdirection) const |
| SiTrackerGSMatchedRecHit2D * | projectOnly (const SiTrackerGSRecHit2D *monoRH, const GeomDet *monoDet, const GluedGeomDet *gluedDet, LocalVector &ldir) const |
| ~GSRecHitMatcher () | |
Definition at line 12 of file GSRecHitMatcher.h.
| typedef std::pair<LocalPoint,LocalPoint> GSRecHitMatcher::StripPosition |
Definition at line 15 of file GSRecHitMatcher.h.
| GSRecHitMatcher::GSRecHitMatcher | ( | ) | [inline] |
Definition at line 17 of file GSRecHitMatcher.h.
{}
| GSRecHitMatcher::~GSRecHitMatcher | ( | ) | [inline] |
Definition at line 18 of file GSRecHitMatcher.h.
{}
| SiTrackerGSMatchedRecHit2D * GSRecHitMatcher::match | ( | const SiTrackerGSRecHit2D * | monoRH, |
| const SiTrackerGSRecHit2D * | stereoRH, | ||
| const GluedGeomDet * | gluedDet, | ||
| LocalVector & | trackdirection | ||
| ) | const |
Definition at line 9 of file GSRecHitMatcher.cc.
References trackerHits::c, alignmentValidation::c1, SiTrackerGSRecHit2D::cluster(), funct::cos(), diffTreeTool::diff, SiTrackerGSRecHit2D::eeId(), error, Exception, TrackingRecHit::geographicalId(), GeomDet::geographicalId(), StripTopology::localPitch(), GSSiTrackerRecHit2DLocalPos::localPosition(), StripTopology::localPosition(), GSSiTrackerRecHit2DLocalPos::localPositionError(), m, PV3DBase< T, PVType, FrameType >::mag2(), Topology::measurementError(), Topology::measurementPosition(), GluedGeomDet::monoDet(), PV3DBase< T, PVType, FrameType >::phi(), position, project(), indexGen::s2, SiTrackerGSRecHit2D::simhitId(), SiTrackerGSRecHit2D::simMultX(), SiTrackerGSRecHit2D::simMultY(), SiTrackerGSRecHit2D::simtrackId(), funct::sin(), GluedGeomDet::stereoDet(), GeomDet::surface(), toLocal(), GeomDetUnit::topology(), PV3DBase< T, PVType, FrameType >::x(), PV2DBase< T, PVType, FrameType >::x(), and xy().
Referenced by SiTrackerGaussianSmearingRecHitConverter::matchHits().
{
// stripdet = mono
// partnerstripdet = stereo
const GeomDetUnit* stripdet = gluedDet->monoDet();
const GeomDetUnit* partnerstripdet = gluedDet->stereoDet();
const StripTopology& topol=(const StripTopology&)stripdet->topology();
LocalPoint position;
// position of the initial and final point of the strip (RPHI cluster) in local strip coordinates
MeasurementPoint RPHIpoint=topol.measurementPosition(monoRH->localPosition());
MeasurementPoint RPHIpointini=MeasurementPoint(RPHIpoint.x(),-0.5);
MeasurementPoint RPHIpointend=MeasurementPoint(RPHIpoint.x(),0.5);
// position of the initial and final point of the strip in local coordinates (mono det)
StripPosition stripmono=StripPosition(topol.localPosition(RPHIpointini),topol.localPosition(RPHIpointend));
if(trackdirection.mag2()<FLT_MIN){// in case of no track hypothesis assume a track from the origin through the center of the strip
LocalPoint lcenterofstrip=monoRH->localPosition();
GlobalPoint gcenterofstrip=(stripdet->surface()).toGlobal(lcenterofstrip);
GlobalVector gtrackdirection=gcenterofstrip-GlobalPoint(0,0,0);
trackdirection=(gluedDet->surface()).toLocal(gtrackdirection);
}
//project mono hit on glued det
StripPosition projectedstripmono=project(stripdet,gluedDet,stripmono,trackdirection);
const StripTopology& partnertopol=(const StripTopology&)partnerstripdet->topology();
//error calculation (the part that depends on mono RH only)
LocalVector RPHIpositiononGluedendvector=projectedstripmono.second-projectedstripmono.first;
double c1=sin(RPHIpositiononGluedendvector.phi());
double s1=-cos(RPHIpositiononGluedendvector.phi());
MeasurementError errormonoRH=topol.measurementError(monoRH->localPosition(),monoRH->localPositionError());
double pitch=topol.localPitch(monoRH->localPosition());
double sigmap12=errormonoRH.uu()*pitch*pitch;
// position of the initial and final point of the strip (STEREO cluster)
MeasurementPoint STEREOpoint=partnertopol.measurementPosition(stereoRH->localPosition());
MeasurementPoint STEREOpointini=MeasurementPoint(STEREOpoint.x(),-0.5);
MeasurementPoint STEREOpointend=MeasurementPoint(STEREOpoint.x(),0.5);
// position of the initial and final point of the strip in local coordinates (stereo det)
StripPosition stripstereo(partnertopol.localPosition(STEREOpointini),partnertopol.localPosition(STEREOpointend));
//project stereo hit on glued det
StripPosition projectedstripstereo=project(partnerstripdet,gluedDet,stripstereo,trackdirection);
//perform the matching
//(x2-x1)(y-y1)=(y2-y1)(x-x1)
AlgebraicMatrix22 m; AlgebraicVector2 c, solution;
m(0,0)=-(projectedstripmono.second.y()-projectedstripmono.first.y()); m(0,1)=(projectedstripmono.second.x()-projectedstripmono.first.x());
m(1,0)=-(projectedstripstereo.second.y()-projectedstripstereo.first.y()); m(1,1)=(projectedstripstereo.second.x()-projectedstripstereo.first.x());
c(0)=m(0,1)*projectedstripmono.first.y()+m(0,0)*projectedstripmono.first.x();
c(1)=m(1,1)*projectedstripstereo.first.y()+m(1,0)*projectedstripstereo.first.x();
m.Invert(); solution = m * c;
position=LocalPoint(solution(0),solution(1));
//
// temporary fix by tommaso
//
LocalError tempError (100,0,100);
// calculate the error
LocalVector stereopositiononGluedendvector=projectedstripstereo.second-projectedstripstereo.first;
double c2=sin(stereopositiononGluedendvector.phi()); double s2=-cos(stereopositiononGluedendvector.phi());
MeasurementError errorstereoRH=partnertopol.measurementError(stereoRH->localPosition(), stereoRH->localPositionError());
pitch=partnertopol.localPitch(stereoRH->localPosition());
double sigmap22=errorstereoRH.uu()*pitch*pitch;
double diff=(c1*s2-c2*s1);
double invdet2=1/(diff*diff);
float xx=invdet2*(sigmap12*s2*s2+sigmap22*s1*s1);
float xy=-invdet2*(sigmap12*c2*s2+sigmap22*c1*s1);
float yy=invdet2*(sigmap12*c2*c2+sigmap22*c1*c1);
LocalError error=LocalError(xx,xy,yy);
// if((gluedDet->surface()).bounds().inside(position,error,3)){
// std::cout<<" ERROR ok "<< std::endl;
// }
// else {
// std::cout<<" ERROR not ok " << std::endl;
// }
//Added by DAO to make sure y positions are zero.
DetId det(monoRH->geographicalId());
if(det.subdetId() > 2) {
SiTrackerGSRecHit2D *adjustedMonoRH = new SiTrackerGSRecHit2D(LocalPoint(monoRH->localPosition().x(),0,0),
monoRH->localPositionError(),
monoRH->geographicalId(),
monoRH->simhitId(),
monoRH->simtrackId(),
monoRH->eeId(),
monoRH->cluster(),
monoRH->simMultX(),
monoRH->simMultY()
);
SiTrackerGSRecHit2D *adjustedStereoRH = new SiTrackerGSRecHit2D(LocalPoint(stereoRH->localPosition().x(),0,0),
stereoRH->localPositionError(),
stereoRH->geographicalId(),
stereoRH->simhitId(),
stereoRH->simtrackId(),
stereoRH->eeId(),
stereoRH->cluster(),
stereoRH->simMultX(),
stereoRH->simMultY()
);
SiTrackerGSMatchedRecHit2D *rV= new SiTrackerGSMatchedRecHit2D(position, error,gluedDet->geographicalId(), monoRH->simhitId(),
monoRH->simtrackId(), monoRH->eeId(), monoRH->cluster(),
monoRH->simMultX(), monoRH->simMultY(),
true, adjustedMonoRH, adjustedStereoRH);
delete adjustedMonoRH;
delete adjustedStereoRH;
return rV;
}
else {
throw cms::Exception("GSRecHitMatcher") << "Matched Pixel!?";
}
}
| GSRecHitMatcher::StripPosition GSRecHitMatcher::project | ( | const GeomDetUnit * | det, |
| const GluedGeomDet * | glueddet, | ||
| const StripPosition & | strip, | ||
| const LocalVector & | trackdirection | ||
| ) | const |
Definition at line 142 of file GSRecHitMatcher.cc.
References GeomDet::surface(), toLocal(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by match().
{
GlobalPoint globalpointini=(det->surface()).toGlobal(strip.first);
GlobalPoint globalpointend=(det->surface()).toGlobal(strip.second);
// position of the initial and final point of the strip in glued local coordinates
LocalPoint positiononGluedini=(glueddet->surface()).toLocal(globalpointini);
LocalPoint positiononGluedend=(glueddet->surface()).toLocal(globalpointend);
//correct the position with the track direction
float scale=-positiononGluedini.z()/trackdirection.z();
LocalPoint projpositiononGluedini= positiononGluedini + scale*trackdirection;
LocalPoint projpositiononGluedend= positiononGluedend + scale*trackdirection;
return StripPosition(projpositiononGluedini,projpositiononGluedend);
}
| SiTrackerGSMatchedRecHit2D * GSRecHitMatcher::projectOnly | ( | const SiTrackerGSRecHit2D * | monoRH, |
| const GeomDet * | monoDet, | ||
| const GluedGeomDet * | gluedDet, | ||
| LocalVector & | ldir | ||
| ) | const |
Definition at line 167 of file GSRecHitMatcher.cc.
References SiTrackerGSRecHit2D::cluster(), delta, Vector3DBase< T, FrameTag >::dot(), SiTrackerGSRecHit2D::eeId(), Exception, TrackingRecHit::geographicalId(), GeomDet::geographicalId(), GSSiTrackerRecHit2DLocalPos::localPosition(), GSSiTrackerRecHit2DLocalPos::localPositionError(), Plane::localZ(), GluedGeomDet::monoDet(), Plane::normalVector(), position, GloballyPositioned< T >::position(), LocalError::rotate(), SiTrackerGSRecHit2D::simhitId(), SiTrackerGSRecHit2D::simMultX(), SiTrackerGSRecHit2D::simMultY(), SiTrackerGSRecHit2D::simtrackId(), GluedGeomDet::stereoDet(), GeomDet::surface(), Surface::toGlobal(), GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by SiTrackerGaussianSmearingRecHitConverter::matchHits().
{
LocalPoint position(monoRH->localPosition().x(), 0.,0.);
const BoundPlane& gluedPlane = gluedDet->surface();
const BoundPlane& hitPlane = monoDet->surface();
double delta = gluedPlane.localZ( hitPlane.position());
LocalPoint lhitPos = gluedPlane.toLocal( monoDet->surface().toGlobal(position ) );
LocalPoint projectedHitPos = lhitPos - ldir * delta/ldir.z();
LocalVector hitXAxis = gluedPlane.toLocal( hitPlane.toGlobal( LocalVector(1,0,0)));
LocalError hitErr = monoRH->localPositionError();
if (gluedPlane.normalVector().dot( hitPlane.normalVector()) < 0) {
// the two planes are inverted, and the correlation element must change sign
hitErr = LocalError( hitErr.xx(), -hitErr.xy(), hitErr.yy());
}
LocalError rotatedError = hitErr.rotate( hitXAxis.x(), hitXAxis.y());
const GeomDetUnit *gluedMonoDet = gluedDet->monoDet();
const GeomDetUnit *gluedStereoDet = gluedDet->stereoDet();
int isMono = 0;
int isStereo = 0;
if(monoDet->geographicalId()==gluedMonoDet->geographicalId()) isMono = 1;
if(monoDet->geographicalId()==gluedStereoDet->geographicalId()) isStereo = 1;
//Added by DAO to make sure y positions are zero and correct Mono or stereo Det is filled.
//Good for debugging.
//std::cout << "The monoDet = " << monoDet->geographicalId() << ". The gluedMonoDet = " << gluedMonoDet->geographicalId() << ". The gluedStereoDet = " << gluedStereoDet->geographicalId()
// << ". isMono = " << isMono << ". isStereo = " << isStereo <<"." << std::endl;
SiTrackerGSRecHit2D *adjustedRH = new SiTrackerGSRecHit2D(LocalPoint(monoRH->localPosition().x(),0,0),
monoRH->localPositionError(),
monoRH->geographicalId(),
monoRH->simhitId(),
monoRH->simtrackId(),
monoRH->eeId(),
monoRH->cluster(),
monoRH->simMultX(),
monoRH->simMultY()
);
//DAO: Not quite sure what to do about the cluster ref, so I will fill it with the monoRH for now...
SiTrackerGSRecHit2D *otherRH = new SiTrackerGSRecHit2D(LocalPoint(-10000,-10000,-10000), LocalError(0,0,0),DetId(000000000), 0,0,0,monoRH->cluster(),0,0);
if ((isMono && isStereo)||(!isMono&&!isStereo)) throw cms::Exception("GSRecHitMatcher") << "Something wrong with DetIds.";
else if (isMono) {
SiTrackerGSMatchedRecHit2D *rV= new SiTrackerGSMatchedRecHit2D(projectedHitPos, rotatedError, gluedDet->geographicalId(),
monoRH->simhitId(), monoRH->simtrackId(), monoRH->eeId(), monoRH->cluster(),
monoRH->simMultX(), monoRH->simMultY(), false, adjustedRH, otherRH);
delete adjustedRH;
delete otherRH;
return rV;
}
else{
SiTrackerGSMatchedRecHit2D *rV=new SiTrackerGSMatchedRecHit2D(projectedHitPos, rotatedError, gluedDet->geographicalId(),
monoRH->simhitId(), monoRH->simtrackId(), monoRH->eeId(), monoRH->cluster(),
monoRH->simMultX(), monoRH->simMultY(), false, otherRH, adjustedRH);
delete adjustedRH;
delete otherRH;
return rV;
}
}
1.7.3