mvapid::MVAAlg Class Reference

#include "MVAAlg.h"

Classes
struct	SortedObj
struct	SumDistance2

Public Member Functions
	MVAAlg (fhicl::ParameterSet const &pset, const art::EDProducer *parentModule)
virtual	~MVAAlg ()
void	GetDetectorEdges ()
void	GetWireNormals ()
void	reconfigure (fhicl::ParameterSet const &p)
void	RunPID (art::Event &evt, std::vector< anab::MVAPIDResult > &result, art::Assns< recob::Track, anab::MVAPIDResult, void > &trackAssns, art::Assns< recob::Shower, anab::MVAPIDResult, void > &showerAssns)

Private Member Functions
int	IsInActiveVol (const TVector3 &pos)
void	PrepareEvent (const art::Event &event)
void	FitAndSortTrack (art::Ptr< recob::Track > track, int &isStoppingReco, SortedObj &sortedObj)
void	SortShower (art::Ptr< recob::Shower > shower, int &isStoppingReco, mvapid::MVAAlg::SortedObj &sortedShower)
void	RunPCA (std::vector< art::Ptr< recob::Hit > > &hits, std::vector< double > &eVals, std::vector< double > &eVecs)
void	_Var_Shape (const SortedObj &track, double &coreHaloRatio, double &concentration, double &conicalness)
double	CalcSegmentdEdxFrac (const SortedObj &track, double start, double end)
double	CalcSegmentdEdxDist (const SortedObj &track, double start, double end)
double	CalcSegmentdEdxDistAtEnd (const mvapid::MVAAlg::SortedObj &track, double distAtEnd)
int	LinFit (const art::Ptr< recob::Track > track, TVector3 &trackPoint, TVector3 &trackDir)
int	LinFitShower (const art::Ptr< recob::Shower > shower, TVector3 &showerPoint, TVector3 &showerDir)

Private Attributes
const calo::CalorimetryAlg	fCaloAlg
double	fEventT0
double	fDetMinX
double	fDetMaxX
double	fDetMinY
double	fDetMaxY
double	fDetMinZ
double	fDetMaxZ
std::map< int, double >	fNormToWiresY
std::map< int, double >	fNormToWiresZ
const art::EDProducer *	fParentModule
std::string	fTrackLabel
std::string	fShowerLabel
std::string	fHitLabel
std::string	fSpacePointLabel
std::string	fTrackingLabel
std::vector< art::Ptr< recob::Track > >	fTracks
std::vector< art::Ptr< recob::Shower > >	fShowers
std::vector< art::Ptr< recob::SpacePoint > >	fSpacePoints
std::vector< art::Ptr< recob::Hit > >	fHits
std::map< art::Ptr< recob::Track >, std::vector< art::Ptr< recob::Hit > > >	fTracksToHits
std::map< art::Ptr< recob::Track >, std::vector< art::Ptr< recob::SpacePoint > > >	fTracksToSpacePoints
std::map< art::Ptr< recob::Shower >, std::vector< art::Ptr< recob::Hit > > >	fShowersToHits
std::map< art::Ptr< recob::Shower >, std::vector< art::Ptr< recob::SpacePoint > > >	fShowersToSpacePoints
std::map< art::Ptr< recob::Hit >, art::Ptr< recob::SpacePoint > >	fHitsToSpacePoints
std::map< art::Ptr< recob::SpacePoint >, art::Ptr< recob::Hit > >	fSpacePointsToHits
anab::MVAPIDResult	fResHolder
TMVA::Reader	fReader
std::vector< std::string >	fMVAMethods
std::vector< std::string >	fWeightFiles
bool	fCheatVertex
TLorentzVector	fVertex4Vect

Detailed Description

Definition at line 47 of file MVAAlg.h.

Constructor & Destructor Documentation

mvapid::MVAAlg::MVAAlg	(	fhicl::ParameterSet const &	pset,
		const art::EDProducer *	parentModule
	)

Definition at line 28 of file MVAAlg.cxx.

References anab::MVAPIDResult::concentration, anab::MVAPIDResult::conicalness, anab::MVAPIDResult::coreHaloRatio, anab::MVAPIDResult::dEdxEnd, anab::MVAPIDResult::dEdxEndRatio, anab::MVAPIDResult::dEdxStart, anab::MVAPIDResult::evalRatio, fCheatVertex, fHitLabel, fMVAMethods, fReader, fResHolder, fShowerLabel, fSpacePointLabel, fTrackingLabel, fTrackLabel, fWeightFiles, and fhicl::ParameterSet::get().

Referenced by mvapid::MVAAlg::SumDistance2::operator()().

                                                                                    :
  fCaloAlg(pset.get<fhicl::ParameterSet>("CalorimetryAlg")), 
  fParentModule(parentModule),fReader(""){
  fHitLabel=pset.get<std::string>("HitLabel");
  fTrackLabel=pset.get<std::string>("TrackLabel");
  fShowerLabel=pset.get<std::string>("ShowerLabel");
  fSpacePointLabel=pset.get<std::string>("SpacePointLabel");
  fTrackingLabel=pset.get<std::string>("TrackingLabel","");

  fCheatVertex=pset.get<bool>("CheatVertex",false);

  fReader.AddVariable("evalRatio",&fResHolder.evalRatio);
  fReader.AddVariable("coreHaloRatio",&fResHolder.coreHaloRatio);
  fReader.AddVariable("concentration",&fResHolder.concentration);
  fReader.AddVariable("conicalness",&fResHolder.conicalness);
  fReader.AddVariable("dEdxStart",&fResHolder.dEdxStart);
  fReader.AddVariable("dEdxEnd",&fResHolder.dEdxEnd);
  fReader.AddVariable("dEdxEndRatio",&fResHolder.dEdxEndRatio);

  fMVAMethods=pset.get<std::vector<std::string> >("MVAMethods");
  std::vector<std::string> weightFileBnames=pset.get<std::vector<std::string> >("WeightFiles");

  cet::search_path searchPath("FW_SEARCH_PATH");
  for(auto fileIter=weightFileBnames.begin();fileIter!=weightFileBnames.end();++fileIter){
    std::string fileWithPath;
    if(!searchPath.find_file(*fileIter,fileWithPath)){
        fWeightFiles.clear();
        fMVAMethods.clear();
        throw cet::exception("MVAPID")<<"Unable to find weight file "<<*fileIter<<" in search path "
                                      <<searchPath.to_string()<<std::endl;
      }
      fWeightFiles.push_back(fileWithPath);
  }

  if(fMVAMethods.size()!=fWeightFiles.size()){
    std::cerr<<"Mismatch in number of MVA methods and weight files!"<<std::endl;
    exit(1);
  }

  for(unsigned int iMethod=0;iMethod!=fMVAMethods.size();++iMethod){
    fReader.BookMVA(fMVAMethods[iMethod], fWeightFiles[iMethod]);
  }
}

mvapid::MVAAlg::~MVAAlg ( )

virtual

Definition at line 72 of file MVAAlg.cxx.

Referenced by mvapid::MVAAlg::SumDistance2::operator()().

{}

Member Function Documentation

void mvapid::MVAAlg::_Var_Shape ( const SortedObj &  track, double &  coreHaloRatio, double &  concentration, double &  conicalness  )

private

Definition at line 539 of file MVAAlg.cxx.

References mvapid::MVAAlg::SortedObj::dir, E, fHitsToSpacePoints, mvapid::MVAAlg::SortedObj::hitMap, mvapid::MVAAlg::SortedObj::length, max, mvapid::MVAAlg::SortedObj::start, and recob::SpacePoint::XYZ().

Referenced by mvapid::MVAAlg::SumDistance2::operator()(), and RunPID().

                                                    {
  
  static const unsigned int conMinHits=3;
  static const double minCharge = 0.1;
  static const double conFracRange=0.2;
  static const double MoliereRadius = 10.1;
  static const double MoliereRadiusFraction = 0.2;


  double totalCharge=0;
  double totalChargeStart=0;
  double totalChargeEnd=0;

  double chargeCore=0;
  double chargeHalo=0;
  double chargeCon=0;
  unsigned int nHits = 0;

  //stuff for conicalness
  double chargeConStart=0;
  double chargeConEnd=0;
  unsigned int nHitsConStart=0;
  unsigned int nHitsConEnd=0;
  

  for(auto hitIter=track.hitMap.begin();hitIter!=track.hitMap.end();++hitIter){
    if(fHitsToSpacePoints.count(hitIter->second)) {
      art::Ptr<recob::SpacePoint> sp=fHitsToSpacePoints.at(hitIter->second);

      double distFromTrackFit = ((TVector3(sp->XYZ()) - track.start).Cross(track.dir)).Mag();
 
      ++nHits;

      if(distFromTrackFit < MoliereRadiusFraction * MoliereRadius)
        chargeCore += hitIter->second->Integral();
      else
        chargeHalo += hitIter->second->Integral();

      totalCharge += hitIter->second->Integral();

      chargeCon += hitIter->second->Integral() / std::max(1.E-2,distFromTrackFit);
      if(hitIter->first/track.length<conFracRange){
        chargeConStart+=distFromTrackFit*distFromTrackFit*hitIter->second->Integral();
        ++nHitsConStart;
        totalChargeStart+=hitIter->second->Integral();
      }
      else if(1.-hitIter->first/track.length<conFracRange){
        chargeConEnd+=distFromTrackFit*distFromTrackFit*hitIter->second->Integral();
        ++nHitsConEnd;
        totalChargeEnd+=hitIter->second->Integral();
      }
    }
  }
  
  coreHaloRatio=chargeHalo/TMath::Max(1.0E-3, chargeCore);
  coreHaloRatio=TMath::Min(100.0, coreHaloRatio);
  concentration=chargeCon/totalCharge;
  if(nHitsConStart>=conMinHits&&nHitsConEnd>=conMinHits&&totalChargeEnd>minCharge&&sqrt(chargeConStart)>minCharge&&totalChargeStart>minCharge){
    conicalness=(sqrt(chargeConEnd)/totalChargeEnd) / (sqrt(chargeConStart)/totalChargeStart);
  }
  else{
    conicalness=1.;
  }
}

double mvapid::MVAAlg::CalcSegmentdEdxDist ( const SortedObj &  track, double  start, double  end  )

private

Definition at line 618 of file MVAAlg.cxx.

References geo::CryostatID::Cryostat, calo::CalorimetryAlg::dEdx_AREA(), dir, mvapid::MVAAlg::SortedObj::dir, fCaloAlg, fEventT0, fNormToWiresY, fNormToWiresZ, mvapid::MVAAlg::SortedObj::hitMap, geo::GeometryCore::Nplanes(), geo::GeometryCore::NTPC(), geo::PlaneID::Plane, geo::TPCID::TPC, recob::Hit::WireID(), and geo::GeometryCore::WirePitch().

Referenced by CalcSegmentdEdxDistAtEnd(), CalcSegmentdEdxFrac(), and mvapid::MVAAlg::SumDistance2::operator()().

                                                                                                  {

  art::ServiceHandle<geo::Geometry> geom;
  
  double totaldEdx=0;
  unsigned int nHits=0;

  //Loop over hits again to calculate average dE/dx and shape variables
  for ( auto hitIter = track.hitMap.begin(); hitIter!=track.hitMap.end(); ++hitIter ){

    if(hitIter->first<start) continue;
    if(hitIter->first>=end) break;

    art::Ptr<recob::Hit> hit=hitIter->second;

    //Pitch to use in dEdx calculation
    double yzPitch = geom->WirePitch(hit->WireID().Plane, hit->WireID().TPC); //pitch not taking into account angle of track or shower
    double xComponent, pitch3D;

    TVector3 dir=track.dir;

    //This assumes equal numbers of TPCs in each cryostat and equal numbers of planes in each TPC
    int planeKey = hit->WireID().Cryostat * geom->NTPC(0) * geom->Nplanes(0, 0) + hit->WireID().TPC * geom->Nplanes(0, 0) + hit->WireID().Plane;

    if(fNormToWiresY.count(planeKey) && fNormToWiresZ.count(planeKey))
      {
        TVector3 normToWires(0.0, fNormToWiresY.at(planeKey), fNormToWiresZ.at(planeKey));
        yzPitch = geom->WirePitch(hit->WireID().Plane, hit->WireID().TPC) / fabs(dir.Dot(normToWires));
      }

    xComponent = yzPitch * dir[0] / sqrt(dir[1] * dir[1] + dir[2] * dir[2]);
    pitch3D = sqrt(xComponent * xComponent + yzPitch * yzPitch);

    double dEdx=fCaloAlg.dEdx_AREA(*hit, pitch3D, fEventT0);
    if( dEdx < 50.){
      ++nHits;
      totaldEdx += dEdx;
    }
  }

  return nHits?totaldEdx/nHits:0;
}

double mvapid::MVAAlg::CalcSegmentdEdxDistAtEnd ( const mvapid::MVAAlg::SortedObj &  track, double  distAtEnd  )

private

Definition at line 612 of file MVAAlg.cxx.

References CalcSegmentdEdxDist(), mvapid::MVAAlg::SortedObj::end, and mvapid::MVAAlg::SortedObj::start.

Referenced by mvapid::MVAAlg::SumDistance2::operator()().

                                                                                                {

  double trackLength=(track.end-track.start).Mag();
  return CalcSegmentdEdxDist(track,trackLength-distAtEnd,trackLength);
}

double mvapid::MVAAlg::CalcSegmentdEdxFrac ( const SortedObj &  track, double  start, double  end  )

private

Definition at line 606 of file MVAAlg.cxx.

References CalcSegmentdEdxDist(), mvapid::MVAAlg::SortedObj::end, and mvapid::MVAAlg::SortedObj::start.

Referenced by mvapid::MVAAlg::SumDistance2::operator()(), and RunPID().

                                                                                                  {

  double trackLength=(track.end-track.start).Mag();
  return CalcSegmentdEdxDist(track,start*trackLength,end*trackLength);
}

void mvapid::MVAAlg::FitAndSortTrack ( art::Ptr< recob::Track >  track, int &  isStoppingReco, SortedObj &  sortedObj  )

private

Definition at line 374 of file MVAAlg.cxx.

References mvapid::MVAAlg::SortedObj::dir, mvapid::MVAAlg::SortedObj::end, recob::Track::End(), fCheatVertex, fHitsToSpacePoints, fTracksToHits, fVertex4Vect, mvapid::MVAAlg::SortedObj::hitMap, hits(), IsInActiveVol(), mvapid::MVAAlg::SortedObj::length, LinFit(), mvapid::MVAAlg::SortedObj::start, track, recob::Track::Vertex(), and recob::SpacePoint::XYZ().

Referenced by mvapid::MVAAlg::SumDistance2::operator()(), and RunPID().

                                                                        {
  
  sortedTrack.hitMap.clear();
  TVector3 trackPoint,trackDir;
  this->LinFit(track,trackPoint,trackDir);

  TVector3 nearestPointStart, nearestPointEnd;

  //For single-particle events can opt to cheat vertex from start of primary trajectory.
  //Ok since in real events it should be possible to identify the true vertex.
  if(fCheatVertex){
    if((track->End()-fVertex4Vect.Vect()).Mag()>(track->Vertex()-fVertex4Vect.Vect()).Mag()){
      nearestPointStart = trackPoint+trackDir*(trackDir.Dot(track->Vertex()-trackPoint)/trackDir.Mag2());
      nearestPointEnd = trackPoint+trackDir*(trackDir.Dot(track->End()-trackPoint)/trackDir.Mag2());
      isStoppingReco = this->IsInActiveVol(track->End());
    }
    else{
      nearestPointStart = trackPoint+trackDir*(trackDir.Dot(track->End()-trackPoint)/trackDir.Mag2());
      nearestPointEnd = trackPoint+trackDir*(trackDir.Dot(track->Vertex()-trackPoint)/trackDir.Mag2());
      isStoppingReco = this->IsInActiveVol(track->Vertex());
      trackDir*=-1.;
    }
  }
  else{
    if(track->End().Z() >= track->Vertex().Z()){ //Otherwise assume particle is forward-going for now...
      nearestPointStart = trackPoint+trackDir*(trackDir.Dot(track->Vertex()-trackPoint)/trackDir.Mag2());
      nearestPointEnd = trackPoint+trackDir*(trackDir.Dot(track->End()-trackPoint)/trackDir.Mag2());
      isStoppingReco = this->IsInActiveVol(track->End());
    }
    else{
      nearestPointStart = trackPoint+trackDir*(trackDir.Dot(track->End()-trackPoint)/trackDir.Mag2());
      nearestPointEnd = trackPoint+trackDir*(trackDir.Dot(track->Vertex()-trackPoint)/trackDir.Mag2());
      isStoppingReco = this->IsInActiveVol(track->Vertex());
    }

    if(trackDir.Z() <= 0){
      trackDir.SetX(-trackDir.X());
      trackDir.SetY(-trackDir.Y());
      trackDir.SetZ(-trackDir.Z());
    }
  }

  sortedTrack.start=nearestPointStart;
  sortedTrack.end=nearestPointEnd;
  sortedTrack.dir=trackDir;
  sortedTrack.length=(nearestPointEnd-nearestPointStart).Mag();

  std::vector<art::Ptr<recob::Hit>> hits=fTracksToHits[track];

  for(auto hitIter=hits.begin();hitIter!=hits.end();++hitIter){
    
    if(!fHitsToSpacePoints.count(*hitIter)) continue;
    art::Ptr<recob::SpacePoint> sp=fHitsToSpacePoints.at(*hitIter);

    TVector3 nearestPoint = trackPoint+trackDir*(trackDir.Dot(TVector3(sp->XYZ())-trackPoint)/trackDir.Mag2());
    double lengthAlongTrack=(nearestPointStart-nearestPoint).Mag();
    sortedTrack.hitMap.insert(std::pair<double,art::Ptr<recob::Hit> >(lengthAlongTrack,*hitIter));
  }
}

void mvapid::MVAAlg::GetDetectorEdges

(

)

Definition at line 89 of file MVAAlg.cxx.

References fDetMaxX, fDetMaxY, fDetMaxZ, fDetMinX, fDetMinY, fDetMinZ, geo::BoxBoundedGeo::MaxX(), geo::BoxBoundedGeo::MaxY(), geo::BoxBoundedGeo::MaxZ(), geo::BoxBoundedGeo::MinX(), geo::BoxBoundedGeo::MinY(), geo::BoxBoundedGeo::MinZ(), geo::GeometryCore::TotalNTPC(), and geo::GeometryCore::TPC().

Referenced by mvapid::MVAPID::beginJob(), and mvapid::MVAAlg::SumDistance2::operator()().

{

  art::ServiceHandle<geo::Geometry> geom;

  fDetMinX = 999999.9;
  fDetMaxX = -999999.9;
  fDetMinY = 999999.9;
  fDetMaxY = -999999.9;
  fDetMinZ = 999999.9;
  fDetMaxZ = -999999.9;

  for(unsigned int t=0; t<geom->TotalNTPC(); t++)
    {
      if(geom->TPC(t).MinX() < fDetMinX)
        fDetMinX = geom->TPC(t).MinX();
      if(geom->TPC(t).MaxX() > fDetMaxX)
        fDetMaxX = geom->TPC(t).MaxX();
      if(geom->TPC(t).MinY() < fDetMinY)
        fDetMinY = geom->TPC(t).MinY();
      if(geom->TPC(t).MaxY() > fDetMaxY)
        fDetMaxY = geom->TPC(t).MaxY();
      if(geom->TPC(t).MinZ() < fDetMinZ)
        fDetMinZ = geom->TPC(t).MinZ();
      if(geom->TPC(t).MaxZ() > fDetMaxZ)
        fDetMaxZ = geom->TPC(t).MaxZ();
    }
}

void mvapid::MVAAlg::GetWireNormals

(

)

Definition at line 118 of file MVAAlg.cxx.

References fNormToWiresY, fNormToWiresZ, geo::GeometryCore::IteratePlanes(), geo::GeometryCore::Nplanes(), and geo::GeometryCore::NTPC().

Referenced by mvapid::MVAPID::beginJob(), and mvapid::MVAAlg::SumDistance2::operator()().

{

  art::ServiceHandle<geo::Geometry> geom;

  fNormToWiresY.clear();
  fNormToWiresZ.clear();

  int planeKey;

  //Get normals to wires for each plane in the detector
  //This assumes equal numbers of TPCs in each cryostat and equal numbers of planes in each TPC
  for (geo::PlaneGeo const& plane: geom->IteratePlanes()) {
    std::string id = std::string(plane.ID());
    int pcryo = id.find("C");
    int ptpc = id.find("T");
    int pplane = id.find("P");
    std::string scryo = id.substr(pcryo+2,2);
    std::string stpc = id.substr(ptpc+2,2);
    std::string splane = id.substr(pplane+2,2);
    int icryo = std::stoi(scryo);
    int itpc = std::stoi(stpc);
    int iplane = std::stoi(splane);
    planeKey = icryo * geom->NTPC(0) * geom->Nplanes(0, 0) + itpc * geom->Nplanes(0, 0) + iplane;  //single index for all planes in detector
    fNormToWiresY.insert(std::make_pair(planeKey, -plane.Wire(0).Direction().Z()));  //y component of normal 
    fNormToWiresZ.insert(std::make_pair(planeKey, plane.Wire(0).Direction().Y()));   //z component of normal

  }
}

int mvapid::MVAAlg::IsInActiveVol ( const TVector3 &  pos )

private

Definition at line 76 of file MVAAlg.cxx.

References fDetMaxX, fDetMaxY, fDetMaxZ, fDetMinX, fDetMinY, and fDetMinZ.

Referenced by FitAndSortTrack(), mvapid::MVAAlg::SumDistance2::operator()(), and SortShower().

{
  const double fiducialDist = 5.0;

  if(pos.X() > (fDetMinX + fiducialDist) && pos.X() < (fDetMaxX - fiducialDist)
     && pos.Y() > (fDetMinY + fiducialDist) && pos.Y() < (fDetMaxY - fiducialDist)
     && pos.Z() > (fDetMinZ + fiducialDist) && pos.Z() < (fDetMaxZ - fiducialDist))
    return 1;
  else
    return 0;

}

int mvapid::MVAAlg::LinFit ( const art::Ptr< recob::Track >  track, TVector3 &  trackPoint, TVector3 &  trackDir  )

private

Definition at line 661 of file MVAAlg.cxx.

References recob::Track::End(), fTracksToSpacePoints, and recob::Track::Vertex().

Referenced by FitAndSortTrack(), and mvapid::MVAAlg::SumDistance2::operator()().

                                                                                              {

  const std::vector<art::Ptr<recob::SpacePoint> >& sp = fTracksToSpacePoints.at(track);

    TGraph2D grFit(1);
    unsigned int iPt=0;
    for(auto spIter=sp.begin();spIter!=sp.end();++spIter){
      TVector3 point=(*spIter)->XYZ();
      grFit.SetPoint(iPt++,point.X(),point.Y(),point.Z());
    }

    //Lift from the ROOT line3Dfit.C tutorial
    ROOT::Fit::Fitter fitter;
    // make the functor object
    mvapid::MVAAlg::SumDistance2 sdist(&grFit);

    ROOT::Math::Functor fcn(sdist,6);

    //Initial fit parameters from track start and end...
    TVector3 trackStart=track->Vertex();
    TVector3 trackEnd=track->End();
    trackDir=(trackEnd-trackStart).Unit();

    TVector3 x0=trackStart-trackDir;
    TVector3 u=trackDir;

    double pStart[6] = {x0.X(),u.X(),x0.Y(),u.Y(),x0.Z(),u.Z()};

    fitter.SetFCN(fcn,pStart);

    bool ok = fitter.FitFCN();
    if (!ok) {
      trackPoint.SetXYZ(x0.X(),x0.Y(),x0.Z());
      trackDir.SetXYZ(u.X(),u.Y(),u.Z());
      trackDir=trackDir.Unit();
      return 1;
    }
    else{
      const ROOT::Fit::FitResult & result = fitter.Result();
      const double * parFit = result.GetParams();
      trackPoint.SetXYZ(parFit[0],parFit[2],parFit[4]);
      trackDir.SetXYZ(parFit[1],parFit[3],parFit[5]);
      trackDir=trackDir.Unit();
      return 0;
    }
}

int mvapid::MVAAlg::LinFitShower ( const art::Ptr< recob::Shower >  shower, TVector3 &  showerPoint, TVector3 &  showerDir  )

private

Definition at line 708 of file MVAAlg.cxx.

References recob::Shower::Direction(), fShowersToSpacePoints, and recob::Shower::ShowerStart().

Referenced by mvapid::MVAAlg::SumDistance2::operator()(), and SortShower().

                                                                                                        {

  const std::vector<art::Ptr<recob::SpacePoint> >& sp = fShowersToSpacePoints.at(shower);

  TGraph2D grFit(1);
  unsigned int iPt=0;
  for(auto spIter=sp.begin();spIter!=sp.end();++spIter){
    TVector3 point=(*spIter)->XYZ();
    grFit.SetPoint(iPt++,point.X(),point.Y(),point.Z());
  }

  //Lift from the ROOT line3Dfit.C tutorial
  ROOT::Fit::Fitter fitter;
  // make the functor object
  mvapid::MVAAlg::SumDistance2 sdist(&grFit);

  ROOT::Math::Functor fcn(sdist,6);

  //Initial fit parameters from shower start and end...
  TVector3 showerStart=shower->ShowerStart();
  showerDir = shower->Direction().Unit();

  TVector3 x0=showerStart-showerDir;
  TVector3 u=showerDir;

  double pStart[6] = {x0.X(),u.X(),x0.Y(),u.Y(),x0.Z(),u.Z()};

  fitter.SetFCN(fcn,pStart);

  bool ok = fitter.FitFCN();
  if (!ok) {
    showerPoint.SetXYZ(x0.X(),x0.Y(),x0.Z());
    showerDir.SetXYZ(u.X(),u.Y(),u.Z());
    showerDir=showerDir.Unit();
    return 1;
  }
  else{
    const ROOT::Fit::FitResult & result = fitter.Result();
    const double * parFit = result.GetParams();
    showerPoint.SetXYZ(parFit[0],parFit[2],parFit[4]);
    showerDir.SetXYZ(parFit[1],parFit[3],parFit[5]);
    showerDir=showerDir.Unit();
    return 0;
  }
}

void mvapid::MVAAlg::PrepareEvent ( const art::Event &  event )

private

Definition at line 272 of file MVAAlg.cxx.

References fCheatVertex, fEventT0, fHitLabel, fHits, fHitsToSpacePoints, fShowerLabel, fShowers, fShowersToHits, fShowersToSpacePoints, fSpacePointLabel, fSpacePoints, fSpacePointsToHits, fTrackingLabel, fTrackLabel, fTracks, fTracksToHits, fTracksToSpacePoints, fVertex4Vect, art::DataViewImpl::getByLabel(), and track.

Referenced by mvapid::MVAAlg::SumDistance2::operator()(), and RunPID().

                                                  {

  fHits.clear();
  fSpacePoints.clear();
  fTracks.clear();
  fShowers.clear();
  fSpacePointsToHits.clear();
  fHitsToSpacePoints.clear();
  fTracksToHits.clear();
  fTracksToSpacePoints.clear();
  fShowersToHits.clear();
  fShowersToSpacePoints.clear();

  auto const* detProp = lar::providerFrom<detinfo::DetectorPropertiesService>();
  fEventT0=detProp->TriggerOffset();

  art::Handle< std::vector<recob::Hit> > hitsHandle;
  evt.getByLabel(fHitLabel, hitsHandle);

  for (unsigned int iHit = 0; iHit < hitsHandle->size(); ++iHit){
    const art::Ptr<recob::Hit> hit(hitsHandle, iHit);
    fHits.push_back(hit);
  }

  art::Handle< std::vector<recob::Track> > tracksHandle;
  evt.getByLabel(fTrackLabel, tracksHandle);

  for (unsigned int iTrack = 0; iTrack < tracksHandle->size(); ++iTrack){
    const art::Ptr<recob::Track> track(tracksHandle, iTrack);
    fTracks.push_back(track);
  }

  art::Handle< std::vector<recob::Shower> > showersHandle;
  evt.getByLabel(fShowerLabel, showersHandle);

  for (unsigned int iShower = 0; iShower < showersHandle->size(); ++iShower){
    const art::Ptr<recob::Shower> shower(showersHandle, iShower);
    fShowers.push_back(shower);
  }

  art::Handle< std::vector<recob::SpacePoint> > spHandle;
  evt.getByLabel(fSpacePointLabel, spHandle);

  for (unsigned int iSP = 0; iSP < spHandle->size(); ++iSP){
    const art::Ptr<recob::SpacePoint> spacePoint(spHandle, iSP);
    fSpacePoints.push_back(spacePoint);
  }

  art::FindManyP<recob::Hit> findTracksToHits(fTracks,evt,fTrackLabel);
  art::FindManyP<recob::Hit> findShowersToHits(fShowers,evt,fShowerLabel);
  art::FindOneP<recob::Hit> findSPToHits(fSpacePoints,evt,fSpacePointLabel);
  
  for(unsigned int iSP = 0; iSP < fSpacePoints.size(); ++iSP)
    {
      const art::Ptr<recob::SpacePoint> spacePoint=fSpacePoints.at(iSP);

      const art::Ptr<recob::Hit> hit = findSPToHits.at(iSP);
      fSpacePointsToHits[spacePoint]=hit;
      fHitsToSpacePoints[hit]=spacePoint;
  }     
  
  for(unsigned int iTrack = 0; iTrack < fTracks.size(); ++iTrack){
    const art::Ptr<recob::Track> track=fTracks.at(iTrack);
    
    const std::vector< art::Ptr<recob::Hit> > trackHits = findTracksToHits.at(iTrack);

    for (unsigned int iHit=0; iHit<trackHits.size(); ++iHit)
    {
      const art::Ptr<recob::Hit> hit = trackHits.at(iHit);          
      fTracksToHits[track].push_back(hit);
      if(fHitsToSpacePoints.count(hit)){
        fTracksToSpacePoints[track].push_back(fHitsToSpacePoints.at(hit));
      }
    }
  }

  for(unsigned int iShower = 0; iShower < fShowers.size(); ++iShower){
    const art::Ptr<recob::Shower> shower=fShowers.at(iShower);
    const std::vector< art::Ptr<recob::Hit> > showerHits = findShowersToHits.at(iShower);

    for (unsigned int iHit=0; iHit<showerHits.size(); ++iHit)
    {
      const art::Ptr<recob::Hit> hit = showerHits.at(iHit);          
      fShowersToHits[shower].push_back(hit);
      if(fHitsToSpacePoints.count(hit)){
        fShowersToSpacePoints[shower].push_back(fHitsToSpacePoints.at(hit));
      }
    }
  }

  if(fCheatVertex){
    art::Handle< std::vector<simb::MCParticle> > partHandle;
    evt.getByLabel(fTrackingLabel, partHandle);

    if(partHandle->size()==0||partHandle->at(0).TrackId()!=1){
      std::cout<<"Error, ID of first track in largeant list is not 0"<<std::endl;
      exit(1);
    } 
    fVertex4Vect=partHandle->at(0).Position();
  }
}

void mvapid::MVAAlg::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 74 of file MVAAlg.cxx.

Referenced by mvapid::MVAAlg::SumDistance2::operator()().

{}

void mvapid::MVAAlg::RunPCA ( std::vector< art::Ptr< recob::Hit > > &  hits, std::vector< double > &  eVals, std::vector< double > &  eVecs  )

private

Definition at line 513 of file MVAAlg.cxx.

References fHitsToSpacePoints, and hits().

Referenced by mvapid::MVAAlg::SumDistance2::operator()(), and RunPID().

                                                                                                              {

  // Define the TPrincipal
  TPrincipal* principal = new TPrincipal(3,"D");
  // Define variables to hold the eigenvalues and eigenvectors
  //const TMatrixD* covar = new TMatrixD();
  //const TVectorD* meanval = new TVectorD();

  for(auto hitIter=hits.begin();hitIter!=hits.end();++hitIter){

    if(fHitsToSpacePoints.count(*hitIter)){
      principal->AddRow( fHitsToSpacePoints.at(*hitIter)->XYZ() );
    }
  }
  
  // PERFORM PCA
  principal->MakePrincipals();
  // GET EIGENVALUES AND EIGENVECTORS
  for(unsigned int i=0;i<3;++i){
    eVals.push_back(principal->GetEigenValues()->GetMatrixArray()[i]);
  }

  for(unsigned int i=0;i<9;++i){
    eVecs.push_back(principal->GetEigenVectors()->GetMatrixArray()[i]);
  }
}

void mvapid::MVAAlg::RunPID	(	art::Event &	evt,
		std::vector< anab::MVAPIDResult > &	result,
		art::Assns< recob::Track, anab::MVAPIDResult, void > &	trackAssns,
		art::Assns< recob::Shower, anab::MVAPIDResult, void > &	showerAssns
	)

Definition at line 148 of file MVAAlg.cxx.

References _Var_Shape(), CalcSegmentdEdxFrac(), anab::MVAPIDResult::concentration, anab::MVAPIDResult::conicalness, anab::MVAPIDResult::coreHaloRatio, util::CreateAssn(), anab::MVAPIDResult::dEdxEnd, anab::MVAPIDResult::dEdxEndRatio, anab::MVAPIDResult::dEdxStart, anab::MVAPIDResult::evalRatio, FitAndSortTrack(), fMVAMethods, fParentModule, fReader, fResHolder, fShowers, fShowersToHits, fTracks, fTracksToHits, mvapid::MVAAlg::SortedObj::hitMap, anab::MVAPIDResult::isStoppingReco, anab::MVAPIDResult::isTrack, anab::MVAPIDResult::length, mvapid::MVAAlg::SortedObj::length, anab::MVAPIDResult::mvaOutput, anab::MVAPIDResult::nSpacePoints, PrepareEvent(), RunPCA(), SortShower(), and anab::MVAPIDResult::trackID.

Referenced by mvapid::MVAAlg::SumDistance2::operator()(), and mvapid::MVAPID::produce().

                                                                                       {

  this->PrepareEvent(evt);

  for(auto trackIter=fTracks.begin();trackIter!=fTracks.end();++trackIter){
    mvapid::MVAAlg::SortedObj sortedObj;

    std::vector<double> eVals,eVecs;
    int isStoppingReco; 
    this->RunPCA(fTracksToHits[*trackIter],eVals,eVecs);
    double evalRatio;
    if(eVals[0] < 0.0001)
      evalRatio=0.0;
    else
      evalRatio=std::sqrt(eVals[1]*eVals[1]+eVals[2]*eVals[2])/eVals[0];
    this->FitAndSortTrack(*trackIter,isStoppingReco,sortedObj);
    double coreHaloRatio,concentration,conicalness;
    this->_Var_Shape(sortedObj,coreHaloRatio,concentration,conicalness);
    double dEdxStart = CalcSegmentdEdxFrac(sortedObj,0.,0.05);
    double dEdxEnd = CalcSegmentdEdxFrac(sortedObj,0.9,1.0);
    double dEdxPenultimate = CalcSegmentdEdxFrac(sortedObj,0.8,0.9);

    /*    
    std::cout<<"coreHaloRatio:   "<<coreHaloRatio<<std::endl;
    std::cout<<"concentration:   "<<concentration<<std::endl;
    std::cout<<"conicalness:     "<<conicalness<<std::endl;
    std::cout<<"dEdxStart: "<<dEdxStart<<std::endl;
    std::cout<<"dEdxEnd: "<<dEdxEnd<<std::endl;
    std::cout<<"dEdxEndRatio: ";
    if(dEdxPenultimate < 0.1)
    std::cout<<"1.0";
    else
    std::cout<<dEdxEnd/dEdxPenultimate;
    std::cout<<std::endl;
    */

    fResHolder.isTrack=1;
    fResHolder.isStoppingReco=isStoppingReco;
    fResHolder.nSpacePoints=sortedObj.hitMap.size();
    fResHolder.trackID=(*trackIter)->ID();
    fResHolder.evalRatio=evalRatio;
    fResHolder.concentration=concentration;
    fResHolder.coreHaloRatio=coreHaloRatio;
    fResHolder.conicalness=conicalness;
    fResHolder.dEdxStart=dEdxStart;
    fResHolder.dEdxEnd=dEdxEnd;
    if(dEdxPenultimate < 0.1)
      fResHolder.dEdxEndRatio=1.0;
    else
      fResHolder.dEdxEndRatio=dEdxEnd/dEdxPenultimate;
    fResHolder.length=sortedObj.length;

    for(auto methodIter=fMVAMethods.begin();methodIter!=fMVAMethods.end();++methodIter){
      fResHolder.mvaOutput[*methodIter]=fReader.EvaluateMVA(*methodIter);
    }
    result.push_back(fResHolder);
    util::CreateAssn(*fParentModule, evt, result, *trackIter, trackAssns);
  }

  for(auto showerIter=fShowers.begin();showerIter!=fShowers.end();++showerIter){
    mvapid::MVAAlg::SortedObj sortedObj;

    std::vector<double> eVals,eVecs;
    int isStoppingReco;

    this->RunPCA(fShowersToHits[*showerIter],eVals,eVecs);

    double evalRatio;
    if(eVals[0]< 0.0001)
      evalRatio=0.0;
    else
      evalRatio=std::sqrt(eVals[1]*eVals[1]+eVals[2]*eVals[2])/eVals[0];

    //this->SortShower(*showerIter,TVector3(eVecs[0],eVecs[3],eVecs[6]),isStoppingReco,sortedObj);
    this->SortShower(*showerIter,isStoppingReco,sortedObj);

    double coreHaloRatio,concentration,conicalness;
    this->_Var_Shape(sortedObj,coreHaloRatio,concentration,conicalness);
    double dEdxStart = CalcSegmentdEdxFrac(sortedObj,0.,0.05);
    double dEdxEnd = CalcSegmentdEdxFrac(sortedObj,0.9,1.0);
    double dEdxPenultimate = CalcSegmentdEdxFrac(sortedObj,0.8,0.9);

    /*    
    std::cout<<"coreHaloRatio:   "<<coreHaloRatio<<std::endl;
    std::cout<<"concentration:   "<<concentration<<std::endl;
    std::cout<<"conicalness:     "<<conicalness<<std::endl;
    std::cout<<"dEdxStart: "<<dEdxStart<<std::endl;
    std::cout<<"dEdxEnd: "<<dEdxEnd<<std::endl;
    std::cout<<"dEdxEndRatio: ";
    if(dEdxPenultimate < 0.1)
    std::cout<<"1.0";
    else
    std::cout<<dEdxEnd/dEdxPenultimate;
    std::cout<<std::endl;
    */

    fResHolder.isTrack=0;
    fResHolder.isStoppingReco=isStoppingReco;
    fResHolder.nSpacePoints=sortedObj.hitMap.size();
    fResHolder.trackID=(*showerIter)->ID()+1000; //For the moment label showers by adding 1000 to ID

    fResHolder.evalRatio=evalRatio;
    fResHolder.concentration=concentration;
    fResHolder.coreHaloRatio=coreHaloRatio;
    fResHolder.conicalness=conicalness;
    fResHolder.dEdxStart=dEdxStart;
    fResHolder.dEdxEnd=dEdxEnd;
    if(dEdxPenultimate < 0.1)
      fResHolder.dEdxEndRatio=1.0;
    else
      fResHolder.dEdxEndRatio=dEdxEnd/dEdxPenultimate;
    fResHolder.length=sortedObj.length;

    for(auto methodIter=fMVAMethods.begin();methodIter!=fMVAMethods.end();++methodIter){
      fResHolder.mvaOutput[*methodIter]=fReader.EvaluateMVA(*methodIter);
    }
    result.push_back(fResHolder);
    util::CreateAssn(*fParentModule, evt, result, *showerIter, showerAssns);
  }

}

void mvapid::MVAAlg::SortShower ( art::Ptr< recob::Shower >  shower, int &  isStoppingReco, mvapid::MVAAlg::SortedObj &  sortedShower  )

private

Definition at line 437 of file MVAAlg.cxx.

References mvapid::MVAAlg::SortedObj::dir, mvapid::MVAAlg::SortedObj::end, fCheatVertex, fHitsToSpacePoints, fShowersToHits, fVertex4Vect, mvapid::MVAAlg::SortedObj::hitMap, hits(), IsInActiveVol(), mvapid::MVAAlg::SortedObj::length, LinFitShower(), recob::Shower::ShowerStart(), mvapid::MVAAlg::SortedObj::start, and recob::SpacePoint::XYZ().

Referenced by mvapid::MVAAlg::SumDistance2::operator()(), and RunPID().

                                                                      {
  sortedShower.hitMap.clear();
  
  std::vector<art::Ptr<recob::Hit>> hits=fShowersToHits[shower];

  TVector3 showerEnd(0, 0, 0); 
  double furthestHitFromStart = -999.9;
  for(auto hitIter=hits.begin();hitIter!=hits.end();++hitIter){

    if(!fHitsToSpacePoints.count(*hitIter)) continue;
    art::Ptr<recob::SpacePoint> sp=fHitsToSpacePoints.at(*hitIter);
    if((TVector3(sp->XYZ()) - shower->ShowerStart()).Mag() > furthestHitFromStart)
      {
        showerEnd = TVector3(sp->XYZ());
        furthestHitFromStart = (TVector3(sp->XYZ()) - shower->ShowerStart()).Mag();
      }
  }

  TVector3 showerPoint,showerDir;
  this->LinFitShower(shower,showerPoint,showerDir);

  TVector3 nearestPointStart, nearestPointEnd;

  //Ensure that shower is fitted in correct direction (assuming for now that particle moves in +z direction)

  if(fCheatVertex){
    if((showerEnd-fVertex4Vect.Vect()).Mag()>(shower->ShowerStart()-fVertex4Vect.Vect()).Mag()){
      nearestPointStart = showerPoint+showerDir*(showerDir.Dot(shower->ShowerStart()-showerPoint)/showerDir.Mag2());
      nearestPointEnd = showerPoint+showerDir*(showerDir.Dot(showerEnd-showerPoint)/showerDir.Mag2());
      isStoppingReco = this->IsInActiveVol(showerEnd);
    }
    else
      {
        nearestPointStart = showerPoint+showerDir*(showerDir.Dot(showerEnd-showerPoint)/showerDir.Mag2());
        nearestPointEnd = showerPoint+showerDir*(showerDir.Dot(shower->ShowerStart()-showerPoint)/showerDir.Mag2());
        isStoppingReco = this->IsInActiveVol(shower->ShowerStart());
        showerDir*=-1.;
      }
  }
  else{
    if(showerEnd.Z() >= shower->ShowerStart().Z()){
      nearestPointStart = showerPoint+showerDir*(showerDir.Dot(shower->ShowerStart()-showerPoint)/showerDir.Mag2());
      nearestPointEnd = showerPoint+showerDir*(showerDir.Dot(showerEnd-showerPoint)/showerDir.Mag2());
      isStoppingReco = this->IsInActiveVol(showerEnd);
    }
    else{
      nearestPointStart = showerPoint+showerDir*(showerDir.Dot(showerEnd-showerPoint)/showerDir.Mag2());
      nearestPointEnd = showerPoint+showerDir*(showerDir.Dot(shower->ShowerStart()-showerPoint)/showerDir.Mag2());
      isStoppingReco = this->IsInActiveVol(shower->ShowerStart());
    }
    
    if(showerDir.Z() <= 0){
      showerDir.SetX(-showerDir.X());
      showerDir.SetY(-showerDir.Y());
      showerDir.SetZ(-showerDir.Z());
    }  
  }

  sortedShower.start=nearestPointStart;
  sortedShower.end=nearestPointEnd;
  //sortedShower.dir=dir;
  sortedShower.dir=showerDir;
  sortedShower.length=(nearestPointEnd-nearestPointStart).Mag();

  for(auto hitIter=hits.begin();hitIter!=hits.end();++hitIter){
    
    if(!fHitsToSpacePoints.count(*hitIter)) continue;
    art::Ptr<recob::SpacePoint> sp=fHitsToSpacePoints.at(*hitIter);

    TVector3 nearestPoint = showerPoint+showerDir*(showerDir.Dot(TVector3(sp->XYZ())-showerPoint)/showerDir.Mag2());
    double lengthAlongShower=(nearestPointStart-nearestPoint).Mag();
    sortedShower.hitMap.insert(std::pair<double,art::Ptr<recob::Hit> >(lengthAlongShower,*hitIter));
  }
  
}

Member Data Documentation

const calo::CalorimetryAlg mvapid::MVAAlg::fCaloAlg

private

Definition at line 125 of file MVAAlg.h.

Referenced by CalcSegmentdEdxDist().

bool mvapid::MVAAlg::fCheatVertex

private

Definition at line 161 of file MVAAlg.h.

Referenced by FitAndSortTrack(), MVAAlg(), PrepareEvent(), and SortShower().

double mvapid::MVAAlg::fDetMaxX

private

Definition at line 129 of file MVAAlg.h.

Referenced by GetDetectorEdges(), and IsInActiveVol().

double mvapid::MVAAlg::fDetMaxY

private

Definition at line 129 of file MVAAlg.h.

Referenced by GetDetectorEdges(), and IsInActiveVol().

double mvapid::MVAAlg::fDetMaxZ

private

Definition at line 129 of file MVAAlg.h.

Referenced by GetDetectorEdges(), and IsInActiveVol().

double mvapid::MVAAlg::fDetMinX

private

Definition at line 129 of file MVAAlg.h.

Referenced by GetDetectorEdges(), and IsInActiveVol().

double mvapid::MVAAlg::fDetMinY

private

Definition at line 129 of file MVAAlg.h.

Referenced by GetDetectorEdges(), and IsInActiveVol().

double mvapid::MVAAlg::fDetMinZ

private

Definition at line 129 of file MVAAlg.h.

Referenced by GetDetectorEdges(), and IsInActiveVol().

double mvapid::MVAAlg::fEventT0

private

Definition at line 127 of file MVAAlg.h.

Referenced by CalcSegmentdEdxDist(), and PrepareEvent().

std::string mvapid::MVAAlg::fHitLabel

private

Definition at line 138 of file MVAAlg.h.

Referenced by MVAAlg(), and PrepareEvent().

std::vector<art::Ptr<recob::Hit> > mvapid::MVAAlg::fHits

private

Definition at line 145 of file MVAAlg.h.

Referenced by PrepareEvent().

std::map<art::Ptr<recob::Hit>,art::Ptr<recob::SpacePoint> > mvapid::MVAAlg::fHitsToSpacePoints

private

Definition at line 151 of file MVAAlg.h.

Referenced by _Var_Shape(), FitAndSortTrack(), PrepareEvent(), RunPCA(), and SortShower().

std::vector<std::string> mvapid::MVAAlg::fMVAMethods

private

Definition at line 158 of file MVAAlg.h.

Referenced by MVAAlg(), and RunPID().

std::map<int,double> mvapid::MVAAlg::fNormToWiresY

private

Definition at line 131 of file MVAAlg.h.

Referenced by CalcSegmentdEdxDist(), and GetWireNormals().

std::map<int,double> mvapid::MVAAlg::fNormToWiresZ

private

Definition at line 132 of file MVAAlg.h.

Referenced by CalcSegmentdEdxDist(), and GetWireNormals().

const art::EDProducer* mvapid::MVAAlg::fParentModule

private

Definition at line 134 of file MVAAlg.h.

Referenced by RunPID().

TMVA::Reader mvapid::MVAAlg::fReader

private

Definition at line 156 of file MVAAlg.h.

Referenced by MVAAlg(), and RunPID().

anab::MVAPIDResult mvapid::MVAAlg::fResHolder

private

Definition at line 154 of file MVAAlg.h.

Referenced by MVAAlg(), and RunPID().

std::string mvapid::MVAAlg::fShowerLabel

private

Definition at line 137 of file MVAAlg.h.

Referenced by MVAAlg(), and PrepareEvent().

std::vector<art::Ptr<recob::Shower> > mvapid::MVAAlg::fShowers

private

Definition at line 143 of file MVAAlg.h.

Referenced by PrepareEvent(), and RunPID().

std::map<art::Ptr<recob::Shower>,std::vector<art::Ptr<recob::Hit> > > mvapid::MVAAlg::fShowersToHits

private

Definition at line 149 of file MVAAlg.h.

Referenced by PrepareEvent(), RunPID(), and SortShower().

std::map<art::Ptr<recob::Shower>,std::vector<art::Ptr<recob::SpacePoint> > > mvapid::MVAAlg::fShowersToSpacePoints

private

Definition at line 150 of file MVAAlg.h.

Referenced by LinFitShower(), and PrepareEvent().

std::string mvapid::MVAAlg::fSpacePointLabel

private

Definition at line 139 of file MVAAlg.h.

Referenced by MVAAlg(), and PrepareEvent().

std::vector<art::Ptr<recob::SpacePoint> > mvapid::MVAAlg::fSpacePoints

private

Definition at line 144 of file MVAAlg.h.

Referenced by PrepareEvent().

std::map<art::Ptr<recob::SpacePoint>,art::Ptr<recob::Hit > > mvapid::MVAAlg::fSpacePointsToHits

private

Definition at line 152 of file MVAAlg.h.

Referenced by PrepareEvent().

std::string mvapid::MVAAlg::fTrackingLabel

private

Definition at line 140 of file MVAAlg.h.

Referenced by MVAAlg(), and PrepareEvent().

std::string mvapid::MVAAlg::fTrackLabel

private

Definition at line 136 of file MVAAlg.h.

Referenced by MVAAlg(), and PrepareEvent().

std::vector<art::Ptr<recob::Track> > mvapid::MVAAlg::fTracks

private

Definition at line 142 of file MVAAlg.h.

Referenced by PrepareEvent(), and RunPID().

std::map<art::Ptr<recob::Track>,std::vector<art::Ptr<recob::Hit> > > mvapid::MVAAlg::fTracksToHits

private

Definition at line 147 of file MVAAlg.h.

Referenced by FitAndSortTrack(), PrepareEvent(), and RunPID().

std::map<art::Ptr<recob::Track>,std::vector<art::Ptr<recob::SpacePoint> > > mvapid::MVAAlg::fTracksToSpacePoints

private

Definition at line 148 of file MVAAlg.h.

Referenced by LinFit(), and PrepareEvent().

TLorentzVector mvapid::MVAAlg::fVertex4Vect

private

Definition at line 163 of file MVAAlg.h.

Referenced by FitAndSortTrack(), PrepareEvent(), and SortShower().

std::vector<std::string> mvapid::MVAAlg::fWeightFiles

private

Definition at line 159 of file MVAAlg.h.

Referenced by MVAAlg().

---

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

• larana/v06_21_01/source/larana/ParticleIdentification/MVAAlg.h
• larana/v06_21_01/source/larana/ParticleIdentification/MVAAlg.cxx