TTrackMC Class Reference

A class to have MC information of TTrack. More...

#include <TTrackMC.h>

List of all members.

Public Member Functions
virtual	~TTrackMC ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
int	hepId (void) const
	returns HEP ID.
const TTrackHEP *const	hep (void) const
	returns a pointer to TTrackHEP.
bool	charge (void) const
	returns charge matching.
double	wireFraction (void) const
	returns wire fraction(F1).
double	wireFractionHEP (void) const
	returns wire fraction(F2).
double	ptFraction (void) const
	returns pt fraction.
double	pzFraction (void) const
	returns pz fraction.
const Hep3Vector &	residual (void) const
	returns residual of momentum.
double	ptResidual (void) const
	returns residual of pt.
double	pzResidual (void) const
	returns residual of pz.
double	ptPull (void) const
	returns pull of pt.
double	pzPull (void) const
	returns pull of pz.
unsigned	state (void) const
	returns state.
unsigned	quality (void) const
	returns quality.
std::string	qualityString (void) const
	returns quality.
void	update (void)
	updates information.
Private Member Functions
	TTrackMC (const TTrack &)
	Constructor.
void	compare (void)
	checks matching of charge.
void	classify (void)
	classifies itself.
Private Attributes
unsigned	_state
unsigned	_quality
const TTrack &	_t
TTrackHEP *	_hep
int	_hepID
double	_wireFraction
double	_wireFractionHEP
bool	_charge
double	_ptFraction
double	_pzFraction
Hep3Vector	_residual
double	_cosOpen
double	_ptResidual
double	_pzResidual
double	_ptPull
double	_pzPull
Friends
class	TrkReco

---

Detailed Description

A class to have MC information of TTrack.

Definition at line 54 of file TTrackMC.h.

---

Constructor & Destructor Documentation

TTrackMC::TTrackMC

(

const TTrack &

)

[private]

Constructor.

Definition at line 46 of file TTrackMC.cxx.

: _t(t),
  _hep(0),
  _hepID(-1),
  _wireFraction(-999.),
  _wireFractionHEP(-999.),
  _charge(false),
  _ptFraction(-999.),
  _pzFraction(-999.),
  _ptResidual(-999.),
  _pzResidual(-999.),
  _ptPull(-999.),
  _pzPull(-999.),
  _state(0),
  _quality(0) {
}

TTrackMC::~TTrackMC

(

)

[virtual]

Destructor.

Definition at line 63 of file TTrackMC.cxx.

                    {
}

---

Member Function Documentation

bool TTrackMC::charge

(

void

)

const [inline]

returns charge matching.

Definition at line 181 of file TTrackMC.h.

References _charge.

Referenced by TTrackManager::saveMCTables().

                           {
    return _charge;
}

void TTrackMC::classify

(

void

)

[private]

classifies itself.

Definition at line 251 of file TTrackMC.cxx.

References _charge, _cosOpen, _hep, _ptFraction, _ptResidual, _pzFraction, _pzResidual, _quality, _state, _wireFraction, TTrackBad, TTrackCharge, TTrackClassified, TTrackGarbage, TTrackGhost, TTrackHep, TTrackMatchingLoose, TTrackMatchingTight, TTrackPt, TTrackPz, and TTrackWire.

                       {
    _state |= TTrackClassified;
    _quality = TTrackGarbage;

    //...HEP matching...
    if (! _hep) return;
    _quality |= TTrackHep;
    if (fabs(_ptFraction) < .1) _quality |= TTrackPt;
    if (fabs(_pzFraction) < .1) _quality |= TTrackPz;
    float momResidual = sqrt(_ptResidual * _ptResidual +
                              _pzResidual * _pzResidual);

    if ((momResidual < 0.100) && (_cosOpen > 0.99)) 
        _quality |= TTrackMatchingLoose;
    if ((momResidual < 0.020) && (_cosOpen > 0.998))
        _quality |= TTrackMatchingTight;

    //...Wire fraction...
    if (_wireFraction < 0.8) return;
    _quality |= TTrackWire;
    _quality |= TTrackCharge;

    //...Charge matching...
    if (! _charge) return;
    _quality |= TTrackBad;

    //...Momentum matching...
    if (_quality & TTrackMatchingLoose)
        _quality |= TTrackGhost;

    //...TTrackGood is set by TrkReco after checking uniqueness...
}

void TTrackMC::compare

(

void

)

[private]

checks matching of charge.

Definition at line 159 of file TTrackMC.cxx.

References _charge, _cosOpen, _hep, _ptFraction, _ptPull, _ptResidual, _pzFraction, _pzPull, _pzResidual, _residual, _t, abs, TTrack::charge(), Helix::Ea(), TTrack::finalHits(), TTrack::helix(), Helix, TMLink::hit(), genRecEmupikp::i, InnerMost(), Helix::kappa(), TMDCWireHit::mc(), Helix::momentum(), TMDCWireHitMC::momentum(), TTrack::p(), Helix::pivot(), TTrack::pt(), TTrackHEP::pType(), TTrack::pz(), and Helix::tanl().

                      {
    if (! _hep) return;

    //...Get charge of HEP particle(This part should be done by LUCHARGE)...
    int id = _hep->pType();
    int aId = abs(id);
    if (aId == 11 || aId == 13 || aId == 15) id *= -1;

    //...Compare charge...
    if ((int) _t.charge() * id > 0) _charge = true;

    //...Get hep mom. at the inner-most hit...
//  AList<TMLink> list = _t.cores();
    AList<TMLink> list = _t.finalHits();
    unsigned n = list.length();
    bool found = false;
    Vector3 pHep;
    Vector3 vHep;
    for (unsigned i = 0; i < n; i++) {
        TMLink * inner = InnerMost(list);
        if (inner->hit()->mc()->hep() == _hep) {
            pHep = inner->hit()->mc()->momentum();
            vHep = inner->hit()->mc()->entrance();
            found = true;
            break;
        }
        list.remove(inner);
    }
    Helix hHep = Helix(vHep, pHep, copysign(1., id));
    hHep.pivot(_t.helix().pivot());
    pHep = hHep.momentum();

    //...For debug...
    if (! found) {
        std::cout << "TTrackMC::compare !!! something wrong with mc hits"
                  << std::endl;

        //...For debug...
//      list = _t.cores();
//      for (unsigned i = 0; i < list.length(); i++) {
//          TMLink * inner = innerMost(list);
//          std::cout << i << " " << inner << ":" << inner->hit()->mc()->hep();
//          std::cout << " " << _hep << std::endl;
//          if (inner->hit()->mc()->hep() == _hep) {
//              pHep = inner->hit()->mc()->momentum();
//              break;
//          }
//          list.remove(inner);
//      }
//      TMLink * t = 0;
//      t->hit();
        //...For debug end...

        return;
    }

    //...Fill caches...
    _residual = _t.p() - pHep;
    _cosOpen = pHep.unit().dot(_t.p().unit());

    //...Compare pt...
    double pt = _t.pt();
    double ptHep = sqrt(pHep.x() * pHep.x() + pHep.y() * pHep.y());
    _ptResidual = pt - ptHep;
    const Helix & h = _t.helix();
    Vector dPt(5, 0);
    dPt[2] = - 1. / (h.kappa() * h.kappa());
    double ptError2 = h.Ea().similarity(dPt);
    if(ptError2<0.0) {
      std::cout << h.kappa() << " " << h.Ea() << " dPt=" << dPt << std::endl;
    }
    double ptError = (ptError2 > 0.) ? sqrt(ptError2) : (DBL_MIN);
    _ptPull = (ptError2 > 0.) ? (_ptResidual) / ptError : (FLT_MAX);
    _ptFraction = (fabs(ptHep)>(FLT_MIN)) ? _ptResidual / ptHep : 0.0;

    //...Compare pz...
    double pz = _t.pz();
    double pzHep = pHep.z();
    _pzResidual = pz - pzHep;
    Vector dPz(5, 0);
    dPz[2] = - h.tanl() / (h.kappa() * h.kappa());
    dPz[4] = 1. / h.kappa();
    double pzError2 = h.Ea().similarity(dPz);
    if(pzError2<0.0) {
      std::cout << h.kappa() << " " << h.Ea() << " dPz=" << dPz << std::endl;
    }
    double pzError = (pzError2 > 0.) ? sqrt(pzError2) : (DBL_MIN);
    _pzPull = (pzError2 > 0.) ? (_pzResidual) / pzError : (FLT_MAX);
    _pzFraction = (abs(pzHep)>FLT_MIN) ? (_pzResidual / pzHep) : (FLT_MAX);
}

void TTrackMC::dump	(	const std::string &	message = std::string(""),
		const std::string &	prefix = std::string("")
	)			const

dumps debug information.

Definition at line 141 of file TTrackMC.cxx.

References _hepID, _ptFraction, _pzFraction, _quality, _state, _t, _wireFraction, _wireFractionHEP, bitDisplay(), TTrack::name(), TTrackBad, TTrackCharge, TTrackGarbage, TTrackGhost, and TTrackGood.

                                                                 {
    std::cout << msg;
    std::cout << _t.name() << ":";
    std::cout << "state=" << _state << ":";
    if (_quality & TTrackGood) std::cout << "good   :";
    else if (_quality & TTrackGhost) std::cout << "ghost  :";
    else if (_quality & TTrackBad) std::cout << "bad    :";
    else if (_quality & TTrackCharge) std::cout << "bad    :";
    else if (_quality & TTrackGarbage) std::cout << "garbage:";
    else std::cout << "classification error:";
    bitDisplay(_quality, 23, 15); std::cout << ":";
    std::cout << _hepID << ":";
    std::cout << _wireFraction << "," << _wireFractionHEP << ":";
    std::cout << _ptFraction << "," << _pzFraction;
    std::cout << std::endl;
}

const TTrackHEP *const TTrackMC::hep

(

void

)

const [inline]

returns a pointer to TTrackHEP.

Definition at line 175 of file TTrackMC.h.

References _hep.

Referenced by TTrackManager::saveMCTables().

                        {
    return _hep;
}

int TTrackMC::hepId

(

void

)

const [inline]

returns HEP ID.

Definition at line 169 of file TTrackMC.h.

References _hepID.

                          {
    return _hepID;
}

double TTrackMC::ptFraction

(

void

)

const [inline]

returns pt fraction.

Definition at line 187 of file TTrackMC.h.

References _ptFraction.

Referenced by TTrackManager::saveMCTables().

                               {
    return _ptFraction;
}

double TTrackMC::ptPull

(

void

)

const [inline]

returns pull of pt.

Definition at line 241 of file TTrackMC.h.

References _ptPull.

                           {
    return _ptPull;
}

double TTrackMC::ptResidual

(

void

)

const [inline]

returns residual of pt.

Definition at line 229 of file TTrackMC.h.

References _ptResidual.

                               {
    return _ptResidual;
}

double TTrackMC::pzFraction

(

void

)

const [inline]

returns pz fraction.

Definition at line 193 of file TTrackMC.h.

References _pzFraction.

Referenced by TTrackManager::saveMCTables().

                               {
    return _pzFraction;
}

double TTrackMC::pzPull

(

void

)

const [inline]

returns pull of pz.

Definition at line 247 of file TTrackMC.h.

References _pzPull.

                           {
    return _pzPull;
}

double TTrackMC::pzResidual

(

void

)

const [inline]

returns residual of pz.

Definition at line 235 of file TTrackMC.h.

References _pzResidual.

                               {
    return _pzResidual;
}

unsigned TTrackMC::quality

(

void

)

const [inline]

returns quality.

Definition at line 217 of file TTrackMC.h.

References _quality.

Referenced by TTrackManager::saveMCTables(), and TrackMCStatus().

                            { 
    return _quality;
} 

std::string TTrackMC::qualityString

(

void

)

const

returns quality.

Definition at line 285 of file TTrackMC.cxx.

References _quality, and TrackMCQualityString().

                                  {
    return TrackMCQualityString(_quality);
}

const Hep3Vector & TTrackMC::residual

(

void

)

const [inline]

returns residual of momentum.

Definition at line 223 of file TTrackMC.h.

References _residual.

                             {
    return _residual;
}

unsigned TTrackMC::state

(

void

)

const [inline]

returns state.

Definition at line 211 of file TTrackMC.h.

References _state.

                          {
    return _state;
}

void TTrackMC::update

(

void

)

updates information.

Definition at line 67 of file TTrackMC.cxx.

References _quality, and _state.

                     {
    _state = 0;
    _quality = 0;
/*
    //...Prepare counters...
    unsigned nHep = TTrackHEP::list().length();
    unsigned nTrk = TrkReco::getTrkReco()->tracks().length();
    unsigned * N1 = (unsigned *) malloc(nHep * sizeof(unsigned));
    float * F1 = (float *) malloc(nHep * sizeof(float));
    unsigned N2 = 0;
    for (unsigned i = 0; i < nHep; i++) {
        N1[i] = 0;
        F1[i] = 0.;
    }
    
    //...Prepare for fraction F1...
//  const AList<TMLink> & cores = _t.cores();
    const AList<TMLink> & cores = _t.finalHits();
    unsigned nCores = cores.length();
    for (unsigned i = 0; i < nCores; i++) {
        TMLink * t = cores[i];
        int hepID = t->hit()->mc()->hep()->id();
        ++N1[hepID];
    }

    //...Calculate fraction F1 and find the best HEP...
    int bestHep = -1;
    TTrackHEP * hep = 0;
    float bestF1 = 0.;
    for (unsigned i = 0; i < nHep; i++) {
        if (nCores) F1[i] = (float) N1[i] / (float) nCores;
        if (F1[i] > bestF1) {
            bestHep = i;
            bestF1 = F1[i];
        }
    }

    //...Check HEP...
    float F2 = 0.;
    if (bestHep != -1) {
        hep = TTrackHEP::list()[bestHep];
        unsigned nAll = 0;
        for (unsigned i = 0; i < hep->hits().length(); i++) {
            const TMDCWireHit * hit = hep->hits()[i]->hit();
            if (! hit) continue;
            if (hit->state() & WireHitInvalidForFit) continue;
            
            ++nAll;
            if (hit->track() == & _t) ++N2;
        }

        //...Calculate fraction F2...
        if (nAll) F2 = (float) N2 / (float) nAll;
    }

    //...Store results...
    _hepID = bestHep;
    _hep = hep;
    _wireFraction = bestF1;
    _wireFractionHEP = F2;

    //...Compare charge and momentum...
    compare();

    //...Classification...
    classify();

    //...Termination...
    free(N1);
    free(F1);
*/
}

double TTrackMC::wireFraction

(

void

)

const [inline]

returns wire fraction(F1).

Definition at line 199 of file TTrackMC.h.

References _wireFraction.

Referenced by TTrackManager::saveMCTables().

                                 {
    return _wireFraction;
}

double TTrackMC::wireFractionHEP

(

void

)

const [inline]

returns wire fraction(F2).

Definition at line 205 of file TTrackMC.h.

References _wireFractionHEP.

Referenced by TTrackManager::saveMCTables().

                                    {
    return _wireFractionHEP;
}

---

Friends And Related Function Documentation

friend class TrkReco [friend]

Definition at line 144 of file TTrackMC.h.

---

Member Data Documentation

bool TTrackMC::_charge [private]

Definition at line 134 of file TTrackMC.h.

Referenced by charge(), classify(), and compare().

double TTrackMC::_cosOpen [private]

Definition at line 138 of file TTrackMC.h.

Referenced by classify(), and compare().

TTrackHEP* TTrackMC::_hep [private]

Definition at line 130 of file TTrackMC.h.

Referenced by classify(), compare(), and hep().

int TTrackMC::_hepID [private]

Definition at line 131 of file TTrackMC.h.

Referenced by dump(), and hepId().

double TTrackMC::_ptFraction [private]

Definition at line 135 of file TTrackMC.h.

Referenced by classify(), compare(), dump(), and ptFraction().

double TTrackMC::_ptPull [private]

Definition at line 141 of file TTrackMC.h.

Referenced by compare(), and ptPull().

double TTrackMC::_ptResidual [private]

Definition at line 139 of file TTrackMC.h.

Referenced by classify(), compare(), and ptResidual().

double TTrackMC::_pzFraction [private]

Definition at line 136 of file TTrackMC.h.

Referenced by classify(), compare(), dump(), and pzFraction().

double TTrackMC::_pzPull [private]

Definition at line 142 of file TTrackMC.h.

Referenced by compare(), and pzPull().

double TTrackMC::_pzResidual [private]

Definition at line 140 of file TTrackMC.h.

Referenced by classify(), compare(), and pzResidual().

unsigned TTrackMC::_quality [private]

Definition at line 127 of file TTrackMC.h.

Referenced by classify(), dump(), quality(), qualityString(), and update().

Hep3Vector TTrackMC::_residual [private]

Definition at line 137 of file TTrackMC.h.

Referenced by compare(), and residual().

unsigned TTrackMC::_state [private]

Definition at line 126 of file TTrackMC.h.

Referenced by classify(), dump(), state(), and update().

const TTrack& TTrackMC::_t [private]

Definition at line 129 of file TTrackMC.h.

Referenced by compare(), and dump().

double TTrackMC::_wireFraction [private]

Definition at line 132 of file TTrackMC.h.

Referenced by classify(), dump(), and wireFraction().

double TTrackMC::_wireFractionHEP [private]

Definition at line 133 of file TTrackMC.h.

Referenced by dump(), and wireFractionHEP().

---