Inheritance diagram for vertex::AggregateVertexAna:

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

	~AggregateVertexAna ()

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

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 ( )

private

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;
   }