cluster::ClusterPCA Class Reference

Inheritance diagram for cluster::ClusterPCA:

Public Member Functions
	ClusterPCA (fhicl::ParameterSet const &pset)
	~ClusterPCA ()

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

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 )

private

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

private

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:

• ClusterPCA_module.cc