sbn::RangePAllPID Class Reference

Inheritance diagram for sbn::RangePAllPID:

Public Member Functions
	RangePAllPID (fhicl::ParameterSet const &p)
	RangePAllPID (RangePAllPID const &)=delete
	RangePAllPID (RangePAllPID &&)=delete
RangePAllPID &	operator= (RangePAllPID const &)=delete
RangePAllPID &	operator= (RangePAllPID &&)=delete
void	produce (art::Event &e) override

Private Attributes
trkf::TrackMomentumCalculator	fRangeCalculator
art::InputTag	fTrackLabel
double	fMuonMass
double	fPionMass

Detailed Description

Definition at line 33 of file RangePAllPID_module.cc.

Constructor & Destructor Documentation

sbn::RangePAllPID::RangePAllPID ( fhicl::ParameterSet const &  p )

explicit

Definition at line 60 of file RangePAllPID_module.cc.

  : EDProducer{p},
    fRangeCalculator(p.get<float>("MinTrackLength", 10.)),
    fTrackLabel(p.get<art::InputTag>("TrackLabel", "pandoraTrack"))
{
  for (unsigned i = 0; i < names.size(); i++) {
    produces<std::vector<sbn::RangeP>>(names[i]);
    produces<art::Assns<recob::Track, sbn::RangeP>>(names[i]); 
  }

  const TDatabasePDG *PDGTable = TDatabasePDG::Instance();
  fMuonMass = PDGTable->GetParticle(13)->Mass();
  fPionMass = PDGTable->GetParticle(211)->Mass();
}

sbn::RangePAllPID::RangePAllPID ( RangePAllPID const &  )

delete

sbn::RangePAllPID::RangePAllPID ( RangePAllPID &&  )

delete

Member Function Documentation

RangePAllPID& sbn::RangePAllPID::operator= ( RangePAllPID const &  )

delete

RangePAllPID& sbn::RangePAllPID::operator= ( RangePAllPID &&  )

delete

void sbn::RangePAllPID::produce ( art::Event &  e )

override

Definition at line 75 of file RangePAllPID_module.cc.

{
  art::Handle<std::vector<recob::Track>> track_handle;
  e.getByLabel(fTrackLabel, track_handle);

  std::vector<art::Ptr<recob::Track>> tracks;
  art::fill_ptr_vector(tracks, track_handle);

  for (unsigned i = 0; i < PIDs.size(); i++) {
    std::unique_ptr<std::vector<sbn::RangeP>> rangecol(new std::vector<sbn::RangeP>);
    std::unique_ptr<art::Assns<recob::Track, sbn::RangeP>> assn(new art::Assns<recob::Track, sbn::RangeP>);

    for (const art::Ptr<recob::Track> track: tracks) {
      sbn::RangeP rangep;
      // Rescale the input and output as described by https://inspirehep.net/literature/1766384 (eq. 6.2)
      if (PIDs[i] == 211) {
        rangep.range_p = fRangeCalculator.GetTrackMomentum(track->Length() * (fMuonMass / fPionMass), 13) * (fPionMass / fMuonMass); 
      }
      else {
        rangep.range_p = fRangeCalculator.GetTrackMomentum(track->Length(), PIDs[i]);
      }
      rangep.trackID = track->ID();
      rangecol->push_back(rangep);
      util::CreateAssn(*this, e, *rangecol, track, *assn, names[i]);
    }

    e.put(std::move(rangecol), names[i]);
    e.put(std::move(assn), names[i]);
  }

}

Member Data Documentation

double sbn::RangePAllPID::fMuonMass

private

Definition at line 52 of file RangePAllPID_module.cc.

double sbn::RangePAllPID::fPionMass

private

Definition at line 53 of file RangePAllPID_module.cc.

trkf::TrackMomentumCalculator sbn::RangePAllPID::fRangeCalculator

private

Definition at line 50 of file RangePAllPID_module.cc.

art::InputTag sbn::RangePAllPID::fTrackLabel

private

Definition at line 51 of file RangePAllPID_module.cc.

---

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

• RangePAllPID_module.cc