|
|
|
#include <TrackInfoProducerAlgorithm.h>
Public Types | |
| typedef TransientTrackingRecHit::ConstRecHitPointer | ConstRecHitPointer |
Public Member Functions | |
| LocalPoint | project (const GeomDet *det, const GeomDet *projdet, LocalPoint position, LocalVector trackdirection) const |
| void | run (const edm::Ref< std::vector< Trajectory > > traj_iterator, reco::TrackRef track, reco::TrackInfo &output, const TrackerGeometry *tracker) |
| TrackInfoProducerAlgorithm (const edm::ParameterSet &conf) | |
| ~TrackInfoProducerAlgorithm () | |
Private Attributes | |
| std::string | backwardPredictedStateTag_ |
| std::string | combinedStateTag_ |
| edm::ParameterSet | conf_ |
| std::string | forwardPredictedStateTag_ |
| std::string | updatedStateTag_ |
Definition at line 20 of file TrackInfoProducerAlgorithm.h.
Definition at line 23 of file TrackInfoProducerAlgorithm.h.
| TrackInfoProducerAlgorithm::TrackInfoProducerAlgorithm | ( | const edm::ParameterSet & | conf | ) | [inline] |
Definition at line 25 of file TrackInfoProducerAlgorithm.h.
:
conf_(conf),
forwardPredictedStateTag_(conf.getParameter<std::string>( "forwardPredictedState" )),
backwardPredictedStateTag_(conf.getParameter<std::string>( "backwardPredictedState" )),
updatedStateTag_(conf.getParameter<std::string>( "updatedState" )),
combinedStateTag_(conf.getParameter<std::string>( "combinedState" )) { }
| TrackInfoProducerAlgorithm::~TrackInfoProducerAlgorithm | ( | ) | [inline] |
Definition at line 32 of file TrackInfoProducerAlgorithm.h.
{}
| LocalPoint TrackInfoProducerAlgorithm::project | ( | const GeomDet * | det, |
| const GeomDet * | projdet, | ||
| LocalPoint | position, | ||
| LocalVector | trackdirection | ||
| ) | const |
Definition at line 168 of file TrackInfoProducerAlgorithm.cc.
References pileupReCalc_HLTpaths::scale, GeomDet::surface(), toLocal(), and PV3DBase< T, PVType, FrameType >::z().
{
GlobalPoint globalpoint=(det->surface()).toGlobal(position);
// position of the initial and final point of the strip in glued local coordinates
LocalPoint projposition=(projdet->surface()).toLocal(globalpoint);
//correct the position with the track direction
float scale=-projposition.z()/trackdirection.z();
projposition+= scale*trackdirection;
return projposition;
}
| void TrackInfoProducerAlgorithm::run | ( | const edm::Ref< std::vector< Trajectory > > | traj_iterator, |
| reco::TrackRef | track, | ||
| reco::TrackInfo & | output, | ||
| const TrackerGeometry * | tracker | ||
| ) |
Definition at line 16 of file TrackInfoProducerAlgorithm.cc.
References reco::BwPredicted, TrajectoryStateCombiner::combine(), reco::Combined, cond::rpcobgas::detid, alignCSCRings::e, reco::FwPredicted, TrackingRecHit::geographicalId(), i, TrackerGeometry::idToDet(), j, mag(), reco::Matched, LocalTrajectoryParameters::momentum(), GluedGeomDet::monoDet(), ProjectedSiStripRecHit2D::originalHit(), PTrajectoryStateOnDet::parameters(), trajectoryStateTransform::persistentState(), pos, LocalTrajectoryParameters::position(), CrabWatch::project, reco::Projected, DetId::rawId(), reco::Single, GluedGeomDet::stereoDet(), GeomDet::toGlobal(), GeomDet::toLocal(), and reco::Updated.
Referenced by TrackInfoProducer::produce().
{
std::vector<TrajectoryMeasurement> measurements =traj_iterator->measurements();
std::vector<TrajectoryMeasurement>::iterator traj_mes_iterator;
//edm::LogInfo("TrackInfoProducer") << "Number of Measurements: "<<measurements.size();
TrackInfo::TrajectoryInfo trajinfo;
int nhit=0;
for(traj_mes_iterator=measurements.begin();traj_mes_iterator!=measurements.end();traj_mes_iterator++){//loop on measurements
TrajectoryStateOnSurface fwdtsos=traj_mes_iterator->forwardPredictedState();
TrajectoryStateOnSurface bwdtsos=traj_mes_iterator->backwardPredictedState();
TrajectoryStateOnSurface updatedtsos=traj_mes_iterator->updatedState();
TrajectoryStateCombiner statecombiner;
TrajectoryStateOnSurface combinedtsos=statecombiner.combine(fwdtsos, bwdtsos);
ConstRecHitPointer ttrh=traj_mes_iterator->recHit();
LocalPoint pos;
if (ttrh->isValid())pos=ttrh->hit()->localPosition() ;
nhit++;
unsigned int detid=ttrh->hit()->geographicalId().rawId();
trackingRecHit_iterator thehit;
TrackingRecHitRef thehitref;
int i=0,j=0;
for (thehit=track->recHitsBegin();thehit!=track->recHitsEnd();thehit++){
i++;
LocalPoint hitpos;
if ((*thehit)->isValid())hitpos=(*thehit)->localPosition();
if((*thehit)->geographicalId().rawId()==detid&&
(hitpos - pos).mag() < 1e-4)
{
thehitref=(*thehit);
j++;
break;
}
}
PTrajectoryStateOnDet const & fwdptsod=trajectoryStateTransform::persistentState( fwdtsos,detid);
PTrajectoryStateOnDet const & bwdptsod=trajectoryStateTransform::persistentState( bwdtsos,detid);
PTrajectoryStateOnDet const & updatedptsod=trajectoryStateTransform::persistentState( updatedtsos,detid);
PTrajectoryStateOnDet const & combinedptsod=trajectoryStateTransform::persistentState( combinedtsos,detid);
const ProjectedSiStripRecHit2D* phit=dynamic_cast<const ProjectedSiStripRecHit2D*>( &*(thehitref));
const SiStripMatchedRecHit2D* matchedhit=dynamic_cast<const SiStripMatchedRecHit2D*>( &*(thehitref));
RecHitType type=Single;
LocalVector monofwd, stereofwd;
LocalVector monobwd, stereobwd;
LocalVector monoco, stereoco;
LocalVector monoup, stereoup;
LocalPoint pmonofwd, pstereofwd;
LocalPoint pmonobwd, pstereobwd;
LocalPoint pmonoco, pstereoco;
LocalPoint pmonoup, pstereoup;
if(matchedhit){
type=Matched;
GluedGeomDet * gdet=(GluedGeomDet *)tracker->idToDet(matchedhit->geographicalId());
GlobalVector gtrkdirfwd=gdet->toGlobal(fwdptsod.parameters().momentum());
GlobalVector gtrkdirbwd=gdet->toGlobal(bwdptsod.parameters().momentum());
GlobalVector gtrkdirup=gdet->toGlobal(updatedptsod.parameters().momentum());
GlobalVector gtrkdirco=gdet->toGlobal(combinedptsod.parameters().momentum());
const GeomDetUnit * monodet=gdet->monoDet();
monofwd=monodet->toLocal(gtrkdirfwd);
monobwd=monodet->toLocal(gtrkdirbwd);
monoup=monodet->toLocal(gtrkdirup);
monoco=monodet->toLocal(gtrkdirco);
pmonofwd=project(gdet,monodet,fwdptsod.parameters().position(),monofwd);
pmonobwd=project(gdet,monodet,bwdptsod.parameters().position(),monobwd);
pmonoup=project(gdet,monodet,updatedptsod.parameters().position(),monoup);
pmonoco=project(gdet,monodet,combinedptsod.parameters().position(),monoco);
const GeomDetUnit * stereodet=gdet->stereoDet();
stereofwd=stereodet->toLocal(gtrkdirfwd);
stereobwd=stereodet->toLocal(gtrkdirbwd);
stereoup=stereodet->toLocal(gtrkdirup);
stereoco=stereodet->toLocal(gtrkdirco);
pstereofwd=project(gdet,stereodet,fwdptsod.parameters().position(),stereofwd);
pstereobwd=project(gdet,stereodet,bwdptsod.parameters().position(),stereobwd);
pstereoup=project(gdet,stereodet,updatedptsod.parameters().position(),stereoup);
pstereoco=project(gdet,stereodet,combinedptsod.parameters().position(),stereoco);
}
else if(phit){
type=Projected;
GluedGeomDet * gdet=(GluedGeomDet *)tracker->idToDet(phit->geographicalId());
GlobalVector gtrkdirfwd=gdet->toGlobal(fwdptsod.parameters().momentum());
GlobalVector gtrkdirbwd=gdet->toGlobal(bwdptsod.parameters().momentum());
GlobalVector gtrkdirup=gdet->toGlobal(updatedptsod.parameters().momentum());
GlobalVector gtrkdirco=gdet->toGlobal(combinedptsod.parameters().momentum());
const SiStripRecHit2D& originalhit=phit->originalHit();
const GeomDetUnit * det;
if(!StripSubdetector(originalhit.geographicalId().rawId()).stereo()){
det=gdet->monoDet();
monofwd= det->toLocal(gtrkdirfwd);
monobwd= det->toLocal(gtrkdirbwd);
monoup= det->toLocal(gtrkdirup);
monoco= det->toLocal(gtrkdirco);
pmonofwd=project(gdet,det,fwdptsod.parameters().position(),monofwd);
pmonobwd=project(gdet,det,bwdptsod.parameters().position(),monobwd);
pmonoup=project(gdet,det,updatedptsod.parameters().position(),monoup);
pmonoco=project(gdet,det,combinedptsod.parameters().position(),monoco);
}
else{
det=gdet->stereoDet();
stereofwd= det->toLocal(gtrkdirfwd);
stereobwd= det->toLocal(gtrkdirbwd);
stereoup= det->toLocal(gtrkdirup);
stereoco= det->toLocal(gtrkdirco);
pstereofwd=project(gdet,det,fwdptsod.parameters().position(),stereofwd);
pstereobwd=project(gdet,det,bwdptsod.parameters().position(),stereobwd);
pstereoup=project(gdet,det,updatedptsod.parameters().position(),stereoup);
pstereoco=project(gdet,det,combinedptsod.parameters().position(),stereoco);
}
}
TrackingRecHitInfo::TrackingStates states;
if(forwardPredictedStateTag_!="") states.insert(std::make_pair(FwPredicted, TrackingStateInfo(std::make_pair(monofwd,stereofwd), std::make_pair(pmonofwd,pstereofwd), fwdptsod)));
if(backwardPredictedStateTag_!="")states.insert(std::make_pair(BwPredicted, TrackingStateInfo(std::make_pair(monobwd,stereobwd), std::make_pair(pmonobwd,pstereobwd), bwdptsod)));
if(updatedStateTag_!="")states.insert(std::make_pair(Updated, TrackingStateInfo(std::make_pair(monoup,stereoup), std::make_pair(pmonoup,pstereoup), updatedptsod)));
if(combinedStateTag_!="")states.insert(std::make_pair(Combined, TrackingStateInfo(std::make_pair(monoco,stereoco), std::make_pair(pmonoco,pstereoco), combinedptsod)));
TrackingRecHitInfo tkRecHitInfo(type, states);
if(j!=0){
trajinfo.insert(std::make_pair(thehitref,tkRecHitInfo));
}
// else edm::LogInfo("TrackInfoProducer") << "RecHit not associated ";
}
//edm::LogInfo("TrackInfoProducer") << "Found "<<nhit<< " hits";
//if(fwdtrajinfo.size()!=nhit) edm::LogInfo("TrackInfoProducer") << "Number of trackinfos "<<fwdtrajinfo.size()<< " doesn't match!";
output=TrackInfo((traj_iterator->seed()),trajinfo);
}
std::string TrackInfoProducerAlgorithm::backwardPredictedStateTag_ [private] |
Definition at line 40 of file TrackInfoProducerAlgorithm.h.
std::string TrackInfoProducerAlgorithm::combinedStateTag_ [private] |
Definition at line 40 of file TrackInfoProducerAlgorithm.h.
Definition at line 39 of file TrackInfoProducerAlgorithm.h.
std::string TrackInfoProducerAlgorithm::forwardPredictedStateTag_ [private] |
Definition at line 40 of file TrackInfoProducerAlgorithm.h.
std::string TrackInfoProducerAlgorithm::updatedStateTag_ [private] |
Definition at line 40 of file TrackInfoProducerAlgorithm.h.
1.7.3