tf::TimingFitAlg Class Reference

Algorithm to determine the direction of tracks using timing. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-23/TimingFit/TimingFitAlg.h"

Public Member Functions
	TimingFitAlg ()
TimingFitResult	HoughFit (const rb::Track *track) const
TimingFitResult	HoughFitPts (const std::vector< Pt > &pts) const
	Helper for HoughFit. May be useful for more advanced users. More...

Detailed Description

Algorithm to determine the direction of tracks using timing.

Definition at line 26 of file TimingFitAlg.h.

Constructor & Destructor Documentation

tf::TimingFitAlg::TimingFitAlg

(

)

Definition at line 31 of file TimingFitAlg.cxx.

  {
  }

Member Function Documentation

TimingFitResult tf::TimingFitAlg::HoughFit

(

const rb::Track *

track

)

const

Definition at line 153 of file TimingFitAlg.cxx.

References rb::Cluster::Cell(), cellMax, dir, rb::Prong::Dir(), Dot(), calib::Calibrator::GetTimeRes(), HoughFitPts(), rb::Cluster::MeanTNS(), rb::Cluster::MeanXYZ(), rb::Cluster::NCell(), elec2geo::pos, gen_hdf5record::pt, rb::Cluster::RecoHit(), rb::RecoHit::T(), confusionMatrixTree::t, getGoodRuns4SAM::t0, rb::RecoHit::X(), rb::RecoHit::Y(), and rb::RecoHit::Z().

Referenced by tf::TimingFit::produce(), and cosrej::CosRej::produce().

  {
    std::vector<Pt> pts;

    art::ServiceHandle<calib::Calibrator> cal;

    const unsigned int cellMax = track->NCell();
    const TVector3 dir = track->Dir();
    // Center Pts around (0,0). Makes natural choice of MINUIT seeding
    const TVector3 origin = track->MeanXYZ();
    const double t0 = track->MeanTNS();
    for(unsigned int cellIdx = 0; cellIdx < cellMax; ++cellIdx){
      const art::Ptr<rb::CellHit> chit = track->Cell(cellIdx);

      const rb::RecoHit rhit = track->RecoHit(chit);
      const double t = rhit.T();

      const TVector3 pos(rhit.X(), rhit.Y(), rhit.Z());

      Pt pt((pos-origin).Dot(dir), t-t0, cal->GetTimeRes(*chit));
      pts.push_back(pt);
    } // end for cellIdx

    return HoughFitPts(pts);
  }

TimingFitResult tf::TimingFitAlg::HoughFitPts

(

const std::vector< Pt > &

pts

)

const

Helper for HoughFit. May be useful for more advanced users.

Definition at line 79 of file TimingFitAlg.cxx.

References tf::TimingFitResult::bestInvSpeed, tf::TimingFitResult::bestScore, tf::TimingFitResult::fwdScore, runNovaSAM::ret, tf::TimingFitResult::revScore, tf::TimingFitErrorHandler(), submit_hadd::u, and tf::TimingFitResult::valid.

Referenced by HoughFit().

  {
    TimingFitResult ret;
    if(pts.size() < 2){
      ret.valid = false;
      return ret;
    }

    gOldHandler = SetErrorHandler(TimingFitErrorHandler);

    HoughScore hs(pts);

    ROOT::Minuit2::MnUserParameters mnPars;
    mnPars.Add("t0", 0, 100);
    mnPars.Add("cinv", 0, 1./30);

    // Best fit
    ROOT::Minuit2::MnMigrad mnMigrad(hs, mnPars);
    ROOT::Minuit2::FunctionMinimum funcMin = mnMigrad();
    ret.bestInvSpeed = mnMigrad.Parameter(1).Value();
    ret.bestScore = -funcMin.Fval()/pts.size();
    ret.valid = funcMin.IsValid();

    // Degugging code to look at fit failures
    /*
    if(!ret.valid){
      TH2F h("", ";t_{0} (ns);1/c (ns/cm)", 100, -300, 300, 100, -.2, +.2);

      for(int x = 0; x < h.GetNbinsX()+2; ++x){
        const double t0 = h.GetXaxis()->GetBinCenter(x);
        for(int y = 0; y < h.GetNbinsY()+2; ++y){
          const double cinv = h.GetYaxis()->GetBinCenter(y);

          h.SetBinContent(x, y, hs.Eval(t0, cinv));
        } // end for y (cinv)
      } // end for x (t0)

      h.Draw("colz");
      gPad->Print("map.eps");

      TGraphErrors* g = new TGraphErrors;
      for(unsigned int i = 0; i < pts.size(); ++i){
        g->SetPoint(g->GetN(), pts[i].L, pts[i].t);
        g->SetPointError(g->GetN()-1, 0, pts[i].err);
      }

      g->SetMarkerStyle(kFullCircle);
      g->Draw("ap");
      gPad->Print("t_vs_z.eps");

      exit(0);
    }
    */

    // Best fit constrained to +c
    mnMigrad.SetValue(0u, 0.);
    mnMigrad.SetValue(1, +1./30);
    mnMigrad.Fix(1);
    funcMin = mnMigrad();
    ret.fwdScore = -funcMin.Fval()/pts.size();

    // Best fit constrained to -c
    mnMigrad.SetValue(0u, 0.);
    mnMigrad.SetValue(1, -1./30);
    mnMigrad.Fix(1);
    funcMin = mnMigrad();
    ret.revScore = -funcMin.Fval()/pts.size();

    SetErrorHandler(gOldHandler);

    return ret;
  }

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-23/TimingFit/TimingFitAlg.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-23/TimingFit/TimingFitAlg.cxx