vertex::FeatureVertexFinderAna Class Reference

Inheritance diagram for vertex::FeatureVertexFinderAna:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer

Public Member Functions
	FeatureVertexFinderAna (fhicl::ParameterSet const &pset)
virtual	~FeatureVertexFinderAna ()
void	analyze (const art::Event &evt)
void	beginJob ()
void	reconfigure (fhicl::ParameterSet const &p)
std::string	workerType () const
bool	modifiesEvent () const
void	registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)
std::string const &	processName () const
bool	wantAllEvents () const
bool	wantEvent (Event const &e)
fhicl::ParameterSetID	selectorConfig () const
art::Handle< art::TriggerResults >	getTriggerResults (Event const &e) const
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)
base_engine_t &	createEngine (seed_t seed)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const
detail::CachedProducts &	cachedProducts ()
void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
void	prepareForJob (fhicl::ParameterSet const &pset)
void	showMissingConsumes () const

Private Attributes
std::string	fLArG4ModuleLabel
std::string	fGenieModuleLabel
std::string	fVertexModuleLabel
std::string	fEndPoint2dModuleLabel
TH1F *	fRun
TH1F *	fEvt
TH1F *	fTruthVtxXPos
TH1F *	fTruthVtxYPos
TH1F *	fTruthVtxZPos
TH1F *	fTruthWireNumberPlane0
TH1F *	fTruthTimeTickPlane0
TH1F *	fTruthWireNumberPlane1
TH1F *	fTruthTimeTickPlane1
TH1F *	fTruthWireNumberPlane2
TH1F *	fTruthTimeTickPlane2
TH1F *	fTruthWireInCmPlane0
TH1F *	fTruthTimeInCmPlane0
TH1F *	fTruthWireInCmPlane1
TH1F *	fTruthTimeInCmPlane1
TH1F *	fTruthWireInCmPlane2
TH1F *	fTruthTimeInCmPlane2
TH1F *	fTwoDNVtxPlane0
TH1F *	fTwoDNVtxPlane1
TH1F *	fTwoDNVtxPlane2
TH1F *	fTwoDWireNumberPlane0
TH1F *	fTwoDTimeTickPlane0
TH1F *	fTwoDWireNumberPlane1
TH1F *	fTwoDTimeTickPlane1
TH1F *	fTwoDWireNumberPlane2
TH1F *	fTwoDTimeTickPlane2
TH1F *	fTwoDWireInCmPlane0
TH1F *	fTwoDTimeInCmPlane0
TH1F *	fTwoDWireInCmPlane1
TH1F *	fTwoDTimeInCmPlane1
TH1F *	fTwoDWireInCmPlane2
TH1F *	fTwoDTimeInCmPlane2
TH1F *	fTwoDStrengthPlane0
TH1F *	fTwoDStrengthPlane1
TH1F *	fTwoDStrengthPlane2
TH1F *	fRecoCheck2dWireNumPlane0
TH1F *	fRecoCheck2dTimeTickPlane0
TH1F *	fRecoCheck2dWireInCmPlane0
TH1F *	fRecoCheck2dTimeInCmPlane0
TH1F *	fRecoCheck2dWireNumPlane1
TH1F *	fRecoCheck2dTimeTickPlane1
TH1F *	fRecoCheck2dWireInCmPlane1
TH1F *	fRecoCheck2dTimeInCmPlane1
TH1F *	fRecoCheck2dWireNumPlane2
TH1F *	fRecoCheck2dTimeTickPlane2
TH1F *	fRecoCheck2dWireInCmPlane2
TH1F *	fRecoCheck2dTimeInCmPlane2
TH1F *	fRecoVtxN3d
TH1F *	fRecoVtxXPos
TH1F *	fRecoVtxYPos
TH1F *	fRecoVtxZPos
TH1F *	fRecoCheck3dVtxX
TH1F *	fRecoCheck3dVtxY
TH1F *	fRecoCheck3dVtxZ
TH2D *	fRecoCheck3dVtxXvsX
TH2D *	fRecoCheck3dVtxYvsY
TH2D *	fRecoCheck3dVtxZvsZ

Detailed Description

Definition at line 80 of file FeatureVertexFinderAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 39 of file EDAnalyzer.h.

using art::EDAnalyzer::WorkerType = WorkerT<EDAnalyzer>

inherited

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

vertex::FeatureVertexFinderAna::FeatureVertexFinderAna ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 174 of file FeatureVertexFinderAna_module.cc.

References reconfigure().

  : EDAnalyzer(pset)
{
  this->reconfigure(pset);
}

vertex::FeatureVertexFinderAna::~FeatureVertexFinderAna ( )

virtual

Definition at line 181 of file FeatureVertexFinderAna_module.cc.

{
}

Member Function Documentation

void vertex::FeatureVertexFinderAna::analyze

(

const art::Event &

evt

)

Definition at line 278 of file FeatureVertexFinderAna_module.cc.

References detinfo::DetectorProperties::ConvertXToTicks(), DEFINE_ART_MODULE, detinfo::DetectorProperties::DriftVelocity(), detinfo::DetectorProperties::Efield(), art::EventID::event(), fEndPoint2dModuleLabel, fEvt, fGenieModuleLabel, fLArG4ModuleLabel, fRecoCheck2dTimeInCmPlane0, fRecoCheck2dTimeInCmPlane1, fRecoCheck2dTimeInCmPlane2, fRecoCheck2dTimeTickPlane0, fRecoCheck2dTimeTickPlane1, fRecoCheck2dTimeTickPlane2, fRecoCheck2dWireInCmPlane0, fRecoCheck2dWireInCmPlane1, fRecoCheck2dWireInCmPlane2, fRecoCheck2dWireNumPlane0, fRecoCheck2dWireNumPlane1, fRecoCheck2dWireNumPlane2, fRecoCheck3dVtxX, fRecoCheck3dVtxXvsX, fRecoCheck3dVtxY, fRecoCheck3dVtxYvsY, fRecoCheck3dVtxZ, fRecoCheck3dVtxZvsZ, fRecoVtxN3d, fRecoVtxXPos, fRecoVtxYPos, fRecoVtxZPos, fRun, fTruthTimeInCmPlane0, fTruthTimeInCmPlane1, fTruthTimeInCmPlane2, fTruthTimeTickPlane0, fTruthTimeTickPlane1, fTruthTimeTickPlane2, fTruthVtxXPos, fTruthVtxYPos, fTruthVtxZPos, fTruthWireInCmPlane0, fTruthWireInCmPlane1, fTruthWireInCmPlane2, fTruthWireNumberPlane0, fTruthWireNumberPlane1, fTruthWireNumberPlane2, fTwoDNVtxPlane0, fTwoDNVtxPlane1, fTwoDNVtxPlane2, fTwoDStrengthPlane0, fTwoDStrengthPlane1, fTwoDStrengthPlane2, fTwoDTimeInCmPlane0, fTwoDTimeInCmPlane1, fTwoDTimeInCmPlane2, fTwoDTimeTickPlane0, fTwoDTimeTickPlane1, fTwoDTimeTickPlane2, fTwoDWireInCmPlane0, fTwoDWireInCmPlane1, fTwoDWireInCmPlane2, fTwoDWireNumberPlane0, fTwoDWireNumberPlane1, fTwoDWireNumberPlane2, fVertexModuleLabel, art::DataViewImpl::getByLabel(), simb::MCTruth::GetParticle(), detinfo::DetectorProperties::GetXTicksOffset(), art::Event::id(), geo::GeometryCore::IteratePlaneIDs(), geo::GeometryCore::NearestWire(), art::PtrVector< T >::push_back(), art::Event::run(), detinfo::DetectorProperties::SamplingRate(), art::PtrVector< T >::size(), detinfo::DetectorProperties::Temperature(), geo::GeometryCore::Views(), simb::MCParticle::Vx(), and geo::GeometryCore::WirePitch().

{
  
  // ###############################
  // ### Filling Run/Event Histo ###
  // ###############################
  fRun->Fill( evt.run() );
  fEvt->Fill( evt.id().event() );
     
  // ##############################################
  // ### Getting the Geant Information Directly ###
  // ##############################################
  art::Handle< std::vector<simb::MCParticle> > mcParticleHandle;
  evt.getByLabel(fLArG4ModuleLabel, mcParticleHandle);

  // #######################################
  // ### Getting MC Truth Info from simb ###
  // #######################################
  art::Handle< std::vector<simb::MCTruth> > mctruthListHandle;
  evt.getByLabel(fGenieModuleLabel,mctruthListHandle); 
  
  // ############################################
  // ### Getting information from BackTrackerService ###
  // ############################################
  art::ServiceHandle<cheat::BackTrackerService> bt_serv;
  
  // ####################################
  // ### Getting Geometry Information ###
  // ####################################
  art::ServiceHandle<geo::Geometry> geom;
  
  // ###################################
  // ### Getting Detector Properties ###
  // ###################################
  const detinfo::DetectorProperties* detp = lar::providerFrom<detinfo::DetectorPropertiesService>();
  
  // ######################################################
  // ### Getting 2d Vertex information (vertex2dHandle) ###
  // ######################################################
  art::Handle< std::vector<recob::EndPoint2D> > vertex2dHandle;
  evt.getByLabel(fEndPoint2dModuleLabel,vertex2dHandle);
  
  // ##################################################
  // ### Getting the 3d Vertex (vertex3dListHandle) ###
  // ##################################################
  art::Handle< std::vector<recob::Vertex> > vertex3dListHandle;
  evt.getByLabel(fVertexModuleLabel,vertex3dListHandle); 
  
  
  
  // ==========================================================
  // === Detector numbers that are useful (to be set later) ===
  // ==========================================================
  //size_t nplanes = 0;
  std::vector<double> WirePitch_CurrentPlane(geom->Views().size(), 0.); //<---Setting the Wire pitch for each plane
                                           // Right now assuming only 3 planes
  // ##############################
  // ### get the wire pitch for each view ###
  // ##############################
  size_t vn = 0;
  for(auto v : geom->Views() ){
    WirePitch_CurrentPlane[vn] = geom->WirePitch(v);
    ++vn;
  }

  // ##################################################
  // ### Calculating the Timetick to CM conversion  ###
  // ##################################################
  float TimeTick      = detp->SamplingRate()/1000.; //<---To get units of microsecond...not nanosec
  float DriftVelocity = detp->DriftVelocity(detp->Efield(),detp->Temperature());
                
  float TimetoCm      = TimeTick*DriftVelocity;




// ===========================================================================================================
// ==================================== Looping over MC information ==========================================
// ===========================================================================================================
 
 
   // ##################################
   // ### Getting MC Truth simb Info ###
   // ##################################
   art::PtrVector<simb::MCTruth> mclist;
   for (unsigned int ii = 0; ii <  mctruthListHandle->size(); ++ii)
        {
        art::Ptr<simb::MCTruth> mctparticle(mctruthListHandle,ii);
        mclist.push_back(mctparticle);
        }
   
   // ##############################################
   // ### Variables for truth vertex information ###
   // ##############################################
   float    truth_vertex[5] = {0.};             //<---Truth x,y,z information
   
   uint32_t VtxWireNum[3] = {0};                //<---Wire number in each plane ( WireNum[plane#] )
   double   VtxTimeTick[3] = {0.};              //<---Time tick in each plane   ( TimeTick[plane#] )
   
   double   VtxWireNum_InCM[3] = {0.};          //<---Wire number in each plane in CM ( WireNum[plane#] )
   double   VtxTimeTick_InCM[3] = {0.};         //<---Time tick in each plane in CM   ( TimeTick[plane#] )
   
   
   // ###################################
   // ### Finding the MC truth vertex ###
   // ###################################
   for( unsigned int i = 0; i < mclist.size(); ++i )
        {
        art::Ptr<simb::MCTruth> mc(mclist[i]);
        simb::MCParticle neut(mc->GetParticle(i)); 
        
        // ############################################
        // ### Filling the vertex x,y,z information ###
        // ############################################
        truth_vertex[0] =neut.Vx();
        truth_vertex[1] =neut.Vy();
        truth_vertex[2] =neut.Vz();
        
        }// end i loop
        
  // ### Filling Histograms ###
  fTruthVtxXPos->Fill( truth_vertex[0] );
  fTruthVtxYPos->Fill( truth_vertex[1] );
  fTruthVtxZPos->Fill( truth_vertex[2] );
  
  // ##############################
  // ### Looping over geo::PlaneIDs ###
  // ##############################
  for(auto const& pid : geom->IteratePlaneIDs()){
    // ############################################################################
    // ### Calculating the nearest wire the vertex corresponds to in each plane ###
    // ############################################################################
    
    try{
      //                  geom->NearestWire(worldLoc[3], Plane#, TPC#, Cyrostat#)
        VtxWireNum[pid.Plane] = geom->NearestWire(truth_vertex,pid.Plane,pid.TPC,pid.Cryostat);
    }
    catch(...){
      mf::LogWarning("FeatureVertexFinderAna") << "Can't find nearest wire";
      continue;
    }
    //             detp->ConvertXToTicks(xpos, plane#, TPC#, Cyrostat#)
    VtxTimeTick[pid.Plane] = (detp->ConvertXToTicks(truth_vertex[0],pid.Plane, pid.TPC, pid.Cryostat)
                          +  detp->GetXTicksOffset(pid.Plane, pid.TPC, pid.Cryostat));
    
    // ======================== Translating each of these in cm =====================
    VtxWireNum_InCM[pid.Plane]  = VtxWireNum[pid.Plane] * WirePitch_CurrentPlane[pid.Plane];
    VtxTimeTick_InCM[pid.Plane] = VtxTimeTick[pid.Plane] * TimetoCm;
  }//<---End pid loop
   
   // ### Filling Histograms ###
   fTruthWireNumberPlane0->Fill( VtxWireNum[0]  );
   fTruthTimeTickPlane0->Fill( VtxTimeTick[0] );
   fTruthWireNumberPlane1->Fill( VtxWireNum[1]  );
   fTruthTimeTickPlane1->Fill( VtxTimeTick[1] );
   fTruthWireNumberPlane2->Fill( VtxWireNum[2]  );
   fTruthTimeTickPlane2->Fill( VtxTimeTick[2] );
   
   fTruthWireInCmPlane0->Fill( VtxWireNum_InCM[0] );
   fTruthTimeInCmPlane0->Fill( VtxTimeTick_InCM[0] );
   fTruthWireInCmPlane1->Fill( VtxWireNum_InCM[1] );
   fTruthTimeInCmPlane1->Fill( VtxTimeTick_InCM[1] );
   fTruthWireInCmPlane2->Fill( VtxWireNum_InCM[2] );
   fTruthTimeInCmPlane2->Fill( VtxTimeTick_InCM[2] );


// ===================================================================================================================
// ==================================== Looping over EndPoint2d information ==========================================
// ===================================================================================================================

   art::PtrVector<recob::EndPoint2D> vert2d;
   
   // Variables for Vertex2d
   
   double Vertex2d_TimeTick[10000] = {0.};      //<---Vertex2d Time Tick for the current plane ( TimeTick[#2d] )
   double Vertex2d_Wire[10000]   = {0.};        //<---Veretx2d Wire # ( Wire[#2d] )
   
   double Vertex2d_TimeTick_InCM[10000] = {0.}; //<---Vertex 2d Time tick in CM ( TimeTick[#2d] )
   double Vertex2d_Wire_InCM[10000]     = {0.}; //<---Veretx2d Wire in CM ( Wire[#2d] )
   int n2dVtx = 0;
   
   int n2dVtxPlane0 = 0, n2dVtxPlane1 = 0, n2dVtxPlane2 = 0;
   
   bool vertexWstrengthplane0 = false , vertexWstrengthplane1 = false;//, vertexWstrengthplane2 = false;   //commented out, Wes, 12.4.13
   
   // ######################################
   // ### Loop over the EndPoint2d List  ###
   // ######################################
   for(size_t ii = 0; ii < vertex2dHandle->size(); ++ii){
     art::Ptr<recob::EndPoint2D> vertex(vertex2dHandle, ii);
     vert2d.push_back(vertex);
   }
   
  // ############################################ 
  // ### If we have 2d vertex, loop over them ###
  // ############################################
  if(vert2d.size() > 0)
        {

        // ##############################
        // ### Looping over geo::PlaneIDs ###
        // ##############################
          for(auto const& pid : geom->IteratePlaneIDs() ){
            for(size_t ww = 0; ww<vert2d.size(); ++ww){
              //std::cout<<"plane = "<<plane<<std::endl;
              //std::cout<<"vert2d[ww]->WireID().Plane = "<<vert2d[ww]->WireID().Plane<<std::endl;
              // Only look at this 2d vertex if it is in the current plane
              if (vert2d[ww]->WireID().planeID() != pid){continue;}
                                        
              Vertex2d_TimeTick[n2dVtx] = vert2d[ww]->DriftTime();
              Vertex2d_Wire[n2dVtx]     = vert2d[ww]->WireID().Wire;
              
              // ======================== Translating each of these in cm =====================
              Vertex2d_Wire_InCM[n2dVtx]           = Vertex2d_Wire[n2dVtx] * WirePitch_CurrentPlane[pid.Plane]; 
              //std::cout<<"Vertex2d_Wire[n2dVtx] = "<<Vertex2d_Wire[n2dVtx]<<" , WirePitch_CurrentPlane[plane] = "<<WirePitch_CurrentPlane[plane]<<std::endl;
              //std::cout<<"Vertex2d_Wire_InCM[n2dVtx] = "<<Vertex2d_Wire_InCM[n2dVtx]<<std::endl;
              Vertex2d_TimeTick_InCM[n2dVtx]     = Vertex2d_TimeTick[n2dVtx] * TimetoCm;
                                        
              // ###########################################################################
              // ### Checking how well we did in reconstructing the vertex (Reco - True) ###
              // ###########################################################################
              
              float RecoCheck_TimeTick = Vertex2d_TimeTick[n2dVtx] - VtxTimeTick[pid.Plane];
              float RecoCheck_WireNum  = Vertex2d_Wire[n2dVtx] - VtxWireNum[pid.Plane];
              
              float RecoCheck_TimeInCm = Vertex2d_TimeTick_InCM[n2dVtx] - VtxTimeTick_InCM[pid.Plane];
              float RecoCheck_WireInCm = Vertex2d_Wire_InCM[n2dVtx] - VtxWireNum_InCM[pid.Plane];
              
                                        
              
              if( vert2d[ww]->Strength() > -1)
                {
                  if(pid.Plane == 0)
                    {
                      vertexWstrengthplane0 = true;
                      //std::cout<<"Filling Plane 0 "<<std::endl;
                      //std::cout<<"Vertex2d_Wire[n2dVtx] = "<<Vertex2d_Wire[n2dVtx]<<std::endl;
                      
                      fTwoDWireNumberPlane0->Fill( Vertex2d_Wire[n2dVtx]  );
                      fTwoDTimeTickPlane0->Fill( Vertex2d_TimeTick[n2dVtx] );
                      fTwoDWireInCmPlane0->Fill( Vertex2d_Wire_InCM[n2dVtx] );
                      fTwoDTimeInCmPlane0->Fill( Vertex2d_TimeTick_InCM[n2dVtx] );
                      
                      fTwoDStrengthPlane0->Fill( vert2d[ww]->Strength() );
                      
                      fRecoCheck2dWireNumPlane0->Fill( RecoCheck_WireNum );
                      fRecoCheck2dTimeTickPlane0->Fill( RecoCheck_TimeTick );
                      fRecoCheck2dWireInCmPlane0->Fill( RecoCheck_WireInCm );
                      fRecoCheck2dTimeInCmPlane0->Fill( RecoCheck_TimeInCm );
                      
                      n2dVtxPlane0++;
                      
                    }//<---End Plane 0
                                
                  if(pid.Plane == 1)
                    {
                      vertexWstrengthplane1 = true;
                      fTwoDWireNumberPlane1->Fill( Vertex2d_Wire[n2dVtx]  );
                      fTwoDTimeTickPlane1->Fill( Vertex2d_TimeTick[n2dVtx] );
                      fTwoDWireInCmPlane1->Fill( Vertex2d_Wire_InCM[n2dVtx] );
                      fTwoDTimeInCmPlane1->Fill( Vertex2d_TimeTick_InCM[n2dVtx] );
                      
                      fTwoDStrengthPlane1->Fill( vert2d[ww]->Strength() );
                      
                      fRecoCheck2dWireNumPlane1->Fill( RecoCheck_WireNum );
                      fRecoCheck2dTimeTickPlane1->Fill( RecoCheck_TimeTick );
                      fRecoCheck2dWireInCmPlane1->Fill( RecoCheck_WireInCm );
                      fRecoCheck2dTimeInCmPlane1->Fill( RecoCheck_TimeInCm );
                      
                      n2dVtxPlane1++;
                      
                    }//<---End Plane 1
                  
                  if(pid.Plane == 2)
                    {
                      //vertexWstrengthplane2 = true;
                      fTwoDWireNumberPlane2->Fill( Vertex2d_Wire[n2dVtx]  );
                      fTwoDTimeTickPlane2->Fill( Vertex2d_TimeTick[n2dVtx] );
                      fTwoDWireInCmPlane2->Fill( Vertex2d_Wire_InCM[n2dVtx] );
                      fTwoDTimeInCmPlane2->Fill( Vertex2d_TimeTick_InCM[n2dVtx] );
                      
                      fTwoDStrengthPlane2->Fill( vert2d[ww]->Strength() );
                      
                      fRecoCheck2dWireNumPlane2->Fill( RecoCheck_WireNum );
                      fRecoCheck2dTimeTickPlane2->Fill( RecoCheck_TimeTick );
                      fRecoCheck2dWireInCmPlane2->Fill( RecoCheck_WireInCm );
                      fRecoCheck2dTimeInCmPlane2->Fill( RecoCheck_TimeInCm );
                      
                      n2dVtxPlane2++;
                      
                      
                    }//<---End Plane 2
                }
              
              
              ++n2dVtx;
            }//<--end ww loop
            
            
          }//<---End plane ID

        }//<--End checking if we have 2d end points
        
  fTwoDNVtxPlane0->Fill( n2dVtxPlane0 );
  fTwoDNVtxPlane1->Fill( n2dVtxPlane1 );
  fTwoDNVtxPlane2->Fill( n2dVtxPlane2 );
  
  // =================================================================================================================
// ==================================== Looping over 3dVertex information ==========================================
// =================================================================================================================
   art::PtrVector<recob::Vertex>     Vertexlist;
   
   double xyz[3] = {0.};

   for(unsigned int ii = 0; ii < vertex3dListHandle->size(); ++ii)
        {
        art::Ptr<recob::Vertex> vertex3d(vertex3dListHandle, ii);
        Vertexlist.push_back(vertex3d); // class member
        
        }//<---End ii loop
        
   if(Vertexlist.size() > 0 && vertexWstrengthplane0 && vertexWstrengthplane1)
        {
        
        fRecoVtxN3d->Fill(Vertexlist.size());
        for(unsigned int ww = 0; ww<Vertexlist.size(); ww++)
                {
                Vertexlist[ww]->XYZ(xyz);
                
                // ###########################################################################
                // ### Checking how well we did in reconstructing the vertex (Reco - True) ###
                // ###########################################################################
   
                // === Finding the Delta X, Y, Z between Reco vtx and truth ===
        
                double DeltaX = xyz[0] - truth_vertex[0] ;
                double DeltaY = xyz[1] - truth_vertex[1] ;
                double DeltaZ = xyz[2] - truth_vertex[2] ;
        
                double DeltaXoverTrueX = DeltaX / truth_vertex[0];
                double DeltaYoverTrueY = DeltaY / truth_vertex[0];
                double DeltaZoverTrueZ = DeltaZ / truth_vertex[0];
   
                fRecoVtxXPos->Fill( xyz[0] );
                fRecoVtxYPos->Fill( xyz[1] );
                fRecoVtxZPos->Fill( xyz[2] );
                
                fRecoCheck3dVtxX->Fill( DeltaX );
                fRecoCheck3dVtxY->Fill( DeltaY );
                fRecoCheck3dVtxZ->Fill( DeltaZ );
   
                fRecoCheck3dVtxXvsX->Fill( xyz[0], DeltaXoverTrueX );
                fRecoCheck3dVtxYvsY->Fill( xyz[1], DeltaYoverTrueY );
                fRecoCheck3dVtxZvsZ->Fill( xyz[2], DeltaZoverTrueZ );
                
                
                
                
                }
        }
   
return;  
}//<---End access the event

void vertex::FeatureVertexFinderAna::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 197 of file FeatureVertexFinderAna_module.cc.

References fEvt, fRecoCheck2dTimeInCmPlane0, fRecoCheck2dTimeInCmPlane1, fRecoCheck2dTimeInCmPlane2, fRecoCheck2dTimeTickPlane0, fRecoCheck2dTimeTickPlane1, fRecoCheck2dTimeTickPlane2, fRecoCheck2dWireInCmPlane0, fRecoCheck2dWireInCmPlane1, fRecoCheck2dWireInCmPlane2, fRecoCheck2dWireNumPlane0, fRecoCheck2dWireNumPlane1, fRecoCheck2dWireNumPlane2, fRecoCheck3dVtxX, fRecoCheck3dVtxXvsX, fRecoCheck3dVtxY, fRecoCheck3dVtxYvsY, fRecoCheck3dVtxZ, fRecoCheck3dVtxZvsZ, fRecoVtxN3d, fRecoVtxXPos, fRecoVtxYPos, fRecoVtxZPos, fRun, fTruthTimeInCmPlane0, fTruthTimeInCmPlane1, fTruthTimeInCmPlane2, fTruthTimeTickPlane0, fTruthTimeTickPlane1, fTruthTimeTickPlane2, fTruthVtxXPos, fTruthVtxYPos, fTruthVtxZPos, fTruthWireInCmPlane0, fTruthWireInCmPlane1, fTruthWireInCmPlane2, fTruthWireNumberPlane0, fTruthWireNumberPlane1, fTruthWireNumberPlane2, fTwoDNVtxPlane0, fTwoDNVtxPlane1, fTwoDNVtxPlane2, fTwoDStrengthPlane0, fTwoDStrengthPlane1, fTwoDStrengthPlane2, fTwoDTimeInCmPlane0, fTwoDTimeInCmPlane1, fTwoDTimeInCmPlane2, fTwoDTimeTickPlane0, fTwoDTimeTickPlane1, fTwoDTimeTickPlane2, fTwoDWireInCmPlane0, fTwoDWireInCmPlane1, fTwoDWireInCmPlane2, fTwoDWireNumberPlane0, fTwoDWireNumberPlane1, fTwoDWireNumberPlane2, and art::TFileDirectory::make().

{
        // get access to the TFile service
        art::ServiceHandle<art::TFileService> tfs;
        // ====================================
        // ==== Outputting TH1F Histograms ====
        // ====================================
        fRun                    = tfs->make<TH1F>("fRun", "Run Number", 1000, 0, 1000);
        fEvt                    = tfs->make<TH1F>("fEvt", "Event Number", 1000, 0, 1000);
        fTruthVtxXPos           = tfs->make<TH1F>("fTruthVtxXPos", "Truth Vertex X Position", 400, 0, 200);
        fTruthVtxYPos           = tfs->make<TH1F>("fTruthVtxYPos", "Truth Vertex Y Position", 400, -100, 100);
        fTruthVtxZPos           = tfs->make<TH1F>("fTruthVtxZPos", "Truth Vertex Z Position", 2000, 0, 1000);
        fTruthWireNumberPlane0  = tfs->make<TH1F>("fTruthWireNumberPlane0", "Truth Wire Number Plane 0", 3000, 0, 3000);
        fTruthTimeTickPlane0    = tfs->make<TH1F>("fTruthTimeTickPlane0", "Truth Time Tick Plane 0", 3200, 0, 3200);
        fTruthWireNumberPlane1  = tfs->make<TH1F>("fTruthWireNumberPlane1", "Truth Wire Number Plane 1", 3000, 0, 3000);
        fTruthTimeTickPlane1    = tfs->make<TH1F>("fTruthTimeTickPlane1", "Truth Time Tick Plane 1", 3200, 0, 3200);
        fTruthWireNumberPlane2  = tfs->make<TH1F>("fTruthWireNumberPlane2", "Truth Wire Number Plane 2", 3000, 0, 3000);
        fTruthTimeTickPlane2    = tfs->make<TH1F>("fTruthTimeTickPlane2", "Truth Time Tick Plane 2", 3200, 0, 3200);
        fTruthWireInCmPlane0    = tfs->make<TH1F>("fTruthWireInCmPlane0", "Truth Wire In CM Plane 0", 2000, 0, 1000);
        fTruthTimeInCmPlane0    = tfs->make<TH1F>("fTruthTimeInCmPlane0", "Truth Time In Cm Plane 0", 2000, 0, 1000);
        fTruthWireInCmPlane1    = tfs->make<TH1F>("fTruthWireInCmPlane1", "Truth Wire In CM Plane 1", 2000, 0, 1000);
        fTruthTimeInCmPlane1    = tfs->make<TH1F>("fTruthTimeInCmPlane1", "Truth Time In Cm Plane 1", 2000, 0, 1000);
        fTruthWireInCmPlane2    = tfs->make<TH1F>("fTruthWireInCmPlane2", "Truth Wire In CM Plane 2", 2000, 0, 1000);
        fTruthTimeInCmPlane2    = tfs->make<TH1F>("fTruthTimeInCmPlane2", "Truth Time In Cm Plane 2", 2000, 0, 1000);
        
        fTwoDNVtxPlane0         = tfs->make<TH1F>("fTwoDNVtxPlane0", "TwoD Number of Verticies Found in Plane 0", 400, 0, 200);
        fTwoDNVtxPlane1         = tfs->make<TH1F>("fTwoDNVtxPlane1", "TwoD Number of Verticies Found in Plane 1", 400, 0, 200);
        fTwoDNVtxPlane2         = tfs->make<TH1F>("fTwoDNVtxPlane2", "TwoD Number of Verticies Found in Plane 2", 400, 0, 200);
        
        fTwoDWireNumberPlane0   = tfs->make<TH1F>("fTwoDWireNumberPlane0", "TwoD Wire Number Plane 0", 3000, 0, 3000);
        fTwoDTimeTickPlane0     = tfs->make<TH1F>("fTwoDTimeTickPlane0", "TwoD Time Tick Plane 0", 3200, 0, 3200);
        fTwoDWireNumberPlane1   = tfs->make<TH1F>("fTwoDWireNumberPlane1", "TwoD Wire Number Plane 1", 3000, 0, 3000);
        fTwoDTimeTickPlane1     = tfs->make<TH1F>("fTwoDTimeTickPlane1", "TwoD Time Tick Plane 1", 3200, 0, 3200);
        fTwoDWireNumberPlane2   = tfs->make<TH1F>("fTwoDWireNumberPlane2", "TwoD Wire Number Plane 2", 3000, 0, 3000);
        fTwoDTimeTickPlane2     = tfs->make<TH1F>("fTwoDTimeTickPlane2", "TwoD Time Tick Plane 2", 3200, 0, 3200);
        fTwoDWireInCmPlane0     = tfs->make<TH1F>("fTwoDWireInCmPlane0", "TwoD Wire In CM Plane 0", 2000, 0, 1000);
        fTwoDTimeInCmPlane0     = tfs->make<TH1F>("fTwoDTimeInCmPlane0", "TwoD Time In Cm Plane 0", 2000, 0, 1000);
        fTwoDWireInCmPlane1     = tfs->make<TH1F>("fTwoDWireInCmPlane1", "TwoD Wire In CM Plane 1", 2000, 0, 1000);
        fTwoDTimeInCmPlane1     = tfs->make<TH1F>("fTwoDTimeInCmPlane1", "TwoD Time In Cm Plane 1", 2000, 0, 1000);
        fTwoDWireInCmPlane2     = tfs->make<TH1F>("fTwoDWireInCmPlane2", "TwoD Wire In CM Plane 2", 2000, 0, 1000);
        fTwoDTimeInCmPlane2     = tfs->make<TH1F>("fTwoDTimeInCmPlane2", "TwoD Time In Cm Plane 2", 2000, 0, 1000);
        
        fTwoDStrengthPlane0     = tfs->make<TH1F>("fTwoDStrengthPlane0", "TwoD Strength Plane 0", 1000, 0, 500);
        fTwoDStrengthPlane1     = tfs->make<TH1F>("fTwoDStrengthPlane1", "TwoD Strength Plane 1", 1000, 0, 500);
        fTwoDStrengthPlane2     = tfs->make<TH1F>("fTwoDStrengthPlane2", "TwoD Strength Plane 2", 1000, 0, 500);
        
        fRecoCheck2dWireNumPlane0       = tfs->make<TH1F>("fRecoCheck2dWireNumPlane0", "Reco Wire Number - True Wire Number Plane 0", 400, -200, 200);
        fRecoCheck2dTimeTickPlane0      = tfs->make<TH1F>("fRecoCheck2dTimeTickPlane0", "Reco Time Tick - True Time Tick Plane 0", 1000, -500, 500);
        fRecoCheck2dWireInCmPlane0      = tfs->make<TH1F>("fRecoCheck2dWireInCmPlane0", "Reco Wire in CM - True Wire in CM Plane 0", 200, -50, 50);     
        fRecoCheck2dTimeInCmPlane0      = tfs->make<TH1F>("fRecoCheck2dTimeInCmPlane0", "Reco Time in CM - True Time in CM Plane 0", 200, -50, 50);
        
        fRecoCheck2dWireNumPlane1       = tfs->make<TH1F>("fRecoCheck2dWireNumPlane1", "Reco Wire Number - True Wire Number Plane 1", 400, -200, 200);
        fRecoCheck2dTimeTickPlane1      = tfs->make<TH1F>("fRecoCheck2dTimeTickPlane1", "Reco Time Tick - True Time Tick Plane 1", 1000, -500, 500);
        fRecoCheck2dWireInCmPlane1      = tfs->make<TH1F>("fRecoCheck2dWireInCmPlane1", "Reco Wire in CM - True Wire in CM Plane 1", 200, -50, 50);     
        fRecoCheck2dTimeInCmPlane1      = tfs->make<TH1F>("fRecoCheck2dTimeInCmPlane1", "Reco Time in CM - True Time in CM Plane 1", 200, -50, 50);
        
        fRecoCheck2dWireNumPlane2       = tfs->make<TH1F>("fRecoCheck2dWireNumPlane2", "Reco Wire Number - True Wire Number Plane 2", 400, -200, 200);
        fRecoCheck2dTimeTickPlane2      = tfs->make<TH1F>("fRecoCheck2dTimeTickPlane2", "Reco Time Tick - True Time Tick Plane 2", 1000, -500, 500);
        fRecoCheck2dWireInCmPlane2      = tfs->make<TH1F>("fRecoCheck2dWireInCmPlane2", "Reco Wire in CM - True Wire in CM Plane 2", 200, -50, 50);     
        fRecoCheck2dTimeInCmPlane2      = tfs->make<TH1F>("fRecoCheck2dTimeInCmPlane2", "Reco Time in CM - True Time in CM Plane 2", 200, -50, 50);
        
        fRecoVtxN3d             = tfs->make<TH1F>("fRecoVtxN3d", "Number of 3d-Reco Verticies", 400, 0, 200);
        fRecoVtxXPos            = tfs->make<TH1F>("fRecoVtxXPos", "Reco Vertex X Position", 400, -10, 200);
        fRecoVtxYPos            = tfs->make<TH1F>("fRecoVtxYPos", "Reco Vertex Y Position", 400, -100, 100);
        fRecoVtxZPos            = tfs->make<TH1F>("fRecoVtxZPos", "Reco Vertex Z Position", 2000, -10, 1000);
        
        fRecoCheck3dVtxX        = tfs->make<TH1F>("fRecoCheck3dVtxX", "Reco X Position - True X Postion", 800, -100, 100);
        fRecoCheck3dVtxY        = tfs->make<TH1F>("fRecoCheck3dVtxY", "Reco Y Position - True Y Postion", 800, -100, 100);
        fRecoCheck3dVtxZ        = tfs->make<TH1F>("fRecoCheck3dVtxZ", "Reco Z Position - True Z Postion", 800, -100, 100);
        
        fRecoCheck3dVtxXvsX     = tfs->make <TH2D>("fRecoCheck3dVtxXvsX", "(Reco X - True X)/True X vs True X", 400, -10, 200, 120, -20, 20);
        fRecoCheck3dVtxYvsY     = tfs->make <TH2D>("fRecoCheck3dVtxYvsY", "(Reco Y - True Y)/True Y vs True Y", 400, -100, 100, 120, -20, 20);
        fRecoCheck3dVtxZvsZ     = tfs->make <TH2D>("fRecoCheck3dVtxZvsZ", "(Reco Z - True Z)/True Z vs True Z", 1000, -10, 1000, 120, -20, 20);
        
        
        return;
}

detail::CachedProducts& art::EventObserverBase::cachedProducts ( )

inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

Referenced by art::EDAnalyzer::doEvent(), and art::OutputModule::doWriteEvent().

  {
    return selectors_;
  }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

{
  return rng()->createEngine(placeholder_schedule_id(), seed);
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make  )

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make);
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make, label_t const &  engine_label  )

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
}

CurrentProcessingContext const * art::EDAnalyzer::currentContext ( ) const

protectedinherited

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

Referenced by art::EDAnalyzer::workerType().

  {
    return current_context_.get();
  }

EngineCreator::seed_t EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset, char const  key[] = "seed", seed_t const  implicit_seed = -1  )

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

{
  auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
  return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
}

art::Handle<art::TriggerResults> art::EventObserverBase::getTriggerResults ( Event const &  e ) const

inlineinherited

Definition at line 61 of file EventObserverBase.h.

References art::detail::CachedProducts::getOneTriggerResults(), and art::EventObserverBase::selectors_.

Referenced by art::OutputModule::doWriteEvent().

  {
    return selectors_.getOneTriggerResults(e);
  }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it )

inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it )

inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

bool art::EventObserverBase::modifiesEvent ( ) const

inlineinherited

Definition at line 25 of file EventObserverBase.h.

  {
    return false;
  }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

Referenced by art::RootOutput::beginSubRun(), art::OutputModule::doBeginRun(), art::OutputModule::doBeginSubRun(), art::OutputModule::doEndRun(), art::OutputModule::doEndSubRun(), art::ProducingService::doPostReadEvent(), art::ProducingService::doPostReadRun(), art::ProducingService::doPostReadSubRun(), art::OutputModule::doWriteEvent(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EventProcessor::readEvent(), art::EventProcessor::readRun(), art::EmptyEvent::readRun_(), art::EventProcessor::readSubRun(), and art::EmptyEvent::readSubRun_().

{
  static Consumer invalid{InvalidTag{}};
  return &invalid;
}

void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

{
  if (!moduleContext_)
    return;

  pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
  for (auto& consumablesPerBranch : consumables_) {
    cet::sort_all(consumablesPerBranch);
  }
}

std::string const& art::EventObserverBase::processName ( ) const

inlineinherited

Definition at line 41 of file EventObserverBase.h.

References art::EventObserverBase::process_name_.

Referenced by art::FileDumperOutput::printPrincipal(), and art::RootOutput::RootOutput().

  {
    return process_name_;
  }

void vertex::FeatureVertexFinderAna::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 187 of file FeatureVertexFinderAna_module.cc.

References fEndPoint2dModuleLabel, fGenieModuleLabel, fLArG4ModuleLabel, fVertexModuleLabel, and fhicl::ParameterSet::get().

Referenced by FeatureVertexFinderAna().

{
        fLArG4ModuleLabel        = p.get< std::string >("LArGeantModuleLabel");
        fGenieModuleLabel        = p.get< std::string >("GenieModuleLabel");
        fVertexModuleLabel       = p.get< std::string >("VertexModuleLabel");
        fEndPoint2dModuleLabel   = p.get< std::string >("EndPoint2dModuleLabel");
        return;
}

void art::EventObserverBase::registerProducts ( MasterProductRegistry &  , ProductDescriptions &  , ModuleDescription const &    )

inlineinherited

Definition at line 33 of file EventObserverBase.h.

  {}

fhicl::ParameterSetID art::EventObserverBase::selectorConfig ( ) const

inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

Referenced by art::RootOutputFile::writeOne().

  {
    return selector_config_id_;
  }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

{
  if (!moduleContext_)
    return;

  // If none of the branches have missing consumes statements, exit early.
  if (std::all_of(cbegin(missingConsumes_),
                  cend(missingConsumes_),
                  [](auto const& perBranch) { return perBranch.empty(); }))
    return;

  constexpr cet::HorizontalRule rule{60};
  mf::LogPrint log{"MTdiagnostics"};
  log << '\n'
      << rule('=') << '\n'
      << "The following consumes (or mayConsume) statements are missing from\n"
      << module_context(moduleDescription_) << '\n'
      << rule('-') << '\n';

  cet::for_all_with_index(
    missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
      for (auto const& pi : perBranch) {
        log << "  "
            << assemble_consumes_statement(static_cast<BranchType>(i), pi)
            << '\n';
      }
    });
  log << rule('=');
}

void art::Consumer::validateConsumedProduct ( BranchType const  bt, ProductInfo const &  pi  )

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

{
  // Early exits if consumes tracking has been disabled or if the
  // consumed product is an allowed consumable.
  if (!moduleContext_)
    return;

  if (cet::binary_search_all(consumables_[bt], pi))
    return;

  if (requireConsumes_) {
    throw Exception(errors::ProductRegistrationFailure,
                    "Consumer: an error occurred during validation of a "
                    "retrieved product\n\n")
      << "The following consumes (or mayConsume) statement is missing from\n"
      << module_context(moduleDescription_) << ":\n\n"
      << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
  }

  missingConsumes_[bt].insert(pi);
}

bool art::EventObserverBase::wantAllEvents ( ) const

inlineinherited

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), art::RootOutput::RootOutput(), and art::OutputWorker::wantAllEvents().

  {
    return wantAllEvents_;
  }

bool art::EventObserverBase::wantEvent ( Event const &  e )

inlineinherited

Definition at line 51 of file EventObserverBase.h.

References art::EventObserverBase::selectors_, and art::detail::CachedProducts::wantEvent().

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), and art::OutputModule::doWriteEvent().

  {
    return selectors_.wantEvent(e);
  }

std::string art::EDAnalyzer::workerType ( ) const

inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

    {
      return "WorkerT<EDAnalyzer>";
    }

Member Data Documentation

std::string vertex::FeatureVertexFinderAna::fEndPoint2dModuleLabel

private

Definition at line 96 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and reconfigure().

TH1F* vertex::FeatureVertexFinderAna::fEvt

private

Definition at line 102 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

std::string vertex::FeatureVertexFinderAna::fGenieModuleLabel

private

Definition at line 94 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and reconfigure().

std::string vertex::FeatureVertexFinderAna::fLArG4ModuleLabel

private

Definition at line 93 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and reconfigure().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dTimeInCmPlane0

private

Definition at line 143 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dTimeInCmPlane1

private

Definition at line 148 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dTimeInCmPlane2

private

Definition at line 153 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dTimeTickPlane0

private

Definition at line 141 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dTimeTickPlane1

private

Definition at line 146 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dTimeTickPlane2

private

Definition at line 151 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dWireInCmPlane0

private

Definition at line 142 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dWireInCmPlane1

private

Definition at line 147 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dWireInCmPlane2

private

Definition at line 152 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dWireNumPlane0

private

Definition at line 140 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dWireNumPlane1

private

Definition at line 145 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck2dWireNumPlane2

private

Definition at line 150 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck3dVtxX

private

Definition at line 160 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH2D* vertex::FeatureVertexFinderAna::fRecoCheck3dVtxXvsX

private

Definition at line 164 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck3dVtxY

private

Definition at line 161 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH2D* vertex::FeatureVertexFinderAna::fRecoCheck3dVtxYvsY

private

Definition at line 165 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoCheck3dVtxZ

private

Definition at line 162 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH2D* vertex::FeatureVertexFinderAna::fRecoCheck3dVtxZvsZ

private

Definition at line 166 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoVtxN3d

private

Definition at line 155 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoVtxXPos

private

Definition at line 156 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoVtxYPos

private

Definition at line 157 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRecoVtxZPos

private

Definition at line 158 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fRun

private

Definition at line 101 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthTimeInCmPlane0

private

Definition at line 113 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthTimeInCmPlane1

private

Definition at line 115 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthTimeInCmPlane2

private

Definition at line 117 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthTimeTickPlane0

private

Definition at line 107 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthTimeTickPlane1

private

Definition at line 109 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthTimeTickPlane2

private

Definition at line 111 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthVtxXPos

private

Definition at line 103 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthVtxYPos

private

Definition at line 104 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthVtxZPos

private

Definition at line 105 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthWireInCmPlane0

private

Definition at line 112 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthWireInCmPlane1

private

Definition at line 114 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthWireInCmPlane2

private

Definition at line 116 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthWireNumberPlane0

private

Definition at line 106 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthWireNumberPlane1

private

Definition at line 108 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTruthWireNumberPlane2

private

Definition at line 110 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDNVtxPlane0

private

Definition at line 119 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDNVtxPlane1

private

Definition at line 120 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDNVtxPlane2

private

Definition at line 121 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDStrengthPlane0

private

Definition at line 136 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDStrengthPlane1

private

Definition at line 137 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDStrengthPlane2

private

Definition at line 138 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDTimeInCmPlane0

private

Definition at line 130 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDTimeInCmPlane1

private

Definition at line 132 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDTimeInCmPlane2

private

Definition at line 134 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDTimeTickPlane0

private

Definition at line 124 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDTimeTickPlane1

private

Definition at line 126 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDTimeTickPlane2

private

Definition at line 128 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDWireInCmPlane0

private

Definition at line 129 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDWireInCmPlane1

private

Definition at line 131 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDWireInCmPlane2

private

Definition at line 133 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDWireNumberPlane0

private

Definition at line 123 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDWireNumberPlane1

private

Definition at line 125 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* vertex::FeatureVertexFinderAna::fTwoDWireNumberPlane2

private

Definition at line 127 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and beginJob().

std::string vertex::FeatureVertexFinderAna::fVertexModuleLabel

private

Definition at line 95 of file FeatureVertexFinderAna_module.cc.

Referenced by analyze(), and reconfigure().

---

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

• larreco/v06_64_02/source/larreco/VertexFinder/FeatureVertexFinderAna_module.cc