#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.

   : 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
   mclists.clear();
 
   if ( fInputModuleLabels.size()==0) {
     evt.getManyByType(mclists);
     // 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);
 
       /**
       std::cout << "+++ mcl " << mcl << "[" << mclists.size() << "] "
                 << "nmc " << nmc << "[" << mclistHandle->size() << "] "
                 << std::endl;
       std::cout << *(mct.get()) << std::endl;
       ***/
     }
 
   }
 
   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.

38 { ; }

Member Function Documentation

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

inline

Definition at line 44 of file MCTruthAndFriendsItr.h.

44 { return thisDk2Nu; }

evgb::MCTruthAndFriendsItr::thisDk2Nu

const bsim::Dk2Nu * thisDk2Nu

Definition: MCTruthAndFriendsItr.h:67

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

inline

Definition at line 42 of file MCTruthAndFriendsItr.h.

42 { return thisGTruth; }

evgb::MCTruthAndFriendsItr::thisGTruth

const simb::GTruth * thisGTruth

Definition: MCTruthAndFriendsItr.h:65

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

inline

Definition at line 47 of file MCTruthAndFriendsItr.h.

47 { return thisLabel; }

evgb::MCTruthAndFriendsItr::thisLabel

std::string thisLabel

Definition: MCTruthAndFriendsItr.h:69

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

inline

Definition at line 43 of file MCTruthAndFriendsItr.h.

43 { return thisMCFlux; }

evgb::MCTruthAndFriendsItr::thisMCFlux

const simb::MCFlux * thisMCFlux

Definition: MCTruthAndFriendsItr.h:66

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

inline

Definition at line 41 of file MCTruthAndFriendsItr.h.

41 { return thisMCTruth; }

evgb::MCTruthAndFriendsItr::thisMCTruth

const simb::MCTruth * thisMCTruth

Definition: MCTruthAndFriendsItr.h:64

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

inline

Definition at line 45 of file MCTruthAndFriendsItr.h.

45 { return thisNuChoice; }

evgb::MCTruthAndFriendsItr::thisNuChoice

const bsim::NuChoice * thisNuChoice

Definition: MCTruthAndFriendsItr.h:68

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 96 of file MCTruthAndFriendsItr.cxx.

 {
   ++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];
 
   /**
   std::cout << "imctruth " << std::setw(3) << imctruth
             << "  [" << indx_handle << "," << indx_within << "]"
             << "  hvMCTruth.isValid() " << hvMCTruth.isValid()
             << " '" << outlabels[indx_handle] << "' "
             << std::endl;
   **/
 
   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.

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

private

Definition at line 53 of file MCTruthAndFriendsItr.h.

int evgb::MCTruthAndFriendsItr::imctruth

private

Definition at line 62 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 57 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 58 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 55 of file MCTruthAndFriendsItr.h.

int evgb::MCTruthAndFriendsItr::nmctruth

private

Definition at line 61 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 59 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 67 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 65 of file MCTruthAndFriendsItr.h.

std::string evgb::MCTruthAndFriendsItr::thisLabel

private

Definition at line 69 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 66 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 64 of file MCTruthAndFriendsItr.h.

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

private

Definition at line 68 of file MCTruthAndFriendsItr.h.

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

nutools/old/EventGeneratorBase/GENIE/MCTruthAndFriendsItr.h
nutools/old/EventGeneratorBase/GENIE/MCTruthAndFriendsItr.cxx

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation