evgb::MCTruthAndFriendsItr Class Reference

#include "MCTruthAndFriendsItr.h"

Public Member Functions
	MCTruthAndFriendsItr (art::Event const &evtIn, std::vector< std::string > const &labels)
virtual	~MCTruthAndFriendsItr ()
bool	Next ()
const simb::MCTruth *	GetMCTruth () const
const simb::GTruth *	GetGTruth () const
const simb::MCFlux *	GetMCFlux () const
const bsim::Dk2Nu *	GetDk2Nu () const
const bsim::NuChoice *	GetNuChoice () const
std::string	GetLabel () const

Private Attributes
art::Event const &	evt
std::vector< std::string > const &	fInputModuleLabels
std::vector< art::Handle< std::vector< simb::MCTruth > > >	mclists
std::set< std::pair< int, int > >	indices
std::set< std::pair< int, int > >::const_iterator	indx_itr
std::vector< std::string >	outlabels
int	nmctruth
int	imctruth
const simb::MCTruth *	thisMCTruth
const simb::GTruth *	thisGTruth
const simb::MCFlux *	thisMCFlux
const bsim::Dk2Nu *	thisDk2Nu
const bsim::NuChoice *	thisNuChoice
std::string	thisLabel

Detailed Description

Definition at line 32 of file MCTruthAndFriendsItr.h.

Constructor & Destructor Documentation

evgb::MCTruthAndFriendsItr::MCTruthAndFriendsItr	(	art::Event const &	evtIn,
		std::vector< std::string > const &	labels
	)

std::cout << "+++ mcl " << mcl << "[" << mclists.size() << "] " << "nmc " << nmc << "[" << mclistHandle->size() << "] " << std::endl; std::cout << *(mct.get()) << std::endl;

Definition at line 15 of file MCTruthAndFriendsItr.cxx.

References evt, fInputModuleLabels, art::ProductRetriever::getByLabel(), art::ProductRetriever::getMany(), indices, indx_itr, art::Handle< T >::isValid(), mclists, nmctruth, outlabels, and art::Handle< T >::provenance().

  : evt(evtIn)
  , fInputModuleLabels(labels)
  , indx_itr(indices.begin())
  , nmctruth(0)
  , imctruth(-1)
  , thisMCTruth(0)
  , thisGTruth(0)
  , thisMCFlux(0)
  , thisDk2Nu(0)
  , thisNuChoice(0)
  , thisLabel("")
{

  // look for any existing MCTruth info in this event
  if ( fInputModuleLabels.size()==0) {
    mclists = evt.getMany<std::vector<simb::MCTruth>>();
    // std::cout << "evt.getManyByType" << std::endl;
  } else {
    mclists.resize(fInputModuleLabels.size());
    for (size_t i=0; i<fInputModuleLabels.size(); ++i) {
      evt.getByLabel(fInputModuleLabels[i],mclists[i]);
      //std::cout << "evt.getByLabel " << fInputModuleLabels[i] << std::endl;
    }
  }

  //std::cerr << "evg::GENIEDumper::analyze got stuff ---------------- "
  //          << evt.id() << std::endl;

  for (size_t mcl = 0; mcl < mclists.size(); ++mcl) {
    art::Handle< std::vector<simb::MCTruth> > mclistHandle = mclists[mcl];
    if ( ! mclistHandle.isValid() ) {
      // std::cerr << "not valid mcl=" <<  mcl << "---------------- " << std::endl;
      continue;
    }
    // processName, moduleLabel, instance (productInstanceName?)
    /*
    std::cout << " mclistHandle processName '"
              << mclistHandle.provenance()->processName()  // top of fcl file
              << "' moduleLabel '"
              << mclistHandle.provenance()->moduleLabel()
              << "' productInstanceName '"
              << mclistHandle.provenance()->productInstanceName()
              << "'"
              << std::endl;
    */

    std::string handleLabel = mclistHandle.provenance()->moduleLabel();
    outlabels.push_back(handleLabel);
    /*
    std::cerr << "mcl=" <<  mcl << " '" << handleLabel << "' ---------------- " << std::endl;
    */

    // loop over mctruths in a list
    for(size_t nmc = 0; nmc < mclistHandle->size(); ++nmc) {
      art::Ptr<simb::MCTruth> mct(mclistHandle, nmc);

      std::pair<int,int> ipair(mcl,nmc);
      indices.insert(ipair);

    }

  }

  indx_itr = indices.begin();
  nmctruth = (int)(indices.size());
  //std::cout << ".... found nmctruth " << nmctruth
  //          << " nlabels " << outlabels.size()
  //          << std::endl;

}

virtual evgb::MCTruthAndFriendsItr::~MCTruthAndFriendsItr ( )

inlinevirtual

Definition at line 38 of file MCTruthAndFriendsItr.h.

References Next().

{ ; }

Member Function Documentation

const bsim::Dk2Nu* evgb::MCTruthAndFriendsItr::GetDk2Nu ( ) const

inline

Definition at line 44 of file MCTruthAndFriendsItr.h.

References thisDk2Nu.

Referenced by evwgh::PPFXCVWeightCalc::GetWeight(), evwgh::PPFXWeightCalc::GetWeight(), evwgh::PPFXTotAbsorpWeightCalc::GetWeight(), evwgh::PPFXThinPionWeightCalc::GetWeight(), evwgh::PPFXThinNucWeightCalc::GetWeight(), evwgh::PPFXThinNucAWeightCalc::GetWeight(), evwgh::PPFXThinNeutronPionWeightCalc::GetWeight(), evwgh::PPFXThinMesonWeightCalc::GetWeight(), evwgh::PPFXThinKaonWeightCalc::GetWeight(), evwgh::PPFXTargAttenWeightCalc::GetWeight(), evwgh::PPFXOtherWeightCalc::GetWeight(), evwgh::PPFXMIPPPionWeightCalc::GetWeight(), and evwgh::PPFXMIPPKaonWeightCalc::GetWeight().

{ return thisDk2Nu;   }

const simb::GTruth* evgb::MCTruthAndFriendsItr::GetGTruth ( ) const

inline

Definition at line 42 of file MCTruthAndFriendsItr.h.

References thisGTruth.

Referenced by evg::GenieOutput::analyze().

{ return thisGTruth;  }

std::string evgb::MCTruthAndFriendsItr::GetLabel ( ) const

inline

Definition at line 47 of file MCTruthAndFriendsItr.h.

References thisLabel.

Referenced by evg::GenieOutput::analyze(), evwgh::PPFXThinMesonWeightCalc::GetWeight(), evwgh::PPFXWeightCalc::GetWeight(), evwgh::PPFXTotAbsorpWeightCalc::GetWeight(), evwgh::PPFXThinPionWeightCalc::GetWeight(), evwgh::PPFXThinNucWeightCalc::GetWeight(), evwgh::PPFXThinNucAWeightCalc::GetWeight(), evwgh::PPFXThinNeutronPionWeightCalc::GetWeight(), evwgh::PPFXMIPPPionWeightCalc::GetWeight(), evwgh::PPFXThinKaonWeightCalc::GetWeight(), evwgh::PPFXTargAttenWeightCalc::GetWeight(), evwgh::PPFXOtherWeightCalc::GetWeight(), evwgh::PPFXCVWeightCalc::GetWeight(), and evwgh::PPFXMIPPKaonWeightCalc::GetWeight().

{ return thisLabel;   } 

const simb::MCFlux* evgb::MCTruthAndFriendsItr::GetMCFlux ( ) const

inline

Definition at line 43 of file MCTruthAndFriendsItr.h.

References thisMCFlux.

Referenced by evg::GenieOutput::analyze().

{ return thisMCFlux;  }

const simb::MCTruth* evgb::MCTruthAndFriendsItr::GetMCTruth ( ) const

inline

Definition at line 41 of file MCTruthAndFriendsItr.h.

References thisMCTruth.

Referenced by evg::GenieOutput::analyze(), evwgh::PPFXThinMesonWeightCalc::GetWeight(), evwgh::PPFXWeightCalc::GetWeight(), evwgh::PPFXTotAbsorpWeightCalc::GetWeight(), evwgh::PPFXThinPionWeightCalc::GetWeight(), evwgh::PPFXThinNucWeightCalc::GetWeight(), evwgh::PPFXThinNucAWeightCalc::GetWeight(), evwgh::PPFXThinNeutronPionWeightCalc::GetWeight(), evwgh::PPFXMIPPPionWeightCalc::GetWeight(), evwgh::PPFXThinKaonWeightCalc::GetWeight(), evwgh::PPFXTargAttenWeightCalc::GetWeight(), evwgh::PPFXOtherWeightCalc::GetWeight(), evwgh::PPFXCVWeightCalc::GetWeight(), and evwgh::PPFXMIPPKaonWeightCalc::GetWeight().

{ return thisMCTruth; }

const bsim::NuChoice* evgb::MCTruthAndFriendsItr::GetNuChoice ( ) const

inline

Definition at line 45 of file MCTruthAndFriendsItr.h.

References thisNuChoice.

{ return thisNuChoice; }

bool evgb::MCTruthAndFriendsItr::Next

(

)

std::cout << "imctruth " << std::setw(3) << imctruth << " [" << indx_handle << "," << indx_within << "]" << " hvMCTruth.isValid() " << hvMCTruth.isValid() << " '" << outlabels[indx_handle] << "' " << std::endl;

Definition at line 94 of file MCTruthAndFriendsItr.cxx.

References evt, art::Ptr< T >::get(), imctruth, indx_itr, mclists, nmctruth, outlabels, thisDk2Nu, thisGTruth, thisLabel, thisMCFlux, thisMCTruth, and thisNuChoice.

Referenced by evg::GenieOutput::analyze(), evwgh::PPFXThinMesonWeightCalc::GetWeight(), evwgh::PPFXMIPPKaonWeightCalc::GetWeight(), evwgh::PPFXWeightCalc::GetWeight(), evwgh::PPFXTotAbsorpWeightCalc::GetWeight(), evwgh::PPFXThinPionWeightCalc::GetWeight(), evwgh::PPFXThinNucWeightCalc::GetWeight(), evwgh::PPFXThinNucAWeightCalc::GetWeight(), evwgh::PPFXOtherWeightCalc::GetWeight(), evwgh::PPFXMIPPPionWeightCalc::GetWeight(), evwgh::PPFXThinKaonWeightCalc::GetWeight(), evwgh::PPFXTargAttenWeightCalc::GetWeight(), evwgh::PPFXThinNeutronPionWeightCalc::GetWeight(), evwgh::PPFXCVWeightCalc::GetWeight(), and ~MCTruthAndFriendsItr().

{
  ++imctruth;   // started at -1, so first call to Next() prepares us for indx=0
  thisMCTruth  = 0;
  thisGTruth   = 0;
  thisMCFlux   = 0;
  thisDk2Nu    = 0;
  thisNuChoice = 0;
  //std::cerr << "Next() called ... moved to imctruth " << imctruth << std::endl;
  if ( imctruth >= nmctruth ) return false;

  size_t indx_handle = (*indx_itr).first;
  size_t indx_within = (*indx_itr).second;

  thisLabel = outlabels[indx_handle];

  art::Handle< std::vector<simb::MCTruth> > hvMCTruth = mclists[indx_handle];

  thisMCTruth = &(hvMCTruth->at(indx_within));

  // inefficient ... should only need to do this bit for every new
  // Handle ...

  //art::FindOneP<recob::Hit> findSPToHits(fSpacePoints,evt,fSpacePointLabel);
  //const art::Ptr<recob::Hit> hit = findSPToHits.at(iSP);

  try {
    art::FindOneP<simb::GTruth> qgtruth(hvMCTruth,evt,outlabels[indx_handle]);
    const art::Ptr<simb::GTruth> gtruthptr = qgtruth.at(indx_within);
    thisGTruth = gtruthptr.get();
  }
  catch (...) {
    // std::cerr << "no GTruth for this handle" << std::endl;
  }

  try {
    art::FindOneP<simb::MCFlux> qmcflux(hvMCTruth,evt,outlabels[indx_handle]);
    const art::Ptr<simb::MCFlux> mcfluxptr = qmcflux.at(indx_within);
    thisMCFlux = mcfluxptr.get();
  }
  catch (...) {
    // std::cerr << "no MCFlux for this handle" << std::endl;
  }

  try {
    art::FindOneP<bsim::Dk2Nu> qdk2nu(hvMCTruth,evt,outlabels[indx_handle]);
    const art::Ptr<bsim::Dk2Nu> dk2nuptr = qdk2nu.at(indx_within);
    thisDk2Nu = dk2nuptr.get();
  }
  catch (...) {
    // std::cerr << "no bsim::Dk2NU for this handle" << std::endl;
  }

  try {
    art::FindOneP<bsim::NuChoice> qnuchoice(hvMCTruth,evt,outlabels[indx_handle]);
    const art::Ptr<bsim::NuChoice> nuchoiceptr = qnuchoice.at(indx_within);
    thisNuChoice = nuchoiceptr.get();
  }
  catch (...) {
    // std::cerr << "no bsim::NuChoice for this handle" << std::endl;
  }

  //std::cerr << "Next() called ... seeing " << thisMCTruth
  //          << " " << thisGTruth << " " << thisMCFlux << std::endl;

  // so user code looks like
  // evgb::MCTruthAndFriendsItr mcitr(evt,labels);
  // while ( mcitr.Next() ) {
  //    const simb::MCTruth* amctruth = mcitr.GetMCTruth();
  //    const simb::GTruth*  agtruth  = mcitr.GetGTruth();
  //...

  /*
  // loop over lists
  try {
    //art::FindOneP<simb::GTruth> QueryG(mclistHandle,evt,handleLabel);
  }
  catch (art::Exception) {
    //std::cerr << "no GTruth for " << handleLabel << std::endl;
  }
  */

  // move the iterator on
  ++indx_itr;
  return true;
}

Member Data Documentation

art::Event const& evgb::MCTruthAndFriendsItr::evt

private

Definition at line 52 of file MCTruthAndFriendsItr.h.

Referenced by MCTruthAndFriendsItr(), and Next().

std::vector<std::string> const& evgb::MCTruthAndFriendsItr::fInputModuleLabels

private

Definition at line 53 of file MCTruthAndFriendsItr.h.

Referenced by MCTruthAndFriendsItr().

int evgb::MCTruthAndFriendsItr::imctruth

private

Definition at line 62 of file MCTruthAndFriendsItr.h.

Referenced by Next().

std::set<std::pair<int,int> > evgb::MCTruthAndFriendsItr::indices

private

Definition at line 57 of file MCTruthAndFriendsItr.h.

Referenced by MCTruthAndFriendsItr().

std::set<std::pair<int,int> >::const_iterator evgb::MCTruthAndFriendsItr::indx_itr

private

Definition at line 58 of file MCTruthAndFriendsItr.h.

Referenced by MCTruthAndFriendsItr(), and Next().

std::vector< art::Handle< std::vector<simb::MCTruth> > > evgb::MCTruthAndFriendsItr::mclists

private

Definition at line 55 of file MCTruthAndFriendsItr.h.

Referenced by MCTruthAndFriendsItr(), and Next().

int evgb::MCTruthAndFriendsItr::nmctruth

private

Definition at line 61 of file MCTruthAndFriendsItr.h.

Referenced by MCTruthAndFriendsItr(), and Next().

std::vector<std::string> evgb::MCTruthAndFriendsItr::outlabels

private

Definition at line 59 of file MCTruthAndFriendsItr.h.

Referenced by MCTruthAndFriendsItr(), and Next().

const bsim::Dk2Nu* evgb::MCTruthAndFriendsItr::thisDk2Nu

private

Definition at line 67 of file MCTruthAndFriendsItr.h.

Referenced by GetDk2Nu(), and Next().

const simb::GTruth* evgb::MCTruthAndFriendsItr::thisGTruth

private

Definition at line 65 of file MCTruthAndFriendsItr.h.

Referenced by GetGTruth(), and Next().

std::string evgb::MCTruthAndFriendsItr::thisLabel

private

Definition at line 69 of file MCTruthAndFriendsItr.h.

Referenced by GetLabel(), and Next().

const simb::MCFlux* evgb::MCTruthAndFriendsItr::thisMCFlux

private

Definition at line 66 of file MCTruthAndFriendsItr.h.

Referenced by GetMCFlux(), and Next().

const simb::MCTruth* evgb::MCTruthAndFriendsItr::thisMCTruth

private

Definition at line 64 of file MCTruthAndFriendsItr.h.

Referenced by GetMCTruth(), and Next().

const bsim::NuChoice* evgb::MCTruthAndFriendsItr::thisNuChoice

private

Definition at line 68 of file MCTruthAndFriendsItr.h.

Referenced by GetNuChoice(), and Next().

---

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

• nugen/EventGeneratorBase/GENIE/MCTruthAndFriendsItr.h
• nugen/EventGeneratorBase/GENIE/MCTruthAndFriendsItr.cxx