evd::SliceNavigator Class Reference

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

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-12-03/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	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)

Public 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

Protected Member Functions
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	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...

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 38 of file SliceNavigator_service.cc.

References fAutoPrintCount, fCurrentSlice, fFullEvent, fInitialSlice, fMaxSliceId, fRun, 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)
  {
    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 58 of file SliceNavigator_service.cc.

References fSlicesToNavigate.

  {
    fSlicesToNavigate.clear();
  }

Member Function Documentation

int evd::SliceNavigator::CurrentSlice

(

)

const

Definition at line 481 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 489 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 71 of file SliceNavigator.h.

References fDimHits.

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

{return fDimHits;}

bool evd::SliceNavigator::DimOthers ( ) const

inline

Definition at line 72 of file SliceNavigator.h.

References fDimOthers.

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 434 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 518 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 583 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 206 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 352 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 371 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 410 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 216 of file SliceNavigator_service.cc.

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

Referenced by 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 253 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 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 233 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 321 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 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();
  }

int evd::SliceNavigator::MaxSlice

(

)

const

Definition at line 512 of file SliceNavigator_service.cc.

References fMaxSliceId.

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

  {
    return fMaxSliceId;
  }

void evd::SliceNavigator::Next

(

)

Definition at line 465 of file SliceNavigator_service.cc.

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

Referenced by 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 84 of file SliceNavigator_service.cc.

Referenced by 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 103 of file SliceNavigator_service.cc.

References om::cerr, art::PtrVector< T >::clear(), CurrentSlice(), allTimeWatchdog::endl, art::Event::event(), exit(), art::Handle< T >::failedToGet(), fAllHits, fAutoPrintCount, fAutoPrintMax, fCellHitInstance, fCellHitLabel, fCurrentSlice, fFilterResults, fFullEvent, shutoffs::filename, fInitialSlice, fMaxSliceId, fRun, fSliceInstance, fSliceLabel, fSlices, fSlicesToNavigate, 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 (fCurrentSlice >= (fAutoPrintMax - 1) && fAutoPrintMax != -1) exit(0);
    }
  }

void evd::SliceNavigator::preBeginRun

(

art::Run const &

r

)

Definition at line 97 of file SliceNavigator_service.cc.

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

Referenced by SliceNavigator().

  {
    fRun = r.run();
  }

void evd::SliceNavigator::Previous

(

)

Definition at line 449 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 64 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 495 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 440 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 669 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

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

protected

Definition at line 683 of file SliceNavigator_service.cc.

Referenced by postProcessEvent().

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

Member Data Documentation

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

protected

Definition at line 121 of file SliceNavigator.h.

Referenced by GetProducts(), and postProcessEvent().

int evd::SliceNavigator::fAutoAdvanceInterval

Wait time in milliseconds.

Definition at line 74 of file SliceNavigator.h.

Referenced by evd::SliceButtonBar::AutoAdvance(), evd::SliceButtonBar::HandleTimer(), and reconfigure().

int evd::SliceNavigator::fAutoPrintCount

Definition at line 78 of file SliceNavigator.h.

Referenced by postProcessEvent(), and SliceNavigator().

int evd::SliceNavigator::fAutoPrintMax

How many events to print (zero = disable printing)

Definition at line 76 of file SliceNavigator.h.

Referenced by postProcessEvent(), and reconfigure().

std::string evd::SliceNavigator::fAutoPrintPattern

Pattern for constructing filename.

Definition at line 75 of file SliceNavigator.h.

Referenced by reconfigure().

std::string evd::SliceNavigator::fCellHitInstance

protected

Definition at line 114 of file SliceNavigator.h.

Referenced by postProcessEvent(), and reconfigure().

std::string evd::SliceNavigator::fCellHitLabel

protected

Definition at line 112 of file SliceNavigator.h.

Referenced by postProcessEvent(), and reconfigure().

int evd::SliceNavigator::fCurrentSlice

protected

Index into fSlicesToNavigate.

Definition at line 116 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 109 of file SliceNavigator.h.

Referenced by DimHits(), and reconfigure().

bool evd::SliceNavigator::fDimOthers

protected

Definition at line 110 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 108 of file SliceNavigator.h.

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

bool evd::SliceNavigator::fFullEvent

protected

Definition at line 115 of file SliceNavigator.h.

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

int evd::SliceNavigator::fInitialSlice

protected

-1 afterwards

Definition at line 123 of file SliceNavigator.h.

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

int evd::SliceNavigator::fMaxSliceId

protected

Definition at line 117 of file SliceNavigator.h.

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

int evd::SliceNavigator::fRun

protected

Definition at line 106 of file SliceNavigator.h.

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

std::string evd::SliceNavigator::fSliceInstance

protected

Definition at line 113 of file SliceNavigator.h.

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

std::string evd::SliceNavigator::fSliceLabel

protected

Definition at line 111 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 120 of file SliceNavigator.h.

Referenced by GetProducts(), and postProcessEvent().

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

protected

Definition at line 118 of file SliceNavigator.h.

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

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-12-03/EventDisplayServices/SliceNavigator.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-12-03/EventDisplayServices/SliceNavigator_service.cc