TBuilderCurl Class Reference

A class to build a Curl track. More...

#include <TBuilderCurl.h>

Inheritance diagram for TBuilderCurl:

Public Member Functions
	TBuilderCurl (const std::string &name)
	Constructor.
virtual	~TBuilderCurl ()
	Destructor.
TTrack *	buildStereo (TTrack &track, const AList< TMLink > &) const
	appends stereo hits to a track.
TTrack *	buildStereoMC (TTrack &track, const AList< TMLink > &) const
void	setParam (const TCurlFinderParameter &)
TTrack *	buildStereo (TTrack &track, const AList< TMLink > &, const AList< TMLink > &) const
bool	buildStereo (TTrack &track, double &dZ, double &tanL) const
Public Member Functions inherited from TBuilder0
	TBuilder0 (const std::string &name)
	Constructor.
	TBuilder0 (const std::string &name, float salvageLevel)
	Constructor with salvage level.
	TBuilder0 (const std::string &name, float stereoZ3, float stereoZ4, float stereoChisq3, float stereoChisq4, float stereoMaxSigma, unsigned fittingCorrections, float salvageLevel)
	Constructor with parameters.
virtual	~TBuilder0 ()
	Destructor.
const std::string &	name (void) const
	returns name.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
const TMSelector &	trackSelector (void) const
	returns a track selector.
TTrack *	buildRphi (const AList< TMLink > &) const
	builds a r/phi track from TMLinks or from Segments.
TTrack *	buildStereo0 (TTrack &track, const AList< TMLink > &) const
	appends stereo hits to a track. (old version)
virtual TTrack *	buildStereo (TTrack &track, const AList< TMLink > &) const
	appends stereo hits to a track.
void	appendClusters (TTrack &track, const AList< TMLink > &) const
	appends TMLinks in a list.
void	salvage (TTrack &track, AList< TMLink > &list) const
	salvages links in a list. Used links will be removed from a list.
virtual int	fit (TTrackBase &) const
	fits a track using a private fitter.
virtual const TMSelector &	trackSelector (const TMSelector &)
	sets a track selector.

Additional Inherited Members
Protected Attributes inherited from TBuilder0
TMSelector	_circleSelector
TMSelector	_trackSelector
TMSelector	_lineSelector
float	_stereoZ3
float	_stereoZ4
float	_stereoChisq3
float	_stereoChisq4
float	_stereoMaxSigma

Detailed Description

A class to build a Curl track.

Definition at line 54 of file TBuilderCurl.h.

Constructor & Destructor Documentation

TBuilderCurl()

TBuilderCurl::TBuilderCurl

(

const std::string &

name

)

Constructor.

Definition at line 59 of file TBuilderCurl.cxx.

: TBuilder0(name),
  _fitter("TBuilderCurl Fitter"){
 
  //jialk
  StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF); 
  if(scmgn!=StatusCode::SUCCESS) { 
    std::cout<< __FILE__<<"Unable to open Magnetic field service"<<std::endl;
  }
  
#if 0
//  if(m_param.ON_CORRECTION){
//    _fitter.sag(true);
//    _fitter.propagation(true);
//    _fitter.tof(true);
//  }
  if (m_param.ON_CORRECTION & 1) _fitter.sag(true);
  if (m_param.ON_CORRECTION & 2) _fitter.propagation(true);
  if (m_param.ON_CORRECTION & 4) _fitter.tof(true);
#endif
}

~TBuilderCurl()

TBuilderCurl::~TBuilderCurl ( )

virtual

Destructor.

Definition at line 81 of file TBuilderCurl.cxx.

                            {
  //  if(m_param.SVD_RECONSTRUCTION){
  //    delete m_svdFinder;
  //    delete m_svdAssociator;
  //  }
}

Member Function Documentation

buildStereo() [1/3]

TTrack * TBuilderCurl::buildStereo ( TTrack &  track, const AList< TMLink > &  stereoList  ) const

virtual

appends stereo hits to a track.

Reimplemented from TBuilder0.

Definition at line 182 of file TBuilderCurl.cxx.

                                                      {
  return NULL;
}

buildStereo() [2/3]

TTrack * TBuilderCurl::buildStereo	(	TTrack &	track,
		const AList< TMLink > &	stereoList,
		const AList< TMLink > &	allStereoList
	)		const

Definition at line 188 of file TBuilderCurl.cxx.

                                                         {
#if (DEBUG_CURL_DUMP+DEBUG_CURL_GNUPLOT+DEBUG_CURL_MC)
  Belle_event_Manager &evtMgr = Belle_event_Manager::get_manager();
  debugMcFlag = 1;
  if(evtMgr.count() != 0 &&
     evtMgr[0].ExpMC() != 2)debugMcFlag = 0;// not MC
#endif
 
  AList<TMLink> list = stereoList;
  unsigned nList = list.length();
 
  //...gets stereo wires from track.links
  for(unsigned i = 0, size = track.links().length(); i < size; ++i){
    unsigned superID = (track.links())[i]->wire()->superLayerId();
    if(superID == 0 || superID == 1 || superID == 5 || 
       superID == 6 || superID == 7 || superID == 8  ){
      int ok = 1;
      for (unsigned j = 0; j < nList; ++j) {
    TMLink * l = list[j];
    if(l->hit()->wire()->id() == (track.links())[i]->hit()->wire()->id()){
      ok = 0; break;
    }
      }
      if(ok == 1)list.append(((track.links())[i]));
    }
    // set LR 2
    (track.links())[i]->leftRight(2);  // returns left-right. 0:left, 1:right, 2:wire
  }
  for(unsigned i = 0, size = list.length(); i < size; ++i){    
    track._links.remove(list[i]);
    // set LR 2
    list[i]->leftRight(2);
  }
 
  if(list.length() < 2 || 
     list.length()+track.nLinks() < 5)return NULL;  
#if DEBUG_CURL_DUMP
  unsigned debug_stereo_counter1 = 0;
  for(unsigned i=0;i<track.links().length();++i){
    unsigned superID = (track.links())[i]->hit()->wire()->superLayerId();
    if(superID == 0 || superID == 1 ||
       superID == 5 || superID == 6 ||
       superID == 7 || superID == 8 )++debug_stereo_counter1;
  }
  std::cout << "(TBuilderCurl)Fitted Track:";
  std::cout << ", A# = " << track.links().length()-debug_stereo_counter1;
  std::cout << ", S# = " << debug_stereo_counter1 << "(==0)";
  std::cout << ", A+S# = " << track.links().length();
  std::cout << ", Cand Stereo Wires # = " << list.length() << std::endl;
  double debugChg = -1.;
  if(track.charge() > 0.)debugChg = 1.;
  if(debugChg > 0.)std::cout << "(TBuilderCurl)  Positive Track" << std::endl;
  else std::cout << "(TBuilderCurl)  Negative Track" << std::endl;
#endif
 
  //...calculates a circle.
  double xc2D;
  double yc2D;
  double r2D;
  AList<TMLink> tmpAxialList = track.links();
  bool err2D = fitWDD(xc2D,yc2D,r2D,tmpAxialList);
 
  //...using SVD   //....Liuqg
  //double dZSVD, tanLSVD;
  //bool initWithSVD = buildStereo(track, dZSVD, tanLSVD);
 
  //...sets arc and z pairs of each stereo hit wire.
  setArcZ(track,list);
  AList<TMLink> removeList;
  for(unsigned i = 0, size = list.length(); i < size; ++i){
    if(list[i]->arcZ(0).x() == -999. && list[i]->arcZ(1).x() == -999.)removeList.append(list[i]);
  }
  list.remove(removeList);
  if(list.length() < 2 || 
     list.length()+track.nLinks() < 5){ // stereo >=2 && axial >= 3
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)  Fail...few wires which can be set Arc-Z # = " 
     << list.length() << std::endl;
#endif    
    return NULL;
  }
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)  Cand Stereo Wires which can be set Arc-Z # = " 
     << list.length() << std::endl;
    plotArcZ(list,0.,0.,0);
#endif 
  
#if 0
  //...separates
  AList<TMLink> layer[18];
  for(unsigned i = 0, size = list.length(); i < size; ++i){
    layer[list[i]->wire()->layer()->axialStereoLayerId()].append(*list[i]);
  }
 
#if DEBUG_CURL_DUMP
  for(int i=0;i<18;++i){
    if(layer[i].length())
      std::cout << "(TBuilderCurl)  Cand Stereo Wires which can be set Arc-Z # = " 
       << layer[i].length()
       << " on " << i << " Layer." << std::endl;
  }
#endif
#endif
 
  // makes a line.
  AList<TMLink> goodL;
  AList<HepPoint3D> goodP;
  double minChi2 = 9999.;
  double goodA   = 9999.;
  double goodB   = 9999.;
  makeLine(track, list, allStereoList, goodL, 
       minChi2, goodA, goodB, goodP);
  HepAListDeleteAll(goodP);
 
#if DEBUG_CURL_DUMP
  std::cout << "(TBuilderCurl)    a = " << goodA << ", b = " << goodB 
        << " (after makeLine-function)" << std::endl;
#endif
 
  // refits
  static TRobustLineFitter rfitter("Can you work well?");
  TMLine newLine0(goodL);
  if(rfitter.fit(newLine0) != 0){
#if DEBUG_CURL_DUMP
      std::cout << "(TBuilderCurl)    Fail...linefitting...fail." << std::endl;
#endif
    return NULL;
  }
 
#if DEBUG_CURL_DUMP
  std::cout << "(TBuilderCurl)    a = " << newLine0.a() << ", b = " << newLine0.b() 
        << " (after robustline-fit)" << std::endl;
#endif
 
  //...appends at last chance
  unsigned nGoodL = goodL.length();
  list.remove(goodL);
  AList<TMLink> goodL2;
  if(m_param.CURL_VERSION == 0){ // default
    if(fabs(newLine0.b()) < 10.){
      appendPoints(list, goodL2,
           newLine0.a(), newLine0.b(), track, m_param.Z_DIFF_FOR_LAST_ATTEND);
    }else{ // in case of bad result of robustLineFitter. (2001/04/04)
      appendPoints(list, goodL2,
           goodA, goodB, track, m_param.Z_DIFF_FOR_LAST_ATTEND);    
    }
  }else{
    // same with b20010409_2122 at least
    appendPoints(list, goodL2,
         newLine0.a(), newLine0.b(), track, m_param.Z_DIFF_FOR_LAST_ATTEND);
  }
  goodL.append(goodL2);
  TMLine newLine(goodL);
 
  //...refits
  if(rfitter.fit(newLine) != 0){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)    Fail...appending and re-fitting...fail." << std::endl;
#endif
    return NULL;
  }
 
#if DEBUG_CURL_DUMP
  std::cout << "(TBuilderCurl)    a = " << newLine.a() << ", b = " << newLine.b() 
        << " (after last-append + re-robustline-fit)" << std::endl;
#endif
  //...makes 3D tracks
  const AList<TMLink> &good = newLine.links();
  track.TTrackBase::append(good);
  if(!check(track)){
#if DEBUG_CURL_DUMP
      std::cout << "(TBuilderCurl)    Fail...checking wire numbers...fail." << std::endl;
#endif
    return NULL;
  }
 
  //...sets initial values
  Vector a = track.helix().a();
  if(err2D){
    double tmpA[3];
    double tmpQ = track.charge() > 0. ? 1. : -1.;
    tmpA[1] = fmod(atan2(tmpQ * yc2D,
             tmpQ * xc2D)
           + 4. * M_PI,
           2. * M_PI);
   // tmpA[2] = tmpQ*Helix::ConstantAlpha/r2D;
   // tmpA[2] = tmpQ*333.564095/r2D;
    tmpA[2] = tmpQ*(333.564095/(-1000*(m_pmgnIMF->getReferField())))/r2D;
//cout<<"TBuilderCurl::tmpA[2]: "<<(333.564095/(-1000*(m_pmgnIMF->getReferField())))<<endl;
    tmpA[0] = xc2D/cos(tmpA[1]) - tmpQ*r2D;
    //std::cout << yc2D/sin(tmpA[1])-tmpQ*r2D << std::endl;
    //std::cout << a[0] << " -0- " << tmpA[0] << std::endl;
    //std::cout << a[1] << " -1- " << tmpA[1] << std::endl;
    //std::cout << a[2] << " -2- " << tmpA[2] << std::endl;
    a[0] = tmpA[0];
    a[1] = tmpA[1];
    a[2] = tmpA[2];
  }
  a[3] = newLine.b();
  a[4] = newLine.a();
//Liuqg
/*  if(initWithSVD){ // use SVD initial values if possible.
    a[3] = dZSVD;
    a[4] = tanLSVD;
  }
*/
  if(m_param.CURL_VERSION == 0){
    if(fabs(a[3]) > 10. && fabs(goodB) < 10.){ // 50cm --> 10cm @ 2001/04/04
      // use initial values when results of RobustFit is bad.
      a[3] = goodB;
      a[4] = goodA;
    }
  }else{
    if(fabs(a[3]) > 50. && fabs(goodB) < 50.){
      a[3] = goodB;
      a[4] = goodA;
      track.TTrackBase::remove(goodL2);
    }
  }
 
  track._helix->a(a);
 
#if DEBUG_CURL_DUMP
  std::cout << "(TBuilderCurl)    Created Line: y = " << newLine.a() 
       << " * x + " << newLine.b() 
       << ", size = " << good.length() << std::endl;
  plotArcZ(const_cast< AList<TMLink>& >(good), newLine.a(), newLine.b(), 0);
#endif
 
  /* std::cout << "1st : "
       << track.helix().a()[0] << " " << track.helix().a()[1] << " "
       << track.helix().a()[2] << " " << track.helix().a()[3] << " "
       << track.helix().a()[4] << std::endl; */
 
#if 1
//  if(m_param.ON_CORRECTION){
//    _fitter.sag();
//    _fitter.propagation();
//    _fitter.tof();
//  }
  if (m_param.ON_CORRECTION & 1) _fitter.sag(true);
  if (m_param.ON_CORRECTION & 2) _fitter.propagation(true);
  if (m_param.ON_CORRECTION & 4) _fitter.tof(true);
#endif
 
  //...fits
  AList<TMLink> bad;
  int err = _fitter.fit(track);
  if (err < 0){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail fitting(0)...error code = " << err << std::endl;
#endif
    return NULL;
  } else if ( fabs(track.helix().a()[3]) > fabs(a[3]) &&
          (fabs(track.helix().a()[3]) > 50.  || // 50cm
           fabs(track.helix().a()[2]) > 100. || // 0.01GeV
           fabs(track.helix().a()[2]) < 0.1) ){ // 10GeV
    // in strange case, set "correction" of wires OFF
    // and then fit with the initial values.
    if(m_param.ON_CORRECTION){
      if (m_param.ON_CORRECTION & 1) _fitter.sag(true);
      if (m_param.ON_CORRECTION & 2) _fitter.propagation(true);
      if (m_param.ON_CORRECTION & 4) _fitter.tof(true);
      //_fitter.sag();
      //_fitter.propagation();
      //_fitter.tof();
      track._helix->a(a);
      //std::cout << "ON --> OFF" << std::endl;
    }else{
      return NULL;
    }
  }
 
  /* std::cout << "2nd : "
       << track.helix().a()[0] << " " << track.helix().a()[1] << " "
       << track.helix().a()[2] << " " << track.helix().a()[3] << " "
       << track.helix().a()[4] << std::endl; */
 
  track.refine(bad, m_param.SELECTOR_MAX_SIGMA * 100.);
  if (!check(track)){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail checking(1)..." << std::endl;
#endif
    return NULL;
  }
  err = _fitter.fit(track);
  if (err < 0){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail fitting(1)...error code = " << err << std::endl;
#endif
    return NULL;
  }
  track.refine(bad, m_param.SELECTOR_MAX_SIGMA * 10.);
  if (!check(track)){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail checking(2)..." << std::endl;
#endif
    return NULL;
  }
  err = _fitter.fit(track);
  if (err < 0){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail fitting(2)...error code = " << err << std::endl;
#endif
    return NULL;
  }
  track.refine(bad, m_param.SELECTOR_MAX_SIGMA);
  if (!check(track)){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail checking(3)..." << std::endl;
#endif
    return NULL;
  }
  err = _fitter.fit(track);
  if (err < 0){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail fitting(3)...error code = " << err << std::endl;
#endif
    return NULL;
  }
 
#if 0
  for(int i=0;i<track.links().length();++i){
    if(track.links()[i]->pull() > 36.){
      std::cout << track.links()[i]->wire()->id() 
       << " :+: " << track.links()[i]->pull() << std::endl;
    }
    if(track.links()[i]->hit()->state() & WireHitInvalidForFit)
      std::cout << "Not Valid For Fit!" << std::endl;
    //if(!(track.links()[i]->hit()->state() & WireHitFittingValid))
    //std::cout << "No-Valid For Fit!" << std::endl;
  }
#endif
  track.refine(bad, m_param.SELECTOR_MAX_SIGMA);
  if (!check(track)){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail checking(4)..." << std::endl;
#endif
    return NULL;
  }
  err = _fitter.fit(track);
  if (err < 0){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)     Fail fitting(4)...error code = " << err << std::endl;
#endif
    return NULL;
  }
 
#if 0    
  for(int i=0;i<track.links().length();++i){
    if(track.links()[i]->pull() > 36.){
      std::cout << track.links()[i]->wire()->id() 
       << " :*: " << track.links()[i]->pull() << std::endl;
    }
    if(track.links()[i]->hit()->state() & WireHitInvalidForFit)
      std::cout << "Not Valid For Fit!" << std::endl;
    //if(!(track.links()[i]->hit()->state() & WireHitFittingValid))
    //std::cout << "No-Valid For Fit!" << std::endl;
  }
#endif
 
  /* std::cout << "3rd : "
       << track.helix().a()[0] << " " << track.helix().a()[1] << " "
       << track.helix().a()[2] << " " << track.helix().a()[3] << " "
       << track.helix().a()[4] << std::endl; */
 
  //...tests
  if (track.nLinks() < 5){ // axial + stereo >= 5
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)      Success fitting, but pre-selection...fail." << std::endl;
#endif
    return NULL;
  }
  if(fabs(track.impact()) > m_param.SELECTOR_MAX_IMPACT ||
     track.pt() < 0.005 || // Pt >= 5MeV
     fabs(track.pz()) > m_param.SELECTOR_STRANGE_PZ){
#if DEBUG_CURL_DUMP
    std::cout << "(TBuilderCurl)      Success fitting, but selection...fail." << std::endl;
    std::cout << "(TBuilderCurl)                           impact = " << track.impact() << std::endl;
    std::cout << "(TBuilderCurl)                           pt     = " << track.pt() << std::endl;
    std::cout << "(TBuilderCurl)                           pz     = " << track.pz() << std::endl;
    std::cout << "(TBuilderCurl)                           dz     = " << track.helix().a()[3] << std::endl;
#endif  
    return NULL;
  }
 
  //...replaces init helix if dz is bad.
  if(fabs(track.helix().a()[3]) > m_param.SELECTOR_REPLACE_DZ &&
     fabs(a[3]) < fabs(track.helix().a()[3])){
    track._helix->a(a);
  }
#if DEBUG_CURL_DUMP
  std::cout << "(TBuilderCurl)      Success Build Stereo!!!" << std::endl;
#endif
  return & track;
}

buildStereo() [3/3]

bool TBuilderCurl::buildStereo	(	TTrack &	track,
		double &	dZ,
		double &	tanL
	)		const

Definition at line 134 of file TBuilderCurl.cxx.

                                  {
/*  
 if(!(m_param.SVD_RECONSTRUCTION))return false;
 
#if 0
  double q = 1.;
  if(track.helix().kappa()<0.)q = -1.;
  double h[5], chisq;
  if(m_svdFinder->recTrk(track.helix().center().x(),
             track.helix().center().y(),
             fabs(track.helix().radius()),
             q,
             0.2,
             h[0], h[1], h[2], h[3], h[4], chisq)){
#if 0
    std::cout << "SVD = "
     << h[0] << " "
     << h[1] << " "
     << h[2] << " "
     << h[3] << " "
     << h[4];
    std::cout << ", chi2 = " << chisq << " pt=" << 1./h[2] << std::endl;
#endif
    dZ   = h[3];
    tanL = h[4];
    return true;
  }else{
    //std::cout << "SVD Fail....." << std::endl;
    return false;
  }
#else
  Helix th(track.helix());
  th.pivot(HepPoint3D(0.,0.,0.));
  AList<TSvdHit> cand;
  if(m_svdAssociator->recTrk(th,dZ,tanL,0.5,-1.0,cand,0.5)){
    //std::cout << "SVD " << dZ << " " << tanL << std::endl;
    return true;
  }else{
    //std::cout << "SVD..." << std::endl;
    return false;
  }
#endif
 */
}

buildStereoMC()

TTrack * TBuilderCurl::buildStereoMC	(	TTrack &	track,
		const AList< TMLink > &	stereoList
	)		const

Definition at line 1050 of file TBuilderCurl.cxx.

                                                                                  {
#if DEBUG_CURL_MC
  AList<TMLink> list = stereoList;
 
  HepPoint3D center = track.helix().center();
  double r  = fabs(track.helix().curv());
  for(unsigned i = 0, size = list.length();
      i < size; ++i){
    HepPoint3D point((list[i]->hit()->mc()->datcdc()->m_xin+
           list[i]->hit()->mc()->datcdc()->m_xout)*0.5,
          (list[i]->hit()->mc()->datcdc()->m_yin+
           list[i]->hit()->mc()->datcdc()->m_yout)*0.5,
          0.);
    double cosdPhi = - center.dot((point - center).unit()) / center.mag();
    double dPhi;
    if(fabs(cosdPhi) < 1.0){
      dPhi = acos(cosdPhi);
    }else if(cosdPhi >= 1.0){
      dPhi = 0.;
    }else{
      dPhi = M_PI;
    }
    list[i]->position(HepPoint3D(r*dPhi,
                  (list[i]->hit()->mc()->datcdc()->m_zin+
                   list[i]->hit()->mc()->datcdc()->m_zout)*0.5,
                  0.));
    /* std::cout << list[i]->wire()->id() << ": " 
     << point.x() << " "
     << point.y() << " "
     << (list[i]->hit()->mc()->datcdc()->m_zin+
     list[i]->hit()->mc()->datcdc()->m_zout)*0.5 << std::endl; */
  }
  //std::cout << "A# = " << track.links().length() << ", S# = " << list.length() << std::endl;
  double xc2D;
  double yc2D;
  double r2D;
  bool err2D = fitWDD(xc2D,yc2D,r2D,track.links()); // axial only
 
  TLine0 newLine(list);
  if(newLine.fit() != 0) return NULL;
  const AList<TMLink> &good = newLine.links();
  track.append(good);
  Vector a = track.helix().a();
  a[3] = newLine.b();
  a[4] = newLine.a();
  if(err2D){
    double tmpA[3];
    double tmpQ = track.charge() > 0. ? 1. : -1.;
    tmpA[1] = fmod(atan2(tmpQ * yc2D,
             tmpQ * xc2D)
           + 4. * M_PI,
           2. * M_PI);
   // tmpA[2] = tmpQ*Helix::ConstantAlpha/r2D;
   // tmpA[2] = tmpQ*333.564095/r2D;
    tmpA[2] =tmpQ*(333.564095/(-1000*(m_pmgnIMF->getReferField())))/r2D;
    tmpA[0] = xc2D/cos(tmpA[1]) - tmpQ*r2D;
    a[0] = tmpA[0];
    a[1] = tmpA[1];
    a[2] = tmpA[2];
  }
  track._helix->a(a);
#if 0
  std::cout << track.helix().a()[0] << " " << track.helix().a()[1] << " "
       << track.helix().a()[2] << " " << track.helix().a()[3] << " "
       << track.helix().a()[4] << std::endl;
#endif
  //...fits
  AList<TMLink> bad;
  int err = _fitter.fit(track);
  if (err < m_param.ERROR_FOR_HELIX_FIT)return NULL;
  track.refine(bad, m_param.SELECTOR_MAX_SIGMA * 100.);
  err = _fitter.fit(track);
  if (err < m_param.ERROR_FOR_HELIX_FIT)return NULL;
  track.refine(bad, m_param.SELECTOR_MAX_SIGMA * 10.);
  err = _fitter.fit(track);
  if (err < m_param.ERROR_FOR_HELIX_FIT)return NULL;
  track.refine(bad, m_param.SELECTOR_MAX_SIGMA);
  err = _fitter.fit(track);
  if (err < m_param.ERROR_FOR_HELIX_FIT)return NULL;
#if 0
  std::cout << track.helix().a()[0] << " " << track.helix().a()[1] << " "
       << track.helix().a()[2] << " " << track.helix().a()[3] << " "
       << track.helix().a()[4] << std::endl;
#endif
  //track._helix->a(a); // re-write
 
  //...checks
  if (track.nLinks() < m_param.SELECTOR_NLINKS)return NULL;
  if (fabs(track.impact()) > m_param.SELECTOR_MAX_IMPACT ||
      track.pt() < m_param.SELECTOR_MIN_PT)return NULL;
  return & track;
#else
  return NULL;
#endif
}

setParam()

void TBuilderCurl::setParam

(

const TCurlFinderParameter &

p

)

Definition at line 89 of file TBuilderCurl.cxx.

                                                    {
  m_param.Z_CUT = p.Z_CUT;
  m_param.Z_DIFF_FOR_LAST_ATTEND = p.Z_DIFF_FOR_LAST_ATTEND;
  m_param.SELECTOR_MAX_SIGMA = p.SELECTOR_MAX_SIGMA;
  m_param.SELECTOR_MAX_IMPACT = p.SELECTOR_MAX_IMPACT;
  m_param.SELECTOR_STRANGE_PZ = p.SELECTOR_STRANGE_PZ;
  m_param.SELECTOR_REPLACE_DZ = p.SELECTOR_REPLACE_DZ;
  m_param.SVD_RECONSTRUCTION = p.SVD_RECONSTRUCTION;
  m_param.MIN_SVD_ELECTRONS = p.MIN_SVD_ELECTRONS;
  m_param.ON_CORRECTION = p.ON_CORRECTION;
  m_param.CURL_VERSION = p.CURL_VERSION;
#if 1
//  if(m_param.ON_CORRECTION){
//    _fitter.sag(true);
//    _fitter.propagation(true);
//    _fitter.tof(true);
//  }
  if (m_param.ON_CORRECTION & 1) _fitter.sag(true);
  if (m_param.ON_CORRECTION & 2) _fitter.propagation(true);
  if (m_param.ON_CORRECTION & 4) _fitter.tof(true);
#endif
  //  if(m_param.SVD_RECONSTRUCTION){
  ////m_svdFinder = new TSvdFinder(-1.*(m_param.MIN_SVD_ELECTRONS),
  ////           m_param.MIN_SVD_ELECTRONS);
  //m_svdAssociator = new TSvdAssociator(-1.*(m_param.MIN_SVD_ELECTRONS),
  //                     m_param.MIN_SVD_ELECTRONS);
  //}
}

Referenced by TCurlFinder::TCurlFinder().

---

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

• TBuilderCurl.h
• TBuilderCurl.cxx