calo::TrackCalorimetryAlg Class Reference

#include "TrackCalorimetryAlg.h"

Classes
struct	HitProperties
struct	HitPropertySorter

Public Types
using	Providers_t = lar::ProviderPack< geo::GeometryCore, detinfo::LArProperties, detinfo::DetectorProperties >

Public Member Functions
	TrackCalorimetryAlg (fhicl::ParameterSet const &p)
void	reconfigure (fhicl::ParameterSet const &p)
void	ExtractCalorimetry (std::vector< recob::Track > const &, std::vector< recob::Hit > const &, std::vector< std::vector< size_t > > const &, std::vector< anab::Calorimetry > &, std::vector< size_t > &, Providers_t providers)

Private Types
typedef std::multiset< HitProperties, HitPropertySorter >	HitPropertiesMultiset_t

Private Member Functions
void	ClearInternalVectors ()
void	ReserveInternalVectors (size_t s)
std::vector< float >	CreatePathLengthFractionVector (recob::Track const &track)
void	AnalyzeHit (recob::Hit const &, recob::Track const &, std::vector< std::pair< geo::WireID, float > > const &, std::vector< float > const &, HitPropertiesMultiset_t &, geo::GeometryCore const &)
bool	IsInvertedTrack (HitPropertiesMultiset_t const &)
void	MakeCalorimetryObject (HitPropertiesMultiset_t const &hpm, recob::Track const &track, size_t const &i_track, std::vector< anab::Calorimetry > &caloVector, std::vector< size_t > &assnTrackCaloVector, geo::PlaneID const &planeID)
void	PrintHitPropertiesMultiset (HitPropertiesMultiset_t const &hpm)

Private Attributes
CalorimetryAlg	caloAlg
unsigned int	fNHitsToDetermineStart

Detailed Description

Definition at line 34 of file TrackCalorimetryAlg.h.

Member Typedef Documentation

typedef std::multiset<HitProperties,HitPropertySorter> calo::TrackCalorimetryAlg::HitPropertiesMultiset_t

private

Definition at line 80 of file TrackCalorimetryAlg.h.

using calo::TrackCalorimetryAlg::Providers_t = lar::ProviderPack< geo::GeometryCore, detinfo::LArProperties, detinfo::DetectorProperties >

Definition at line 38 of file TrackCalorimetryAlg.h.

Constructor & Destructor Documentation

calo::TrackCalorimetryAlg::TrackCalorimetryAlg

(

fhicl::ParameterSet const &

p

)

Title: Track Calorimetry Algorithim Class Author: Wes Ketchum (wketc.nosp@m.hum@.nosp@m.lanl..nosp@m.gov), based on code in the Calorimetry_module

Description: Algorithm that produces a calorimetry object given a track Input: recob::Track, Assn<recob::Spacepoint,recob::Track>, Assn<recob::Hit,recob::Track> Output: anab::Calorimetry, (and Assn<anab::Calorimetry,recob::Track>)

Definition at line 16 of file TrackCalorimetryAlg.cxx.

References reconfigure().

                                                                      :
  caloAlg(p.get<fhicl::ParameterSet>("CalorimetryAlg"))
{
  this->reconfigure(p);
}

Member Function Documentation

void calo::TrackCalorimetryAlg::AnalyzeHit ( recob::Hit const &  hit, recob::Track const &  track, std::vector< std::pair< geo::WireID, float > > const &  traj_points_in_plane, std::vector< float > const &  path_length_fraction_vec, HitPropertiesMultiset_t &  HitPropertiesMultiset, geo::GeometryCore const &  geom  )

private

Definition at line 119 of file TrackCalorimetryAlg.cxx.

References caloAlg, calo::CalorimetryAlg::dEdx_AREA(), recob::Hit::Integral(), recob::Track::LocationAtPoint(), geo::PlaneID::Plane, lar::util::TrackPitchInView(), geo::GeometryCore::View(), and recob::Hit::WireID().

Referenced by ExtractCalorimetry(), and ReserveInternalVectors().

                                                                       {

  size_t traj_iter = std::distance(traj_points_in_plane.begin(),
                                   std::min_element(traj_points_in_plane.begin(),
                                                    traj_points_in_plane.end(),
                                                    dist_projected(hit,geom)));
  float pitch = lar::util::TrackPitchInView(track, geom.View(hit.WireID().Plane),traj_iter);

  HitPropertiesMultiset.emplace(hit.Integral(),
                                hit.Integral()/pitch,
                                caloAlg.dEdx_AREA(hit,pitch),
                                pitch,
                                track.LocationAtPoint(traj_iter),
                                path_length_fraction_vec[traj_iter]);
}

void calo::TrackCalorimetryAlg::ClearInternalVectors ( )

inlineprivate

Definition at line 83 of file TrackCalorimetryAlg.h.

Referenced by ExtractCalorimetry().

{}

std::vector< float > calo::TrackCalorimetryAlg::CreatePathLengthFractionVector ( recob::Track const &  track )

private

Definition at line 105 of file TrackCalorimetryAlg.cxx.

References recob::Track::Length(), recob::Track::LocationAtPoint(), and recob::Track::NumberTrajectoryPoints().

Referenced by ExtractCalorimetry(), and ReserveInternalVectors().

                                                                                               {

  std::vector<float> trk_path_length_frac_vec(track.NumberTrajectoryPoints());

  float cumulative_path_length=0;
  const float total_path_length = track.Length();
  for(size_t i_trj=1; i_trj<track.NumberTrajectoryPoints(); i_trj++){    
    cumulative_path_length+=(track.LocationAtPoint(i_trj)-track.LocationAtPoint(i_trj-1)).Mag();
    trk_path_length_frac_vec[i_trj]=cumulative_path_length/total_path_length;
  }

  return trk_path_length_frac_vec;
}

void calo::TrackCalorimetryAlg::ExtractCalorimetry	(	std::vector< recob::Track > const &	trackVector,
		std::vector< recob::Hit > const &	hitVector,
		std::vector< std::vector< size_t > > const &	hit_indices_per_track,
		std::vector< anab::Calorimetry > &	caloVector,
		std::vector< size_t > &	assnTrackCaloVector,
		Providers_t	providers
	)

Definition at line 27 of file TrackCalorimetryAlg.cxx.

References AnalyzeHit(), ClearInternalVectors(), CreatePathLengthFractionVector(), lar::ProviderPack< Providers >::get(), recob::Track::LocationAtPoint(), MakeCalorimetryObject(), recob::Track::NumberTrajectoryPoints(), ReserveInternalVectors(), and track.

Referenced by calo::TrackCalorimetry::produce().

{
  auto const& geom = *(providers.get<geo::GeometryCore>());
//  auto const& larp = *(providers.get<detinfo::LArProperties>());
  auto const& detprop = *(providers.get<detinfo::DetectorProperties>());
  
  //loop over the track list
  for(size_t i_track=0; i_track<trackVector.size(); i_track++){

    recob::Track const& track = trackVector[i_track];
    std::vector<float> path_length_fraction_vec(CreatePathLengthFractionVector(track));

    //sort hits into each plane
    std::vector< std::vector<size_t> > hit_indices_per_plane(geom.Nplanes());
    for(auto const& i_hit : hit_indices_per_track[i_track])
      hit_indices_per_plane[hitVector[i_hit].WireID().Plane].push_back(i_hit);
    
    //loop over the planes
    for(size_t i_plane=0; i_plane<geom.Nplanes(); i_plane++){

      ClearInternalVectors();
      ReserveInternalVectors(hit_indices_per_plane[i_plane].size());

      //project down the track into wire/tick space for this plane
      std::vector< std::pair<geo::WireID,float> > traj_points_in_plane(track.NumberTrajectoryPoints());
      for(size_t i_trjpt=0; i_trjpt<track.NumberTrajectoryPoints(); i_trjpt++){
        double x_pos = track.LocationAtPoint(i_trjpt).X();
        float tick = detprop.ConvertXToTicks(x_pos,(int)i_plane,0,0);
        traj_points_in_plane[i_trjpt] = std::make_pair(geom.NearestWireID(track.LocationAtPoint(i_trjpt),i_plane),
                                                       tick);
      }
      
      HitPropertiesMultiset_t HitPropertiesMultiset;
      //now loop through hits
      for(auto const& i_hit : hit_indices_per_plane[i_plane])
        AnalyzeHit(hitVector[i_hit],
                   track,
                   traj_points_in_plane,
                   path_length_fraction_vec,
                   HitPropertiesMultiset,
                   geom);


      //PrintHitPropertiesMultiset(HitPropertiesMultiset);
      geo::PlaneID planeID(0,0,i_plane);
      MakeCalorimetryObject(HitPropertiesMultiset, track, i_track, caloVector, assnTrackCaloVector, planeID);

    }//end loop over planes

  }//end loop over tracks

}//end ExtractCalorimetry

bool calo::TrackCalorimetryAlg::IsInvertedTrack ( HitPropertiesMultiset_t const &  hpm )

private

Definition at line 140 of file TrackCalorimetryAlg.cxx.

References fNHitsToDetermineStart.

Referenced by MakeCalorimetryObject(), PrintHitPropertiesMultiset(), and ReserveInternalVectors().

                                                                               {

  if(hpm.size() <= fNHitsToDetermineStart) return false;
  
  float charge_start=0,charge_end=0;
  unsigned int counter=0;
  for(HitPropertiesMultiset_t::iterator it_hpm=hpm.begin(); 
      it_hpm!=hpm.end();
      it_hpm++)
    {
      charge_start += (*it_hpm).charge;
      counter++;
      if(counter==fNHitsToDetermineStart) break;
    }

  counter=0;
  for(HitPropertiesMultiset_t::reverse_iterator it_hpm=hpm.rbegin(); 
      it_hpm!=hpm.rend();
      it_hpm++)
    {
      charge_end += (*it_hpm).charge;
      counter++;
      if(counter==fNHitsToDetermineStart) break;
    }

  return (charge_start > charge_end);
}

void calo::TrackCalorimetryAlg::MakeCalorimetryObject ( HitPropertiesMultiset_t const &  hpm, recob::Track const &  track, size_t const &  i_track, std::vector< anab::Calorimetry > &  caloVector, std::vector< size_t > &  assnTrackCaloVector, geo::PlaneID const &  planeID  )

private

Definition at line 168 of file TrackCalorimetryAlg.cxx.

References IsInvertedTrack(), and recob::Track::Length().

Referenced by ExtractCalorimetry(), and ReserveInternalVectors().

                                                                                {
  size_t n_hits = hpm.size();
  std::vector<double> dEdxVector,dQdxVector,resRangeVector,deadWireVector,pitchVector;
  std::vector<TVector3> XYZVector;

  dEdxVector.reserve(n_hits);
  dQdxVector.reserve(n_hits);
  resRangeVector.reserve(n_hits);
  deadWireVector.reserve(n_hits);
  pitchVector.reserve(n_hits);
  XYZVector.reserve(n_hits);

  double kinetic_energy=0,track_length=track.Length();
  if(IsInvertedTrack(hpm)){

    for(HitPropertiesMultiset_t::iterator it_hpm=hpm.begin(); 
        it_hpm!=hpm.end();
        it_hpm++)
      {
        dEdxVector.push_back((*it_hpm).dEdx);
        dQdxVector.push_back((*it_hpm).dQdx);
        resRangeVector.push_back((*it_hpm).path_fraction*track_length);
        pitchVector.push_back((*it_hpm).pitch);
        XYZVector.push_back((*it_hpm).xyz);
        kinetic_energy += dEdxVector.back()*pitchVector.back();
      }
        
  }
  else{

    for(HitPropertiesMultiset_t::reverse_iterator it_hpm=hpm.rbegin(); 
        it_hpm!=hpm.rend();
        it_hpm++)
      {
        dEdxVector.push_back((*it_hpm).dEdx);
        dQdxVector.push_back((*it_hpm).dQdx);
        resRangeVector.push_back((1-(*it_hpm).path_fraction)*track_length);
        pitchVector.push_back((*it_hpm).pitch);
        XYZVector.push_back((*it_hpm).xyz);
        kinetic_energy += dEdxVector.back()*pitchVector.back();
      }

  }

  caloVector.emplace_back(kinetic_energy,
                          dEdxVector,
                          dQdxVector,
                          resRangeVector,
                          deadWireVector,
                          track_length,
                          pitchVector,
                          XYZVector,
                          planeID);
  assnTrackCaloVector.emplace_back(i_track);
}

void calo::TrackCalorimetryAlg::PrintHitPropertiesMultiset ( HitPropertiesMultiset_t const &  hpm )

private

Definition at line 229 of file TrackCalorimetryAlg.cxx.

References IsInvertedTrack().

Referenced by ReserveInternalVectors().

                                                                                          {

  for(auto const& hit : hpm)
    hit.Print();

  std::cout << "Inverted? " << IsInvertedTrack(hpm) << std::endl;

}

void calo::TrackCalorimetryAlg::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 22 of file TrackCalorimetryAlg.cxx.

References caloAlg, fNHitsToDetermineStart, fhicl::ParameterSet::get(), and calo::CalorimetryAlg::reconfigure().

Referenced by calo::TrackCalorimetry::reconfigure(), and TrackCalorimetryAlg().

                                                                   {
  caloAlg.reconfigure(p.get<fhicl::ParameterSet>("CalorimetryAlg"));
  fNHitsToDetermineStart = p.get<unsigned int>("NHitsToDetermineStart",3);
}

void calo::TrackCalorimetryAlg::ReserveInternalVectors ( size_t  s )

inlineprivate

Definition at line 84 of file TrackCalorimetryAlg.h.

References AnalyzeHit(), CreatePathLengthFractionVector(), IsInvertedTrack(), MakeCalorimetryObject(), PrintHitPropertiesMultiset(), track, and lar::dump::vector().

Referenced by ExtractCalorimetry().

{}

Member Data Documentation

CalorimetryAlg calo::TrackCalorimetryAlg::caloAlg

private

Definition at line 53 of file TrackCalorimetryAlg.h.

Referenced by AnalyzeHit(), and reconfigure().

unsigned int calo::TrackCalorimetryAlg::fNHitsToDetermineStart

private

Definition at line 54 of file TrackCalorimetryAlg.h.

Referenced by IsInvertedTrack(), and reconfigure().

---

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

• larreco/v06_64_02/source/larreco/Calorimetry/TrackCalorimetryAlg.h
• larreco/v06_64_02/source/larreco/Calorimetry/TrackCalorimetryAlg.cxx