|
|
|
#include <InOutConversionSeedFinder.h>
Definition at line 39 of file InOutConversionSeedFinder.h.
typedef FreeTrajectoryState InOutConversionSeedFinder::FTS [private] |
Definition at line 44 of file InOutConversionSeedFinder.h.
typedef TrajectoryStateOnSurface InOutConversionSeedFinder::TSOS [private] |
Definition at line 45 of file InOutConversionSeedFinder.h.
| InOutConversionSeedFinder::InOutConversionSeedFinder | ( | const edm::ParameterSet & | config | ) |
Definition at line 19 of file InOutConversionSeedFinder.cc.
References conf_, edm::ParameterSet::getParameter(), maxNumberOfInOutSeedsPerInputTrack_, the2ndHitdphi_, the2ndHitdzConst_, and the2ndHitdznSigma_.
: ConversionSeedFinder( conf ), conf_(conf) { // std::cout << " InOutConversionSeedFinder CTOR " << "\n"; maxNumberOfInOutSeedsPerInputTrack_ = conf_.getParameter<int>("maxNumOfSeedsInOut"); //the2ndHitdphi_ = 0.008; the2ndHitdphi_ = 0.01; the2ndHitdzConst_ = 5.; the2ndHitdznSigma_ = 2.; }
| InOutConversionSeedFinder::~InOutConversionSeedFinder | ( | ) | [virtual] |
Definition at line 37 of file InOutConversionSeedFinder.cc.
{
//std::cout << " InOutConversionSeedFinder DTOR " << "\n";
}
| void InOutConversionSeedFinder::completeSeed | ( | const TrajectoryMeasurement & | m1, |
| FreeTrajectoryState & | fts, | ||
| const Propagator * | propagator, | ||
| int | ilayer | ||
| ) | const [private] |
} / end debug
Definition at line 591 of file InOutConversionSeedFinder.cc.
References Reference_intrackfit_cff::barrel, createSeed(), ConversionSeedFinder::getMeasurementTracker(), i, ConversionSeedFinder::layerList(), DetLayer::location(), LayerMeasurements::measurements(), ConversionSeedFinder::printLayer(), TrajectoryMeasurement::recHit(), mathSSE::sqrt(), GeometricSearchDet::surface(), the2ndHitdphi_, the2ndHitdzConst_, and the2ndHitdznSigma_.
Referenced by findSeeds().
{
//std::cout<< "InOutConversionSeedFinder::completeSeed ilayer " << ilayer << "\n";
// A seed is made from 2 Trajectory Measuremennts. The 1st is the input
// argument m1. This routine looks for the 2nd measurement in layer ilayer
// Begin by making a new much stricter MeasurementEstimator based on the
// position errors of the 1st hit.
printLayer(ilayer);
MeasurementEstimator * newEstimator;
std::vector<const DetLayer*> allLayers=layerList();
const DetLayer * layer = allLayers[ilayer];
// if ( layer->location() == GeomDetEnumerators::barrel ) {const BarrelDetLayer * barrelLayer = dynamic_cast<const BarrelDetLayer*>(layer);
// //std::cout << "InOutConversionSeedFinder::completeSeed **** Barrel on layer " << ilayer << " R= " << barrelLayer->specificSurface().radius() << "\n";
// } else {
// const ForwardDetLayer * forwardLayer = dynamic_cast<const ForwardDetLayer*>(layer);
// //std::cout << "InOutConversionSeedFinder::completeSeed **** Forw on layer " << ilayer << " Z= " << forwardLayer->specificSurface().position().z() << "\n";
if (layer->location() == GeomDetEnumerators::barrel ) {
float dz = sqrt(the2ndHitdznSigma_*the2ndHitdznSigma_*m1.recHit()->globalPositionError().czz()
+ the2ndHitdzConst_*the2ndHitdzConst_);
newEstimator = new ConversionBarrelEstimator(-the2ndHitdphi_, the2ndHitdphi_, -dz, dz);
}
else {
float m1dr = sqrt(m1.recHit()->localPositionError().yy());
float dr = sqrt(the2ndHitdznSigma_*the2ndHitdznSigma_*m1dr*m1dr
+ the2ndHitdzConst_*the2ndHitdznSigma_);
newEstimator = new ConversionForwardEstimator(-the2ndHitdphi_, the2ndHitdphi_, dr);
}
//std::cout << "InOutConversionSeedFinder::completeSeed fts For the TSOS " << fts << "\n";
TSOS tsos(fts, layer->surface() );
if ( !tsos.isValid() ) {
//std::cout << "InOutConversionSeedFinder::completeSeed TSOS is not valid " << "\n";
}
//std::cout << "InOutConversionSeedFinder::completeSeed TSOS " << tsos << "\n";
//std::cout << "InOutConversionSeedFinder::completeSeed propagationDirection " << int(propagator->propagationDirection() ) << "\n";
//std::cout << "InOutConversionSeedFinder::completeSeed pointer to estimator " << newEstimator << "\n";
LayerMeasurements theLayerMeasurements_(this->getMeasurementTracker() );
std::vector<TrajectoryMeasurement> measurements = theLayerMeasurements_.measurements( *layer, tsos, *propagator, *newEstimator);
//std::cout << "InOutConversionSeedFinder::completeSeed Found " << measurements.size() << " second hits " << "\n";
delete newEstimator;
for(unsigned int i = 0; i < measurements.size(); ++i) {
if( measurements[i].recHit()->isValid() ) {
createSeed(m1, measurements[i]);
}
}
}
| void InOutConversionSeedFinder::createSeed | ( | const TrajectoryMeasurement & | m1, |
| const TrajectoryMeasurement & | m2 | ||
| ) | const [private] |
Definition at line 664 of file InOutConversionSeedFinder.cc.
References alongMomentum, TrajectoryStateOnSurface::curvilinearError(), benchmark_cfg::errors, TrajectoryMeasurement::estimate(), TrajectoryStateOnSurface::freeTrajectoryState(), TrajectoryStateOnSurface::isValid(), TrajectoryMeasurement::layer(), maxNumberOfInOutSeedsPerInputTrack_, nSeedsPerInputTrack_, trajectoryStateTransform::persistentState(), TrajectoryMeasurement::predictedState(), Propagator::propagate(), edm::OwnVector< T, P >::push_back(), TrajectoryMeasurement::recHit(), ConversionSeedFinder::thePropagatorAlongMomentum_, ConversionSeedFinder::theSeeds_, ConversionSeedFinder::theUpdator_, track2Charge_, and track2InitialMomentum_.
Referenced by completeSeed().
{
//std::cout << "InOutConversionSeedFinder::createSeed " << "\n";
if ( m1.predictedState().isValid() ) {
GlobalTrajectoryParameters newgtp( m1.recHit()->globalPosition(), track2InitialMomentum_, track2Charge_, &(*theMF_) );
CurvilinearTrajectoryError errors = m1.predictedState().curvilinearError();
FreeTrajectoryState fts(newgtp, errors);
TrajectoryStateOnSurface state1 = thePropagatorAlongMomentum_->propagate(fts, m1.recHit()->det()->surface());
/*
//std::cout << "hit surface " << m1.recHit()->det()->surface().position() << "\n";
//std::cout << "prop to " << typeid( m1.recHit()->det()->surface() ).name() <<"\n";
//std::cout << "prop to first hit " << state1 << "\n";
//std::cout << "update to " << m1.recHit()->globalPosition() << "\n";
*/
if ( state1.isValid() ) {
TrajectoryStateOnSurface updatedState1 = theUpdator_.update(state1, *m1.recHit() );
if ( updatedState1.isValid() ) {
TrajectoryStateOnSurface state2 = thePropagatorAlongMomentum_->propagate(*updatedState1.freeTrajectoryState(), m2.recHit()->det()->surface());
if ( state2.isValid() ) {
TrajectoryStateOnSurface updatedState2 = theUpdator_.update(state2, *m2.recHit() );
TrajectoryMeasurement meas1(state1, updatedState1, m1.recHit() , m1.estimate(), m1.layer());
TrajectoryMeasurement meas2(state2, updatedState2, m2.recHit() , m2.estimate(), m2.layer());
edm::OwnVector<TrackingRecHit> myHits;
myHits.push_back(meas1.recHit()->hit()->clone());
myHits.push_back(meas2.recHit()->hit()->clone());
//std::cout << "InOutConversionSeedFinder::createSeed new seed " << "\n";
if ( nSeedsPerInputTrack_ >= maxNumberOfInOutSeedsPerInputTrack_ ) return;
PTrajectoryStateOnDet const & ptsod= trajectoryStateTransform::persistentState(state2, meas2.recHit()->hit()->geographicalId().rawId() );
//std::cout << " InOutConversionSeedFinder::createSeed New seed parameters " << state2 << "\n";
theSeeds_.push_back(TrajectorySeed( ptsod, myHits, alongMomentum ));
nSeedsPerInputTrack_++;
//std::cout << "InOutConversionSeedFinder::createSeed New seed hit 1 R " << m1.recHit()->globalPosition().perp() << "\n";
//std::cout << "InOutConversionSeedFinder::createSeed New seed hit 2 R " << m2.recHit()->globalPosition().perp() << "\n";
}
}
}
}
}
| void InOutConversionSeedFinder::fillClusterSeeds | ( | ) | const [private, virtual] |
Definition at line 66 of file InOutConversionSeedFinder.cc.
References GeomDetEnumerators::barrel, TrajectoryStateOnSurface::freeTrajectoryState(), i, TrackingGeometry::idToDet(), TrajectoryStateOnSurface::isValid(), ConversionSeedFinder::layerList(), myPointer, nSeedsPerInputTrack_, ConversionSeedFinder::printLayer(), Propagator::propagate(), TrajectoryMeasurement::recHit(), startSeed(), GeometricSearchDet::surface(), theOutInTracks_, ConversionSeedFinder::thePropagatorOppositeToMomentum_, ConversionSeedFinder::theTrackerGeom_, and TrajectoryMeasurement::updatedState().
Referenced by makeSeeds().
{
std::vector<Trajectory>::const_iterator outInTrackItr;
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds outInTracks_.size " << theOutInTracks_.size() << "\n";
//Start looking for seeds for both of the 2 best tracks from the inward tracking
/*
for(outInTrackItr = theOutInTracks_.begin(); outInTrackItr != theOutInTracks_.end(); ++outInTrackItr) {
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds out in input track hits " << (*outInTrackItr).foundHits() << "\n";
DetId tmpId = DetId( (*outInTrackItr).seed().startingState().detId());
const GeomDet* tmpDet = this->getMeasurementTracker()->geomTracker()->idToDet( tmpId );
GlobalVector gv = tmpDet->surface().toGlobal( (*outInTrackItr).seed().startingState().parameters().momentum() );
//std::cout << " InOutConversionSeedFinder::fillClusterSeed was built from seed position " <<gv << " charge " << (*outInTrackItr).seed().startingState().parameters().charge() << "\n";
Trajectory::DataContainer m= outInTrackItr->measurements();
int nHit=0;
for (Trajectory::DataContainer::iterator itm = m.begin(); itm != m.end(); ++itm) {
if ( itm->recHit()->isValid() ) {
nHit++;
//std::cout << nHit << ") Valid RecHit global position " << itm->recHit()->globalPosition() << " R " << itm->recHit()->globalPosition().perp() << " phi " << itm->recHit()->globalPosition().phi() << " eta " << itm->recHit()->globalPosition().eta() << "\n";
}
}
}
*/
//Start looking for seeds for both of the 2 best tracks from the inward tracking
for(outInTrackItr = theOutInTracks_.begin(); outInTrackItr != theOutInTracks_.end(); ++outInTrackItr) {
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds out in input track hits " << (*outInTrackItr).foundHits() << "\n";
nSeedsPerInputTrack_=0;
//Find the first valid hit of the track
// Measurements are ordered according to the direction in which the trajectories were built
std::vector<TrajectoryMeasurement> measurements = (*outInTrackItr).measurements();
std::vector<const DetLayer*> allLayers=layerList();
//std::cout << " InOutConversionSeedFinder::fill clusterSeed allLayers.size " << allLayers.size() << "\n";
for(unsigned int i = 0; i < allLayers.size(); ++i) {
//std::cout << " allLayers " << allLayers[i] << "\n";
printLayer(i);
}
std::vector<const DetLayer*> myLayers;
myLayers.clear();
std::vector<TrajectoryMeasurement>::reverse_iterator measurementItr;
std::vector<TrajectoryMeasurement*> myItr;
// TrajectoryMeasurement* myPointer=0;
myPointer=0;
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds measurements.size " << measurements.size() <<"\n";
for(measurementItr = measurements.rbegin() ; measurementItr != measurements.rend(); ++measurementItr) {
if( (*measurementItr).recHit()->isValid()) {
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds measurement on layer " << measurementItr->layer() << " " <<&(*measurementItr) << " position " << measurementItr->recHit()->globalPosition() << " R " << sqrt( measurementItr->recHit()->globalPosition().x()*measurementItr->recHit()->globalPosition().x() + measurementItr->recHit()->globalPosition().y()*measurementItr->recHit()->globalPosition().y() ) << " Z " << measurementItr->recHit()->globalPosition().z() << " phi " << measurementItr->recHit()->globalPosition().phi() << "\n";
myLayers.push_back( measurementItr->layer() ) ;
myItr.push_back( &(*measurementItr) );
}
}
//std::cout << " InOutConversionSeedFinder::fillClusterSeed myLayers.size " << myLayers.size() << "\n";
// for( unsigned int i = 0; i < myLayers.size(); ++i) {
// }
if ( myItr.size()==0 ) {
//std::cout << "HORRENDOUS ERROR! No meas on track!" << "\n";
}
unsigned int ilayer;
for(ilayer = 0; ilayer < allLayers.size(); ++ilayer) {
//std::cout << " allLayers in the search loop " << allLayers[ilayer] << " " << myLayers[0] << "\n";
if ( allLayers[ilayer] == myLayers[0]) {
myPointer=myItr[0];
//std::cout << " allLayers in the search loop allLayers[ilayer] == myLayers[0]) " << allLayers[ilayer] << " " << myLayers[0] << " myPointer " << myPointer << "\n";
//std::cout << "Layer " << ilayer << " contains the first valid measurement " << "\n";
printLayer(ilayer);
if ( (myLayers[0])->location() == GeomDetEnumerators::barrel ) {
// const BarrelDetLayer * barrelLayer = dynamic_cast<const BarrelDetLayer*>(myLayers[0]);
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds **** firstHit found in Barrel on layer " << ilayer << " R= " << barrelLayer->specificSurface().radius() << "\n";
} else {
//const ForwardDetLayer * forwardLayer = dynamic_cast<const ForwardDetLayer*>(myLayers[0]);
//std::cout << " InOutwardConversionSeedFinder::fillClusterSeeds **** firstHit found in Forw on layer " << ilayer << " Z= " << forwardLayer->specificSurface().position().z() << "\n";
}
break;
} else if ( allLayers[ilayer] == myLayers[1] ) {
myPointer=myItr[1];
//std::cout << " allLayers in the search loop allLayers[ilayer] == myLayers[1]) " << allLayers[ilayer] << " " << myLayers[1] << " myPointer " << myPointer << "\n";
//std::cout << "Layer " << ilayer << " contains the first innermost valid measurement " << "\n";
if ( (myLayers[1])->location() == GeomDetEnumerators::barrel ) {
// const BarrelDetLayer * barrelLayer = dynamic_cast<const BarrelDetLayer*>(myLayers[1]);
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds **** 2ndHit found in Barrel on layer " << ilayer << " R= " << barrelLayer->specificSurface().radius() << "\n";
} else {
//const ForwardDetLayer * forwardLayer = dynamic_cast<const ForwardDetLayer*>(myLayers[1]);
//std::cout << " InOutwardConversionSeedFinder::fillClusterSeeds **** 2ndHitfound on forw layer " << ilayer << " Z= " << forwardLayer->specificSurface().position().z() << "\n";
}
break;
}
}
if(ilayer == allLayers.size()) {
//std::cout << "InOutConversionSeedFinder::fillClusterSeeds ERROR could not find layer on list" << "\n";
return;
}
//PropagatorWithMaterial reversePropagator(oppositeToMomentum, 0.000511, &(*theMF_) );
FreeTrajectoryState * fts = myPointer->updatedState().freeTrajectoryState();
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds First FTS charge " << fts->charge() << " Position " << fts->position() << " momentum " << fts->momentum() << " R " << sqrt(fts->position().x()*fts->position().x() + fts->position().y()* fts->position().y() ) << " Z " << fts->position().z() << " phi " << fts->position().phi() << " fts parameters " << fts->parameters() << "\n";
while (ilayer > 0) {
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds looking for 2nd seed from layer " << ilayer << "\n";
// if ( (allLayers[ilayer])->location() == GeomDetEnumerators::barrel ) {const BarrelDetLayer * barrelLayer = dynamic_cast<const BarrelDetLayer*>(allLayers[ilayer]);
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds **** Barrel on layer " << ilayer << " R= " << barrelLayer->specificSurface().radius() << "\n";
// } else {
//const ForwardDetLayer * forwardLayer = dynamic_cast<const ForwardDetLayer*>(allLayers[ilayer]);
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds **** Forw on layer " << ilayer << " Z= " << forwardLayer->specificSurface().position().z() << "\n";
// }
const DetLayer * previousLayer = allLayers[ilayer];
TrajectoryStateOnSurface stateAtPreviousLayer;
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds previousLayer->surface() position before " <<allLayers[ilayer] << " " << previousLayer->surface().position() << " layer location " << previousLayer->location() << "\n";
// Propagate to the previous layer
// The present layer is actually included in the loop so that a partner can be searched for
// Applying the propagator to the same layer does not do any harm. It simply does nothing
// const Propagator& newProp=thePropagatorOppositeToMomentum_;
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds reversepropagator direction " << thePropagatorOppositeToMomentum_->propagationDirection() << "\n";
if (ilayer-1>0) {
if ( allLayers[ilayer] == myLayers[0] ) {
//std::cout << " innermost hit R " << myPointer->recHit()->globalPosition().perp() << " Z " << myPointer->recHit()->globalPosition().z() << " phi " <<myPointer->recHit()->globalPosition().phi() << "\n";
//std::cout << " surface R " << theTrackerGeom_->idToDet( myPointer->recHit() ->geographicalId())->surface().position().perp() << " Z " << theTrackerGeom_->idToDet( myPointer->recHit() ->geographicalId())->surface().position().z() << " phi " << theTrackerGeom_->idToDet( myPointer->recHit() ->geographicalId())->surface().position().phi() << "\n";
stateAtPreviousLayer= thePropagatorOppositeToMomentum_->propagate(*fts, theTrackerGeom_->idToDet( myPointer->recHit() ->geographicalId())->surface() );
} else {
stateAtPreviousLayer= thePropagatorOppositeToMomentum_->propagate(*fts, previousLayer->surface() );
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds previousLayer->surface() position after " << previousLayer->surface().position() << " layer location " << previousLayer->location() << "\n";
}
} else if ( ilayer-1==0) {
//stateAtPreviousLayer= thePropagatorOppositeToMomentum_->propagate(*fts, theTrackerGeom_->idToDet( myPointer->recHit() ->geographicalId())->surface() );
stateAtPreviousLayer= thePropagatorOppositeToMomentum_->propagate(*fts, previousLayer->surface() );
}
if(!stateAtPreviousLayer.isValid()) {
//std::cout << "InOutConversionSeedFinder::fillClusterSeeds ERROR:could not propagate back to layer " << ilayer << "\n";
} else {
//std::cout << "InOutConversionSeedFinder::fillClusterSeeds stateAtPreviousLayer is valid. Propagating back to layer " << ilayer << "\n";
//std::cout << "InOutConversionSeedFinder::fillClusterSeeds stateAtPreviousLayer R " << stateAtPreviousLayer.globalPosition().perp() << " Z " << stateAtPreviousLayer.globalPosition().z() << " phi " << stateAtPreviousLayer.globalPosition().phi() << "\n";
startSeed(fts, stateAtPreviousLayer, -1, ilayer );
}
--ilayer;
}
if ( ilayer == 0) {
// if ( (allLayers[ilayer])->location() == GeomDetEnumerators::barrel ) {const BarrelDetLayer * barrelLayer = dynamic_cast<const BarrelDetLayer*>(allLayers[ilayer]);
// //std::cout << " InOutConversionSeedFinder::fillClusterSeeds **** Barrel on layer " << ilayer << " R= " << barrelLayer->specificSurface().radius() << "\n";
// } else {
//const ForwardDetLayer * forwardLayer = dynamic_cast<const ForwardDetLayer*>(allLayers[ilayer]);
//std::cout << " InOutConversionSeedFinder::fillClusterSeeds **** Forw on layer " << ilayer << " Z= " << forwardLayer->specificSurface().position().z() << "\n";
// }
const DetLayer * previousLayer = allLayers[ilayer];
TrajectoryStateOnSurface stateAtPreviousLayer;
stateAtPreviousLayer= thePropagatorOppositeToMomentum_->propagate(*fts, previousLayer->surface() );
if(!stateAtPreviousLayer.isValid()) {
//std::cout << "InOutConversionSeedFinder::fillClusterSeeds ERROR:could not propagate back to layer " << ilayer << "\n";
} else {
//std::cout << "InOutConversionSeedFinder::fillClusterSeeds stateAtPreviousLayer is valid. Propagating back to layer " << ilayer << "\n";
//std::cout << "InOutConversionSeedFinder::fillClusterSeeds stateAtPreviousLayer R " << stateAtPreviousLayer.globalPosition().perp() << " Z " << stateAtPreviousLayer.globalPosition().z() << " phi " << stateAtPreviousLayer.globalPosition().phi() << "\n";
startSeed(fts, stateAtPreviousLayer, -1, ilayer );
}
}
} // End loop over Out In tracks
}
| void InOutConversionSeedFinder::findSeeds | ( | const TrajectoryStateOnSurface & | startingState, |
| float | signedpt, | ||
| unsigned int | startingLayer | ||
| ) | const [private, virtual] |
Rememeber that this alwyas give back at least one dummy-innvalid it which prevents from everything getting stopped
Definition at line 408 of file InOutConversionSeedFinder.cc.
References Reference_intrackfit_cff::barrel, completeSeed(), dir, ConversionSeedFinder::getMeasurementTracker(), TrajectoryStateOnSurface::globalDirection(), TrajectoryStateOnSurface::globalPosition(), ConversionSeedFinder::layerList(), DetLayer::location(), m, myPointer, reco::CaloCluster::position(), TrajectoryMeasurement::recHit(), ConversionFastHelix::stateAtVertex(), GeometricSearchDet::surface(), theFirstMeasurements_, ConversionSeedFinder::thePropagatorAlongMomentum_, theSecondBC_, track2InitialMomentum_, x, detailsBasic3DVector::y, and z.
Referenced by startSeed().
{
std::vector<const DetLayer*> allLayers=layerList();
//std::cout << "InOutConversionSeedFinder::findSeeds starting forward propagation from startingLayer " << startingLayer << "\n";
// create error matrix
AlgebraicSymMatrix55 m = AlgebraicMatrixID();
m(0,0) = 0.1; m(1,1) = 0.0001 ; m(2,2) = 0.0001 ;
m(3,3) = 0.0001 ; m(4,4) = 0.001;
// Make an FTS consistent with the start point, start direction and curvature
FreeTrajectoryState fts(GlobalTrajectoryParameters(startingState.globalPosition(),
startingState.globalDirection(),
double(transverseCurvature), 0, &(*theMF_) ),
CurvilinearTrajectoryError(m));
if (fts.momentum().mag2() == 0){
edm::LogWarning("FailedToInitiateSeeding")<< " initial FTS has a zero momentum, probably because of the zero field. ";
return;
}
//std::cout << " InOutConversionSeedFinder::findSeeds startingState R "<< startingState.globalPosition().perp() << " Z " << startingState.globalPosition().z() << " phi " << startingState.globalPosition().phi() << " position " << startingState.globalPosition() << "\n";
//std::cout << " InOutConversionSeedFinder::findSeeds Initial FTS charge " << fts.charge() << " curvature " << transverseCurvature << "\n";
//std::cout << " InOutConversionSeedFinder::findSeeds Initial FTS parameters " << fts << "\n";
//float dphi = 0.01;
float dphi = 0.03;
float zrange = 5.;
for( unsigned int ilayer = startingLayer; ilayer <= startingLayer+1 && (ilayer < allLayers.size()-2); ++ilayer) {
const DetLayer * layer = allLayers[ilayer];
// if ( layer->location() == GeomDetEnumerators::barrel ) {const BarrelDetLayer * barrelLayer = dynamic_cast<const BarrelDetLayer*>(layer);
// //std::cout << "InOutConversionSeedFinder::findSeeds **** Barrel on layer " << ilayer << " R= " << barrelLayer->specificSurface().radius() << "\n";
// } else {
// const ForwardDetLayer * forwardLayer = dynamic_cast<const ForwardDetLayer*>(layer);
// }
// // end debug
MeasurementEstimator * newEstimator=0;
if (layer->location() == GeomDetEnumerators::barrel ) {
//std::cout << "InOutConversionSeedFinder::findSeeds Barrel ilayer " << ilayer << "\n";
newEstimator = new ConversionBarrelEstimator(-dphi, dphi, -zrange, zrange);
}
else {
//std::cout << "InOutConversionSeedFinder::findSeeds Forward ilayer " << ilayer << "\n";
newEstimator = new ConversionForwardEstimator(-dphi, dphi, 15.);
}
theFirstMeasurements_.clear();
// Get measurements compatible with the FTS and Estimator
TSOS tsos(fts, layer->surface() );
//std::cout << "InOutConversionSeedFinder::findSeed propagationDirection " << int(thePropagatorAlongMomentum_->propagationDirection() ) << "\n";
LayerMeasurements theLayerMeasurements_(this->getMeasurementTracker() );
theFirstMeasurements_ = theLayerMeasurements_.measurements( *layer, tsos, *thePropagatorAlongMomentum_, *newEstimator);
delete newEstimator;
//std::cout << "InOutConversionSeedFinder::findSeeds Found " << theFirstMeasurements_.size() << " first hits" << "\n";
if ( theFirstMeasurements_.size() == 1 ) { // only dummy hit found: start finding the seed from the innermost hit of the OutIn track
GlobalPoint bcPos((theSecondBC_.position()).x(),(theSecondBC_.position()).y(),(theSecondBC_.position()).z());
GlobalVector dir = startingState.globalDirection();
GlobalPoint back1mm = myPointer->recHit()->globalPosition();
back1mm -= dir.unit()*0.1;
//std::cout << " InOutConversionSeedFinder:::findSeeds going to make the helix using back1mm " << back1mm << "\n";
ConversionFastHelix helix(bcPos, myPointer->recHit()->globalPosition(), back1mm, &(*theMF_));
helix.stateAtVertex();
//std::cout << " InOutConversionSeedFinder:::findSeeds helix status " <<helix.isValid() << std::endl;
if ( !helix.isValid() ) continue;
track2InitialMomentum_= helix.stateAtVertex().momentum();
// Make a new FTS
FreeTrajectoryState newfts(GlobalTrajectoryParameters(
myPointer->recHit()->globalPosition(), startingState.globalDirection(),
helix.stateAtVertex().transverseCurvature(), 0, &(*theMF_)),
CurvilinearTrajectoryError(m));
completeSeed(*myPointer, newfts, thePropagatorAlongMomentum_, ilayer+1);
completeSeed(*myPointer, newfts, thePropagatorAlongMomentum_, ilayer+2);
} else {
//Loop over compatible hits
int mea=0;
for(std::vector<TrajectoryMeasurement>::iterator tmItr = theFirstMeasurements_.begin(); tmItr !=theFirstMeasurements_.end(); ++tmItr) {
mea++;
if (tmItr->recHit()->isValid() ) {
// Make a new helix as in fillClusterSeeds() but using the hit position
//std::cout << "InOutConversionSeedFinder::findSeeds hit R " << tmItr->recHit()->globalPosition().perp() << " Z " << tmItr->recHit()->globalPosition().z() << " " << tmItr->recHit()->globalPosition() << "\n";
GlobalPoint bcPos((theSecondBC_.position()).x(),(theSecondBC_.position()).y(),(theSecondBC_.position()).z());
GlobalVector dir = startingState.globalDirection();
GlobalPoint back1mm = tmItr->recHit()->globalPosition();
back1mm -= dir.unit()*0.1;
//std::cout << " InOutConversionSeedFinder:::findSeeds going to make the helix using back1mm " << back1mm << "\n";
ConversionFastHelix helix(bcPos, tmItr->recHit()->globalPosition(), back1mm, &(*theMF_));
helix.stateAtVertex();
//std::cout << " InOutConversionSeedFinder:::findSeeds helix status " <<helix.isValid() << std::endl;
if ( !helix.isValid() ) continue;
track2InitialMomentum_= helix.stateAtVertex().momentum();
//std::cout << "InOutConversionSeedFinder::findSeeds Updated estimatedPt = " << helix.stateAtVertex().momentum().perp() << " curvature " << helix.stateAtVertex().transverseCurvature() << "\n";
// << ", bcet = " << theBc->Et()
// << ", estimatedPt/bcet = " << estimatedPt/theBc->Et() << endl;
// Make a new FTS
FreeTrajectoryState newfts(GlobalTrajectoryParameters(
tmItr->recHit()->globalPosition(), startingState.globalDirection(),
helix.stateAtVertex().transverseCurvature(), 0, &(*theMF_)),
CurvilinearTrajectoryError(m));
//std::cout << "InOutConversionSeedFinder::findSeeds new FTS charge " << newfts.charge() << "\n";
/*
// Soome diagnostic output
// may be useful - comparission of the basic cluster position
// with the ecal impact position of the track
TrajectoryStateOnSurface stateAtECAL
= forwardPropagator.propagate(newfts, ECALSurfaces::barrel());
if (!stateAtECAL.isValid() || abs(stateAtECAL.globalPosition().eta())>1.479) {
if (startingState.globalDirection().eta() > 0.) {
stateAtECAL = forwardPropagator.propagate(newfts,
ECALSurfaces::positiveEtaEndcap());
} else {
stateAtECAL = forwardPropagator.propagate(newfts,
ECALSurfaces::negativeEtaEndcap());
}
}
GlobalPoint ecalImpactPosition = stateAtECAL.isValid() ? stateAtECAL.globalPosition() : GlobalPoint(0.,0.,0.);
cout << "Projected fts positon at ECAL surface: " <<
ecalImpactPosition << " bc position: " << theBc->Position() << endl;
*/
completeSeed(*tmItr, newfts, thePropagatorAlongMomentum_, ilayer+1);
completeSeed(*tmItr, newfts, thePropagatorAlongMomentum_, ilayer+2);
}
}
}
}
}
| std::vector< const reco::CaloCluster * > InOutConversionSeedFinder::getSecondCaloClusters | ( | const GlobalPoint & | conversionPosition, |
| float | charge | ||
| ) | const [private] |
Definition at line 365 of file InOutConversionSeedFinder.cc.
References bcCollection_, i, PV3DBase< T, PVType, FrameType >::phi(), and query::result.
Referenced by startSeed().
{
std::vector<const reco::CaloCluster*> result;
//std::cout << "InOutConversionSeedFinder::getSecondCaloClusters" << "\n";
Geom::Phi<float> theConvPhi(conversionPosition.phi() );
for (unsigned i = 0; i < bcCollection_->size(); ++i ) {
Geom::Phi<float> theBcPhi( bcCollection_->ptrAt(i)->position().phi() );
//std::cout<< "InOutConversionSeedFinder::getSecondCaloClusters BC energy " << bcCollection_->ptrAt(i)->energy() << " Calo cluster phi " << theBcPhi << " " << bcCollection_->ptrAt(i)->position().phi()<< " theConvPhi " << theConvPhi << "\n";
// Require phi of cluster to be consistent with the conversion
// position and the track charge
if (fabs(theBcPhi-theConvPhi ) < .5 &&
((charge<0 && theBcPhi-theConvPhi >-.5) ||
(charge>0 && theBcPhi-theConvPhi <.5))){
//result.push_back(&(*bcItr));
result.push_back(&(*(bcCollection_->ptrAt(i)) ));
}
}
return result;
}
| void InOutConversionSeedFinder::makeSeeds | ( | const edm::Handle< edm::View< reco::CaloCluster > > & | allBc | ) | const [virtual] |
Implements ConversionSeedFinder.
Definition at line 43 of file InOutConversionSeedFinder.cc.
References bcCollection_, fillClusterSeeds(), ConversionSeedFinder::findLayers(), inputTracks_, theFirstMeasurements_, theOutInTracks_, and ConversionSeedFinder::theSeeds_.
Referenced by ConversionTrackCandidateProducer::buildCollections().
{
//std::cout << " InOutConversionSeedFinder::makeSeeds() " << "\n";
theSeeds_.clear();
//std::cout << " Check Calo cluster collection size " << allBC->size() << "\n";
bcCollection_= allBC;
findLayers();
fillClusterSeeds();
//std::cout << "Built vector of seeds of size " << theSeeds_.size() << "\n" ;
theOutInTracks_.clear();
inputTracks_.clear();
theFirstMeasurements_.clear();
}
| void InOutConversionSeedFinder::setTracks | ( | std::vector< Trajectory > const & | in | ) | [inline] |
Definition at line 63 of file InOutConversionSeedFinder.h.
References recoMuon::in, and theOutInTracks_.
Referenced by ConversionTrackCandidateProducer::buildCollections().
{ theOutInTracks_ = in;}
| void InOutConversionSeedFinder::startSeed | ( | FreeTrajectoryState * | fts, |
| const TrajectoryStateOnSurface & | stateAtPreviousLayer, | ||
| int | charge, | ||
| int | layer | ||
| ) | const [private] |
Definition at line 311 of file InOutConversionSeedFinder.cc.
References FreeTrajectoryState::charge(), dir, findSeeds(), getSecondCaloClusters(), TrajectoryStateOnSurface::globalDirection(), TrajectoryStateOnSurface::globalPosition(), PV3DBase< T, PVType, FrameType >::phi(), reco::CaloCluster::position(), ConversionFastHelix::stateAtVertex(), theSecondBC_, track2Charge_, TrajectoryStateOnSurface::transverseCurvature(), x, detailsBasic3DVector::y, and z.
Referenced by fillClusterSeeds().
{
//std::cout << "InOutConversionSeedFinder::startSeed ilayer " << ilayer << "\n";
// Get a list of basic clusters that are consistent with a track
// starting at the assumed conversion point with opp. charge to the
// inward track. Loop over these basic clusters.
track2Charge_ = charge*fts->charge();
std::vector<const reco::CaloCluster*> bcVec;
//std::cout << "InOutConversionSeedFinder::startSeed charge assumed for the in-out track " << track2Charge_ << "\n";
Geom::Phi<float> theConvPhi( stateAtPreviousLayer.globalPosition().phi());
//std::cout << "InOutConversionSeedFinder::startSeed stateAtPreviousLayer phi " << stateAtPreviousLayer.globalPosition().phi() << " R " << stateAtPreviousLayer.globalPosition().perp() << " Z " << stateAtPreviousLayer.globalPosition().z() << "\n";
bcVec = getSecondCaloClusters(stateAtPreviousLayer.globalPosition(),track2Charge_);
std::vector<const reco::CaloCluster*>::iterator bcItr;
//std::cout << "InOutConversionSeedFinder::startSeed bcVec.size " << bcVec.size() << "\n";
// debug
// for(bcItr = bcVec.begin(); bcItr != bcVec.end(); ++bcItr) {
// //std::cout << "InOutConversionSeedFinder::startSeed list of bc eta " << (*bcItr)->position().eta() << " phi " << (*bcItr)->position().phi() << " x " << (*bcItr)->position().x() << " y " << (*bcItr)->position().y() << " z " << (*bcItr)->position().z() << "\n";
// }
for(bcItr = bcVec.begin(); bcItr != bcVec.end(); ++bcItr) {
theSecondBC_ = **bcItr;
GlobalPoint bcPos((theSecondBC_.position()).x(),
(theSecondBC_.position()).y(),
(theSecondBC_.position()).z());
//std::cout << "InOutConversionSeedFinder::startSeed for bc position x " << bcPos.x() << " y " << bcPos.y() << " z " << bcPos.z() << " eta " << bcPos.eta() << " phi " << bcPos.phi() << "\n";
GlobalVector dir = stateAtPreviousLayer.globalDirection();
GlobalPoint back1mm = stateAtPreviousLayer.globalPosition();
//std::cout << "InOutConversionSeedFinder::startSeed stateAtPreviousLayer.globalPosition() " << back1mm << "\n";
back1mm -= dir.unit()*0.1;
//std::cout << " InOutConversionSeedFinder:::startSeed going to make the helix using back1mm " << back1mm <<"\n";
ConversionFastHelix helix(bcPos, stateAtPreviousLayer.globalPosition(), back1mm, &(*theMF_));
helix.stateAtVertex();
//std::cout << " InOutConversionSeedFinder:::startSeed helix status " <<helix.isValid() << std::endl;
if ( !helix.isValid() ) continue;
findSeeds(stateAtPreviousLayer, helix.stateAtVertex().transverseCurvature(), ilayer);
}
}
edm::Handle<edm::View<reco::CaloCluster> > InOutConversionSeedFinder::bcCollection_ [mutable, private] |
Definition at line 96 of file InOutConversionSeedFinder.h.
Referenced by getSecondCaloClusters(), and makeSeeds().
Reimplemented from ConversionSeedFinder.
Definition at line 68 of file InOutConversionSeedFinder.h.
Referenced by InOutConversionSeedFinder().
std::vector<Trajectory> InOutConversionSeedFinder::inputTracks_ [mutable, private] |
Definition at line 91 of file InOutConversionSeedFinder.h.
Referenced by makeSeeds().
Definition at line 86 of file InOutConversionSeedFinder.h.
Referenced by createSeed(), and InOutConversionSeedFinder().
TrajectoryMeasurement* InOutConversionSeedFinder::myPointer [mutable, private] |
Definition at line 89 of file InOutConversionSeedFinder.h.
Referenced by fillClusterSeeds(), and findSeeds().
int InOutConversionSeedFinder::nSeedsPerInputTrack_ [mutable, private] |
Definition at line 85 of file InOutConversionSeedFinder.h.
Referenced by createSeed(), and fillClusterSeeds().
float InOutConversionSeedFinder::the2ndHitdphi_ [private] |
Definition at line 80 of file InOutConversionSeedFinder.h.
Referenced by completeSeed(), and InOutConversionSeedFinder().
float InOutConversionSeedFinder::the2ndHitdzConst_ [private] |
Definition at line 81 of file InOutConversionSeedFinder.h.
Referenced by completeSeed(), and InOutConversionSeedFinder().
float InOutConversionSeedFinder::the2ndHitdznSigma_ [private] |
Definition at line 82 of file InOutConversionSeedFinder.h.
Referenced by completeSeed(), and InOutConversionSeedFinder().
std::vector<TrajectoryMeasurement> InOutConversionSeedFinder::theFirstMeasurements_ [mutable, private] |
Definition at line 93 of file InOutConversionSeedFinder.h.
Referenced by findSeeds(), and makeSeeds().
std::vector<Trajectory> InOutConversionSeedFinder::theOutInTracks_ [mutable, private] |
Definition at line 92 of file InOutConversionSeedFinder.h.
Referenced by fillClusterSeeds(), makeSeeds(), and setTracks().
reco::CaloCluster InOutConversionSeedFinder::theSecondBC_ [mutable, private] |
Definition at line 95 of file InOutConversionSeedFinder.h.
Referenced by findSeeds(), and startSeed().
int InOutConversionSeedFinder::track2Charge_ [mutable, private] |
Definition at line 83 of file InOutConversionSeedFinder.h.
Referenced by createSeed(), and startSeed().
GlobalVector InOutConversionSeedFinder::track2InitialMomentum_ [mutable, private] |
Definition at line 84 of file InOutConversionSeedFinder.h.
Referenced by createSeed(), and findSeeds().
1.7.3