Inheritance diagram for vertex::AggregateVertexAna:

Public Member Functions
	AggregateVertexAna (fhicl::ParameterSet const &pset)

	~AggregateVertexAna ()

Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDAnalyzer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept

	Analyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	Analyzer (Table< Config > const &config)

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Observer
	~Observer () noexcept

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

fhicl::ParameterSetID	selectorConfig () const

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Member Functions
void	analyze (const art::Event &evt)

void	beginJob ()

Private Attributes
TH1F *	HnTrksVtx

TH1F *	HnVtxes

TH1F *	HVtxSep

TH2F *	HVtxRZ

std::string	fDBScanModuleLabel

std::string	fHoughModuleLabel

std::string	fHitModuleLabel

std::string	fTrack3DModuleLabel

std::string	fEndPointModuleLabel

std::string	fVertexModuleLabel

art::PtrVector< recob::Hit >	fhitlist

art::PtrVector< recob::EndPoint2D >	feplist

art::PtrVector< recob::Track >	ftracklist

art::PtrVector< recob::Vertex >	fVertexlist

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

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

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 34 of file AggregateVertexAna_module.cc.

Constructor & Destructor Documentation

vertex::AggregateVertexAna::AggregateVertexAna ( fhicl::ParameterSet const & pset )

explicit

Definition at line 71 of file AggregateVertexAna_module.cc.

     : EDAnalyzer(pset)
     , fHitModuleLabel     (pset.get< std::string >("FFFTHitModuleLabel"))
     , fTrack3DModuleLabel (pset.get< std::string >("Track3DModuleLabel"))
     , fEndPointModuleLabel(pset.get< std::string >("EndPointModuleLabel"))
     , fVertexModuleLabel  (pset.get< std::string >("VertexModuleLabel"))
   {
 
 
   }

vertex::AggregateVertexAna::~AggregateVertexAna ( )

Definition at line 83 of file AggregateVertexAna_module.cc.

84 {

85 }

Member Function Documentation

void vertex::AggregateVertexAna::analyze ( const art::Event & evt )

private

Definition at line 101 of file AggregateVertexAna_module.cc.

   {
     art::Handle< std::vector<recob::Hit> > hitListHandle;
     evt.getByLabel(fHitModuleLabel,hitListHandle);
     for(unsigned int ii = 0; ii < hitListHandle->size(); ++ii){
       art::Ptr<recob::Hit> hit(hitListHandle, ii);
       fhitlist.push_back(hit); // class member
     }
 
     art::Handle< std::vector<recob::EndPoint2D> > epListHandle;
     evt.getByLabel(fEndPointModuleLabel,epListHandle);
     for(unsigned int ii = 0; ii < epListHandle->size(); ++ii){
       art::Ptr<recob::EndPoint2D> ep(epListHandle, ii);
       feplist.push_back(ep); // class member
     }
 
     art::Handle< std::vector<recob::Track> > trackListHandle;
     evt.getByLabel(fTrack3DModuleLabel,trackListHandle);
     for(unsigned int ii = 0; ii < trackListHandle->size(); ++ii){
       art::Ptr<recob::Track> track(trackListHandle, ii);
       ftracklist.push_back(track); // class member
     }
 
     art::Handle< std::vector<recob::Vertex> > vertexListHandle;
     evt.getByLabel(fVertexModuleLabel,vertexListHandle);
     for(unsigned int ii = 0; ii < vertexListHandle->size(); ++ii){
       art::Ptr<recob::Vertex> vertex(vertexListHandle, ii);
       fVertexlist.push_back(vertex); // class member
     }
 
     HnVtxes->Fill(feplist.size(),1);
 
     art::FindManyP<recob::Track> fmt(vertexListHandle, evt, fVertexModuleLabel);
     art::FindManyP<recob::Hit>   fmh(vertexListHandle, evt, fVertexModuleLabel);
 
     for(size_t v1 = 0; v1 < fVertexlist.size(); ++v1)  {
 
       std::vector< art::Ptr<recob::Track> > tvlist = fmt.at(v1);
 
       HnTrksVtx->Fill(tvlist.size(),1);
 
       if(tvlist.size() < 1) continue;
 
       std::vector< art::Ptr<recob::Hit> > hitvlist = fmh.at(v1);
 
       // Hits no longer has XYZ() method. To get 3d hit position info I'm going to have to
       // loop on all the SpacePoints and loop on all Hits from there till it matches
       // one from this vertex. This affects the two Fill() efforts below. EC, 19-Nov-2010.
       art::PtrVector<recob::Hit>::const_iterator hitv = hitvlist.begin();
 
       for(size_t v2 = v1+1; v2 < fVertexlist.size(); ++v2){
 
         std::vector< art::Ptr<recob::Hit> > hitvlist2 = fmh.at(v2);
 
         std::vector< art::Ptr<recob::Hit> >::const_iterator hitv2 = hitvlist2.begin();
 
         // These two whiles should be each precisely one iteration long.
         while( hitv != hitvlist.end() ){
           while( hitv2 != hitvlist2.end() ){
             TVector3 dist;
             mf::LogInfo("AggregateVertexAna") << "AggregateVertexAna: dist is " << dist.Mag() << ".";
             HVtxSep->Fill(dist.Mag(),1);
             hitv2++;
           }
           hitv++;
         }// end loop over hitv entries
 
       }// end loop over v2
     }// end loop over v1
 
     return;
   }// end analyze

void vertex::AggregateVertexAna::beginJob ( )

privatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 88 of file AggregateVertexAna_module.cc.

   {
     art::ServiceHandle<art::TFileService const> tfs;
 
     HnVtxes   = tfs->make<TH1F>("Num Vertices","Num Vertices",8,-0.5,7.5);
     HVtxSep   = tfs->make<TH1F>("Vertices spacing","Vertices spacing",20,0.001,5.0);
     HVtxRZ    = tfs->make<TH2F>("Vtx in RZ","Vtx in RZ",20,-50.0,+50.0,20,0.0,50.0);
     HnTrksVtx = tfs->make<TH1F>("Tracks per vtx","Tracks per vtx",8,-0.5,7.5);
 
     return;
   }