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
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
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
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.
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
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
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
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.
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
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
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
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.)
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.)
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.)

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 InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

Member Function Documentation

axialHits() [1/3]

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 482 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;
}

axialHits() [2/3]

const AList< TMDCWireHit > & TMDC::axialHits

(

unsigned

mask = 0

)

const

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

axialHits() [3/3]

const AList< TMDCWireHit > & TMDC::axialHits

(

unsigned

mask = 0

)

const

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

axialStereoLayerId() [1/3]

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

static

Definition at line 939 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 layerId(), and TMDCLayer::TMDCLayer().

axialStereoLayerId() [2/3]

unsigned TMDC::axialStereoLayerId ( const MdcGeoLayer * const )

static

axialStereoLayerId() [3/3]

unsigned TMDC::axialStereoLayerId ( const MdcGeoLayer * const )

static

badHits() [1/3]

const AList< TMDCWireHit > & TMDC::badHits

(

void

)

returns bad hits(finding invalid hits).

Definition at line 503 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;
}

badHits() [2/3]

const AList< TMDCWireHit > & TMDC::badHits

(

void

)

returns bad hits(finding invalid hits).

badHits() [3/3]

const AList< TMDCWireHit > & TMDC::badHits

(

void

)

returns bad hits(finding invalid hits).

cdcVersion() [1/3]

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

inline

returns MDC version.

Definition at line 269 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return std::string( _cdcVersion ); }

Referenced by dump(), getTMDC(), and getTMDC().

cdcVersion() [2/3]

std::string TMDC::cdcVersion

(

void

)

const

returns MDC version.

cdcVersion() [3/3]

std::string TMDC::cdcVersion

(

void

)

const

returns MDC version.

classification() [1/3]

void TMDC::classification

(

void

)

classify hits.

Definition at line 411 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().

classification() [2/3]

void TMDC::classification

(

void

)

classify hits.

classification() [3/3]

void TMDC::classification

(

void

)

classify hits.

clear() [1/3]

void TMDC::clear ( void )

virtual

clears all TMDC information.

Reimplemented from TUpdater.

Definition at line 247 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();
}

clear() [2/3]

void TMDC::clear ( void )

virtual

clears all TMDC information.

Reimplemented from TUpdater.

clear() [3/3]

void TMDC::clear ( void )

virtual

clears all TMDC information.

Reimplemented from TUpdater.

debugLevel() [1/6]

int TMDC::debugLevel ( int a )

inline

sets debug level.

Definition at line 240 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _debugLevel = a; }

debugLevel() [2/6]

int TMDC::debugLevel

(

int

)

sets debug level.

debugLevel() [3/6]

int TMDC::debugLevel

(

int

)

sets debug level.

debugLevel() [4/6]

int TMDC::debugLevel ( void ) const

inline

returns debug level.

Definition at line 238 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _debugLevel; }

debugLevel() [5/6]

int TMDC::debugLevel

(

void

)

const

returns debug level.

debugLevel() [6/6]

int TMDC::debugLevel

(

void

)

const

returns debug level.

driftDistance() [1/3]

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 975 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()};
  double            p[3]   = { tp.x(), tp.y(), tp.z() };
  const HepPoint3D& onWire = l.positionOnWire();
  // maqm float x[3] = {onWire.x(), onWire.y(), onWire.z()};
  double x[3] = { onWire.x(), onWire.y(), onWire.z() };
  float  time = h.reccdc()->tdc + t0Offset - tof;
  float  dist = 1.e6;// add initialize 25-05-15
  float  edist = 1.e9;// add initialize 25-05-15
  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 );
  }
}

driftDistance() [2/3]

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)

driftDistance() [3/3]

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)

dump() [1/3]

void TMDC::dump

(

const std::string &

message

)

const

dumps debug information.

Definition at line 157 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().

dump() [2/3]

void TMDC::dump

(

const std::string &

message

)

const

dumps debug information.

dump() [3/3]

void TMDC::dump

(

const std::string &

message

)

const

dumps debug information.

fastClear() [1/3]

void TMDC::fastClear

(

void

)

clears TMDC information.

Definition at line 261 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().

fastClear() [2/3]

void TMDC::fastClear

(

void

)

clears TMDC information.

fastClear() [3/3]

void TMDC::fastClear

(

void

)

clears TMDC information.

fudgeFactor() [1/6]

float TMDC::fudgeFactor ( float a )

inline

sets/returns fudge factor for drift time error.

Definition at line 267 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _fudgeFactor = a; }

fudgeFactor() [2/6]

float TMDC::fudgeFactor

(

float

)

sets/returns fudge factor for drift time error.

fudgeFactor() [3/6]

float TMDC::fudgeFactor

(

float

)

sets/returns fudge factor for drift time error.

fudgeFactor() [4/6]

float TMDC::fudgeFactor ( void ) const

inline

returns fudge factor for drift time error.

Definition at line 265 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _fudgeFactor; }

fudgeFactor() [5/6]

float TMDC::fudgeFactor

(

void

)

const

returns fudge factor for drift time error.

fudgeFactor() [6/6]

float TMDC::fudgeFactor

(

void

)

const

returns fudge factor for drift time error.

getTMDC() [1/6]

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

static

returns TMDC object.

Definition at line 79 of file TMDC.cxx.

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

getTMDC() [2/6]

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

static

returns TMDC object.

getTMDC() [3/6]

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

static

returns TMDC object.

getTMDC() [4/6]

TMDC * TMDC::getTMDC ( void )

static

Definition at line 84 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().

getTMDC() [5/6]

TMDC * TMDC::getTMDC ( void )

static

getTMDC() [6/6]

TMDC * TMDC::getTMDC ( void )

static

hits() [1/3]

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 496 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;
}

hits() [2/3]

const AList< TMDCWireHit > & TMDC::hits

(

unsigned

mask = 0

)

const

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

hits() [3/3]

const AList< TMDCWireHit > & TMDC::hits

(

unsigned

mask = 0

)

const

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

hitsMC() [1/3]

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

inline

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

Definition at line 250 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _hitsMC; }

Referenced by TPerfectFinder::doit().

hitsMC() [2/3]

const AList< TMDCWireHitMC > & TMDC::hitsMC

(

void

)

const

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

hitsMC() [3/3]

const AList< TMDCWireHitMC > & TMDC::hitsMC

(

void

)

const

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

layer() [1/3]

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 248 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _layers[id]; }

Referenced by wire().

layer() [2/3]

const TMDCLayer *const TMDC::layer

(

unsigned

id

)

const

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

layer() [3/3]

const TMDCLayer *const TMDC::layer

(

unsigned

id

)

const

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

layerId() [1/9]

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

static

Definition at line 618 of file TMDC.cxx.

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

layerId() [2/9]

unsigned TMDC::layerId ( const MdcGeoWire * const )

static

layerId() [3/9]

unsigned TMDC::layerId ( const MdcGeoWire * const )

static

layerId() [4/9]

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

static

Definition at line 951 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() [5/9]

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

static

layerId() [6/9]

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

static

layerId() [7/9]

unsigned TMDC::layerId ( unsigned wireId )

static

Definition at line 556 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().

layerId() [8/9]

unsigned TMDC::layerId ( unsigned wireId )

static

layerId() [9/9]

unsigned TMDC::layerId ( unsigned wireId )

static

localId() [1/6]

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

static

Definition at line 705 of file TMDC.cxx.

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

localId() [2/6]

unsigned TMDC::localId ( const MdcGeoWire * const )

static

localId() [3/6]

unsigned TMDC::localId ( const MdcGeoWire * const )

static

localId() [4/6]

unsigned TMDC::localId ( unsigned wireId )

static

Definition at line 632 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().

localId() [5/6]

unsigned TMDC::localId ( unsigned wireId )

static

localId() [6/6]

unsigned TMDC::localId ( unsigned wireId )

static

localLayerId() [1/9]

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

static

Definition at line 868 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/9]

unsigned TMDC::localLayerId ( const MdcGeoLayer * const )

static

localLayerId() [3/9]

unsigned TMDC::localLayerId ( const MdcGeoLayer * const )

static

localLayerId() [4/9]

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

static

Definition at line 853 of file TMDC.cxx.

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

localLayerId() [5/9]

unsigned TMDC::localLayerId ( const MdcGeoWire * const )

static

localLayerId() [6/9]

unsigned TMDC::localLayerId ( const MdcGeoWire * const )

static

localLayerId() [7/9]

unsigned TMDC::localLayerId ( unsigned wireId )

static

Definition at line 784 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().

localLayerId() [8/9]

unsigned TMDC::localLayerId ( unsigned wireId )

static

localLayerId() [9/9]

unsigned TMDC::localLayerId ( unsigned wireId )

static

name() [1/3]

std::string TMDC::name

(

void

)

const

returns name.

Definition at line 73 of file TMDC.cxx.

{ return "TMDC"; }

Referenced by dump().

name() [2/3]

std::string TMDC::name

(

void

)

const

returns name.

name() [3/3]

std::string TMDC::name

(

void

)

const

returns name.

nLayers() [1/3]

unsigned TMDC::nLayers ( void ) const

inline

Definition at line 275 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _nLayers; }

nLayers() [2/3]

unsigned TMDC::nLayers

(

void

)

const

nLayers() [3/3]

unsigned TMDC::nLayers

(

void

)

const

nLocalLayer() [1/3]

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 242 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

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

nLocalLayer() [2/3]

unsigned TMDC::nLocalLayer

(

unsigned

superLayerId

)

const

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

nLocalLayer() [3/3]

unsigned TMDC::nLocalLayer

(

unsigned

superLayerId

)

const

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

nSuperLayers() [1/3]

unsigned TMDC::nSuperLayers ( void ) const

inline

Definition at line 273 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _nSuperLayers; }

nSuperLayers() [2/3]

unsigned TMDC::nSuperLayers

(

void

)

const

nSuperLayers() [3/3]

unsigned TMDC::nSuperLayers

(

void

)

const

nWires() [1/3]

unsigned TMDC::nWires ( void ) const

inline

Definition at line 271 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _nWires; }

nWires() [2/3]

unsigned TMDC::nWires

(

void

)

const

nWires() [3/3]

unsigned TMDC::nWires

(

void

)

const

phi() [1/3]

float TMDC::phi ( float phi )

inlinestatic

Definition at line 258 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

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

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

phi() [2/3]

float TMDC::phi ( float phi )

static

phi() [3/3]

float TMDC::phi ( float phi )

static

stereoHits() [1/3]

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 489 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;
}

stereoHits() [2/3]

const AList< TMDCWireHit > & TMDC::stereoHits

(

unsigned

mask = 0

)

const

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

stereoHits() [3/3]

const AList< TMDCWireHit > & TMDC::stereoHits

(

unsigned

mask = 0

)

const

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

superLayer() [1/3]

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 252 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

                                                                         {
  return _superLayers[id];
}

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

superLayer() [2/3]

const AList< TMDCLayer > *const TMDC::superLayer

(

unsigned

id

)

const

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

superLayer() [3/3]

const AList< TMDCLayer > *const TMDC::superLayer

(

unsigned

id

)

const

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

superLayerId() [1/9]

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

static

Definition at line 762 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/9]

unsigned TMDC::superLayerId ( const MdcGeoLayer * const )

static

superLayerId() [3/9]

unsigned TMDC::superLayerId ( const MdcGeoLayer * const )

static

superLayerId() [4/9]

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

static

Definition at line 749 of file TMDC.cxx.

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

superLayerId() [5/9]

unsigned TMDC::superLayerId ( const MdcGeoWire * const )

static

superLayerId() [6/9]

unsigned TMDC::superLayerId ( const MdcGeoWire * const )

static

superLayerId() [7/9]

unsigned TMDC::superLayerId ( unsigned wireId )

static

Definition at line 719 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(), TMDCWire::superLayerId(), TMDCLayer::TMDCLayer(), and wireName().

superLayerId() [8/9]

unsigned TMDC::superLayerId ( unsigned wireId )

static

superLayerId() [9/9]

unsigned TMDC::superLayerId ( unsigned wireId )

static

update() [1/3]

void TMDC::update

(

bool

mcAnalysis = true

)

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

Definition at line 277 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();
}

update() [2/3]

void TMDC::update

(

bool

mcAnalysis = true

)

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

update() [3/3]

void TMDC::update

(

bool

mcAnalysis = true

)

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

updateMC() [1/3]

void TMDC::updateMC

(

void

)

updates TMDC information for MC.

Definition at line 336 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().

updateMC() [2/3]

void TMDC::updateMC

(

void

)

updates TMDC information for MC.

updateMC() [3/3]

void TMDC::updateMC

(

void

)

updates TMDC information for MC.

version() [1/3]

std::string TMDC::version

(

void

)

const

returns version.

Definition at line 75 of file TMDC.cxx.

{ return "2.18"; }

Referenced by dump(), and getTMDC().

version() [2/3]

std::string TMDC::version

(

void

)

const

returns version.

version() [3/3]

std::string TMDC::version

(

void

)

const

returns version.

wire() [1/12]

const TMDCWire *const TMDC::wire

(

const HepPoint3D &

point

)

const

returns a pointer to a wire.

Definition at line 216 of file TMDC.cxx.

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

wire() [2/12]

const TMDCWire *const TMDC::wire

(

const HepPoint3D &

point

)

const

returns a pointer to a wire.

wire() [3/12]

const TMDCWire *const TMDC::wire

(

const HepPoint3D &

point

)

const

returns a pointer to a wire.

wire() [4/12]

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

returns a pointer to a wire.

Definition at line 222 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() [5/12]

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

returns a pointer to a wire.

wire() [6/12]

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

returns a pointer to a wire.

wire() [7/12]

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 210 of file TMDC.cxx.

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

wire() [8/12]

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.

wire() [9/12]

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.

wire() [10/12]

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 256 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h.

{ return _wires[id]; }

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

wire() [11/12]

const TMDCWire *const TMDC::wire

(

unsigned

wireId

)

const

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

wire() [12/12]

const TMDCWire *const TMDC::wire

(

unsigned

wireId

)

const

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

wireName() [1/9]

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

static

Definition at line 542 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/9]

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

static

wireName() [3/9]

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

static

wireName() [4/9]

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

static

Definition at line 967 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() [5/9]

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

static

wireName() [6/9]

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

static

wireName() [7/9]

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

static

Definition at line 536 of file TMDC.cxx.

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

Referenced by wireName(), and wireName().

wireName() [8/9]

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

static

wireName() [9/9]

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

static

---

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

• InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TMDC.h
• InstallArea/x86_64-el9-gcc13-opt/include/TrkReco/TMDC.h
• Reconstruction/TrkReco/include/TrkReco/TMDC.h
• TMDC.cxx