Inheritance diagram for TrackHitAna::TrackHitAna:

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

virtual	~TrackHitAna ()

void	beginJob ()

void	endJob ()

void	beginRun (const art::Run &run)

void	reconfigure (fhicl::ParameterSet const &pset)

void	analyze (const art::Event &evt)

Private Member Functions
int	traversePFParticleHierarchy (art::Handle< std::vector< recob::PFParticle >> &, size_t, const art::FindManyP< recob::Track > &, const art::FindManyP< recob::Vertex > &, int &, int &) const

double	length (const recob::Track *track)

double	projectedLength (const recob::Track *track)

Private Attributes
std::vector< art::InputTag >	fHitProducerLabelVec

std::vector< art::InputTag >	fWireProducerLabelVec

std::vector< art::InputTag >	fPFParticleProducerLabelVec

std::vector< art::InputTag >	fTrackProducerLabelVec

int	fEvent

int	fRun

int	fSubRun

int	fNumEvents

std::vector< std::unique_ptr < IHitHistogramTool > >	fHitHistogramToolVec

std::vector< std::vector < double > >	fChannelPedVec

const geo::GeometryCore *	fGeometry

const lariov::DetPedestalProvider &	fPedestalRetrievalAlg
	Keep track of an instance to the pedestal retrieval alg. More...

Detailed Description

Definition at line 76 of file TrackHitAna_module.cc.

Constructor & Destructor Documentation

TrackHitAna::TrackHitAna::TrackHitAna ( fhicl::ParameterSet const & pset )

explicit

Definition at line 147 of file TrackHitAna_module.cc.

     : EDAnalyzer(parameterSet),
       fPedestalRetrievalAlg(*lar::providerFrom<lariov::DetPedestalService>())
 
 {
     fGeometry = lar::providerFrom<geo::Geometry>();
 
     // Read in the parameters from the .fcl file.
     this->reconfigure(parameterSet);
 }

TrackHitAna::TrackHitAna::~TrackHitAna ( )

virtual

Definition at line 160 of file TrackHitAna_module.cc.

161 {}

Member Function Documentation

void TrackHitAna::TrackHitAna::analyze ( const art::Event & evt )

Definition at line 213 of file TrackHitAna_module.cc.

 {
     // Start by fetching some basic event information for our n-tuple.
     fEvent  = event.id().event();
     fRun    = event.run();
     fSubRun = event.subRun();
 
     fNumEvents++;
     
     for(const auto& hitLabel : fHitProducerLabelVec)
     {
         // Make a pass through all hits to make contrasting plots
         art::Handle< std::vector<recob::Hit> > hitHandle;
         event.getByLabel(hitLabel, hitHandle);
         
         if (hitHandle.isValid())
         {
             IHitHistogramTool::HitPtrVec allHitVec;
             art::fill_ptr_vector(allHitVec, hitHandle);
             
             for(auto& hitHistTool : fHitHistogramToolVec) hitHistTool->fillHistograms(allHitVec);
         }
     }
 
     return;
 }

void TrackHitAna::TrackHitAna::beginJob ( )

Definition at line 164 of file TrackHitAna_module.cc.

 {
     // Get the detector length, to determine the maximum bin edge of one
     // of the histograms.
     //double detectorLength = fGeometry->DetLength();
 
     // Access ART's TFileService, which will handle creating and writing
     // histograms and n-tuples for us.
     art::ServiceHandle<art::TFileService> tfs;
 
     // The arguments to 'make<whatever>' are the same as those passed
     // to the 'whatever' constructor, provided 'whatever' is a ROOT
     // class that TFileService recognizes.
     for (auto& hitHistTool : fHitHistogramToolVec) hitHistTool->initializeHists(tfs, "AllHits");
 
     // zero out the event counter
     fNumEvents = 0;
 }

void TrackHitAna::TrackHitAna::beginRun ( const art::Run & run )

Definition at line 184 of file TrackHitAna_module.cc.

 {
     // How to convert from number of electrons to GeV.  The ultimate
     // source of this conversion factor is
     // ${LARSIM_DIR}/include/SimpleTypesAndConstants/PhysicalConstants.h.
 //    art::ServiceHandle<sim::LArG4Parameters> larParameters;
 //    fElectronsToGeV = 1./larParameters->GeVToElectrons();
 }

void TrackHitAna::TrackHitAna::endJob ( )

Definition at line 240 of file TrackHitAna_module.cc.

 {
     // Make a call to normalize histograms if so desired
     for(auto& hitHistTool : fHitHistogramToolVec) hitHistTool->endJob(fNumEvents);
 
     return;
 }

double TrackHitAna::TrackHitAna::length ( const recob::Track * track )

private

Definition at line 250 of file TrackHitAna_module.cc.

 {
     double   result(0.);
 /*
     TVector3 disp(track->LocationAtPoint(0));
     TVector3 lastPoint(track->LocationAtPoint(0));
     TVector3 lastDir(0.,0.,0.);
     int      n(track->NumberTrajectoryPoints());
 
     for(int i = 1; i < n; ++i)
     {
         const TVector3& pos = track->LocationAtPoint(i);
 
         TVector3 trajDir = pos - lastPoint;
 
         if (trajDir.Mag2()) trajDir.SetMag(1.);
 
 //        if (lastDir.Dot(trajDir) >= 0.)
 //        {
             disp   -= pos;
             result += disp.Mag();
             disp    = pos;
 //        }
 
         lastPoint = pos;
         lastDir   = trajDir;
     }
 */
     return result;
 }

double TrackHitAna::TrackHitAna::projectedLength ( const recob::Track * track )

private

Definition at line 283 of file TrackHitAna_module.cc.

 {
     double   result(0.);
 /*
     TVector3 lastPoint(track->LocationAtPoint(0));
     TVector3 lastDir(track->DirectionAtPoint(0));
     int      n(track->NumberTrajectoryPoints());
 
     for(int i = 1; i < n; ++i)
     {
         const TVector3& newPoint = track->LocationAtPoint(i);
 
         TVector3 lastToNewPoint = newPoint - lastPoint;
         double   arcLenToDoca   = lastDir.Dot(lastToNewPoint);
 
         result    += arcLenToDoca;
         lastPoint  = lastPoint + arcLenToDoca * lastDir;
         lastDir    = track->DirectionAtPoint(i);
     }
 */
     return result;
 }

void TrackHitAna::TrackHitAna::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 194 of file TrackHitAna_module.cc.

 {
     // Read parameters from the .fcl file. The names in the arguments
     // to p.get<TYPE> must match names in the .fcl file.
     fHitProducerLabelVec        = p.get< std::vector<art::InputTag> >("HitModuleLabel",          std::vector<art::InputTag>() = {"gauss"});
     fPFParticleProducerLabelVec = p.get< std::vector<art::InputTag> >("PFParticleProducerLabel", std::vector<art::InputTag>() = {"cluster3d"});
     fTrackProducerLabelVec      = p.get< std::vector<art::InputTag> >("TrackProducerLabel",      std::vector<art::InputTag>() = {"trackkalmanhit"});
     fWireProducerLabelVec       = p.get< std::vector<art::InputTag> >("WireProducerLabel",       std::vector<art::InputTag>() = {"caldata"});
     
     // Implement the tools for handling the responses
     const std::vector<fhicl::ParameterSet>& hitHistogramToolVec = p.get<std::vector<fhicl::ParameterSet>>("HitHistogramToolList");
     
     for(auto& hitHistogramTool : hitHistogramToolVec)
         fHitHistogramToolVec.push_back(art::make_tool<IHitHistogramTool>(hitHistogramTool));
 
     return;
 }

int TrackHitAna::TrackHitAna::traversePFParticleHierarchy	(	art::Handle< std::vector< recob::PFParticle >> &	pfParticleHandle,
		size_t	pfParticleIdx,
		const art::FindManyP< recob::Track > &	trackAssns,
		const art::FindManyP< recob::Vertex > &	vertexAssns,
		int &	nTracks,
		int &	nVertices
	)		const

private

Definition at line 306 of file TrackHitAna_module.cc.

 {
     // So far no daughters...
     int nDaughters(0);
 /*
     // Get pointer to PFParticle
     art::Ptr<recob::PFParticle> pfParticle(pfParticleHandle,pfParticleIdx);
 
     // Recover tracks/vertices associated to this PFParticle
     std::vector<art::Ptr<recob::Track>>  pfPartTrackVec  = trackAssns.at(pfParticle.key());
     std::vector<art::Ptr<recob::Vertex>> pfPartVertexVec = vertexAssns.at(pfParticle.key());
 
     nTracks    += pfPartTrackVec.size();
     nVertices  += pfPartVertexVec.size();
     nDaughters += pfParticle->Daughters().size();
 
     for(auto& daughterIdx : pfParticle->Daughters())
     {
         nDaughters += traversePFParticleHierarchy(pfParticleHandle, daughterIdx, trackAssns, vertexAssns, nTracks, nVertices);
     }
 */
     return nDaughters;
 }

Member Data Documentation

std::vector<std::vector<double> > TrackHitAna::TrackHitAna::fChannelPedVec

private

Definition at line 133 of file TrackHitAna_module.cc.

int TrackHitAna::TrackHitAna::fEvent

private

Definition at line 125 of file TrackHitAna_module.cc.

const geo::GeometryCore* TrackHitAna::TrackHitAna::fGeometry

private

Definition at line 136 of file TrackHitAna_module.cc.

std::vector<std::unique_ptr<IHitHistogramTool> > TrackHitAna::TrackHitAna::fHitHistogramToolVec

private

Definition at line 131 of file TrackHitAna_module.cc.

std::vector<art::InputTag> TrackHitAna::TrackHitAna::fHitProducerLabelVec

private

Definition at line 119 of file TrackHitAna_module.cc.

int TrackHitAna::TrackHitAna::fNumEvents

private

Definition at line 128 of file TrackHitAna_module.cc.

const lariov::DetPedestalProvider& TrackHitAna::TrackHitAna::fPedestalRetrievalAlg

private

Keep track of an instance to the pedestal retrieval alg.

Definition at line 137 of file TrackHitAna_module.cc.

std::vector<art::InputTag> TrackHitAna::TrackHitAna::fPFParticleProducerLabelVec

private

Definition at line 121 of file TrackHitAna_module.cc.

int TrackHitAna::TrackHitAna::fRun

private

Definition at line 126 of file TrackHitAna_module.cc.

int TrackHitAna::TrackHitAna::fSubRun

private

Definition at line 127 of file TrackHitAna_module.cc.

std::vector<art::InputTag> TrackHitAna::TrackHitAna::fTrackProducerLabelVec

private

Definition at line 122 of file TrackHitAna_module.cc.

std::vector<art::InputTag> TrackHitAna::TrackHitAna::fWireProducerLabelVec

private

Definition at line 120 of file TrackHitAna_module.cc.

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

TrackHitAna_module.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation