|
|
|
#include "RecoTracker/CkfPattern/interface/GroupedCkfTrajectoryBuilder.h"#include "RecoTracker/CkfPattern/interface/TrajectorySegmentBuilder.h"#include "TrackingTools/DetLayers/interface/DetLayer.h"#include "TrackingTools/PatternTools/interface/TrajMeasLessEstim.h"#include "TrackingTools/KalmanUpdators/interface/KFUpdator.h"#include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"#include "TrackingTools/TrackFitters/interface/KFTrajectoryFitter.h"#include "RecoTracker/CkfPattern/interface/GroupedTrajCandLess.h"#include "TrackingTools/TrajectoryFiltering/interface/RegionalTrajectoryFilter.h"#include "TrackingTools/PatternTools/interface/TempTrajectory.h"#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"#include "TrackingTools/MeasurementDet/interface/LayerMeasurements.h"#include "TrackingTools/Records/interface/TrackingComponentsRecord.h"#include "TrackingTools/DetLayers/interface/DetGroup.h"#include "RecoTracker/TkDetLayers/interface/GeometricSearchTracker.h"#include "TrackingTools/Records/interface/TransientRecHitRecord.h"#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHitBuilder.h"#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"#include "TrackingTools/PatternTools/interface/TrajectoryStateUpdator.h"#include "TrackingTools/PatternTools/interface/Trajectory.h"#include "TrackingTools/TrajectoryState/interface/BasicSingleTrajectoryState.h"#include "FWCore/MessageLogger/interface/MessageLogger.h"#include "TrackingTools/PatternTools/interface/TransverseImpactPointExtrapolator.h"#include "TrackingTools/TrajectoryCleaning/interface/TrajectoryCleanerBySharedHits.h"#include <algorithm>#include "TrackingTools/DetLayers/interface/BarrelDetLayer.h"#include "TrackingTools/DetLayers/interface/ForwardDetLayer.h"Go to the source code of this file.
Functions | |
| std::string | whatIsTheNextStep (TempTrajectory &traj, std::pair< TrajectoryStateOnSurface, std::vector< const DetLayer * > > &stateAndLayers) |
| std::string | whatIsTheStateToUse (TrajectoryStateOnSurface &initial, TrajectoryStateOnSurface &stateToUse, const DetLayer *l) |
| std::string whatIsTheNextStep | ( | TempTrajectory & | traj, |
| std::pair< TrajectoryStateOnSurface, std::vector< const DetLayer * > > & | stateAndLayers | ||
| ) |
Definition at line 278 of file GroupedCkfTrajectoryBuilder.cc.
References TempTrajectory::lastLayer(), GloballyPositioned< T >::phi(), GloballyPositioned< T >::position(), Cylinder::radius(), ForwardDetLayer::specificSurface(), and BarrelDetLayer::specificSurface().
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer().
{
std::stringstream buffer;
vector<const DetLayer*> & nl = stateAndLayers.second;
#include "TrackingTools/DetLayers/interface/BarrelDetLayer.h"
#include "TrackingTools/DetLayers/interface/ForwardDetLayer.h"
//B.M. TkLayerName layerName;
//B.M. buffer << "Started from " << layerName(traj.lastLayer())
const BarrelDetLayer* sbdl = dynamic_cast<const BarrelDetLayer*>(traj.lastLayer());
const ForwardDetLayer* sfdl = dynamic_cast<const ForwardDetLayer*>(traj.lastLayer());
if (sbdl) {
buffer << "Started from " << traj.lastLayer() << " r=" << sbdl->specificSurface().radius()
<< " phi=" << sbdl->specificSurface().phi() << endl;
} else if (sfdl) {
buffer << "Started from " << traj.lastLayer() << " z " << sfdl->specificSurface().position().z()
<< " phi " << sfdl->specificSurface().phi() << endl;
}
buffer << "Trying to go to";
for ( vector<const DetLayer*>::iterator il=nl.begin();
il!=nl.end(); il++){
//B.M. buffer << " " << layerName(*il) << " " << *il << endl;
const BarrelDetLayer* bdl = dynamic_cast<const BarrelDetLayer*>(*il);
const ForwardDetLayer* fdl = dynamic_cast<const ForwardDetLayer*>(*il);
if (bdl) buffer << " r " << bdl->specificSurface().radius() << endl;
if (fdl) buffer << " z " << fdl->specificSurface().position().z() << endl;
//buffer << " " << *il << endl;
}
return buffer.str();
}
| std::string whatIsTheStateToUse | ( | TrajectoryStateOnSurface & | initial, |
| TrajectoryStateOnSurface & | stateToUse, | ||
| const DetLayer * | l | ||
| ) |
Definition at line 308 of file GroupedCkfTrajectoryBuilder.cc.
References TrajectoryStateOnSurface::charge(), TrajectoryStateOnSurface::curvilinearError(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), i, CurvilinearTrajectoryError::matrix(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), mathSSE::sqrt(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer().
{
std::stringstream buffer;
buffer << "GCTB: starting from "
<< " r / phi / z = " << stateToUse.globalPosition().perp()
<< " / " << stateToUse.globalPosition().phi()
<< " / " << stateToUse.globalPosition().z()
<< " , pt / phi / pz /charge = "
<< stateToUse.globalMomentum().perp() << " / "
<< stateToUse.globalMomentum().phi() << " / "
<< stateToUse.globalMomentum().z() << " / "
<< stateToUse.charge()
<< " for layer at "<< l << endl;
buffer << " errors:";
for ( int i=0; i<5; i++ ) buffer << " " << sqrt(stateToUse.curvilinearError().matrix()(i,i));
buffer << endl;
//buffer << "GCTB: starting from r / phi / z = " << initial.globalPosition().perp()
//<< " / " << initial.globalPosition().phi()
//<< " / " << initial.globalPosition().z() << " , pt / pz = "
//<< initial.globalMomentum().perp() << " / "
//<< initial.globalMomentum().z() << " for layer at "
//<< l << endl;
//buffer << " errors:";
//for ( int i=0; i<5; i++ ) buffer << " " << sqrt(initial.curvilinearError().matrix()(i,i));
//buffer << endl;
return buffer.str();
}
1.7.3