cluster::ClusterPCA Class Reference

Inheritance diagram for cluster::ClusterPCA:

Public Member Functions
	ClusterPCA (fhicl::ParameterSet const &pset)
	~ClusterPCA ()
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	PerformClusterPCA (const std::vector< art::Ptr< recob::Hit > > &HitsThisCluster, double *PrincDirectionWT, double &PrincValue, double &TotalCharge, bool NormPC)
void	analyze (art::Event const &evt)
void	beginJob ()

Private Attributes
std::string	fClusterModuleLabel
bool	fNormPC
TTree *	fTree
Int_t	fView
Float_t	fPrincDirW
Float_t	fPrincDirT
Float_t	fPrincValue
Float_t	fTotalCharge
Float_t	fNHits

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 38 of file ClusterPCA_module.cc.

Constructor & Destructor Documentation

cluster::ClusterPCA::ClusterPCA ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 75 of file ClusterPCA_module.cc.

    : EDAnalyzer(pset)
    , fClusterModuleLabel(pset.get<std::string>("ClusterModuleLabel"))
    , fNormPC            (pset.get<bool>("NormPC"))
  {
  }

cluster::ClusterPCA::~ClusterPCA

(

)

Definition at line 83 of file ClusterPCA_module.cc.

  {
  }

Member Function Documentation

void cluster::ClusterPCA::analyze ( art::Event const &  evt )

privatevirtual

Implements art::EDAnalyzer.

Definition at line 102 of file ClusterPCA_module.cc.

  {
    // Get a Handle for the input Cluster object(s).
    art::Handle< std::vector<recob::Cluster> > clusterVecHandle;
    evt.getByLabel(fClusterModuleLabel,clusterVecHandle);

    constexpr size_t nViews = 3;

    // one vector of cluster index in the original collection, per view;
    // to support all present views, replace with a dynamically growing array
    std::array<std::vector<size_t>, nViews> ClsIndex;

    // loop over the input Clusters
    for(size_t i = 0; i < clusterVecHandle->size(); ++i){

      //get a art::Ptr to each Cluster
      art::Ptr<recob::Cluster> cl(clusterVecHandle, i);

      switch(cl->View()){
      case geo::kU :
        fView=0;
        break;
      case geo::kV :
        fView=1;
        break;
      case geo::kZ :
        fView=2;
        break;
      default :
        mf::LogError("ClusterPCA")
          << "Hit belongs to an unsupported view (#" << cl->View() << ")";
        break;
      }// end switch on view

      ClsIndex[fView].push_back(i);
    }// end loop over input clusters

    art::FindManyP<recob::Hit> fm(clusterVecHandle, evt, fClusterModuleLabel);
    for(fView = 0; fView < (int) nViews; ++fView){

      const std::vector<size_t>& ClsIndices = ClsIndex[fView];

      for(size_t c = 0; c < ClsIndices.size(); ++c){

        // find the hits associated with the current cluster
        const std::vector<art::Ptr<recob::Hit>>& ptrvs = fm.at(ClsIndices[c]);

        double PrincDir[2], PrincValue=0;
        double TotalCharge=0;

        PerformClusterPCA( ptrvs, PrincDir, PrincValue, TotalCharge, fNormPC);

        fPrincDirW   = PrincDir[0];
        fPrincDirT   = PrincDir[1];
        fPrincValue  = PrincValue;
        fTotalCharge = TotalCharge;
        fNHits       = ptrvs.size();

        fTree->Fill();

      }// end loop over first cluster iterator
     }// end loop over planes

    return;

  } // ClusterPCA::analyze()

void cluster::ClusterPCA::beginJob ( )

privatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 89 of file ClusterPCA_module.cc.

  {
    art::ServiceHandle<art::TFileService const> tfs;
    fTree = tfs->make<TTree>("PCATree","PCATree");
    fTree->Branch("View",      &fView,      "View/I");
    fTree->Branch("PrincDirW", &fPrincDirW, "PrincDirW/F");
    fTree->Branch("PrincDirT", &fPrincDirT, "PrincDirT/F");
    fTree->Branch("PrincValue",&fPrincValue, "PrincValue/F");
    fTree->Branch("TotalCharge",&fTotalCharge, "TotalCharge/F");
    fTree->Branch("NHits",      &fNHits,     "fNHits/F");
  }

void cluster::ClusterPCA::PerformClusterPCA ( const std::vector< art::Ptr< recob::Hit > > &  HitsThisCluster, double *  PrincDirectionWT, double &  PrincValue, double &  TotalCharge, bool  NormPC  )

private

Definition at line 188 of file ClusterPCA_module.cc.

  {

  double Center[2] = {0,0};
  TotalCharge = 0;

  for(auto itHit = HitsThisCluster.begin(); itHit!=HitsThisCluster.end(); ++itHit)
    {
      Center[0] += (*itHit)->WireID().Wire;
      Center[1] += (*itHit)->PeakTime();
      TotalCharge += (*itHit)->Integral();
    }

  Center[0] /= float(HitsThisCluster.size());
  Center[1] /= float(HitsThisCluster.size());


  double WireTime[2];

  std::string OptionString;

  if(NormPC==false)
    OptionString = "D";
  else
    OptionString = "ND";

  TPrincipal  pc(2, OptionString.c_str());


  for(auto itHit = HitsThisCluster.begin(); itHit!=HitsThisCluster.end(); ++itHit)
    {
      WireTime[0] = (*itHit)->WireID().Wire - Center[0];
      WireTime[1] = (*itHit)->PeakTime()    - Center[1];

      pc.AddRow(WireTime);

    }

  pc.MakePrincipals();

  PrincipalEigenvalue = (*pc.GetEigenValues())[0];

  for(size_t n=0; n!=2; ++n)
    {
      PrincipalDirection[n]=        (*pc.GetEigenVectors())[0][n];
    }


  // Comment this out if you want to shut it up
  pc.Print("MSEV");

  pc.Clear();
  return;
}

Member Data Documentation

std::string cluster::ClusterPCA::fClusterModuleLabel

private

Definition at line 52 of file ClusterPCA_module.cc.

Float_t cluster::ClusterPCA::fNHits

private

Definition at line 62 of file ClusterPCA_module.cc.

bool cluster::ClusterPCA::fNormPC

private

Definition at line 53 of file ClusterPCA_module.cc.

Float_t cluster::ClusterPCA::fPrincDirT

private

Definition at line 59 of file ClusterPCA_module.cc.

Float_t cluster::ClusterPCA::fPrincDirW

private

Definition at line 58 of file ClusterPCA_module.cc.

Float_t cluster::ClusterPCA::fPrincValue

private

Definition at line 60 of file ClusterPCA_module.cc.

Float_t cluster::ClusterPCA::fTotalCharge

private

Definition at line 61 of file ClusterPCA_module.cc.

TTree* cluster::ClusterPCA::fTree

private

Definition at line 55 of file ClusterPCA_module.cc.

Int_t cluster::ClusterPCA::fView

private

Definition at line 57 of file ClusterPCA_module.cc.

---

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

• larreco/larreco/ClusterFinder/ClusterPCA_module.cc