---

cms::HITrackVertexMaker Class Reference

#include <RecoHIMuon/HiMuTracking/interface/HITrackVertexMaker.h>

List of all members.

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
edm::InputTag	L2candTag_
edm::ESHandle< MagneticField >	magfield
edm::ESHandle< MeasurementTracker >	measurementTrackerHandle
edm::InputTag	primaryVertexTag
edm::ESHandle< Propagator >	propagatorAlongHandle
edm::ESHandle< Propagator >	propagatorOppositeHandle
edm::ParameterSet	pset_
edm::ESHandle < TransientTrackingRecHitBuilder >	recHitBuilderHandle
edm::InputTag	rphirecHitsTag
HICMeasurementEstimator *	theEstimator
MinPtTrajectoryFilter *	theMinPtFilter
const NavigationSchool *	theNavigationSchool
HICTrajectoryBuilder *	theTrajectoryBuilder
edm::ESHandle < GeometricSearchTracker >	tracker
edm::ESHandle < TrajectoryStateUpdator >	updatorHandle

---

Detailed Description

Definition at line 62 of file HITrackVertexMaker.h.

---

Constructor & Destructor Documentation

cms::HITrackVertexMaker::HITrackVertexMaker	(	const edm::ParameterSet &	ps1,
		const edm::EventSetup &	es1
	)			[explicit]

Definition at line 104 of file HITrackVertexMaker.cc.

References builderName, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), edm::EventSetup::get(), edm::ParameterSet::getParameter(), L2candTag_, magfield, measurementTrackerHandle, primaryVertexTag, edm::ESHandle< T >::product(), propagatorAlongHandle, propagatorOppositeHandle, pset_, ptmin, recHitBuilderHandle, rphirecHitsTag, theEstimator, theMinPtFilter, theTrajectoryBuilder, tracker, and updatorHandle.

{

   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;

//
// Initializetion from Records
    es1.get<TrackerRecoGeometryRecord>().get( tracker );
    es1.get<IdealMagneticFieldRecord>().get(magfield);
    es1.get<CkfComponentsRecord>().get("",measurementTrackerHandle);
    es1.get<TransientRecHitRecord>().get(builderName,recHitBuilderHandle); 
       
    double theChiSquareCut = 500.;
    double nsig = 3.;
    int theLowMult = 1;
    theEstimator = new HICMeasurementEstimator(&(*tracker), &(*magfield), theChiSquareCut, nsig);
    std::string updatorName = "KFUpdator";
    std::string propagatorAlongName    = "PropagatorWithMaterial";
    std::string propagatorOppositeName = "PropagatorWithMaterialOpposite";
    es1.get<TrackingComponentsRecord>().get(propagatorAlongName,propagatorAlongHandle);
    es1.get<TrackingComponentsRecord>().get(propagatorOppositeName,propagatorOppositeHandle);
    es1.get<TrackingComponentsRecord>().get(updatorName,updatorHandle);
    double ptmin=1.;
    theMinPtFilter = new MinPtTrajectoryFilter(ptmin); 
    theTrajectoryBuilder = new HICTrajectoryBuilder(        pset_, es1,
                                                     updatorHandle.product(),
                                                     propagatorAlongHandle.product(),
                                                     propagatorOppositeHandle.product(),
                                                     theEstimator,
                                                     recHitBuilderHandle.product(),
                                                     measurementTrackerHandle.product(),
                                                     theMinPtFilter);
#ifdef DEBUG
    std::cout<<" HICTrajectoryBuilder constructed "<<std::endl;
#endif
}

cms::HITrackVertexMaker::~HITrackVertexMaker

(

)

Definition at line 153 of file HITrackVertexMaker.cc.

{
} 

---

Member Function Documentation

bool cms::HITrackVertexMaker::produceTracks	(	const edm::Event &	e1,
		const edm::EventSetup &	es1,
		HICConst *	theHICConst,
		FmpConst *	theFmpConst
	)

Definition at line 157 of file HITrackVertexMaker.cc.

References alongMomentum, builderName, FreeTrajectoryState::charge(), GenMuonPlsPt100GeV_cfg::cout, cms::HICFTSfromL1orL2::createFTSfromL2(), FreeTrajectoryState::curvilinearError(), lat::endl(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::freeTrajectoryState(), edm::EventSetup::get(), edm::Event::getByLabel(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), TransientVertex::hasRefittedTracks(), insideOut, TrajectoryStateOnSurface::isValid(), TrajectoryStateClosestToBeamLine::isValid(), TransientVertex::isValid(), it, L2candTag_, magfield, matrix, measurementTrackerHandle, FreeTrajectoryState::momentum(), VarParsing::mult, TransientVertex::normalisedChiSquared(), oppositeToMomentum, TransientVertex::originalTracks(), outsideIn, p, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), TransientVertex::position(), FreeTrajectoryState::position(), primaryVertexTag, cms::DiMuonSeedGeneratorHIC::produce(), edm::ESHandle< T >::product(), recHitBuilderHandle, TransientVertex::refittedTracks(), rphirecHitsTag, edm::second(), HICMeasurementEstimator::setHICConst(), HICMeasurementEstimator::setMult(), HICTrajectoryBuilder::settracker(), cms::FmpConst::setVertex(), cms::HICConst::setVertex(), cms::HICTkOuterStartingLayerFinder::startingLayers(), state, theEstimator, theNavigationSchool, theTrajectoryBuilder, TransientVertex::totalChiSquared(), tracker, tracks, TrajectoryStateClosestToBeamLine::trackStateAtPCA(), HICTrajectoryBuilder::trajectories(), 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(), and cms::TestMuL1L2FilterSTA::filter().

{
    
   bool dimuon = false;
   edm::Handle<reco::VertexCollection> vertexcands;
   e1.getByLabel (primaryVertexTag,vertexcands);
   int iv = 0;

//   cout<<" Number of vertices primary  "<<vertexcands->size()<<endl;

   if(vertexcands->size()<1) return dimuon;

   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 "<<theHICConst->zvert<<endl;
 
//
// Get recHit builder
//

  es1.get<TransientRecHitRecord>().get(builderName,recHitBuilderHandle);

//
// Get measurement tracker
//
//  
 
  es1.get<CkfComponentsRecord>().get("",measurementTrackerHandle);
    
#ifdef DEBUG  
  std::cout<<" Before first tracker update "<<std::endl;
#endif

  measurementTrackerHandle->update(e1);
  
#ifdef DEBUG
  std::cout<<" After first tracker update "<<std::endl;
#endif
//
// Get L1 muon info
//

  vector<L1MuGMTExtendedCand> excall;

/*  
try{  
edm::Handle<L1MuGMTReadoutCollection> gmtrc1_handle;
e1.getByLabel("gmtDigis",gmtrc1_handle);

L1MuGMTReadoutCollection const* gmtrc1 = gmtrc1_handle.product();


vector<L1MuGMTReadoutRecord> gmtrecords = gmtrc1->getRecords();
vector<L1MuGMTReadoutRecord>::const_iterator igmtrr1;
std::cout<<" Number of :L1 "<<gmtrecords.size()<<std::endl;


for(igmtrr1=gmtrecords.begin(); igmtrr1!=gmtrecords.end(); igmtrr1++)
{

vector<L1MuRegionalCand>::const_iterator iter11;
vector<L1MuRegionalCand> rmc1;

vector<L1MuGMTExtendedCand>::const_iterator gmt_iter1; 
vector<L1MuGMTExtendedCand> exc1 = igmtrr1->getGMTCands();

std::cout<<" Number of :L1MuGMTExtendedCand "<<exc1.size()<<std::endl;
int igmt1 = 0;
for(gmt_iter1 = exc1.begin(); gmt_iter1!=exc1.end(); gmt_iter1++)
        {
        double etal1=(*gmt_iter1).etaValue();
        double phil1=(*gmt_iter1).phiValue();
        int chargel1=(*gmt_iter1).charge();
        double ptl1=(*gmt_iter1).ptValue();
        cout<<"etal1= "<<etal1<<"phil1= "<<phil1<<"ptl1= "<<ptl1<<"chargel1= "<<chargel1<<endl;
        igmt1++;
        excall.push_back(*gmt_iter1);

        } 
}

} 
        catch (cms::Exception& e) { // can't find it!
              throw e;
       }

*/

   edm::Handle<RecoChargedCandidateCollection> L2mucands;
   //edm::Handle<TrackCollection> L2mucands;
   e1.getByLabel (L2candTag_,L2mucands);
   RecoChargedCandidateCollection::const_iterator L2cand1;
   RecoChargedCandidateCollection::const_iterator L2cand2;
      
#ifdef DEBUG
   cout<<" Number of muon candidates "<<L2mucands->size()<<endl;
#endif

   if( L2mucands->size() < 2 ) return dimuon;
   
// For trajectory builder   
   int  theLowMult = 1;
   theEstimator->setHICConst(theHICConst);
   theEstimator->setMult(theLowMult);
   theNavigationSchool = new HICSimpleNavigationSchool(&(*tracker), &(*magfield));
   theTrajectoryBuilder->settracker(measurementTrackerHandle.product());
//============================   
    FastMuPropagator* theFmp = new FastMuPropagator(&(*magfield),theFmpConst); 
    StateOnTrackerBound state(theFmp); 
    TrajectoryStateOnSurface tsos;
    
#ifdef DEBUG   
//   for (cand1=L2mucands->begin(); cand1!=L2mucands->end(); cand1++) {
//     reco::TrackRef mytr = (*cand1).track(); 
//      std::cout<<" Inner position "<<(*mytr).innerPosition().x()<<" "<<(*mytr).innerPosition().y()<<" "<<(*mytr).innerPosition().z()<<std::endl;
//
//   } 
#endif 
  
    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> allTraj;
   
   int numofseeds = 0; 
   for(vector<FreeTrajectoryState>::iterator ifts=theFts.begin(); ifts!=theFts.end(); ifts++)
   {
#ifdef DEBUG
     cout<<" cycle on Muon Trajectory State " <<(*ifts).parameters().position().perp()<<
                                          " " <<(*ifts).parameters().position().z()   <<endl;
#endif
     tsos=state((*ifts));
     
     int iseedn = 0; 
     
     if(tsos.isValid())
     {
     
#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 ) continue;

        map<DetLayer*, DiMuonSeedGeneratorHIC::SeedContainer> seedmap = Seed.produce(e1 ,es1, (*ftsnew), tsos, (*ifts), 
                                               recHitBuilderHandle.product(), measurementTrackerHandle.product(), &seedlayers);
        
        int ilnum = 0;
          for( vector<DetLayer*>::iterator it = seedlayers.begin(); it != seedlayers.end(); it++)
          {
            ilnum++;
            vector<Trajectory> theTmpTrajectories0;
            DiMuonSeedGeneratorHIC::SeedContainer seeds = (*seedmap.find(*it)).second;
    // set the correct navigation
            NavigationSetter setter( *theNavigationSchool);
#ifdef DEBUG
       std::cout<<" Layer "<<ilnum<<" Number of seeds in layer "<<seeds.size()<<std::endl;          
#endif      
            if(seeds.size() > 0) {
       for(DiMuonSeedGeneratorHIC::SeedContainer::iterator iseed=seeds.begin();iseed!=seeds.end();iseed++)
       {
         std::vector<TrajectoryMeasurement> theV = (*iseed).measurements();
#ifdef DEBUG
         std::cout<< "Layer " <<ilnum<<" Seed number "<<iseedn<<"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) {
                                            allTraj.insert(allTraj.end(),theTmpTrajectories.begin(),theTmpTrajectories.end());
                                            theTmpTrajectories0.insert(theTmpTrajectories0.end(),theTmpTrajectories.begin(),theTmpTrajectories.end());
                                        }
        if(theTmpTrajectories0.size() > 0) { 
#ifdef DEBUG
          std::cout<<"We found trajectories for at least one seed "<<std::endl; 
#endif
          break;
        }
                                        
        iseedn++;
            
       } // seeds           
            } // there are seeds/layer
            
//          std::cout<<" Size of seed container "<<seeds.size()<<" size of trajectories "<<theTmpTrajectories0.size()<<std::endl;  
            
            if( theTmpTrajectories0.size() > 0 ) break;
            
          } // Detlayer cycle  
        
     } // tsos. isvalid
     if(iseedn > 0) numofseeds++;
   } // Muon Free trajectory state

#ifdef DEBUG
   cout<<" Number of muons trajectories in MUON "<<theFts.size()<<" Number of seeds "<<numofseeds+1<<endl;
#endif
        
//
// start fitting procedure
//
#ifdef DEBUG
    if(allTraj.size()>0 ) {std::cout<<" Event reconstruction finished with "<<allTraj.size()<<std::endl;} else {std::cout<<" Event reconstruction::no tracks found"<<std::endl;}
#endif
    if(allTraj.size()<2)  return dimuon;

    edm::ESHandle<GlobalTrackingGeometry> globTkGeomHandle;
    es1.get<GlobalTrackingGeometryRecord>().get(globTkGeomHandle);

    vector<reco::Track> thePositiveTracks;
    vector<reco::Track> theNegativeTracks; 
    vector<reco::TrackRef> thePositiveTrackRefs;
    vector<reco::TrackRef> theNegativeTrackRefs;
    vector<reco::TransientTrack> thePositiveTransTracks;
    vector<reco::TransientTrack> theNegativeTransTracks;

    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;

    edm::ESHandle<TrajectoryFitter> theFitterTrack;
    edm::ESHandle<TrajectorySmoother> theSmootherTrack;
    edm::ESHandle<Propagator> thePropagatorTrack;
    es1.get<TrackingComponentsRecord>().get("KFFitterForRefitInsideOut",theFitterTrack);
    es1.get<TrackingComponentsRecord>().get("KFSmootherForRefitInsideOut",theSmootherTrack);
    es1.get<TrackingComponentsRecord>().get("SmartPropagatorAny",thePropagatorTrack);
    enum RefitDirection{insideOut,outsideIn,undetermined};


    for(vector<Trajectory>::iterator it = allTraj.begin(); it!= allTraj.end(); it++)
    {
//
// Refit the trajectory
//
// Refit the trajectory before putting into the last result

  Trajectory::ConstRecHitContainer recHitsForReFit = (*it).recHits();
  
  PTrajectoryStateOnDet garbage1;
  edm::OwnVector<TrackingRecHit> garbage2;
  TrajectoryStateOnSurface firstTSOS = (*it).firstMeasurement().updatedState();

//   cout<<" firstTSOS "<<firstTSOS.freeTrajectoryState()->position().perp()<<" "<<firstTSOS.freeTrajectoryState()->position().z()<<endl;
//  PropagationDirection propDir =
//    (firstTSOS.globalPosition().basicVector().dot(firstTSOS.globalMomentum().basicVector())>0) ? alongMomentum : oppositeToMomentum;

  PropagationDirection propDir = oppositeToMomentum;
  
//  RefitDirection theRefitDirection = insideOut;
//  propDir = alongMomentum;
//  RefitDirection theRefitDirection = outsideIn;

//  if(propDir == alongMomentum && theRefitDirection == outsideIn)  {cout<<" oppositeToMomentum "<<endl; propDir=oppositeToMomentum;}
//  if(propDir == oppositeToMomentum && theRefitDirection == insideOut) {cout<<" alongMomentum "<<endl; propDir=alongMomentum;}

  TrajectorySeed seed(garbage1,garbage2,propDir);
  vector<Trajectory> trajectories = theFitterTrack->fit(seed,recHitsForReFit,firstTSOS);

//  vector<Trajectory> trajectories = theFitterTrack->fit(*it);

//  if(trajectories.empty()) {cout<<" No refit is done "<<endl;};
  
//  Trajectory trajectoryBW = trajectories.front();
//  vector<Trajectory> trajectoriesSM = theSmoother->trajectories(trajectoryBW);
//  if(trajectoriesSM.size()<1) continue;

    TrajectoryStateOnSurface innertsos;
    if (it->direction() == alongMomentum) {
   //    cout<<" Direction is along momentum "<<endl;
      innertsos = it->firstMeasurement().updatedState();
    } else {
      innertsos = it->lastMeasurement().updatedState();
   //   cout<<" Direction is opposite to momentum "<<endl;
      
    }
   //  cout<<" Position of the innermost point "<<innertsos.freeTrajectoryState()->position().perp()<<" "<<innertsos.freeTrajectoryState()->position().z()<<endl;
    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(it->chiSquared(),
                   it->ndof(), 
                   pos, mom, tscbl.trackStateAtPCA().charge(), tscbl.trackStateAtPCA().curvilinearError(),Algo);
    TransientTrack tmpTk( theTrack, &(*magfield), globTkGeomHandle );
    if( tscbl.trackStateAtPCA().charge() > 0)
    {
      thePositiveTracks.push_back( theTrack );
      thePositiveTransTracks.push_back( tmpTk );
    }
      else
      {
        theNegativeTracks.push_back( theTrack );
        theNegativeTransTracks.push_back( tmpTk );
      }
    } // Traj

   if( thePositiveTracks.size() < 1 || theNegativeTracks.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 = thePositiveTransTracks.begin(); iplus != thePositiveTransTracks.end(); iplus++)
   {
     theTwoTransTracks.clear();
     theTwoTransTracks.push_back(*iplus);
     for(vector<reco::TransientTrack>::iterator iminus = theNegativeTransTracks.begin(); iminus != theNegativeTransTracks.end(); iminus++)
     {
       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();
     vector<TransientTrack> refittedTrax;
     if( theRecoVertex.hasRefittedTracks() ) {
          refittedTrax = theRecoVertex.refittedTracks();
     }

     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);
     } // iminus
  } // iplus
  
       if( theVertexContainer.size() < 1 ) { 
#ifdef DEBUG
        std::cout<<" No vertex found in event "<<std::endl; 
#endif
       return dimuon; 
       } else {
          // cout<<" Event vertex is selected "<<endl;
       }


    dimuon = true;
    return dimuon;
} 

---

Member Data Documentation

std::string cms::HITrackVertexMaker::builderName [private]

Definition at line 87 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

edm::InputTag cms::HITrackVertexMaker::L2candTag_ [private]

Definition at line 82 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

edm::ESHandle<MagneticField> cms::HITrackVertexMaker::magfield [private]

Definition at line 89 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

edm::ESHandle<MeasurementTracker> cms::HITrackVertexMaker::measurementTrackerHandle [private]

Definition at line 91 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

edm::InputTag cms::HITrackVertexMaker::primaryVertexTag [private]

Definition at line 84 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

edm::ESHandle<Propagator> cms::HITrackVertexMaker::propagatorAlongHandle [private]

Definition at line 98 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker().

edm::ESHandle<Propagator> cms::HITrackVertexMaker::propagatorOppositeHandle [private]

Definition at line 99 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker().

edm::ParameterSet cms::HITrackVertexMaker::pset_ [private]

Definition at line 86 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker().

edm::ESHandle<TransientTrackingRecHitBuilder> cms::HITrackVertexMaker::recHitBuilderHandle [private]

Definition at line 90 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

edm::InputTag cms::HITrackVertexMaker::rphirecHitsTag [private]

Definition at line 83 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

HICMeasurementEstimator* cms::HITrackVertexMaker::theEstimator [private]

Definition at line 96 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

MinPtTrajectoryFilter* cms::HITrackVertexMaker::theMinPtFilter [private]

Definition at line 95 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker().

const NavigationSchool* cms::HITrackVertexMaker::theNavigationSchool [private]

Definition at line 93 of file HITrackVertexMaker.h.

Referenced by produceTracks().

HICTrajectoryBuilder* cms::HITrackVertexMaker::theTrajectoryBuilder [private]

Definition at line 94 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

edm::ESHandle<GeometricSearchTracker> cms::HITrackVertexMaker::tracker [private]

Definition at line 92 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker(), and produceTracks().

edm::ESHandle<TrajectoryStateUpdator> cms::HITrackVertexMaker::updatorHandle [private]

Definition at line 100 of file HITrackVertexMaker.h.

Referenced by HITrackVertexMaker().

---

The documentation for this class was generated from the following files:

• RecoHIMuon/HiMuTracking/interface/HITrackVertexMaker.h
• RecoHIMuon/HiMuTracking/src/HITrackVertexMaker.cc

---