protoana::ProtoDUNEBeamCuts Class Reference

#include <ProtoDUNEBeamCuts.h>

Public Member Functions
	ProtoDUNEBeamCuts ()
	ProtoDUNEBeamCuts (const fhicl::ParameterSet &pset)
bool	IsBeamlike (const recob::Track &, const art::Event &, std::string)
bool	IsBeamlike (const recob::Shower &, const art::Event &, std::string)

Private Member Functions
BeamVals	GetMCBeam (const art::Event &)
BeamVals	GetDataBeam (const art::Event &)

Private Attributes
std::map< std::string, fhicl::ParameterSet >	DataCuts
std::map< std::string, fhicl::ParameterSet >	MCCuts
std::vector< std::string >	valid_momenta

Detailed Description

Definition at line 22 of file ProtoDUNEBeamCuts.h.

Constructor & Destructor Documentation

protoana::ProtoDUNEBeamCuts::ProtoDUNEBeamCuts ( )

inline

Definition at line 25 of file ProtoDUNEBeamCuts.h.

{};

protoana::ProtoDUNEBeamCuts::ProtoDUNEBeamCuts

(

const fhicl::ParameterSet &

pset

)

Definition at line 5 of file ProtoDUNEBeamCuts.cxx.

                                                                            {


  std::vector< fhicl::ParameterSet > DataCuts_sets = pset.get< std::vector< fhicl::ParameterSet > >
                                                             ("DataCuts");
  std::vector< fhicl::ParameterSet > MCCuts_sets   = pset.get< std::vector< fhicl::ParameterSet > >
                                                             ("MCCuts");

  //Go through the data cuts 
  for( size_t i = 0; i < DataCuts_sets.size(); ++i ){
    std::string momentum = DataCuts_sets[i].get<std::string>("Momentum");
    if( std::find( valid_momenta.begin(), valid_momenta.end(), momentum ) == valid_momenta.end() ){
      std::exception e;
      throw e;
    }

    if( DataCuts.find( momentum ) != DataCuts.end() ){
      std::cerr << "Error: Attempting to duplicate a cut parameter set" << std::endl;
      std::exception e;
      throw e;
    }
     
    DataCuts[momentum] = DataCuts_sets[i];
  }


  //Now go through the MC cuts
  for( size_t i = 0; i < MCCuts_sets.size(); ++i ){
    
    std::string momentum = DataCuts_sets[i].get<std::string>("Momentum");
    if( std::find( valid_momenta.begin(), valid_momenta.end(), momentum ) == valid_momenta.end() ){
      std::exception e;
      throw e;
    }

    if( MCCuts.find( momentum ) != MCCuts.end() ){
      std::cerr << "Error: Attempting to duplicate a cut parameter set" << std::endl;
      std::exception e;
      throw e;
    }
     
    MCCuts[momentum] = MCCuts_sets[i];
  }


}

Member Function Documentation

protoana::BeamVals protoana::ProtoDUNEBeamCuts::GetDataBeam ( const art::Event &  evt )

private

Definition at line 206 of file ProtoDUNEBeamCuts.cxx.

                                                                             {
  std::vector<art::Ptr<beam::ProtoDUNEBeamEvent>> beamVec;
  auto beamHandle = evt.getValidHandle< std::vector< beam::ProtoDUNEBeamEvent > >("beamevent");
                        
  BeamVals result;
  result.Valid = false;

  if( beamHandle.isValid()){
    art::fill_ptr_vector(beamVec, beamHandle);
  }
  //Should just have one
  const beam::ProtoDUNEBeamEvent & beamEvent = *(beamVec.at(0));

  const std::vector< recob::Track > & beamTracks = beamEvent.GetBeamTracks();
  if( beamTracks.size() == 0 ){
    std::cout << "Warning: no tracks associated to beam data" << std::endl;
    return result;  
  }
  else if( beamTracks.size() > 1 ){
    std::cout << "Warning: mutiple tracks associated to beam data" << std::endl;
    return result;  
  }
  
  result.Valid = true;
  result.X = beamTracks.at(0).Trajectory().End().X();
  result.Y = beamTracks.at(0).Trajectory().End().Y();

  result.DirX = beamTracks.at(0).EndDirection().X();
  result.DirY = beamTracks.at(0).EndDirection().Y();
  result.DirZ = beamTracks.at(0).EndDirection().Z();
  return result;
}

protoana::BeamVals protoana::ProtoDUNEBeamCuts::GetMCBeam ( const art::Event &  evt )

private

Definition at line 181 of file ProtoDUNEBeamCuts.cxx.

                                                                           {

  protoana::ProtoDUNETruthUtils truthUtil;
  auto mcTruths = evt.getValidHandle< std::vector< simb::MCTruth > >("generator");
  const simb::MCParticle* true_beam_particle = truthUtil.GetGeantGoodParticle((*mcTruths)[0],evt);

  BeamVals result;
  result.Valid = false;

  if( !true_beam_particle ){
    std::cout << "No true beam particle" << std::endl;       
    return result;
  }

  result.Valid = true;
  result.DirX = true_beam_particle->Px() / true_beam_particle->P();
  result.DirY = true_beam_particle->Py() / true_beam_particle->P();
  result.DirZ = true_beam_particle->Pz() / true_beam_particle->P();

  //Project the beam to Z = 0
  result.X = true_beam_particle->Position(0).X() + (-1.*true_beam_particle->Position(0).Z())*(result.DirX / result.DirZ);
  result.Y = true_beam_particle->Position(0).Y() + (-1.*true_beam_particle->Position(0).Z())*(result.DirY / result.DirZ);
  return result;
}

bool protoana::ProtoDUNEBeamCuts::IsBeamlike	(	const recob::Track &	,
		const art::Event &	,
		std::string
	)

bool protoana::ProtoDUNEBeamCuts::IsBeamlike	(	const recob::Shower &	,
		const art::Event &	,
		std::string
	)

Member Data Documentation

std::map< std::string, fhicl::ParameterSet > protoana::ProtoDUNEBeamCuts::DataCuts

private

Definition at line 35 of file ProtoDUNEBeamCuts.h.

std::map< std::string, fhicl::ParameterSet > protoana::ProtoDUNEBeamCuts::MCCuts

private

Definition at line 36 of file ProtoDUNEBeamCuts.h.

std::vector< std::string > protoana::ProtoDUNEBeamCuts::valid_momenta

private

Initial value:

= {
        ".3", ".5", "1", "2", "3", "6", "7"
      }

Definition at line 38 of file ProtoDUNEBeamCuts.h.

---

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

• protoduneana/protoduneana/Utilities/ProtoDUNEBeamCuts.h
• protoduneana/protoduneana/Utilities/ProtoDUNEBeamCuts.cxx