filt::GenNuFilter Class Reference

Inheritance diagram for filt::GenNuFilter:

Public Member Functions
	GenNuFilter (fhicl::ParameterSet const &pset)
virtual bool	filter (art::Event &e) override
void	reconfigure (fhicl::ParameterSet const &pset)
virtual void	beginJob () override

Private Attributes
bool	fVtxInTPC
std::vector< int >	fLepPDGs
bool	fCC
bool	fNC

Detailed Description

Definition at line 15 of file GenNuFilter_module.cc.

Constructor & Destructor Documentation

filt::GenNuFilter::GenNuFilter ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 32 of file GenNuFilter_module.cc.

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

Member Function Documentation

void filt::GenNuFilter::beginJob ( )

overridevirtual

Definition at line 85 of file GenNuFilter_module.cc.

                             {
  }

bool filt::GenNuFilter::filter ( art::Event &  e )

overridevirtual

Definition at line 48 of file GenNuFilter_module.cc.

                                      {
    //std::vector< art::Handle< std::vector<simb::MCTruth> > > mclists;
    //e.getManyByType(mclists);
    auto mclists = e.getMany< std::vector<simb::MCTruth> >();
    for (unsigned int i = 0; i < mclists.size() ; i++){
      for (unsigned int j = 0; j < mclists[i]->size(); j++){
        //Should have the truth record for the event now
        const art::Ptr<simb::MCTruth> mc_truth(mclists[i],j);
        if(mc_truth->Origin() == simb::kBeamNeutrino){
          simb::MCNeutrino nu = mc_truth->GetNeutrino();
          // Check vertex in TPC
          double vtxX = nu.Nu().Vx();
          double vtxY = nu.Nu().Vy();
          double vtxZ = nu.Nu().Vz();
          TVector3 vtx(vtxX, vtxY, vtxZ);
          if(fVtxInTPC && !RecoUtils::IsInsideTPC(vtx, 0)) continue;
          // Check if CC or NC specified
          if(!(fCC && fNC)){
            if(fCC && nu.CCNC() == simb::kNC) continue;
            if(fNC && nu.CCNC() == simb::kCC) continue;
          }
          // Check which flavour of interaction
          int lepPdg = nu.Lepton().PdgCode();
          if(fLepPDGs.size() > 0){
            if(fLepPDGs[0] != 0){
              if(std::find(fLepPDGs.begin(), fLepPDGs.end(), lepPdg) == fLepPDGs.end()) continue;
            }
          }
          return true;
        }
      }
    }

    return false;
  }

void filt::GenNuFilter::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 40 of file GenNuFilter_module.cc.

                                                            {
    fVtxInTPC = pset.get<bool>("VtxInTPC");
    fLepPDGs = pset.get<std::vector<int>>("LepPDGs");
    fCC = pset.get<bool>("CC");
    fNC = pset.get<bool>("NC");
  }

Member Data Documentation

bool filt::GenNuFilter::fCC

private

Definition at line 26 of file GenNuFilter_module.cc.

std::vector<int> filt::GenNuFilter::fLepPDGs

private

Definition at line 25 of file GenNuFilter_module.cc.

bool filt::GenNuFilter::fNC

private

Definition at line 27 of file GenNuFilter_module.cc.

bool filt::GenNuFilter::fVtxInTPC

private

Definition at line 24 of file GenNuFilter_module.cc.

---

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

• GenNuFilter_module.cc