TMDC Class Reference

#include <TMDC.h>

Inheritance diagram for TMDC:

Public Member Functions
std::string	name (void) const
	returns name.
std::string	version (void) const
	returns version.
std::string	cdcVersion (void) const
	returns MDC version.
void	dump (const std::string &message) const
	dumps debug information.
int	debugLevel (void) const
	returns debug level.
float	fudgeFactor (void) const
	returns fudge factor for drift time error.
void	clear (void)
	clears all TMDC information.
void	fastClear (void)
	clears TMDC information.
void	update (bool mcAnalysis=true)
	updates TMDC information. clear() is called in this function.
int	debugLevel (int)
	sets debug level.
float	fudgeFactor (float)
	sets/returns fudge factor for drift time error.
const TMDCWire *const	wire (unsigned wireId) const
	returns a pointer to a wire. 0 will be returned if 'wireId' is invalid.
const TMDCWire *const	wire (unsigned layerId, int localId) const
	returns a pointer to a wire. 'localId' can be negative. 0 will be returned if 'layerId' is invalid.
const TMDCWire *const	wire (const HepPoint3D &point) const
	returns a pointer to a wire.
const TMDCWire *const	wire (float r, float phi) const
	returns a pointer to a wire.
const TMDCLayer *const	layer (unsigned id) const
	returns a pointer to a layer. 0 will be returned if 'id' is invalid.
const AList< TMDCLayer > *const	superLayer (unsigned id) const
	returns a pointer to a super-layer. 0 will be returned if 'id' is invalid.
unsigned	nLocalLayer (unsigned superLayerId) const
	returns # of wire layers in a super layer. 0 will be returned if 'superLayerId' is invalid.
const AList< TMDCWireHit > &	hits (unsigned mask=0) const
	returns a list of TMDCWireHit. 'update()' must be called before calling this function.
const AList< TMDCWireHit > &	axialHits (unsigned mask=0) const
	returns a list of axial hits. 'update()' must be called before calling this function.
const AList< TMDCWireHit > &	stereoHits (unsigned mask=0) const
	returns a list of stereo hits. 'update()' must be called before calling this function.
const AList< TMDCWireHitMC > &	hitsMC (void) const
	returns a list of TMDCWireHitMC. 'updateMC()' must be called before calling this function.
const AList< TMDCWireHit > &	badHits (void)
	returns bad hits(finding invalid hits).
unsigned	nWires (void) const
unsigned	nSuperLayers (void) const
unsigned	nLayers (void) const
void	classification (void)
	classify hits.
void	updateMC (void)
	updates TMDC information for MC.

Static Public Member Functions
static TMDC *	getTMDC (const std::string &cdcVersion)
	returns TMDC object.
static TMDC *	getTMDC (void)
static std::string	wireName (unsigned wireId)
static std::string	wireName (const MdcGeoWire *const)
static std::string	wireName (const MdcRec_wirhit &)
static unsigned	layerId (unsigned wireId)
static unsigned	layerId (const MdcGeoWire *const)
static unsigned	localId (unsigned wireId)
static unsigned	localId (const MdcGeoWire *const)
static unsigned	superLayerId (unsigned wireId)
static unsigned	superLayerId (const MdcGeoWire *const)
static unsigned	superLayerId (const MdcGeoLayer *const)
static unsigned	localLayerId (unsigned wireId)
static unsigned	localLayerId (const MdcGeoWire *const)
static unsigned	localLayerId (const MdcGeoLayer *const)
static unsigned	axialStereoLayerId (const MdcGeoLayer *const)
static unsigned	layerId (unsigned axialStereo, unsigned axialStereoLayerId)
static float	phi (float phi)
static void	driftDistance (TMLink &link, const TTrack &track, unsigned correctionFlag=0, float T0Offset=0.)
	calculates corrected drift time. correctionFlag(bit 0:tof, 1:T0 offset, 2:propagation delay, 3:tan(lambda) correction)

Detailed Description

A class to represent MDC.

The instance of TMDC is a singleton. 'TMDC::getTMDC()' gives you a pointer to access the TMDC. Geometrical information is initialized automatically. Before accessing hit information, user has to call 'update()' to initialize hit information event by event.

Definition at line 61 of file TMDC.h.

Member Function Documentation

axialHits()

const AList< TMDCWireHit > & TMDC::axialHits

(

unsigned

mask = 0

)

const

returns a list of axial hits. 'update()' must be called before calling this function.

Definition at line 518 of file TMDC.cxx.

                                   {
    if (! mask) return _axialHits;
    else if (mask == WireHitFindingValid) return _axialHits;
    std::cout << "TMDC::axialHits !!! unsupported mask given" << std::endl;
    return _axialHits;
}

Referenced by TrkReco::execute().

axialStereoLayerId()

unsigned TMDC::axialStereoLayerId ( const MdcGeoLayer * const l )

static

Definition at line 991 of file TMDC.cxx.

                                                    {
    if (! l) return 9999;
    unsigned id = l->Id();
 
    if (id < 8) return id;
    else if (id < 20) return id-8;
    else if (id < 36) return id-12;
    else if (id < 43) return id-24;
 
    return 9999;    
}

Referenced by TMDCLayer::TMDCLayer().

badHits()

const AList< TMDCWireHit > & TMDC::badHits

(

void

)

returns bad hits(finding invalid hits).

Definition at line 542 of file TMDC.cxx.

                  {
    if (! updated()) update();
    if (_badHits.length()) return _badHits;
 
    //...Loop over RECMDC_WIRHIT bank...
//    unsigned nReccdc = BsCouTab(RECMDC_WIRHIT);
    unsigned nReccdc = MdcRecWirhitCol::getMdcRecWirhitCol()->size();
    for (unsigned i = 0; i < nReccdc; i++) {
//  struct reccdc_wirhit * h =
//      (struct reccdc_wirhit *)
//      BsGetEnt(RECMDC_WIRHIT, i + 1, BBS_No_Index);
    MdcRec_wirhit* h = &(*MdcRecWirhitCol::getMdcRecWirhitCol())[i];
    
    //...Check validity...
    if (h->stat & WireHitFindingValid) continue;
 
    //...Obtain a pointer to GEOMDC...
//  geocdc_wire * g =
//      (geocdc_wire *) BsGetEnt(GEOMDC_WIRE, h->m_geo, BBS_No_Index);
    const MdcGeoWire* g = h->geo;
    
    //...Get a pointer to a TMDCWire...
//  TMDCWire * w = _wires[g->m_ID - 1];
    TMDCWire * w = _wires[g->Id()];
    
    //...Create TMDCWireHit...
    _badHits.append(new TMDCWireHit(w, h, _fudgeFactor));
    }
 
    return _badHits;
}

cdcVersion()

std::string TMDC::cdcVersion ( void ) const

inline

returns MDC version.

Definition at line 304 of file TMDC.h.

                           {
    return std::string(_cdcVersion);
}

Referenced by dump().

classification()

void TMDC::classification

(

void

)

classify hits.

Definition at line 445 of file TMDC.cxx.

                         {
    unsigned n = _hits.length();
 
    for (unsigned i = 0; i < n; i++) {
    TMDCWireHit * h = _hits[i];
    const TMDCWire * w = h->wire();
    unsigned state = h->state();
 
    //...Cache pointers to a neighbor...
    const TMDCWire * neighbor[6];
    for (unsigned j = 0; j < 6; j++) neighbor[j] = w->neighbor(j);
 
    //output to test ...
//  cout<<w->localId()<<endl;
//  if(neighbor[0]) cout<<"inner local:"<<neighbor[0]->localId()<<", "<<neighbor[1]->localId()<<endl;
//  if(neighbor[5]) cout<<"outer local:"<<neighbor[4]->localId()<<", "<<neighbor[5]->localId()<<endl;   
        
    //...Decide hit pattern...
    unsigned pattern = 0;
    for (unsigned j = 0; j < 6; j++) {
        if (neighbor[j])
        if (neighbor[j]->hit())
            pattern += (1 << j);
    }
    state |= (pattern << WireHitNeighborHit);
 
    //...Check isolation...
    const TMDCWireHit * hr1 = neighbor[2]->hit();
    const TMDCWireHit * hl1 = neighbor[3]->hit();
    if ((hr1 == 0) && (hl1 == 0)) {
        state |= WireHitIsolated;
    }
    else {
        const TMDCWireHit * hr2 = neighbor[2]->neighbor(2)->hit();
        const TMDCWireHit * hl2 = neighbor[3]->neighbor(3)->hit();
        if ((hr2 == 0) && (hr1 != 0) && (hl1 == 0) ||
        (hl2 == 0) && (hl1 != 0) && (hr1 == 0))
        state |= WireHitIsolated;
    }
 
    //...Check continuation...
    unsigned superLayer = w->superLayerId();
    bool previous = false;
    bool next = false;
    if (neighbor[0] == 0) previous = true;
    else {
        if ((neighbor[0]->hit()) || neighbor[1]->hit())
        previous = true;
    }
    if (neighbor[5] == 0) next = true;
    else {
        if ((neighbor[4]->hit()) || neighbor[5]->hit())
        next = true;
    }
    // if (previous && next) state |= WireHitContinuous;
    if (previous || next) state |= WireHitContinuous;
 
    //...Solve LR locally...
    if ((pattern == 34) || (pattern == 42) ||
        (pattern == 40) || (pattern == 10) ||
        (pattern == 35) || (pattern == 50))
        state |= WireHitPatternRight;
    else if ((pattern == 17) || (pattern == 21) ||
         (pattern == 20) || (pattern ==  5) ||
         (pattern == 19) || (pattern == 49))
        state |= WireHitPatternLeft;
 
    //...Store it...
    h->state(state);
    }
}

Referenced by update().

clear()

void TMDC::clear ( void )

virtual

clears all TMDC information.

Reimplemented from TUpdater.

Definition at line 282 of file TMDC.cxx.

                {
    unsigned i = 0;
    while (TMDCWire * w = _wires[i++])
    w->clear();
 
    _hitWires.removeAll();
    _axialHits.removeAll();
    _stereoHits.removeAll();
    HepAListDeleteAll(_hits);
    HepAListDeleteAll(_hitsMC);
    HepAListDeleteAll(_badHits);
 
    TUpdater::clear();
}

debugLevel() [1/2]

int TMDC::debugLevel ( int a )

inline

sets debug level.

Definition at line 245 of file TMDC.h.

                      {
    return _debugLevel = a;
}

debugLevel() [2/2]

int TMDC::debugLevel ( void ) const

inline

returns debug level.

Definition at line 239 of file TMDC.h.

                           {
    return _debugLevel;
}

driftDistance()

void TMDC::driftDistance ( TMLink & link, const TTrack & track, unsigned correctionFlag = 0, float T0Offset = 0. )

static

calculates corrected drift time. correctionFlag(bit 0:tof, 1:T0 offset, 2:propagation delay, 3:tan(lambda) correction)

Definition at line 1029 of file TMDC.cxx.

                            {
 
    //...No correction...
    if (flag == 0) {
    if (l.hit()) {
        l.drift(0, l.hit()->drift(0));
        l.drift(1, l.hit()->drift(1));
        l.dDrift(0, l.hit()->dDrift(0));
        l.dDrift(1, l.hit()->dDrift(1));
    }
    else {
        l.drift(0, 0.);
        l.drift(1, 0.);
        l.dDrift(0, 0.);
        l.dDrift(1, 0.);
    }
 
    return;
    }
 
    //...TOF correction...
    float tof = 0.;
    if (flag && 1) {
    int imass = 3;
    float tl = t.helix().a()[4];
    float f = sqrt(1. + tl * tl);
    float s = fabs(t.helix().curv()) * fabs(l.dPhi()) * f;
    float p = f / fabs(t.helix().a()[2]);
//zsl   calcdc_tof2_(& imass, & p, & s, & tof); 
    }
 
    //...T0 correction....
    if (! (flag && 2)) t0Offset = 0.;
 
    //...Propagation corrections...
    const TMDCWireHit & h = * l.hit();
    int wire = h.wire()->id();
    HepVector3D tp = t.helix().momentum(l.dPhi());
    float p[3] = {tp.x(), tp.y(), tp.z()};
    const HepPoint3D & onWire = l.positionOnWire();
    float x[3] = {onWire.x(), onWire.y(), onWire.z()};
    float time = h.reccdc()->tdc + t0Offset - tof;
    float dist;
    float edist;
    int prop = (flag & 4);
 
    //...Calculation with left side...
    int side = -1;
    if (side == 0) side = -1;
//zsl    calcdc_driftdist_(& prop,
//            & wire,
//            & side,
//            p,
//            x,
//            & time,
//            & dist,
//            & edist);
    l.drift(0, dist);
    l.dDrift(0, edist);
 
    //...Calculation with left side...
    side = 1;
//zsl    calcdc_driftdist_(& prop,
//            & wire,
//            & side,
//            p,
//            x,
//            & time,
//            & dist,
//            & edist);
    l.drift(1, dist);
    l.dDrift(1, edist);
 
    //...tan(lambda) correction...
    if (flag && 8) {
    float tanl = abs(p[2]) / tp.perp();
    float c;
    if ((tanl >= 0.0) && (tanl < 0.5))      c = -0.48 * tanl + 1.3;
    else if ((tanl >= 0.5) && (tanl < 1.0)) c = -0.28 * tanl + 1.2;
    else if ((tanl >= 1.0) && (tanl < 1.5)) c = -0.16 * tanl + 1.08;
    else                                    c =  0.84;
 
    l.dDrift(0, l.dDrift(0) * c);
    l.dDrift(1, l.dDrift(1) * c);
    }
}

dump()

void TMDC::dump

(

const std::string &

message

)

const

dumps debug information.

Definition at line 181 of file TMDC.cxx.

                                      {
//      if (   msg.index("name")    != -1
//      || msg.index("version") != -1
//      || msg.index("detail")  != -1
//      || msg == "") {
  if (   msg.find("name")    != std::string::npos
    || msg.find("version") != std::string::npos
    || msg.find("detail")  != std::string::npos
    || msg == "") {
    std::cout << name() << "(" << version() << ") ";
    }
    if (msg.find("detail") != std::string::npos || msg.find("state") != std::string::npos) {
    std::cout << "Debug Level=" << _debugLevel;
    }
    std::cout << std::endl;
 
    std::string tab("        ");
 
    if (msg == "" || msg.find("geometry") != std::string::npos) {
 
    //...Get information..."
    unsigned nLayer = _layers.length();
    std::cout << "    version    : " << version() << std::endl;
    std::cout << "    cdc version: " << cdcVersion() << std::endl;
    std::cout << "    # of wires : " << _wires.length() << std::endl;
    std::cout << "    # of layers: " << nLayer << std::endl;
    std::cout << "    super layer information" << std::endl;
    std::cout << "        # of super layers = " << _superLayers.length()
          << std::endl;
 
    std::cout << "    layer information" << std::endl;
    for (unsigned i = 0; i < _nLayers; i++)
        _layers[i]->dump("", tab);
 
    std::cout << "    wire information" << std::endl;
    for (unsigned i = 0; i < _nWires; i++)
        (_wires[i])->dump("neighbor", tab);
 
    return;
    }
    if (msg.find("hits") != std::string::npos) {
    std::cout << "    hits     : " << _hits.length() << std::endl;
    for (unsigned i = 0; i < _hits.length(); i++) {
        _hits[i]->dump("state mc", tab);
    }
    }
    if (msg.find("axialHits") != std::string::npos) {
    std::cout << "    hits     : " << _axialHits.length() << std::endl;
    for (unsigned i = 0; i < _axialHits.length(); i++) {
        _axialHits[i]->dump("state mc", tab);
    }
    }
    if (msg.find("stereoHits") != std::string::npos) {
    std::cout << "    hits     : " << _stereoHits.length() << std::endl;
    for (unsigned i = 0; i < _stereoHits.length(); i++) {
        _stereoHits[i]->dump("state mc", tab);
    }
    }
}

Referenced by dump().

fastClear()

void TMDC::fastClear

(

void

)

clears TMDC information.

Definition at line 298 of file TMDC.cxx.

                    {
    unsigned i = 0;
    while (TMDCWire * w = _hitWires[i++])
    w->clear();
    i = 0;
    while (TMDCWireHit * h = _badHits[i++])
    ((TMDCWire *) h->wire())->clear();
 
    _hitWires.removeAll();
    _axialHits.removeAll();
    _stereoHits.removeAll();
    HepAListDeleteAll(_hits);
    HepAListDeleteAll(_hitsMC);
    HepAListDeleteAll(_badHits);
 
    TUpdater::clear();
}

Referenced by TrkReco::finalize(), and update().

fudgeFactor() [1/2]

float TMDC::fudgeFactor ( float a )

inline

sets/returns fudge factor for drift time error.

Definition at line 298 of file TMDC.h.

                         {
    return _fudgeFactor = a;
}

fudgeFactor() [2/2]

float TMDC::fudgeFactor ( void ) const

inline

returns fudge factor for drift time error.

Definition at line 292 of file TMDC.h.

                            {
    return _fudgeFactor;
}

getTMDC() [1/2]

TMDC * TMDC::getTMDC ( const std::string & cdcVersion )

static

returns TMDC object.

Definition at line 89 of file TMDC.cxx.

                                       {
    if (! _cdc) _cdc = new TMDC(version);
    return _cdc;
}

getTMDC() [2/2]

TMDC * TMDC::getTMDC ( void )

static

Definition at line 95 of file TMDC.cxx.

                  {
    if (! _cdc) _cdc = new TMDC("1.0");
    return _cdc;
}

Referenced by TBuilder0::buildStereo(), TConformalFinder::conformal2Wire(), TMDCTsf::createTsf(), TPerfectFinder::doit(), and TMDCWire::neighbor().

hits()

const AList< TMDCWireHit > & TMDC::hits

(

unsigned

mask = 0

)

const

returns a list of TMDCWireHit. 'update()' must be called before calling this function.

Definition at line 534 of file TMDC.cxx.

                              {
    if (! mask) return _hits;
    else if (mask == WireHitFindingValid) return _hits;
    std::cout << "TMDC::hits !!! unsupported mask given" << std::endl;
    return _hits;
}

Referenced by TrkReco::execute().

hitsMC()

const AList< TMDCWireHitMC > & TMDC::hitsMC ( void ) const

inline

returns a list of TMDCWireHitMC. 'updateMC()' must be called before calling this function.

Definition at line 265 of file TMDC.h.

                       {
    return _hitsMC;
}

Referenced by TPerfectFinder::doit().

layer()

const TMDCLayer *const TMDC::layer ( unsigned id ) const

inline

returns a pointer to a layer. 0 will be returned if 'id' is invalid.

Definition at line 259 of file TMDC.h.

                             {
    return _layers[id];
}

Referenced by TMDCTsf::createTsf(), TMDCWire::neighbor(), and wire().

layerId() [1/3]

unsigned TMDC::layerId ( const MdcGeoWire * const w )

static

Definition at line 661 of file TMDC.cxx.

                                        {
    if (! w) return 9999;
    //zsl    unsigned id = w->Id();
    //    return layerId(id);
    return w->Layer();
}

layerId() [2/3]

unsigned TMDC::layerId ( unsigned axialStereo, unsigned axialStereoLayerId )

static

Definition at line 1004 of file TMDC.cxx.

                                      {
 
    //...Axial case...
    if (as == 0) {
    if (id < 12) return id + 8;
    else if (id < 19) return id + 24;
    return 9999;
    }
 
    //...Stereo case...
    if (id < 8) return id;
    else if (id < 24) return id + 12;
    return 9999;
}

layerId() [3/3]

unsigned TMDC::layerId ( unsigned wireId )

static

Definition at line 598 of file TMDC.cxx.

                         {
//    int length = _wires.length();
//    if(id < length) return _wires[id]->geocdc()->Layer();
//    return 9999;    
 
    if (id < 40) return 0;
    else if (id < 84) return 1;
    else if (id < 132) return 2;
    else if (id < 188) return 3;
 
    else if (id < 252) return 4;
    else if (id < 324) return 5;
    else if (id < 404) return 6;
    else if (id < 484) return 7;
 
    else if (id < 560) return 8;
    else if (id < 636) return 9;
    else if (id < 724) return 10;
    else if (id < 812) return 11;
 
    else if (id < 912) return 12;
    else if (id < 1012) return 13;
    else if (id < 1124) return 14;
    else if (id < 1236) return 15;
 
    else if (id < 1364) return 16;
    else if (id < 1492) return 17;
    else if (id < 1632) return 18;
    else if (id < 1772) return 19;
 
    else if (id < 1932) return 20;
    else if (id < 2092) return 21;
    else if (id < 2252) return 22;
    else if (id < 2412) return 23;
 
    else if (id < 2604) return 24;
    else if (id < 2796) return 25;
    else if (id < 2988) return 26;
    else if (id < 3180) return 27;
 
    else if (id < 3388) return 28;
    else if (id < 3596) return 29;
    else if (id < 3804) return 30;
    else if (id < 4012) return 31;
 
    else if (id < 4252) return 32;
    else if (id < 4492) return 33;
    else if (id < 4732) return 34;
    else if (id < 4972) return 35;
 
    else if (id < 5228) return 36;
    else if (id < 5484) return 37;
    else if (id < 5740) return 38;
    else if (id < 5996) return 39;
 
    else if (id < 6284) return 40;
    else if (id < 6572) return 41;
    else if (id < 6860) return 42;
 
    return 9999;    
}

Referenced by localLayerId(), wire(), and wireName().

localId() [1/2]

unsigned TMDC::localId ( const MdcGeoWire * const w )

static

Definition at line 750 of file TMDC.cxx.

                                        {
    if (! w) return 9999;
//    unsigned id = w->Id();
//    return localId(id);
    return w->Cell();    
}

localId() [2/2]

unsigned TMDC::localId ( unsigned wireId )

static

Definition at line 676 of file TMDC.cxx.

                         {
    if (id < 40) return id;
    else if (id < 84) return id-40;
    else if (id < 132) return id-84;
    else if (id < 188) return id-132;
 
    else if (id < 252) return id-188;
    else if (id < 324) return id-252;
    else if (id < 404) return id-324;
    else if (id < 484) return id-404;
 
    else if (id < 560) return id-484;
    else if (id < 636) return id-560;
    else if (id < 724) return id-636;
    else if (id < 812) return id-724;
 
    else if (id < 912) return id-812;
    else if (id < 1012) return id-912;
    else if (id < 1124) return id-1012;
    else if (id < 1236) return id-1124;
 
    else if (id < 1364) return id-1236;
    else if (id < 1492) return id-1364;
    else if (id < 1632) return id-1492;
    else if (id < 1772) return id-1632;
 
    else if (id < 1932) return id-1772;
    else if (id < 2092) return id-1932;
    else if (id < 2252) return id-2092;
    else if (id < 2412) return id-2252;
 
    else if (id < 2604) return id-2412;
    else if (id < 2796) return id-2604;
    else if (id < 2988) return id-2796;
    else if (id < 3180) return id-2988;
 
    else if (id < 3388) return id-3180;
    else if (id < 3596) return id-3388;
    else if (id < 3804) return id-3596;
    else if (id < 4012) return id-3804;
 
    else if (id < 4252) return id-4012;
    else if (id < 4492) return id-4252;
    else if (id < 4732) return id-4492;
    else if (id < 4972) return id-4732;
 
    else if (id < 5228) return id-4972;
    else if (id < 5484) return id-5228;
    else if (id < 5740) return id-5484;
    else if (id < 5996) return id-5740;
 
    else if (id < 6284) return id-5996;
    else if (id < 6572) return id-6284;
    else if (id < 6860) return id-6572;
 
    return 9999;    
 
/*    if (id < 384) return id % 64;
    else if (id < 624) return (id - 384) % 80;
    else if (id < 1200) return (id - 624) % 96;
    else if (id < 1584) return (id - 1200) % 128;
    else if (id < 2304) return (id - 1584) % 144;
    else if (id < 2944) return (id - 2304) % 160;
    else if (id < 3904) return (id - 2944) % 176;
    else if (id < 4736) return (id - 3904) % 208;
    else if (id < 5936) return (id - 4736) % 240;
    else if (id < 6960) return (id - 5936) % 256;
    else if (id < 8400) return (id - 6960) % 288;
 
    return 9999;
*/
}

Referenced by wire(), wire(), and wireName().

localLayerId() [1/3]

unsigned TMDC::localLayerId ( const MdcGeoLayer * const l )

static

Definition at line 919 of file TMDC.cxx.

                                              {
    if (! l) return 9999;
    unsigned id = l->Id();
    if (id < 44) return id%4;      
    return 9999;
    
/*    if (id < 1) return 0;
    else if (id < 2) return 1;
    else if (id < 3) return 2;
    else if (id < 4) return 3;
    else if (id < 5) return 4;
    else if (id < 6) return 5;
 
    else if (id < 7) return 0;
    else if (id < 8) return 1;
    else if (id < 9) return 2;
 
    else if (id < 10) return 0;
    else if (id < 11) return 1;
    else if (id < 12) return 2;
    else if (id < 13) return 3;
    else if (id < 14) return 4;
    else if (id < 15) return 5;
 
    else if (id < 16) return 0;
    else if (id < 17) return 1;
    else if (id < 18) return 2;
 
    else if (id < 19) return 0;
    else if (id < 20) return 1;
    else if (id < 21) return 2;
    else if (id < 22) return 3;
    else if (id < 23) return 4;
 
    else if (id < 24) return 0;
    else if (id < 25) return 1;
    else if (id < 26) return 2;
    else if (id < 27) return 3;
 
    else if (id < 28) return 0;
    else if (id < 29) return 1;
    else if (id < 30) return 2;
    else if (id < 31) return 3;
    else if (id < 32) return 4;
 
    else if (id < 33) return 0;
    else if (id < 34) return 1;
    else if (id < 35) return 2;
    else if (id < 36) return 3;
 
    else if (id < 37) return 0;
    else if (id < 38) return 1;
    else if (id < 39) return 2;
    else if (id < 40) return 3;
    else if (id < 41) return 4;
 
    else if (id < 42) return 0;
    else if (id < 43) return 1;
    else if (id < 44) return 2;
    else if (id < 45) return 3;
 
    else if (id < 46) return 0;
    else if (id < 47) return 1;
    else if (id < 48) return 2;
    else if (id < 49) return 3;
    else if (id < 50) return 4;
 
    return 9999;    
*/
}

localLayerId() [2/3]

unsigned TMDC::localLayerId ( const MdcGeoWire * const w )

static

Definition at line 903 of file TMDC.cxx.

                                             {
    if (! w) return 9999;
    return w->Layer()%4;
    
//    unsigned id = w->Id();
//    return localLayerId(id);
}

localLayerId() [3/3]

unsigned TMDC::localLayerId ( unsigned wireId )

static

Definition at line 833 of file TMDC.cxx.

                              {
    if(id < 6860) return layerId(id)%4;
    return 9999;
    
/*    if (id < 64) return 0;
    else if (id < 128) return 1;
    else if (id < 192) return 2;
    else if (id < 256) return 3;
    else if (id < 320) return 4;
    else if (id < 384) return 5;
 
    else if (id < 464) return 0;
    else if (id < 544) return 1;
    else if (id < 624) return 2;
 
    else if (id < 720) return 0;
    else if (id < 816) return 1;
    else if (id < 912) return 2;
    else if (id < 1008) return 3;
    else if (id < 1104) return 4;
    else if (id < 1200) return 5;
 
    else if (id < 1328) return 0;
    else if (id < 1456) return 1;
    else if (id < 1584) return 2;
 
    else if (id < 1728) return 0;
    else if (id < 1872) return 1;
    else if (id < 2016) return 2;
    else if (id < 2160) return 3;
    else if (id < 2304) return 4;
 
    else if (id < 2464) return 0;
    else if (id < 2624) return 1;
    else if (id < 2784) return 2;
    else if (id < 2944) return 3;
 
    else if (id < 3136) return 0;
    else if (id < 3328) return 1;
    else if (id < 3520) return 2;
    else if (id < 3712) return 3;
    else if (id < 3904) return 4;
 
    else if (id < 4112) return 0;
    else if (id < 4320) return 1;
    else if (id < 4528) return 2;
    else if (id < 4736) return 3;
 
    else if (id < 4976) return 0;
    else if (id < 5216) return 1;
    else if (id < 5456) return 2;
    else if (id < 5696) return 3;
    else if (id < 5936) return 4;
 
    else if (id < 6192) return 0;
    else if (id < 6448) return 1;
    else if (id < 6704) return 2;
    else if (id < 6960) return 3;
 
    else if (id < 7248) return 0;
    else if (id < 7536) return 1;
    else if (id < 7824) return 2;
    else if (id < 8112) return 3;
    else if (id < 8400) return 4;
 
    return 9999;    
*/
}

Referenced by TMDCLayer::TMDCLayer().

name()

std::string TMDC::name

(

void

)

const

returns name.

Definition at line 76 of file TMDC.cxx.

                     {
    return "TMDC";
}

Referenced by dump().

nLayers()

unsigned TMDC::nLayers ( void ) const

inline

Definition at line 322 of file TMDC.h.

                        {
    return _nLayers;
}

Referenced by TMDCWire::neighbor().

nLocalLayer()

unsigned TMDC::nLocalLayer ( unsigned superLayerId ) const

inline

returns # of wire layers in a super layer. 0 will be returned if 'superLayerId' is invalid.

Definition at line 251 of file TMDC.h.

                                             {
    AList<TMDCLayer> * superLayer = _superLayers[superLayerId];
    if (! superLayer) return 0;
    return superLayer->length();
}

nSuperLayers()

unsigned TMDC::nSuperLayers ( void ) const

inline

Definition at line 316 of file TMDC.h.

                             {
    return _nSuperLayers;
}

nWires()

unsigned TMDC::nWires ( void ) const

inline

Definition at line 310 of file TMDC.h.

                       {
    return _nWires;
}

phi()

float TMDC::phi ( float phi )

inlinestatic

Definition at line 283 of file TMDC.h.

                 {
    static const float pi2 = 2. * M_PI;
    a = std::fmod(a, pi2);
    if (a < 0.) return a + pi2;
    return a;
}

Referenced by wire(), and wire().

stereoHits()

const AList< TMDCWireHit > & TMDC::stereoHits

(

unsigned

mask = 0

)

const

returns a list of stereo hits. 'update()' must be called before calling this function.

Definition at line 526 of file TMDC.cxx.

                                    {
    if (! mask) return _stereoHits;
    else if (mask == WireHitFindingValid) return _stereoHits;
    std::cout << "TMDC::stereoHits !!! unsupported mask given" << std::endl;
    return _stereoHits;
}

Referenced by TrkReco::execute().

superLayer()

const AList< TMDCLayer > *const TMDC::superLayer ( unsigned id ) const

inline

returns a pointer to a super-layer. 0 will be returned if 'id' is invalid.

Definition at line 271 of file TMDC.h.

                                  {
    return _superLayers[id];
}

Referenced by TBuilder0::buildStereo(), classification(), and nLocalLayer().

superLayerId() [1/3]

unsigned TMDC::superLayerId ( const MdcGeoLayer * const l )

static

Definition at line 810 of file TMDC.cxx.

                                              {
    if (! l) return 9999;
    unsigned id = l->Id();
 
    if(id < 44) return id/4;
    return 9999;     
/*    if (id < 6) return 0;
    else if (id < 9) return 1;
    else if (id < 15) return 2;
    else if (id < 18) return 3;
    else if (id < 23) return 4;
    else if (id < 27) return 5;
    else if (id < 32) return 6;
    else if (id < 36) return 7;
    else if (id < 41) return 8;
    else if (id < 45) return 9;
    else if (id < 50) return 10;
 
    return 9999;
*/
}

superLayerId() [2/3]

unsigned TMDC::superLayerId ( const MdcGeoWire * const w )

static

Definition at line 796 of file TMDC.cxx.

                                             {
    if (! w) return 9999;
    unsigned id = w->Id();
    return superLayerId(id);
}

superLayerId() [3/3]

unsigned TMDC::superLayerId ( unsigned wireId )

static

Definition at line 765 of file TMDC.cxx.

                              {
    if (id < 188) return 0;
    else if (id < 484) return 1;
    else if (id < 812) return 2;
    else if (id < 1236) return 3;
    else if (id < 1772) return 4;
    else if (id < 2412) return 5;
    else if (id < 3180) return 6;
    else if (id < 4012) return 7;
    else if (id < 4972) return 8;
    else if (id < 5996) return 9;
    else if (id < 6860) return 10;
 
    return 9999;    
/*    if (id < 384) return 0;
    else if (id < 624) return 1;
    else if (id < 1200) return 2;
    else if (id < 1584) return 3;
    else if (id < 2304) return 4;
    else if (id < 2944) return 5;
    else if (id < 3904) return 6;
    else if (id < 4736) return 7;
    else if (id < 5936) return 8;
    else if (id < 6960) return 9;
    else if (id < 8400) return 10;
 
    return 9999;
*/
}

Referenced by nLocalLayer(), superLayerId(), TMDCLayer::TMDCLayer(), and wireName().

update()

void TMDC::update

(

bool

mcAnalysis = true

)

updates TMDC information. clear() is called in this function.

Definition at line 317 of file TMDC.cxx.

                            {
    //...Already updated?...
//zsl    if (updated()) return;
 
    //...Clear old information...
    fastClear();
 
    //...Loop over RECMDC_WIRHIT bank...
//    unsigned nReccdc = BsCouTab(RECMDC_WIRHIT);
     unsigned nReccdc = MdcRecWirhitCol::getMdcRecWirhitCol()->size();
#ifdef TRKRECO_DEBUG
    cout<<"size of MdcRecWirhit:"<<nReccdc<<endl;   
#endif
    for (unsigned i = 0; i < nReccdc; i++) {
//  struct reccdc_wirhit * h =
//      (struct reccdc_wirhit *)
//      BsGetEnt(RECMDC_WIRHIT, i + 1, BBS_No_Index);
    MdcRec_wirhit* h = &(*MdcRecWirhitCol::getMdcRecWirhitCol())[i]; 
    //...Check validity...
//zsl   if (! (h->stat & WireHitFindingValid)) continue;
 
    //...Obtain a pointer to GEOMDC...
//  struct geocdc_wire * g =
//      (struct geocdc_wire *)
//      BsGetEnt(GEOMDC_WIRE, h->m_geo, BBS_No_Index);
    const MdcGeoWire* g = h->geo;
//  if (_newCdc)
//      if (g->m_ID < 256)
//      continue;
 
    //...Get a pointer to a TMDCWire...
//  TMDCWire * w = _wires[g->m_ID - 1];
        TMDCWire * w = _wires[g->Id()];   
    _hitWires.append(w);
 
    //...Create TMDCWireHit...
    TMDCWireHit * hit = new TMDCWireHit(w, h, _fudgeFactor);
    _hits.append(hit);
 
    //... update the wirehit pointer in the TMDCWire  (by zang shilei)
        w->hit(hit);
       //if (w->hit()) cout<<"test .......... right~~~~~"<<endl; 
    
    //...Axial or stereo...
    if (w->axial()) _axialHits.append(hit);
    else            _stereoHits.append(hit);
    }
        
    //...Hit classfication...
    classification();
 
    //...MC information...
    if (mcAnalysis) updateMC();
 
    //...Update information...
    TUpdater::update();
}

Referenced by TrkReco::execute().

updateMC()

void TMDC::updateMC

(

void

)

updates TMDC information for MC.

Definition at line 376 of file TMDC.cxx.

                   {
 
    //...Create TTrackHEP...
    TTrackHEP::update();
 
    //...Loop over DATMDC_MCWIRHIT bank...
    unsigned n = 0;
//    for (unsigned i = 0; i < BsCouTab(DATMDC_MCWIRHIT); i++) {
    for (unsigned i = 0; i < MdcDatMcwirhitCol::getMdcDatMcwirhitCol()->size(); i++) {   
//        struct datcdc_mcwirhit * h =
//      (struct datcdc_mcwirhit *)
//      BsGetEnt(DATMDC_MCWIRHIT, i + 1, BBS_No_Index);
    MdcDat_mcwirhit* h = &(*MdcDatMcwirhitCol::getMdcDatMcwirhitCol())[i];
    //...Get a pointer to RECMDC_WIRHIT...
//  reccdc_wirhit * whp =
//      (reccdc_wirhit *) BsGetEnt(RECMDC_WIRHIT, h->m_dat, BBS_No_Index);
    MdcRec_wirhit* whp = h->dat->rec;
    //...Get TrkReco objects...
    TMDCWireHit * wh = 0;
    TMDCWire * w = 0;
    if (whp) {
        if (whp->stat & WireHitFindingValid) {
        unsigned n = _hits.length();
        unsigned j = (whp->id < n) ? whp->id : n;
        while (j) {
            --j;
            if (_hits[j]->reccdc() == whp) {
            wh = _hits[j];
            w = _wires[wh->wire()->id()];
            break;
            }
        }
        }
    }
    if (! w) {
//      geocdc_wire * g =
//      (geocdc_wire *) BsGetEnt(GEOMDC_WIRE, h->m_geo, BBS_No_Index);
//      w = _wires[g->m_ID - 1];
      const MdcGeoWire* g = h->geo;
    }
 
    //...Create TMDCWireHitMC...
    TMDCWireHitMC * hit = new TMDCWireHitMC(w, wh, h);
    _hitsMC.append(hit);
    w->hit(hit);
    if (wh) wh->mc(hit);
 
    //...TTrackHEP...
//  TTrackHEP * hep = TTrackHEP::list()[h->m_hep - 1];
        TTrackHEP * hep = TTrackHEP::list()[h->hep->id];     
    hit->_hep = hep;
    if (hep) hep->_hits.append(hit);
    else {
        std::cout << "TMDC::updateMC !!! mission impossible" << std::endl;
        std::cout << "                   This error will cause TrkReco crush";
        std::cout << std::endl;
#ifdef TRKRECO_DEBUG_DETAIL
        //cout << "    h->m_hep, h->m_hep -1 = " << h->m_hep;
        //std::cout << ", " << h->m_hep - 1 << std::endl;
        std::cout << "    TTrackHEP list length = ";
        std::cout << TTrackHEP::list().length() << std::endl;
//      BsShwDat(GEN_HEPEVT);
//      BsShwDat(DATMDC_MCWIRHIT);
#endif
    }
    }
}

Referenced by update().

version()

std::string TMDC::version

(

void

)

const

returns version.

Definition at line 81 of file TMDC.cxx.

                        {
    return "2.18";
}

Referenced by dump(), and getTMDC().

wire() [1/4]

const TMDCWire *const TMDC::wire

(

const HepPoint3D &

point

)

const

returns a pointer to a wire.

Definition at line 250 of file TMDC.cxx.

                                     {
    float r = p.mag();
    float phi = p.phi();
    return wire(r, phi);
}

wire() [2/4]

const TMDCWire *const TMDC::wire	(	float	r,
		float	phi ) const

returns a pointer to a wire.

Definition at line 257 of file TMDC.cxx.

                                 {
 
    std::cout << "r,phi = " << r << "," << p << std::endl;
 
    unsigned id = 25;
    bool ok = false;
    const TMDCLayer * l;
    while (! ok) {
    l = layer(id);
    if (! l) return NULL;
    
    const MdcGeoLayer * geo = l->geocdc();
    if (geo->Radius()/10 + geo->RCSiz2()/10 < r) ++id;
    else if (geo->Radius()/10 - geo->RCSiz1()/10 > r) --id;
    else ok = true;
    }
//    float dPhi = 2. * M_PI / float(l->nWires());
//    if (l->geocdc()->m_offset > 0.) p -= dPhi / 2.;
    float dPhi = 2. * M_PI / float(l->nWires());   
    p -= l->geocdc()->Offset()/dPhi; 
    unsigned localId = unsigned(phi(p) / dPhi);
    return l->wire(localId);
}

wire() [3/4]

const TMDCWire *const TMDC::wire	(	unsigned	layerId,
		int	localId ) const

returns a pointer to a wire. 'localId' can be negative. 0 will be returned if 'layerId' is invalid.

Definition at line 242 of file TMDC.cxx.

                                              {
    if (layerId < _nLayers)
    return _layers[layerId]->wire(localId);
 
    return NULL;
}

wire() [4/4]

const TMDCWire *const TMDC::wire ( unsigned wireId ) const

inline

returns a pointer to a wire. 0 will be returned if 'wireId' is invalid.

Definition at line 277 of file TMDC.h.

                            {
    return _wires[id];
}

Referenced by TConformalFinder::conformal2Wire(), driftDistance(), TMDCWire::neighbor(), and wire().

wireName() [1/3]

std::string TMDC::wireName ( const MdcGeoWire * const w )

static

Definition at line 583 of file TMDC.cxx.

                                         {
    if (! w) return std::string("no such a wire");
//    unsigned id = w->m_ID - 1;
    unsigned id = w->Id();
    return wireName(id);
}

wireName() [2/3]

std::string TMDC::wireName ( const MdcRec_wirhit & h )

static

Definition at line 1020 of file TMDC.cxx.

                                      {
//    geocdc_wire * g = (geocdc_wire *) BsGetEnt(GEOMDC_WIRE,
//                         h.m_geo,
//                         BBS_No_Index);
      const MdcGeoWire* g = h.geo;
      return wireName(g);
}

wireName() [3/3]

std::string TMDC::wireName ( unsigned wireId )

static

Definition at line 575 of file TMDC.cxx.

                              {
    if (superLayerId(wireId) % 2)
    return itostring(layerId(wireId)) + "=" +
        itostring(localId(wireId));
    return itostring(layerId(wireId)) + "-" + itostring(localId(wireId));
}

Referenced by wireName(), and wireName().

---

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

• TMDC.h
• TMDC.cxx