TPerfectFinder Class Reference

A class to find tracks using MC info. More...

#include <TPerfectFinder.h>

Inheritance diagram for TPerfectFinder:

List of all members.

Public Member Functions
	TPerfectFinder (int perfectFitting, float maxSigma, float maxSigmaStereo, unsigned fittingFlag)
	Constructor.
virtual	~TPerfectFinder ()
	Destructor.
std::string	name (void) const
	returns name.
std::string	version (void) const
	returns version.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
void	clear (void)
	clear internal information.
bool	perfectFitting (bool)
	gets perfect momentum from MC info.
int	doit (const AList< TMDCWireHit > &axialHits, const AList< TMDCWireHit > &stereoHits, AList< TTrack > &tracks, AList< TTrack > &tracks2D)
	finds tracks.
virtual int	debugLevel (void) const
	returns debug level.
virtual int	debugLevel (int)
	sets debug level.
virtual bool	doStereo (bool)
	sets flag to reconstruct 3D.
virtual bool	doSalvage (bool)
	sets flag to salvage hits.
Private Member Functions
float	charge (int pType) const
Private Attributes
bool	_perfectFitting
AList< TMLink >	_links
float	_maxSigma
float	_maxSigmaStereo
unsigned	_fittingFlag
TBuilder	_builder
THelixFitter	_fitter

---

Detailed Description

A class to find tracks using MC info.

Definition at line 28 of file TPerfectFinder.h.

---

Constructor & Destructor Documentation

TPerfectFinder::TPerfectFinder	(	int	perfectFitting,
		float	maxSigma,
		float	maxSigmaStereo,
		unsigned	fittingFlag
	)

Constructor.

Definition at line 21 of file TPerfectFinder.cxx.

: _perfectFitting(perfectFitting),
  _builder("conformal builder",
           maxSigma,
           maxSigmaStereo,
           0,
           0,
           0,
           fittingFlag),
  _fitter("helix fitter") {
}

TPerfectFinder::~TPerfectFinder

(

)

[virtual]

Destructor.

Definition at line 36 of file TPerfectFinder.cxx.

                                {
}

---

Member Function Documentation

float TPerfectFinder::charge

(

int

pType

)

const [private]

Definition at line 209 of file TPerfectFinder.cxx.

Referenced by doit().

                                      {
    float chg;

    if (pType == 11) chg = -1;
    else if (pType == -11) chg = 1;
    else if (pType == 13) chg = -1;
    else if (pType == -13) chg = 1;
    else if (pType == 211) chg = 1;
    else if (pType == -211) chg = -1;
    else if (pType == 321) chg = 1;
    else if (pType == -321) chg = -1;
    else if (pType == 2212) chg = 1;
    else if (pType == -2212) chg = -1;
    else {
        std::cout << "TPerfectFinder !!! charge of particle type=";
        std::cout << pType << " is unknown" << std::endl;
    }

    return chg;
}

void TPerfectFinder::clear

(

void

)

[virtual]

clear internal information.

Implements TFinderBase.

Definition at line 204 of file TPerfectFinder.cxx.

References _links.

                          {
    HepAListDeleteAll(_links);
}

int TFinderBase::debugLevel

(

int

)

[inline, virtual, inherited]

sets debug level.

Definition at line 96 of file TFinderBase.h.

References TFinderBase::_debugLevel.

                             {
    return _debugLevel = a;
}

int TFinderBase::debugLevel

(

void

)

const [inline, virtual, inherited]

returns debug level.

Definition at line 90 of file TFinderBase.h.

References TFinderBase::_debugLevel.

Referenced by TConformalFinder0::doit(), TConformalFinder::doit(), TrkReco::initialize(), TConformalFinder0::specialFinding(), and TConformalFinder0::standardFinding().

                                  {
    return _debugLevel;
}

int TPerfectFinder::doit	(	const AList< TMDCWireHit > &	axialHits,
		const AList< TMDCWireHit > &	stereoHits,
		AList< TTrack > &	tracks,
		AList< TTrack > &	tracks2D
	)			[virtual]

finds tracks.

Implements TFinderBase.

Definition at line 51 of file TPerfectFinder.cxx.

References _fitter, _links, _perfectFitting, charge(), THelixFitter::fit(), TMDC::getTMDC(), TMDC::hitsMC(), genRecEmupikp::i, ganga-rec::j, TMLink::leftRight(), Helix::pivot(), t(), TMLink::wire(), and TMDCWire::xyPosition().

                                               {

    //...Preparations...
    static const HepPoint3D dummy(0, 0, 0);
    CAList<TMDCWireHit> hits;
    hits.append(axialHits);
    hits.append(stereoHits);

    //...Make a list of HEP track...
    CAList<TTrackHEP> heps;
    const unsigned nHits = hits.length();
    for (unsigned i = 0; i < nHits; i++) {
        const TMDCWireHitMC * const mc = hits[i]->mc();
        if (! mc) continue;
        const TTrackHEP * const hep = mc->hep();

        if (! hep) continue;
        heps.append(hep);
    }
    heps.purge();

    //...Create tracks...
    const AList<TMDCWireHitMC> & mcHits = TMDC::getTMDC()->hitsMC();
    const unsigned nHitsMC = mcHits.length();
    const unsigned nHeps = heps.length();
    for (unsigned i = 0; i < nHeps; i++) {
        const TTrackHEP & hep = * heps[i];
        const float chg = charge(hep.pType());
        AList<TMLink> hepLinks;
        HepPoint3D posIn, posOut;
        HepVector3D momIn, momOut;
        float rMin = 99999;
        float rMax = 0.;
        TMLink * linkIn;
        unsigned lastLayer = 0;
        for (unsigned j = 0; j < mcHits.length(); j++) {
            const TMDCWireHitMC * const mc = mcHits[j];
            if (! mc) continue;
            if (! mc->hit()) continue;
            if (& hep != mc->hep()) continue;
            if (! hits.hasMember(mc->hit())) continue;

            //...Remove hits by curl back...(assuming order of mc hits)
            if (mc->wire()->layerId() < lastLayer) break;
            lastLayer = mc->wire()->layerId();
            HepPoint3D ent = mc->entrance();
            HepVector3D dir = mc->direction();
            ent.setZ(0.);
            dir.setZ(0.);
            if ((ent.unit()).dot(dir.unit()) < 0.5) continue;

            //...Hit...
            TMLink * l = new TMLink(0, mc->hit(), dummy);
            l->leftRight(mc->leftRight());
            hepLinks.append(l);
            _links.append(l);

            //...Check r to get MC mom...
            const float r = ent.mag();
            if (r < rMin) {
                posIn = mc->entrance();
                momIn = mc->momentum();
                rMin = r;
                linkIn = l;
            }
            if (r > rMax) {
                posOut = mc->entrance();
                momOut = mc->momentum();
                rMax = r;
            }
        }
        if (_links.length() == 0) continue;

        //...Do perfect fitting...
        TTrack * t = 0;
//          Helix h(posIn, momIn, chg);
        Helix h(posOut, momOut, chg);
//          h.pivot(posIn);
        h.pivot(linkIn->wire()->xyPosition());
        t = new TTrack(h);
        t->append(hepLinks);
        t->assign(0);

//          Helix hX(posOut, momOut, chg);
        
//          std::cout << "    gen@cdc i = " << h.a() << std::endl;
//          std::cout << "       pivot  = " << h.pivot() << std::endl;
//          std::cout << "       mom    = " << h.momentum() << std::endl;
//          std::cout << "    gen@cdc o = " << hX.a() << std::endl;
//          std::cout << "       pivot  = " << hX.pivot() << std::endl;
//          std::cout << "       mom    = " << hX.momentum() << std::endl;
//          std::cout << "    momOut    = " << momOut << std::endl;
//          std::cout << "    posOut    = " << posOut << std::endl;
//          std::cout << "    ptOut     = " << momOut.perp() << std::endl;
//          std::cout << "    costhOut  = " << momOut.z() / momOut.mag();
//          std::cout << std::endl;

        if (_perfectFitting == 0) {
            std::cout << "special test in perfect finder" << std::endl;
            AList<TMLink> tmp;
            for (unsigned i = 0; i < t->nLinks(); i++) {
                bool rem = false;
                if (t->links()[i]->wire()->name() == "0-56")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "2-56")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "5-57")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "6=68")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "7=69")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "8=68")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "10-85")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "20-126")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "39-210")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "41=221")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "43=221")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "46-251")
                    rem = true;
                else if (t->links()[i]->wire()->name() == "48-251")
                    rem = true;
                if (rem) {
                    std::cout << t->links()[i]->wire()->name() << " removed" << std::endl;
                    tmp.append(* t->links()[i]);
                }
            }
            t->remove(tmp);
            t->dump("detail","@lastHit");
            Helix & h = (Helix &) t->helix();
            h.pivot(HepVector3D(0., 0., 0.));
            t->dump("detail","@origin ");
            _fitter.fit(* t);
            t->dump("detail","fitted  ");
            static const HepVector3D p0(1.226, -1.025, 0.120);
            std::cout << "Pdif mag=" << (t->p() - p0).mag() << std::endl;
        }

        tracks.append(t);
    }
    return 0;
}

bool TFinderBase::doSalvage

(

bool

)

[inline, virtual, inherited]

sets flag to salvage hits.

Reimplemented in TConformalFinder0.

Definition at line 108 of file TFinderBase.h.

Referenced by TrkReco::initialize().

                             {
    return a;
}

bool TFinderBase::doStereo

(

bool

)

[inline, virtual, inherited]

sets flag to reconstruct 3D.

Reimplemented in TConformalFinder0.

Definition at line 102 of file TFinderBase.h.

Referenced by TrkReco::initialize().

                            {
    return a;
}

void TPerfectFinder::dump	(	const std::string &	message = std::string(""),
		const std::string &	prefix = std::string("")
	)			const [virtual]

dumps debug information.

Reimplemented from TFinderBase.

Definition at line 45 of file TPerfectFinder.cxx.

References TFinderBase::dump().

                                                                       {
    std::cout << pre;
    TFinderBase::dump(msg);
}

std::string TPerfectFinder::name

(

void

)

const [inline, virtual]

returns name.

Implements TFinderBase.

Definition at line 90 of file TPerfectFinder.h.

                               {
    return "Perfect Finder";
}

bool TPerfectFinder::perfectFitting

(

bool

)

[inline]

gets perfect momentum from MC info.

Definition at line 96 of file TPerfectFinder.h.

References _perfectFitting.

                                     {
    return _perfectFitting = a;
}

std::string TPerfectFinder::version

(

void

)

const [virtual]

returns version.

Implements TFinderBase.

Definition at line 40 of file TPerfectFinder.cxx.

                                  {
    return "2.04";
}

---

Member Data Documentation

TBuilder TPerfectFinder::_builder [private]

Definition at line 73 of file TPerfectFinder.h.

THelixFitter TPerfectFinder::_fitter [private]

Definition at line 74 of file TPerfectFinder.h.

Referenced by doit().

unsigned TPerfectFinder::_fittingFlag [private]

Definition at line 72 of file TPerfectFinder.h.

AList<TMLink> TPerfectFinder::_links [private]

Definition at line 69 of file TPerfectFinder.h.

Referenced by clear(), and doit().

float TPerfectFinder::_maxSigma [private]

Definition at line 70 of file TPerfectFinder.h.

float TPerfectFinder::_maxSigmaStereo [private]

Definition at line 71 of file TPerfectFinder.h.

bool TPerfectFinder::_perfectFitting [private]

Definition at line 68 of file TPerfectFinder.h.

Referenced by doit(), and perfectFitting().

---