|
|
|
#include <TrackTransformerForGlobalCosmicMuons.h>
Definition at line 39 of file TrackTransformerForGlobalCosmicMuons.h.
| TrackTransformerForGlobalCosmicMuons::TrackTransformerForGlobalCosmicMuons | ( | const edm::ParameterSet & | parameterSet | ) |
Constructor.
Definition at line 42 of file TrackTransformerForGlobalCosmicMuons.cc.
References edm::ParameterSet::getParameter().
{
theTrackerRecHitBuilderName = parameterSet.getParameter<string>("TrackerRecHitBuilder");
theMuonRecHitBuilderName = parameterSet.getParameter<string>("MuonRecHitBuilder");
theRPCInTheFit = parameterSet.getParameter<bool>("RefitRPCHits");
theCacheId_TC = theCacheId_GTG = theCacheId_MG = theCacheId_TRH = 0;
theSkipStationDT = parameterSet.getParameter<int>("SkipStationDT");
theSkipStationCSC = parameterSet.getParameter<int>("SkipStationCSC");
theSkipWheelDT = parameterSet.getParameter<int>("SkipWheelDT");
theTrackerSkipSystem = parameterSet.getParameter<int>("TrackerSkipSystem");
theTrackerSkipSection = parameterSet.getParameter<int>("TrackerSkipSection");//layer, wheel, or disk depending on the system
}
| TrackTransformerForGlobalCosmicMuons::~TrackTransformerForGlobalCosmicMuons | ( | ) | [virtual] |
| ESHandle< TrajectoryFitter > TrackTransformerForGlobalCosmicMuons::fitter | ( | bool | up | ) | const |
the refitter used to refit the reco::Track
Definition at line 665 of file TkLasBeamFitter.cc.
References asciidump::at, TkLasBeamFitter::atFunction(), TkLasBeamFitter::fitBeamSplitters_, TkLasBeamFitter::gBSparam, TkLasBeamFitter::gHitZprime, fitWZ::gMinuit, TkLasBeam::isAlignmentTube(), TkLasBeam::isTecInternal(), TkLasBeamFitter::tecMinusFunction(), and TkLasBeamFitter::tecPlusFunction().
{
TGraphErrors *lasData = new TGraphErrors(gHitZprime.size(),
&(gHitZprime[0]), &(hitPhi[0]),
&(hitZprimeError[0]), &(hitPhiError[0]));
// do fit (R = entire range)
if(beam.isTecInternal(1)){
TF1 tecPlus("tecPlus", tecPlusFunction, zMin, zMax, nFitParams );
tecPlus.SetParameter( 1, 0 ); // slope
tecPlus.SetParameter( nFitParams - 1, 0 ); // BS
lasData->Fit("tecPlus", "R"); // "R", "RV" or "RQ"
}
else if(beam.isTecInternal(-1)){
TF1 tecMinus("tecMinus", tecMinusFunction, zMin, zMax, nFitParams );
tecMinus.SetParameter( 1, 0 ); // slope
tecMinus.SetParameter( nFitParams - 1, 0 ); // BS
lasData->Fit("tecMinus", "R");
}
else{
TF1 at("at", atFunction, zMin, zMax, nFitParams );
at.SetParameter( 1, 0 ); // slope
at.SetParameter( nFitParams - 1, 0 ); // BS
lasData->Fit("at","R");
}
// get values and errors for offset and slope
gMinuit->GetParameter(0, offset, offsetError);
gMinuit->GetParameter(1, slope, slopeError);
// additional AT parameters
// define param errors that are not used later
double bsAngleParamError(0.), phiAtMinusParamError(0.),
phiAtPlusParamError(0.), atThetaSplitParamError(0.);
if(beam.isAlignmentTube()){
gMinuit->GetParameter(2, phiAtMinusParam, phiAtMinusParamError);
gMinuit->GetParameter(3, phiAtPlusParam, phiAtPlusParamError);
gMinuit->GetParameter(4, atThetaSplitParam, atThetaSplitParamError);
}
// get Beam Splitter parameters
if(fitBeamSplitters_){
if(beam.isAlignmentTube() && hitsAtTecPlus == 0){
bsAngleParam = gBSparam;
}
else{
gMinuit->GetParameter( nFitParams - 1 , bsAngleParam, bsAngleParamError);
}
}
else{
bsAngleParam = gBSparam;
}
// fill covariance matrix
vector<double> vec( covMatrix.num_col() * covMatrix.num_col() );
gMinuit->mnemat( &vec[0], covMatrix.num_col() );
for(int col = 0; col < covMatrix.num_col(); col++){
for(int row = 0; row < covMatrix.num_col(); row++){
covMatrix[col][row] = vec[row + covMatrix.num_col()*col];
}
}
// compute correlation between parameters
// double corr01 = covMatrix[1][0]/(offsetError*slopeError);
delete lasData;
}
| TransientTrackingRecHit::ConstRecHitContainer TrackTransformerForGlobalCosmicMuons::getTransientRecHits | ( | const reco::TransientTrack & | track | ) | const |
Definition at line 108 of file TrackTransformerForGlobalCosmicMuons.cc.
References filterCSVwithJSON::copy, CSC(), CSCDetId, DetId::det(), PXFDetId::disk(), GeomDetEnumerators::DT, TIBDetId::layer(), TOBDetId::layer(), PXBDetId::layer(), LogTrace, DetId::Muon, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, reco::TransientTrack::recHitsBegin(), reco::TransientTrack::recHitsEnd(), dedefs::RPC, RPCDetId, DetId::subdetId(), sistripsummary::TEC, sistripsummary::TIB, sistripsummary::TID, sistripsummary::TOB, align::Tracker, DetId::Tracker, TIDDetId::wheel(), and TECDetId::wheel().
{
TransientTrackingRecHit::ConstRecHitContainer tkHits;
TransientTrackingRecHit::ConstRecHitContainer staHits;
for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit) {
if((*hit)->isValid()) {
if ( (*hit)->geographicalId().det() == DetId::Tracker && TrackerKeep((*hit)->geographicalId())) {
tkHits.push_back(theTrackerRecHitBuilder->build(&**hit));
} else if ( (*hit)->geographicalId().det() == DetId::Muon && MuonKeep((*hit)->geographicalId())){
if( (*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit){
LogTrace("Reco|TrackingTools|TrackTransformer") << "RPC Rec Hit discarged";
continue;
}
staHits.push_back(theMuonRecHitBuilder->build(&**hit));
}
}
}
if(staHits.empty()) return staHits;
/*
bool up = staHits.front()->globalPosition().y()>0 ? true : false;
if(up){
reverse(staHits.begin(),staHits.end());
reverse(tkHits.begin(),tkHits.end());
}
*/
copy(staHits.begin(),staHits.end(),back_inserter(tkHits));
for(TransientTrackingRecHit::ConstRecHitContainer::const_iterator hit = tkHits.begin();
hit !=tkHits.end(); ++hit){
DetId hitId = (*hit)->geographicalId();
GlobalPoint glbpoint = trackingGeometry()->idToDet(hitId)->position();
if(hitId.det() == DetId::Tracker) {
if (hitId.subdetId() == StripSubdetector::TIB )
LogTrace("TrackFitters") << glbpoint << " I am TIB " << TIBDetId(hitId).layer();
else if (hitId.subdetId() == StripSubdetector::TOB )
LogTrace("TrackFitters") << glbpoint << " I am TOB " << TOBDetId(hitId).layer();
else if (hitId.subdetId() == StripSubdetector::TEC )
LogTrace("TrackFitters") << glbpoint << " I am TEC " << TECDetId(hitId).wheel();
else if (hitId.subdetId() == StripSubdetector::TID )
LogTrace("TrackFitters") << glbpoint << " I am TID " << TIDDetId(hitId).wheel();
else if (hitId.subdetId() == StripSubdetector::TID )
LogTrace("TrackFitters") << glbpoint << " I am TID " << TIDDetId(hitId).wheel();
else if (hitId.subdetId() == (int) PixelSubdetector::PixelBarrel )
LogTrace("TrackFitters") << glbpoint << " I am PixBar " << PXBDetId(hitId).layer();
else if (hitId.subdetId() == (int) PixelSubdetector::PixelEndcap )
LogTrace("TrackFitters") << glbpoint << " I am PixFwd " << PXFDetId(hitId).disk();
else
LogTrace("TrackFitters") << " UNKNOWN TRACKER HIT TYPE ";
} else if(hitId.det() == DetId::Muon) {
if(hitId.subdetId() == MuonSubdetId::DT)
LogTrace("TrackFitters") << glbpoint << " I am DT " << DTWireId(hitId);
else if (hitId.subdetId() == MuonSubdetId::CSC )
LogTrace("TrackFitters") << glbpoint << " I am CSC " << CSCDetId(hitId);
else if (hitId.subdetId() == MuonSubdetId::RPC )
LogTrace("TrackFitters") << glbpoint << " I am RPC " << RPCDetId(hitId);
else
LogTrace("TrackFitters") << " UNKNOWN MUON HIT TYPE ";
} else
LogTrace("TrackFitters") << " UNKNOWN HIT TYPE ";
}
return tkHits;
}
| const MagneticField* TrackTransformerForGlobalCosmicMuons::magneticField | ( | ) | const [inline] |
the magnetic field
Definition at line 55 of file TrackTransformerForGlobalCosmicMuons.h.
References theMGField.
{return &*theMGField;}
| bool TrackTransformerForGlobalCosmicMuons::MuonKeep | ( | DetId | id | ) | const |
check (via options) if this is a muon rec hit for removal
Definition at line 306 of file TrackTransformerForGlobalCosmicMuons.cc.
References MuonSubdetId::CSC, MuonSubdetId::DT, DetId::Muon, MuonSubdetId::RPC, relativeConstraints::station, RPCDetId::station(), DTChamberId::station(), CSCDetId::station(), and DTChamberId::wheel().
{
if (id.det() != DetId::Muon) return false;
if (theSkipStationDT < 0 && theSkipStationCSC < 0) return true;
int station = -999;
int wheel = -999;
bool isRPC = false;
bool isDT = false;
bool isCSC = false;
if ( id.subdetId() == MuonSubdetId::DT ) {
DTChamberId did(id.rawId());
station = did.station();
wheel = did.wheel();
isDT = true;
} else if ( id.subdetId() == MuonSubdetId::CSC ) {
CSCDetId did(id.rawId());
station = did.station();
isCSC = true;
} else if ( id.subdetId() == MuonSubdetId::RPC ) {
RPCDetId rpcid(id.rawId());
station = rpcid.station();
isRPC = true;
}
if (isRPC && (station == theSkipStationCSC || station == theSkipStationDT)) return false;
if (isDT && station == theSkipStationDT ) return false;
if (isCSC && station == theSkipStationCSC ) return false;
if (isDT && theSkipWheelDT > -998 && wheel == theSkipWheelDT) return false;
return true;
}
| ESHandle< Propagator > TrackTransformerForGlobalCosmicMuons::propagator | ( | bool | up | ) | const [private] |
Definition at line 132 of file MuonServiceProxy.cc.
{
propagators::const_iterator prop = thePropagators.find(propagatorName);
if (prop == thePropagators.end()){
LogError("Muon|RecoMuon|MuonServiceProxy")
<< "MuonServiceProxy: propagator with name: "<< propagatorName <<" not found! Please load it in the MuonServiceProxy.cff";
return ESHandle<Propagator>(0);
}
return prop->second;
}
| void TrackTransformerForGlobalCosmicMuons::setServices | ( | const edm::EventSetup & | setup | ) | [virtual] |
set the services needed by the TrackTransformer
Implements TrackTransformerBase.
Definition at line 61 of file TrackTransformerForGlobalCosmicMuons.cc.
References edm::EventSetup::get(), LogTrace, and metname.
{
const std::string metname = "Reco|TrackingTools|TrackTransformer";
setup.get<TrajectoryFitter::Record>().get("KFFitterForRefitInsideOut",theFitterIO);
setup.get<TrajectoryFitter::Record>().get("KFSmootherForRefitInsideOut",theSmootherIO);
setup.get<TrajectoryFitter::Record>().get("KFFitterForRefitOutsideIn",theFitterOI);
setup.get<TrajectoryFitter::Record>().get("KFSmootherForRefitOutsideIn",theSmootherOI);
unsigned long long newCacheId_TC = setup.get<TrackingComponentsRecord>().cacheIdentifier();
if ( newCacheId_TC != theCacheId_TC ){
LogTrace(metname) << "Tracking Component changed!";
theCacheId_TC = newCacheId_TC;
setup.get<TrackingComponentsRecord>().get("SmartPropagatorRK",thePropagatorIO);
setup.get<TrackingComponentsRecord>().get("SmartPropagatorRKOpposite",thePropagatorOI);
}
// Global Tracking Geometry
unsigned long long newCacheId_GTG = setup.get<GlobalTrackingGeometryRecord>().cacheIdentifier();
if ( newCacheId_GTG != theCacheId_GTG ) {
LogTrace(metname) << "GlobalTrackingGeometry changed!";
theCacheId_GTG = newCacheId_GTG;
setup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);
}
// Magfield Field
unsigned long long newCacheId_MG = setup.get<IdealMagneticFieldRecord>().cacheIdentifier();
if ( newCacheId_MG != theCacheId_MG ) {
LogTrace(metname) << "Magnetic Field changed!";
theCacheId_MG = newCacheId_MG;
setup.get<IdealMagneticFieldRecord>().get(theMGField);
}
// Transient Rechit Builders
unsigned long long newCacheId_TRH = setup.get<TransientRecHitRecord>().cacheIdentifier();
if ( newCacheId_TRH != theCacheId_TRH ) {
theCacheId_TRH = newCacheId_TRH;
LogTrace(metname) << "TransientRecHitRecord changed!";
setup.get<TransientRecHitRecord>().get(theTrackerRecHitBuilderName,theTrackerRecHitBuilder);
setup.get<TransientRecHitRecord>().get(theMuonRecHitBuilderName,theMuonRecHitBuilder);
}
}
| ESHandle< TrajectorySmoother > TrackTransformerForGlobalCosmicMuons::smoother | ( | bool | up | ) | const |
the smoother used to smooth the trajectory which came from the refitting step
Definition at line 186 of file TrackTransformerForGlobalCosmicMuons.cc.
{
if(up) return theSmootherOI;
else return theSmootherIO;
}
| bool TrackTransformerForGlobalCosmicMuons::TrackerKeep | ( | DetId | id | ) | const |
check (via options) if this is a tracker rec hit for removal
Definition at line 256 of file TrackTransformerForGlobalCosmicMuons.cc.
References PXFDetId::disk(), TIBDetId::layer(), TOBDetId::layer(), PXBDetId::layer(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, DetId::Tracker, TIDDetId::wheel(), and TECDetId::wheel().
{
bool retVal = true;
if (id.det() != DetId::Tracker ) return false;
if (theTrackerSkipSystem < 0 ) return true;
int layer = -999;
int disk = -999;
int wheel = -999;
if ( id.subdetId() == theTrackerSkipSystem){
if (theTrackerSkipSystem == PixelSubdetector::PixelBarrel) {
PXBDetId did(id.rawId());
layer = did.layer();
}
if (theTrackerSkipSystem == StripSubdetector::TIB) {
TIBDetId did(id.rawId());
layer = did.layer();
}
if (theTrackerSkipSystem == StripSubdetector::TOB) {
TOBDetId did(id.rawId());
layer = did.layer();
}
if (theTrackerSkipSystem == PixelSubdetector::PixelEndcap) {
PXFDetId did(id.rawId());
disk = did.disk();
}
if (theTrackerSkipSystem == StripSubdetector::TID) {
TIDDetId did(id.rawId());
wheel = did.wheel();
}
if (theTrackerSkipSystem == StripSubdetector::TEC) {
TECDetId did(id.rawId());
wheel = did.wheel();
}
}
if (theTrackerSkipSection > -998 && layer == theTrackerSkipSection) retVal = false;
if (theTrackerSkipSection > -998 && disk == theTrackerSkipSection) retVal = false;
if (theTrackerSkipSection > -998 && wheel == theTrackerSkipSection) retVal = false;
return retVal;
}
| edm::ESHandle<GlobalTrackingGeometry> TrackTransformerForGlobalCosmicMuons::trackingGeometry | ( | ) | const [inline] |
the tracking geometry
Definition at line 58 of file TrackTransformerForGlobalCosmicMuons.h.
References theTrackingGeometry.
{return theTrackingGeometry;}
| vector< Trajectory > TrackTransformerForGlobalCosmicMuons::transform | ( | const reco::Track & | tr | ) | const [virtual] |
Convert a reco::Track into Trajectory.
Convert Tracks into Trajectories.
Implements TrackTransformerBase.
Definition at line 199 of file TrackTransformerForGlobalCosmicMuons.cc.
References alongMomentum, reco::TransientTrack::innermostMeasurementState(), LogTrace, metname, oppositeToMomentum, reco::TransientTrack::outermostMeasurementState(), LargeD0_PixelPairStep_cff::propagator, and reco::TransientTrack::track().
{
const std::string metname = "Reco|TrackingTools|TrackTransformer";
reco::TransientTrack track(tr,magneticField(),trackingGeometry());
// Build the transient Rechits
TransientTrackingRecHit::ConstRecHitContainer recHitsForReFit = getTransientRecHits(track);
if(recHitsForReFit.size() < 2) return vector<Trajectory>();
bool up = recHitsForReFit.back()->globalPosition().y()>0 ? true : false;
LogTrace(metname) << "Up ? " << up;
PropagationDirection propagationDirection = up ? oppositeToMomentum : alongMomentum;
TrajectoryStateOnSurface firstTSOS = up ? track.outermostMeasurementState() : track.innermostMeasurementState();
unsigned int innerId = up ? track.track().outerDetId() : track.track().innerDetId();
LogTrace(metname) << "Prop Dir: " << propagationDirection << " FirstId " << innerId << " firstTSOS " << firstTSOS;
TrajectorySeed seed(PTrajectoryStateOnDet(),TrajectorySeed::recHitContainer(),propagationDirection);
if(recHitsForReFit.front()->geographicalId() != DetId(innerId)){
LogTrace(metname)<<"Propagation occurring"<<endl;
firstTSOS = propagator(up)->propagate(firstTSOS, recHitsForReFit.front()->det()->surface());
LogTrace(metname)<<"Final destination: " << recHitsForReFit.front()->det()->surface().position() << endl;
if(!firstTSOS.isValid()){
LogTrace(metname)<<"Propagation error!"<<endl;
return vector<Trajectory>();
}
}
vector<Trajectory> trajectories = fitter(up)->fit(seed,recHitsForReFit,firstTSOS);
if(trajectories.empty()){
LogTrace(metname)<<"No Track refitted!"<<endl;
return vector<Trajectory>();
}
Trajectory trajectoryBW = trajectories.front();
vector<Trajectory> trajectoriesSM = smoother(up)->trajectories(trajectoryBW);
if(trajectoriesSM.empty()){
LogTrace(metname)<<"No Track smoothed!"<<endl;
return vector<Trajectory>();
}
return trajectoriesSM;
}
unsigned long long TrackTransformerForGlobalCosmicMuons::theCacheId_GTG [private] |
Definition at line 88 of file TrackTransformerForGlobalCosmicMuons.h.
unsigned long long TrackTransformerForGlobalCosmicMuons::theCacheId_MG [private] |
Definition at line 89 of file TrackTransformerForGlobalCosmicMuons.h.
unsigned long long TrackTransformerForGlobalCosmicMuons::theCacheId_TC [private] |
Definition at line 87 of file TrackTransformerForGlobalCosmicMuons.h.
unsigned long long TrackTransformerForGlobalCosmicMuons::theCacheId_TRH [private] |
Definition at line 90 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 102 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 103 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 100 of file TrackTransformerForGlobalCosmicMuons.h.
Referenced by magneticField().
edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformerForGlobalCosmicMuons::theMuonRecHitBuilder [private] |
Definition at line 112 of file TrackTransformerForGlobalCosmicMuons.h.
std::string TrackTransformerForGlobalCosmicMuons::theMuonRecHitBuilderName [private] |
Definition at line 111 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 81 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 82 of file TrackTransformerForGlobalCosmicMuons.h.
bool TrackTransformerForGlobalCosmicMuons::theRPCInTheFit [private] |
Definition at line 92 of file TrackTransformerForGlobalCosmicMuons.h.
int TrackTransformerForGlobalCosmicMuons::theSkipStationCSC [private] |
Definition at line 94 of file TrackTransformerForGlobalCosmicMuons.h.
int TrackTransformerForGlobalCosmicMuons::theSkipStationDT [private] |
Definition at line 93 of file TrackTransformerForGlobalCosmicMuons.h.
int TrackTransformerForGlobalCosmicMuons::theSkipWheelDT [private] |
Definition at line 95 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 105 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 106 of file TrackTransformerForGlobalCosmicMuons.h.
edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformerForGlobalCosmicMuons::theTrackerRecHitBuilder [private] |
Definition at line 109 of file TrackTransformerForGlobalCosmicMuons.h.
std::string TrackTransformerForGlobalCosmicMuons::theTrackerRecHitBuilderName [private] |
Definition at line 108 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 97 of file TrackTransformerForGlobalCosmicMuons.h.
Definition at line 96 of file TrackTransformerForGlobalCosmicMuons.h.
edm::ESHandle<GlobalTrackingGeometry> TrackTransformerForGlobalCosmicMuons::theTrackingGeometry [private] |
Definition at line 99 of file TrackTransformerForGlobalCosmicMuons.h.
Referenced by trackingGeometry().
1.7.3