BCReco Class Reference

#include <BCReco.h>

Inheritance diagram for BCReco:

List of all members.

Public Member Functions
	BCReco (const char *version)
	~BCReco ()
JobCResult	Reco (EDMEventHandle &evt)
JobCResult	Ana (const EDMEventHandle &evt)
void	Update (const CfgConfig &c)
Private Member Functions
void	DoCrossings (float z1, float z2, float z3)
void	HitAssoc ()
void	TrkFind ()
void	ViewMatch ()
void	MakeHistos ()
	......................................................................
void	MatchSU ()
void	MatchTV ()
void	TrkFit ()
void	FitTMinuit ()
void	FitLSQ ()
void	CleanUpTrk ()
void	FillNt ()
void	Clear ()
Private Attributes
std::vector< Hit * >	fBCHitList [3][4]
std::list< AssocHit_t * >	fBCAHitList [3][4]
std::list< CandPair_t * >	fCandTrk [4]
std::list< CandTrajectory_t * >	fCandTraj
std::list< TrackFit_t * >	fTrackFit
int	fNTrack
TTree *	fNtTree
	Ntuple tree.
BCRecoNt *	fNtuple
	Branch for the ntuple tree.
TMinuit *	fMinuit
	MINUIT interface.
TH1F *	fTrkFindHisto [4]
	Predicted wire - hit wire in TrkFind.
TH1F *	fViewMatchHisto [3][4]
	Predicted wire - hit wire in ViewMatch.
TH1F *	fnTrackHisto
	Number of tracks.
TH1F *	fx1Histo
	x1
TH1F *	fy1Histo
	y1
TH1F *	fx3Histo
	x3
TH1F *	fy3Histo
	y3
TH1F *	fdxdzHisto
	dx/dz
TH1F *	fdydzHisto
	dy/dz
TH1F *	fdsdzHisto
	ds/dz
TH1F *	ftargetxHisto
	x at target
TH1F *	ftargetyHisto
	y at target
TH2F *	ftargetxyHisto
	x:y at target
TH2F *	fBC3xyHisto
	x:y
TH2F *	fBC3xdxdzHisto
	x:dx/dz
TH2F *	fBC3ydydzHisto
	y:dy/dz
TH1F *	fChi2Histo
	chi2
TH1F *	fRChi2Histo
	Reduced chi2.
TH1F *	fFinalTrkHisto
	Distance(Fit track - hit wire).
int	fDebug
int	fPrintLevel
int	fN_TrkFind
int	fN_ViewMatch
int	fMaxItr
int	fNWireItr
float	fWidth_TrkFit
float	fWidthTrks
int	fFitFunc
int	fMaxNWirePerChamber
int	fMaxTotalNWire
Classes
struct	AssocHit_t
struct	CandPair_t
struct	CandTrajectory_t
struct	TrackFit_t

---

Detailed Description

Definition at line 68 of file BCReco.h.

---

Constructor & Destructor Documentation

BCReco::BCReco

(

const char *

version

)

Definition at line 171 of file BCReco.cxx.

References fDebug, fMinuit, fPrintLevel, fTrkFindHisto, fViewMatchHisto, gsCotAng, gsDelX, gsFCN(), gsXWire1, kBC1, kNPlane, MakeHistos(), and JobCModule::SetCfgVersion().

                                  : 
  JobCModule("BCReco"),
  fNTrack(0),
  fNtTree(0),
  fNtuple(new BCRecoNt),
  fMinuit(new TMinuit(4)),
  fnTrackHisto(0),
  fx1Histo(0),
  fy1Histo(0),
  fx3Histo(0),
  fy3Histo(0),
  fdxdzHisto(0),
  fdydzHisto(0),
  fdsdzHisto(0),
  ftargetxHisto(0),
  ftargetyHisto(0),
  ftargetxyHisto(0),
  fBC3xyHisto(0),
  fBC3xdxdzHisto(0),
  fBC3ydydzHisto(0),
  fChi2Histo(0),
  fRChi2Histo(0),
  fFinalTrkHisto(0),
  fDebug(1),
  fPrintLevel(-1),
  fN_TrkFind(5),
  fN_ViewMatch(3),
  fMaxItr(30),
  fNWireItr(3),
  fWidth_TrkFit(0.5),
  fWidthTrks(0.5),
  fFitFunc(0),
  fMaxNWirePerChamber(10),
  fMaxTotalNWire(25)
{
  // Do this first to load default configuration
  this->SetCfgVersion(version);

  for (int iplane=0; iplane < GDCGeom::kNPlane; ++iplane) {
    gsCotAng[iplane] = 1./
      tan(M_PI/2.-GDCGeom::Id(GDCGeom::kBC1).GetWrAng(iplane+1));
    gsDelX[iplane] = GDCGeom::Id(GDCGeom::kBC1).GetWrDx(iplane+1);
    gsXWire1[iplane] = gsDelX[iplane]*
      (GDCGeom::Id(GDCGeom::kBC1).GetMaxWire(iplane+1)/2.-0.5);
  }

  for (int iplane=0; iplane < GDCGeom::kNPlane; ++iplane) {
    fTrkFindHisto[iplane] = 0;
    for (int ichamber=0; ichamber<3; ++ichamber)
      fViewMatchHisto[ichamber][iplane] = 0;
  }

  fMinuit->SetFCN(gsFCN);
  fMinuit->SetPrintLevel(fPrintLevel);

  if (fDebug) MakeHistos();
}

BCReco::~BCReco

(

)

Definition at line 314 of file BCReco.cxx.

References fBC3xdxdzHisto, fBC3xyHisto, fBC3ydydzHisto, fChi2Histo, fdsdzHisto, fdxdzHisto, fdydzHisto, fFinalTrkHisto, fnTrackHisto, fNtTree, fNtuple, fRChi2Histo, ftargetxHisto, ftargetxyHisto, ftargetyHisto, fTrkFindHisto, fViewMatchHisto, fx1Histo, fx3Histo, fy1Histo, fy3Histo, and Write().

{
  for (int iplane=0; iplane<4; ++iplane) {
    if (fTrkFindHisto[iplane]) { 
      fTrkFindHisto[iplane]->Write(); 
      delete fTrkFindHisto[iplane]; 
    }
    for (int ichamber=0; ichamber<3; ++ichamber)
      if (fViewMatchHisto[ichamber][iplane]) { 
    fViewMatchHisto[ichamber][iplane]->Write(); 
    delete fViewMatchHisto[ichamber][iplane]; 
      }
  }

  if (fnTrackHisto) {
    fnTrackHisto->Write();
    delete fnTrackHisto;
  }
  if (fx1Histo) {
    fx1Histo->Write();
    delete fx1Histo;
  }
  if (fy1Histo) {
    fy1Histo->Write();
    delete fy1Histo;
  }
  if (fx3Histo) {
    fx3Histo->Write();
    delete fx3Histo;
  }
  if (fy3Histo) {
    fy3Histo->Write();
    delete fy3Histo;
  }
  if (fdxdzHisto) {
    fdxdzHisto->Write();
    delete fdxdzHisto;
  }
  if (fdydzHisto) {
    fdydzHisto->Write();
    delete fdydzHisto;
  }
  if (fdsdzHisto) {
    fdsdzHisto->Write();
    delete fdsdzHisto;
  }
  if (ftargetxHisto) {
    ftargetxHisto->Write();
    delete ftargetxHisto;
  }
  if (ftargetyHisto) {
    ftargetyHisto->Write();
    delete ftargetyHisto;
  }
  if (ftargetxyHisto) {
    ftargetxyHisto->Write();
    delete ftargetxyHisto;
  }
  if (fBC3xyHisto) {
    fBC3xyHisto->Write();
    delete fBC3xyHisto;
  }
  if (fBC3xdxdzHisto) {
    fBC3xdxdzHisto->Write();
    delete fBC3xdxdzHisto;
  }
  if (fBC3ydydzHisto) {
    fBC3ydydzHisto->Write();
    delete fBC3ydydzHisto;
  }
  if (fChi2Histo) {
    fChi2Histo->Write();
    delete fChi2Histo;
  }
  if (fRChi2Histo) {
    fRChi2Histo->Write();
    delete fRChi2Histo;
  }
  if (fFinalTrkHisto) {
    fFinalTrkHisto->Write();
    delete fFinalTrkHisto;
  }

  if (fNtTree) {
    fNtTree->Write();
    delete fNtTree;
    delete fNtuple;
  }
}

---

Member Function Documentation

JobCResult BCReco::Ana

(

const EDMEventHandle &

evt

)

[virtual]

Reimplemented from JobCModule.

Definition at line 651 of file BCReco.cxx.

References JobCModule::kPassed.

{
  return JobCModule::kPassed;
}

void BCReco::CleanUpTrk

(

)

[private]

Definition at line 1429 of file BCReco.cxx.

References BCReco::TrackFit_t::chi2, fNTrack, fTrackFit, fWidthTrks, BCReco::TrackFit_t::x1, BCReco::TrackFit_t::x3, BCReco::TrackFit_t::y1, and BCReco::TrackFit_t::y3.

Referenced by Reco().

{
  //====================================================================
  // A clean up step is performed. This step loops over all trajectory
  // candidates and compares the parameters, if the parameters fall
  // within some window (fWidthTrks) the two candidates are declared to
  // be the same and the candidate with the best chi-squired is chosen.
  //====================================================================
  
  std::list<TrackFit_t*>::iterator trackFitItr(fTrackFit.begin());
  std::list<TrackFit_t*>::iterator trackFitItrEnd(fTrackFit.end());  
  std::list<TrackFit_t*>::iterator trackFitItr1;
  for ( ; trackFitItr != trackFitItrEnd; ++trackFitItr) {
    trackFitItr1 = trackFitItr;
    for ( trackFitItr1++; trackFitItr1 != trackFitItrEnd; 
      ++trackFitItr1) {
      
      TrackFit_t* t1 = *trackFitItr;
      TrackFit_t* t2 = *trackFitItr1;
      if ( fabs(t1->x1 - t2->x1) < fWidthTrks &&
       fabs(t1->y1 - t2->y1) < fWidthTrks &&
       fabs(t1->x3 - t2->x3) < fWidthTrks &&
       fabs(t1->y3 - t2->y3) < fWidthTrks ) {
    // One of the tracks needs to be deleted

    fNTrack--;
    if (t1->chi2 <= t2->chi2) {
      delete t2;
      std::list<TrackFit_t*>::iterator saveItr(trackFitItr1);
      saveItr--;
      fTrackFit.erase(trackFitItr1);
      trackFitItr1 = saveItr;
    }
    else {
      delete t1;
      std::list<TrackFit_t*>::iterator saveItr(trackFitItr);
      saveItr++;
      fTrackFit.erase(trackFitItr);
      if (saveItr != fTrackFit.begin()) saveItr--;
      trackFitItr = saveItr;
      break;
    }
      }
    }
  }

  if (fNTrack>1) {
    trackFitItr = fTrackFit.begin();
    trackFitItrEnd = fTrackFit.end();
    for ( ; trackFitItr != trackFitItrEnd; ++trackFitItr) {
      (*trackFitItr)->goldenTrk = 3;
    }
  }

  if (fDebug) fnTrackHisto->Fill(fNTrack);

  // Fill histograms for dx/dz && dy/dz
  if (fDebug && fNTrack == 1) {
    float x1, y1;
    float x3, y3;
    float z1, z3;
    float dx, dy, dz;

    z1 = (GDCGeom::Id(GDCGeom::kBC1).GetZPos(2) 
      + GDCGeom::Id(GDCGeom::kBC1).GetZPos(3))/2.;
    z3 = (GDCGeom::Id(GDCGeom::kBC3).GetZPos(2) 
      + GDCGeom::Id(GDCGeom::kBC3).GetZPos(3))/2.;

    TrackFit_t* fitTrack = *fTrackFit.begin();
    x1 = fitTrack->x1;
    y1 = fitTrack->y1;
    x3 = fitTrack->x3;
    y3 = fitTrack->y3;

    dx = x3 - x1;
    dy = y3 - y1;
    dz = z3 - z1;

    fx1Histo->Fill(x1);
    fy1Histo->Fill(y1);
    fx3Histo->Fill(x3);
    fy3Histo->Fill(y3);

    fdxdzHisto->Fill(dx/dz*1000);
    fdydzHisto->Fill(dy/dz*1000);
    fdsdzHisto->Fill(sqrt(dx/dz*dx/dz + dy/dz*dy/dz + 1));
    fChi2Histo->Fill(fitTrack->chi2);
    fRChi2Histo->Fill(fitTrack->chi2/fitTrack->nHit);

    ftargetxHisto->Fill(x3+dx/dz*2000); // x position at target
    ftargetyHisto->Fill(y3+dy/dz*2000); // y position at target
    ftargetxyHisto->Fill(x3+dx/dz*2000, y3+dy/dz*2000); // position at target
    fBC3xyHisto->Fill(x3, y3);
    fBC3xdxdzHisto->Fill(x3, dx/dz*1000);
    fBC3ydydzHisto->Fill(y3, dy/dz*1000);

    for (int i=0; i<GDCGeom::kNBCPlane; ++i) 
      fFinalTrkHisto->Fill(gsResidual[i]);
  }

  if (fDebug>1) {
    trackFitItr=fTrackFit.begin();
    trackFitItrEnd = fTrackFit.end();  
    for (; trackFitItr != trackFitItrEnd; ++trackFitItr) {
      std::cout << "Track Fit " 
        << "x1 = " << (*trackFitItr)->x1 << " "
        << "y1 = " << (*trackFitItr)->y1 << " "
        << "x3 = " << (*trackFitItr)->x3 << " "
        << "y3 = " << (*trackFitItr)->y3 << std::endl;
    }
  }
}

void BCReco::Clear

(

)

[private]

Definition at line 1612 of file BCReco.cxx.

References fBCAHitList, fBCHitList, fCandTraj, fCandTrk, fTrackFit, and kNBC.

Referenced by Reco().

{
//======================================================================
// Delete all objects that got stored in vectors and clear all vectors
//======================================================================
  fNTrack = 0;

  for (int i = 0; i < GDCGeom::kNBC; ++i) {
    for (int j = 0; j < GDCGeom::kNPlane; ++j) {
      std::vector<Hit*>::iterator itr(fBCHitList[i][j].begin());
      std::vector<Hit*>::iterator itrEnd(fBCHitList[i][j].end());
      for ( ; itr != itrEnd; ++itr) delete (*itr);
      fBCHitList[i][j].clear();

      std::list<AssocHit_t*>::iterator itrA(fBCAHitList[i][j].begin());
      std::list<AssocHit_t*>::iterator itrAEnd(fBCAHitList[i][j].end());
      for ( ; itrA != itrAEnd; ++itrA) delete (*itrA);
      fBCAHitList[i][j].clear();
    }
  }

  for (int i = 0; i < GDCGeom::kNPlane; ++i) {
    std::list<CandPair_t*>::iterator itr(fCandTrk[i].begin());
    std::list<CandPair_t*>::iterator itrEnd(fCandTrk[i].end());
    for ( ; itr != itrEnd; ++itr) delete (*itr);
    fCandTrk[i].clear();
  }

  std::list<CandTrajectory_t*>::iterator itrC(fCandTraj.begin());
  std::list<CandTrajectory_t*>::iterator itrCEnd(fCandTraj.end());
  for ( ; itrC != itrCEnd; ++itrC) delete (*itrC);
  fCandTraj.clear();

  std::list<TrackFit_t*>::iterator itrT(fTrackFit.begin());
  std::list<TrackFit_t*>::iterator itrTEnd(fTrackFit.end());
  for ( ; itrT != itrTEnd; ++itrT) delete (*itrT);
  fTrackFit.clear();
}

void BCReco::DoCrossings	(	float	z1,
		float	z2,
		float	z3
	)			[private]

Definition at line 534 of file BCReco.cxx.

References BCReco::TrackFit_t::chi2, DistanceWireLine(), fBCHitList, fNTrack, fTrackFit, fWidth_TrkFit, BCCrossing::GetCrossing(), BCReco::TrackFit_t::goldenTrk, gsResidual, gsWire, kBC1, BCReco::TrackFit_t::nHit, BCReco::TrackFit_t::wire, x, BCReco::TrackFit_t::x1, BCReco::TrackFit_t::x3, x3, y, BCReco::TrackFit_t::y1, y1, BCReco::TrackFit_t::y3, and y3.

Referenced by Reco().

{
  int nTmp = 0;
  float x1, x3, y1, y3;
  float x[3],y[3];
  int w[3] = {0,0,0};
  for (int ichamber=0; ichamber<3; ichamber++) {
    std::list<Crossing*> crossing(0);
    BCCrossing c;
  
    bool hasCrossing = c.GetCrossing(ichamber, fBCHitList[ichamber],
                     crossing);
    
    if (hasCrossing) {
      if (crossing.size() == 1) {
      Crossing* cross = crossing.front();
      x[ichamber] = cross->x();
      y[ichamber] = cross->y();
      w[ichamber] = cross->w();
      nTmp++;
      }
    }
  }
  
  if (nTmp>=2) {
    if (w[0] == 0) {
      x1 = x[1] + (x[2]-x[1])/(z3-z2)*(z1-z2);
      y1 = y[1] + (y[2]-y[1])/(z3-z2)*(z1-z2);
      x3 = x[2];
      y3 = y[2];
    }
    else if (w[1] == 0) {
      x1 = x[0];
      y1 = y[0];
      x3 = x[2];
      y3 = y[2];
    }
    else if (w[2] == 0) {
      x1 = x[0];
      y1 = y[0];
      x3 = x[0] + (x[1]-x[0])/(z2-z1)*(z3-z1);
      y3 = y[0] + (y[1]-y[0])/(z2-z1)*(z3-z1);  
    }
    else if (w[1] <= w[0] && w[1] <= w[2]) {
      x1 = x[0];
      y1 = y[0];
      x3 = x[2];
      y3 = y[2];
    }
    else if (w[0] < w[1] && w[0] <= w[2]) {
      x1 = x[1] + (x[2]-x[1])/(z3-z2)*(z1-z2);
      y1 = y[1] + (y[2]-y[1])/(z3-z2)*(z1-z2);
      x3 = x[2];
      y3 = y[2];
    }
    else if (w[2] < w[0] && w[2] < w[1]) {
      x1 = x[0];
      y1 = y[0];
      x3 = x[0] + (x[1]-x[0])/(z2-z1)*(z3-z1);
      y3 = y[0] + (y[1]-y[0])/(z2-z1)*(z3-z1);
    }

    // loop thru all raw hits and find the wire closest to the track
    for (int ichamber=0; ichamber<3; ichamber++)
      for (int iplane=0; iplane<4; iplane++) {
    std::vector<Hit*>::iterator 
      hitItr(fBCHitList[ichamber][iplane].begin());
    std::vector<Hit*>::iterator 
      hitItrEnd(fBCHitList[ichamber][iplane].end());
    
    double tmpD = 9999;
    gsWire[iplane+ichamber*4] = -1;
    gsResidual[iplane+ichamber*4] = 9999.;
    for ( ; hitItr != hitItrEnd; ++hitItr) {
      int wire = (*hitItr)->wire();
      double d = DistanceWireLine(iplane+ichamber*4,wire,x1,y1,x3,y3);
      
      if (d<tmpD && d<fWidth_TrkFit) {
        tmpD = d;
        gsWire[iplane+ichamber*4] = wire;
        gsResidual[iplane+ichamber*4] = d;
      }
    }
      }
    
    TrackFit_t* trkfit = new TrackFit_t;
    int        nHit = 0;
    double     chi2 = 0.0;
    
    trkfit->x1        = x1;
    trkfit->y1        = y1;
    trkfit->x3        = x3;
    trkfit->y3        = y3;
    
    for (int i=0; i<GDCGeom::kNBCPlane; i++) {
      trkfit->wire[i]    = gsWire[i];
      if (gsWire[i]>0) {
    nHit++;
    double d = DistanceWireLine(i,gsWire[i],trkfit->x1,trkfit->y1,
                    trkfit->x3, trkfit->y3);
    chi2 =+ d*d;
      }
    }
    trkfit->chi2 = chi2/nHit
      /(GDCGeom::Id(GDCGeom::kBC1).GetWrSep(1)
    *GDCGeom::Id(GDCGeom::kBC1).GetWrSep(1)/12);
    
    trkfit->nHit      = nHit;
    trkfit->goldenTrk = 2;
    
    fTrackFit.push_back(trkfit);
    fNTrack++;
  }
}

void BCReco::FillNt

(

)

[private]

Definition at line 1544 of file BCReco.cxx.

References BCRecoNt::chi2, fNtTree, fNtuple, fTrackFit, BCRecoNt::goldenTrk, BCRecoNt::nHit, BCRecoNt::ntrack, BCRecoNt::wire, BCRecoNt::x1, BCRecoNt::x3, BCRecoNt::y1, and BCRecoNt::y3.

Referenced by Reco().

{
//======================================================================
// Store track information in the event
// ONLY if total number of tracks < 10 !!!
//======================================================================

  std::list<TrackFit_t*>::iterator trackFitItr(fTrackFit.begin());
  std::list<TrackFit_t*>::iterator trackFitItrEnd(fTrackFit.end());

  int iTrack = 0;

  for ( ; trackFitItr != trackFitItrEnd; ++trackFitItr) {
    fNtuple->nHit[iTrack]      = (*trackFitItr)->nHit;
    fNtuple->x1[iTrack]        = (*trackFitItr)->x1;
    fNtuple->y1[iTrack]        = (*trackFitItr)->y1;
    fNtuple->x3[iTrack]        = (*trackFitItr)->x3;
    fNtuple->y3[iTrack]        = (*trackFitItr)->y3;
    fNtuple->chi2[iTrack]      = (*trackFitItr)->chi2;
    fNtuple->goldenTrk[iTrack] = (*trackFitItr)->goldenTrk;
    for (int i=0; i<GDCGeom::kNBCPlane; i++) {
      fNtuple->wire[iTrack][i] = (*trackFitItr)->wire[i];
    }

    ++iTrack;
  }
  
  fNtuple->ntrack = iTrack;
  
  for (int j=iTrack; j<10; j++) {
    fNtuple->nHit[j]      = 0;
    fNtuple->x1[j]        = 9999.;
    fNtuple->y1[j]        = 9999.;
    fNtuple->x3[j]        = 9999.;
    fNtuple->y3[j]        = 9999.;
    fNtuple->chi2[j]      = 9999.;
    fNtuple->goldenTrk[j] = 9999;
    for (int i=0; i<GDCGeom::kNBCPlane; i++) fNtuple->wire[j][i] = -1;
  }
  
  if (iTrack==0) {
    fNtuple->nHit[iTrack] = 0;
    fNtuple->chi2[iTrack] = 9999.;
    for (int i=0; i<GDCGeom::kNBCPlane; i++) fNtuple->wire[iTrack][i] = -1;
  }
  
  if (fNtTree) fNtTree->Fill();
}

void BCReco::FitLSQ

(

)

[private]

Definition at line 1237 of file BCReco.cxx.

References alpha, DistanceWireLine(), fBCHitList, fCandTraj, fNTrack, fNWireItr, fTrackFit, fWidth_TrkFit, gsResidual, gsWire, kBC1, kBC2, kBC3, kNBC, p, BCReco::CandTrajectory_t::x1, BCReco::CandTrajectory_t::x3, x3, BCReco::CandTrajectory_t::y1, y1, BCReco::CandTrajectory_t::y3, and y3.

Referenced by TrkFit().

{
  // Takes the candidate wire pairs from view pairings and performs a
  // least-squared fit with the 12 planes of wires to four parameters.

  std::list<CandTrajectory_t*>::iterator candTrajItr(fCandTraj.begin());
  std::list<CandTrajectory_t*>::iterator candTrajItrEnd(fCandTraj.end());

  double x1, y1;  // space points in BC1
  double x3, y3;  // space points in BC3

  bool noFit = false;
  
  const GDCGeom* geom[GDCGeom::kNBC] = {
    &GDCGeom::Id(GDCGeom::kBC1),
    &GDCGeom::Id(GDCGeom::kBC2),
    &GDCGeom::Id(GDCGeom::kBC3)
  };


  static float z1 = (geom[0]->GetZPos(2) + 
             geom[0]->GetZPos(3))/2.;  // middle pos. of chamber 1
  static float z3 = (geom[2]->GetZPos(2) + 
             geom[2]->GetZPos(3))/2.;  // middle pos. of chamber 3 
  
  float sep[GDCGeom::kNBCPlane]; // wire spacing (cm)
  float d1w[GDCGeom::kNBCPlane]; // distance of wire 1 from the origin (cm)
  float c[GDCGeom::kNBCPlane][4]; // geometric constants (angles of wires?)
  double coordinate[4];  // dummy array to store init. and final pos. of trk
  
  // Loop over every trajectory and fit it, store results into
  // fTrackFit list
  for ( ; candTrajItr != candTrajItrEnd; ++candTrajItr) {
    TrackFit_t* trackFit = new TrackFit_t;
    CandTrajectory_t* candTraj = *candTrajItr;
    ++fNTrack;
    
    for (int nItr=0; nItr<fNWireItr; nItr++) {
      x1 = candTraj->x1;
      y1 = candTraj->y1;
      x3 = candTraj->x3;
      y3 = candTraj->y3;
      
      // loop thru all raw hits and find the wire closest to the track
      for (int ichamber=0; ichamber<GDCGeom::kNBC; ichamber++) {
    for (int iplane=0; iplane<GDCGeom::kNPlane; iplane++) {
      int planeIndx = iplane + ichamber * GDCGeom::kNPlane;

      std::vector<Hit*>::iterator 
        hitItr(fBCHitList[ichamber][iplane].begin());
      std::vector<Hit*>::iterator 
        hitItrEnd(fBCHitList[ichamber][iplane].end());
      
      double tmpD = 9999;
      gsWire[planeIndx] = -1;
      gsResidual[planeIndx] = 9999.;
      for ( ; hitItr != hitItrEnd; ++hitItr) {
        int wire = (*hitItr)->wire();
        double d = DistanceWireLine(planeIndx,wire,x1,y1,x3,y3);
        
        if (d<tmpD && d<fWidth_TrkFit) {
          tmpD = d;
          gsWire[planeIndx]    = wire;
          gsResidual[planeIndx] = d;
        }
      }
    }
      }    
  
      /*
      // loop thru associated hits and ...????
      for (int i=0; i<GDCGeom::kNBCPlane; i++) {
      if (candTraj->ahit[i]!=NULL) {
      if (candTraj->ahit[i]->twoHits == false)
      gsWire[i] = int(candTraj->ahit[i]->aHitWire);
      else {
      int hit1 = candTraj->ahit[i]->hit1->wire();
      int hit2 = candTraj->ahit[i]->hit2->wire();
      double d1 = 9999.;
      double d2 = 9999.;
      if (hit1>0) d1 = DistanceWireLine(i,hit1,x1,y1,x3,y3);
      if (hit2>0) d2 = DistanceWireLine(i,hit2,x1,y1,x3,y3);
      if (d1<d2)
      gsWire[i] = hit1;
      else
      gsWire[i] = hit2;
      }
      }
      else
      gsWire[i] = -1;
      }
      */
      for (int ichamber = 0; ichamber < GDCGeom::kNBC; ichamber++) {
    for (int iplane=0; iplane < GDCGeom::kNPlane; iplane++) {
      int planeIndx = iplane + (ichamber * 4);

      float alpha = geom[ichamber]->GetWrAng(iplane+1) - M_PI/2.;
      
      sep[planeIndx] = geom[ichamber]->GetWrSep(iplane+1);
      
      d1w[planeIndx] = geom[ichamber]->GetD1W(iplane+1);
      
      float sa = sin(alpha);
      float ca = cos(alpha);
      
      float z = geom[ichamber]->GetZPos(iplane+1);
      c[planeIndx][0] = (z3 - z)*ca/(z1 - z3)/sep[planeIndx];
      c[planeIndx][1] = (z3 - z)*sa/(z1 - z3)/sep[planeIndx];
      c[planeIndx][2] = (z - z1)*ca/(z1 - z3)/sep[planeIndx];
      c[planeIndx][3] = (z - z1)*sa/(z1 - z3)/sep[planeIndx];
    }
      }
      
      int weight[12];
      for (int i=0; i < GDCGeom::kNBCPlane; i++) {
    if (gsWire[i]>0) weight[i] = 1;
    else             weight[i] = 0;
      }
      
      double array[16] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
      double wire[4] = {0,0,0,0};
      for (int i = 0; i < GDCGeom::kNBCPlane; i++) {
    for (int j = 0; j < 16; ++j) {
      array[j] += weight[i] * c[i][j/4] * c[i][j%4];
    }

    for (int j = 0; j < 4; ++j) {
      // Subtract 1. since wire number starts from 1
      wire[j] += weight[i]*c[i][j]*
        (gsWire[i] - 1. - d1w[i]/sep[i]);
    }
      }
      
      TMatrixD cc(4,4);
      cc.SetMatrixArray(array);
      
      double det;
      cc.Invert(&det);

      if (det == 0.) {
    noFit = true;
    break;
      }

      TMatrixD w(4,1);
      w.SetMatrixArray(wire);
      
      TMatrixD p(4,1);
      p = cc*w;
      
      p.GetMatrix2Array(coordinate);
      
      candTraj->x1 = coordinate[0];
      candTraj->y1 = coordinate[1];
      candTraj->x3 = coordinate[2];
      candTraj->y3 = coordinate[3];
    }
    
    if (noFit == true) break;
    
    trackFit->x1 = coordinate[0];
    trackFit->y1 = coordinate[1];
    trackFit->x3 = coordinate[2];
    trackFit->y3 = coordinate[3];
    
    int nHit = 0;
    for (int i=0; i<GDCGeom::kNBCPlane; i++) {
      trackFit->wire[i] = gsWire[i];
      if (gsWire[i]>0) nHit++;
    }
    trackFit->nHit = nHit;
    
    // Calculate chi**2
    double chi2 = 0.0;
    
    for (int i=0; i < GDCGeom::kNBCPlane; ++i) {
      if (gsWire[i] > 0) {
    double d = DistanceWireLine(i,gsWire[i],trackFit->x1,trackFit->y1,
                    trackFit->x3, trackFit->y3);
    chi2 =+ d*d;
      }
    }
    trackFit->chi2 = chi2/nHit
      /(geom[0]->GetWrSep(1) * geom[0]->GetWrSep(1)/12);
    
    trackFit->goldenTrk = 1;
    fTrackFit.push_back(trackFit);
  }
}

void BCReco::FitTMinuit

(

)

[private]

Definition at line 1137 of file BCReco.cxx.

References BCReco::TrackFit_t::chi2, DistanceWireLine(), fBCHitList, fCandTraj, fMaxItr, fMinuit, fNTrack, fTrackFit, fWidth_TrkFit, BCReco::TrackFit_t::goldenTrk, gsResidual, gsWire, kNBC, BCReco::TrackFit_t::nHit, BCReco::TrackFit_t::wire, BCReco::TrackFit_t::x1, BCReco::TrackFit_t::x3, x3, BCReco::TrackFit_t::y1, y1, BCReco::TrackFit_t::y3, and y3.

Referenced by TrkFit().

{
  std::list<CandTrajectory_t*>::iterator candTrajItr(fCandTraj.begin());
  std::list<CandTrajectory_t*>::iterator candTrajItrEnd(fCandTraj.end());

  double x1, y1;  // space points in BC1
  double x3, y3;  // space points in BC3

  int ierr;
  double arglist[10];
  double x1err, y1err;
  double x3err, y3err;
  double chi2, fedm, errdef;
  int npari, nparx;
  
  for ( ; candTrajItr != candTrajItrEnd; ++candTrajItr) {
    CandTrajectory_t* candTraj = *candTrajItr;
    TrackFit_t* trackFit = new TrackFit_t;
    ++fNTrack;

    x1 = candTraj->x1;
    y1 = candTraj->y1;
    x3 = candTraj->x3;
    y3 = candTraj->y3;
    
    for (int itr=0; itr<fMaxItr; ++itr) {
      ierr = 0;
      arglist[0] = 1.;
      fMinuit->mnexcm("SET ERR", arglist, 1, ierr);
      fMinuit->mnparm(0, "x1", x1, 0.01, 0, 0, ierr);
      fMinuit->mnparm(1, "y1", y1, 0.01, 0, 0, ierr);
      fMinuit->mnparm(2, "x3", x3, 0.01, 0, 0, ierr);
      fMinuit->mnparm(3, "y3", y3, 0.01, 0, 0, ierr);
      
      arglist[0] = 10.;
      arglist[1] = 1.0;
      fMinuit->mnexcm("MIGRAD", arglist, 2, ierr);
      
      fMinuit->GetParameter(0, x1, x1err);
      fMinuit->GetParameter(1, y1, y1err);
      fMinuit->GetParameter(2, x3, x3err);
      fMinuit->GetParameter(3, y3, y3err);
      
      fMinuit->mnstat(chi2, fedm, errdef, npari, nparx, ierr);
    }
    
    // loop thru all raw hits and find the wire closest to the track
    for (int ichamber=0; ichamber<GDCGeom::kNBC; ichamber++)
      for (int iplane=0; iplane<GDCGeom::kNPlane; iplane++) {
    std::vector<Hit*>::iterator 
      hitItr(fBCHitList[ichamber][iplane].begin());
    std::vector<Hit*>::iterator 
      hitItrEnd(fBCHitList[ichamber][iplane].end());
    
    double tmpD = 9999;
    gsWire[iplane+ichamber*4] = -1;
    gsResidual[iplane+ichamber*4] = 9999.;
    for ( ; hitItr != hitItrEnd; ++hitItr) {
      int wire = (*hitItr)->wire();
      double d = DistanceWireLine(iplane+ichamber*4,wire,x1,y1,x3,y3);
      
      if (d<tmpD && d<fWidth_TrkFit) {
        tmpD = d;
        gsWire[iplane+ichamber*4]    = wire;
        gsResidual[iplane+ichamber*4] = d;
      }
    }
      }
    
    /*
      for (int i=0; i<GDCGeom::kNBCPlane; i++)
      if (candTraj->ahit[i] != NULL) {
      gsWire[i]  = int(candTraj->ahit[i]->aHitWire);
      gsResidual[i] = DistanceWireLine(i,gsWire[i],x1,y1,x3,y3);
      }
      else {
      gsWire[i] = -1;
      gsResidual[i] = 9999.;
      }
    */
    
    int nHit = 0;
    for (int i=0; i<GDCGeom::kNBCPlane; ++i)
      if (gsWire[i]>0) nHit++;
    
    trackFit->nHit = nHit;
    for (int i=0; i<GDCGeom::kNBCPlane; i++)
      trackFit->wire[i] = gsWire[i];
    trackFit->x1   = x1;
    trackFit->y1   = y1;
    trackFit->x3   = x3;
    trackFit->y3   = y3;
    trackFit->chi2 = chi2;
    trackFit->goldenTrk = 1;
    fTrackFit.push_back(trackFit);
  }
}

void BCReco::HitAssoc

(

)

[private]

Definition at line 658 of file BCReco.cxx.

References BCReco::AssocHit_t::aHitWire, fBCAHitList, fBCHitList, fDebug, BCReco::AssocHit_t::hit1, BCReco::AssocHit_t::hit2, kNBC, kNPlane, and BCReco::AssocHit_t::twoHits.

Referenced by Reco().

{
//======================================================================
// Hit Association: reduce the number of first and last chamber
// combinations.  The adjacent hits within each plane are found and
// grouped together.
// If w_{i+1} = w_i+1, then only w_i+0.5 is entered into the list.
//======================================================================
  for (int ichamber=0; ichamber<GDCGeom::kNBC; ++ichamber)
    for (int iplane=0; iplane<GDCGeom::kNPlane; ++iplane) {
      std::vector<Hit*>::iterator wireItr
    (fBCHitList[ichamber][iplane].begin());
      std::vector<Hit*>::iterator wireItrEnd
    (fBCHitList[ichamber][iplane].end());

      AssocHit_t* aHitInfo = 0;
      for ( ; wireItr != wireItrEnd; ++wireItr) {
    Hit* thisHit = *wireItr;

    if (aHitInfo && 
        ((aHitInfo->aHitWire + 1) == thisHit->wire())) {
      // This wire is already has a neighbor in the list, so 
      // we just modify the previous entry
      assert(!fBCAHitList[ichamber][iplane].empty());

      aHitInfo->aHitWire += 0.5;
      aHitInfo->hit2      = thisHit;
      aHitInfo->twoHits   = true;
    }
    else {
      aHitInfo = new AssocHit_t;
      aHitInfo->aHitWire = thisHit->wire();
      aHitInfo->hit1     = thisHit;
      aHitInfo->hit2     = thisHit;
      aHitInfo->twoHits  = false;
      fBCAHitList[ichamber][iplane].push_back(aHitInfo);      
    }
      }
    }

  if (fDebug>1) {
    std::cout << "Hit Association" << std::endl;
    for (int ichamber=0; ichamber<GDCGeom::kNBC; ++ichamber)
      for (int iplane=0; iplane<GDCGeom::kNPlane; ++iplane) {
    std::list<AssocHit_t*>::iterator wireItr
      (fBCAHitList[ichamber][iplane].begin());
    std::list<AssocHit_t*>::iterator wireItrEnd
      (fBCAHitList[ichamber][iplane].end());
    for ( ; wireItr != wireItrEnd; ++wireItr) {
      AssocHit_t* thisHit = (*wireItr);
      std::cout << "Chamber " << ichamber << "  "
            << "Plane " << iplane << "  "
            << thisHit->aHitWire << " "
            << thisHit->hit1->wire() << " "
            << thisHit->hit2->wire() << " "
            << thisHit->hit1->tdc() << " "
            << thisHit->hit2->tdc() << " "
            << thisHit->twoHits << " "
            << std::endl;
    }
      }
  }
}

void BCReco::MakeHistos

(

)

[private]

......................................................................

Definition at line 231 of file BCReco.cxx.

References fBC3xdxdzHisto, fBC3xyHisto, fBC3ydydzHisto, fChi2Histo, fdsdzHisto, fdxdzHisto, fdydzHisto, fFinalTrkHisto, fnTrackHisto, fNtTree, fNtuple, fRChi2Histo, ftargetxHisto, ftargetxyHisto, ftargetyHisto, fTrkFindHisto, fViewMatchHisto, fx1Histo, fx3Histo, fy1Histo, and fy3Histo.

Referenced by BCReco().

{
  fNtTree = new TTree("fBCnt","BC Track Fit Ntuple");
  fNtTree->Branch("fBCnt", "BCRecoNt", &fNtuple);
  fTrkFindHisto[0] = new TH1F("fTrkFindHisto_u",
                  "Predicted wire - Hit wire (u-view)",
                  640, -320., 320.);
  fTrkFindHisto[1] = new TH1F("fTrkFindHisto_v",
                  "Predicted wire - Hit wire (v-view)",
                  640, -320., 320.);
  fTrkFindHisto[2] = new TH1F("fTrkFindHisto_s",
                  "Predicted wire - Hit wire (s-view)",
                  640, -320., 320.);
  fTrkFindHisto[3] = new TH1F("fTrkFindHisto_t",
                  "Predicted wire - Hit wire (t-view)",
                  640, -320., 320.);

  fViewMatchHisto[0][0] = new TH1F("fViewMatchHisto_ch1_u",
                   "Predicted wire - Hit wire (u-view)",
                   640, -320., 320.);
  fViewMatchHisto[0][1] = new TH1F("fViewMatchHisto_ch1_v",
                   "Predicted wire - Hit wire (v-view)",
                   640, -320., 320.);
  fViewMatchHisto[0][2] = new TH1F("fViewMatchHisto_ch1_s",
                   "Predicted wire - Hit wire (s-view)",
                   640, -320., 320.);
  fViewMatchHisto[0][3] = new TH1F("fViewMatchHisto_ch1_t",
                   "Predicted wire - Hit wire (t-view)",
                   640, -320., 320.);

  fViewMatchHisto[1][0] = new TH1F("fViewMatchHisto_ch2_u",
                   "Predicted wire - Hit wire (u-view)",
                   640, -320., 320.);
  fViewMatchHisto[1][1] = new TH1F("fViewMatchHisto_ch2_v",
                   "Predicted wire - Hit wire (v-view)",
                   640, -320., 320.);
  fViewMatchHisto[1][2] = new TH1F("fViewMatchHisto_ch2_s",
                   "Predicted wire - Hit wire (s-view)",
                   640, -320., 320.);
  fViewMatchHisto[1][3] = new TH1F("fViewMatchHisto_ch2_t",
                   "Predicted wire - Hit wire (t-view)",
                   640, -320., 320.);

  fViewMatchHisto[2][0] = new TH1F("fViewMatchHisto_ch3_u",
                   "Predicted wire - Hit wire (u-view)",
                   640, -320., 320.);
  fViewMatchHisto[2][1] = new TH1F("fViewMatchHisto_ch3_v",
                   "Predicted wire - Hit wire (v-view)",
                   640, -320., 320.);
  fViewMatchHisto[2][2] = new TH1F("fViewMatchHisto_ch3_s",
                   "Predicted wire - Hit wire (s-view)",
                   640, -320., 320.);
  fViewMatchHisto[2][3] = new TH1F("fViewMatchHisto_ch3_t",
                   "Predicted wire - Hit wire (t-view)",
                   640, -320., 320.);

  fnTrackHisto = new TH1F("fnTrackHisto", "Number of Track", 10, 0, 10);
  fx1Histo = new TH1F("fx1Histo","x1", 200, -10., 10.);
  fy1Histo = new TH1F("fy1Histo","y1", 200, -10., 10.);
  fx3Histo = new TH1F("fx3Histo","x3", 200, -10., 10.);
  fy3Histo = new TH1F("fy3Histo","y3", 200, -10., 10.);
  fdxdzHisto = new TH1F("fdxdzHisto","dx/dz", 200, -5., 5.);
  fdydzHisto = new TH1F("fdydzHisto","dy/dz", 200, -5., 5.);
  fdsdzHisto = new TH1F("fdsdzHisto","ds/dz", 200, 0.9999, 1.0001);
  ftargetxHisto = new TH1F("ftargetxHisto","x at target", 200, -10., 10.);
  ftargetyHisto = new TH1F("ftargetyHisto","y at target", 200, -10., 10.);
  ftargetxyHisto = new TH2F("ftargetxyHisto", "y:x", 200, -10.,10., 
                200,-10.,10.);
  fBC3xyHisto = new TH2F("fBC3xyHisto", "y:x", 
             200, -10., 10., 200,-10.,10.);
  fBC3xdxdzHisto = new TH2F("fBC3xdxdzHisto", "x:dxdz", 
                200, -10.,10., 200,-5.,5.);
  fBC3ydydzHisto = new TH2F("fBC3ydydzHisto", "y:dydz", 
                200, -10.,10., 200,-5.,5.);
  fChi2Histo = new TH1F("fChi2Histo","Chi2", 4000, 0., 2000.);
  fRChi2Histo = new TH1F("fRChi2Histo","Reduced Chi2", 400, 0., 200.);
  fFinalTrkHisto = new TH1F("fFinalTrkHisto",
                "Distance(Fit Track - Hit wire)", 
                100, 0., 0.5);
}

void BCReco::MatchSU

(

)

[private]

Definition at line 837 of file BCReco.cxx.

References BCReco::CandTrajectory_t::ahit, BCReco::AssocHit_t::aHitWire, fBCAHitList, fCandTraj, fCandTrk, fDebug, fN_ViewMatch, fViewMatchHisto, gsCotAng, gsDelX, gsXWire1, BCReco::AssocHit_t::hit1, BCReco::AssocHit_t::hit2, kBC1, kBC3, kNBC, NULL, x, BCReco::CandTrajectory_t::x1, BCReco::CandTrajectory_t::x3, x3, y, BCReco::CandTrajectory_t::y1, y1, BCReco::CandTrajectory_t::y3, and y3.

Referenced by ViewMatch().

{  
  float xu, xs, yu, ys;      // x- or y-intercepts of a wire in u & s views
  float x1, x3, y1, y3, z1, z3;
  float x, y, z;
  float w;                   // predicted wire number

  // u-s view pairings
  z1 = (GDCGeom::Id(GDCGeom::kBC1).GetZPos(1)
    + GDCGeom::Id(GDCGeom::kBC1).GetZPos(3))/2.;
  z3 = (GDCGeom::Id(GDCGeom::kBC3).GetZPos(1)
    + GDCGeom::Id(GDCGeom::kBC3).GetZPos(3))/2.;

  std::list<CandPair_t*>::iterator uCandTrkItr(fCandTrk[0].begin());
  std::list<CandPair_t*>::iterator uCandTrkItrEnd(fCandTrk[0].end());
  std::list<CandPair_t*>::iterator sCandTrkItr(fCandTrk[2].begin());
  std::list<CandPair_t*>::iterator sCandTrkItrEnd(fCandTrk[2].end());

  for ( ; uCandTrkItr != uCandTrkItrEnd; ++uCandTrkItr) {
    CandPair_t* pair0 = *uCandTrkItr;

    for (sCandTrkItr=fCandTrk[2].begin();
     sCandTrkItr != sCandTrkItrEnd; ++sCandTrkItr) {
      CandPair_t* pair2 = *sCandTrkItr;

      xu = gsXWire1[0] - gsDelX[0]*(pair0->hit1->aHitWire-1);
      xs = gsXWire1[2] - gsDelX[2]*(pair2->hit1->aHitWire-1);
      yu = -xu*gsCotAng[0];
      ys = -xs*gsCotAng[2];
      
      x1 = (ys-yu)/(gsCotAng[0] - gsCotAng[2]);
      y1 = (ys*gsCotAng[0]-yu*gsCotAng[2])/(gsCotAng[0] - gsCotAng[2]);

      xu = gsXWire1[0] - gsDelX[0]*(pair0->hit3->aHitWire-1);
      xs = gsXWire1[2] - gsDelX[2]*(pair2->hit3->aHitWire-1);
      yu = -xu*gsCotAng[0];
      ys = -xs*gsCotAng[2];

      x3 = (ys-yu)/(gsCotAng[0] - gsCotAng[2]);
      y3 = (ys*gsCotAng[0]-yu*gsCotAng[2])/(gsCotAng[0] - gsCotAng[2]);

      CandTrajectory_t* candTraj = new CandTrajectory_t;

      candTraj->x1  = x1;
      candTraj->y1  = y1;
      candTraj->x3  = x3;
      candTraj->y3  = y3;

      candTraj->ahit[0]  = pair0->hit1;
      candTraj->ahit[2]  = pair2->hit1;
      candTraj->ahit[4]  = pair0->hit2;
      candTraj->ahit[6]  = pair2->hit2;
      candTraj->ahit[8]  = pair0->hit3;
      candTraj->ahit[10] = pair2->hit3;
      candTraj->ahit[1]  = NULL;
      candTraj->ahit[3]  = NULL;
      candTraj->ahit[5]  = NULL;
      candTraj->ahit[7]  = NULL;
      candTraj->ahit[9]  = NULL;
      candTraj->ahit[11] = NULL;

      int num = 0;

      std::list<AssocHit_t*>::iterator wItr;
      std::list<AssocHit_t*>::iterator wItrEnd;

      for (int ichamber=0; ichamber< GDCGeom::kNBC; ++ichamber) {
    const GDCGeom& geomID = GDCGeom::Id(ichamber);
    for (int iplane=1; iplane<4; iplane=iplane+2) {
      bool hasHit = false;

      wItr = fBCAHitList[ichamber][iplane].begin();
      wItrEnd = fBCAHitList[ichamber][iplane].end();

      z = geomID.GetZPos(iplane+1); 
      x = x1 + (x3-x1)/(z3-z1)*(z-z1);
      y = y1 + (y3-y1)/(z3-z1)*(z-z1);

      w = (y/gsCotAng[iplane]-x+gsXWire1[iplane])/gsDelX[iplane]+1;

      float minD = 320.;
      float diff = 320.;
      for ( ; wItr != wItrEnd; ++wItr) {
        AssocHit_t* thisHit = *wItr;

        if ((w-fN_ViewMatch)<= thisHit->aHitWire 
        && (w+fN_ViewMatch)>= thisHit->aHitWire) {
          hasHit = true;
          if (diff < fabs(w-thisHit->aHitWire) 
          || wItr==fBCAHitList[ichamber][iplane].begin()) {
        diff = fabs(w-thisHit->aHitWire);

        if (iplane==1)
          switch (ichamber) {
          case 0: candTraj->ahit[1] = thisHit; break;
          case 1: candTraj->ahit[5] = thisHit; break;
          case 2: candTraj->ahit[9] = thisHit; break;
          }
        else
          switch (ichamber) {
          case 0: candTraj->ahit[3]  = thisHit; break;
          case 1: candTraj->ahit[7]  = thisHit; break;
          case 2: candTraj->ahit[11] = thisHit; break;
          }
          }
        }
        if (fabs(w-thisHit->aHitWire) < minD) {
          minD = w-thisHit->aHitWire;
        }
      }
      if (fViewMatchHisto[ichamber][iplane] && minD!=320.)
        fViewMatchHisto[ichamber][iplane]->Fill(minD);

      if (!hasHit) ++num;
      if (!fDebug && num>1) break;
    }
      }

      if (num<=1) {
    fCandTraj.push_back(candTraj);

    if (fDebug>2) {
      std::cout << "View Match 1 " 
            << "x1 = " << candTraj->x1 << " "
            << "y1 = " << candTraj->y1 << " "
            << "x3 = " << candTraj->x3 << " "
            << "y3 = " << candTraj->y3 << std::endl;
      for (int i=0; i<GDCGeom::kNBCPlane; i++)
        if (candTraj->ahit[i]!=NULL)
          std::cout << candTraj->ahit[i]->aHitWire << " ";
      std::cout << std::endl;
    }
      }
      else delete candTraj;
    }
  }
}

void BCReco::MatchTV

(

)

[private]

Definition at line 977 of file BCReco.cxx.

References BCReco::CandTrajectory_t::ahit, BCReco::AssocHit_t::aHitWire, fBCAHitList, fCandTraj, fCandTrk, fDebug, fN_ViewMatch, fViewMatchHisto, gsCotAng, gsDelX, gsXWire1, BCReco::AssocHit_t::hit1, BCReco::AssocHit_t::hit2, kBC1, kBC3, kNBC, NULL, x, BCReco::CandTrajectory_t::x1, BCReco::CandTrajectory_t::x3, x3, y, BCReco::CandTrajectory_t::y1, y1, BCReco::CandTrajectory_t::y3, and y3.

Referenced by ViewMatch().

                     {
  float xt, xv, yt, yv;      // x- or y-intercepts of a wire in t & u views
  float x1, x3, y1, y3, z1, z3;
  float x, y, z;
  float w;                   // predicted wire number

  // v-t view pairings
  z1 = (GDCGeom::Id(GDCGeom::kBC1).GetZPos(2)
    + GDCGeom::Id(GDCGeom::kBC1).GetZPos(4))/2.;
  z3 = (GDCGeom::Id(GDCGeom::kBC3).GetZPos(2)
    + GDCGeom::Id(GDCGeom::kBC3).GetZPos(4))/2.;

  std::list<CandPair_t*>::iterator vCandTrkItr(fCandTrk[1].begin());
  std::list<CandPair_t*>::iterator vCandTrkItrEnd(fCandTrk[1].end());
  std::list<CandPair_t*>::iterator tCandTrkItr(fCandTrk[3].begin());
  std::list<CandPair_t*>::iterator tCandTrkItrEnd(fCandTrk[3].end());

  for ( ; vCandTrkItr != vCandTrkItrEnd; ++vCandTrkItr) {
    CandPair_t* pair1 = *vCandTrkItr;

    for (tCandTrkItr=fCandTrk[3].begin();
     tCandTrkItr != tCandTrkItrEnd; ++tCandTrkItr) {
      CandPair_t* pair3 = *tCandTrkItr;

      xv = gsXWire1[1] - gsDelX[1]*(pair1->hit1->aHitWire-1);
      xt = gsXWire1[3] - gsDelX[3]*(pair3->hit1->aHitWire-1);
      yv = -xv*gsCotAng[1];
      yt = -xt*gsCotAng[3];
      
      x1 = (yt-yv)/(gsCotAng[1] - gsCotAng[3]);
      y1 = (yt*gsCotAng[1]-yv*gsCotAng[3])/(gsCotAng[1] - gsCotAng[3]);

      xv = gsXWire1[1] - gsDelX[1]*(pair1->hit3->aHitWire-1);
      xt = gsXWire1[3] - gsDelX[3]*(pair3->hit3->aHitWire-1);
      yv = -xv*gsCotAng[1];
      yt = -xt*gsCotAng[3];

      x3 = (yt-yv)/(gsCotAng[1] - gsCotAng[3]);
      y3 = (yt*gsCotAng[1]-yv*gsCotAng[3])/(gsCotAng[1] - gsCotAng[3]);
      
      CandTrajectory_t* candTraj = new CandTrajectory_t;

      candTraj->x1  = x1;
      candTraj->y1  = y1;
      candTraj->x3  = x3;
      candTraj->y3  = y3;

      candTraj->ahit[0]  = NULL;
      candTraj->ahit[2]  = NULL;
      candTraj->ahit[4]  = NULL;
      candTraj->ahit[6]  = NULL;
      candTraj->ahit[8]  = NULL;
      candTraj->ahit[10] = NULL;
      candTraj->ahit[1]  = pair1->hit1;
      candTraj->ahit[3]  = pair3->hit1;
      candTraj->ahit[5]  = pair1->hit2;
      candTraj->ahit[7]  = pair3->hit2;
      candTraj->ahit[9]  = pair1->hit3;
      candTraj->ahit[11] = pair3->hit3;

      int num = 0;

      std::list<AssocHit_t*>::iterator wItr;
      std::list<AssocHit_t*>::iterator wItrEnd;

      for (int ichamber=0; ichamber< GDCGeom::kNBC; ++ichamber) {
    const GDCGeom& geomID = GDCGeom::Id(ichamber);
    for (int iplane=0; iplane<GDCGeom::kNPlane; iplane=iplane+2) {
      bool hasHit = false;

      z = geomID.GetZPos(iplane+1); 
      x = x1 + (x3-x1)/(z3-z1)*(z-z1);
      y = y1 + (y3-y1)/(z3-z1)*(z-z1);

      w = (y/gsCotAng[iplane]-x+gsXWire1[iplane])/gsDelX[iplane]+1;

      float minD = 320.;
      float diff = 320;

      wItr = fBCAHitList[ichamber][iplane].begin();
      wItrEnd = fBCAHitList[ichamber][iplane].end();
      for ( ; wItr != wItrEnd; ++wItr) {
        AssocHit_t* thisHit = *wItr;
        if ((w-fN_ViewMatch) <= thisHit->aHitWire
        &&(w+fN_ViewMatch) >= thisHit->aHitWire) {        
          hasHit = true;
          if (diff < fabs(w-thisHit->aHitWire)
          || wItr==fBCAHitList[ichamber][iplane].begin()) {
        diff = fabs(w-thisHit->aHitWire);
        if (iplane==0)
          switch (ichamber) {
          case 0: candTraj->ahit[0] = thisHit; break;
          case 1: candTraj->ahit[4] = thisHit; break;
          case 2: candTraj->ahit[8] = thisHit; break;
          }
        else
          switch (ichamber) {
          case 0: candTraj->ahit[2]  = thisHit; break;
          case 1: candTraj->ahit[6]  = thisHit; break;
          case 2: candTraj->ahit[10] = thisHit; break;
          }
          }
        }
        if (fabs(w-thisHit->aHitWire) < minD) {
          minD = w-thisHit->aHitWire;
        }
      }
      if (fViewMatchHisto[ichamber][iplane] && minD!=320.)
        fViewMatchHisto[ichamber][iplane]->Fill(minD);

      if (!hasHit) ++num;
      if (!fDebug && num>1) break;
    }
      }
 
      if (num<=1) {
    fCandTraj.push_back(candTraj);

    if (fDebug>2) {
      std::cout << "View Match 2 " 
            << "x1 = " << candTraj->x1 << " "
            << "y1 = " << candTraj->y1 << " "
            << "x3 = " << candTraj->x3 << " "
            << "y3 = " << candTraj->y3 << std::endl;
      for (int i=0; i<GDCGeom::kNBCPlane; i++)
        if (candTraj->ahit[i]!=NULL)
          std::cout << candTraj->ahit[i]->aHitWire << " ";
      std::cout << std::endl;
    }
      }
      else delete candTraj;
    }
  }
}

JobCResult BCReco::Reco

(

EDMEventHandle &

evt

)

[virtual]

Reimplemented from JobCModule.

Definition at line 406 of file BCReco.cxx.

References EDMEventHandle::Cal(), CleanUpTrk(), Clear(), DoCrossings(), EDMEventHeader::Event(), BCRecoNt::evt, fBCHitList, fCandTraj, FillNt(), fMaxNWirePerChamber, fMaxTotalNWire, fNTrack, fNtuple, fTrackFit, EDMDataBucket::Get(), EDMDataBucket::GetFolder(), EDMEventHandle::Header(), HitAssoc(), kBC1, kBC2, kBC3, JobCModule::kFailed, kNBC, kNPlane, JobCModule::kPassed, EDMDataBucket::MakeFolder(), EDMDataBucket::Put(), EDMEventHandle::Reco(), EDMEventHeader::Run(), BCRecoNt::run, BCRecoNt::sub, EDMEventHeader::SubRun(), TrkFind(), TrkFit(), ViewMatch(), x3, y1, and y3.

{

  if (evt.Reco().GetFolder("./bc")==0) {
    evt.Reco().MakeFolder("./bc");
  }
  if (evt.Reco().GetFolder("./bc/linefit")==0) {
    evt.Reco().MakeFolder("./bc/linefit");
  }

  // Record header information to ntuple
  fNtuple->run = evt.Header().Run();
  fNtuple->sub = evt.Header().SubRun();
  fNtuple->evt = evt.Header().Event();

  // Extract wire hits from the event
  const char* dir1 = "./bc1";
  const char* dir2 = "./bc2";
  const char* dir3 = "./bc3";

  int nWire = 0;

  // Use Calibrated data
  std::vector<const DCWire*> bc1digit;
  std::vector<const DCWire*> bc2digit;
  std::vector<const DCWire*> bc3digit;

  try { 
    evt.Cal().Get(dir1, bc1digit); 
    if ((int) bc1digit.size() > fMaxNWirePerChamber) return kFailed;
    nWire += bc1digit.size();
  }
  catch (EDMException e) { return JobCModule::kFailed; }
  try { 
    evt.Cal().Get(dir2, bc2digit); 
    if ((int) bc2digit.size() > fMaxNWirePerChamber) return kFailed;
    nWire += bc2digit.size();
  }
  catch (EDMException e) { return JobCModule::kFailed; }
  try { 
    evt.Cal().Get(dir3, bc3digit); 
    if ((int) bc3digit.size() > fMaxNWirePerChamber) return kFailed;
    nWire += bc3digit.size();
  }
  catch (EDMException e) { return JobCModule::kFailed; }

  if (nWire > fMaxTotalNWire) return kFailed;
  
  for (int i=0; i<int(bc1digit.size()); ++i) {
    short t = bc1digit[i]->GetTime();
    if (t>-20 && t<40) {
      Hit* hit = new Hit(bc1digit[i]->GetWire(), t);
      fBCHitList[0][bc1digit[i]->GetPlane()-1].push_back(hit);
    }
  }

  for (int i=0; i<int(bc2digit.size()); ++i) {
    short t = bc2digit[i]->GetTime();
    if (t>-20 && t<40) {
      Hit* hit = new Hit(bc2digit[i]->GetWire(), t);
      fBCHitList[1][bc2digit[i]->GetPlane()-1].push_back(hit);
    }
  }

  for (int i=0; i<int(bc3digit.size()); ++i) {
    short t = bc3digit[i]->GetTime();
    if (t>-20 && t<40) {
      Hit* hit = new Hit(bc3digit[i]->GetWire(), t);
      fBCHitList[2][bc3digit[i]->GetPlane()-1].push_back(hit);
    }
  }

  for (int ichamber=0; ichamber<GDCGeom::kNBC; ++ichamber) {
    for (int iplane=0; iplane<GDCGeom::kNPlane; ++iplane) {
      sort(fBCHitList[ichamber][iplane].begin(), 
       fBCHitList[ichamber][iplane].end(),
       HitComp());
    }
  }

  HitAssoc();
  TrkFind();
  ViewMatch();
  TrkFit();
  CleanUpTrk();

  float z1, z2, z3;
  z1 = (GDCGeom::Id(GDCGeom::kBC1).GetZPos(1)
    + GDCGeom::Id(GDCGeom::kBC1).GetZPos(3))/2.;
  z2 = (GDCGeom::Id(GDCGeom::kBC2).GetZPos(1)
    + GDCGeom::Id(GDCGeom::kBC2).GetZPos(3))/2.;
  z3 = (GDCGeom::Id(GDCGeom::kBC3).GetZPos(1)
    + GDCGeom::Id(GDCGeom::kBC3).GetZPos(3))/2.;

  if (fCandTraj.empty()) DoCrossings(z1, z2, z3);

  // Fill the ntuple
  if (fNTrack<=10) FillNt();

  // Store the line fit in the event
  float x1, y1;
  float x3, y3;
  int   wire[12];

  std::list<TrackFit_t*>::iterator trackFitItr(fTrackFit.begin());
  std::list<TrackFit_t*>::iterator trackFitItrEnd(fTrackFit.end());

  for ( ; trackFitItr != trackFitItrEnd; ++trackFitItr) {
    x1 = (*trackFitItr)->x1;
    y1 = (*trackFitItr)->y1;
    x3 = (*trackFitItr)->x3;
    y3 = (*trackFitItr)->y3;

    for (int i=0; i<12; i++)
      wire[i] = (*trackFitItr)->wire[i];

    BCLine l(x1, y1, z1, x3, y3, z3, wire, (*trackFitItr)->goldenTrk);

    evt.Reco().Put(l, "./bc/linefit");
  }

  Clear();

  return JobCModule::kPassed;
}

void BCReco::TrkFind

(

)

[private]

Definition at line 724 of file BCReco.cxx.

References c1, fBCAHitList, fCandTrk, fN_TrkFind, fTrkFindHisto, kBC1, kBC2, kBC3, and kNPlane.

Referenced by Reco().

{
//======================================================================
// Track Finding
// - Loops through all pairs of associated hits in BC1 and BC3.
// - Search a road +-n wires wide (fN_TrkFind) around the prediction w2
//   in the list of hits BC2
// - Pairs passing this test are kept in track candidate list (fCandTrk)
//======================================================================
  for (int iplane=0; iplane<GDCGeom::kNPlane; ++iplane) {
    float z1, z2, z3;
    float c1, c2;      // geometric constants

    z1 = GDCGeom::Id(GDCGeom::kBC1).GetZPos(iplane+1);
    z2 = GDCGeom::Id(GDCGeom::kBC2).GetZPos(iplane+1);
    z3 = GDCGeom::Id(GDCGeom::kBC3).GetZPos(iplane+1);
    c1 = (z3-z2)/(z3-z1);  
    c2 = (z2-z1)/(z3-z1);

    std::list<AssocHit_t*>::iterator w1Itr(fBCAHitList[0][iplane].begin());
    std::list<AssocHit_t*>::iterator w1ItrEnd(fBCAHitList[0][iplane].end());

    for ( ; w1Itr != w1ItrEnd; ++w1Itr) {
      AssocHit_t* hit0 = *w1Itr;

      std::list<AssocHit_t*>::iterator w3Itr(fBCAHitList[2][iplane].begin());
      std::list<AssocHit_t*>::iterator w3ItrEnd(fBCAHitList[2][iplane].end());
    
      for ( ; w3Itr != w3ItrEnd; ++w3Itr) {
    AssocHit_t* hit2 = *w3Itr;
    
    float w2;      // predicted wire in BC2
    w2 = c1*hit0->aHitWire + c2*hit2->aHitWire;

    std::list<AssocHit_t*>::iterator 
      w2Itr(fBCAHitList[1][iplane].begin());
    std::list<AssocHit_t*>::iterator 
      w2ItrEnd(fBCAHitList[1][iplane].end());

    float minD = 320.;
    for ( ; w2Itr != w2ItrEnd; ++w2Itr) {
      AssocHit_t* hit1 = *w2Itr;

      if ((w2-fN_TrkFind)<= hit1->aHitWire &&
          (w2+fN_TrkFind)>= hit1->aHitWire) {
        CandPair_t* candPair = new CandPair_t;  // pairs of candidates
        candPair->hit1 = hit0;
        candPair->hit2 = hit1;
        candPair->hit3 = hit2;
        fCandTrk[iplane].push_back(candPair);
      }
      if (fabs(w2 - hit1->aHitWire) < minD) {
        minD = w2 - hit1->aHitWire;
      }
    }
    if (fTrkFindHisto[iplane] && minD!=320.) {
      fTrkFindHisto[iplane]->Fill(minD);
    }
      }
    }
  }

  if (fDebug>1) {
    std::cout << "Track Find" << std::endl;
    for (int iplane=0; iplane<GDCGeom::kNPlane; ++iplane) {
      std::list<CandPair_t*>::iterator candTrkItr(fCandTrk[iplane].begin());
      std::list<CandPair_t*>::iterator candTrkItrEnd(fCandTrk[iplane].end());

      for ( ; candTrkItr != candTrkItrEnd; ++candTrkItr) {
    CandPair_t* thisPair = (*candTrkItr);
    std::cout << "Plane " << iplane << " "
          << thisPair->hit1->aHitWire << " "
          << thisPair->hit2->aHitWire << " "
          << thisPair->hit3->aHitWire
          << std::endl;
      }
    } 
  }
}

void BCReco::TrkFit

(

)

[private]

Definition at line 1114 of file BCReco.cxx.

References fFitFunc, FitLSQ(), FitTMinuit(), and fNTrack.

Referenced by Reco().

{
  fNTrack = 0;

  switch (fFitFunc) {
  case 0: // TMinuit
    FitTMinuit();
    return;

  case 1: // LSQ
    FitLSQ();
    return;

  default:
    std::cerr << "Invalid BCReco fit function " << fFitFunc
          << std::endl;
    exit(-1);
    break;
  }
}

void BCReco::Update

(

const CfgConfig &

c

)

[virtual]

Implements CfgObserver.

Definition at line 1595 of file BCReco.cxx.

References fDebug, fFitFunc, fMaxItr, fMaxNWirePerChamber, fMaxTotalNWire, fN_TrkFind, fN_ViewMatch, fNWireItr, fPrintLevel, fWidth_TrkFit, and fWidthTrks.

{
  c("Debug").       Get(fDebug);
  c("PrintLevel").  Get(fPrintLevel);
  c("N_TrkFind").   Get(fN_TrkFind);
  c("N_ViewMatch"). Get(fN_ViewMatch);
  c("MaxItr").      Get(fMaxItr);
  c("NWireItr").    Get(fNWireItr);
  c("Width_TrkFit").Get(fWidth_TrkFit);
  c("WidthTrks").   Get(fWidthTrks);
  c("FitFunc").     Get(fFitFunc);
  c("MaxNWirePerChamber").Get(fMaxNWirePerChamber);
  c("MaxTotalNWire").Get(fMaxTotalNWire);
}

void BCReco::ViewMatch

(

)

[private]

Definition at line 806 of file BCReco.cxx.

References ct, fCandTraj, fDebug, MatchSU(), MatchTV(), and NULL.

Referenced by Reco().

{
//======================================================================
// View Matching
// - Loop over track candidate list.
// - Predict wire indices in complementary views.
// - Form straight line 3D trajectory.
//======================================================================
  MatchSU();
  MatchTV();

  if (fDebug>1) {
    std::cout << "View Match Trajectory " << std::endl;
    std::list<CandTrajectory_t*>::iterator candTrajItr(fCandTraj.begin());
    std::list<CandTrajectory_t*>::iterator candTrajItrEnd(fCandTraj.end());
    for ( ; candTrajItr != candTrajItrEnd; ++candTrajItr) {
      CandTrajectory_t* ct = *candTrajItr;
      std::cout << "x1 = " << ct->x1 << " "
        << "y1 = " << ct->y1 << " "
        << "x3 = " << ct->x3 << " "
        << "y3 = " << ct->y3 << std::endl;
      for (int i=0; i<GDCGeom::kNBCPlane; i++)
    if (ct->ahit[i]!=NULL)
      std::cout << ct->ahit[i]->aHitWire << " ";
      std::cout << std::endl;
    }
  }
}

---

Member Data Documentation

TH2F* BCReco::fBC3xdxdzHisto [private]

x:dx/dz

Definition at line 153 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH2F* BCReco::fBC3xyHisto [private]

x:y

Definition at line 152 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH2F* BCReco::fBC3ydydzHisto [private]

y:dy/dz

Definition at line 154 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

std::list<AssocHit_t*> BCReco::fBCAHitList[3][4] [private]

Definition at line 129 of file BCReco.h.

Referenced by Clear(), HitAssoc(), MatchSU(), MatchTV(), and TrkFind().

std::vector<Hit*> BCReco::fBCHitList[3][4] [private]

Definition at line 128 of file BCReco.h.

Referenced by Clear(), DoCrossings(), FitLSQ(), FitTMinuit(), HitAssoc(), and Reco().

std::list<CandTrajectory_t*> BCReco::fCandTraj [private]

Definition at line 131 of file BCReco.h.

Referenced by Clear(), FitLSQ(), FitTMinuit(), MatchSU(), MatchTV(), Reco(), and ViewMatch().

std::list<CandPair_t*> BCReco::fCandTrk[4] [private]

Definition at line 130 of file BCReco.h.

Referenced by Clear(), MatchSU(), MatchTV(), and TrkFind().

TH1F* BCReco::fChi2Histo [private]

chi2

Definition at line 155 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

int BCReco::fDebug [private]

Definition at line 161 of file BCReco.h.

Referenced by BCReco(), HitAssoc(), MatchSU(), MatchTV(), Update(), and ViewMatch().

TH1F* BCReco::fdsdzHisto [private]

ds/dz

Definition at line 148 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH1F* BCReco::fdxdzHisto [private]

dx/dz

Definition at line 146 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH1F* BCReco::fdydzHisto [private]

dy/dz

Definition at line 147 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH1F* BCReco::fFinalTrkHisto [private]

Distance(Fit track - hit wire).

Definition at line 157 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

int BCReco::fFitFunc [private]

Definition at line 170 of file BCReco.h.

Referenced by TrkFit(), and Update().

int BCReco::fMaxItr [private]

Definition at line 165 of file BCReco.h.

Referenced by FitTMinuit(), and Update().

int BCReco::fMaxNWirePerChamber [private]

Definition at line 171 of file BCReco.h.

Referenced by Reco(), and Update().

int BCReco::fMaxTotalNWire [private]

Definition at line 172 of file BCReco.h.

Referenced by Reco(), and Update().

TMinuit* BCReco::fMinuit [private]

MINUIT interface.

Definition at line 138 of file BCReco.h.

Referenced by BCReco(), and FitTMinuit().

int BCReco::fN_TrkFind [private]

Definition at line 163 of file BCReco.h.

Referenced by TrkFind(), and Update().

int BCReco::fN_ViewMatch [private]

Definition at line 164 of file BCReco.h.

Referenced by MatchSU(), MatchTV(), and Update().

int BCReco::fNTrack [private]

Definition at line 133 of file BCReco.h.

Referenced by CleanUpTrk(), DoCrossings(), FitLSQ(), FitTMinuit(), Reco(), and TrkFit().

TH1F* BCReco::fnTrackHisto [private]

Number of tracks.

Definition at line 141 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TTree* BCReco::fNtTree [private]

Ntuple tree.

Definition at line 135 of file BCReco.h.

Referenced by FillNt(), MakeHistos(), and ~BCReco().

BCRecoNt* BCReco::fNtuple [private]

Branch for the ntuple tree.

Definition at line 136 of file BCReco.h.

Referenced by FillNt(), MakeHistos(), Reco(), and ~BCReco().

int BCReco::fNWireItr [private]

Definition at line 166 of file BCReco.h.

Referenced by FitLSQ(), and Update().

int BCReco::fPrintLevel [private]

Definition at line 162 of file BCReco.h.

Referenced by BCReco(), and Update().

TH1F* BCReco::fRChi2Histo [private]

Reduced chi2.

Definition at line 156 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH1F* BCReco::ftargetxHisto [private]

x at target

Definition at line 149 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH2F* BCReco::ftargetxyHisto [private]

x:y at target

Definition at line 151 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH1F* BCReco::ftargetyHisto [private]

y at target

Definition at line 150 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

std::list<TrackFit_t*> BCReco::fTrackFit [private]

Definition at line 132 of file BCReco.h.

Referenced by CleanUpTrk(), Clear(), DoCrossings(), FillNt(), FitLSQ(), FitTMinuit(), and Reco().

TH1F* BCReco::fTrkFindHisto[4] [private]

Predicted wire - hit wire in TrkFind.

Definition at line 139 of file BCReco.h.

Referenced by BCReco(), MakeHistos(), TrkFind(), and ~BCReco().

TH1F* BCReco::fViewMatchHisto[3][4] [private]

Predicted wire - hit wire in ViewMatch.

Definition at line 140 of file BCReco.h.

Referenced by BCReco(), MakeHistos(), MatchSU(), MatchTV(), and ~BCReco().

float BCReco::fWidth_TrkFit [private]

Definition at line 168 of file BCReco.h.

Referenced by DoCrossings(), FitLSQ(), FitTMinuit(), and Update().

float BCReco::fWidthTrks [private]

Definition at line 169 of file BCReco.h.

Referenced by CleanUpTrk(), and Update().

TH1F* BCReco::fx1Histo [private]

x1

Definition at line 142 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH1F* BCReco::fx3Histo [private]

x3

Definition at line 144 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH1F* BCReco::fy1Histo [private]

y1

Definition at line 143 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

TH1F* BCReco::fy3Histo [private]

y3

Definition at line 145 of file BCReco.h.

Referenced by MakeHistos(), and ~BCReco().

---

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

• BCReco.h
• BCReco.cxx

---