MCReco Class Reference

Inheritance diagram for MCReco:

Public Member Functions
	MCReco (fhicl::ParameterSet const &p)
void	produce (art::Event &e) override
template<typename T >
void	MakeMCEdep (art::Event &evt)

Private Attributes
art::InputTag	fMCParticleLabel
art::InputTag	fMCParticleLiteLabel
art::InputTag	fSimChannelLabel
bool	fUseSimEnergyDeposit
bool	fUseSimEnergyDepositLite
bool	fIncludeDroppedParticles
::sim::MCRecoPart	fPart
::sim::MCRecoEdep	fEdep
::sim::MCShowerRecoAlg	fMCSAlg
::sim::MCTrackRecoAlg	fMCTAlg

Detailed Description

Definition at line 28 of file MCReco_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 52 of file MCReco_module.cc.

  : EDProducer{pset}
  , fPart   (pset.get< fhicl::ParameterSet >("MCRecoPart"))
  , fEdep   (pset.get< fhicl::ParameterSet >("MCRecoEdep"))
  , fMCSAlg (pset.get< fhicl::ParameterSet >("MCShowerRecoAlg"))
  , fMCTAlg (pset.get< fhicl::ParameterSet >("MCTrackRecoAlg"))
{

  //for backwards compatibility to using the "G4ModName" label...
  if(!(pset.get_if_present<art::InputTag>("MCParticleLabel",fMCParticleLabel) &&
       pset.get_if_present<art::InputTag>("SimChannelLabel",fSimChannelLabel)) ){

    mf::LogWarning("MCReco_module") << "USING DEPRECATED G4ModName CONFIG IN MCRECO_MODULE"
            << "\nUse 'MCParticleLabel' and 'SimChannelLabel' instead.";

    fMCParticleLabel = pset.get<art::InputTag>("G4ModName","largeant");
    fMCParticleLiteLabel = pset.get<art::InputTag>("G4ModName","largeant");
    fSimChannelLabel = pset.get<art::InputTag>("G4ModName","largeant");
  }
  else {
    fMCParticleLiteLabel = pset.get<art::InputTag>("MCParticleLiteLabel", "largeant");
  }

  fUseSimEnergyDeposit = pset.get<bool>("UseSimEnergyDeposit",false);
  fUseSimEnergyDepositLite = pset.get<bool>("UseSimEnergyDepositLite",false);
  fIncludeDroppedParticles = pset.get<bool>("IncludeDroppedParticles",false);

  if (fUseSimEnergyDepositLite && fUseSimEnergyDeposit) {
    mf::LogWarning("MCReco_module") << "Asked to use both SimEnergyDeposit and SimEnergyDepositLite - will use SimEnergyDeposit.";
  }

  produces< std::vector< sim::MCShower> >();
  produces< std::vector< sim::MCTrack>  >();
  // Call appropriate produces<>() functions here.

//MCReco::~MCReco()
//{
  // Clean up dynamic memory and other resources here.
//}
}

Member Function Documentation

template<typename T >

void MCReco::MakeMCEdep

(

art::Event &

evt

)

Definition at line 142 of file MCReco_module.cc.

                                                           {
  // Retrieve T
  auto const& sed_array = *evt.getValidHandle<std::vector<T>>(fSimChannelLabel);
  fEdep.MakeMCEdep(sed_array);
}

void MCReco::produce ( art::Event &  e )

override

Definition at line 93 of file MCReco_module.cc.

{
//  std::unique_ptr< std::vector<sim::MCTrack> > outTrackArray(new std::vector<sim::MCTrack>);

  // Retrieve mcparticles
  art::Handle<std::vector<simb::MCParticle> > mcpHandle;
  evt.getByLabel(fMCParticleLabel,mcpHandle);
  if(!mcpHandle.isValid()) throw cet::exception(__FUNCTION__) << "Failed to retrieve simb::MCParticle";;

  // Find associations
  art::FindOneP<simb::MCTruth> ass(mcpHandle, evt, fMCParticleLabel);
  std::vector<simb::Origin_t> orig_array;
  orig_array.reserve(mcpHandle->size());
  for(size_t i=0; i<mcpHandle->size(); ++i) {
    const art::Ptr<simb::MCTruth> &mct = ass.at(i);
    orig_array.push_back(mct->Origin());
  }

  const std::vector<simb::MCParticle>& mcp_array(*mcpHandle);

  if (fIncludeDroppedParticles) {
    auto const& mcmp_array = *evt.getValidHandle<std::vector<sim::MCParticleLite>>(fMCParticleLiteLabel);
    fPart.AddParticles(mcp_array, orig_array, mcmp_array);
  } // end if fIncludeDroppedParticles
  else {
    fPart.AddParticles(mcp_array,orig_array);
  }

  // change implemented by David Caratelli to allow for MCRECO to run without SimChannels and using
  // SimEnergyDeposits instead
  if (fUseSimEnergyDeposit == true) {
    MakeMCEdep<sim::SimEnergyDeposit>(evt);
  }
  // change implemented by Laura Domine to allow for MCRECO to run with SimEnergyDepositLite
  else if (fUseSimEnergyDepositLite == true) {
    MakeMCEdep<sim::SimEnergyDepositLite>(evt);
  }
  else {
    MakeMCEdep<sim::SimChannel>(evt);
  }

  //Add MCShowers and MCTracks to the event
  evt.put(fMCSAlg.Reconstruct(fPart,fEdep));
  evt.put(fMCTAlg.Reconstruct(fPart,fEdep));

  fEdep.Clear();
  fPart.clear();
}

Member Data Documentation

::sim::MCRecoEdep MCReco::fEdep

private

Definition at line 47 of file MCReco_module.cc.

bool MCReco::fIncludeDroppedParticles

private

Definition at line 44 of file MCReco_module.cc.

art::InputTag MCReco::fMCParticleLabel

private

Definition at line 39 of file MCReco_module.cc.

art::InputTag MCReco::fMCParticleLiteLabel

private

Definition at line 40 of file MCReco_module.cc.

::sim::MCShowerRecoAlg MCReco::fMCSAlg

private

Definition at line 48 of file MCReco_module.cc.

::sim::MCTrackRecoAlg MCReco::fMCTAlg

private

Definition at line 49 of file MCReco_module.cc.

::sim::MCRecoPart MCReco::fPart

private

Definition at line 46 of file MCReco_module.cc.

art::InputTag MCReco::fSimChannelLabel

private

Definition at line 41 of file MCReco_module.cc.

bool MCReco::fUseSimEnergyDeposit

private

Definition at line 42 of file MCReco_module.cc.

bool MCReco::fUseSimEnergyDepositLite

private

Definition at line 43 of file MCReco_module.cc.

---

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

• MCReco_module.cc