quad::EvalVtx Class Reference

Inheritance diagram for quad::EvalVtx:

Public Member Functions
	EvalVtx (const fhicl::ParameterSet &pset)
void	analyze (const art::Event &evt) override

Private Attributes
std::string	fTruthLabel
std::vector< std::string >	fVertexLabels
TTree *	fTree
int	gEvt
double	gTrueX
double	gTrueY
double	gTrueZ
std::vector< double >	gRecoX
std::vector< double >	gRecoY
std::vector< double >	gRecoZ

Detailed Description

Definition at line 21 of file EvalVtx_module.cc.

Constructor & Destructor Documentation

quad::EvalVtx::EvalVtx ( const fhicl::ParameterSet &  pset )

explicit

Definition at line 43 of file EvalVtx_module.cc.

                                              :
  EDAnalyzer(pset),
  fTruthLabel(pset.get<std::string>("TruthLabel")),
  fVertexLabels(pset.get<std::vector<std::string>>("VertexLabels"))
{
  art::ServiceHandle<art::TFileService> tfs;

  fTree = tfs->make<TTree>("vtxs", "vtxs");

  fTree->Branch("evt", &gEvt);
  fTree->Branch("true_x", &gTrueX);
  fTree->Branch("true_y", &gTrueY);
  fTree->Branch("true_z", &gTrueZ);

  gRecoX.resize(fVertexLabels.size());
  gRecoY.resize(fVertexLabels.size());
  gRecoZ.resize(fVertexLabels.size());

  for(unsigned int i = 0; i < fVertexLabels.size(); ++i){
    const std::string& l = fVertexLabels[i];
    fTree->Branch((l+"_x").c_str(), &gRecoX[i]);
    fTree->Branch((l+"_y").c_str(), &gRecoY[i]);
    fTree->Branch((l+"_z").c_str(), &gRecoZ[i]);
  }
}

Member Function Documentation

void quad::EvalVtx::analyze ( const art::Event &  evt )

override

Definition at line 111 of file EvalVtx_module.cc.

{
  const auto& truths = *evt.getValidHandle<std::vector<simb::MCTruth>>(fTruthLabel);
  if(truths.empty()) return;

  const simb::MCParticle& nu = truths[0].GetNeutrino().Nu();

  gEvt = evt.event();

  gTrueX = nu.Vx();
  gTrueY = nu.Vy();
  gTrueZ = nu.Vz();

  for(unsigned int i = 0; i < fVertexLabels.size(); ++i){
    // NB - vertices returned by these functions can be invalid (default
    // constructed) in case no vertex exists. We stil have to write them, since
    // another algorithm may have made a valid vertex for this
    // event. Downstream plotting code should be aware and treat (0,0,0)
    // vertices specially.
    const recob::Vertex vtx = (fVertexLabels[i] == "pandora") ? GetVtxByAssns(fVertexLabels[i], evt) : GetFirstVertex(fVertexLabels[i], evt);
    const recob::Vertex::Point_t reco_vtx = vtx.position();

    gRecoX[i] = reco_vtx.x();
    gRecoY[i] = reco_vtx.y();
    gRecoZ[i] = reco_vtx.z();
  } // end for i

  fTree->Fill();
}

Member Data Documentation

TTree* quad::EvalVtx::fTree

private

Definition at line 32 of file EvalVtx_module.cc.

std::string quad::EvalVtx::fTruthLabel

private

Definition at line 29 of file EvalVtx_module.cc.

std::vector<std::string> quad::EvalVtx::fVertexLabels

private

Definition at line 30 of file EvalVtx_module.cc.

int quad::EvalVtx::gEvt

private

Definition at line 35 of file EvalVtx_module.cc.

std::vector<double> quad::EvalVtx::gRecoX

private

Definition at line 37 of file EvalVtx_module.cc.

std::vector<double> quad::EvalVtx::gRecoY

private

Definition at line 37 of file EvalVtx_module.cc.

std::vector<double> quad::EvalVtx::gRecoZ

private

Definition at line 37 of file EvalVtx_module.cc.

double quad::EvalVtx::gTrueX

private

Definition at line 36 of file EvalVtx_module.cc.

double quad::EvalVtx::gTrueY

private

Definition at line 36 of file EvalVtx_module.cc.

double quad::EvalVtx::gTrueZ

private

Definition at line 36 of file EvalVtx_module.cc.

---

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

• EvalVtx_module.cc