TPerfectFinder Class Reference

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

#include <TPerfectFinder.h>

Inheritance diagram for TPerfectFinder:

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.
	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.
Public Member Functions inherited from TFinderBase
	TFinderBase ()
	Constructor.
virtual	~TFinderBase ()
	Destructor.
virtual std::string	name (void) const =0
	returns name.
virtual std::string	version (void) const =0
	returns version.
virtual void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual int	debugLevel (void) const
	returns debug level.
virtual void	clear (void)=0
	clear internal information.
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.
virtual int	doit (const AList< TMDCWireHit > &axialHits, const AList< TMDCWireHit > &stereoHits, AList< TTrack > &tracks3D, AList< TTrack > &tracks2D)=0
	finds tracks. 'hits' are used to reconstruct. 'tracks' can be used for both inputs and outputs. Return value = (0, +, -) means (success, warning, fatal error).
	TFinderBase ()
	Constructor.
virtual	~TFinderBase ()
	Destructor.
virtual std::string	name (void) const =0
	returns name.
virtual std::string	version (void) const =0
	returns version.
virtual void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual int	debugLevel (void) const
	returns debug level.
virtual void	clear (void)=0
	clear internal information.
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.
virtual int	doit (const AList< TMDCWireHit > &axialHits, const AList< TMDCWireHit > &stereoHits, AList< TTrack > &tracks3D, AList< TTrack > &tracks2D)=0
	finds tracks. 'hits' are used to reconstruct. 'tracks' can be used for both inputs and outputs. Return value = (0, +, -) means (success, warning, fatal error).

Detailed Description

A class to find tracks using MC info.

Definition at line 28 of file InstallArea/include/TrkReco/TrkReco/TPerfectFinder.h.

Constructor & Destructor Documentation

TPerfectFinder() [1/2]

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() [1/2]

TPerfectFinder::~TPerfectFinder ( )

virtual

Destructor.

Definition at line 36 of file TPerfectFinder.cxx.

                                {
}

TPerfectFinder() [2/2]

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

Constructor.

~TPerfectFinder() [2/2]

virtual TPerfectFinder::~TPerfectFinder ( )

virtual

Destructor.

Member Function Documentation

clear() [1/2]

void TPerfectFinder::clear ( void  )

virtual

clear internal information.

Implements TFinderBase.

Definition at line 204 of file TPerfectFinder.cxx.

                          {
    HepAListDeleteAll(_links);
}

clear() [2/2]

void TPerfectFinder::clear ( void  )

virtual

clear internal information.

Implements TFinderBase.

doit() [1/2]

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.

                                       {
 
    //...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;
}

doit() [2/2]

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

virtual

finds tracks.

Implements TFinderBase.

dump() [1/2]

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.

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

dump() [2/2]

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

virtual

dumps debug information.

Reimplemented from TFinderBase.

name() [1/2]

std::string TPerfectFinder::name ( void  ) const

inlinevirtual

returns name.

Implements TFinderBase.

Definition at line 90 of file InstallArea/include/TrkReco/TrkReco/TPerfectFinder.h.

                               {
    return "Perfect Finder";
}

name() [2/2]

std::string TPerfectFinder::name ( void  ) const

virtual

returns name.

Implements TFinderBase.

perfectFitting() [1/2]

bool TPerfectFinder::perfectFitting ( bool  a )

inline

gets perfect momentum from MC info.

Definition at line 96 of file InstallArea/include/TrkReco/TrkReco/TPerfectFinder.h.

                                     {
    return _perfectFitting = a;
}

perfectFitting() [2/2]

bool TPerfectFinder::perfectFitting

(

bool

)

gets perfect momentum from MC info.

version() [1/2]

std::string TPerfectFinder::version ( void  ) const

virtual

returns version.

Implements TFinderBase.

Definition at line 40 of file TPerfectFinder.cxx.

                                  {
    return "2.04";
}

version() [2/2]

std::string TPerfectFinder::version ( void  ) const

virtual

returns version.

Implements TFinderBase.

---

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

• InstallArea/include/TrkReco/TrkReco/TPerfectFinder.h
• Reconstruction/TrkReco/TrkReco-00-08-59-patch4-slc6tag/TrkReco/TPerfectFinder.h
• TPerfectFinder.cxx