|
|
|
#include <HITrackVertexMaker.h>
Public Member Functions | |
| HITrackVertexMaker (const edm::ParameterSet &, const edm::EventSetup &es1) | |
| bool | produceTracks (const edm::Event &, const edm::EventSetup &, HICConst *, FmpConst *) |
| ~HITrackVertexMaker () | |
Private Attributes | |
| std::string | builderName |
| int | eventCount |
| edm::InputTag | L2candTag_ |
| edm::InputTag | primaryVertexTag |
| edm::ParameterSet | pset_ |
| edm::InputTag | rphirecHitsTag |
| std::vector< const NavigationSchool * > | theNavigationSchoolV |
Definition at line 62 of file HITrackVertexMaker.h.
| cms::HITrackVertexMaker::HITrackVertexMaker | ( | const edm::ParameterSet & | ps1, |
| const edm::EventSetup & | es1 | ||
| ) | [explicit] |
Definition at line 110 of file HITrackVertexMaker.cc.
References gather_cfg::cout, and edm::ParameterSet::getParameter().
{
L2candTag_ = ps1.getParameter< edm::InputTag > ("L2CandTag");
rphirecHitsTag = ps1.getParameter< edm::InputTag > ("rphiRecHits");
builderName = ps1.getParameter< std::string > ("TTRHBuilder");
primaryVertexTag = ps1.getParameter< edm::InputTag > ("PrimaryVertexTag");
#ifdef DEBUG
std::cout<<" Start HI TrackVertexMaker constructor "<<std::endl;
#endif
pset_ = ps1;
std::cout<<" No initialization "<<std::endl;
eventCount = 0;
#ifdef DEBUG
std::cout<<" HICTrajectoryBuilder constructed "<<std::endl;
#endif
}
| cms::HITrackVertexMaker::~HITrackVertexMaker | ( | ) |
Definition at line 132 of file HITrackVertexMaker.cc.
{
// std::cout<<" HITrackVertexMaker::destructor "<<std::endl;
}
| bool cms::HITrackVertexMaker::produceTracks | ( | const edm::Event & | e1, |
| const edm::EventSetup & | es1, | ||
| HICConst * | theHICConst, | ||
| FmpConst * | theFmpConst | ||
| ) |
Definition at line 137 of file HITrackVertexMaker.cc.
References alongMomentum, ClosestApproachInRPhi::calculate(), FreeTrajectoryState::charge(), gather_cfg::cout, cms::HICFTSfromL1orL2::createFTSfromL2(), FreeTrajectoryState::curvilinearError(), Reference_intrackfit_cff::endcap, first, TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::freeTrajectoryState(), edm::EventSetup::get(), edm::Event::getByLabel(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), i, iseed, TrajectoryStateOnSurface::isValid(), TrajectoryStateClosestToBeamLine::isValid(), TransientVertex::isValid(), j, HLT_ES_cff::magfield, makeMuonMisalignmentScenario::matrix, FreeTrajectoryState::momentum(), VarParsing::mult, TransientVertex::normalisedChiSquared(), TransientVertex::originalTracks(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), pos, TransientVertex::position(), FreeTrajectoryState::position(), cms::DiMuonSeedGeneratorHIC::produce(), edm::ESHandle< T >::product(), ptmin, edm::second(), HICMeasurementEstimator::setHICConst(), HICMeasurementEstimator::setMult(), cms::FmpConst::setVertex(), cms::HICConst::setVertex(), cms::HICTkOuterStartingLayerFinder::startingLayers(), evf::utils::state, TransientVertex::totalChiSquared(), patCandidatesForDimuonsSequences_cff::tracker, testEve_cfg::tracks, TrajectoryStateClosestToBeamLine::trackStateAtPCA(), reco::TrackBase::undefAlgorithm, v, KalmanVertexFitter::vertex(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and cms::HICConst::zvert.
Referenced by cms::TestMuL1L2Filter::filter(), cms::TestMuL1L2FilterSTA::filter(), and cms::HLTHIMuL1L2L3Filter::hltFilter().
{
bool dimuon = false;
edm::Handle<RecoChargedCandidateCollection> L2mucands;
// edm::Handle<TrackCollection> L2mucands;
e1.getByLabel (L2candTag_,L2mucands);
#ifdef DEBUG
cout<<" Number of muon candidates "<<L2mucands->size()<<endl;
if( L2mucands->size() < 2 ) cout<<"L2 muon failed"<<endl;
#endif
if( L2mucands->size() < 2 ) return dimuon;
#ifdef DEBUG
cout<<"L2 muon accepted"<<endl;
#endif
// std::cout<<" Just do nothing for L3 but initiate ESHandles "<<std::endl;
edm::Handle<reco::VertexCollection> vertexcands;
e1.getByLabel (primaryVertexTag,vertexcands);
#ifdef DEBUG
cout<<" Number of vertices primary "<<vertexcands->size()<<endl;
if(vertexcands->size()<1) cout<<" Primary vertex failed "<<endl;
#endif
if(vertexcands->size()<1) return dimuon;
#ifdef DEBUG_COUNT
cout<<" Accepted for L3 propagation "<<endl;
#endif
int iv = 0;
for (reco::VertexCollection::const_iterator ipvertex=vertexcands->begin();ipvertex!=vertexcands->end();ipvertex++)
{
// cout<<" Vertex position from pixels "<<(*ipvertex).position().z()<<endl;
if (iv == 0) {theHICConst->setVertex((*ipvertex).position().z()); theFmpConst->setVertex((*ipvertex).position().z());}
iv++;
}
// cout << " Vertex is set to (found by pixel finder)"<<theHICConst->zvert<<endl;
eventCount++;
// ============================ Event accepted for L3
// Initialization from Records
//
std::string updatorName = "KFUpdator";
std::string propagatorAlongName = "PropagatorWithMaterial";
std::string propagatorOppositeName = "PropagatorWithMaterialOpposite";
double theChiSquareCut = 500.;
double nsig = 3.;
double ptmin=1.;
edm::ESHandle<MagneticField> magfield;
edm::ESHandle<TransientTrackingRecHitBuilder> recHitBuilderHandle;
edm::ESHandle<MeasurementTracker> measurementTrackerHandle;
edm::ESHandle<GeometricSearchTracker> tracker;
edm::ESHandle<Propagator> propagatorAlongHandle;
edm::ESHandle<Propagator> propagatorOppositeHandle;
edm::ESHandle<TrajectoryStateUpdator> updatorHandle;
es1.get<TrackerRecoGeometryRecord>().get( tracker );
es1.get<IdealMagneticFieldRecord>().get(magfield);
es1.get<CkfComponentsRecord>().get("",measurementTrackerHandle);
es1.get<TransientRecHitRecord>().get(builderName,recHitBuilderHandle);
es1.get<TrackingComponentsRecord>().get(propagatorAlongName,propagatorAlongHandle);
es1.get<TrackingComponentsRecord>().get(propagatorOppositeName,propagatorOppositeHandle);
es1.get<TrackingComponentsRecord>().get(updatorName,updatorHandle);
// Initialization of navigation school
if(eventCount == 1) {
int lost=-1;
int lostf=-1;
theNavigationSchoolV.push_back(new HICSimpleNavigationSchool(&(*tracker), &(*magfield), lost, lostf));
for(int i=11;i>-1;i--) {
theNavigationSchoolV.push_back(new HICSimpleNavigationSchool(&(*tracker), &(*magfield), i, lostf));
}
lost=-1;
for(int i=12;i>-1;i--) {
theNavigationSchoolV.push_back(new HICSimpleNavigationSchool(&(*tracker), &(*magfield), lost, i));
}
} // initialization at the first event
MinPtTrajectoryFilter theMinPtFilter(ptmin);
HICMeasurementEstimator theEstimator(&(*tracker), &(*magfield), theChiSquareCut, nsig);
HICTrajectoryBuilder theTrajectoryBuilder( pset_, es1,
updatorHandle.product(),
propagatorAlongHandle.product(),
propagatorOppositeHandle.product(),
&theEstimator,
recHitBuilderHandle.product(),
measurementTrackerHandle.product(),
&theMinPtFilter);
measurementTrackerHandle->update(e1);
#ifdef DEBUG
std::cout<<" After first tracker update "<<std::endl;
#endif
// For trajectory builder
int theLowMult = 1;
theEstimator.setHICConst(theHICConst);
theEstimator.setMult(theLowMult);
theTrajectoryBuilder.settracker(measurementTrackerHandle.product());
//============================
// FastMuPropagator* theFmp = new FastMuPropagator(&(*magfield),theFmpConst);
// StateOnTrackerBound state(theFmp);
FastMuPropagator theFmp(&(*magfield),theFmpConst);
StateOnTrackerBound state(&theFmp);
TrajectoryStateOnSurface tsos;
HICFTSfromL1orL2 vFts(&(*magfield));
int NumOfSigma=4;
HICTkOuterStartingLayerFinder TkOSLF(NumOfSigma, &(*magfield), &(*tracker), theHICConst);
int mult = 1;
DiMuonSeedGeneratorHIC Seed(rphirecHitsTag,&(*magfield),&(*tracker), theHICConst, builderName, mult);
// vector<FreeTrajectoryState> theFts = vFts.createFTSfromStandAlone((*mucands));
vector<FreeTrajectoryState> theFts = vFts.createFTSfromL2((*L2mucands));
#ifdef DEBUG
cout<<" Size of the freeTS "<<theFts.size()<<endl;
#endif
DiMuonSeedGeneratorHIC::SeedContainer myseeds;
map<DetLayer*, DiMuonSeedGeneratorHIC::SeedContainer> seedmap;
vector<Trajectory> theTmpTrajectories0;
// map< FreeTrajectoryState*, Trajectory>
vector<FreeTrajectoryState*> theFoundFts;
map<FreeTrajectoryState*, vector<Trajectory> > theMapFtsTraj;
for(vector<FreeTrajectoryState>::iterator ifts=theFts.begin(); ifts!=theFts.end(); ifts++)
{
theTmpTrajectories0.clear();
#ifdef DEBUG
cout<<" cycle on Muon Trajectory State " <<(*ifts).parameters().position().perp()<<
" " <<(*ifts).parameters().position().z() <<endl;
#endif
tsos=state((*ifts));
if(tsos.isValid())
{
// vector<Trajectory> theTmpTrajectories0;
#ifdef DEBUG
cout<<" Position "<<tsos.globalPosition().perp()<<" "<<tsos.globalPosition().phi()<<
" "<<tsos.globalPosition().z()<<" "<<tsos.globalMomentum().perp()<<endl;
#endif
// Start to find starting layers
FreeTrajectoryState* ftsnew=tsos.freeTrajectoryState();
vector<DetLayer*> seedlayers = TkOSLF.startingLayers((*ftsnew));
#ifdef DEBUG
std::cout<<" The size of the starting layers "<<seedlayers.size()<<std::endl;
#endif
if( seedlayers.size() == 0 ) {
#ifdef DEBUG
cout<<" Starting layers failed for muon candidate "<<endl;
#endif
continue;
}
map<DetLayer*, DiMuonSeedGeneratorHIC::SeedContainer> seedmap = Seed.produce(e1 ,es1,
(*ftsnew), tsos, (*ifts),
recHitBuilderHandle.product(),
measurementTrackerHandle.product(),
&seedlayers);
for( vector<DetLayer*>::iterator it = seedlayers.begin(); it != seedlayers.end(); it++)
{
DiMuonSeedGeneratorHIC::SeedContainer seeds = (*seedmap.find(*it)).second;
#ifdef DEBUG
std::cout<<" Layer::Position z "<<(**it).surface().position().z()<<
" Number of seeds in layer "<<seeds.size()<<std::endl;
#endif
if(seeds.size() == 0) {
#ifdef DEBUG
std::cout<<" No seeds are found: do not continue with this Detlayer "<<std::endl;
#endif
continue;
}
// set the first navigation (all layers without gap)
NavigationSetter setter( *(theNavigationSchoolV[0]));
for(DiMuonSeedGeneratorHIC::SeedContainer::iterator iseed=seeds.begin();
iseed!=seeds.end();iseed++)
{
std::vector<TrajectoryMeasurement> theV = (*iseed).measurements();
#ifdef DEBUG
std::cout<< "RecHIT Layer position r "<<
theV[0].recHit()->globalPosition().perp()<<" phi "<<
theV[0].recHit()->globalPosition().phi()<<" z "<<
theV[0].recHit()->globalPosition().z()<<" momentum "<<
theV[0].updatedState().freeTrajectoryState()->parameters().momentum().perp()<<" "<<
theV[0].updatedState().freeTrajectoryState()->parameters().momentum().z()<<std::endl;
#endif
vector<Trajectory> theTmpTrajectories = theTrajectoryBuilder.trajectories(*iseed);
#ifdef DEBUG
cout<<" Number of found trajectories "<<theTmpTrajectories.size()<<endl;
#endif
if(theTmpTrajectories.size()>0) {
theTmpTrajectories0.insert(theTmpTrajectories0.end(),
theTmpTrajectories.begin(),
theTmpTrajectories.end());
#ifdef DEBUG
std::cout<<"We found trajectories for at least one seed "<<theTmpTrajectories0.size()<<std::endl;
#endif
break;
} // There are trajectories
} // seeds
if( theTmpTrajectories0.size() > 0 ) {
#ifdef DEBUG
std::cout<<"We found trajectories for at least one seed "<<theTmpTrajectories0.size()<<"break layer cycle "<<std::endl;
#endif
break;
}
} // seedlayer
if( theTmpTrajectories0.size() > 0 ) {
#ifdef DEBUG
std::cout<<"We found trajectories for at least one seed "<<theTmpTrajectories0.size()
<<"continue"<<std::endl;
#endif
theFoundFts.push_back(&(*ifts));
theMapFtsTraj[&(*ifts)] = theTmpTrajectories0;
continue;
}
#ifdef DEBUG
std::cout<<"No trajectories for this FTS "<<theTmpTrajectories0.size()<<
"try second path"<<std::endl;
#endif
// No trajectories found for this Muon FTS although seeds exist.
// Try to allow another trajectory map with lost 1 layer (only for barrel layers)
for( vector<DetLayer*>::iterator it = seedlayers.begin(); it != seedlayers.end(); it++)
{
DiMuonSeedGeneratorHIC::SeedContainer seeds = (*seedmap.find(*it)).second;
if(seeds.size() == 0) {
#ifdef DEBUG
std::cout<<" Second path: No seeds are found: do not continue with this Detlayer "<<std::endl;
#endif
continue;
}
if((**it).location() == GeomDetEnumerators::endcap) {
for(unsigned int j=13; j<theNavigationSchoolV.size();j++){
NavigationSetter setter( *(theNavigationSchoolV[j]));
for(DiMuonSeedGeneratorHIC::SeedContainer::iterator iseed=seeds.begin();
iseed!=seeds.end();iseed++)
{
std::vector<TrajectoryMeasurement> theV = (*iseed).measurements();
vector<Trajectory> theTmpTrajectories = theTrajectoryBuilder.trajectories(*iseed);
if(theTmpTrajectories.size()>0) {
theTmpTrajectories0.insert(theTmpTrajectories0.end(),
theTmpTrajectories.begin(),
theTmpTrajectories.end());
#ifdef DEBUG
std::cout<<"Second path: We found trajectories for at least one seed in barrel layer "<<
theTmpTrajectories0.size()<<std::endl;
#endif
break;
} // There are trajectories
} // seeds
if( theTmpTrajectories0.size() > 0 ) {
#ifdef DEBUG
std::cout<<"Second path: no trajectories: we try next barrel layer "<<
theTmpTrajectories0.size()<<std::endl;
#endif
break;
}
} // navigation maps
} else {
for(int j=1; j<13;j++){
NavigationSetter setter(*(theNavigationSchoolV[j]));
for(DiMuonSeedGeneratorHIC::SeedContainer::iterator iseed=seeds.begin();
iseed!=seeds.end();iseed++)
{
std::vector<TrajectoryMeasurement> theV = (*iseed).measurements();
vector<Trajectory> theTmpTrajectories = theTrajectoryBuilder.trajectories(*iseed);
if(theTmpTrajectories.size()>0) {
theTmpTrajectories0.insert(theTmpTrajectories0.end(),
theTmpTrajectories.begin(),
theTmpTrajectories.end());
#ifdef DEBUG
std::cout<<"Second path: We found trajectories for at least one seed in barrel layer "<<
theTmpTrajectories0.size()<<std::endl;
#endif
break;
} // There are trajectories
} // seeds
if( theTmpTrajectories0.size() > 0 ) {
#ifdef DEBUG
std::cout<<
"Second path: We found trajectories for at least one seed in barrel/endcap layer"<<
theTmpTrajectories0.size()<<std::endl;
#endif
break;
}
} // navigation maps
} // barrel or endcap seed
if( theTmpTrajectories0.size() > 0 ) {
#ifdef DEBUG
std::cout<<
"Second path: We found trajectories for at least one seed in barrel/endcap layer "
<<
theTmpTrajectories0.size()<<std::endl;
#endif
theFoundFts.push_back(&(*ifts));
theMapFtsTraj[&(*ifts)] = theTmpTrajectories0;
break;
}
} // seedlayer
} // tsos. isvalid
} // Muon Free trajectory state
//
// start fitting procedure
//
#ifdef DEBUG
if(theFoundFts.size()>0 ) {
std::cout<<" Event reconstruction finished with "<<theFoundFts.size()
<<std::endl;}
else {std::cout<<" Event reconstruction::no tracks found"<<std::endl;}
#endif
if(theFoundFts.size()<2) return dimuon;
// Look for vertex constraints
edm::ESHandle<GlobalTrackingGeometry> globTkGeomHandle;
es1.get<GlobalTrackingGeometryRecord>().get(globTkGeomHandle);
reco::BeamSpot::CovarianceMatrix matrix;
matrix(2,2) = 0.001;
matrix(3,3) = 0.001;
reco::BeamSpot bs( reco::BeamSpot::Point(0., 0., theHICConst->zvert),
0.1,
0.,
0.,
0.,
matrix
);
reco::TrackBase::TrackAlgorithm Algo = reco::TrackBase::undefAlgorithm;
// For trajectory refitting
vector<reco::Track> firstTrack;
vector<reco::TransientTrack> firstTransTracks;
vector<reco::TrackRef> firstTrackRefs;
vector<reco::Track> secondTrack;
vector<reco::TransientTrack> secondTransTracks;
vector<reco::TrackRef> secondTrackRefs;
for(vector<FreeTrajectoryState*>::iterator it = theFoundFts.begin(); it!= theFoundFts.end(); it++)
{
vector<Trajectory> first = (*theMapFtsTraj.find(*it)).second;
for(vector<Trajectory>::iterator im=first.begin();im!=first.end(); im++) {
TrajectoryStateOnSurface innertsos;
if (im->direction() == alongMomentum) {
innertsos = im->firstMeasurement().updatedState();
} else {
innertsos = im->lastMeasurement().updatedState();
}
// CMSSW31X
TSCBLBuilderNoMaterial tscblBuilder;
// CMSSW22X
//TrajectoryStateClosestToBeamLineBuilder tscblBuilder;
TrajectoryStateClosestToBeamLine tscbl = tscblBuilder(*(innertsos.freeState()),bs);
if (tscbl.isValid()==false) {
//cout<<" false track "<<endl;
continue;
}
GlobalPoint v = tscbl.trackStateAtPCA().position();
math::XYZPoint pos( v.x(), v.y(), v.z() );
// cout<<" Position of track close to vertex "<<v.perp()<<" "<<v.z()<<" Primary vertex "<<theHICConst->zvert<<
// " charge "<<tscbl.trackStateAtPCA().charge()<<endl;
if(v.perp() > 0.1 ) continue;
if(fabs(v.z() - theHICConst->zvert ) > 0.4 ) continue;
GlobalVector p = tscbl.trackStateAtPCA().momentum();
math::XYZVector mom( p.x(), p.y(), p.z() );
Track theTrack(im->chiSquared(),
im->ndof(),
pos, mom, tscbl.trackStateAtPCA().charge(),
tscbl.trackStateAtPCA().curvilinearError(),Algo);
TransientTrack tmpTk( theTrack, &(*magfield), globTkGeomHandle );
firstTrack.push_back( theTrack );
firstTransTracks.push_back( tmpTk );
}
if( firstTrack.size() == 0 ) continue;
for(vector<FreeTrajectoryState*>::iterator jt = it+1; jt!= theFoundFts.end(); jt++)
{
vector<Trajectory> second = (*theMapFtsTraj.find(*jt)).second;
// cout<<" Number of trajectories first "<<first.size()<<" second "<<second.size()<<endl;
for(vector<Trajectory>::iterator im=second.begin();im!=second.end(); im++) {
TrajectoryStateOnSurface innertsos;
if (im->direction() == alongMomentum) {
innertsos = im->firstMeasurement().updatedState();
} else {
innertsos = im->lastMeasurement().updatedState();
}
TSCBLBuilderNoMaterial tscblBuilder;
//TrajectoryStateClosestToBeamLineBuilder tscblBuilder;
TrajectoryStateClosestToBeamLine tscbl = tscblBuilder(*(innertsos.freeState()),bs);
if (tscbl.isValid()==false) {
//cout<<" false track "<<endl;
continue;
}
GlobalPoint v = tscbl.trackStateAtPCA().position();
math::XYZPoint pos( v.x(), v.y(), v.z() );
// cout<<" Position of track close to vertex "<<v.perp()<<" "<<v.z()<<" Primary vertex "<<theHICConst->zvert<<
// " charge "<<tscbl.trackStateAtPCA().charge()<<endl;
if(v.perp() > 0.1 ) continue;
if(fabs(v.z() - theHICConst->zvert ) > 0.4 ) continue;
GlobalVector p = tscbl.trackStateAtPCA().momentum();
math::XYZVector mom( p.x(), p.y(), p.z() );
Track theTrack(im->chiSquared(),
im->ndof(),
pos, mom, tscbl.trackStateAtPCA().charge(),
tscbl.trackStateAtPCA().curvilinearError(),Algo);
TransientTrack tmpTk( theTrack, &(*magfield), globTkGeomHandle );
secondTrack.push_back( theTrack );
secondTransTracks.push_back( tmpTk );
}
if( secondTrack.size() == 0 ) continue;
} // FTS
} // FTS
if( firstTrack.size() < 1 || secondTrack.size() < 1 ){
#ifdef DEBUG
cout<<" No enough tracks to get vertex "<<endl;
#endif
return dimuon;
}
bool useRefTrax=true;
KalmanVertexFitter theFitter(useRefTrax);
TransientVertex theRecoVertex;
//
// Create possible two particles vertices
//
vector<reco::TransientTrack> theTwoTransTracks;
vector<TransientVertex> theVertexContainer;
for(vector<reco::TransientTrack>::iterator iplus = firstTransTracks.begin();
iplus != firstTransTracks.end(); iplus++)
{
for(vector<reco::TransientTrack>::iterator iminus = secondTransTracks.begin();
iminus != secondTransTracks.end(); iminus++)
{
// To chech CAIR error before the vertex fitting
TwoTrackMinimumDistance ttmd;
bool CAIR_ST = false;
GlobalTrajectoryParameters sta = (*iminus).impactPointState().globalParameters();
GlobalTrajectoryParameters stb = (*iplus).impactPointState().globalParameters();
ClosestApproachInRPhi theIniAlgo;
CAIR_ST = theIniAlgo.calculate( sta, stb );
//cout<<"%%%%% CAIR_ST : "<<CAIR_ST<<" %%%%%"<<endl;
if(CAIR_ST == 0) continue;
theTwoTransTracks.clear();
theTwoTransTracks.push_back(*iplus);
theTwoTransTracks.push_back(*iminus);
theRecoVertex = theFitter.vertex(theTwoTransTracks);
if( !theRecoVertex.isValid() ) {
continue;
}
// cout<<" Vertex is found "<<endl;
// cout<<" Chi2 = "<<theRecoVertex.totalChiSquared()<<
// " r= "<<theRecoVertex.position().perp()<<
// " z= "<<theRecoVertex.position().z()<<endl;
// Additional cuts
if ( theRecoVertex.totalChiSquared() > 0.0002 ) {
// cout<<" Vertex is failed with Chi2 : "<<theRecoVertex.totalChiSquared()<<endl;
continue;
}
if ( theRecoVertex.position().perp() > 0.08 ) {
// cout<<" Vertex is failed with r position : "<<theRecoVertex.position().perp()<<endl;
continue;
}
if ( fabs(theRecoVertex.position().z()-theHICConst->zvert) > 0.06 ) {
// cout<<" Vertex is failed with z position : "<<theRecoVertex.position().z()<<endl;
continue;
}
double quality = theRecoVertex.normalisedChiSquared();
std::vector<reco::TransientTrack> tracks = theRecoVertex.originalTracks();
for (std::vector<reco::TransientTrack>::iterator ivb = tracks.begin(); ivb != tracks.end(); ivb++)
{
quality = quality * (*ivb).chi2() /(*ivb).ndof();
}
if( quality > 70. ) {
// cout<<" Vertex failed quality cut "<<quality<<endl;
continue;
}
theVertexContainer.push_back(theRecoVertex);
dimuon = true;
break;
} // iminus
if(dimuon) break;
} // iplus
return dimuon;
}
std::string cms::HITrackVertexMaker::builderName [private] |
Definition at line 86 of file HITrackVertexMaker.h.
int cms::HITrackVertexMaker::eventCount [private] |
Definition at line 80 of file HITrackVertexMaker.h.
Definition at line 81 of file HITrackVertexMaker.h.
Definition at line 83 of file HITrackVertexMaker.h.
Definition at line 85 of file HITrackVertexMaker.h.
Definition at line 82 of file HITrackVertexMaker.h.
std::vector<const NavigationSchool*> cms::HITrackVertexMaker::theNavigationSchoolV [private] |
Definition at line 87 of file HITrackVertexMaker.h.
1.7.3