evd::SliceNavigator Class Reference

Look up products dependent on the current slice. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/EventDisplayServices/SliceNavigator.h"

Inheritance diagram for evd::SliceNavigator:

Public Member Functions
	SliceNavigator (const fhicl::ParameterSet &pset, art::ActivityRegistry &reg)
	~SliceNavigator ()
void	reconfigure (const fhicl::ParameterSet &pset) override
void	preBeginRun (art::Run const &r)
void	postProcessEvent (const art::Event &evt)
template<class T >
size_t	GetProducts (const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
template<class T >
size_t	GetProducts (const art::Event &evt, const std::string &whichTrack, const std::string &whichMichel, const std::string &instance, std::vector< const T * > &prods)
bool	FullEventMode () const
void	SetFullEventMode (bool mode)
void	Previous ()
void	Next ()
int	CurrentSlice () const
int	CurrentSliceIndex () const
void	AutoAdvance ()
void	StopAutoAdvance ()
bool	AutoAdvancing ()
void	SetCurrentSlice (int slice)
int	MaxSlice () const
std::vector< const rb::Cluster * >	SlicesToNavigate (const art::Event &evt)
void	GetBounds (const art::Event &evt, double &xmin, double &xmax, double &ymin, double &ymax, double &zmin, double &zmax, double trimfrac=0)
void	GetBox (const art::Event &evt, double &xmin, double &xmax, double &ymin, double &ymax, double &zmin, double &zmax, double trimfrac=0)
bool	DimHits () const
bool	DimOthers () const
template<>
size_t	GetProductsImplFull (const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const rb::Cluster * > &clusts)
template<>
size_t	GetProducts (const art::Event &evt, const std::string &whichTrack, const std::string &whichMichel, const std::string &instance, std::vector< const rb::Track * > &tracks)
template<>
size_t	GetProducts (const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const rb::Cluster * > &clusts)
template<>
size_t	GetProducts (const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const simb::MCTruth * > &truths)

Protected Member Functions
Bool_t	HandleTimer (TTimer *t)
template<class T >
size_t	GetProductsImpl (const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
template<class T >
size_t	GetProductsImplFull (const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
template<class T >
size_t	GetProductsImplPart (const art::Event &evt, const std::string &which, const std::string &instance, std::vector< const T * > &prods)
void	ParseModuleInstanceString (std::string &mod, std::string &inst)
std::string	stringReplace (std::string &s, const std::string &toReplace, const std::string &replaceWith)

Protected Attributes
int	fAutoAdvanceInterval
	Wait time in milliseconds. More...
std::string	fAutoPrintPattern
	Pattern for constructing filename. More...
int	fAutoPrintMax
	How many events to print (zero = disable printing) More...
int	fAutoPrintCount
int	fRun
bool	fFilterResults
	Only navigate slices that aren't filtered out. More...
bool	fDimHits
bool	fDimOthers
std::string	fSliceLabel
std::string	fCellHitLabel
std::string	fSliceInstance
std::string	fCellHitInstance
bool	fFullEvent
int	fCurrentSlice
	Index into fSlicesToNavigate. More...
int	fMaxSliceId
std::vector< int >	fSlicesToNavigate
std::vector< art::Ptr< rb::Cluster > >	fSlices
art::PtrVector< rb::CellHit >	fAllHits
int	fInitialSlice
	-1 afterwards More...
TTimer *	fTimer

Private Member Functions
void	do_reconfigure (fhicl::ParameterSet const &pset)

Detailed Description

Look up products dependent on the current slice.

Definition at line 23 of file SliceNavigator.h.

Constructor & Destructor Documentation

evd::SliceNavigator::SliceNavigator ( const fhicl::ParameterSet &  pset, art::ActivityRegistry &  reg  )

explicit

Definition at line 40 of file SliceNavigator_service.cc.

References fAutoPrintCount, fCurrentSlice, fFullEvent, fInitialSlice, fMaxSliceId, fRun, fTimer, evdb::kDRAWING_SERVICE, evdb::ServiceTable::OverrideCategory(), postProcessEvent(), preBeginRun(), and reconfigure().

    : evdb::Reconfigurable{pset},
      fAutoPrintCount(0),
      fRun(0),
      fFullEvent(true),
      fCurrentSlice(0),
      fMaxSliceId(0),
      fInitialSlice(-1),
      fTimer(0)
  {
    reconfigure(pset);
    
    evdb::ServiceTable::OverrideCategory("SliceNavigator",
                                         evdb::kDRAWING_SERVICE);

    reg.sPostProcessEvent.watch(this, &SliceNavigator::postProcessEvent);
    reg.sPreBeginRun.watch(this, &SliceNavigator::preBeginRun);
  }

evd::SliceNavigator::~SliceNavigator

(

)

Definition at line 61 of file SliceNavigator_service.cc.

References fSlicesToNavigate, and fTimer.

  {
    fSlicesToNavigate.clear();
    delete fTimer;
  }

Member Function Documentation

void evd::SliceNavigator::AutoAdvance

(

)

Definition at line 494 of file SliceNavigator_service.cc.

References fAutoAdvanceInterval, and fTimer.

Referenced by evd::SliceButtonBar::AutoAdvance().

  {
    delete fTimer;
    fTimer = new TTimer;
    fTimer->SetObject(this);
    // One-shot timer, will be re-primed by postProcessEvent()
    fTimer->Start(fAutoAdvanceInterval, true);
  }

bool evd::SliceNavigator::AutoAdvancing

(

)

Definition at line 513 of file SliceNavigator_service.cc.

References fTimer.

Referenced by evd::SliceButtonBar::AutoAdvance().

  {
    return fTimer != 0;
  }

int evd::SliceNavigator::CurrentSlice

(

)

const

Definition at line 519 of file SliceNavigator_service.cc.

References fCurrentSlice, and fSlicesToNavigate.

Referenced by postProcessEvent(), evd::SliceButtonBar::SliceButtonBar(), and evd::SliceButtonBar::Update().

  {
    if(fSlicesToNavigate.empty() || fCurrentSlice < 0) return 0;

    return fSlicesToNavigate [fCurrentSlice ];
  }

int evd::SliceNavigator::CurrentSliceIndex

(

)

const

Definition at line 527 of file SliceNavigator_service.cc.

References fCurrentSlice.

Referenced by evd::RawDataDrawer::RawDigit2D(), evd::RawDataDrawer::RawDigit3D(), and evd::RawDataDrawer::RawDigitHardware().

  {
    return fCurrentSlice;
  }

bool evd::SliceNavigator::DimHits ( ) const

inline

Definition at line 75 of file SliceNavigator.h.

References fDimHits.

Referenced by evd::RecoBaseDrawer::CellHit2D().

{return fDimHits;}

bool evd::SliceNavigator::DimOthers ( ) const

inline

Definition at line 76 of file SliceNavigator.h.

References fDimOthers, GetProductsImpl(), GetProductsImplFull(), GetProductsImplPart(), HandleTimer(), APDHVSetting::mod, ParseModuleInstanceString(), string, stringReplace(), and confusionMatrixTree::t.

Referenced by evd::RawDataDrawer::RawDigit2D(), evd::RawDataDrawer::RawDigit3D(), and evd::RawDataDrawer::RawDigitHardware().

{return fDimOthers;}

void evdb::Reconfigurable::do_reconfigure ( fhicl::ParameterSet const &  pset )

inlineinherited

Definition at line 17 of file Reconfigurable.h.

{ reconfigure(pset); }

bool evd::SliceNavigator::FullEventMode

(

)

const

Definition at line 447 of file SliceNavigator_service.cc.

References fFullEvent.

Referenced by evd::RawDataDrawer::RawDigit2D(), evd::RawDataDrawer::RawDigit3D(), evd::RawDataDrawer::RawDigitHardware(), and evd::SliceButtonBar::SliceButtonBar().

  {
    return fFullEvent;
  }

void evd::SliceNavigator::GetBounds	(	const art::Event &	evt,
		double &	xmin,
		double &	xmax,
		double &	ymin,
		double &	ymax,
		double &	zmin,
		double &	zmax,
		double	trimfrac = 0
	)

Definition at line 556 of file SliceNavigator_service.cc.

References rb::CellHit::Cell(), rb::Cluster::Cell(), om::cerr, geo::GeometryBase::DetHalfHeight(), geo::GeometryBase::DetHalfWidth(), geo::GeometryBase::DetLength(), allTimeWatchdog::endl, art::Handle< T >::failedToGet(), fCurrentSlice, fSliceInstance, fSliceLabel, fSlicesToNavigate, geom(), art::DataViewImpl::getByLabel(), MECModelEnuComparisons::i, geo::kX, rb::Cluster::MaxXYZ(), rb::Cluster::MinXYZ(), rb::Cluster::NCell(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), and rb::CellHit::View().

Referenced by evd::TZProjPad::AutoZoomSlice().

  {
    // Fallback to the whole detector in case we don't get values below
    art::ServiceHandle<geo::Geometry> geom;
    xmin = -geom->DetHalfWidth();
    xmax = +geom->DetHalfWidth();
    ymin = -geom->DetHalfHeight();
    ymax = +geom->DetHalfHeight();
    zmin = 0;
    zmax = geom->DetLength();

    art::Handle<std::vector<rb::Cluster>> slices;
    evt.getByLabel( art::InputTag(fSliceLabel, fSliceInstance), slices);

    if(slices.failedToGet()){
      std::cerr << "Couldn't find any slices under label "
                << fSliceLabel << std::endl;
      return;
    }

    if(fCurrentSlice >= 0 && fCurrentSlice < int(slices->size())){
      const rb::Cluster& slice = (*slices)[fSlicesToNavigate[fCurrentSlice]];

      if(trimfrac == 0){
        const TVector3 r0 = slice.MinXYZ();
        const TVector3 r1 = slice.MaxXYZ();
        xmin = r0.X(); xmax = r1.X();
        ymin = r0.Y(); ymax = r1.Y();
        zmin = r0.Z(); zmax = r1.Z();
        return;
      }

      std::vector<double> xs, ys, zs;

      for(unsigned int i = 0; i < slice.NCell(); ++i){
        const art::Ptr<rb::CellHit>& chit = slice.Cell(i);
        double xyz[3];
        geom->Plane(chit->Plane())->Cell(chit->Cell())->GetCenter(xyz);
        zs.push_back(xyz[2]);
        if(chit->View() == geo::kX)
          xs.push_back(xyz[0]);
        else
          ys.push_back(xyz[1]);
      }

      std::sort(xs.begin(), xs.end());
      std::sort(ys.begin(), ys.end());
      std::sort(zs.begin(), zs.end());

      xmin = xs[xs.size()*trimfrac];
      ymin = ys[ys.size()*trimfrac];
      zmin = zs[zs.size()*trimfrac];

      // For small events we want to pick up the last hit, not the N-1th hit
      xmax = xs[(xs.size()-1)*(1-trimfrac)+1];
      ymax = ys[(ys.size()-1)*(1-trimfrac)+1];
      zmax = zs[(zs.size()-1)*(1-trimfrac)+1];
    }
  }

void evd::SliceNavigator::GetBox	(	const art::Event &	evt,
		double &	xmin,
		double &	xmax,
		double &	ymin,
		double &	ymax,
		double &	zmin,
		double &	zmax,
		double	trimfrac = 0
	)

Definition at line 621 of file SliceNavigator_service.cc.

References rb::CellHit::Cell(), rb::Cluster::Cell(), om::cerr, geo::GeometryBase::DetHalfHeight(), geo::GeometryBase::DetHalfWidth(), geo::GeometryBase::DetLength(), allTimeWatchdog::endl, art::Handle< T >::failedToGet(), fCurrentSlice, evd::GeometryDrawingOptions::fSliceBoxSize, fSliceInstance, fSliceLabel, fSlicesToNavigate, geom(), art::DataViewImpl::getByLabel(), MECModelEnuComparisons::i, geo::kX, rb::Cluster::MaxXYZ(), rb::Cluster::MinXYZ(), rb::Cluster::NCell(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), nd_projection_maker::sizes, rb::CellHit::View(), make_mec_shifts_plots::xmin, and ymin.

Referenced by evd::TZProjPad::AutoZoomBox().

  {
    // Fallback to the whole detector in case we don't get values below
    art::ServiceHandle<geo::Geometry> geom;
    xmin = -geom->DetHalfWidth();
    xmax = +geom->DetHalfWidth();
    ymin = -geom->DetHalfHeight();
    ymax = +geom->DetHalfHeight();
    zmin = 0;
    zmax = geom->DetLength();

    // Use Fiducial dimensions if smaller than slice
    art::ServiceHandle<evd::GeometryDrawingOptions> drawopt;

    float xbox=0;
    float ybox=0;
    float zbox=0;
    if ( (drawopt->fSliceBoxSize).size() == 3 ) {
      std::vector< float > sizes = drawopt->fSliceBoxSize;
      xbox = sizes[0];
      ybox = sizes[1];
      zbox = sizes[2];
    }
    art::Handle<std::vector<rb::Cluster>> slices;
    evt.getByLabel( art::InputTag(fSliceLabel, fSliceInstance), slices);

    if(slices.failedToGet()){
      std::cerr << "Couldn't find any slices under label "
                << fSliceLabel << std::endl;
      return;
    }

    if(fCurrentSlice >= 0 && fCurrentSlice < int(slices->size())){
      const rb::Cluster& slice = (*slices)[fSlicesToNavigate[fCurrentSlice]];

      if(trimfrac == 0){
        const TVector3 r0 = slice.MinXYZ();
        const TVector3 r1 = slice.MaxXYZ();
        xmin = r0.X(); xmax = r1.X();
        ymin = r0.Y(); ymax = r1.Y();
        zmin = r0.Z(); zmax = r1.Z();
        return;
      }

      std::vector<double> xs, ys, zs;

      for(unsigned int i = 0; i < slice.NCell(); ++i){
        const art::Ptr<rb::CellHit>& chit = slice.Cell(i);
        double xyz[3];
        geom->Plane(chit->Plane())->Cell(chit->Cell())->GetCenter(xyz);
        zs.push_back(xyz[2]);
        if(chit->View() == geo::kX)
          xs.push_back(xyz[0]);
        else
          ys.push_back(xyz[1]);
      }
      std::sort(xs.begin(), xs.end());
      std::sort(ys.begin(), ys.end());
      std::sort(zs.begin(), zs.end());

      xmin = xs[xs.size()*trimfrac];
      ymin = ys[ys.size()*trimfrac];
      zmin = zs[zs.size()*trimfrac];

      // For small events we want to pick up the last hit, not the N-1th hit
      xmax = xs[(xs.size()-1)*(1-trimfrac)+1];
      ymax = ys[(ys.size()-1)*(1-trimfrac)+1];
      zmax = zs[(zs.size()-1)*(1-trimfrac)+1];

      double xmid=(xmax+xmin)/2;
      double ymid=(ymax+ymin)/2;
      double zmid=(zmax+zmin)/2;
      xmax= ( xbox==0 ) ? xmax : xmid+xbox/2;
      ymax= ( ybox==0 ) ? ymax : ymid+ybox/2;
      zmax= ( zbox==0 ) ? zmax : zmid+zbox/2;
      xmin= ( xbox==0 ) ? xmin : xmid-xbox/2;
      ymin= ( ybox==0 ) ? ymin : ymid-ybox/2;
      zmin= ( zbox==0 ) ? zmin : zmid-zbox/2;
    }
  }

template<class T >

size_t evd::SliceNavigator::GetProducts	(	const art::Event &	evt,
		const std::string &	which,
		const std::string &	instance,
		std::vector< const T * > &	prods
	)

Definition at line 219 of file SliceNavigator_service.cc.

References GetProductsImpl().

Referenced by evd::RecoBaseDrawer::CellHit2D(), evd::RecoBaseDrawer::CellHit3D(), evd::RecoBaseDrawer::Cluster2D(), evd::RecoBaseDrawer::FillTQHisto(), evd::SimulationDrawer::GetMCTruth(), GetProducts(), GetProductsImplPart(), evd::RecoBaseDrawer::HoughResult2D(), evd::RecoBaseDrawer::OfflineChans2D(), postProcessEvent(), evd::RecoBaseDrawer::Prong2D(), evd::RecoBaseDrawer::Prong3D(), evd::RecoBaseDrawer::Track2D(), and evd::RecoBaseDrawer::Vertex2D().

  {
    return GetProductsImpl(evt, which, instance, prods);
  }

template<class T >

size_t evd::SliceNavigator::GetProducts	(	const art::Event &	evt,
		const std::string &	whichTrack,
		const std::string &	whichMichel,
		const std::string &	instance,
		std::vector< const T * > &	prods
	)

template<>

size_t evd::SliceNavigator::GetProducts	(	const art::Event &	evt,
		const std::string &	whichTrack,
		const std::string &	whichMichel,
		const std::string &	instance,
		std::vector< const rb::Track * > &	tracks
	)

Definition at line 365 of file SliceNavigator_service.cc.

References plot_validation_datamc::c, GetProducts(), and GetProductsImpl().

  {
    std::vector<const me::TrkME*> clusts;
    GetProductsImpl(evt, whichTrack, instance, tracks);
    GetProductsImpl(evt, whichMichel, instance, clusts);

    for(const me::TrkME* c: clusts)
      tracks.push_back(c->ParentTrk().get());

    return tracks.size();
  }

template<>

size_t evd::SliceNavigator::GetProducts	(	const art::Event &	evt,
		const std::string &	which,
		const std::string &	instance,
		std::vector< const rb::Cluster * > &	clusts
	)

Definition at line 384 of file SliceNavigator_service.cc.

References om::cerr, allTimeWatchdog::endl, art::Handle< T >::failedToGet(), fCurrentSlice, fFilterResults, fFullEvent, fSliceLabel, fSlicesToNavigate, art::DataViewImpl::getByLabel(), GetProducts(), GetProductsImplFull(), compare_h5_caf::idx, and rb::IsFiltered().

  {
    // If it's the actual slices, look them up by index.
    // Otherwise we'll get nothing
    // (not possible to associate Cluster with Cluster)
    if(which == fSliceLabel && instance == ""){
      art::Handle<std::vector<rb::Cluster>> slices;
      evt.getByLabel(fSliceLabel, slices);
      if(slices.failedToGet()){
        std::cerr << "Couldn't find any slices under label "
                  << fSliceLabel << std::endl;
        return 0;
      }
      if( fFullEvent ){
        int nSli = slices->size();
        for( int iSli=0; iSli < nSli; ++iSli){
          if ( fFilterResults && rb::IsFiltered(evt,slices, iSli ) )
          continue;
          clusts.push_back((art::Ptr<rb::Cluster>(slices, iSli).get())) ;
        }
      }// if full event
      else{
        int idx = fSlicesToNavigate[ fCurrentSlice];
        if( fCurrentSlice >=0 )
        clusts.push_back( (art::Ptr<rb::Cluster>(slices, idx)).get() );
      }
    }
    else
      GetProductsImplFull(evt, which, instance, clusts);

    return clusts.size();
  }

template<>

size_t evd::SliceNavigator::GetProducts	(	const art::Event &	evt,
		const std::string &	which,
		const std::string &	instance,
		std::vector< const simb::MCTruth * > &	truths
	)

Definition at line 423 of file SliceNavigator_service.cc.

References fAllHits, fCurrentSlice, fFullEvent, fSlices, fSlicesToNavigate, GetProductsImpl(), and cheat::BackTracker::SliceToMCTruth().

  {
    truths.clear();

    if(fFullEvent){
      GetProductsImpl(evt, which, instance, truths);
    }
    else{
      art::ServiceHandle<cheat::BackTracker> bt;

      std::vector<cheat::NeutrinoEffPur> eps =
        bt->SliceToMCTruth(fSlices[fSlicesToNavigate[fCurrentSlice]],
                           fAllHits);

      if(!eps.empty()) truths = {eps.front().neutrinoInt.get()};
    }

    return truths.size();
  }

template<class T >

size_t evd::SliceNavigator::GetProductsImpl ( const art::Event &  evt, const std::string &  which, const std::string &  instance, std::vector< const T * > &  prods  )

protected

Definition at line 229 of file SliceNavigator_service.cc.

References fFullEvent, GetProductsImplFull(), and GetProductsImplPart().

Referenced by DimOthers(), and GetProducts().

  {
    if(fFullEvent)
      return GetProductsImplFull(evt, which, instance, prods);
    else
      return GetProductsImplPart(evt, which, instance, prods);
  }

template<class T >

size_t evd::SliceNavigator::GetProductsImplFull ( const art::Event &  evt, const std::string &  which, const std::string &  instance, std::vector< const T * > &  prods  )

protected

Definition at line 266 of file SliceNavigator_service.cc.

References om::cerr, distance(), e, allTimeWatchdog::endl, fFilterResults, art::fill_ptr_vector(), fSliceInstance, fSliceLabel, art::DataViewImpl::getByLabel(), GetProductsImplPart(), MECModelEnuComparisons::i, rb::IsFiltered(), T, and APDHVSetting::temp.

Referenced by DimOthers(), GetProducts(), GetProductsImpl(), and GetProductsImplFull().

  {

    art::Handle<std::vector<T> > tempHandle;

    std::vector<art::Ptr<T> > temp;

    try{
      evt.getByLabel(art::InputTag(which, instance), tempHandle);
      art::fill_ptr_vector( temp,tempHandle);
    }
    catch(cet::exception& e){
      std::cerr << "Failed to get products of type " << typeid(T).name()
                << " under label \"" << which << "\": "
                << e << std::endl;
    }

    std::vector<art::Ptr<T> > tempcopy;
    // If fFilterResults is turned on, we don't want to
    // draw anything from a slice that was filtered out
    // So each product has to be mapped to its associated
    // slice cluster and check if the slice was filtered out.
    if( fFilterResults &&
        typeid(T) != typeid(rb::CellHit) &&
        typeid(T) != typeid(simb::MCTruth)){

      art::Handle<std::vector<rb::Cluster> > slicesHandle;
      evt.getByLabel( art::InputTag(fSliceLabel, fSliceInstance), slicesHandle);
      std::vector<art::Ptr<rb::Cluster> >    slices;
      art::fill_ptr_vector( slices, slicesHandle);
      art::FindOneP<rb::Cluster> fot(temp, evt, which);

      if(!fot.isValid()){
        std::cerr << "Nothing of type " << typeid(T).name()
        << " with label \"" << which << "\" and instance \"" << instance
        << "\" associated with slices" << std::endl;
        return 0;
      }

      int nT = temp.size();
      for( int i = 0; i < nT; i++){

        art::Ptr<rb::Cluster> slice = fot.at( i );
        std::vector< art::Ptr< rb::Cluster> >::iterator slicePos =
        std::find(slices.begin(), slices.end(), slice) ;
        int iSlice = std::distance( slices.begin(), slicePos );
        
        if( fFilterResults && rb::IsFiltered(evt,slicesHandle, iSlice ) )
        continue;
        tempcopy.push_back( temp[i] );
        
      }// end loop over temp
    }
    else
      tempcopy.swap( temp );

    for( art::Ptr<T> p: tempcopy)
      prods.push_back( p.get());

    return prods.size();
  }

template<>

size_t evd::SliceNavigator::GetProductsImplFull	(	const art::Event &	evt,
		const std::string &	which,
		const std::string &	instance,
		std::vector< const rb::Cluster * > &	clusts
	)

Definition at line 246 of file SliceNavigator_service.cc.

References om::cerr, e, allTimeWatchdog::endl, GetProductsImplFull(), art::DataViewImpl::getView(), and APDHVSetting::temp.

  {
    std::vector<const rb::Cluster*> temp(clusts);
    try{
      evt.getView(which, instance, temp);
      temp.swap(clusts);
    }
    catch(cet::exception& e){
      std::cerr << "Failed to get products of type " << typeid(rb::Cluster).name()
                << " under label \"" << which << "\": "
                << e << std::endl;
    }
    return clusts.size();
  }

template<class T >

size_t evd::SliceNavigator::GetProductsImplPart ( const art::Event &  evt, const std::string &  which, const std::string &  instance, std::vector< const T * > &  prods  )

protected

Definition at line 334 of file SliceNavigator_service.cc.

References om::cerr, allTimeWatchdog::endl, art::Handle< T >::failedToGet(), fCurrentSlice, PandAna.Demos.pi0_spectra::fmt, fSliceInstance, fSliceLabel, fSlicesToNavigate, art::DataViewImpl::getByLabel(), GetProducts(), and T.

Referenced by DimOthers(), GetProductsImpl(), and GetProductsImplFull().

  {
    art::Handle<std::vector<rb::Cluster>> slices;
    evt.getByLabel(art::InputTag(fSliceLabel, fSliceInstance), slices);

    if(slices.failedToGet()){
      std::cerr << "Couldn't find any slices under label "
                << fSliceLabel << std::endl;
      return 0;
    }

    art::FindManyP<T> fmt(slices, evt, art::InputTag(which, instance));
    if(!fmt.isValid()){
      std::cerr << "Nothing of type " << typeid(T).name()
                << " with label \"" << which << "\" and instance \"" << instance
                << "\" associated with slices" << std::endl;
      return 0;
    }

    std::vector<art::Ptr<T>> prodPtrs = fmt.at( fSlicesToNavigate[fCurrentSlice] );
    prods.clear();
    for(art::Ptr<T> p: prodPtrs) prods.push_back(p.get());
    return prods.size();
  }

Bool_t evd::SliceNavigator::HandleTimer ( TTimer *  t )

protected

Definition at line 729 of file SliceNavigator_service.cc.

References Next().

Referenced by DimOthers().

  {
    // This is the auto-advance timer
    Next();
    return true;
  }

int evd::SliceNavigator::MaxSlice

(

)

const

Definition at line 550 of file SliceNavigator_service.cc.

References fMaxSliceId.

Referenced by evd::SliceButtonBar::Update().

  {
    return fMaxSliceId;
  }

void evd::SliceNavigator::Next

(

)

Definition at line 478 of file SliceNavigator_service.cc.

References fCurrentSlice, fSlicesToNavigate, evdb::kNEXT_EVENT, evdb::kRELOAD_EVENT, and evdb::NavState::Set().

Referenced by HandleTimer(), and evd::SliceButtonBar::Next().

  {
    ++fCurrentSlice;

    if(fCurrentSlice >= int(fSlicesToNavigate.size())){
      fCurrentSlice = -2;
      fSlicesToNavigate.clear();
      evdb::NavState::Set(evdb::kNEXT_EVENT);
    }
    else{
      // Go to the slice requested
      evdb::NavState::Set(evdb::kRELOAD_EVENT);
    }
  }

void evd::SliceNavigator::ParseModuleInstanceString ( std::string &  mod, std::string &  inst  )

protected

Definition at line 88 of file SliceNavigator_service.cc.

Referenced by DimOthers(), and reconfigure().

  {
    const size_t splitPos = mod.find_first_of('_');
    if(splitPos == std::string::npos){
      inst = "";
      return;
    }
    inst = mod.substr(splitPos+1);
    mod.resize(splitPos);
  }

void evd::SliceNavigator::postProcessEvent

(

const art::Event &

evt

)

Definition at line 107 of file SliceNavigator_service.cc.

References om::cerr, art::PtrVector< T >::clear(), om::cout, CurrentSlice(), allTimeWatchdog::endl, art::Event::event(), exit(), art::Handle< T >::failedToGet(), fAllHits, fAutoAdvanceInterval, fAutoPrintCount, fAutoPrintMax, fAutoPrintPattern, fCellHitInstance, fCellHitLabel, fCurrentSlice, fFilterResults, fFullEvent, shutoffs::filename, fInitialSlice, fMaxSliceId, fRun, fSliceInstance, fSliceLabel, fSlices, fSlicesToNavigate, fTimer, art::DataViewImpl::getByLabel(), evdb::Printable::GetPrintables(), GetProducts(), MECModelEnuComparisons::i, rb::IsFiltered(), it, evdb::kNEXT_EVENT, evdb::kPREV_EVENT, evdb::Printable::Print(), nd_projection_maker::ps, art::PtrVector< T >::push_back(), evdb::NavState::Set(), string, stringReplace(), and art::to_string().

Referenced by SliceNavigator().

  {
    // Clear the slices to navigate vector
    fSlicesToNavigate.clear();

    art::Handle<std::vector<rb::Cluster>> slices;
    
    art::InputTag sliceTag(fSliceLabel, fSliceInstance);
    
    evt.getByLabel( sliceTag, slices);

    if(slices.failedToGet()){
      std::cerr << "Couldn't find any slices under label "
                << fSliceLabel << std::endl;
      fMaxSliceId = 0;
      return;
    }

    int nSli = slices->size();
    for( int iSli = 0; iSli < nSli; ++iSli){
      // if we only want to look at slices that
      // survived filtering, throw the rest away
      if(!fFullEvent && fFilterResults && rb::IsFiltered(evt,slices,iSli) )
      continue;
      if(!fFullEvent && (*slices)[iSli].IsNoise())
      continue;

      fSlicesToNavigate.push_back( iSli );
      if(iSli == fInitialSlice){
        fCurrentSlice = fSlicesToNavigate.size()-1;
        fInitialSlice = -1; // Don't do this again
      }
    }
    if (fSlicesToNavigate.size() > 0) fMaxSliceId = fSlicesToNavigate.back();

    // Protecting against segfaults:  If we've advanced to the next event
    // with fCurrentSlice set to the max number of slices, and the next event
    // has fewer slices than the previous, we need to make sure we don't 
    // suffer from an index out of range problem.
    if( (fCurrentSlice + 1 > int(fSlicesToNavigate.size())) && !fSlicesToNavigate.empty() ) 
      fCurrentSlice = fSlicesToNavigate.size() - 1;

    if(fCurrentSlice == -1){
      if(!fSlicesToNavigate.empty())
        // We pressed the prev slice button
        fCurrentSlice = fSlicesToNavigate.size()-1;
      else
        evdb::NavState::Set(evdb::kPREV_EVENT);
    }

    if(fCurrentSlice == -2 ){
      if(!fSlicesToNavigate.empty())
        // We pressed the next button, but there
        // no more slices left in the previous event
        fCurrentSlice = 0;
      else
        evdb::NavState::Set(evdb::kNEXT_EVENT);
    }

    art::Handle<std::vector<rb::Cluster>> slicecol;
    evt.getByLabel(sliceTag, slicecol);
    fSlices.clear();
    for(unsigned int i = 0; i < slicecol->size(); ++i)
      fSlices.push_back(art::Ptr<rb::Cluster>(slicecol, i));

    art::Handle<std::vector<rb::CellHit>> hitcol;
    art::InputTag cellTag(fCellHitLabel, fCellHitInstance);
    evt.getByLabel(cellTag, hitcol);
    fAllHits.clear();
    for(unsigned int i = 0; i < hitcol->size(); ++i)
      fAllHits.push_back(art::Ptr<rb::CellHit>(hitcol, i));
      
    // Auto-printing
    if (fAutoPrintMax > 0 || fAutoPrintMax == -1){
      ++fAutoPrintCount;
      std::map<std::string, evdb::Printable*>& ps = evdb::Printable::GetPrintables();
      for(std::map<std::string, evdb::Printable*>::iterator it = ps.begin(); it != ps.end(); ++it){
        evdb::Printable* p = it->second;
        
        std::string filename = fAutoPrintPattern;
        // std::string filename = "evd-slice_r%r_e%e_s%s.gif";
        
        // Insert instances of run variable into filename
        if (filename.find("%r") != std::string::npos)
          stringReplace(filename, "%r", std::to_string(fRun));
        
        // Insert instances of event variable into filename
        if (filename.find("%e") != std::string::npos)
          stringReplace(filename, "%e", std::to_string(evt.event()));
        
        // Insert instances of slice variable into filename
        if (filename.find("%s") != std::string::npos)
          stringReplace(filename, "%s", std::to_string(this->CurrentSlice()));
        
        // Print file
        p->Print(TString(filename));
      }
      if (fAutoPrintCount >= (fAutoPrintMax - 1) && fAutoPrintMax != -1)
      {
        std::cout << "Hit autoprint max and exiting." << std::endl;
        exit(0);
      }
    }

    if(fTimer){
      // Auto-advance is active. Fire off one more shot
      fTimer->Start(fAutoAdvanceInterval, true);
    }
  }

void evd::SliceNavigator::preBeginRun

(

art::Run const &

r

)

Definition at line 101 of file SliceNavigator_service.cc.

References fRun, and art::Run::run().

Referenced by SliceNavigator().

  {
    fRun = r.run();
  }

void evd::SliceNavigator::Previous

(

)

Definition at line 462 of file SliceNavigator_service.cc.

References fCurrentSlice, fSlicesToNavigate, evdb::kPREV_EVENT, evdb::kRELOAD_EVENT, and evdb::NavState::Set().

Referenced by evd::SliceButtonBar::Previous().

  {
    --fCurrentSlice;

    if(fCurrentSlice < 0){
      fCurrentSlice = -1;
      evdb::NavState::Set(evdb::kPREV_EVENT);
      fSlicesToNavigate.clear();
    }
    else{
      // Go to the slice requested
      evdb::NavState::Set(evdb::kRELOAD_EVENT);
    }
  }

void evd::SliceNavigator::reconfigure ( const fhicl::ParameterSet &  pset )

override

Definition at line 68 of file SliceNavigator_service.cc.

References fAutoAdvanceInterval, fAutoPrintMax, fAutoPrintPattern, fCellHitInstance, fCellHitLabel, fDimHits, fDimOthers, fFilterResults, fFullEvent, fInitialSlice, fSliceInstance, fSliceLabel, fhicl::ParameterSet::get(), ParseModuleInstanceString(), and string.

Referenced by SliceNavigator().

  {
    fSliceLabel          = pset.get<std::string>("SliceLabel");
    fCellHitLabel        = pset.get<std::string>("CellHitLabel");
    fAutoAdvanceInterval = pset.get<int>("AutoAdvanceInterval");
    fAutoPrintPattern    = pset.get<std::string>("AutoPrintPattern");
    fAutoPrintMax        = pset.get<int>("AutoPrintMax");
    fFilterResults       = (pset.get<int>("SliceOpt.val") & 1);
    fDimHits             = (pset.get<int>("SliceOpt.val") & 2);
    fDimOthers           = (pset.get<int>("SliceOpt.val") & 4);

    fFullEvent = !pset.get<bool>("InitiallyEnabled");
    fInitialSlice = pset.get<int>("InitialSlice");

    this->ParseModuleInstanceString(fSliceLabel, fSliceInstance);
    this->ParseModuleInstanceString(fCellHitLabel, fCellHitInstance);
    
  }

void evd::SliceNavigator::SetCurrentSlice

(

int

slice

)

Definition at line 533 of file SliceNavigator_service.cc.

References fCurrentSlice, fSlicesToNavigate, it, evdb::kRELOAD_EVENT, and evdb::NavState::Set().

Referenced by evd::SliceButtonBar::GoTo().

  {
    auto it = std::find(fSlicesToNavigate.begin(),
                        fSlicesToNavigate.end(),
                        slice);
    if(it != fSlicesToNavigate.end()){
      fCurrentSlice = it-fSlicesToNavigate.begin();
      evdb::NavState::Set(evdb::kRELOAD_EVENT);
    }
    else{
      // What else can we do?
      fCurrentSlice = 0;
      evdb::NavState::Set(evdb::kRELOAD_EVENT);
    }
  }

void evd::SliceNavigator::SetFullEventMode

(

bool

mode

)

Definition at line 453 of file SliceNavigator_service.cc.

References fFullEvent, evdb::kRELOAD_EVENT, submit_nova_art::mode, and evdb::NavState::Set().

Referenced by evd::SliceButtonBar::Enable().

  {
    if(mode == fFullEvent) return;

    fFullEvent = mode;
    evdb::NavState::Set(evdb::kRELOAD_EVENT);
  }

std::vector< const rb::Cluster * > evd::SliceNavigator::SlicesToNavigate

(

const art::Event &

evt

)

Definition at line 707 of file SliceNavigator_service.cc.

References art::Handle< T >::failedToGet(), fSliceInstance, fSliceLabel, fSlicesToNavigate, art::DataViewImpl::getByLabel(), and MECModelEnuComparisons::i.

Referenced by evd::RecoBaseDrawer::CellHit2D(), evd::RawDataDrawer::RawDigit2D(), evd::RawDataDrawer::RawDigit3D(), and evd::RawDataDrawer::RawDigitHardware().

  {
    std::vector<const rb::Cluster*> clusts;
    art::Handle<std::vector<rb::Cluster>> slices;
    evt.getByLabel(art::InputTag(fSliceLabel, fSliceInstance), slices);
    if(slices.failedToGet()) return {};

    for(unsigned int i = 0; i < fSlicesToNavigate.size(); ++i)
      clusts.push_back(&(*slices)[fSlicesToNavigate[i]]);

    return clusts;
  }// End of SlicesToNavigate

void evd::SliceNavigator::StopAutoAdvance

(

)

Definition at line 504 of file SliceNavigator_service.cc.

References fTimer.

Referenced by evd::SliceButtonBar::AutoAdvance().

  {
    if(!fTimer) return;
    fTimer->Stop();
    delete fTimer;
    fTimer = 0;
  }

std::string evd::SliceNavigator::stringReplace ( std::string &  s, const std::string &  toReplace, const std::string &  replaceWith  )

protected

Definition at line 721 of file SliceNavigator_service.cc.

Referenced by DimOthers(), and postProcessEvent().

  {
      return(s.replace(s.find(toReplace), toReplace.length(), replaceWith));
  }

Member Data Documentation

art::PtrVector<rb::CellHit> evd::SliceNavigator::fAllHits

protected

Definition at line 127 of file SliceNavigator.h.

Referenced by GetProducts(), and postProcessEvent().

int evd::SliceNavigator::fAutoAdvanceInterval

protected

Wait time in milliseconds.

Definition at line 106 of file SliceNavigator.h.

Referenced by AutoAdvance(), postProcessEvent(), and reconfigure().

int evd::SliceNavigator::fAutoPrintCount

protected

Definition at line 110 of file SliceNavigator.h.

Referenced by postProcessEvent(), and SliceNavigator().

int evd::SliceNavigator::fAutoPrintMax

protected

How many events to print (zero = disable printing)

Definition at line 108 of file SliceNavigator.h.

Referenced by postProcessEvent(), and reconfigure().

std::string evd::SliceNavigator::fAutoPrintPattern

protected

Pattern for constructing filename.

Definition at line 107 of file SliceNavigator.h.

Referenced by postProcessEvent(), and reconfigure().

std::string evd::SliceNavigator::fCellHitInstance

protected

Definition at line 120 of file SliceNavigator.h.

Referenced by postProcessEvent(), and reconfigure().

std::string evd::SliceNavigator::fCellHitLabel

protected

Definition at line 118 of file SliceNavigator.h.

Referenced by postProcessEvent(), and reconfigure().

int evd::SliceNavigator::fCurrentSlice

protected

Index into fSlicesToNavigate.

Definition at line 122 of file SliceNavigator.h.

Referenced by CurrentSlice(), CurrentSliceIndex(), GetBounds(), GetBox(), GetProducts(), GetProductsImplPart(), Next(), postProcessEvent(), Previous(), SetCurrentSlice(), and SliceNavigator().

bool evd::SliceNavigator::fDimHits

protected

Definition at line 115 of file SliceNavigator.h.

Referenced by DimHits(), and reconfigure().

bool evd::SliceNavigator::fDimOthers

protected

Definition at line 116 of file SliceNavigator.h.

Referenced by DimOthers(), and reconfigure().

bool evd::SliceNavigator::fFilterResults

protected

Only navigate slices that aren't filtered out.

Definition at line 114 of file SliceNavigator.h.

Referenced by GetProducts(), GetProductsImplFull(), postProcessEvent(), and reconfigure().

bool evd::SliceNavigator::fFullEvent

protected

Definition at line 121 of file SliceNavigator.h.

Referenced by FullEventMode(), GetProducts(), GetProductsImpl(), postProcessEvent(), reconfigure(), SetFullEventMode(), and SliceNavigator().

int evd::SliceNavigator::fInitialSlice

protected

-1 afterwards

Definition at line 129 of file SliceNavigator.h.

Referenced by postProcessEvent(), reconfigure(), and SliceNavigator().

int evd::SliceNavigator::fMaxSliceId

protected

Definition at line 123 of file SliceNavigator.h.

Referenced by MaxSlice(), postProcessEvent(), and SliceNavigator().

int evd::SliceNavigator::fRun

protected

Definition at line 112 of file SliceNavigator.h.

Referenced by postProcessEvent(), preBeginRun(), and SliceNavigator().

std::string evd::SliceNavigator::fSliceInstance

protected

Definition at line 119 of file SliceNavigator.h.

Referenced by GetBounds(), GetBox(), GetProductsImplFull(), GetProductsImplPart(), postProcessEvent(), reconfigure(), and SlicesToNavigate().

std::string evd::SliceNavigator::fSliceLabel

protected

Definition at line 117 of file SliceNavigator.h.

Referenced by GetBounds(), GetBox(), GetProducts(), GetProductsImplFull(), GetProductsImplPart(), postProcessEvent(), reconfigure(), and SlicesToNavigate().

std::vector<art::Ptr<rb::Cluster> > evd::SliceNavigator::fSlices

protected

Definition at line 126 of file SliceNavigator.h.

Referenced by GetProducts(), and postProcessEvent().

std::vector<int> evd::SliceNavigator::fSlicesToNavigate

protected

Definition at line 124 of file SliceNavigator.h.

Referenced by CurrentSlice(), GetBounds(), GetBox(), GetProducts(), GetProductsImplPart(), Next(), postProcessEvent(), Previous(), SetCurrentSlice(), SlicesToNavigate(), and ~SliceNavigator().

TTimer* evd::SliceNavigator::fTimer

protected

Definition at line 131 of file SliceNavigator.h.

Referenced by AutoAdvance(), AutoAdvancing(), postProcessEvent(), SliceNavigator(), StopAutoAdvance(), and ~SliceNavigator().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/EventDisplayServices/SliceNavigator.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/EventDisplayServices/SliceNavigator_service.cc