FluxGeoFilter Class Reference

Inheritance diagram for FluxGeoFilter:

Public Member Functions
	FluxGeoFilter (fhicl::ParameterSet const &p)
	FluxGeoFilter (FluxGeoFilter const &)=delete
	FluxGeoFilter (FluxGeoFilter &&)=delete
FluxGeoFilter &	operator= (FluxGeoFilter const &)=delete
FluxGeoFilter &	operator= (FluxGeoFilter &&)=delete
bool	filter (art::Event &e) override

Private Attributes
std::set< std::string >	fVolFlux

Detailed Description

Definition at line 35 of file FluxGeoFilter_module.cc.

Constructor & Destructor Documentation

FluxGeoFilter::FluxGeoFilter ( fhicl::ParameterSet const &  p )

explicit

Definition at line 55 of file FluxGeoFilter_module.cc.

  : EDFilter{p}  // ,
  // More initializers here.
{
  std::cout<<"Configuring flux filter."<<std::endl;
  std::vector<std::string> tmp=p.get<std::vector<std::string> >("volumes");
  fVolFlux.insert(tmp.begin(),tmp.end());

  std::cout<<"Filtering flux through volumes: "<<std::endl;
  for (auto s: fVolFlux) std::cout<<"\t"<<s<<std::endl;
}

FluxGeoFilter::FluxGeoFilter ( FluxGeoFilter const &  )

delete

FluxGeoFilter::FluxGeoFilter ( FluxGeoFilter &&  )

delete

Member Function Documentation

bool FluxGeoFilter::filter ( art::Event &  e )

override

Definition at line 67 of file FluxGeoFilter_module.cc.

{

  // Implementation of required member function here.
  bool result=false; //false filters out, true passes

  //check if neutrino goes through volTPCActive
  geo::GeometryCore const* geom = lar::providerFrom<geo::Geometry>();
  TGeoManager* rgeo=geom->ROOTGeoManager();

  // art::Handle< std::vector<simb::MCFlux> > mcFluxHandle;
  // e.getByLabel("flux",mcFluxHandle);
  // std::vector<simb::MCFlux> const& fluxlist = *mcFluxHandle;

  art::Handle< std::vector<simb::MCTruth> > mctruthHandle;
  e.getByLabel("flux",mctruthHandle);
  std::vector<simb::MCTruth> const& mclist = *mctruthHandle;

  for(unsigned int inu = 0; inu < mclist.size(); inu++){
    simb::MCParticle nu = mclist[inu].GetNeutrino().Nu();
    rgeo->SetCurrentPoint(nu.Vx(),nu.Vy(),nu.Vz());
    rgeo->SetCurrentDirection(nu.Px(),nu.Py(),nu.Pz());
    TGeoNode* node=rgeo->FindNode();
    while (node) {
      std::string volname=node->GetVolume()->GetName();
      rgeo->FindNextBoundary();
      node=gGeoManager->Step();
      if (fVolFlux.find(volname)!=fVolFlux.end()) {
        result=true;
        break;
      }
    }
  }
  return result;
}

FluxGeoFilter& FluxGeoFilter::operator= ( FluxGeoFilter const &  )

delete

FluxGeoFilter& FluxGeoFilter::operator= ( FluxGeoFilter &&  )

delete

Member Data Documentation

std::set<std::string> FluxGeoFilter::fVolFlux

private

Definition at line 52 of file FluxGeoFilter_module.cc.

---

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

• FluxGeoFilter_module.cc