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

Inheritance diagram for KVFTest:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

virtual void	endJob ()

	KVFTest (const edm::ParameterSet &)

	~KVFTest ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

edm::EDGetTokenT < reco::TrackCollection >	token_tracks

edm::EDGetTokenT < TrackingParticleCollection >	token_TrackTruth

edm::EDGetTokenT < TrackingVertexCollection >	token_VertexTruth

SimpleVertexTree *	tree

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 37 of file KVFTest.h.

Constructor & Destructor Documentation

KVFTest::KVFTest ( const edm::ParameterSet & iConfig )

explicit

Definition at line 24 of file KVFTest.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), kvfPSet, outputFile_, rootFile_, AlCaHLTBitMon_QueryRunRegistry::string, token_tracks, token_TrackTruth, and token_VertexTruth.

   : theConfig(iConfig), associatorForParamAtPca(0), tree(0)
 {
   token_tracks = consumes<TrackCollection>(iConfig.getParameter<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";
 
   token_TrackTruth = consumes<TrackingParticleCollection>(edm::InputTag("trackingtruth", "TrackTruth"));
   token_VertexTruth = consumes<TrackingVertexCollection>(edm::InputTag("trackingtruth", "VertexTruth"));
 
 }

KVFTest::~KVFTest ( )

Definition at line 40 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 57 of file KVFTest.cc.

References TrackAssociatorByChi2::associateRecoToSim(), associatorForParamAtPca, gather_cfg::cout, SimpleVertexTree::fill(), edm::EventSetup::get(), edm::Event::getByToken(), getSimVertex(), edm::EventBase::id(), edm::HandleBase::isValid(), TransientVertex::position(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), token_tracks, token_TrackTruth, tree, and KalmanVertexFitter::vertex().

 {
   if ( associatorForParamAtPca==0 ) {
     edm::ESHandle<TrackAssociatorBase> theAssociatorForParamAtPca;
     iSetup.get<TrackAssociatorRecord>().get("TrackAssociatorByChi2",theAssociatorForParamAtPca);
     associatorForParamAtPca = (TrackAssociatorByChi2 *) theAssociatorForParamAtPca.product();
 
     tree = new SimpleVertexTree("VertexFitter", associatorForParamAtPca);
   }
 
 
 
   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.getByToken(token_tracks, 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";
     std::cout << "got " << (*tks).size() << " tracks " << std::endl;
     
     // Transform Track to TransientTrack
 
     //get the builder:
     edm::ESHandle<TransientTrackBuilder> theB;
     iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",theB);
     //do the conversion:
     std::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.getByToken(token_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 ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 44 of file KVFTest.cc.

44 {

45 }

void KVFTest::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 48 of file KVFTest.cc.

References tree.

                      {
   delete tree;
 }

TrackingVertex KVFTest::getSimVertex ( const edm::Event & iEvent ) const

private

Definition at line 125 of file KVFTest.cc.

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

Referenced by analyze().

 {
    // get the simulated vertices
   edm::Handle<TrackingVertexCollection>  TVCollectionH ;
   iEvent.getByToken(token_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());
 }