---

KVFTest Class Reference

This is a very simple test analyzer mean to test the KalmanVertexFitter. More...

#include <RecoVertex/KVFTest/src/KVFTest.cc>

Inheritance diagram for KVFTest:

List of all members.

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob (edm::EventSetup const &)
virtual void	endJob ()
	KVFTest (const edm::ParameterSet &)
	~KVFTest ()
Private Member Functions
TrackingVertex	getSimVertex (const edm::Event &iEvent) const
Private Attributes
TrackAssociatorByChi2 *	associatorForParamAtPca
edm::ParameterSet	kvfPSet
std::string	outputFile_
TFile *	rootFile_
edm::ParameterSet	theConfig
std::string	trackLabel_
SimpleVertexTree *	tree

---

Detailed Description

This is a very simple test analyzer mean to test the KalmanVertexFitter.

Description: steers tracker primary vertex reconstruction and storage.

Implementation: <Notes on="" implementation>="">

Definition at line 42 of file KVFTest.h.

---

Constructor & Destructor Documentation

KVFTest::KVFTest

(

const edm::ParameterSet &

iConfig

)

[explicit]

Definition at line 27 of file KVFTest.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), kvfPSet, outputFile_, rootFile_, and trackLabel_.

  : theConfig(iConfig)
{
  trackLabel_ = iConfig.getParameter<std::string>("TrackLabel");
  outputFile_ = iConfig.getUntrackedParameter<std::string>("outputFile");
  kvfPSet = iConfig.getParameter<edm::ParameterSet>("KVFParameters");
  rootFile_ = TFile::Open(outputFile_.c_str(),"RECREATE"); 
  edm::LogInfo("RecoVertex/KVFTest") 
    << "Initializing KVF TEST analyser  - Output file: " << outputFile_ <<"\n";
}

KVFTest::~KVFTest

(

)

Definition at line 39 of file KVFTest.cc.

References rootFile_.

                  {
  delete rootFile_;
}

---

Member Function Documentation

void KVFTest::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)			[virtual]

Implements edm::EDAnalyzer.

Definition at line 61 of file KVFTest.cc.

References TrackAssociatorByChi2::associateRecoToSim(), associatorForParamAtPca, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), SimpleVertexTree::fill(), edm::EventSetup::get(), edm::Event::getByLabel(), getSimVertex(), edm::Event::id(), edm::Handle< T >::isValid(), TransientVertex::position(), edm::Handle< T >::product(), sv, trackLabel_, tree, and KalmanVertexFitter::vertex().

{
  edm::LogInfo("RecoVertex/KVFTest") 
    << "Reconstructing event number: " << iEvent.id() << "\n";
    
  // get RECO tracks from the event
  // `tks` can be used as a ptr to a reco::TrackCollection
  edm::Handle<edm::View<reco::Track> > tks;
  iEvent.getByLabel(trackLabel_, tks);
  if (!tks.isValid()) {
    edm::LogInfo("RecoVertex/KVFTest") 
      << "Exception during event number: " << iEvent.id()
      << "\n";
  } else {
    edm::LogInfo("RecoVertex/KVFTest") 
      << "Found: " << (*tks).size() << " reconstructed tracks" << "\n";
    cout << "got " << (*tks).size() << " tracks " << endl;
    
    // Transform Track to TransientTrack

    //get the builder:
    edm::ESHandle<TransientTrackBuilder> theB;
    iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",theB);
    //do the conversion:
    vector<TransientTrack> t_tks = (*theB).build(tks);

    edm::LogInfo("RecoVertex/KVFTest") 
      << "Found: " << t_tks.size() << " reconstructed tracks" << "\n";
    
    // Call the KalmanVertexFitter if more than 1 track
    if (t_tks.size() > 1) {
      //      KalmanVertexFitter kvf(kvfPSet);
      KalmanVertexFitter kvf(true);
      TransientVertex tv = kvf.vertex(t_tks);

      std::cout << "Position: " << Vertex::Point(tv.position()) << "\n";

      // For the analysis: compare to your SimVertex
      TrackingVertex sv = getSimVertex(iEvent);
      edm::Handle<TrackingParticleCollection>  TPCollectionH ;
      iEvent.getByLabel("trackingtruth","TrackTruth",TPCollectionH);
      if (!TPCollectionH.isValid()) {
        edm::LogInfo("RecoVertex/KVFTest") 
          << "Exception during event number: " << iEvent.id() 
          << "\n";
      } else {
        const TrackingParticleCollection tPC = *(TPCollectionH.product());
        reco::RecoToSimCollection recSimColl=associatorForParamAtPca->associateRecoToSim(tks,
                                                                                         TPCollectionH,
                                                                                         &iEvent);    
        tree->fill(tv, &sv, &recSimColl);
      }
    }
  }  
}

void KVFTest::beginJob

(

edm::EventSetup const &

setup

)

[virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 43 of file KVFTest.cc.

References associatorForParamAtPca, edm::EventSetup::get(), edm::ESHandle< T >::product(), and tree.

                                                {
  edm::ESHandle<TrackAssociatorBase> theAssociatorForParamAtPca;
  setup.get<TrackAssociatorRecord>().get("TrackAssociatorByChi2",theAssociatorForParamAtPca);
  associatorForParamAtPca = (TrackAssociatorByChi2 *) theAssociatorForParamAtPca.product();

  tree = new SimpleVertexTree("VertexFitter", associatorForParamAtPca);
}

void KVFTest::endJob

(

void

)

[virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 52 of file KVFTest.cc.

References tree.

                     {
  delete tree;
}

TrackingVertex KVFTest::getSimVertex

(

const edm::Event &

iEvent

)

const [private]

Definition at line 119 of file KVFTest.cc.

References edm::Event::getByLabel(), and edm::Handle< T >::product().

Referenced by analyze().

{
   // get the simulated vertices
  edm::Handle<TrackingVertexCollection>  TVCollectionH ;
  iEvent.getByLabel("trackingtruth","VertexTruth",TVCollectionH);
  const TrackingVertexCollection tPC = *(TVCollectionH.product());

//    Handle<edm::SimVertexContainer> simVtcs;
//    iEvent.getByLabel("g4SimHits", simVtcs);
//    std::cout << "SimVertex " << simVtcs->size() << std::endl;
//    for(edm::SimVertexContainer::const_iterator v=simVtcs->begin();
//        v!=simVtcs->end(); ++v){
//      std::cout << "simvtx "
//             << v->position().x() << " "
//             << v->position().y() << " "
//             << v->position().z() << " "
//             << v->parentIndex() << " "
//             << v->noParent() << " "
//               << std::endl;
//    }
   return *(tPC.begin());
}

---

Member Data Documentation

TrackAssociatorByChi2* KVFTest::associatorForParamAtPca [private]

Definition at line 58 of file KVFTest.h.

Referenced by analyze(), and beginJob().

edm::ParameterSet KVFTest::kvfPSet [private]

Definition at line 57 of file KVFTest.h.

Referenced by KVFTest().

std::string KVFTest::outputFile_ [private]

Definition at line 62 of file KVFTest.h.

Referenced by KVFTest().

TFile* KVFTest::rootFile_ [private]

Definition at line 60 of file KVFTest.h.

Referenced by KVFTest(), and ~KVFTest().

edm::ParameterSet KVFTest::theConfig [private]

Definition at line 56 of file KVFTest.h.

std::string KVFTest::trackLabel_ [private]

Definition at line 63 of file KVFTest.h.

Referenced by analyze(), and KVFTest().

SimpleVertexTree* KVFTest::tree [private]

Definition at line 59 of file KVFTest.h.

Referenced by analyze(), beginJob(), and endJob().

---

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

• RecoVertex/KalmanVertexFit/plugins/KVFTest.h
• RecoVertex/KalmanVertexFit/plugins/KVFTest.cc

---