#include <TofTrack.h>

Public Member Functions
	TofTrack ()

	~TofTrack ()

int	trackId () const

int	tofTrackId () const

int	id1 () const

int	id2 () const

unsigned int	barrel () const

ExtTrackCase	hitCase () const

double	p () const

double	path () const

double	path1 () const

double	path2 () const

double	zrhit1 () const

double	zrhit2 () const

double	errzr1 () const

double	errzr2 () const

int	kal (unsigned int i) const

double	zrhit1 (unsigned int i) const

double	zrhit2 (unsigned int i) const

double	errzr1 (unsigned int i) const

double	errzr2 (unsigned int i) const

std::vector< TofData * >	tofData1 () const

std::vector< TofData * >	tofData2 () const

int	size1 () const

int	size2 () const

int	size3 () const

int	size4 () const

double	theta1 () const

double	theta2 () const

double	phi1 () const

double	phi2 () const

unsigned int	quality1 () const

unsigned int	quality2 () const

unsigned int	quality () const

int	tofId1 () const

int	tofId2 () const

double	ph11 () const

double	ph12 () const

double	ph21 () const

double	ph22 () const

double	ph1 () const

double	ph2 () const

double	ph () const

double	tof11 (unsigned int i) const

double	tof12 (unsigned int i) const

double	tof21 (unsigned int i) const

double	tof22 (unsigned int i) const

double	tof1 (unsigned int i) const

double	tof2 (unsigned int i) const

double	tof (unsigned int i) const

double	qch1 () const

double	qch2 () const

double	qch3 () const

double	qch4 () const

double	adc1 () const

double	adc2 () const

double	adc3 () const

double	adc4 () const

double	tdc1 () const

double	tdc2 () const

double	tdc3 () const

double	tdc4 () const

double	texpInner (unsigned int i) const

double	texpOuter (unsigned int i) const

double	texp (unsigned int i) const

double	ztdc1 () const

double	ztdc2 () const

double	zadc1 () const

double	zadc2 () const

double	estime () const

double	tdiff1 () const

double	tdiff2 () const

int	t0Stat () const

unsigned int	flag () const

bool	isNoHit () const

void	setQuality (int qual)

void	setQuality1 (int qual1)

void	setQuality2 (int qual2)

void	setFlag (unsigned int flag)

void	setExtTrack (RecExtTrack *extTrack, int kal[5], double t0, int t0Stat)

void	getMultiHit (TofTrack *&)

void	setTofData (TofDataMap tofDataMap)

void	tofDataAnalysis (TofData *tof, unsigned int iflag)

void	getTofData (TofData *tof, unsigned int iflag)

void	getTofDataEast (TofData *tof, unsigned int iflag)

void	getTofDataWest (TofData *tof, unsigned int iflag)

void	getTofDataNohit (unsigned int iflag)

void	match (bool forCalibration, std::vector< int > deadId, std::vector< TofTrack * > *&tofTrackVec)

void	findTofDataBarrel (std::vector< TofData * > tofDataVec1, std::vector< TofData * > tofDataVec2, double zrhit, unsigned int iflag, std::vector< TofTrack * > *&tofTrackVec)

TofData *	chooseTofData (std::vector< TofData * > tofDataVec, double zrhit)

TofData *	compareTofData (TofData tofData1, TofData tofData2, double zrhit)

void	findTofDataEndcap (std::vector< TofData * > tofDataVec1, std::vector< TofData * > tofDataVec2, double zr1[5])

TofData *	chooseTofDataEndcap (std::vector< TofData * > tofDataVec, double zr1[5])

TofData *	compareTofDataEndcap (TofData tofData1, TofData tofData2)

void	tofDataStudy ()

void	setCalibration ()

void	convert2RecTofTrackCol (RecTofTrackCol *recTofTrackCol)

void	buildRecTofTrack (RecTofTrack *track, int layerorend)

void	setRecTofTrack (RecTofTrack *track, int layerorend)

void	convert2RecBTofCalHitCol (int runNumber, int eventNumber, RecBTofCalHitCol *btofCalHitCol, std::string calibData)

void	convert2RecETofCalHitCol (int runNumber, int eventNumber, RecETofCalHitCol *etofCalHitCol, std::string calibData)

void	qualityAnalysis ()

Detailed Description

Definition at line 20 of file TofTrack.h.

Constructor & Destructor Documentation

◆ TofTrack()

TofTrack::TofTrack ( )

Definition at line 14 of file TofTrack.cxx.

                   {
  m_trackId = -1;
  m_tofTrackId = -1;
  m_id1 = -99;
  m_id2 = -99;
  m_hitCase = NoHit;
  m_momentum = -99.0;
  m_path  = 0.0;
  m_path1 = 0.0;
  m_path2 = 0.0;
  m_zrhit1 = 9999.0;
  m_errzr1 = 9999.0;
  m_zrhit2 = 9999.0;
  m_errzr2 = 9999.0;
  for( unsigned int i=0; i<5; i++ ) {
    m_kal[i]  = -1;
    m_zr1[i]  = 9999.0;
    m_zr2[i]  = 9999.0;
    m_ezr1[i] = 9999.0;
    m_ezr2[i] = 9999.0;
  }
  if( m_tofData1.size()>0 ) {
    std::vector<TofData*>::iterator iter1 = m_tofData1.begin();
    m_tofData1.clear();
  }
  if( m_tofData2.size()>0 ) {
    std::vector<TofData*>::iterator iter2 = m_tofData2.begin();
    m_tofData2.clear();
  }
  if( m_tofData3.size()>0 ) {
    std::vector<TofData*>::iterator iter3 = m_tofData3.begin();
    m_tofData3.clear();
  }
  if( m_tofData4.size()>0 ) {
    std::vector<TofData*>::iterator iter4 = m_tofData4.begin();
    m_tofData4.clear();
  }
  m_quality1 = 0;
  m_quality2 = 0;
  m_quality  = 10;
 
  m_tofId1 = -99;
  m_tofId2 = -99;
 
  m_ph11  = -999.0;
  m_ph12  = -999.0;
  m_ph21  = -999.0;
  m_ph22  = -999.0;
  m_ph1   = -999.0;
  m_ph2   = -999.0;
  m_ph    = -999.0;
 
  for( unsigned int i=0; i<5; i++ ) {
    m_tof11[i]  = 0.0;
    m_tof12[i]  = 0.0;
    m_tof21[i]  = 0.0;
    m_tof22[i]  = 0.0;
    m_tof1[i]   = 0.0;
    m_tof2[i]   = 0.0;
    m_tof[i]    = 0.0;
  }
 
  m_sigma11 = -999.0;
  m_sigma12 = -999.0;
  m_sigma21 = -999.0;
  m_sigma22 = -999.0;
  m_sigma1  = -999.0;
  m_sigma2  = -999.0;
  m_sigma   = -999.0;
 
  m_qch1  = -999.0;
  m_qch2  = -999.0;
  m_qch3  = -999.0;
  m_qch4  = -999.0;
  m_adc1  = -999.0;
  m_adc2  = -999.0;
  m_adc3  = -999.0;
  m_adc4  = -999.0;
  m_tdc1  = -999.0;
  m_tdc2  = -999.0;
  m_tdc3  = -999.0;
  m_tdc4  = -999.0;
 
  for( unsigned int i=0; i<5; i++ ) {
    m_texpInner[i] = -999.0;
    m_texpOuter[i] = -999.0;
    m_texp[i] = -999.0;
  }
 
  m_ztdc1 = -999.0;
  m_ztdc2 = -999.0;
  m_zadc1 = -999.0;
  m_zadc2 = -999.0;
 
  m_estime = -999.0;
  m_tdiff1  = -999.0;
  m_tdiff2  = -999.0;
 
  m_flag  = 0;
 
  return;
}

◆ ~TofTrack()

TofTrack::~TofTrack ( )

Definition at line 118 of file TofTrack.cxx.

                    {
  if( m_tofData1.size()>0 ) {
    std::vector<TofData*>::iterator iter1 = m_tofData1.begin();
    m_tofData1.clear();
  }
  if( m_tofData2.size()>0 ) {
    std::vector<TofData*>::iterator iter2 = m_tofData2.begin();
    m_tofData2.clear();
  }
  if( m_tofData3.size()>0 ) {
    std::vector<TofData*>::iterator iter3 = m_tofData3.begin();
    m_tofData3.clear();
  }
  if( m_tofData4.size()>0 ) {
    std::vector<TofData*>::iterator iter4 = m_tofData4.begin();
    m_tofData4.clear();
  }
  return;
}

Member Function Documentation

◆ adc1()

double TofTrack::adc1 ( ) const

inline

Definition at line 81 of file TofTrack.h.

81{ return m_adc1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ adc2()

double TofTrack::adc2 ( ) const

inline

Definition at line 82 of file TofTrack.h.

82{ return m_adc2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ adc3()

double TofTrack::adc3 ( ) const

inline

Definition at line 83 of file TofTrack.h.

83{ return m_adc3; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ adc4()

double TofTrack::adc4 ( ) const

inline

Definition at line 84 of file TofTrack.h.

84{ return m_adc4; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ barrel()

unsigned int TofTrack::barrel ( ) const

inline

Definition at line 29 of file TofTrack.h.

29{ return m_barrel; }

Referenced by convert2RecTofTrackCol(), TofCheckDigi::Fill_TofTrack(), match(), and setCalibration().

◆ buildRecTofTrack()

void TofTrack::buildRecTofTrack	(	RecTofTrack *	track,
		int	layerorend
	)

Definition at line 1558 of file TofTrack.cxx.

                                                                    {
 
  track->setTofTrackID( m_tofTrackId );
  track->setTrackID( m_trackId );
 
  track->setErrTof( 0.0 );
  track->setBeta( 0.0 );
 
  double sigma[6];
  for( int i=0; i<6; i++ ) {
    sigma[i]   = 0.0;
  }
  track->setSigma( sigma );
  track->setQuality( m_quality );
  track->setT0( m_estime );
  track->setErrT0( 0.0 );
  track->setPhi( 9999.0 );
  track->setErrPhi( 9999.0 );
  track->setEnergy( 9999.0 );
  track->setErrEnergy( 9999.0 );
 
  if( ( layerorend == 11 ) || ( layerorend == 12 ) || ( layerorend == 1 ) ) {
    track->setTofID( m_tofId1 );
    track->setPath( m_path1 );
    track->setZrHit( m_zrhit1 );
    track->setErrZ( m_errzr1 );
    track->setTexp( m_texpInner );
 
    setRecTofTrack( track, layerorend );
  }
 
  if( ( layerorend==21 ) || ( layerorend==22 ) || ( layerorend==2 ) ) {
    track->setTofID( m_tofId2 );
    track->setPath( m_path2 );
    track->setZrHit( m_zrhit2 );
    track->setErrZ( m_errzr2 );
    track->setTexp( m_texpOuter );
 
    setRecTofTrack( track, layerorend );
  }
 
  if( layerorend==0 ) {
    track->setTofID( m_tofId1 );
    track->setPath( m_path2 );
    track->setZrHit( m_zrhit2 );
    track->setErrZ( m_errzr2 );
    track->setTexp( m_texp );
 
    setRecTofTrack( track, layerorend );
  }
 
  if( layerorend == 3 ) {
    track->setTofID( m_id1 );
    track->setPath( m_path1 );
    track->setZrHit( m_zrhit1 );
    track->setErrZ( m_errzr1 );
    track->setTexp( m_texpInner );
  }
 
  return;
}

Referenced by convert2RecTofTrackCol().

◆ chooseTofData()

TofData * TofTrack::chooseTofData	(	std::vector< TofData * >	tofDataVec,
		double	zrhit
	)

Definition at line 780 of file TofTrack.cxx.

                                                                               {
  if( tofDataVec.size() == 0 ) {
    cout << "TofRec::TofTrack::ChooseTofData: Size of TofData Vector is Zero!" << endl;
    return 0;
  }
  std::vector<TofData*>::iterator igood = tofDataVec.begin();
  if( tofDataVec.size() > 1 ) {
    double deltaZ = 1000.0;
    std::vector<TofData*>::iterator iter  = tofDataVec.begin();
    // ZTDC compare
    for( ; iter != tofDataVec.end(); iter++ ) {
      if( ( (*iter)->quality() & 0x5 ) == 0x5 ) {
    if( abs( (*iter)->ztdc() - zrhit ) < deltaZ ) {
      deltaZ = abs( (*iter)->ztdc() - zrhit );
      igood = iter;
    }
      }
    }
    // ZADC compare
    if( deltaZ > 999.0 ) {
      iter  = tofDataVec.begin();
      for( ; iter != tofDataVec.end(); iter++ ) {
    if( ( (*iter)->quality() & 0xa ) == 0xa ) {
      if( abs( (*iter)->zadc() - zrhit ) < deltaZ ) {
        deltaZ = abs( (*iter)->zadc() - zrhit );
        igood = iter;
      }
    }
      }
    }
    // Max Q
    if( deltaZ > 999.0 ) {
      unsigned int ibad = 0xf0;
      iter  = tofDataVec.begin();
      for( ; iter != tofDataVec.end(); iter++ ) {
    if( ( (*iter)->quality() & 0xf0 ) < ibad ) {
      igood = iter;
      ibad = ( (*iter)->quality() & 0xf0 );
    }
    else if( ( (*iter)->quality() & 0xf0 ) == ibad ) {
      if( ( (*iter)->adc1() + (*iter)->adc2() ) > ( (*igood)->adc1() + (*igood)->adc2() ) ) {
        igood = iter;
        ibad = ( (*iter)->quality() & 0xf0 );
      }
    }
      }
    }
  }
 
  return (*igood);
}

Referenced by findTofDataBarrel().

◆ chooseTofDataEndcap()

TofData * TofTrack::chooseTofDataEndcap	(	std::vector< TofData * >	tofDataVec,
		double	zr1[5]
	)

Definition at line 940 of file TofTrack.cxx.

                                                                                      {
  if( tofDataVec.size() == 0 ) {
    cout << "TofRec::TofTrack::ChooseTofData: Size of TofData Vector is Zero!" << endl;
    return 0;
  }
  std::vector<TofData*>::iterator igood = tofDataVec.begin();
  if( tofDataVec.size() > 1 ) {
    bool multihit = false;
    std::vector<TofData*>::iterator iter  = tofDataVec.begin();
    for( ; iter != tofDataVec.end(); iter++ ) {
      if( (*iter)->qtimes1()>1 ) { multihit = true; }
    }
    iter = tofDataVec.begin();
    if( multihit ) {
      double tcorr     = -999.0;
      double deltaTMin = 999.0;
      for( ; iter != tofDataVec.end(); iter++ ) {
    tcorr = tofCaliSvc->ETime( (*iter)->adc(), (*iter)->tdc()-m_estime, zr1[2], (*iter)->tofId() );
    for( unsigned int i=0; i<5; i++ ) {
      if( abs(tcorr-m_texpInner[i]) < deltaTMin ) {
        deltaTMin = abs(tcorr-m_texpInner[i]);
        igood     = iter;
      }
    }
      }
    }
    else {
      double maxQ = 0.0;
      for( ; iter != tofDataVec.end(); iter++ ) {
    if( (*iter)->adc() > maxQ ) {
      maxQ  = (*iter)->adc();
      igood = iter;
    }
      }
    }
  }
  return (*igood);
} 

Referenced by findTofDataEndcap().

◆ compareTofData()

TofData * TofTrack::compareTofData	(	TofData *	tofData1,
		TofData *	tofData2,
		double	zrhit
	)

Definition at line 837 of file TofTrack.cxx.

                                                                                      {
  TofData* tof = tofData1;
  // ZTDC compare
  if( abs(tofData1->ztdc() - zrhit ) > abs(tofData2->ztdc() - zrhit ) ) {
    // SingleEnd/NoT/NoQ compare
    if( ( tofData1->quality() & 0xf0 ) >= ( tofData1->quality() & 0xf0 ) ) {
      // QDC compare
      //    if( ( tofData1->adc1() + tofData1->adc2() ) < ( tofData2->adc1() + tofData2->adc2() ) ) {
      tof = tofData2;
      //    }
    }
    //    }
  }
 
  return tof;
}

Referenced by findTofDataBarrel().

◆ compareTofDataEndcap()

TofData * TofTrack::compareTofDataEndcap	(	TofData *	tofData1,
		TofData *	tofData2
	)

Definition at line 983 of file TofTrack.cxx.

                                                                              {
  TofData* tof = tofData1;
  if( tof->adc() < tofData2->adc() ) {
    tof = tofData2;
  }
  return tof;
}

Referenced by findTofDataEndcap().

◆ convert2RecBTofCalHitCol()

void TofTrack::convert2RecBTofCalHitCol	(	int	runNumber,
		int	eventNumber,
		RecBTofCalHitCol *	btofCalHitCol,
		std::string	calibData
	)

Definition at line 1705 of file TofTrack.cxx.

                                                                                                                              {
 
  if( ( m_quality1 & 0x800 ) == 0x800 ) {
 
    RecBTofCalHit* ahit = new RecBTofCalHit;
    ahit->setRun( runNumber );
    ahit->setEvent( eventNumber );
    ahit->setMod( m_tofId1 );
    ahit->setQual( m_hitCase );
 
    for( int i=0; i<5; i++ ) {
      ahit->setTpred( i, m_texpInner[i] );
    }
    if( calibData == "Dimu" ) {
      ahit->setTpred( m_texpInner[1] );
      ahit->setZHit( m_zr1[1] );
      ahit->setdZHit( m_ezr1[1] );
    }
    else {
      ahit->setTpred( m_texpInner[0] );
      ahit->setZHit( m_zr1[0] );
      ahit->setdZHit( m_ezr1[0] );
    }
    
    ahit->setTdc1( m_tdc1-m_estime );
    ahit->setTdc2( m_tdc2-m_estime );
    ahit->setAdc1( m_adc1 );
    ahit->setAdc2( m_adc2 );
    //  ahit->setZHit( m_zrhit1 );
    //  ahit->setdZHit( m_errzr1 );
    ahit->setDeltaPhi( m_estime );
    ahit->setsinTheta( m_theta1 );
    ahit->setP( m_momentum );
    ahit->setQ( m_ph1 );
    ahit->setPath( m_path1 );
 
    btofCalHitCol->push_back( ahit );
 
    if( ( m_quality2 & 0x800 ) == 0x800 ) {
 
      RecBTofCalHit* bhit = new RecBTofCalHit;
      bhit->setRun( runNumber );
      bhit->setEvent( eventNumber );
      bhit->setMod( m_tofId2 );
      bhit->setQual( m_hitCase );
      
      for( int i=0; i<5; i++ ) {
    bhit->setTpred( i, m_texpOuter[i] );
      }
      if( calibData == "Dimu" ) {
    bhit->setTpred( m_texpOuter[1] );
    bhit->setZHit( m_zr2[1] );
    bhit->setdZHit( m_ezr2[1] );
      }
      else {
    bhit->setTpred( m_texpOuter[0] );
    bhit->setZHit( m_zr2[0] );
    bhit->setdZHit( m_ezr2[0] );
      }
 
      bhit->setTdc1( m_tdc3-m_estime );
      bhit->setTdc2( m_tdc4-m_estime );
      bhit->setAdc1( m_adc3 );
      bhit->setAdc2( m_adc4 );
      //  bhit->setZHit( m_zrhit2 );
      //  bhit->setdZHit( m_errzr2 );
      bhit->setDeltaPhi( m_estime );
      bhit->setsinTheta( m_theta2 );
      bhit->setP( m_momentum );
      bhit->setQ( m_ph2 );
      bhit->setPath( m_path2 );
 
      ahit->setnext(bhit);
 
      btofCalHitCol->push_back( bhit );
 
    }
  }
 
  return;
}

◆ convert2RecETofCalHitCol()

void TofTrack::convert2RecETofCalHitCol	(	int	runNumber,
		int	eventNumber,
		RecETofCalHitCol *	etofCalHitCol,
		std::string	calibData
	)

Definition at line 1788 of file TofTrack.cxx.

                                                                                                                              {
 
  if( ( m_quality1 & 0x800 ) != 0x800 ) return;
 
  RecETofCalHit* chit = new RecETofCalHit;
  chit->setRun( runNumber );
  chit->setEvent( eventNumber );
  chit->setMod( m_tofId1 );
  chit->setQual( m_hitCase );
 
  for( int i=0; i<5; i++ ) {
    chit->setTpred( i, m_texpInner[i] );
  }
  if( calibData == "Dimu" ) {
    chit->setTpred( m_texpInner[1] );
    chit->setRHit( m_zr1[1] );
    chit->setdRHit( m_ezr1[1] );
  }
  else {
    chit->setTpred( m_texpInner[0] );
    chit->setRHit( m_zr1[0] );
    chit->setdRHit( m_ezr1[0] );
  }
 
  chit->setTdc( m_tdc1-m_estime );
  chit->setAdc( m_adc1 );
  //  chit->setRHit( m_zrhit1 );
  //  chit->setdRHit( m_errzr1 );
  chit->setDeltaPhi( m_estime );
  chit->setcosTheta( m_theta1 );
  chit->setQ( m_ph1 );
  chit->setP( m_momentum );
  chit->setPath( m_path1 );
 
  etofCalHitCol->push_back( chit );
 
  return;
}

◆ convert2RecTofTrackCol()

void TofTrack::convert2RecTofTrackCol ( RecTofTrackCol * recTofTrackCol )

Definition at line 1344 of file TofTrack.cxx.

                                                                      {
 
  bool barrel = ( ( m_hitCase == InnerLayer ) || ( m_hitCase == OuterLayer ) || ( m_hitCase == DoubleLayer ) );
 
  bool innerEast = ( ( m_quality1 & 0xc ) == 0xc );
  bool innerWest = ( ( m_quality1 & 0x3 ) == 0x3 );
  bool outerEast = ( ( m_quality2 & 0xc ) == 0xc );
  bool outerWest = ( ( m_quality2 & 0x3 ) == 0x3 );
 
  if( barrel ) {
 
    if( innerEast ) {
      RecTofTrack* atrack11 = new RecTofTrack;
      buildRecTofTrack( atrack11, 11 );            // innerlayer east readout
      TofHitStatus* hitStatus11 = new TofHitStatus;
      if( innerWest ) {
    hitStatus11->setBarrelReadout( 1, true );  // innerlayer east readout
      }
      else {
    if( m_hitCase == InnerLayer ) {
      hitStatus11->setBarrelCluster( 11 );     // innerlayer east cluster
    }
    else if( m_hitCase == DoubleLayer ) {
      if( outerEast && outerWest ) {
        hitStatus11->setBarrelCounter( 11 );   // innerlayer east counter
      }
      else {
        hitStatus11->setBarrelCluster( 11 );   // innerlayer east cluster
      }
    }
    else {
      cout << "TofRec::TofTrack:convert2RecTofTrackCol: 11- Impossible!" << endl;
    }
      }
      atrack11->setStatus( hitStatus11->value() );
      delete hitStatus11;
      recTofTrackCol->push_back( atrack11 );
    }
 
    if( innerWest ) {
      RecTofTrack* atrack12 = new RecTofTrack;
      buildRecTofTrack( atrack12, 12 );            // innerlayer west readout
      TofHitStatus* hitStatus12 = new TofHitStatus;
      if( innerEast ) {
    hitStatus12->setBarrelReadout( 1, false ); // innerlayer west
      }
      else {
    if( m_hitCase == InnerLayer ) {
      hitStatus12->setBarrelCluster( 12 );     // innerlayer west cluster
    }
    else if( m_hitCase == DoubleLayer ) {
      if( outerEast && outerWest ) {
        hitStatus12->setBarrelCounter( 12 );   // innerlayer west counter
      }
      else {
        hitStatus12->setBarrelCluster( 12 );   // innerlayer west cluster
      }
    }
    else {
      cout << "TofRec::TofTrack:convert2RecTofTrackCol: 12- Impossible!" << endl;
    }
      }
      atrack12->setStatus( hitStatus12->value() );
      delete hitStatus12;
      recTofTrackCol->push_back( atrack12 );
    }
 
    if( innerEast && innerWest ) {
      RecTofTrack* atrack1  = new RecTofTrack;
      buildRecTofTrack( atrack1, 1 );             // innerlayer counter
      TofHitStatus* hitStatus1 = new TofHitStatus;
      if( m_hitCase == InnerLayer ) {
    hitStatus1->setBarrelCluster( 1 );        // innerlayer cluster and counter
      }
      else if( m_hitCase == DoubleLayer ) {
    if( outerEast && outerWest ) {
      hitStatus1->setBarrelCounter( 1 );      // innerlayer counter
    }
    else {
      hitStatus1->setBarrelCluster( 1 );      // innerlayer cluster and counter
    }
      }
      else {
    cout << "TofRec::TofTrack:convert2RecTofTrackCol: 1- Impossible!" << endl;
      }
      atrack1->setStatus( hitStatus1->value() );
      delete hitStatus1;
      recTofTrackCol->push_back( atrack1 );
    }
 
    if( outerEast ) {
      RecTofTrack* atrack21 = new RecTofTrack;
      buildRecTofTrack( atrack21, 21 );           // outerlayer east readout
      TofHitStatus* hitStatus21 = new TofHitStatus;
      if( outerWest ) {
    hitStatus21->setBarrelReadout( 2, true ); // outerlayer east readout
      }
      else {
    if( m_hitCase == OuterLayer ) {
      hitStatus21->setBarrelCluster( 21 );    // outerlayer east cluster
    }
    else if( m_hitCase == DoubleLayer ) {
      if( innerEast || innerWest ) {
        hitStatus21->setBarrelCounter( 21 );  // outerlayer east counter
      }
      //      else {
      //        hitStatus21->setBarrelCluster( 21 );  // outerlayer east cluster
      //      }
    }
    else {
      cout << "TofRec::TofTrack:convert2RecTofTrackCol: 21- Impossible!" << endl;
    }
      }
      atrack21->setStatus( hitStatus21->value() );
      delete hitStatus21;
      recTofTrackCol->push_back( atrack21 );
    }
 
    if( outerWest ) {
      RecTofTrack* atrack22 = new RecTofTrack;
      buildRecTofTrack( atrack22, 22 );            // outerlayer west readout
      TofHitStatus* hitStatus22 = new TofHitStatus;
      if( outerEast ) {
    hitStatus22->setBarrelReadout( 2, false ); // outerlayer west readout
      }
      else {
    if( m_hitCase == OuterLayer ) {
      hitStatus22->setBarrelCluster( 22 );     // outerlayer west cluster
    }
    else if( m_hitCase == DoubleLayer ) {
      if( innerEast || innerWest ) {
        hitStatus22->setBarrelCounter( 22 );   // outerlayer west counter
      }
      //      else {
      //        hitStatus22->setBarrelCluster( 22 );   // outerlayer west cluster
      //      }
    }
    else {
      cout << "TofRec::TofTrack:convert2RecTofTrackCol: 22- Impossible!" << endl;
    }
      }
      atrack22->setStatus( hitStatus22->value() );
      delete hitStatus22;
      recTofTrackCol->push_back( atrack22 );
    }
 
    if( outerEast && outerWest ) {
      RecTofTrack* atrack2  = new RecTofTrack;
      buildRecTofTrack( atrack2, 2 );             // outerlayer counter
      TofHitStatus* hitStatus2 = new TofHitStatus;
      if( m_hitCase == OuterLayer ) {
    hitStatus2->setBarrelCluster( 2 );        // outerlayer cluster and counter
      }
      else if( m_hitCase == DoubleLayer ) {
    if( innerEast && innerWest ) {
      hitStatus2->setBarrelCounter( 2 );      // outerlayer counter
    }
    else {
      hitStatus2->setBarrelCluster( 2 );      // outerlayer cluster and counter
    }
      }
      else {
    cout << "TofRec::TofTrack:convert2RecTofTrackCol: 2- Impossible!" << endl;
      }
      atrack2->setStatus( hitStatus2->value() );
      delete hitStatus2;
      recTofTrackCol->push_back( atrack2 );
    }
 
    if( innerEast && innerWest && outerEast && outerWest ) {
      RecTofTrack* atrack  = new RecTofTrack;
      buildRecTofTrack( atrack, 0 );              // doublelayer cluster
      TofHitStatus* hitStatus = new TofHitStatus;
      hitStatus->setBarrelCluster( 3 );           // doublelayer cluster
      atrack->setStatus( hitStatus->value() );
      delete hitStatus;
      recTofTrackCol->push_back( atrack );
    }
 
  }
 
  if( ( m_hitCase == EastEndcap ) || ( m_hitCase == WestEndcap ) ) {
    RecTofTrack* atrack = new RecTofTrack;
    buildRecTofTrack( atrack, 11 );               // eastendcap counter
    TofHitStatus* hitStatus = new TofHitStatus;
    if( m_hitCase == EastEndcap ) {
      hitStatus->setEndcapCluster( true );        // east endcap cluster counter readout
    }
    else if( m_hitCase == WestEndcap ) {
      hitStatus->setEndcapCluster( false );       // west endcap cluster counter readout
    }
    else {
      cout << "TofRec::TofTrack:convert2RecTofTrackCol: endcap- Impossible!" << endl;
    }
    atrack->setStatus( hitStatus->value() );
    delete hitStatus;
    recTofTrackCol->push_back( atrack );
  }
 
  if( m_hitCase == NoHit ) {
    RecTofTrack* atrack = new RecTofTrack;
    buildRecTofTrack( atrack, 3 );            // no hit
    TofHitStatus* hitStatus = new TofHitStatus;
    hitStatus->setNoHit();                    // no hit
    atrack->setStatus( hitStatus->value() );
    delete hitStatus;
    recTofTrackCol->push_back( atrack );
  }
 
  return;
}

◆ errzr1() [1/2]

double TofTrack::errzr1 ( ) const

inline

Definition at line 37 of file TofTrack.h.

37{ return m_errzr1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ errzr1() [2/2]

double TofTrack::errzr1 ( unsigned int i ) const

inline

Definition at line 42 of file TofTrack.h.

42{ return m_ezr1[i]; }

◆ errzr2() [1/2]

double TofTrack::errzr2 ( ) const

inline

Definition at line 38 of file TofTrack.h.

38{ return m_errzr2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ errzr2() [2/2]

double TofTrack::errzr2 ( unsigned int i ) const

inline

Definition at line 43 of file TofTrack.h.

43{ return m_ezr2[i]; }

◆ estime()

double TofTrack::estime ( ) const

inline

Definition at line 99 of file TofTrack.h.

99{ return m_estime; }

◆ findTofDataBarrel()

void TofTrack::findTofDataBarrel	(	std::vector< TofData * >	tofDataVec1,
		std::vector< TofData * >	tofDataVec2,
		double	zrhit,
		unsigned int	iflag,
		std::vector< TofTrack * > *&	tofTrackVec
	)

Definition at line 609 of file TofTrack.cxx.

                                                                                                                                                                         {
 
  unsigned int qual = 0xf;
  TofData* tof = 0;
  if( tofDataVec2.size() == 0 ) {
    if( tofDataVec1.size() == 0 ) {
      qual = 0;
    }
    else if( tofDataVec1.size() == 1 ) {
      std::vector<TofData*>::iterator iter1 = tofDataVec1.begin();
      tof = (*iter1);
      qual = 0x1;
    }
    else if( tofDataVec1.size() > 1 ) {
      tof= chooseTofData( tofDataVec1, zrhit ); 
      qual = 0x2;
    }
    else {
      cout << "TofRec::TofTrack::findTofDataBarrel: 1- Impossible!" << endl;
    }
  }
  else if( ( tofDataVec2.size() == 1 ) ) {
    if( tofDataVec1.size() == 0 ) {
      std::vector<TofData*>::iterator iter2 = tofDataVec2.begin();
      tof = (*iter2);
      qual = 0x4;
    }
    else if( tofDataVec1.size() == 1 ) {
      std::vector<TofData*>::iterator iter1 = tofDataVec1.begin();
      if( ((*iter1)->quality()&0x1ff)==0x01f && abs((*iter1)->ztdc()-zrhit)<ztdc_Cut ) {
    tof = (*iter1);
      }
      else {
    std::vector<TofData*>::iterator iter2 = tofDataVec2.begin();
        tof = compareTofData( (*iter1), (*iter2), zrhit );
      }
      qual = 0x5;
    }
    else if( tofDataVec1.size() > 1 ) {
      TofData* tofData1 = chooseTofData( tofDataVec1, zrhit ); 
      if( (tofData1->quality()&0x1ff)==0x01f && abs(tofData1->ztdc()-zrhit)<ztdc_Cut ) {
    tof = tofData1;
      }
      else {
    std::vector<TofData*>::iterator iter2 = tofDataVec2.begin();
    tof = compareTofData( tofData1, (*iter2), zrhit );
      }
      qual = 0x6;
    }
    else {
      cout << "TofRec::TofTrack::findTofDataBarrel: 2- Impossible!" << endl;
    }
  }
  else if( ( tofDataVec2.size() > 1 ) ) {
    if( tofDataVec1.size() == 0 ) {
      tof = chooseTofData( tofDataVec2, zrhit );
      qual = 0x8;
    }
    else if( tofDataVec1.size() == 1 ) {
      std::vector<TofData*>::iterator iter1 = tofDataVec1.begin();
      if( ((*iter1)->quality()&0x1ff)==0x01f && abs((*iter1)->ztdc()-zrhit)<ztdc_Cut ) {
    tof = (*iter1);
      }
      else {
    TofData* tofData2 = chooseTofData( tofDataVec2, zrhit ); 
    tof = compareTofData( (*iter1), tofData2, zrhit );
      }
      qual = 0x9;
    }
    else if( tofDataVec1.size() > 1 ) {
      TofData* tofData1 = chooseTofData( tofDataVec1, zrhit ); 
      if( (tofData1->quality()&0x1ff)==0x01f && abs(tofData1->ztdc()-zrhit)<ztdc_Cut ) {
    tof = tofData1;
      }
      else {
    TofData* tofData2 = chooseTofData( tofDataVec2, zrhit ); 
    tof = compareTofData( tofData1, tofData2, zrhit );
      }
      qual = 0xa;
    }
    else {
      cout << "TofRec::TofTrack::findTofDataBarrel: 3- Impossible!" << endl;
    }
  }
 
  if( qual != 0 ) {
    if( !(tof->used()) ) {
      getTofData( tof, iflag );
    }
    else {
      bool z1=false,  z2=false;
      bool zc1=false, zc2=false;
      TofTrack* track=0;
      if( iflag==1 ) {
    z1  = ( abs( m_zrhit1 - tof->ztdc() ) < ztdc_Cut );
    zc1 = ( m_zrhit1 > tof->ztdc() );
    std::vector<TofTrack*>::iterator iter = tofTrackVec->begin(); 
    for( ; iter!=tofTrackVec->end(); iter++ ) {
      if( tof->tofId()==(*iter)->tofId1() ) {
        track = (*iter);
        z2    = ( abs( (*iter)->zrhit1() - tof->ztdc() ) < ztdc_Cut );
        zc2   = ( (*iter)->zrhit1() > tof->ztdc() );
      }
    }
      }
      else if( iflag==2 ) {
    z1  = ( abs( m_zrhit2 - tof->ztdc() ) < ztdc_Cut );
    zc1 = ( m_zrhit2 > tof->ztdc() );
    std::vector<TofTrack*>::iterator iter = tofTrackVec->begin(); 
    for( ; iter!=tofTrackVec->end(); iter++ ) {
      if( tof->tofId()==(*iter)->tofId2() ) {
        track = (*iter);
        z2    = ( abs( (*iter)->zrhit2() - tof->ztdc() ) < ztdc_Cut );
        zc2   = ( (*iter)->zrhit2() > tof->ztdc() );
      }
    }
      }
      
      if( ( z1 && z2 )||( (!z1) && (!z2) ) ) {
    if( zc1 && !zc2 ) {
      getTofDataEast( tof, iflag );
      track->getTofDataWest( tof, iflag );
    }
    else if( !zc1 && zc2 ) {
      getTofDataWest( tof, iflag );
      track->getTofDataEast( tof, iflag );
    }
      }
      else if( z1 && !z2 ) {
    getTofData( tof, iflag );
    track->getTofDataNohit( iflag );
    cout << iflag << " " << track->quality1() << " " << track->quality2() << endl;
      }
      else if( !z1 && z2 ) {
    qual = 0;
      }
    }
  }
 
  if( qual == 0 ) {
    if( ( iflag == 1 ) || ( iflag == 3 ) ) {
      m_quality1 = ( m_quality1 | 0x300 );
    }
    else if( iflag == 2 ) {
      m_quality2 = ( m_quality2 | 0x300 );
    }
    else {
      cout << "TofRec::TofTrack::findTofDataBarrel: the 1- IFLAG is Out of Range!" << endl;
    }
  }
  else {
    qual = ( qual << 12 );
    if( ( iflag == 1 ) || ( iflag == 3 ) ) {
      m_quality1 = ( m_quality1 | qual );
    }
    else if( iflag == 2 ) {
      m_quality2 = ( m_quality2 | qual );
    }
    else {
      cout << "TofRec::TofTrack::findTofDataBarrel: the 2- IFLAG is Out of Range!" << endl;
    }
  }
 
  return;
}

Referenced by match().

◆ findTofDataEndcap()

void TofTrack::findTofDataEndcap	(	std::vector< TofData * >	tofDataVec1,
		std::vector< TofData * >	tofDataVec2,
		double	zr1[5]
	)

Definition at line 859 of file TofTrack.cxx.

                                                                                                                  {
 
  unsigned int iflag = 3;
  unsigned int qual = 0xf;
 
  if( tofDataVec2.size() == 0 ) {
    if( tofDataVec1.size() == 0 ) {
      qual = 0;
    }
    else if( tofDataVec1.size() == 1 ) {
      std::vector<TofData*>::iterator iter1 = tofDataVec1.begin();
      getTofData( (*iter1), iflag );
      qual = 0x1;
    }
    else if( tofDataVec1.size() > 1 ) {
      getTofData( chooseTofDataEndcap( tofDataVec1, zr1 ), iflag ); 
      qual = 0x2;
    }
    else {
      cout << "TofRec::TofTrack::findTofDataEndcap: 1- Impossible!" << endl;
    }
  }
  else if( ( tofDataVec2.size() == 1 ) ) {
    if( tofDataVec1.size() == 0 ) {
      std::vector<TofData*>::iterator iter2 = tofDataVec2.begin();
      getTofData( (*iter2), iflag );
      qual = 0x4;
    }
    else if( tofDataVec1.size() == 1 ) {
      std::vector<TofData*>::iterator iter1 = tofDataVec1.begin();
      std::vector<TofData*>::iterator iter2 = tofDataVec2.begin();
      getTofData( compareTofDataEndcap( (*iter1), (*iter2) ), iflag );
      qual = 0x5;
    }
    else if( tofDataVec1.size() > 1 ) {
      TofData* tofData1 = chooseTofDataEndcap( tofDataVec1, zr1 ); 
      std::vector<TofData*>::iterator iter2 = tofDataVec2.begin();
      getTofData( compareTofDataEndcap( tofData1, (*iter2) ), iflag );
      qual = 0x6;
    }
    else {
      cout << "TofRec::TofTrack::findTofDataBarrel: 2- Impossible!" << endl;
    }
  }
  else if( ( tofDataVec2.size() > 1 ) ) {
    if( tofDataVec1.size() == 0 ) {
      getTofData( chooseTofDataEndcap( tofDataVec2, zr1 ), iflag ); 
      qual = 0x8;
    }
    else if( tofDataVec1.size() == 1 ) {
      std::vector<TofData*>::iterator iter1 = tofDataVec1.begin();
      TofData* tofData2 = chooseTofDataEndcap( tofDataVec2, zr1 ); 
      getTofData( compareTofDataEndcap( (*iter1), tofData2 ), iflag );
      qual = 0x9;
    }
    else if( tofDataVec1.size() > 1 ) {
      TofData* tofData1 = chooseTofDataEndcap( tofDataVec1, zr1 ); 
      TofData* tofData2 = chooseTofDataEndcap( tofDataVec2, zr1 ); 
      getTofData( compareTofDataEndcap( tofData1, tofData2 ), iflag );
      qual = 0xa;
    }
    else {
      cout << "TofRec::TofTrack::findTofDataBarrel: 3- Impossible!" << endl;
    }
  }
 
  if( qual == 0 ) {
    m_quality1 = ( m_quality1 | 0x300 );
  }
  else {
    qual = ( qual << 12 );
    m_quality1 = ( m_quality1 | qual );
  }
 
  return;
}

Referenced by match().

◆ flag()

unsigned int TofTrack::flag ( ) const

inline

Definition at line 105 of file TofTrack.h.

105{ return m_flag; }

Referenced by setFlag().

◆ getMultiHit()

void TofTrack::getMultiHit ( TofTrack *& track )

Definition at line 244 of file TofTrack.cxx.

                                             {
  if( m_hitCase == InnerLayer || m_hitCase == OuterLayer || m_hitCase == DoubleLayer ) {
 
    if( ( m_hitCase==InnerLayer || m_hitCase==DoubleLayer ) && ( track->hitCase()==InnerLayer || track->hitCase()==DoubleLayer ) ) {
      if( ( abs(m_id1-track->id1())<=1 ) || ( m_id1==0 && track->id1()==87 ) || ( m_id1==87 && track->id1()==0 ) ) {
    track->setQuality1( ( track->quality1() | 0x400 ) );
    m_quality1 = ( m_quality1 | 0x400 );    
      }
    }
 
    if( ( m_hitCase==OuterLayer || m_hitCase==DoubleLayer ) && ( track->hitCase()==OuterLayer || track->hitCase()==DoubleLayer ) ) {
      if( ( abs(m_id2-track->id2())<=1 ) || ( m_id2==88 && track->id2()==175 ) || ( m_id2==175 && track->id2()==88 ) ) {
    track->setQuality2( ( track->quality2() | 0x400 ) );
    m_quality2 = ( m_quality2 | 0x400 );    
      }
    }
 
  }
  else if( m_hitCase == EastEndcap ) {
    if( track->hitCase()==EastEndcap ) {
      if( ( abs(m_id1-track->id1())<=1 ) || ( m_id1==0 && track->id1()==47 ) || ( m_id1==47 && track->id1()==0 ) ) {
    track->setQuality1( ( track->quality1() | 0x400 ) );
    m_quality1 = ( m_quality1 | 0x400 );    
      } 
    }
  }
  else if( m_hitCase == WestEndcap ) {
    if( track->hitCase()==WestEndcap ) {
      if( ( abs(m_id1-track->id1())<=1 ) || ( m_id1==48 && track->id1()==95 ) || ( m_id1==95 && track->id1()==48 ) ) {
    track->setQuality1( ( track->quality1() | 0x400 ) );
    m_quality1 = ( m_quality1 | 0x400 );    
      } 
    }
  }
 
  return;
}

Referenced by TofRec::execute().

◆ getTofData()

void TofTrack::getTofData	(	TofData *	tof,
		unsigned int	iflag
	)

Definition at line 995 of file TofTrack.cxx.

                                                            {
 
  if( iflag == 1 ) { 
    m_tofId1   = tof->tofId();
    m_qch1     = tof->adcChannelEast();
    m_adc1     = tof->adc1();
    m_tdc1     = tof->tdc1();
    m_qch2     = tof->adcChannelWest();
    m_adc2     = tof->adc2();
    m_tdc2     = tof->tdc2();
    m_ztdc1    = tof->ztdc();
    m_zadc1    = tof->zadc();
    m_quality1 = ( m_quality1 | ( 0x1f & tof->quality() ) );
    if( ( ( tof->quality() & 0x5 ) != 0x5 ) || ( ( ( tof->quality() & 0x5 ) == 0x5 ) && ( abs( tof->ztdc() - m_zrhit1 ) > ztdc_Cut ) ) ) { 
      m_quality1 = ( m_quality1 | 0x100 );
    }
    if( ( ( tof->quality() & 0xa ) != 0xa ) || ( ( ( tof->quality() & 0xa ) == 0xa ) && ( abs( tof->zadc() - m_zrhit1 ) > zadc_Cut ) ) ) { 
      m_quality1 = ( m_quality1 | 0x200 );
    }
  }
  else if( iflag == 2 ) {
    m_tofId2   = tof->tofId();
    m_qch3     = tof->adcChannelEast();
    m_adc3     = tof->adc1();
    m_tdc3     = tof->tdc1();
    m_qch4     = tof->adcChannelWest();
    m_adc4     = tof->adc2();
    m_tdc4     = tof->tdc2();
    m_ztdc2    = tof->ztdc();
    m_zadc2    = tof->zadc();
    m_quality2 = ( m_quality2 | ( 0x1f & tof->quality() ) );
    if( ( ( tof->quality() & 0x5 ) != 0x5 ) || ( ( ( tof->quality() & 0x5 ) == 0x5 ) && ( abs( tof->ztdc() - m_zrhit2 ) > ztdc_Cut ) ) ) { 
      m_quality2 = ( m_quality2 | 0x100 );
    }
    if( ( ( tof->quality() & 0xa ) != 0xa ) || ( ( ( tof->quality() & 0xa ) == 0xa ) && ( abs( tof->zadc() - m_zrhit2 ) > zadc_Cut ) ) ) { 
      m_quality2 = ( m_quality2 | 0x200 );
    }
  }
  else if( iflag == 3 ) {
    m_tofId1   = tof->tofId();
    m_qch1     = tof->adcChannel();
    m_adc1     = tof->adc();
    m_tdc1     = tof->tdc();
    m_quality1 = ( m_quality1 | ( 0x1f & tof->quality() ) );
    m_quality1 = ( m_quality1 | 0x300 );
  }
  else {
    cout << "TofRec::TofTrack::getTofData: Flag which sign the Barrel/Endcap or Inner/Outer is wrong! Please check it!" << endl;
  }
  tof->setUsed();
  return;
}

Referenced by findTofDataBarrel(), and findTofDataEndcap().

◆ getTofDataEast()

void TofTrack::getTofDataEast	(	TofData *	tof,
		unsigned int	iflag
	)

Definition at line 1052 of file TofTrack.cxx.

                                                                {
 
  if( iflag == 1 ) { 
    m_tofId1   = tof->tofId();
    m_qch1     = tof->adcChannelEast();
    m_adc1     = tof->adc1();
    m_tdc1     = tof->tdc1();
    m_qch2     = -999.0;
    m_adc2     = -999.0;
    m_tdc2     = -999.0;
    m_ztdc1    = tof->ztdc();
    m_zadc1    = tof->zadc();
    m_quality1 = ( ( m_quality1 & 0xfffffff0 ) | ( 0x1c & tof->quality() ) );
    if( ( ( tof->quality() & 0x5 ) != 0x5 ) || ( ( ( tof->quality() & 0x5 ) == 0x5 ) && ( abs( tof->ztdc() - m_zrhit1 ) > ztdc_Cut ) ) ) { 
      m_quality1 = ( m_quality1 | 0x100 );
    }
    if( ( ( tof->quality() & 0xa ) != 0xa ) || ( ( ( tof->quality() & 0xa ) == 0xa ) && ( abs( tof->zadc() - m_zrhit1 ) > zadc_Cut ) ) ) { 
      m_quality1 = ( m_quality1 | 0x200 );
    }
  }
  else if( iflag == 2 ) {
    m_tofId2   = tof->tofId();
    m_qch3     = tof->adcChannelEast();
    m_adc3     = tof->adc1();
    m_tdc3     = tof->tdc1();
    m_qch4     = -999.0;
    m_adc4     = -999.0;
    m_tdc4     = -999.0;
    m_ztdc2    = tof->ztdc();
    m_zadc2    = tof->zadc();
    m_quality2 = ( ( m_quality2 & 0xfffffff0 ) | ( 0x1c & tof->quality() ) );
    if( ( ( tof->quality() & 0x5 ) != 0x5 ) || ( ( ( tof->quality() & 0x5 ) == 0x5 ) && ( abs( tof->ztdc() - m_zrhit2 ) > ztdc_Cut ) ) ) { 
      m_quality2 = ( m_quality2 | 0x100 );
    }
    if( ( ( tof->quality() & 0xa ) != 0xa ) || ( ( ( tof->quality() & 0xa ) == 0xa ) && ( abs( tof->zadc() - m_zrhit2 ) > zadc_Cut ) ) ) { 
      m_quality2 = ( m_quality2 | 0x200 );
    }
  }
  else {
    cout << "TofRec::TofTrack::getTofDataEast: Flag which sign the Barrel/Endcap or Inner/Outer is wrong! Please check it!" << endl;
  }
  tof->setUsed();
  return;
}

Referenced by findTofDataBarrel().

◆ getTofDataNohit()

void TofTrack::getTofDataNohit ( unsigned int iflag )

Definition at line 1150 of file TofTrack.cxx.

                                                   {
 
  if( iflag == 1 ) { 
    m_tofId1   = -99;
    m_qch1     = -999.0;
    m_adc1     = -999.0;
    m_tdc1     = -999.0;
    m_qch2     = -999.0;
    m_adc2     = -999.0;
    m_tdc2     = -999.0;
    m_ztdc1    = -999.0;
    m_zadc1    = -999.0;
    m_quality1 = ( m_quality1 & 0x700 );
    if( m_hitCase == InnerLayer ) { m_hitCase = NoHit; }
    else if( m_hitCase == DoubleLayer ) { m_hitCase = OuterLayer; }
  }
  else if( iflag == 2 ) {
    m_tofId2   = -99;
    m_qch3     = -999.0;
    m_adc3     = -999.0;
    m_tdc3     = -999.0;
    m_qch4     = -999.0;
    m_adc4     = -999.0;
    m_tdc4     = -999.0;
    m_ztdc2    = -999.0;
    m_zadc2    = -999.0;
    m_quality2 = ( m_quality2 & 0x700 );
    if( m_hitCase == OuterLayer ) { m_hitCase = NoHit; }
    else if( m_hitCase == DoubleLayer ) { m_hitCase = InnerLayer; }
  }
  else {
    cout << "TofRec::TofTrack::getTofData: Flag which sign the Barrel/Endcap or Inner/Outer is wrong! Please check it!" << endl;
  }
 
  return;
}

Referenced by findTofDataBarrel().

◆ getTofDataWest()

void TofTrack::getTofDataWest	(	TofData *	tof,
		unsigned int	iflag
	)

Definition at line 1101 of file TofTrack.cxx.

                                                                {
 
  if( iflag == 1 ) { 
    m_tofId1   = tof->tofId();
    m_qch1     = -999.0;
    m_adc1     = -999.0;
    m_tdc1     = -999.0;
    m_qch2     = tof->adcChannelWest();
    m_adc2     = tof->adc2();
    m_tdc2     = tof->tdc2();
    m_ztdc1    = tof->ztdc();
    m_zadc1    = tof->zadc();
    m_quality1 = ( ( m_quality1 & 0xfffffff0 ) | ( 0x13 & tof->quality() ) );
    if( ( ( tof->quality() & 0x5 ) != 0x5 ) || ( ( ( tof->quality() & 0x5 ) == 0x5 ) && ( abs( tof->ztdc() - m_zrhit1 ) > ztdc_Cut ) ) ) { 
      m_quality1 = ( m_quality1 | 0x100 );
    }
    if( ( ( tof->quality() & 0xa ) != 0xa ) || ( ( ( tof->quality() & 0xa ) == 0xa ) && ( abs( tof->zadc() - m_zrhit1 ) > zadc_Cut ) ) ) { 
      m_quality1 = ( m_quality1 | 0x200 );
    }
  }
  else if( iflag == 2 ) {
    m_tofId2   = tof->tofId();
    m_qch3     = -999.0;
    m_adc3     = -999.0;
    m_tdc3     = -999.0;
    m_qch4     = tof->adcChannelWest();
    m_adc4     = tof->adc2();
    m_tdc4     = tof->tdc2();
    m_ztdc2    = tof->ztdc();
    m_zadc2    = tof->zadc();
    m_quality2 = ( ( m_quality2 & 0xfffffff0 ) | ( 0x13 & tof->quality() ) );
    if( ( ( tof->quality() & 0x5 ) != 0x5 ) || ( ( ( tof->quality() & 0x5 ) == 0x5 ) && ( abs( tof->ztdc() - m_zrhit2 ) > ztdc_Cut ) ) ) { 
      m_quality2 = ( m_quality2 | 0x100 );
    }
    if( ( ( tof->quality() & 0xa ) != 0xa ) || ( ( ( tof->quality() & 0xa ) == 0xa ) && ( abs( tof->zadc() - m_zrhit2 ) > zadc_Cut ) ) ) { 
      m_quality2 = ( m_quality2 | 0x200 );
    }
  }
  else {
    cout << "TofRec::TofTrack::getTofDataWest: Flag which sign the Barrel/Endcap or Inner/Outer is wrong! Please check it!" << endl;
  }
  tof->setUsed();
  return;
}

Referenced by findTofDataBarrel().

◆ hitCase()

ExtTrackCase TofTrack::hitCase ( ) const

inline

Definition at line 30 of file TofTrack.h.

30{ return m_hitCase; }

Referenced by TofCheckDigi::Fill_TofTrack(), getMultiHit(), TofCount::setTrack1(), TofCount::setTrack2(), TofCount::setTrack3(), and TofCount::setTrack4().

◆ id1()

int TofTrack::id1 ( ) const

inline

Definition at line 27 of file TofTrack.h.

27{ return m_id1; }

Referenced by TofCheckDigi::Fill_TofTrack(), and getMultiHit().

◆ id2()

int TofTrack::id2 ( ) const

inline

Definition at line 28 of file TofTrack.h.

28{ return m_id2; }

Referenced by TofCheckDigi::Fill_TofTrack(), and getMultiHit().

◆ isNoHit()

bool TofTrack::isNoHit ( ) const

inline

Definition at line 108 of file TofTrack.h.

108{ return m_hitCase==NoHit; }

◆ kal()

int TofTrack::kal ( unsigned int i ) const

inline

Definition at line 39 of file TofTrack.h.

39{ return m_kal[i]; }

Referenced by TofCheckDigi::Fill_TofTrack(), and setExtTrack().

◆ match()

void TofTrack::match	(	bool	forCalibration,
		std::vector< int >	deadId,
		std::vector< TofTrack * > *&	tofTrackVec
	)

Definition at line 493 of file TofTrack.cxx.

                                                                                                     {
 
  if( m_hitCase == NoHit ) return;
 
  if( m_hitCase == InnerLayer ) {
    findTofDataBarrel( m_tofData1, m_tofData2, m_zrhit1, 1, tofTrackVec );
    if( ( m_quality1 & 0x10 ) == 0 ) { m_hitCase = NoHit; }
  }
  else if( m_hitCase == OuterLayer ) {
    findTofDataBarrel( m_tofData3, m_tofData4, m_zrhit2, 2, tofTrackVec );
    if( ( m_quality2 & 0x10 ) == 0 ) { m_hitCase = NoHit; }
  }
  else if( m_hitCase == DoubleLayer ) {
    findTofDataBarrel( m_tofData1, m_tofData2, m_zrhit1, 1, tofTrackVec );
    if( ( m_quality1 & 0x10 ) == 0 ) { m_hitCase = OuterLayer; }
    findTofDataBarrel( m_tofData3, m_tofData4, m_zrhit2, 2, tofTrackVec );
    if( ( m_quality2 & 0x10 ) == 0 ) {
      if( m_hitCase == DoubleLayer ) {
    m_hitCase = InnerLayer;
      }
      else if( m_hitCase == OuterLayer ) {
    m_hitCase = NoHit;
      }
      else {
    cout << "TofRec::TofTrack::match: 2- Impossible!" << endl;
      }
    }
  }
  else if( ( m_hitCase == EastEndcap ) || ( m_hitCase == WestEndcap ) ) {
    findTofDataEndcap( m_tofData1, m_tofData2, m_zr1 );
  }
  else {
    cout << "TofRec::TofTrack::match: 1- Impossible!" << endl; 
  }
 
  if( forCalibration ) {
    // set Data Sample for Calibration, double layer, only one hit for counter, T and Q.
    if( m_hitCase == DoubleLayer ) {
      if( ( ( m_quality1 & 0xf ) == 0xf ) && ( ( m_quality2 & 0xf ) == 0xf ) ) {
    m_quality1 = ( m_quality1 | 0x800 );        // Calibration Sample
    m_quality2 = ( m_quality2 | 0x800 );        // Calibration Sample
      }
      else {
    std::vector<int>::iterator iter = deadId.begin();
    for( ; iter != deadId.end(); iter++ ) {
      Identifier iden = Identifier(*iter);
      int barrel = TofID::barrel_ec(iden);
      int layer  = TofID::layer(iden);
      int tofId  = TofID::phi_module(iden);
      int east   = TofID::end(iden);
      if( barrel == 1 ) {
        if( layer==0 ) {
          if( m_tofId1 == tofId ) {
        if( ( m_quality2 & 0xf ) == 0xf ) {
          if( ( ( east == 0 ) && ( ( m_quality1 & 0xf ) == 0x3 ) ) || ( ( east == 1 ) && ( ( m_quality1 & 0xf ) == 0xc ) ) ) {
            m_quality1 = ( m_quality1 | 0x800 );
            m_quality2 = ( m_quality2 | 0x800 );
 
          }
        }
          }
        }
        else if( layer == 1 ) {
          if( m_tofId2 == (tofId+88) ) {
        if( ( m_quality1 & 0xf ) == 0xf ) {
          if( ( ( east == 0 ) && ( ( m_quality2 & 0xf ) == 0x3 ) ) || ( ( east == 1 ) && ( ( m_quality2 & 0xf ) == 0xc ) ) ) {
            m_quality1 = ( m_quality1 | 0x800 );
            m_quality2 = ( m_quality2 | 0x800 );
          }
        }
          }
        }
      }
    }
      }
    }
    // set Data Sample for Calibration, only one hit for counter, T and Q.
    else if( m_hitCase == InnerLayer ) {
      if( ( m_quality1 & 0xf ) == 0xf ) {
    m_quality1 = ( m_quality1 | 0x800 );        // Calibration Sample
      }
      else {
    std::vector<int>::iterator iter = deadId.begin();
    for( ; iter != deadId.end(); iter++ ) {
      Identifier iden = Identifier(*iter);
      int barrel = TofID::barrel_ec(iden);
      int layer  = TofID::layer(iden);
      int tofId  = TofID::phi_module(iden);
      int east   = TofID::end(iden);
      if( barrel == 1 ) {
        if( layer==0 ) {
          if( m_tofId1 == tofId ) {
        if( ( ( east == 0 ) && ( ( m_quality1 & 0xf ) == 0x3 ) ) || ( ( east == 1 ) && ( ( m_quality1 & 0xf ) == 0xc ) ) ) {
          m_quality1 = ( m_quality1 | 0x800 );
        }
          }
        }
      }
    }
      }
    }
 
    // set Data Sample for Calibration, only one hit for counter, T and Q.
    if( ( ( m_hitCase == EastEndcap ) || ( m_hitCase == WestEndcap ) ) && ( ( m_quality1 & 0xf ) == 0xc ) ) {
      m_quality1 = ( m_quality1 | 0x800 );        // Calibration Sample
    }
  }
 
  return;
}

◆ p()

double TofTrack::p ( ) const

inline

Definition at line 31 of file TofTrack.h.

31{ return m_momentum; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ path()

double TofTrack::path ( ) const

inline

Definition at line 32 of file TofTrack.h.

32{ return m_path; }

◆ path1()

double TofTrack::path1 ( ) const

inline

Definition at line 33 of file TofTrack.h.

33{ return m_path1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ path2()

double TofTrack::path2 ( ) const

inline

Definition at line 34 of file TofTrack.h.

34{ return m_path2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ ph()

double TofTrack::ph ( ) const

inline

Definition at line 67 of file TofTrack.h.

67{ return m_ph; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ ph1()

double TofTrack::ph1 ( ) const

inline

Definition at line 65 of file TofTrack.h.

65{ return m_ph1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ ph11()

double TofTrack::ph11 ( ) const

inline

Definition at line 61 of file TofTrack.h.

61{ return m_ph11; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ ph12()

double TofTrack::ph12 ( ) const

inline

Definition at line 62 of file TofTrack.h.

62{ return m_ph12; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ ph2()

double TofTrack::ph2 ( ) const

inline

Definition at line 66 of file TofTrack.h.

66{ return m_ph2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ ph21()

double TofTrack::ph21 ( ) const

inline

Definition at line 63 of file TofTrack.h.

63{ return m_ph21; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ ph22()

double TofTrack::ph22 ( ) const

inline

Definition at line 64 of file TofTrack.h.

64{ return m_ph22; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ phi1()

double TofTrack::phi1 ( ) const

inline

Definition at line 52 of file TofTrack.h.

52{ return m_phi1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ phi2()

double TofTrack::phi2 ( ) const

inline

Definition at line 53 of file TofTrack.h.

53{ return m_phi2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ qch1()

double TofTrack::qch1 ( ) const

inline

Definition at line 77 of file TofTrack.h.

77{ return m_qch1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ qch2()

double TofTrack::qch2 ( ) const

inline

Definition at line 78 of file TofTrack.h.

78{ return m_qch2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ qch3()

double TofTrack::qch3 ( ) const

inline

Definition at line 79 of file TofTrack.h.

79{ return m_qch3; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ qch4()

double TofTrack::qch4 ( ) const

inline

Definition at line 80 of file TofTrack.h.

80{ return m_qch4; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ quality()

unsigned int TofTrack::quality ( ) const

inline

Definition at line 56 of file TofTrack.h.

56{ return m_quality; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ quality1()

unsigned int TofTrack::quality1 ( ) const

inline

Definition at line 54 of file TofTrack.h.

54{ return m_quality1; }

Referenced by TofCheckDigi::Fill_TofTrack(), findTofDataBarrel(), and getMultiHit().

◆ quality2()

unsigned int TofTrack::quality2 ( ) const

inline

Definition at line 55 of file TofTrack.h.

55{ return m_quality2; }

Referenced by TofCheckDigi::Fill_TofTrack(), findTofDataBarrel(), and getMultiHit().

◆ qualityAnalysis()

void TofTrack::qualityAnalysis ( )

◆ setCalibration()

void TofTrack::setCalibration ( )

Definition at line 1193 of file TofTrack.cxx.

                              {
 
  bool barrel = ( ( m_hitCase == InnerLayer ) || ( m_hitCase == OuterLayer ) || ( m_hitCase == DoubleLayer ) );
  bool endcap = ( ( m_hitCase == EastEndcap ) || ( m_hitCase == WestEndcap ) );
 
  bool innerEast  = ( ( m_quality1 & 0xc ) == 0xc );
  bool innerWest  = ( ( m_quality1 & 0x3 ) == 0x3 );
  bool outerEast  = ( ( m_quality2 & 0xc ) == 0xc );
  bool outerWest  = ( ( m_quality2 & 0x3 ) == 0x3 );
  bool innerLayer = ( ( m_quality1 & 0xf ) == 0xf );
  bool outerLayer = ( ( m_quality2 & 0xf ) == 0xf );
 
  bool endcapData = ( ( m_quality1 & 0xc ) == 0xc );
 
  if( m_hitCase == DoubleLayer ) {
    for( unsigned int i=0; i<5; i++ ) {
      m_texp[i] = tofCaliSvc->BTimeCluster( m_texpInner[i], m_texpOuter[i], m_zr1[i], m_zr2[i], m_tofId1, m_tofId2 );
    }
    m_path    =  tofCaliSvc->BTimeCluster( m_path1, m_path2, m_zrhit1, m_zrhit2, m_tofId1, m_tofId2 );
  }
 
  if( barrel ) {
    if( innerEast ) {
      for( unsigned int i=0; i<5; i++ ) {
    m_tof11[i] = tofCaliSvc->BTime1( m_adc1, m_tdc1-m_estime, m_zr1[i], m_tofId1, m_estime );
      }
      m_sigma11 = tofCaliSvc->BSigma1( m_zrhit1, m_tofId1 );
      m_ph11 = m_adc1;
    }
    
    if( innerWest ) {
      for( unsigned int i=0; i<5; i++ ) {
    m_tof12[i] = tofCaliSvc->BTime2( m_adc2, m_tdc2-m_estime, m_zr1[i], m_tofId1, m_estime );
      }
      m_sigma12 = tofCaliSvc->BSigma2( m_zrhit1, m_tofId1 );
      m_ph12 = m_adc2;
    }
 
    if( innerLayer ) {
      for( unsigned int i=0; i<5; i++ ) {
    m_tof1[i]  = tofCaliSvc->BTimeCounter( m_tof11[i], m_tof12[i], m_zr1[i], m_tofId1 );
      }
      m_sigma1  = tofCaliSvc->BSigmaCounter( m_zrhit1, m_tofId1 );
      m_ph1  = tofCaliSvc->BPulseHeight( m_adc1, m_adc2, m_zrhit1, m_theta1, m_tofId1 );
      /*
      if( m_hitCase == InnerLayer ) {
    for( unsigned int i=0; i<5; i++ ) {
      m_tof[i] = m_tof1[i];
    }
    m_sigma   = m_sigma1;
      }
      */
    }
 
    if( outerEast ) {
      for( unsigned int i=0; i<5; i++ ) {
    m_tof21[i] = tofCaliSvc->BTime1( m_adc3, m_tdc3-m_estime, m_zr2[i], m_tofId2, m_estime );
      }
      m_sigma21 = tofCaliSvc->BSigma1( m_zrhit2, m_tofId2 );
      m_ph21 = m_adc3;
    }
 
    if( outerWest ) {
      for( unsigned int i=0; i<5; i++ ) {
    m_tof22[i] = tofCaliSvc->BTime2( m_adc4, m_tdc4-m_estime, m_zr2[i], m_tofId2, m_estime );
      }
      m_sigma22 = tofCaliSvc->BSigma2( m_zrhit2, m_tofId2 );
      m_ph22 = m_adc4;
    }
 
    if( outerLayer ) {
      for( unsigned int i=0; i<5; i++ ) {
    m_tof2[i]  = tofCaliSvc->BTimeCounter( m_tof21[i], m_tof22[i], m_zr2[i], m_tofId2 );
      }
      m_sigma2  = tofCaliSvc->BSigmaCounter( m_zrhit2, m_tofId2 );
      m_ph2  = tofCaliSvc->BPulseHeight( m_adc3, m_adc4, m_zrhit2, m_theta2, m_tofId2 );
      /*
      if( m_hitCase == OuterLayer ) {
    for( unsigned int i=0; i<5; i++ ) {
      m_tof[i] = m_tof2[i];
    }
    m_sigma   = m_sigma2;
      }
      */
    }
 
    if( innerLayer && outerLayer ) {
      for( unsigned int i=0; i<5; i++ ) {
    m_tof[i] = tofCaliSvc->BTimeCluster( m_tof1[i], m_tof2[i], m_zr1[i], m_zr2[i], m_tofId1, m_tofId2 );
      }
      m_sigma = tofCaliSvc->BSigmaCluster( m_zrhit1, m_zrhit2, m_tofId1, m_tofId2 );
      m_ph    = tofCaliSvc->BTimeCluster( m_ph1, m_ph2, m_zrhit1, m_zrhit2, m_tofId1, m_tofId2 );
    }
  }
 
  if( endcap ) {
    if( endcapData ) {
      for( unsigned int i=0; i<5; i++ ) {
    m_tof11[i] = tofCaliSvc->ETime( m_adc1, m_tdc1-m_estime, m_zr1[i], m_tofId1 );
      }
      m_sigma11 = tofCaliSvc->ESigma( m_zrhit1, m_tofId1 );
      m_ph11  = tofCaliSvc->EPulseHeight( m_adc1, m_zrhit1, m_theta1, m_tofId1 );
      m_quality = 1;
      if( (m_quality1&0xa000)!=0 ) { m_quality = 4; }
    }
  }
 
  // set Quality of Barrel TOF
  if( barrel ) {
 
    // double layer
    if( innerLayer && outerLayer ) {
      m_quality = 1;
    }
    else {
      // single layer
      if( innerLayer || outerLayer ) {
    m_quality = 2;
      }
      else {
    // single-end of one layer
        if( innerEast || innerWest || outerEast || outerWest ) {
      m_quality = 3;
    }
      }
    }
 
    // multi-hit
    if( ( (m_quality1&0xa000)!=0 ) || ( (m_quality2&0xa000)!=0 ) ) {
      m_quality = m_quality + 3;
    }
 
    // ztdc and extrapolated zhit is not matched
    if( ( (m_quality1&0x100)==0x100 ) || ( (m_quality2&0x100)==0x100 ) ) {
      if( ( m_quality == 1 ) || ( m_quality == 4 ) ) { m_quality = 7; }
      else if( ( m_quality == 2 ) || ( m_quality == 5 ) ) { m_quality = 8; }
      else if( ( m_quality == 3 ) || ( m_quality == 6 ) ) { m_quality = 9; }
      else {
    cout << "TofRec::TofTrack::setCalibration: Impossible!" << endl;
      }
    }
 
  }
 
  return;
}

◆ setExtTrack()

void TofTrack::setExtTrack	(	RecExtTrack *	extTrack,
		int	kal[5],
		double	t0,
		int	t0Stat
	)

Definition at line 144 of file TofTrack.cxx.

                                                                                     {
 
  m_estime = t0;
  m_t0Stat = t0Stat;
 
  int tofId1 = extTrack->tof1VolumeNumber();
  int tofId2 = extTrack->tof2VolumeNumber();
  
  m_barrel = 3;
  if( tofId1>=0 && tofId1<=87 ) {
    m_id1 = tofId1;
    m_barrel = 1;
    m_hitCase = InnerLayer;
  }
  else if( tofId1>=176 && tofId1<=223 ) {
    m_id1 = tofId1 - 176 + 48;
    m_barrel = 2;
    m_hitCase = WestEndcap;
  }
  else if( tofId1>=224 && tofId1<=271 ) { 
    m_id1 = tofId1 - 176 - 48;
    m_barrel = 0;
    m_hitCase = EastEndcap;
  }
  else {
    m_barrel = 3;
    m_hitCase = NoHit;
  }
 
  if( tofId2>=88 && tofId2<=175 ) {
    m_id2 = tofId2;
    m_barrel = 1;
    if( m_hitCase == InnerLayer ) {
      m_hitCase = DoubleLayer;
    }
    else if( m_hitCase==NoHit ) {
      m_hitCase = OuterLayer;
    }
  }
 
  m_trackId  = extTrack->trackId();
 
  m_momentum = extTrack->tof1Momentum().r();
  if( m_hitCase == InnerLayer || m_hitCase == DoubleLayer ) {
    for( unsigned int i=0; i<5; i++ ) {
      m_texpInner[i] = extTrack->tof1(i);
    }
    m_path1  = extTrack->tof1Path();
    m_zrhit1 = extTrack->tof1Position().z();
    m_errzr1 = extTrack->tof1PosSigmaAlongZ();
    m_theta1 = extTrack->tof1Momentum().rho()/extTrack->tof1Momentum().r();
    m_phi1   = extTrack->tof1Position().phi();
    for( unsigned int i=0; i<5; i++ ) {
      m_zr1[i]  = extTrack->tof1Position(i).z();
      m_ezr1[i] = extTrack->tof1PosSigmaAlongZ(i); 
    }
  }
  else if( m_hitCase == EastEndcap || m_hitCase == WestEndcap ) {
    for( unsigned int i=0; i<5; i++ ) {
      m_texpInner[i] = extTrack->tof1(i);
    }
    m_path1  = extTrack->tof1Path();
    m_zrhit1 = extTrack->tof1Position().rho();
    m_errzr1 = sqrt( extTrack->tof1PosSigmaAlongX()*extTrack->tof1PosSigmaAlongX() + extTrack->tof1PosSigmaAlongY()*extTrack->tof1PosSigmaAlongY() );
    m_theta1 = extTrack->tof1Momentum().z()/extTrack->tof1Momentum().r();
    m_phi1   = extTrack->tof1Position().phi();
    for( unsigned int i=0; i<5; i++ ) {
      m_zr1[i]  = extTrack->tof1Position(i).rho();
      m_ezr1[i] = sqrt( extTrack->tof1PosSigmaAlongX(i)*extTrack->tof1PosSigmaAlongX(i) + extTrack->tof1PosSigmaAlongY(i)*extTrack->tof1PosSigmaAlongY(i) ); 
    }
  }
 
  if( m_hitCase == OuterLayer || m_hitCase == DoubleLayer ) {
    for( unsigned int i=0; i<5; i++ ) {
      m_texpOuter[i] = extTrack->tof2(i);
    }
    m_path2  = extTrack->tof2Path();
    m_zrhit2 = extTrack->tof2Position().z();
    m_errzr2 = extTrack->tof2PosSigmaAlongZ();
    m_theta2 = extTrack->tof2Momentum().rho()/extTrack->tof2Momentum().r();
    m_phi2   = extTrack->tof2Position().phi();
    for( unsigned int i=0; i<5; i++ ) {
      m_zr2[i]  = extTrack->tof2Position(i).z();
      m_ezr2[i] = extTrack->tof2PosSigmaAlongZ(i); 
    }
  }
 
  if( m_hitCase == NoHit ) { m_quality = 11; }
 
  for( unsigned int i=0; i<5; i++ ) {
    m_kal[i] = kal[i];
  }
 
  return;
}

Referenced by TofRec::execute().

◆ setFlag()

void TofTrack::setFlag ( unsigned int flag )

inline

Definition at line 112 of file TofTrack.h.

112{ m_flag = flag; }

TofTrack::flag

unsigned int flag() const

Definition: TofTrack.h:105

◆ setQuality()

void TofTrack::setQuality ( int qual )

inline

Definition at line 109 of file TofTrack.h.

109{ m_quality = qual; }

◆ setQuality1()

void TofTrack::setQuality1 ( int qual1 )

inline

Definition at line 110 of file TofTrack.h.

110{ m_quality1 = qual1; }

Referenced by getMultiHit().

◆ setQuality2()

void TofTrack::setQuality2 ( int qual2 )

inline

Definition at line 111 of file TofTrack.h.

111{ m_quality2 = qual2; }

Referenced by getMultiHit().

◆ setRecTofTrack()

void TofTrack::setRecTofTrack	(	RecTofTrack *	track,
		int	layerorend
	)

Definition at line 1621 of file TofTrack.cxx.

                                                                  {
 
  double toffset[6];
  for( unsigned int i=0; i<6; i++ ) {
    toffset[i] = 0.0;
  }
 
  if( layerorend == 0 ) {   // cluster or double layer hit
    track->setPh( m_ph );
    track->setTof( m_tof[0] );
    track->setSigmaElectron( m_sigma );
    for( unsigned int i=0; i<5; i++ ) {
      toffset[i] = m_tof[0] - m_tof[i];
    }
    track->setToffset( toffset );
    track->setBeta( m_path/m_tof[0]/30.0 );
  }
  else if( layerorend == 1 ) { // inner layer
    track->setPh( m_ph1 );
    track->setTof( m_tof1[0] );
    track->setSigmaElectron( m_sigma1 );
    for( unsigned int i=0; i<5; i++ ) {
      toffset[i] = m_tof1[0] - m_tof1[i];
    }
    track->setToffset( toffset );
    track->setBeta( m_path1/m_tof1[0]/30.0 );
  }
  else if( layerorend == 2 ) { // outer layer
    track->setPh( m_ph2 );
    track->setTof( m_tof2[0] );
    track->setSigmaElectron( m_sigma2 );
    for( unsigned int i=0; i<5; i++ ) {
      toffset[i] = m_tof2[0] - m_tof2[i];
    }
    track->setToffset( toffset );
    track->setBeta( m_path2/m_tof2[0]/30.0 );
  }
  else if( layerorend == 11 ) { // inner layer east end readout
    track->setPh( m_ph11 );
    track->setTof( m_tof11[0] );
    track->setSigmaElectron( m_sigma11 );
    for( unsigned int i=0; i<5; i++ ) {
      toffset[i] = m_tof11[0] - m_tof11[i];
    }
    track->setToffset( toffset );
    track->setBeta( m_path1/m_tof11[0]/30.0 );
  }
  else if( layerorend == 12 ) { // inner layer west end readout
    track->setPh( m_ph12 );
    track->setTof( m_tof12[0] );
    track->setSigmaElectron( m_sigma12 );
    for( unsigned int i=0; i<5; i++ ) {
      toffset[i] = m_tof12[0] - m_tof12[i];
    }
    track->setToffset( toffset );
    track->setBeta( m_path1/m_tof12[0]/30.0 );
  }
  else if( layerorend == 21 ) { // outer layer east end readout
    track->setPh( m_ph21 );
    track->setTof( m_tof21[0] );
    track->setSigmaElectron( m_sigma21 );
    for( unsigned int i=0; i<5; i++ ) {
      toffset[i] = m_tof21[0] - m_tof21[i];
    }
    track->setToffset( toffset );
    track->setBeta( m_path2/m_tof21[0]/30.0 );
  }
  else if( layerorend == 22 ) { // outer layer west end readout
    track->setPh( m_ph22 );
    track->setTof( m_tof22[0] );
    track->setSigmaElectron( m_sigma22 );
    for( unsigned int i=0; i<5; i++ ) {
      toffset[i] = m_tof22[0] - m_tof22[i];
    }
    track->setToffset( toffset );
    track->setBeta( m_path2/m_tof22[0]/30.0 );
  }
  else{
    cout << "TofRec TofTrack::SetRecTofTrack  layerorend = " << layerorend << endl;
  }
  return;
}

Referenced by buildRecTofTrack().

◆ setTofData()

void TofTrack::setTofData ( TofDataMap tofDataMap )

Definition at line 288 of file TofTrack.cxx.

                                                 {
 
  if( m_hitCase == NoHit ) return;
 
  unsigned int identify0 = 0x0000c000;
  unsigned int identify1 = 0x0000c000;
  unsigned int identify2 = 0x0000c000;
  unsigned int count0 = 0;
  unsigned int count1 = 0;
  unsigned int count2 = 0;
 
  if( ( ( m_hitCase == InnerLayer ) || ( m_hitCase == DoubleLayer ) ) && ( m_id1 > -1 ) ) {
    int tofid0 = m_id1;
    identify0 = TofID::getIntID( 1, 0, tofid0, 0 );
    count0 = tofDataMap.count( identify0 );
    int tofid1 = tofid0 - 1;
    if( tofid1 == -1 ) { tofid1 = 87; }
    identify1 = TofID::getIntID( 1, 0, tofid1, 0 );
    count1 = tofDataMap.count( identify1 );
    int tofid2 = tofid0 + 1;
    if( tofid2 == 88 ) { tofid2 = 0; }
    identify2 = TofID::getIntID( 1, 0, tofid2, 0 );
    count2 = tofDataMap.count( identify2 );
  }
 
  if( ( m_hitCase == EastEndcap ) && ( m_id1 > -1 ) ) {
    int tofid0 = m_id1;
    identify0 = TofID::getIntID( 0, 0, tofid0, 0 );
    count0 = tofDataMap.count( identify0 );
    int tofid1 = tofid0 - 1;
    if( tofid1 == -1 ) { tofid1 = 47; }
    identify1 = TofID::getIntID( 0, 0, tofid1, 0 );
    count1 = tofDataMap.count( identify1 );
    int tofid2 = tofid0 + 1;
    if( tofid2 == 48 ) { tofid2 = 0; }
    identify2 = TofID::getIntID( 0, 0, tofid2, 0 );
    count2 = tofDataMap.count( identify2 );
  }
 
  if( ( m_hitCase == WestEndcap ) && ( m_id1 > -1 ) ) {
    int tofid0 = m_id1 - 48;
    identify0 = TofID::getIntID( 2, 0, tofid0, 0 );
    count0 = tofDataMap.count( identify0 );
    int tofid1 = tofid0 - 1;
    if( tofid1 == -1 ) { tofid1 = 47; }
    identify1 = TofID::getIntID( 2, 0, tofid1, 0 );
    count1 = tofDataMap.count( identify1 );
    int tofid2 = tofid0 + 1;
    if( tofid2 == 48 ) { tofid2 = 0; }
    identify2 = TofID::getIntID( 2, 0, tofid2, 0 );
    count2 = tofDataMap.count( identify2 );
  }
 
  if( count0 > 0 ) {
    pair< IterTofDataMap, IterTofDataMap > range = tofDataMap.equal_range( identify0 );
    IterTofDataMap iter = range.first;
    for( unsigned int i = 0; i < count0; i++, iter++ ) {
      tofDataAnalysis( (*iter).second, 1 );
    }
  }
  if( count1 > 0 ) {
    pair< IterTofDataMap, IterTofDataMap > range = tofDataMap.equal_range( identify1 );
    IterTofDataMap iter = range.first;
    for( unsigned int i = 0; i < count1; i++, iter++ ) {
      tofDataAnalysis( (*iter).second, 2 );
    }
  }
  if( count2 > 0 ) {
    pair< IterTofDataMap, IterTofDataMap > range = tofDataMap.equal_range( identify2 );
    IterTofDataMap iter = range.first;
    for( unsigned int i = 0; i < count2; i++, iter++ ) {
      tofDataAnalysis( (*iter).second, 2 );
    }
  }
 
  if( ( count0 == 0 ) && ( count1 == 0 ) && ( count2 == 0 ) ) {
    if( m_hitCase != DoubleLayer ) {
      m_hitCase = NoHit; 
      m_quality = 12;
    }
    else {
      m_hitCase = OuterLayer;
    }
  }
 
  if( ( ( m_hitCase == OuterLayer ) || ( m_hitCase == DoubleLayer ) ) && ( m_id2 > 87 ) ) {
    int tofid0 = m_id2 - 88;
    identify0 = TofID::getIntID( 1, 1, tofid0, 0 );
    count0 = tofDataMap.count( identify0 );
    int tofid1 = tofid0 - 1;
    if( tofid1 == -1 ) { tofid1 = 87; }
    identify1 = TofID::getIntID( 1, 1, tofid1, 0 );
    count1 = tofDataMap.count( identify1 );
    int tofid2 = tofid0 + 1;
    if( tofid2 == 88 ) { tofid2 = 0; }
    identify2 = TofID::getIntID( 1, 1, tofid2, 0 );
    count2 = tofDataMap.count( identify2 );
 
    if( count0 > 0 ) {
      pair< IterTofDataMap, IterTofDataMap > range = tofDataMap.equal_range( identify0 );
      IterTofDataMap iter = range.first;
      for( unsigned int i = 0; i < count0; i++, iter++ ) {
    tofDataAnalysis( (*iter).second, 3 );
      }
    }
    if( count1 > 0 ) {
      pair< IterTofDataMap, IterTofDataMap > range = tofDataMap.equal_range( identify1 );
      IterTofDataMap iter = range.first;
      for( unsigned int i = 0; i < count1; i++, iter++ ) {
    tofDataAnalysis( (*iter).second, 4 );
      }
    }
    if( count2 > 0 ) {
      pair< IterTofDataMap, IterTofDataMap > range = tofDataMap.equal_range( identify2 );
      IterTofDataMap iter = range.first;
      for( unsigned int i = 0; i < count2; i++, iter++ ) {
    tofDataAnalysis( (*iter).second, 4 );
      }
    }
 
    if( ( count0 == 0 ) && ( count1 == 0 ) && ( count2 == 0 ) ) {
      if( m_hitCase != DoubleLayer ) {
    m_hitCase = NoHit; 
    m_quality = 12;
      }
      else {
    m_hitCase = InnerLayer;
      }
    }
  }
 
  return;
}

◆ size1()

int TofTrack::size1 ( ) const

inline

Definition at line 46 of file TofTrack.h.

46{ return m_tofData1.size(); }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ size2()

int TofTrack::size2 ( ) const

inline

Definition at line 47 of file TofTrack.h.

47{ return m_tofData2.size(); }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ size3()

int TofTrack::size3 ( ) const

inline

Definition at line 48 of file TofTrack.h.

48{ return m_tofData3.size(); }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ size4()

int TofTrack::size4 ( ) const

inline

Definition at line 49 of file TofTrack.h.

49{ return m_tofData4.size(); }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ t0Stat()

int TofTrack::t0Stat ( ) const

inline

Definition at line 103 of file TofTrack.h.

103{ return m_t0Stat; }

Referenced by setExtTrack().

◆ tdc1()

double TofTrack::tdc1 ( ) const

inline

Definition at line 85 of file TofTrack.h.

85{ return m_tdc1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tdc2()

double TofTrack::tdc2 ( ) const

inline

Definition at line 86 of file TofTrack.h.

86{ return m_tdc2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tdc3()

double TofTrack::tdc3 ( ) const

inline

Definition at line 87 of file TofTrack.h.

87{ return m_tdc3; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tdc4()

double TofTrack::tdc4 ( ) const

inline

Definition at line 88 of file TofTrack.h.

88{ return m_tdc4; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tdiff1()

double TofTrack::tdiff1 ( ) const

inline

Definition at line 100 of file TofTrack.h.

100{ return m_tdiff1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tdiff2()

double TofTrack::tdiff2 ( ) const

inline

Definition at line 101 of file TofTrack.h.

101{ return m_tdiff2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ texp()

double TofTrack::texp ( unsigned int i ) const

inline

Definition at line 92 of file TofTrack.h.

92{ return m_texp[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ texpInner()

double TofTrack::texpInner ( unsigned int i ) const

inline

Definition at line 90 of file TofTrack.h.

90{ return m_texpInner[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ texpOuter()

double TofTrack::texpOuter ( unsigned int i ) const

inline

Definition at line 91 of file TofTrack.h.

91{ return m_texpOuter[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ theta1()

double TofTrack::theta1 ( ) const

inline

Definition at line 50 of file TofTrack.h.

50{ return m_theta1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ theta2()

double TofTrack::theta2 ( ) const

inline

Definition at line 51 of file TofTrack.h.

51{ return m_theta2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tof()

double TofTrack::tof ( unsigned int i ) const

inline

Definition at line 75 of file TofTrack.h.

75{ return m_tof[i]; }

Referenced by compareTofData(), compareTofDataEndcap(), TofCheckDigi::Fill_TofTrack(), findTofDataBarrel(), getTofData(), getTofDataEast(), getTofDataWest(), and tofDataAnalysis().

◆ tof1()

double TofTrack::tof1 ( unsigned int i ) const

inline

Definition at line 73 of file TofTrack.h.

73{ return m_tof1[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tof11()

double TofTrack::tof11 ( unsigned int i ) const

inline

Definition at line 69 of file TofTrack.h.

69{ return m_tof11[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tof12()

double TofTrack::tof12 ( unsigned int i ) const

inline

Definition at line 70 of file TofTrack.h.

70{ return m_tof12[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tof2()

double TofTrack::tof2 ( unsigned int i ) const

inline

Definition at line 74 of file TofTrack.h.

74{ return m_tof2[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tof21()

double TofTrack::tof21 ( unsigned int i ) const

inline

Definition at line 71 of file TofTrack.h.

71{ return m_tof21[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tof22()

double TofTrack::tof22 ( unsigned int i ) const

inline

Definition at line 72 of file TofTrack.h.

72{ return m_tof22[i]; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ tofData1()

std::vector< TofData * > TofTrack::tofData1 ( ) const

inline

Definition at line 44 of file TofTrack.h.

44{ return m_tofData1; }

Referenced by compareTofData(), compareTofDataEndcap(), findTofDataBarrel(), and findTofDataEndcap().

◆ tofData2()

std::vector< TofData * > TofTrack::tofData2 ( ) const

inline

Definition at line 45 of file TofTrack.h.

45{ return m_tofData2; }

Referenced by compareTofData(), compareTofDataEndcap(), findTofDataBarrel(), and findTofDataEndcap().

◆ tofDataAnalysis()

void TofTrack::tofDataAnalysis	(	TofData *	tof,
		unsigned int	iflag
	)

Definition at line 426 of file TofTrack.cxx.

                                                                 {
 
  unsigned int qual = tof->quality();
 
  if( ( qual & 0x10 ) == 0 ) {
 
    //    qual = ( qual | 0x310 );    // zadc, ztdc unmatched, and track matched
    qual = ( qual | 0x10 );    // zadc, ztdc unmatched, and track matched
 
    if( tof->barrel() ) {       // Barrel, Endcap has been done
      if( ( tof->quality() == 0x7 ) || ( tof->quality() == 0xd ) ) {
    qual = ( qual | 0x20 );                       // lost one Q
      }
      
      if( ( tof->quality() == 0xb ) || ( tof->quality() == 0xe ) ) {
    qual = ( qual | 0x40 );                       // lost one T
      }
      
      if( ( tof->quality() == 0x3 ) || ( tof->quality() == 0xc ) ) {
    qual = ( qual | 0x80 );                       // single end
      }
 
      if( ( tof->quality() & 0x5 ) == 0x5 ) {
    double ztdc = tofCaliSvc->ZTDC( tof->tdc1(), tof->tdc2(), tof->tofId() );
    tof->setZTdc( ztdc );
    /*
    if( ( ( (iflag==1) || (iflag==2) ) && ( abs(ztdc-m_zrhit1)<ztdc_Cut ) ) || ( ( (iflag==3) || (iflag==4) ) && ( abs(ztdc-m_zrhit2)<ztdc_Cut ) ) ) { 
      qual = ( qual & 0xeffffeff );               // Z_tdc un-matched.
    }
    else {
      qual = ( qual | 0x100 );
    }
    */
      }
 
      if( ( tof->quality() & 0xa ) == 0xa ) {
    double zadc = tofCaliSvc->ZADC( tof->adc1(), tof->adc2(), tof->tofId() );
    tof->setZAdc( zadc );
    /*
    if( ( ( (iflag==1) || (iflag==2) ) && ( abs(zadc-m_zrhit1)<zadc_Cut ) ) || ( ( (iflag==3) || (iflag==4) ) && ( abs(zadc-m_zrhit2)<zadc_Cut ) ) ) { 
      qual = ( qual & 0xeffffbff );               // Z_adc un-matched.
    }
    else {
      qual = ( qual | 0x200 );
    }
    */
      }
    }
 
    tof->setQuality( qual );
  }
 
  if( iflag == 1 )      m_tofData1.push_back( tof );
  else if( iflag == 2 ) m_tofData2.push_back( tof );
  else if( iflag == 3 ) m_tofData3.push_back( tof );
  else if( iflag == 4 ) m_tofData4.push_back( tof );
  else {
    cout << "TofRec::TofTrack::TofDataAnalylsis: the Flag should be 1-4, out of the Range!" << endl;
  }
 
  return;
}

Referenced by setTofData().

◆ tofDataStudy()

void TofTrack::tofDataStudy ( )

◆ tofId1()

int TofTrack::tofId1 ( ) const

inline

Definition at line 58 of file TofTrack.h.

58{ return m_tofId1; }

Referenced by TofCheckDigi::Fill_TofTrack(), and setExtTrack().

◆ tofId2()

int TofTrack::tofId2 ( ) const

inline

Definition at line 59 of file TofTrack.h.

59{ return m_tofId2; }

Referenced by TofCheckDigi::Fill_TofTrack(), and setExtTrack().

◆ tofTrackId()

int TofTrack::tofTrackId ( ) const

inline

Definition at line 26 of file TofTrack.h.

26{ return m_tofTrackId; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ trackId()

int TofTrack::trackId ( ) const

inline

Definition at line 25 of file TofTrack.h.

25{ return m_trackId; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ zadc1()

double TofTrack::zadc1 ( ) const

inline

Definition at line 96 of file TofTrack.h.

96{ return m_zadc1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ zadc2()

double TofTrack::zadc2 ( ) const

inline

Definition at line 97 of file TofTrack.h.

97{ return m_zadc2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ zrhit1() [1/2]

double TofTrack::zrhit1 ( ) const

inline

Definition at line 35 of file TofTrack.h.

35{ return m_zrhit1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ zrhit1() [2/2]

double TofTrack::zrhit1 ( unsigned int i ) const

inline

Definition at line 40 of file TofTrack.h.

40{ return m_zr1[i]; }

◆ zrhit2() [1/2]

double TofTrack::zrhit2 ( ) const

inline

Definition at line 36 of file TofTrack.h.

36{ return m_zrhit2; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ zrhit2() [2/2]

double TofTrack::zrhit2 ( unsigned int i ) const

inline

Definition at line 41 of file TofTrack.h.

41{ return m_zr2[i]; }

◆ ztdc1()

double TofTrack::ztdc1 ( ) const

inline

Definition at line 94 of file TofTrack.h.

94{ return m_ztdc1; }

Referenced by TofCheckDigi::Fill_TofTrack().

◆ ztdc2()

double TofTrack::ztdc2 ( ) const

inline

Definition at line 95 of file TofTrack.h.

95{ return m_ztdc2; }

Referenced by TofCheckDigi::Fill_TofTrack().

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ TofTrack()

◆ ~TofTrack()

Member Function Documentation

◆ adc1()

◆ adc2()

◆ adc3()

◆ adc4()

◆ barrel()

◆ buildRecTofTrack()

◆ chooseTofData()

◆ chooseTofDataEndcap()

◆ compareTofData()

◆ compareTofDataEndcap()

◆ convert2RecBTofCalHitCol()

◆ convert2RecETofCalHitCol()

◆ convert2RecTofTrackCol()

◆ errzr1() [1/2]

◆ errzr1() [2/2]

◆ errzr2() [1/2]

◆ errzr2() [2/2]

◆ estime()

◆ findTofDataBarrel()

◆ findTofDataEndcap()

◆ flag()

◆ getMultiHit()

◆ getTofData()

◆ getTofDataEast()

◆ getTofDataNohit()

◆ getTofDataWest()

◆ hitCase()

◆ id1()

◆ id2()

◆ isNoHit()

◆ kal()

◆ match()

◆ p()

◆ path()

◆ path1()

◆ path2()

◆ ph()

◆ ph1()

◆ ph11()

◆ ph12()

◆ ph2()

◆ ph21()

◆ ph22()

◆ phi1()

◆ phi2()

◆ qch1()

◆ qch2()

◆ qch3()

◆ qch4()

◆ quality()

◆ quality1()

◆ quality2()

◆ qualityAnalysis()

◆ setCalibration()

◆ setExtTrack()

◆ setFlag()

◆ setQuality()

◆ setQuality1()

◆ setQuality2()

◆ setRecTofTrack()

◆ setTofData()

◆ size1()

◆ size2()

◆ size3()

◆ size4()

◆ t0Stat()

◆ tdc1()

◆ tdc2()

◆ tdc3()

◆ tdc4()

◆ tdiff1()

◆ tdiff2()

◆ texp()

◆ texpInner()