sbn::TrackStoppingChi2Alg Class Reference

#include <TrackStoppingChi2Alg.h>

Public Member Functions
	TrackStoppingChi2Alg (fhicl::ParameterSet const &p)
StoppingChi2Fit	RunFit (const std::vector< float > &dEdxVec, const std::vector< float > &resRangeVec) const
StoppingChi2Fit	RunFit (const anab::Calorimetry &calo) const
StoppingChi2Fit	RunFitForCosmicID (const anab::Calorimetry &calo) const

Private Attributes
const float	fFitRange
const float	fMaxdEdx
const unsigned int	fMinHits
sbn::TPCGeoAlg	fTpcGeo

Detailed Description

Definition at line 22 of file TrackStoppingChi2Alg.h.

Constructor & Destructor Documentation

sbn::TrackStoppingChi2Alg::TrackStoppingChi2Alg

(

fhicl::ParameterSet const &

p

)

Definition at line 7 of file TrackStoppingChi2Alg.cc.

                                                                        :
  fFitRange(p.get<float>("FitRange"))
  , fMaxdEdx(p.get<float>("MaxdEdx"))
  , fMinHits(p.get<unsigned int>("MinHits"))
{
}

Member Function Documentation

sbn::StoppingChi2Fit sbn::TrackStoppingChi2Alg::RunFit	(	const std::vector< float > &	dEdxVec,
		const std::vector< float > &	resRangeVec
	)		const

Definition at line 14 of file TrackStoppingChi2Alg.cc.

{
  if (dEdxVec.size() != resRangeVec.size())
    throw cet::exception("TrackStoppingChi2Alg") << "dEdx and Res Range do not have same length: " << dEdxVec.size() << " and " << resRangeVec.size() << std::endl;

  if (dEdxVec.size() < fMinHits)
    return StoppingChi2Fit();

  const auto graph(std::make_unique<TGraph>(dEdxVec.size(), &resRangeVec[0], &dEdxVec[0]));

  // Try and fit a flat polynomial
  graph->Fit("pol0", "Q");
  const TF1* polFit = graph->GetFunction("pol0");
  const float pol0Chi2(polFit ? polFit->GetChisquare() : -5.f);
  const float pol0Fit(polFit ? polFit->GetParameter(0) : -5.f);

  // Try to fit an exponential
  graph->Fit("expo", "Q");
  const TF1* expFit = graph->GetFunction("expo");
  const float expChi2(expFit ? expFit->GetChisquare() : -5.f);

  return StoppingChi2Fit(pol0Chi2, expChi2, pol0Fit);
}

sbn::StoppingChi2Fit sbn::TrackStoppingChi2Alg::RunFit

(

const anab::Calorimetry &

calo

)

const

Definition at line 38 of file TrackStoppingChi2Alg.cc.

{

  std::vector<float> dEdxVec, resRangeVec;
  // Fill the dEdx vs res range vectors, ignoring the first/last points
  for (size_t i = 1; i < calo.dEdx().size() - 1; i++) {
    const float thisdEdx(calo.dEdx()[i]);
    const float thisResRange(calo.ResidualRange()[i]);
    if (thisResRange > fFitRange || thisdEdx > fMaxdEdx)
      continue;

    dEdxVec.push_back(thisdEdx);
    resRangeVec.push_back(thisResRange);
  }

  return this->RunFit(dEdxVec, resRangeVec);
}

sbn::StoppingChi2Fit sbn::TrackStoppingChi2Alg::RunFitForCosmicID

(

const anab::Calorimetry &

calo

)

const

Definition at line 56 of file TrackStoppingChi2Alg.cc.

{
  if(calo.XYZ().size() == 0)
    return StoppingChi2Fit();

  geo::Point_t start(calo.XYZ().front());
  geo::Point_t end(calo.XYZ().back());

  std::vector<float> dEdxVec, resRangeVec;
  // Fill the dEdx vs res range vectors, ignoring the first/last points
  for (size_t i = 1; i < calo.dEdx().size() - 1; i++) {
    const float thisdEdx(calo.dEdx()[i]);
    const float thisResRange(calo.ResidualRange()[i]);
    if (thisResRange > fFitRange || thisdEdx > fMaxdEdx)
      continue;

    dEdxVec.push_back(thisdEdx);
    resRangeVec.push_back(thisResRange);
  }

  if(fTpcGeo.MinDistToWall(start) > fTpcGeo.MinDistToWall(end))
    {
      std::reverse(dEdxVec.begin(), dEdxVec.end());
      std::reverse(resRangeVec.begin(), resRangeVec.end());
    }

  return this->RunFit(dEdxVec, resRangeVec);
}

Member Data Documentation

const float sbn::TrackStoppingChi2Alg::fFitRange

private

Definition at line 38 of file TrackStoppingChi2Alg.h.

const float sbn::TrackStoppingChi2Alg::fMaxdEdx

private

Definition at line 38 of file TrackStoppingChi2Alg.h.

const unsigned int sbn::TrackStoppingChi2Alg::fMinHits

private

Definition at line 39 of file TrackStoppingChi2Alg.h.

sbn::TPCGeoAlg sbn::TrackStoppingChi2Alg::fTpcGeo

private

Definition at line 41 of file TrackStoppingChi2Alg.h.

---

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

• TrackStoppingChi2Alg.h
• TrackStoppingChi2Alg.cc