DTBtiChip Class Reference

#include <DTBtiChip.h>

Public Member Functions
void	add_digi (int cell, const DTDigi *digi)
	Add a digi to the DTBtiChip. More...
void	add_digi_clock (int cell, int clock_digi)
	Add a clock digi to the DTBtiChip. More...
void	addTrig (int step, DTBtiTrig *btitrig)
void	clear ()
	delete hits and triggers More...
GlobalPoint	CMSPosition () const
	Position in CMS frame. More...
DTConfigBti *	config () const
	Configuration set. More...
	DTBtiChip (DTBtiCard *card, DTTrigGeom *geom, int supl, int n, DTConfigBti *_config)
	original constructor More...
	DTBtiChip (const DTBtiChip &bti)
	Copy constructor. More...
DTTrigGeom *	geom () const
	Return trigger geometry. More...
std::vector< const DTDigi * >	get_CellDigis (int cell)
	get digi vector - SV 28/XI/02 More...
DTBtiId	id () const
	Return the DTBtiChip Id. More...
void	init_clock ()
LocalPoint	localPosition () const
	Position in chamber frame (x is the one of first traco in slave plane) More...
int	nCellHit () const
	Number of cells with hits. More...
int	nTrig (int step) const
	Number of triggers found. More...
int	number () const
	Return DTBtiChip number. More...
DTBtiChip &	operator= (const DTBtiChip &bti)
	Assignment operator. More...
void	run ()
	Run DTBtiChip algorithm. More...
int	sector () const
	Return sector number. More...
int	station () const
	Return station number. More...
int	superlayer () const
	Return superlayer. More...
DTBtiTrig *	trigger (int step, unsigned n) const
	Return the requested trigger. More...
DTBtiTrigData	triggerData (int step, unsigned n) const
	Return the data part of the requested trigger. More...
std::vector< DTBtiTrig * >	trigList (int step) const
int	wheel () const
	Return wheel number. More...
	~DTBtiChip ()
	Destructor. More...

Private Member Functions
void	acceptMask (BitArray< 80 > *BitArrPtr, int k, int accep)
void	computeEqs ()
void	computeSums ()
int	currentIntStep () const
int	currentStep () const
void	doLTS ()
void	eraseTrigger (int step, unsigned n)
void	findTrig ()
void	init ()
int	keepTrig (const int eq, const int acp, const int code)
int	keepTrigPatt (int flag, const int eq, const int pattType, int hlflag)
bool	matchEq (float eqA, float eqB, int AC)
int	reSum (int a, int b)
int	reSum23 (int a, int b)
void	reSumSet ()
void	setSnap ()
int	store (const int eq, const int code, const int K, const int X, float KeqAB=0., float KeqBC=0., float KeqCD=0., float KeqAC=0., float KeqBD=0., float KeqAD=0.)
void	sum (const int s, const int a, const int b)
void	tick ()

Private Attributes
int	_busyStart_clock [9]
DTBtiCard *	_card
DTConfigBti *	_config
int	_curStep
float	_difs [25]
std::vector< const DTDigi * >	_digis [9]
std::vector< int >	_digis_clock [9]
DTTrigGeom *	_geom
std::vector< DTBtiHit * >	_hits [9]
DTBtiId	_id
float	_JTR [32][3]
float	_Keq [32][6]
float	_KTR [32][2]
int	_MaxKAcc
int	_MaxKcenterTraco
int	_MaxKleftTraco
int	_MaxKrightTraco
int	_MinKAcc
int	_MinKcenterTraco
int	_MinKleftTraco
int	_MinKrightTraco
int	_nStepUsedHits
float	_sums [25]
DTBtiHit *	_thisStepUsedHit [9]
int	_thisStepUsedTimes [9]
std::vector< DTBtiTrig * >	_trigs [DTConfig::NSTEPL-DTConfig::NSTEPF+1]
float	_Xeq [32][2]
float	_XeqAB_patt0
float	_XeqAC_patt0
float	_XeqBD_patt0
float	_XeqCD_patt0
int	init_done
int	RE23
int	RE43
int	reSumAr [3][5]
int	reSumAr23 [3][5]
int	ST
int	ST2
int	ST23
int	ST3
int	ST4
int	ST43
int	ST5
int	ST7

Detailed Description

Implementation of DTBtiChip trigger algorithm Internally uses DTBtiHit to store muon digis

Date:

2010/01/21 10:22:12

Revision:

1.8

Author

S. Vanini

Definition at line 48 of file DTBtiChip.h.

Constructor & Destructor Documentation

DTBtiChip::DTBtiChip	(	DTBtiCard *	card,
		DTTrigGeom *	geom,
		int	supl,
		int	n,
		DTConfigBti *	_config
	)

original constructor

new constructor with configuration

Definition at line 48 of file DTBtiChip.cc.

References _busyStart_clock, _card, _digis, _geom, _hits, _id, _MaxKAcc, _MinKAcc, _trigs, DTTrigGeom::cellH(), DTTrigGeom::cellPitch(), DTConfigBti::CH(), DTConfigBti::CL(), CMSPosition(), DTTrigGeom::CMSPosition(), config(), gather_cfg::cout, DTConfigBti::DEADpar(), debug, i, init_done, DTConfigBti::KAccTheta(), n, DTConfig::NSTEPF, DTConfig::NSTEPL, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), reSumSet(), setSnap(), evf::utils::sid, DTConfigBti::ST(), DTTrigGeom::statId(), DTBtiId::superlayer(), funct::tan(), theta(), DTBtiCard::useAcceptParamFlag(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                         : _card(card), _geom(geom), _config(conf) {


 // original constructor
 setSnap();
 reSumSet(); 

 // Debugging...
  if(config()->debug()>2){
    cout << "DTBtiChip constructor called for BTI number " << n;
    cout << " in superlayer " << supl << endl;
  }

  // reserve the appropriate amount of space for vectors
  int i=0;
  for(i=0;i<DTConfig::NSTEPL - DTConfig::NSTEPF;i++) {
    _trigs[i].reserve(2);
  }

  for(i=0;i<9;i++) {
    _digis[i].reserve(10);
    _hits[i].reserve(10);
  }

  //SV wire dead time init
  int DEAD = config()->DEADpar();
  for(int cell=1; cell<=9; cell++){
    _busyStart_clock[cell-1] = - DEAD -1;
  }

  // Identifier
  DTChamberId sid = _geom->statId();
  _id = DTBtiId(sid, supl, n);

  //if(config()->trigSetupGeom() == 0){
    // set K acceptance in theta view for this BTI: 6 bit resolution....  
    _MinKAcc = 0;
    _MaxKAcc = 63;
    
/*     DTBtiId _id1 = DTBtiId(sid,supl,1);
       
       cout <<"superlayer" << _id.superlayer()<< "BTI1   " <<  _id1.bti()  << " BTICur " << _id.bti()<< endl;
       cout <<endl;
       GlobalPoint gp1 = _geom->CMSPosition(_id1);
       cout << "pos of BTI "<<  _id1.bti()  << gp1 <<endl;
       // K of tracks from vertex
       GlobalPoint gp = CMSPosition();
       cout << "pos of BTI" << _id.bti()  << gp <<endl;    
       cout << endl ; */
       
       
//     // theta bti acceptance cut is in bti chip (no traco in theta!)
//     // acceptance from orca geom: bti theta angle in CMS frame +-2 in K units 
//     if(_id.superlayer()==2){
//       float distp2 = (int)(2*_geom->cellH()*config()->ST()/_geom->cellPitch());
//       float K0 = config()->ST();

// /*      DTBtiId _id1 = DTBtiId(sid,supl,1);
      
//       cout << "BTI1   " <<  _id1.bti() << endl;
//       cout << "BTICur " << _id.bti() <<endl;
//       GlobalPoint gp1 = _geom->CMSPosition(_id1);
//       cout << "pos of BTI 1 " << gp1 <<endl;*/
      
//       // K of tracks from vertex
//       GlobalPoint gp = CMSPosition();
//       if(config()->debug()>3){
//         cout << "Position: R=" << gp.perp() << "cm, Phi=" << gp.phi()*180/3.14159;
//         cout << " deg, Z=" << gp.z() << " cm" << endl;
//       }
//       // CB TEST WITH NEW GEOMETRY
//       // new geometry: modified wrt old due to specularity of theta SLs (still to understand on wheel zero) 19/06/06
//       float theta;
//       if (_id.wheel()==0){ 
//  if(_id.sector()%4>1) theta = atan( gp.z()/gp.perp() );
//  else theta = atan( -gp.z()/gp.perp() );
//       }
//       else theta = atan( fabs(gp.z())/gp.perp() );
//       // .11 =TAN(6.3 deg) ==> k=2 (e' ancora vero? forse questa parte va aggiornata sena ripassare per gli angoli) 19/6/06
//       float thetamin = theta-config()->KAccTheta()*0.055;
//       float thetamax = theta+config()->KAccTheta()*0.055;

//       float fktmin = tan(thetamin)*distp2 + K0;
//       int ktmin = (fktmin>0) ? (int)(fktmin+0.5) : (int)(fktmin-0.5);
//       float fktmax = tan(thetamax)*distp2 + K0;
//       int ktmax = (fktmax>0) ? (int)(fktmax+0.5) : (int)(fktmax-0.5);
// //      float fkbti = -gp.z()/gp.perp()*distp2;
// //      int kbti = (fkbti>0) ? (int)(fkbti+0.5) : (int)(fkbti-0.5);
// //      // K acceptance to point to vertex
// //      int ktmin = kbti-config()->KAccTheta();  // minimum
// //      int ktmax = kbti+config()->KAccTheta();  // maximum
//       if(ktmin>_MinKAcc)_MinKAcc=ktmin;
//       if(ktmax<_MaxKAcc)_MaxKAcc=ktmax;
//     }

//     // debugging
//     if(config()->debug()>2){
//       cout << "CMS position:" << CMSPosition() << endl;
//       cout << "K acceptance:" << _MinKAcc << "," << _MaxKAcc << endl;
//     }
//     // end debugging
// theta bti acceptance cut is in bti chip (no traco in theta!)
    // acceptance is determined about BTI angle wrt vertex with programmable value 
    if(_id.superlayer()==2){
      // 091105 SV theta bti trigger angular acceptance in CMSSW is computed from geometry 
      // (theta of the bti) +- a given tolerance config()->KAccTheta(): tolerance NOT in 
      // hardware configuration. The hw tolerance is given in the system and the 
      // overall acceptance is computed "before" data (CH,CL) is given to the MC
      // or written in the DB. No way to "extract" the tolerance from DB yet.

      if(_card->useAcceptParamFlag()==0){

    //float distp2 = (int)(2*_geom->cellH()*config()->ST()/_geom->cellPitch());   SV fix 17/III/07
        float distp2 = 2*_geom->cellH()*config()->ST()/_geom->cellPitch();
        float K0 = config()->ST();

        // position of BTI 1 and of current one
        DTBtiId _id1 = DTBtiId(sid,supl,1);
        GlobalPoint gp1 = _geom->CMSPosition(_id1); 
        GlobalPoint gp = CMSPosition();
        if(config()->debug()>3){
            cout << "Position: R=" << gp.perp() << "cm, Phi=" << gp.phi()*180/3.14159;
            cout << " deg, Z=" << gp.z() << " cm" << endl;
        }
    // new geometry: modified wrt old due to specularity of theta SLs --> fixed 6/9/06 
        float theta;
        if(gp1.z() < 0.0) 
            theta = atan( -(gp.z())/gp.perp() );                
        else 
            theta = atan( (gp.z())/gp.perp() );

    // set BTI acceptance window : fixed wrt ORCA on 6/9/06  
        float fktmin = tan(theta)*distp2+K0 ;
    int ktmin = static_cast<int>(fktmin)-config()->KAccTheta();
        float fktmax = tan(theta)*distp2+K0+1;
    int ktmax = static_cast<int>(fktmax)+config()->KAccTheta();
        if(ktmin>_MinKAcc)
        _MinKAcc=ktmin;
        if(ktmax<_MaxKAcc)
        _MaxKAcc=ktmax;
      }
      // 091105 SV acceptance is taken simply from CH, CL parameters
      else {
    _MinKAcc = config()->CL();
    _MaxKAcc = config()->CH();
      } 

      // debugging
      if(config()->debug()>2){
        cout << "CMS position:" << CMSPosition() << endl;
        cout << "K acceptance (theta view):" << _MinKAcc << "," << _MaxKAcc  << endl;
      }// end debugging
     
    }//end theta acceptance computation 

  //}// end if trigSetupGeom=0


  //SV flag for initialization....
  init_done = 0; 

}

DTBtiChip::DTBtiChip

(

const DTBtiChip &

bti

)

Copy constructor.

Definition at line 212 of file DTBtiChip.cc.

References _difs, _digis, _hits, _JTR, _Keq, _KTR, _MaxKAcc, _MinKAcc, _sums, _thisStepUsedHit, _thisStepUsedTimes, _trigs, _Xeq, i, j, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, p1, reSumSet(), and setSnap().

                                         :
  _geom(bti._geom), _id(bti._id),
  _curStep(bti._curStep), _nStepUsedHits(bti._nStepUsedHits){

  setSnap();
  reSumSet();

  int i=0;
  for(i=0;i<DTConfig::NSTEPL - DTConfig::NSTEPF;i++) {
    _trigs[i].reserve(2);
    vector<DTBtiTrig*>::const_iterator p;
    for(p=bti._trigs[i].begin();p<bti._trigs[i].end();p++){
      _trigs[i].push_back(*p);
    }
  }
  for(i=0;i<9;i++) {
    _digis[i].reserve(10);
    vector<const DTDigi*>::const_iterator p;
    for(p=bti._digis[i].begin();p<bti._digis[i].end();p++){
      _digis[i].push_back(*p);
    }
    _hits[i].reserve(10);
    vector<DTBtiHit*>::const_iterator p1;
    for(p1=bti._hits[i].begin();p1<bti._hits[i].end();p1++){
      _hits[i].push_back(*p1);
    }
    _thisStepUsedTimes[i] = bti._thisStepUsedTimes[i];
    _thisStepUsedHit[i] = _thisStepUsedHit[i];
  }
  for(i=0;i<25;i++){
    _sums[i] = bti._sums[i];
    _difs[i] = bti._difs[i];
  }
  for(i=0;i<26;i++){
    int j = 0;
    for(j=0;j<6;j++){
      _Keq[i][j] = bti._Keq[i][j];
    }
    for(j=0;j<3;j++){
      _JTR[i][j] = bti._JTR[i][j];
    }
    for(j=0;j<2;j++){
      _Xeq[i][j] = bti._Xeq[i][j];
      _KTR[i][j] = bti._KTR[i][j];
    }
  }
  _MinKAcc = bti._MinKAcc;
  _MaxKAcc = bti._MaxKAcc;



}

DTBtiChip::~DTBtiChip

(

)

Destructor.

Definition at line 267 of file DTBtiChip.cc.

References clear().

                     {
  clear();
}

Member Function Documentation

void DTBtiChip::acceptMask ( BitArray< 80 > *  BitArrPtr, int  k, int  accep  )

private

Definition at line 678 of file DTBtiChipSelTrig.cc.

References HDQMDatabaseProducer::config, gather_cfg::cout, debug, and BitArray< N >::set().

{
   if(k>=0&&k<78){
     if(config()->debug()>4)
       cout<<"DTBtiChip::acceptMask ->  Setting acceptance for k="<<k<<endl;

     if(accep==0){ //code 00
       if(k>=1) 
         BitArrPtr->set(k-1);
       BitArrPtr->set(k);
       BitArrPtr->set(k+1);
     }
     if(accep==1){ //code 01
       BitArrPtr->set(k);
       BitArrPtr->set(k+1);
       BitArrPtr->set(k+2);
     }
     if(accep==2){ //code 10
       if(k>1)
         BitArrPtr->set(k-2);
       if(k>=1)
         BitArrPtr->set(k-1);
       BitArrPtr->set(k);
     }
     if(accep==3){ //code 11
       if(k>1)
         BitArrPtr->set(k-2);
       if(k>=1)
         BitArrPtr->set(k-1);
       BitArrPtr->set(k);
       BitArrPtr->set(k+1);
       BitArrPtr->set(k+2);
     }
   }


  if(config()->debug()>4)
    cout<<"DTBtiChip::acceptMask ->  END "<<endl;

}

void DTBtiChip::add_digi	(	int	cell,
		const DTDigi *	digi
	)

Add a digi to the DTBtiChip.

Definition at line 327 of file DTBtiChip.cc.

References _busyStart_clock, _digis, _geom, _id, DTBtiHit::_stepTimeTdc, DTBtiId::bti(), config(), DTDigi::countsTDC(), gather_cfg::cout, DTConfigBti::DEADpar(), debug, DTTrigGeom::nCell(), DTDigi::print(), and superlayer().

Referenced by DTBtiCard::loadBTI().

                                                {

  if(_id.bti()<1 || _id.bti() >_geom->nCell(superlayer()))return;
  if(cell<1 || cell>9){
    cout << "DTBtiChip::add_digi : wrong cell number: " << cell;
    cout << ". Digi not added!" << endl;
    return;
  }

  int DEAD = config()->DEADpar();
  float stepTimeTdc = DTBtiHit::_stepTimeTdc;


  if( int(digi->countsTDC()/stepTimeTdc) - _busyStart_clock[cell-1] > DEAD ){
    _busyStart_clock[cell-1] = int(digi->countsTDC()/stepTimeTdc);
    _digis[cell-1].push_back(digi);

    // debugging
    if(config()->debug()>1){
    cout << "DTBtiChip::add_digi: DTBtiChip # " <<_id.bti() <<
    " cell " << cell << " --> drift time (tdc units)= " << digi->countsTDC()<< endl;
    digi->print();
    }
  }
  else {
  // debugging
  if(config()->debug()>1)
    cout << "DTBtiChip::add_digi: DTBtiChip # " <<_id.bti() <<
     " cell " << cell << " in dead time -> digi not added! " << endl;
  }

}

void DTBtiChip::add_digi_clock	(	int	cell,
		int	clock_digi
	)

Add a clock digi to the DTBtiChip.

Definition at line 362 of file DTBtiChip.cc.

References _busyStart_clock, _digis_clock, config(), gather_cfg::cout, DTConfigBti::DEADpar(), debug, and number().

                                            {

  if(cell<1 || cell>9){
    cout << "DTBtiChip::add_digi_clock : wrong cell number: " << cell;
    cout << ". Digi not added!" << endl;
    return;
  }

  int DEAD = config()->DEADpar();

  if( digi - _busyStart_clock[cell-1] > DEAD ){
    _busyStart_clock[cell-1] = digi;
    _digis_clock[cell-1].push_back(digi);
    // debugging
    if(config()->debug()>1)
      cout << "DTBtiChip::add_digi_clock: DTBtiChip # " <<number() <<
      " cell " << cell << " --> clock time = " << digi << endl;
  }
  else{
  // debugging
  if(config()->debug()>1)
    cout << "DTBtiChip::add_digi_clock: DTBtiChip # " << number() <<
     " cell " << cell << " in dead time -> digi not added! " << endl;
  }
}

void DTBtiChip::addTrig	(	int	step,
		DTBtiTrig *	btitrig
	)

Add a DTBtiChip trigger (normally used by DTBtiChip itself - may be used for debugging by other classes)

Definition at line 441 of file DTBtiChip.cc.

References _trigs, config(), gather_cfg::cout, debug, DTConfig::NSTEPF, DTConfig::NSTEPL, and relval_parameters_module::step.

Referenced by store(), and DTBtiCard::storeTrigger().

                                               { 
  if(step>=DTConfig::NSTEPF&&step<=DTConfig::NSTEPL){
    if(config()->debug()>3) 
      cout << "DTBtiChip: adding trigger..." <<endl;
    _trigs[step-DTConfig::NSTEPF].push_back(btitrig);
  } else {
    if(config()->debug()>3){    
      cout << "DTBtiChip::addTrig: step " << step ;
      cout << " outside range. Trigger not added" << endl;
    }
  } 
}

void DTBtiChip::clear

(

void

)

delete hits and triggers

Definition at line 522 of file DTBtiChip.cc.

References _digis, _digis_clock, _hits, _trigs, begin, trackerHits::c, config(), gather_cfg::cout, debug, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, and p1.

Referenced by ~DTBtiChip().

                 {

  if(config()->debug()>3)
    cout << "DTBtiChip::clear()" << endl;

  for(int c=0;c<9;c++) {

    _digis[c].clear();
    _digis_clock[c].clear();

    vector<DTBtiHit*>::iterator p;
    for(p=_hits[c].begin();p<_hits[c].end();p++){
      delete (*p);
    }
    _hits[c].clear();
  }

  vector<DTBtiTrig*>::iterator p1;
  for(int is=0;is<DTConfig::NSTEPL-DTConfig::NSTEPF+1;is++){
    for(p1=_trigs[is].begin();p1<_trigs[is].end();p1++){
      delete (*p1);
    }
    _trigs[is].clear();
  }
}

GlobalPoint DTBtiChip::CMSPosition ( ) const

inline

Position in CMS frame.

Definition at line 98 of file DTBtiChip.h.

References _geom, _id, and DTTrigGeom::CMSPosition().

Referenced by DTBtiChip().

{ return _geom->CMSPosition(_id); }

void DTBtiChip::computeEqs ( )

private

Definition at line 145 of file DTBtiChipEquations.cc.

References begin, HDQMDatabaseProducer::config, gather_cfg::cout, i, j, and findQualityFiles::size.

Referenced by run().

                     {
  //
  // Compute all K and X equations of DTBtiChip algorithm
  // NB now Keq=Kman

  float K0 = config()->ST();
  //cout <<"K0="<<K0<<endl;

  //enabled patterns....
  int PTMS[32];
  for(int i=0; i<32; i++){
    PTMS[i] = config()->PTMSflag(i);
  }

    int i;
  // redundant patterns added by Sara Vanini
  i=0;  //  1324A  --> 1L3L2R4L
  if(PTMS[i] ){
  _Keq[i][0] =  _difs[1]          + 2.*K0;                  //eq. AB
  _Keq[i][1] = -_sums[5]          + 2.*K0;                  //eq. BC
  _Keq[i][2] =  _sums[6];                           //eq. CD
  _Keq[i][3] = -(_sums[0]/2.)  + 2.*K0 + 0.01;              //eq. AC
  _Keq[i][4] =  (_difs[8]/2.)  +    K0 + 0.01 ;             //eq. BD
  //_Keq[i][5] =  (_difs[2]/3.) + 4.*K0/3. + 0.51;                 //eq. AD
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(_difs[2]),3.)) ) + 
                int( (double(_difs[2])/3.) );                        //eq. AD

  //patt 0 always uses Xbc, because Xad could be negative...
  //when wire B is missing,  TshiftB=0
  //when hit B is gone out of shift register, Tshift=K0+1
  float _difs_p0 = _difs[5]; 

  float TshiftB = 0;
  float TshiftC = 0;
  if(_thisStepUsedHit[3-1]==0){
    if(_hits[3-1].size()==0 )
      TshiftB = 0;
    if(_hits[3-1].size()>0 && (*(_hits[3-1].begin()))->clockTime()<=-K0 )
      TshiftB = K0+1; 
  }
  else
   TshiftB = _thisStepUsedHit[3-1]->jtrig();  

  if(_thisStepUsedHit[2-1]==0){
    if(_hits[2-1].size()==0 )
      TshiftC = 0;
    if(_hits[2-1].size()>0 && (*(_hits[2-1].begin()))->clockTime()<=-K0 )
      TshiftC = K0+1; 
  }
  else
    TshiftC = _thisStepUsedHit[2-1]->jtrig();  

  _difs_p0 = (float)(TshiftB - TshiftC);

// Zotto's
  _XeqAB_patt0 = (_sums[1] - K0) / 4.;                       //eq. AB
  _Xeq[i][0]   = (_difs_p0 + K0) / 4.;                       //eq. BC
  _XeqCD_patt0 = (_difs[6] + K0) / 4.;                        //eq. CD
  _XeqAC_patt0 = -(_difs[0])   / 4.;                         //eq. AC
  _XeqBD_patt0 =  (_sums[8])   / 4.;                         //eq. BD
  _Xeq[i][1] = (_sums[2] - K0) / 4.;                         //eq. AD
  
/*
// my eq
  _XeqAB_patt0 = (_sums[1] - K0) / 4.;                       //eq. AB
  _Xeq[i][0]   = (_difs[5] + K0) / 4.;                       //eq. BC
  _XeqCD_patt0 = (_difs[6] + K0) / 4.;                        //eq. CD
  _XeqAC_patt0 = -(_difs[0])   / 4.;                         //eq. AC
  _XeqBD_patt0 =  (_sums[8])   / 4.;                         //eq. BD
  _Xeq[i][1] = (_sums[2] - K0) / 4.;                         //eq. AD
*/
   
  }
 

  i=1;  //  1324B  --> 1L3L2R4R
  if(PTMS[i] ){
  _Keq[i][0] =  _difs[1]        + 2.*K0;
  _Keq[i][1] = -(_sums[5])  + 2.*K0;
  _Keq[i][2] = -_difs[6]        + 2.*K0;
  _Keq[i][3] = -(_sums[0]/2.) + 2.*K0 + 0.01;
  _Keq[i][4] = -(_sums[8]/2.) + 2.*K0 + 0.01;
  //_Keq[i][5] = -(_sums[2]/3.)     + 2.*K0 + 0.51;
  _Keq[i][5] = ST2 +
               reSum23( 0 , int(fmod(double(-_sums[2]),3.)) ) +
               int( (double(-_sums[2])/3.) );   

  _Xeq[i][0] = ( _difs[5] + K0) / 4.;
  _Xeq[i][1] = (-_difs[2] + K0) / 4.;
  }
 
  i=2;  //  1324C  --> 1R3L2R4L
  if(PTMS[i] ){
  _Keq[i][0] =  _sums[1];
  _Keq[i][1] = -_sums[5]       + 2.*K0;
  _Keq[i][2] =  _sums[6];
  _Keq[i][3] = -(_difs[0]/2.) +   K0 + 0.01;
  _Keq[i][4] =  (_difs[8]/2.) +   K0 + 0.01;
  //_Keq[i][5] =  (_sums[2]/3.)    + 2.*K0/3. + 0.51;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(_sums[2]),3.)) ) + 
                int( (double(_sums[2])/3.) );                        //eq. AD
/*  
cout << "ST23 =" << ST23 << endl;
cout << "RE23 =" << RE23 << endl;
cout << "fmod(double(_sums[2]),3.) =" << fmod(double(_sums[2]),3.) << endl;
cout << "reSum23 = "<< reSum23( RE23 , fmod(double(_sums[2]),3.) ) << endl;
cout << "double(_sums[2])/3.="<<double(_sums[2])/3.<< endl;
cout << "int('') = " << int( (double(_sums[2])/3.) ) << endl;
*/
  _Xeq[i][0] = (_difs[5] + K0) / 4.;
  _Xeq[i][1] = (_difs[2] + K0) / 4.;
  }
 
  i=3;  //  1324D  --> 1R3L2R4R
  if(PTMS[i] ){
  _Keq[i][0] =  _sums[1];
  _Keq[i][1] = -_sums[5]      + 2.*K0;
  _Keq[i][2] = -_difs[6]      + 2.*K0;
  _Keq[i][3] = -(_difs[0]/2.) +    K0 + 0.01;
  _Keq[i][4] = -(_sums[8]/2.) + 2.*K0 + 0.01;
  //_Keq[i][5] = -(_difs[2]/3.) + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(-_difs[2]),3.)) ) + 
                int( (double(-_difs[2])/3.) );                        //eq. AD


  _Xeq[i][0] = ( _difs[5] +   K0) / 4.;
  _Xeq[i][1] = (-_sums[2] + 3.*K0) / 4.;
  }
 

  i=4;  //  i = DTBtiChipEQMAP->index("1364A");  --> 1L3L6L4R
  if(PTMS[i] ){
  _Keq[i][0] =  _difs[1]        + 2.*K0;
  _Keq[i][1] =  _difs[10]       + 2.*K0;
  _Keq[i][2] = -(_sums[13]) + 2.*K0;
  _Keq[i][3] =  (_difs[3]/2.)   + 2.*K0 + 0.01;
  _Keq[i][4] = -(_sums[8]/2.)   + 2.*K0 + 0.01;
  //_Keq[i][5] = -(_sums[2]/3.)   + 2.*K0 + 0.51;
  _Keq[i][5] = ST2 +
               reSum23( 0 , int(fmod(double(-_sums[2]),3.)) ) +
               int( (double(-_sums[2])/3.) );   

  _Xeq[i][0] = ( _sums[10] + K0) / 4.;
  _Xeq[i][1] = (-_difs[2]  + K0) / 4.;
  }
 
  i=5;  //  i = DTBtiChipEQMAP->index("1364B");  --> 1R3L6L4R
  if(PTMS[i] ){
  _Keq[i][0] =  _sums[1];
  _Keq[i][1] =  _difs[10]      + 2.*K0;
  _Keq[i][2] = -_sums[13]      + 2.*K0;
  _Keq[i][3] =  (_sums[3]/2.)  +   K0 + 0.01;
  _Keq[i][4] = -(_sums[8]/2.)  + 2.*K0 + 0.01;
  //_Keq[i][5] = -(_difs[2]/3.)  + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(-_difs[2]),3.)) ) + 
                int( (double(-_difs[2])/3.) );                        //eq. AD

  _Xeq[i][0] = ( _sums[10]+   K0) / 4.;
  _Xeq[i][1] = (-_sums[2] + 3.*K0) / 4.;
  }

  i=6;  //  i = DTBtiChipEQMAP->index("1364C");  --> 1R3R6L4R
  if(PTMS[i] ){
  _Keq[i][0] = -_difs[1]       + 2.*K0;
  _Keq[i][1] =  _sums[10];
  _Keq[i][2] = -_sums[13]      + 2.*K0;
  _Keq[i][3] =  (_sums[3]/2.)  +    K0 + 0.01;
  _Keq[i][4] = -(_difs[8]/2.)  +    K0 + 0.01;
  //_Keq[i][5] = -(_difs[2]/3.)  + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(-_difs[2]),3.)) ) + 
                int( (double(-_difs[2])/3.) );                        //eq. AD

  _Xeq[i][0] = ( _difs[10]+ 3.*K0) / 4.;
  _Xeq[i][1] = (-_sums[2] + 3.*K0) / 4.;
  }
 
  i=7;  //  i = DTBtiChipEQMAP->index("1368A");  --> 1R3R6L8L
  if(PTMS[i] ){
  _Keq[i][0] = -_difs[1]        + 2.*K0;
  _Keq[i][1] =  (_sums[10]);
  _Keq[i][2] =  _difs[20]       + 2.*K0;
  _Keq[i][3] =  (_sums[3]/2.)   +    K0 + 0.01;
  _Keq[i][4] =  (_sums[11]/2.)  +    K0 + 0.01;
  //_Keq[i][5] =  (_sums[4]/3.)   + 4.*K0/3. + 0.51;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(_sums[4]),3.)) ) + 
                int( (double(_sums[4])/3.) );                        //eq. AD

  _Xeq[i][0] =  (_difs[10]+ 3.*K0) / 4.;
  _Xeq[i][1] =  (_difs[4] + 3.*K0) / 4.;
  }
 
  i=8;  //  i = DTBtiChipEQMAP->index("1368B");  --> 1R3R6R8L
  if(PTMS[i] ){
  _Keq[i][0] = -_difs[1]          + 2.*K0;
  _Keq[i][1] = -_difs[10]         + 2.*K0;
  _Keq[i][2] =  (_sums[20]);
  _Keq[i][3] = -(_difs[3]/2.)     + 2.*K0 + 0.01;
  _Keq[i][4] =  (_sums[11]/2.)    +    K0 + 0.01;
  //_Keq[i][5] =  (_sums[4]/3.)     + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(_sums[4]),3.)) ) + 
                int( (double(_sums[4])/3.) );                        //eq. AD


  _Xeq[i][0] = (-_sums[10]+ 5.*K0) / 4.;
  _Xeq[i][1] = ( _difs[4] + 3.*K0) / 4.;
  }

  i=9;  //  i = DTBtiChipEQMAP->index("1368C");  --> 1R3L6L8L
  if(PTMS[i] ){
  _Keq[i][0] =  (_sums[1]);
  _Keq[i][1] =  _difs[10]        + 2.*K0;
  _Keq[i][2] =  _difs[20]        + 2.*K0;
  _Keq[i][3] =  (_sums[3]/2.)    +    K0 + 0.01;
  _Keq[i][4] =  (_difs[11]/2.)   + 2.*K0 + 0.01;
  //_Keq[i][5] =  (_sums[4]/3.)    + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(_sums[4]),3.)) ) + 
                int( (double(_sums[4])/3.) );                        //eq. AD


  _Xeq[i][0] =  (_sums[10]+  K0) / 4.;
  _Xeq[i][1] =  (_difs[4] + 3.*K0) / 4.;
  }
  
  i=10;  //  i = DTBtiChipEQMAP->index("5324A");  --> 5L3L2R4L
  if(PTMS[i] ){
  _Keq[i][0] = -_difs[9];
  _Keq[i][1] = -_sums[5]       + 2.*K0;
  _Keq[i][2] =  _sums[6];
  _Keq[i][3] = -(_sums[7]/2.)  +    K0 + 0.01;
  _Keq[i][4] =  (_difs[8]/2.)  +    K0 + 0.01;
  //_Keq[i][5] = -(_difs[12]/3.) + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(-_difs[12]),3.)) ) + 
                int( (double(-_difs[12])/3.) );                        //eq. AD

  _Xeq[i][0] =  (_difs[5] + K0) / 4.;
  _Xeq[i][1] =  (_sums[12]+ K0) / 4.;
  }

  i=11;  //  i = DTBtiChipEQMAP->index("5324B");  --> 5L3R2R4L
  if(PTMS[i] ){
  _Keq[i][0] = -_sums[9]       + 2.*K0;
  _Keq[i][1] =  _difs[5];
  _Keq[i][2] =  _sums[6];
  _Keq[i][3] = -(_sums[7]/2.)  +    K0 + 0.01;
  _Keq[i][4] =  (_sums[8]/2.   + 0.01);
  //_Keq[i][5] = -(_difs[12]/3.) + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(-_difs[12]),3.)) ) + 
                int( (double(-_difs[12])/3.) );                        //eq. AD

  _Xeq[i][0] = (-_sums[5] + 3.*K0) / 4.;
  _Xeq[i][1] = ( _sums[12]+   K0) / 4.;
  }

  i=12;  //  i = DTBtiChipEQMAP->index("5324C");  --> 5R3R2R4L
  if(PTMS[i] ){
  _Keq[i][0] =  _difs[9];
  _Keq[i][1] =  _difs[5];
  _Keq[i][2] =  (_sums[6]);
  _Keq[i][3] =  (_difs[7]/2.) + 0.01;
  _Keq[i][4] =  (_sums[8]/2.)  + 0.01;
  _Keq[i][5] =  (_sums[12]/3.) + 0.51;

  _Xeq[i][0] = (-_sums[5] + 3. * K0) / 4.;
  _Xeq[i][1] = (-_difs[12]+ 3. * K0) / 4.;
  }

  i=13;  //  i = DTBtiChipEQMAP->index("5364A");  --> 5L3R6L4L
  if(PTMS[i] ){
  _Keq[i][0] = -_sums[9]        + 2.*K0;
  _Keq[i][1] =  _sums[10];
  _Keq[i][2] = -_difs[13];
  _Keq[i][3] =  (_difs[16]/2.)  +    K0 + 0.01;
  _Keq[i][4] =  (_sums[8]/2.    + 0.01);
  //_Keq[i][5] = -(_difs[12]/3.)  +  2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(-_difs[12]),3.)) ) + 
                int( (double(-_difs[12])/3.) );                        //eq. AD

  _Xeq[i][0] =  (_difs[10]+ 3.*K0) / 4.;
  _Xeq[i][1] =  (_sums[12]+   K0) / 4.;
  }
 
  i=14;  //  i = DTBtiChipEQMAP->index("5364B");  --> 5L3R6L4R
  if(PTMS[i] ){
  _Keq[i][0] = -_sums[9]        + 2.*K0;
  _Keq[i][1] =  _sums[10];
  _Keq[i][2] = -_sums[13]       + 2.*K0;
  _Keq[i][3] =  (_difs[16]/2.)  +    K0 + 0.01;
  _Keq[i][4] = -(_difs[8]/2.)   +    K0 + 0.01;
  //_Keq[i][5] = -(_sums[12]/3.)  + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(-_sums[12]),3.)) ) + 
                int( (double(-_sums[12])/3.) ); 

  _Xeq[i][0] =  (_difs[10]+ 3.*K0) / 4.;
  _Xeq[i][1] =  (_difs[12]+ 3.*K0) / 4.;
  }
 
  i=15;  //  i = DTBtiChipEQMAP->index("5364C");  --> 5R3R6L4L
  if(PTMS[i] ){
  _Keq[i][0] =  _difs[9];
  _Keq[i][1] =  (_sums[10]);
  _Keq[i][2] = -_difs[13];
  _Keq[i][3] =  (_sums[16]/2.  + 0.01);
  _Keq[i][4] =  (_sums[8]/2.   + 0.01);
  _Keq[i][5] =  (_sums[12]/3.) + 0.51;

  _Xeq[i][0] = ( _difs[10]+ 3.*K0) / 4.;
  _Xeq[i][1] = (-_difs[12]+ 3.*K0) / 4.;
  }
 
  i=16;  //  i = DTBtiChipEQMAP->index("5364D");  --> 5R3R6L4R
  if(PTMS[i] ){
  _Keq[i][0] =  _difs[9];
  _Keq[i][1] =  _sums[10];
  _Keq[i][2] = -_sums[13]      + 2.*K0;
  _Keq[i][3] =  (_sums[16]/2.  + 0.01);
  _Keq[i][4] = -(_difs[8]/2.)  +    K0 + 0.01;
  //_Keq[i][5] =  (_difs[12]/3.) + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(_difs[12]),3.)) ) + 
                int( (double(_difs[12])/3.) );                        //eq. AD

  _Xeq[i][0] = ( _difs[10]+ 3.*K0) / 4.;
  _Xeq[i][1] = (-_sums[12]+ 5.*K0) / 4.;
  }
 
  i=17;  //  i = DTBtiChipEQMAP->index("5368A");  --> 5L3R6L8L
  if(PTMS[i] ){
  _Keq[i][0] = -_sums[9]        + 2.*K0;
  _Keq[i][1] =  _sums[10];
  _Keq[i][2] =  _difs[20]       + 2.*K0;
  _Keq[i][3] =  (_difs[16]/2.)  +    K0 + 0.01;
  _Keq[i][4] =  (_sums[11]/2.)  +    K0 + 0.01;
  //_Keq[i][5] =  (_difs[18]/3.)  + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(_difs[18]),3.)) ) + 
                int( (double(_difs[18])/3.) );                        //eq. AD

  _Xeq[i][0] =  (_difs[10]+ 3.*K0) / 4.;
  _Xeq[i][1] =  (_sums[18]+ 3.*K0) / 4.;
  }
 
  i=18;  //  i = DTBtiChipEQMAP->index("5368B");  --> 5L3R6R8L
  if(PTMS[i] ){
  _Keq[i][0] = -_sums[9]        + 2.*K0;
  _Keq[i][1] = -_difs[10]       + 2.*K0;
  _Keq[i][2] =  _sums[20];
  _Keq[i][3] = -(_sums[16]/2.)  + 2.*K0 + 0.01;
  _Keq[i][4] =  (_sums[11]/2.)  +    K0 + 0.01;
  //_Keq[i][5] =  (_difs[18]/3.)  + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(_difs[18]),3.)) ) + 
                int( (double(_difs[18])/3.) );                        //eq. AD

  _Xeq[i][0] = (-_sums[10]+ 5.*K0) / 4.;
  _Xeq[i][1] = ( _sums[18]+ 3.*K0) / 4.;
  }
 
  i=19;  //  i = DTBtiChipEQMAP->index("5368C");  --> 5L3R6R8R
  if(PTMS[i] ){
  _Keq[i][0] = -(_sums[9])  + 2.*K0;
  _Keq[i][1] = -_difs[10]       + 2.*K0;
  _Keq[i][2] = -_difs[20]       + 2.*K0;
  _Keq[i][3] = -(_sums[16]/2.)  + 2.*K0 + 0.01;
  _Keq[i][4] = -(_difs[11]/2.)  + 2.*K0 + 0.01;
  //_Keq[i][5] = -(_sums[18]/3.)  + 2.*K0 + 0.51;
  _Keq[i][5] = ST2 +
               reSum23( 0 , int(fmod(double(-_sums[18]),3.)) ) +
               int( (double(-_sums[18])/3.) );   

  _Xeq[i][0] = (-_sums[10]+ 5.*K0) / 4.;
  _Xeq[i][1] = (-_difs[18]+ 5.*K0) / 4.;
  }
 
  i=20;  //  i = DTBtiChipEQMAP->index("5764A");  --> 5R7L6L4R
  if(PTMS[i] ){
  _Keq[i][0] =  _sums[17];
  _Keq[i][1] = -_difs[19];
  _Keq[i][2] = -_sums[13]       + 2.*K0;
  _Keq[i][3] =  (_sums[16]/2.   + 0.01);
  _Keq[i][4] = -(_sums[14]/2.)  +    K0 + 0.01;
  //_Keq[i][5] =  (_difs[12]/3.)  + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(_difs[12]),3.)) ) + 
                int( (double(_difs[12])/3.) );                        //eq. AD

  _Xeq[i][0] = ( _sums[19]+ 3.*K0) / 4.;
  _Xeq[i][1] = (-_sums[12]+ 5.*K0) / 4.;
  }
 
  i=21;  //  i = DTBtiChipEQMAP->index("5764B");  --> 5R7L6R4R
  if(PTMS[i] ){
  _Keq[i][0] =  _sums[17];
  _Keq[i][1] = -_sums[19]      + 2.*K0;
  _Keq[i][2] =  _difs[13];
  _Keq[i][3] = -(_difs[16]/2.) +    K0 + 0.01;
  _Keq[i][4] = -(_sums[14]/2.) +    K0 + 0.01;
  //_Keq[i][5] =  (_difs[12]/3.) + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(_difs[12]),3.)) ) + 
                int( (double(_difs[12])/3.) );                        //eq. AD

  _Xeq[i][0] = ( _difs[19]+ 5.*K0) / 4.;
  _Xeq[i][1] = (-_sums[12]+ 5.*K0) / 4.;
  }
 
  i=22;  //  i = DTBtiChipEQMAP->index("5764C");  --> 5R7L6L4L
  if(PTMS[i] ){
  _Keq[i][0] =  (_sums[17]);
  _Keq[i][1] = -_difs[19];
  _Keq[i][2] = -_difs[13];
  _Keq[i][3] =  (_sums[16]/2.  + 0.01);
  _Keq[i][4] = -(_difs[14]/2.) + 0.01;
  _Keq[i][5] =  (_sums[12]/3.) + 0.51;

  _Xeq[i][0] = ( _sums[19]+ 3.*K0) / 4.;
  _Xeq[i][1] = (-_difs[12]+ 3.*K0) / 4.;
  }
 
  i=23;  //  i = DTBtiChipEQMAP->index("9764A");  --> 9L7L6L4R
  if(PTMS[i] ){
  _Keq[i][0] = -_difs[23];
  _Keq[i][1] = -_difs[19];
  _Keq[i][2] = -(_sums[13]) + 2.*K0;
  _Keq[i][3] = -(_difs[21]/2.) + 0.01;
  _Keq[i][4] = -(_sums[14]/2.)  +    K0 + 0.01;
  //_Keq[i][5] = -(_sums[15]/3.)  + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(-_sums[15]),3.)) ) + 
                int( (double(-_sums[15])/3.) );                        //eq. AD

  _Xeq[i][0] =  (_sums[19]+ 3.*K0) / 4.;
  _Xeq[i][1] =  (_difs[15]+ 5.*K0) / 4.;
  }
 
  i=24;  //  i = DTBtiChipEQMAP->index("9764B");  --> 9L7L6R4R
  if(PTMS[i] ){
  _Keq[i][0] = -_difs[23];
  _Keq[i][1] = -(_sums[19])   + 2.*K0;
  _Keq[i][2] =  _difs[13];
  _Keq[i][3] = -(_sums[21]/2.)    +    K0 + 0.01;
  _Keq[i][4] = -(_sums[14]/2.)    +    K0 + 0.01;
  //_Keq[i][5] = -(_sums[15]/3.)    + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(-_sums[15]),3.)) ) + 
                int( (double(-_sums[15])/3.) );                        //eq. AD

  _Xeq[i][0] =  (_difs[19]+ 5.*K0) / 4.;
  _Xeq[i][1] =  (_difs[15]+ 5.*K0) / 4.;
  }
 
  i=25;  //  i = DTBtiChipEQMAP->index("9764C");  --> 9L7R6R4R
  if(PTMS[i] ){
  _Keq[i][0] = -(_sums[23])   + 2.*K0;
  _Keq[i][1] =  _difs[19];
  _Keq[i][2] =  _difs[13];
  _Keq[i][3] = -(_sums[21]/2.)    +   K0 + 0.01;
  _Keq[i][4] =  (_difs[14]/2.) + 0.01;
  //_Keq[i][5] = -(_sums[15]/3.)    + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(-_sums[15]),3.)) ) + 
                int( (double(-_sums[15])/3.) );                        //eq. AD

  _Xeq[i][0] = (-_sums[19]+ 7.*K0) / 4.;
  _Xeq[i][1] = ( _difs[15]+ 5.*K0) / 4.;
  }

  i=26;  //  int i = DTBtiChipEQMAP->index("5768A") --> 5L7L6R8L
  if(PTMS[i] ){
  _Keq[i][0] =  _difs[17]       + 2.*K0;
  _Keq[i][1] = -_sums[19]       + 2.*K0;
  _Keq[i][2] =  _sums[20];
  _Keq[i][3] = -(_sums[16]/2.)  + 2.*K0 + 0.01;
  _Keq[i][4] =  (_difs[22]/2.)  +    K0 + 0.01;
  //_Keq[i][5] =  (_difs[18]/3.)  + 4.*K0/3.;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(_difs[18]),3.)) ) + 
                int( (double(_difs[18])/3.) );                        //eq. AD

  _Xeq[i][0] =  (_difs[19] + 5.*K0) / 4.;
  _Xeq[i][1] =  (_sums[18] + 3.*K0) / 4.;
  }
 
  i=27;  //  i = DTBtiChipEQMAP->index("5768B");  --> 5L7L6R8R
  if(PTMS[i] ){
  _Keq[i][0] =  _difs[17]       + 2.*K0;
  _Keq[i][1] = -(_sums[19]) + 2.*K0;
  _Keq[i][2] = -_difs[20]       + 2.*K0;
  _Keq[i][3] = -(_sums[16]/2.)  + 2.*K0 + 0.01;
  _Keq[i][4] = -(_sums[22]/2.)  + 2.*K0 + 0.01;
  //_Keq[i][5] = -(_sums[18]/3.)  + 2.*K0 + 0.51;
  _Keq[i][5] = ST2 +
               reSum23( 0 , int(fmod(double(-_sums[18]),3.)) ) +
               int( (double(-_sums[18])/3.) );   
  
  _Xeq[i][0] = ( _difs[19] + 5.*K0) / 4.;
  _Xeq[i][1] = (-_difs[18] + 5.*K0) / 4.;
  }
 
  i=28;  //  i = DTBtiChipEQMAP->index("5768C");  --> 5R7L6R8L
  if(PTMS[i] ){
  _Keq[i][0] =  _sums[17];
  _Keq[i][1] = -_sums[19]       + 2.*K0;
  _Keq[i][2] =  _sums[20];
  _Keq[i][3] = -(_difs[16]/2.)  +   K0 + 0.01;
  _Keq[i][4] =  (_difs[22]/2.)  +   K0 + 0.01;
  //_Keq[i][5] =  (_sums[18]/3.)  + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(_sums[18]),3.)) ) + 
                int( (double(_sums[18])/3.) );                        //eq. AD

  _Xeq[i][0] =  (_difs[19] + 5.*K0) / 4.;
  _Xeq[i][1] =  (_difs[18] + 5.*K0) / 4.;
  }

  i=29;  //  i = DTBtiChipEQMAP->index("5768D");  --> 5R7L6R8R
  if(PTMS[i] ){
  _Keq[i][0] =  _sums[17];
  _Keq[i][1] = -_sums[19]       + 2.*K0;
  _Keq[i][2] = -_difs[20]       + 2.*K0;
  _Keq[i][3] = -(_difs[16]/2.)  +    K0 + 0.01;
  _Keq[i][4] = -(_sums[22]/2.)  + 2.*K0 + 0.01;
  //_Keq[i][5] = -(_difs[18]/3.)  + 4.*K0/3. ;
  _Keq[i][5] =  ST43 + 
                reSum( RE43 , int(fmod(double(-_difs[18]),3.)) ) + 
                int( (double(-_difs[18])/3.) );                        //eq. AD

  _Xeq[i][0] = ( _difs[19] + 5.*K0) / 4.;
  _Xeq[i][1] = (-_sums[18] + 7.*K0) / 4.;
  }

  i=30;  //  9768A  --> 9L7L6R8L
  if(PTMS[i] ){
  _Keq[i][0] = -_difs[23];
  _Keq[i][1] = -_sums[19]       + 2.*K0;
  _Keq[i][2] =  _sums[20];
  _Keq[i][3] = -(_sums[21]/2.)  +    K0 + 0.01;
  _Keq[i][4] =  (_difs[22]/2.)  +    K0 + 0.01;
  //_Keq[i][5] = -(_difs[24]/3.)  + 2.*K0/3. ;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(-_difs[24]),3.)) ) + 
                int( (double(-_difs[24])/3.) );                        //eq. AD

  _Xeq[i][0] = (_difs[19] + 5.*K0) / 4.;
  _Xeq[i][1] = (_sums[24] + 5.*K0) / 4.;
  }
 
  i=31;  //  9768B  --> 9L7R6R8L
  if(PTMS[i] ){
  _Keq[i][0] = -_sums[23]       + 2.*K0;
  _Keq[i][1] =  _difs[19];
  _Keq[i][2] =  _sums[20];
  _Keq[i][3] = -(_sums[21]/2.)  +   K0 + 0.01;
  _Keq[i][4] =  (_sums[22]/2.) + 0.01;
  //_Keq[i][5] = -(_difs[24]/3.)  + 2.*K0/3.;
  _Keq[i][5] =  ST23 + 
                reSum23( RE23 , int(fmod(double(-_difs[24]),3.)) ) + 
                int( (double(-_difs[24])/3.) );                        //eq. AD

  _Xeq[i][0] = (-_sums[19] + 7.*K0) / 4.;
  _Xeq[i][1] = ( _sums[24] + 5.*K0) / 4.;
  }
 
  // debugging
  if(config()->debug()>3){
    cout << endl << " Step: " << currentStep() << endl;
    for(i=0;i<32;i++){
      if(PTMS[i] ){
        cout << "K Equation " << i << " --> ";
        int j=0;
        for(j=0;j<6;j++){
      cout << _Keq[i][j] << " ";
        }
        cout << endl;
        cout << "X Equation " << i << " --> ";
        for(j=0;j<2;j++){
      cout << _Xeq[i][j] << " ";
        }
        if( i==0 ){
          cout << _XeqAB_patt0 << " " << _XeqCD_patt0 << " ";
          cout << _XeqAC_patt0 << " " << _XeqBD_patt0 << " ";
        }
        cout << endl;
      }
    }
  }
  // end debugging
  
}

void DTBtiChip::computeSums ( )

private

Definition at line 39 of file DTBtiChipEquations.cc.

References HDQMDatabaseProducer::config, gather_cfg::cout, and debug.

Referenced by run().

                      {
  //
  // compute all sums and diffs
  // nothing more than a table...
  // numbering convention here is the same as the fortran version:
  //   sum number in the range [1,25]
  //   cell numbers in the range [1,9]
  //     --> sum(int,int,int) decreases the indices by 1
  // sum (sum_number, first_cell, second_cell)
  //
  
  if(config()->debug()>3){
    cout << "DTBtiChip::computeSums called" << endl; 
  }

  sum( 1,2,1);
  sum( 2,3,1);
  sum( 3,4,1);
  sum( 4,6,1);
  sum( 5,8,1);
  sum( 6,3,2);
  sum( 7,4,2);
  sum( 8,5,2);
  sum( 9,4,3);
  sum(10,5,3);
  sum(11,6,3);
  sum(12,8,3);
  sum(13,5,4);
  sum(14,6,4);
  sum(15,7,4);
  sum(16,9,4);
  sum(17,6,5);
  sum(18,7,5);
  sum(19,8,5);
  sum(20,7,6);
  sum(21,8,6);
  sum(22,9,6);
  sum(23,8,7);
  sum(24,9,7);
  sum(25,9,8);
}

DTConfigBti* DTBtiChip::config ( void  ) const

inline

Configuration set.

testing DTConfigBti

Definition at line 119 of file DTBtiChip.h.

References _config.

Referenced by add_digi(), add_digi_clock(), addTrig(), clear(), doLTS(), DTBtiChip(), init(), init_clock(), nCellHit(), run(), setSnap(), store(), and tick().

{ return _config; }

int DTBtiChip::currentIntStep ( ) const

inlineprivate

Definition at line 147 of file DTBtiChip.h.

References _curStep.

Referenced by run(), and tick().

{ return _curStep; } // int. step (80MHz)

int DTBtiChip::currentStep ( ) const

inlineprivate

Definition at line 146 of file DTBtiChip.h.

References _curStep.

Referenced by run(), and store().

{ return (int)(((float)(_curStep)+0.5)/2); }

void DTBtiChip::doLTS ( )

private

Definition at line 741 of file DTBtiChip.cc.

References DTBtiTrig::code(), config(), gather_cfg::cout, debug, DTConfigBti::debug(), eraseTrigger(), DTConfigBti::LTS(), DTConfig::NSTEPF, DTConfig::NSTEPL, nTrig(), DTConfigBti::SET(), and trigger().

Referenced by run().

                 {
 
  if(config()->debug()>2)
    cout<<"Do LTS"<<endl;
  int lts = config()->LTS();
  int nbxlts = config()->SET();

  // Do LTS only on the requested SL
  //if (superlayer()==2 && lts==1) return;
  //if (superlayer()!=2 && lts==2) return;
  //new DTConfig: do LTS only is LTS!=0  --> somewhat redundant !
  if (lts==0) return;

  // loop on steps
  for(int is=DTConfig::NSTEPF; is<=DTConfig::NSTEPL; is++) {
    if(nTrig(is)>0) { // non empty step
      if( trigger(is,1)->code()==8 ) { // HTRIG at this step
        // do LTS on nbxLTS[superlayer] following steps
        for(int js=is+1;(js<=is+nbxlts&&js<=DTConfig::NSTEPL);js++){
          if(nTrig(js)>0) { // non empty step
            DTBtiTrig* tr = trigger(js,1);
            if( tr->code()<8 && (lts==1 || lts==3)) {
              if(config()->debug()>3)
                cout<<"LTS: erasing trigger!"<<endl; 
              eraseTrigger(js,1); // delete trigger
            }
          }
        }
        // do LTS on previous step
        if(is>DTConfig::NSTEPF && nTrig(is-1)>0) { // non empty step
          DTBtiTrig* tr = trigger(is-1,1);
          if( tr->code()<8 && (lts==2 || lts==3) ) {
            if(config()->debug()>3)
                cout<<"LTS: erasing trigger!"<<endl;                                
            eraseTrigger(is-1,1); // delete trigger
          }
        }

      }
    }
  }
}

void DTBtiChip::eraseTrigger ( int  step, unsigned  n  )

private

Definition at line 507 of file DTBtiChip.cc.

References _trigs, gather_cfg::cout, n, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, findQualityFiles::size, and relval_parameters_module::step.

Referenced by doLTS().

                                            {
  if(step<DTConfig::NSTEPF||step>DTConfig::NSTEPL){
    cout << "DTBtiChip::eraseTrigger: step out of range: " << step ;
    cout << " trigger not deleted!" << endl;
  } 
  if(n<1 || n>_trigs[step-DTConfig::NSTEPF].size()) {
    cout << "DTBtiChip::trigger: requested trigger does not exist: " << n;
    cout << " trigger not deleted!" << endl;
  }
  vector<DTBtiTrig*>::iterator p = _trigs[step-DTConfig::NSTEPF].begin()+n-1;
  if(&(*p))delete (*p);
  _trigs[step-DTConfig::NSTEPF].erase(p);
}

void DTBtiChip::findTrig ( )

private

Definition at line 40 of file DTBtiChipSelTrig.cc.

References HDQMDatabaseProducer::config, gather_cfg::cout, debug, and i.

Referenced by run().

                   {

  if(config()->debug()>3){
    cout << "DTBtiChip::findTrig called" << endl; 
  }

  //pattern type:       1 = normal
  //            2 = 1L1
  //            3 = 11L
  //            4 = L11
  //            5 = 1H1
  //            6 = 11H
  //            7 = H11

//  int accpatB = config()->AccPattB(); //default +-1
//  int accpatA = config()->AccPattA(); //default +-2
//  int tiKes = config()->XON();
  int RON = config()->RONflag();  //default 1, redundant patterns enabled
  int PTMS[32];
  for(int i=0; i<32; i++){
    PTMS[i] = config()->PTMSflag(i);
  }

  //Triggers (same order as manual):
  for(int hl=0;hl<2;hl++){  //hl=0 high   hl=1 low
  if( RON==1 ){
    if( keepTrigPatt(PTMS[0],0,1,hl) ) return;   // 1324A --> 0 - 1L3L2R4L: nor
  }
  if( keepTrigPatt(PTMS[1],1,2,hl) ) return;   // 1324B --> 1 - 1L3L2R4R: 1L1
  if( keepTrigPatt(PTMS[2],2,1,hl) ) return;   // 1324C --> 2 - 1R3L2R4L: nor
  if( keepTrigPatt(PTMS[3],3,1,hl) ) return;   // 1324D --> 3 - 1R3L2R4R: nor
  if( keepTrigPatt(PTMS[4],4,3,hl) ) return;   // 1364A -->  4 - 1L3L6L4R: 11L
  if( keepTrigPatt(PTMS[5],5,1,hl) ) return;   // 1364B -->  5 - 1R3L6L4R: nor
  if( keepTrigPatt(PTMS[6],6,1,hl) ) return;   // 1364C -->  6 - 1R3R6L4R: nor
  if( keepTrigPatt(PTMS[7],7,5,hl) ) return;   // 1368A -->  7 - 1R3R6L8L: 1H1
  if( keepTrigPatt(PTMS[8],8,6,hl) ) return;   // 1368B -->  8 - 1R3R6R8L: 11H
  if( keepTrigPatt(PTMS[9],9,7,hl) ) return;   // 1368C -->  9 - 1R3L6L8L: H11
  if( keepTrigPatt(PTMS[10],10,1,hl) ) return;   // 5324A --> 10 - 5L3L2R4L: nor
  if( keepTrigPatt(PTMS[11],11,1,hl) ) return;   // 5324B --> 11 - 5L3R2R4L: nor
  if( keepTrigPatt(PTMS[12],12,6,hl) ) return;   // 5324C --> 12 - 5R3R2R4L: 11H
  if( keepTrigPatt(PTMS[13],13,1,hl) ) return;   // 5364A --> 13 - 5L3R6L4L: nor
  if( keepTrigPatt(PTMS[14],14,1,hl) ) return;   // 5364B --> 14 - 5L3R6L4R: nor
  if( keepTrigPatt(PTMS[15],15,5,hl) ) return;   // 5364C --> 15 - 5R3R6L4L: 1H1
  if( keepTrigPatt(PTMS[16],16,1,hl) ) return;   // 5364D --> 16 - 5R3R6L4R: nor
  if( keepTrigPatt(PTMS[17],17,1,hl) ) return;   // 5368A --> 17 - 5L3R6L8L: nor
  if( keepTrigPatt(PTMS[18],18,1,hl) ) return;   // 5368B --> 18 - 5L3R6R8L: nor
  if( keepTrigPatt(PTMS[19],19,4,hl) ) return;   // 5368C --> 19 - 5L3R6R8R: L11
  if( keepTrigPatt(PTMS[20],20,1,hl) ) return;   // 5764A --> 20 - 5R7L6L4R: nor
  if( keepTrigPatt(PTMS[21],21,1,hl) ) return;   // 5764B --> 21 - 5R7L6R4R: nor
  if( keepTrigPatt(PTMS[22],22,7,hl) ) return;   // 5764C --> 22 - 5R7L6L4L: H11
  if( keepTrigPatt(PTMS[23],23,3,hl) ) return;   // 9764A --> 23 - 9L7L6L4R: 11L
  if( keepTrigPatt(PTMS[24],24,2,hl) ) return;   // 9764B --> 24 - 9L7L6R4R: 1L1
  if( keepTrigPatt(PTMS[25],25,4,hl) ) return;   // 9764C --> 25 - 9L7R6R4R: L11
  if( keepTrigPatt(PTMS[26],26,1,hl) ) return;   // 5768A -->  26 - 5L7L6R8L: nor   
  if( RON==1 ){
    if( keepTrigPatt(PTMS[27],27,2,hl) ) return;   // 5768B --> 27 - 5L7L6R8R: 1L1
    if( keepTrigPatt(PTMS[28],28,1,hl) ) return;   // 5768C --> 28 - 5R7L6R8L: nor
    if( keepTrigPatt(PTMS[29],29,1,hl) ) return;   // 5768D --> 29 - 5R7L6R8R: nor
    if( keepTrigPatt(PTMS[30],30,1,hl) ) return;   // 9768A --> 30 - 9L7L6R8L: nor
    if( keepTrigPatt(PTMS[31],31,1,hl) ) return;   // 9768B --> 31 - 9L7R6R8L: nor 
  }
 
  }//end h/l loop



/*
  for(int hl=0;hl<2;hl++){  //hl=0 high   hl=1 low
  if( keepTrigPatt(PTMS[0],0,1,hl) ) return;   // 5768A -->  0 - 5L7L6R8L: nor   
  if( RON==1 ){
    if( keepTrigPatt(PTMS[1],1,2,hl) ) return;   // 5768B -->  1 - 5L7L6R8R: 1L1
    if( keepTrigPatt(PTMS[2],2,1,hl) ) return;   // 5768C -->  2 - 5R7L6R8L: nor
    if( keepTrigPatt(PTMS[3],3,1,hl) ) return;   // 5768D -->  3 - 5R7L6R8R: nor
  }
  if( keepTrigPatt(PTMS[4],4,3,hl) ) return;   // 1364A -->  4 - 1L3L6L4R: 11L
  if( keepTrigPatt(PTMS[5],5,1,hl) ) return;   // 1364B -->  5 - 1R3L6L4R: nor
  if( keepTrigPatt(PTMS[6],6,1,hl) ) return;   // 1364C -->  6 - 1R3R6L4R: nor
  if( keepTrigPatt(PTMS[7],7,5,hl) ) return;   // 1368A -->  7 - 1R3R6L8L: 1H1
  if( keepTrigPatt(PTMS[8],8,6,hl) ) return;   // 1368B -->  8 - 1R3R6R8L: 11H
  if( keepTrigPatt(PTMS[9],9,7,hl) ) return;   // 1368C -->  9 - 1R3L6L8L: H11
  if( keepTrigPatt(PTMS[10],10,1,hl) ) return;   // 5324A --> 10 - 5L3L2R4L: nor
  if( keepTrigPatt(PTMS[11],11,1,hl) ) return;   // 5324B --> 11 - 5L3R2R4L: nor
  if( keepTrigPatt(PTMS[12],12,6,hl) ) return;   // 5324C --> 12 - 5R3R2R4L: 11H
  if( keepTrigPatt(PTMS[13],13,1,hl) ) return;   // 5364A --> 13 - 5L3R6L4L: nor
  if( keepTrigPatt(PTMS[14],14,1,hl) ) return;   // 5364B --> 14 - 5L3R6L4R: nor
  if( keepTrigPatt(PTMS[15],15,5,hl) ) return;   // 5364C --> 15 - 5R3R6L4L: 1H1
  if( keepTrigPatt(PTMS[16],16,1,hl) ) return;   // 5364D --> 16 - 5R3R6L4R: nor
  if( keepTrigPatt(PTMS[17],17,1,hl) ) return;   // 5368A --> 17 - 5L3R6L8L: nor
  if( keepTrigPatt(PTMS[18],18,1,hl) ) return;   // 5368B --> 18 - 5L3R6R8L: nor
  if( keepTrigPatt(PTMS[19],19,4,hl) ) return;   // 5368C --> 19 - 5L3R6R8R: L11
  if( keepTrigPatt(PTMS[20],20,1,hl) ) return;   // 5764A --> 20 - 5R7L6L4R: nor
  if( keepTrigPatt(PTMS[21],21,1,hl) ) return;   // 5764B --> 21 - 5R7L6R4R: nor
  if( keepTrigPatt(PTMS[22],22,7,hl) ) return;   // 5764C --> 22 - 5R7L6L4L: H11
  if( keepTrigPatt(PTMS[23],23,3,hl) ) return;   // 9764A --> 23 - 9L7L6L4R: 11L
  if( keepTrigPatt(PTMS[24],24,2,hl) ) return;   // 9764B --> 24 - 9L7L6R4R: 1L1
  if( keepTrigPatt(PTMS[25],25,4,hl) ) return;   // 9764C --> 25 - 9L7R6R4R: L11
  if( RON==1 ){
    if( keepTrigPatt(PTMS[26],26,1,hl) ) return;   // 1324A --> 26 - 1L3L2R4L: nor
  }
  if( keepTrigPatt(PTMS[27],27,2,hl) ) return;   // 1324B --> 27 - 1L3L2R4R: 1L1
  if( keepTrigPatt(PTMS[28],28,1,hl) ) return;   // 1324C --> 28 - 1R3L2R4L: nor
  if( keepTrigPatt(PTMS[29],29,1,hl) ) return;   // 1324D --> 29 - 1R3L2R4R: nor
  if( RON==1 ){
    if( keepTrigPatt(PTMS[30],30,1,hl) ) return;   // 9768A --> 30 - 9L7L6R8L: nor
  }
  if( keepTrigPatt(PTMS[31],31,1,hl) ) return;   // 9768B --> 31 - 9L7R6R8L: nor 

  }//end h/l loop

*/
  
/* 
  // High level triggers:
  if( keepTrig( 1,accpatB,8) ) return;   // 5768B -->  1 - acc. patt. B
  if( keepTrig( 2,accpatA,8) ) return;   // 5768C -->  2 - acc. patt. A
  if( keepTrig( 3,accpatA,8) ) return;   // 5768D -->  3 - acc. patt. A
  if( keepTrig( 4,accpatB,8) ) return;   // 1364A -->  4 - acc. patt. B
  if( keepTrig( 5,accpatA,8) ) return;   // 1364B -->  5 - acc. patt. A
  if( keepTrig( 6,accpatA,8) ) return;   // 1364C -->  6 - acc. patt. A
  if( keepTrig( 7,accpatB,8) ) return;   // 1368A -->  7 - acc. patt. B
  if( keepTrig( 8,accpatB,8) ) return;   // 1368B -->  8 - acc. patt. B
  if( keepTrig( 9,accpatB,8) ) return;   // 1368C -->  9 - acc. patt. B
  if( keepTrig(10,accpatA,8) ) return;   // 5324A --> 10 - acc. patt. A
  if( keepTrig(11,accpatA,8) ) return;   // 5324B --> 11 - acc. patt. A
  if( keepTrig(12,accpatB,8) ) return;   // 5324C --> 12 - acc. patt. B
  if( keepTrig(13,accpatA,8) ) return;   // 5364A --> 13 - acc. patt. A
  if( keepTrig(14,accpatA,8) ) return;   // 5364B --> 14 - acc. patt. A
  if( keepTrig(15,accpatB,8) ) return;   // 5364C --> 15 - acc. patt. B
  if( keepTrig(16,accpatA,8) ) return;   // 5364D --> 16 - acc. patt. A
  if( keepTrig(17,accpatA,8) ) return;   // 5368A --> 17 - acc. patt. A
  if( keepTrig(18,accpatA,8) ) return;   // 5368B --> 18 - acc. patt. A
  if( keepTrig(19,accpatB,8) ) return;   // 5368C --> 19 - acc. patt. B
  if( keepTrig(20,accpatA,8) ) return;   // 5764A --> 20 - acc. patt. A
  if( keepTrig(21,accpatA,8) ) return;   // 5764B --> 21 - acc. patt. A
  if( keepTrig(22,accpatB,8) ) return;   // 5764C --> 22 - acc. patt. B
  if( keepTrig(23,accpatB,8) ) return;   // 9764A --> 23 - acc. patt. B
  if( keepTrig(24,accpatB,8) ) return;   // 9764B --> 24 - acc. patt. B
  if( keepTrig(25,accpatB,8) ) return;   // 9764C --> 25 - acc. patt. B
  if( keepTrig( 0,accpatA,8) ) return;   // 5768A -->  0 - acc. patt. A
  */ 
 /* 
  // Low level triggers -B
  if( keepTrig( 1,accpatB,2) ) return;   // 5768B -->  1 - acc. patt. B
  if( keepTrig( 2,accpatA,2) ) return;   // 5768C -->  2 - acc. patt. A
  if(tiKes==1) {
    if( keepTrig( 3,accpatA,2) ) return; // 5768D -->  3 - acc. patt. A
  }
  if( keepTrig( 4,accpatB,2) ) return;   // 1364A -->  4 - acc. patt. B
  if( keepTrig( 5,accpatA,2) ) return;   // 1364B -->  5 - acc. patt. A
  if( keepTrig( 6,accpatA,2) ) return;   // 1364C -->  6 - acc. patt. A
  if( keepTrig( 7,accpatB,2) ) return;   // 1368A -->  7 - acc. patt. B
  if( keepTrig( 8,accpatB,2) ) return;   // 1368B -->  8 - acc. patt. B
  if( keepTrig( 9,accpatB,2) ) return;   // 1368C -->  9 - acc. patt. B
  if( keepTrig(10,accpatA,2) ) return;   // 5324A --> 10 - acc. patt. A
  if( keepTrig(11,accpatA,2) ) return;   // 5324B --> 11 - acc. patt. A
  if( keepTrig(12,accpatB,2) ) return;   // 5324C --> 12 - acc. patt. B
  if(tiKes==1) {
    if( keepTrig(13,accpatA,2) ) return; // 5364A --> 13 - acc. patt. A
  }
  if( keepTrig(14,accpatA,2) ) return;   // 5364B --> 14 - acc. patt. A
  if( keepTrig(15,accpatB,2) ) return;   // 5364C --> 15 - acc. patt. B
  if( keepTrig(16,accpatA,2) ) return;   // 5364D --> 16 - acc. patt. A
  if(tiKes==1) {
    if( keepTrig(17,accpatA,2) ) return; // 5368A --> 17 - acc. patt. A
  }
  if( keepTrig(18,accpatA,2) ) return;   // 5368B --> 18 - acc. patt. A
  if( keepTrig(19,accpatB,2) ) return;   // 5368C --> 19 - acc. patt. B
  if( keepTrig(20,accpatA,2) ) return;   // 5764A --> 20 - acc. patt. A
  if(tiKes==1) {
    if( keepTrig(21,accpatA,2) ) return; // 5764B --> 21 - acc. patt. A
  }
  if( keepTrig(22,accpatB,2) ) return;   // 5764C --> 22 - acc. patt. B
  if( keepTrig(23,accpatB,2) ) return;   // 9764A --> 23 - acc. patt. B
  if( keepTrig(24,accpatB,2) ) return;   // 9764B --> 24 - acc. patt. B
  if( keepTrig(25,accpatB,2) ) return;   // 9764C --> 25 - acc. patt. B
  if( keepTrig( 0,accpatA,2) ) return;   // 5768A -->  0 - acc. patt. A

  // Low level triggers -C
  if( keepTrig( 1,accpatB,3) ) return;   // 5768B -->  1 - acc. patt. B
  if( keepTrig( 2,accpatA,3) ) return;   // 5768C -->  2 - acc. patt. A
  if( keepTrig( 3,accpatA,3) ) return;   // 5768D -->  3 - acc. patt. A
  if( keepTrig( 4,accpatB,3) ) return;   // 1364A -->  4 - acc. patt. B
  if( keepTrig( 5,accpatA,3) ) return;   // 1364B -->  5 - acc. patt. A
  if(tiKes==1) {
    if( keepTrig( 6,accpatA,3) ) return; // 1364C -->  6 - acc. patt. A
  }
  if( keepTrig( 7,accpatB,3) ) return;   // 1368A -->  7 - acc. patt. B
  if( keepTrig( 8,accpatB,3) ) return;   // 1368B -->  8 - acc. patt. B
  if( keepTrig( 9,accpatB,3) ) return;   // 1368C -->  9 - acc. patt. B
  if(tiKes==1) {
    if( keepTrig(10,accpatA,3) ) return; // 5324A --> 10 - acc. patt. A
  }
  if( keepTrig(11,accpatA,3) ) return;   // 5324B --> 11 - acc. patt. A
  if( keepTrig(12,accpatB,3) ) return;   // 5324C --> 12 - acc. patt. B
  if( keepTrig(13,accpatA,3) ) return;   // 5364A --> 13 - acc. patt. A
  if( keepTrig(14,accpatA,3) ) return;   // 5364B --> 14 - acc. patt. A
  if( keepTrig(15,accpatB,3) ) return;   // 5364C --> 15 - acc. patt. B
  if(tiKes==1) {
    if( keepTrig(16,accpatA,3) ) return; // 5364D --> 16 - acc. patt. A
  }
  if( keepTrig(17,accpatA,3) ) return;   // 5368A --> 17 - acc. patt. A
  if( keepTrig(18,accpatA,3) ) return;   // 5368B --> 18 - acc. patt. A
  if( keepTrig(19,accpatB,3) ) return;   // 5368C --> 19 - acc. patt. B
  if( keepTrig(20,accpatA,3) ) return;   // 5764A --> 20 - acc. patt. A
  if( keepTrig(21,accpatA,3) ) return;   // 5764B --> 21 - acc. patt. A
  if( keepTrig(22,accpatB,3) ) return;   // 5764C --> 22 - acc. patt. B
  if( keepTrig(23,accpatB,3) ) return;   // 9764A --> 23 - acc. patt. B
  if( keepTrig(24,accpatB,3) ) return;   // 9764B --> 24 - acc. patt. B
  if( keepTrig(25,accpatB,3) ) return;   // 9764C --> 25 - acc. patt. B
  if(tiKes==1) {
    if( keepTrig( 0,accpatA,3) ) return; // 5768A -->  0 - acc. patt. A
  }

  // Low level triggers -A
  if( keepTrig( 1,accpatB,1) ) return;   // 5768B -->  1 - acc. patt. B
  if( keepTrig( 2,accpatA,1) ) return;   // 5768C -->  2 - acc. patt. A
  if( keepTrig( 3,accpatA,1) ) return;   // 5768D -->  3 - acc. patt. A
  if( keepTrig( 4,accpatB,1) ) return;   // 1364A -->  4 - acc. patt. B
  if( keepTrig( 5,accpatA,1) ) return;   // 1364B -->  5 - acc. patt. A
  if( keepTrig( 6,accpatA,1) ) return;   // 1364C -->  6 - acc. patt. A
  if( keepTrig( 7,accpatB,1) ) return;   // 1368A -->  7 - acc. patt. B
  if( keepTrig( 8,accpatB,1) ) return;   // 1368B -->  8 - acc. patt. B
  if(tiKes==1) {
    if( keepTrig( 9,accpatB,1) ) return; // 1368C -->  9 - acc. patt. B
  }
  if( keepTrig(10,accpatA,1) ) return;   // 5324A --> 10 - acc. patt. A
  if( keepTrig(11,accpatA,1) ) return;   // 5324B --> 11 - acc. patt. A
  if( keepTrig(12,accpatB,1) ) return;   // 5324C --> 12 - acc. patt. B
  if( keepTrig(13,accpatA,1) ) return;   // 5364A --> 13 - acc. patt. A
  if( keepTrig(14,accpatA,1) ) return;   // 5364B --> 14 - acc. patt. A
  if( keepTrig(15,accpatB,1) ) return;   // 5364C --> 15 - acc. patt. B
  if( keepTrig(16,accpatA,1) ) return;   // 5364D --> 16 - acc. patt. A
  if( keepTrig(17,accpatA,1) ) return;   // 5368A --> 17 - acc. patt. A
  if( keepTrig(18,accpatA,1) ) return;   // 5368B --> 18 - acc. patt. A
  if(tiKes==1) {
    if( keepTrig(19,accpatB,1) ) return; // 5368C --> 19 - acc. patt. B
  }
  if( keepTrig(20,accpatA,1) ) return;   // 5764A --> 20 - acc. patt. A
  if( keepTrig(21,accpatA,1) ) return;   // 5764B --> 21 - acc. patt. A
  if(tiKes==1) {
    if( keepTrig(22,accpatB,1) ) return; // 5764C --> 22 - acc. patt. B
  }
  if( keepTrig(23,accpatB,1) ) return;   // 9764A --> 23 - acc. patt. B
  if( keepTrig(24,accpatB,1) ) return;   // 9764B --> 24 - acc. patt. B
  if(tiKes==1) {
    if( keepTrig(25,accpatB,1) ) return; // 9764C --> 25 - acc. patt. B
  }
  if( keepTrig( 0,accpatA,1) ) return;   // 5768A -->  0 - acc. patt. A

  // Low level triggers -D
  if( keepTrig( 0,accpatA,4) ) return;   // 5768A -->  0 - acc. patt. A
  if( keepTrig( 1,accpatB,4) ) return;   // 5768B -->  1 - acc. patt. B
  if( keepTrig( 2,accpatA,4) ) return;   // 5768C -->  2 - acc. patt. A
  if( keepTrig( 3,accpatA,4) ) return;   // 5768D -->  3 - acc. patt. A
  if(tiKes==1) {
    if( keepTrig( 4,accpatB,4) ) return; // 1364A -->  4 - acc. patt. B
  }
  if( keepTrig( 5,accpatA,4) ) return;   // 1364B -->  5 - acc. patt. A
  if( keepTrig( 6,accpatA,4) ) return;   // 1364C -->  6 - acc. patt. A
  if( keepTrig( 7,accpatB,4) ) return;   // 1368A -->  7 - acc. patt. B
  if(tiKes==1) {
    if( keepTrig( 8,accpatB,4) ) return; // 1368B -->  8 - acc. patt. B
  }
  if( keepTrig( 9,accpatB,4) ) return;   // 1368C -->  9 - acc. patt. B
  if( keepTrig(10,accpatA,4) ) return;   // 5324A --> 10 - acc. patt. A
  if( keepTrig(11,accpatA,4) ) return;   // 5324B --> 11 - acc. patt. A
  if(tiKes==1) {
    if( keepTrig(12,accpatB,4) ) return; // 5324C --> 12 - acc. patt. B
  }
  if( keepTrig(13,accpatA,4) ) return;   // 5364A --> 13 - acc. patt. A
  if( keepTrig(14,accpatA,4) ) return;   // 5364B --> 14 - acc. patt. A
  if( keepTrig(15,accpatB,4) ) return;   // 5364C --> 15 - acc. patt. B
  if( keepTrig(16,accpatA,4) ) return;   // 5364D --> 16 - acc. patt. A
  if( keepTrig(17,accpatA,4) ) return;   // 5368A --> 17 - acc. patt. A
  if( keepTrig(18,accpatA,4) ) return;   // 5368B --> 18 - acc. patt. A
  if( keepTrig(19,accpatB,4) ) return;   // 5368C --> 19 - acc. patt. B
  if( keepTrig(20,accpatA,4) ) return;   // 5764A --> 20 - acc. patt. A
  if( keepTrig(21,accpatA,4) ) return;   // 5764B --> 21 - acc. patt. A
  if( keepTrig(22,accpatB,4) ) return;   // 5764C --> 22 - acc. patt. B
  if(tiKes==1) {
    if( keepTrig(23,accpatB,4) ) return; // 9764A --> 23 - acc. patt. B
  }
  if( keepTrig(24,accpatB,4) ) return;   // 9764B --> 24 - acc. patt. B
  if( keepTrig(25,accpatB,4) ) return;   // 9764C --> 25 - acc. patt. B
*/
}

DTTrigGeom* DTBtiChip::geom ( ) const

inline

Return trigger geometry.

Definition at line 123 of file DTBtiChip.h.

References _geom.

{ return _geom; }

std::vector<const DTDigi*> DTBtiChip::get_CellDigis ( int  cell )

inline

get digi vector - SV 28/XI/02

Definition at line 74 of file DTBtiChip.h.

References _digis.

{ return _digis[cell];}

DTBtiId DTBtiChip::id ( void  ) const

inline

Return the DTBtiChip Id.

Definition at line 126 of file DTBtiChip.h.

References _id.

Referenced by DTBtiTrig::DTBtiTrig().

{ return _id; }

void DTBtiChip::init ( void  )

private

Definition at line 549 of file DTBtiChip.cc.

References _curStep, _difs, _digis, _hits, _sums, _thisStepUsedHit, begin, config(), gather_cfg::cout, debug, i, AlCaHLTBitMon_ParallelJobs::p, p1, and DTConfigBti::WENflag().

Referenced by run().

                {

  if(config()->debug()>3)
    cout << "DTBtiChip::init() -> initializing bti chip" << endl;

  _curStep=0;
  for(int i=0;i<25;i++) {
    _sums[i] = 1000;
    _difs[i] = 1000;
  } 

  for(int cell=0;cell<9;cell++) {
    int WEN = config()->WENflag(cell+1);
    if( WEN==1 ){
      _thisStepUsedHit[cell]=0;
      vector<const DTDigi*>::const_iterator p;
      for(p=_digis[cell].begin();p<_digis[cell].end();p++){
        DTBtiHit* hit = new DTBtiHit(*p,config());
        //int clockTime = (int)(fabs(((*p)->time()+config()->SetupTime())/12.5));
        //DTBtiHit* hit = new DTBtiHit(clockTime,config());
        _hits[cell].push_back(hit);
      }

      //debugging
      if(config()->debug()>2){
        vector<DTBtiHit*>::const_iterator p1;
        for(p1=_hits[cell].begin();p1<_hits[cell].end();p1++){
            cout << " Filling hit in cell " << cell+1;
                if((*p1)->curTime()!=4000) 
                  cout << " raw time in trigger: " << (*p1)->curTime() << endl;
                cout << " time (clock units): " << (*p1)->clockTime() << endl; 
        }
      }
      // end debugging
    }
  }
}

void DTBtiChip::init_clock

(

)

Definition at line 589 of file DTBtiChip.cc.

References _curStep, _difs, _digis_clock, _hits, _sums, _thisStepUsedHit, begin, config(), gather_cfg::cout, debug, i, init_done, p1, and DTConfigBti::WENflag().

                      {

  if(config()->debug()>3)
    cout << "DTBtiChip::init_clock() -> initializing bti chip" << endl;

  init_done = 1;
  _curStep=0;

  for(int i=0;i<25;i++) {
    _sums[i] = 1000;
    _difs[i] = 1000;
  } 

  for(int cell=0;cell<9;cell++) {
    int WEN = config()->WENflag(cell+1);
    if( WEN==1 ){
      _thisStepUsedHit[cell]=0;
    for(unsigned int i=0; i<_digis_clock[cell].size(); i++){
      const int clockTime = (_digis_clock[cell])[i];
      DTBtiHit* hit = new DTBtiHit(clockTime,config());
      _hits[cell].push_back(hit);
    }
    
    //debugging
    if(config()->debug()>2){
      vector<DTBtiHit*>::const_iterator p1;
      for(p1=_hits[cell].begin();p1<_hits[cell].end();p1++){
        cout << " Filling hit in cell " << cell+1;
        if((*p1)->curTime()!=4000) 
          cout << " time: " << (*p1)->curTime() << endl;
        else
          cout << " time (clock units): " << (*p1)->clockTime() << endl; 
        }
      }
      // end debugging
    }
  }
}

int DTBtiChip::keepTrig ( const int  eq, const int  acp, const int  code  )

private

Definition at line 722 of file DTBtiChipSelTrig.cc.

References abs, HDQMDatabaseProducer::config, gather_cfg::cout, debug, BitArray< N >::element(), i, BitArray< N >::print(), and BitArray< N >::set().

                                                               {

  if(config()->debug()>4){
    cout << "DTBtiChip::keepTrig called with arguments: ";
    cout << eq << ", " << acp << ", " << code << endl; 
  }

  int const KCen = 40; // Arrays will start from 0 --> use 40 instead of 41
  BitArray<80> val0, val1, val2, val3, val4, val5;
  int mm[6];

  int i=0;
  for(i=0;i<6;i++){
    mm[i]=0;
    int mk = (int)(2*_Keq[eq][i]);
    if(abs(mk) > config()->KCut() )
        continue;
    mm[i]=(int)(_Keq[eq][i]*2)+KCen;
    //    if(_Keq[eq][i]<0){
    //      mm[i]=_Keq[eq][i]*2-0.5+KCen;
    //    } else {
    //      mm[i]=_Keq[eq][i]*2+KCen;
    //    }
  }

  if(mm[0]>0 && (code==8 || code==3 || code==4) ){
    val0.set(mm[0]-1);
    val0.set(mm[0]);
    val0.set(mm[0]+1);
    if(acp==2 && (code==8 || code==3) ) {
      val0.set(mm[0]-2);
      val0.set(mm[0]+2);
    }
  }

  if(mm[1]>0 && (code==8 || code==1 || code==4) ){
    val1.set(mm[1]-1);
    val1.set(mm[1]);
    val1.set(mm[1]+1);
    if(acp==2 && code==8 ) {
      val1.set(mm[1]-2);
      val1.set(mm[1]+2);
    }
  }

  if(mm[2]>0 && (code==8 || code==1 || code==2) ){
    val2.set(mm[2]-1);
    val2.set(mm[2]);
    val2.set(mm[2]+1);
    if(acp==2 && (code==8 || code==2) ) {
      val2.set(mm[2]-2);
      val2.set(mm[2]+2);
    }
  }

  if(mm[3]>0 && (code==8 || code==2 || code==4) ){
    val3.set(mm[3]-1);
    val3.set(mm[3]);
    val3.set(mm[3]+1);
  }

  if(mm[4]>0 && (code==8 || code==1 || code==3) ){
    val4.set(mm[4]-1);
    val4.set(mm[4]);
    val4.set(mm[4]+1);
  }

  if(mm[5]>0 && (code==8 || code==2 || code==3) ){
    val5.set(mm[5]);
  }

  // debugging
  if(config()->debug()>4){
    cout << " dump of val arrays: " << endl;
    //    cout << val0.to_string() << endl;
    //    cout << val1.to_string() << endl;
    //    cout << val2.to_string() << endl;
    //    cout << val3.to_string() << endl;
    //    cout << val4.to_string() << endl;
    //    cout << val5.to_string() << endl;
    val0.print();
    cout << endl;
    val1.print();
    cout << endl;
    val2.print();
    cout << endl;
    val3.print();
    cout << endl;
    val4.print();
    cout << endl;
    val5.print();
    cout << endl;
  }
  // end debugging

  int KMax = 0;
  int LKMax = -1;
  for(i=0;i<80;i++){
    int val = val0.element(i)+val1.element(i)+val2.element(i)+
              val3.element(i)+val4.element(i)+val5.element(i);
    //    int val = val0.test(i)+val1.test(i)+val2.test(i)+
    //              val3.test(i)+val4.test(i)+val5.test(i);
    if(val>KMax) {
      KMax=val;
      LKMax=i;
    }
  }
  
  // Note that all bits in val are shifted by one w.r.t. FORTRAN version
  // The output K will be the same because of the different value of Kcen

  if        (KMax==6 && code==8) {
    store(eq,8,LKMax-KCen,int(_Xeq[eq][1]));
    return 1;
  } else if (KMax==3 && code!=8) {
    if(code==1 || code==4) {
      store(eq,code,LKMax-KCen,int(_Xeq[eq][0]));
    } else {
      store(eq,code,LKMax-KCen,int(_Xeq[eq][1]));
    }
    return 1;
  }
  return 0; 
  // return value is ITFL of FORTRAN version
}

int DTBtiChip::keepTrigPatt ( int  flag, const int  eq, const int  pattType, int  hlflag  )

private

Definition at line 327 of file DTBtiChipSelTrig.cc.

References abs, HDQMDatabaseProducer::config, gather_cfg::cout, debug, BitArray< N >::element(), Vispa.Plugins.EdmBrowser.EdmDataAccessor::eq(), i, BitArray< N >::print(), BitArray< N >::reset(), BitArray< N >::set(), and X.

                                                                                      {
  //if pattern is not enabled, return
  if(flag==0)
    return  0;

  int AC1 = config()->AccPattAC1(); //default 0
  int AC2 = config()->AccPattAC2(); //default 3
  int ACH = config()->AccPattACH(); //default 1
  int ACL = config()->AccPattACL(); //default 2
  int tiKes = config()->XON();

  if(config()->debug()>4){
    cout << "DTBtiChip::keepTrigPatt called with arguments: ";
    cout << eq << ", " << pattType << ", " << hlflag << endl; 
    cout<<"AC1,AC2,ACH,ACL="<<AC1<<" "<<AC2<<" "<<ACH<<" "<<ACL<<endl;
  }

  BitArray<80> val0, val1, val2, val3, val4, val5;
  int mm[6];

  //NB _Keq = (Kman - K0)/2   where K0=tmax*2  --->OBSOLETE, now Keq=Kman
  //int K0 = 2 * (config()->ST()/2.); 
  int K0 = int (config()->ST());

  int i=0;
  float Keqfloat[6] = {0,0,0,0,0,0};
  for(i=0;i<6;i++){
    mm[i] = -1;
    int mk = (int)(_Keq[eq][i] - K0);
    //if(abs(mk) > config()->KCut(_id.superlayer()-1))continue;
    if(abs(mk) > 2*K0)  continue;
    Keqfloat[i]=_Keq[eq][i];
    mm[i]=(int)(_Keq[eq][i]);
    //    if(_Keq[eq][i]<0){
    //      mm[i]=_Keq[eq][i]*2-0.5+KCen;
    //    } else {
    //      mm[i]=_Keq[eq][i]*2+KCen;
    //    }
  }

  switch(pattType)
    {
      case 1: //normal pattern
    //if(hlflag==1 && (eq!=2 && eq!=14 && eq!=28) ){  //test for L
/*  if(hlflag==1){  //test for L
          acceptMask(&val0,mm[0],AC1); //eqAB
          acceptMask(&val1,mm[1],AC1); //eqBC
          acceptMask(&val2,mm[2],AC1); //eqCD
        }
        else
*/
        {
          acceptMask(&val0,mm[0],AC2); //eqAB
          acceptMask(&val1,mm[1],AC2); //eqBC
          acceptMask(&val2,mm[2],AC2); //eqCD
        }
    break;

      case 2: //1L1 pattern
    acceptMask(&val0,mm[0],AC1); //eqAB
        acceptMask(&val1,mm[1],ACL); //eqBC
        acceptMask(&val2,mm[2],AC1); //eqCD
    break;

      case 3: //11L pattern
    acceptMask(&val0,mm[0],AC1); //eqAB
        acceptMask(&val1,mm[1],AC1); //eqBC
        acceptMask(&val2,mm[2],ACL); //eqCD
    break;

      case 4: //L11 pattern
    acceptMask(&val0,mm[0],ACL); //eqAB
        acceptMask(&val1,mm[1],AC1); //eqBC
        acceptMask(&val2,mm[2],AC1); //eqCD
    break;

      case 5: //1H1 pattern
    acceptMask(&val0,mm[0],AC1); //eqAB
        acceptMask(&val1,mm[1],ACH); //eqBC
        acceptMask(&val2,mm[2],AC1); //eqCD
    break;

      case 6: //11H pattern
    acceptMask(&val0,mm[0],AC1); //eqAB
        acceptMask(&val1,mm[1],AC1); //eqBC
        acceptMask(&val2,mm[2],ACH); //eqCD
    break;

      case 7: //H11 pattern
    acceptMask(&val0,mm[0],ACH); //eqAB
        acceptMask(&val1,mm[1],AC1); //eqBC
        acceptMask(&val2,mm[2],AC1); //eqCD
    break;
 
      default:
    acceptMask(&val0,mm[0],AC2); //eqAB
        acceptMask(&val1,mm[1],AC2); //eqBC
        acceptMask(&val2,mm[2],AC2); //eqCD
    break;
    }   

  //eq. AC and BD acceptance are always +-1 ->code 00
  int acc = 0;
  acceptMask(&val3,mm[3],acc); //eq.AC
  acceptMask(&val4,mm[4],acc); //eq.BD

  //eq. AD is the reference value!
  if(mm[5]>0){
    val5.set(mm[5]);
  }

  // debugging: print() method prints from last to first bit!
  if(config()->debug()>4){
    cout << " dump of val arrays: " << endl;
    //    cout << val0.to_string() << endl;
    //    cout << val1.to_string() << endl;
    //    cout << val2.to_string() << endl;
    //    cout << val3.to_string() << endl;
    //    cout << val4.to_string() << endl;
    //    cout << val5.to_string() << endl;
    val0.print(); 
    cout << endl;
    val1.print();
    cout << endl;
    val2.print();
    cout << endl;
    val3.print();
    cout << endl;
    val4.print();
    cout << endl;
    val5.print();
    cout << endl;
   }
  // end debugging

  //search for High trigger:
  if(hlflag==0){
    int code = 0;
    int KMax = 0;
    int LKMax = -1;
    for(i=0;i<80;i++){
      int val = val0.element(i)+val1.element(i)+val2.element(i)+
                val3.element(i)+val4.element(i)+val5.element(i);
      if(val>KMax) {
        KMax=val;
        LKMax=i;
      }
    }
  
    //SV: K value is stored in 6 bits, so K=64->0, etc
    if(LKMax>63)
      LKMax-=64;

    if(KMax==6) {
      code=8;
      int X;
      if( eq==0 )  //store Xbc only for patt 0, else store Xad
        X=int(_Xeq[eq][0]);
      else
        X=int(_Xeq[eq][1]);
      store(eq,code,LKMax,X,Keqfloat[0],Keqfloat[1],
           Keqfloat[2],Keqfloat[3],Keqfloat[4],Keqfloat[5]);
      return 1;
    }
    return 0;
  } //end H 

  //search for Low triggers:
  if(hlflag==1){
    int code = 0;
    int RON = config()->RONflag();


    //hit in B is missing
    if(config()->debug()>4)
      cout << "SelTrig: searching low-B" << endl;
    int LKMax = -1;
    for(i=0;i<80;i++){
      int val = val2.element(i)+val3.element(i)+val5.element(i);
      if(val==3){ //Ref. is eqAD
        code=2;
        LKMax=i;
        int storefg = 1; 

        //SV - XON time-ind.Keq suppr. XON=0 do not abilitate patterns 
        if(tiKes==0) {
          if(eq==3 || eq==13 || eq==17 || eq==21 || eq==29){  
            if(config()->debug()>3)
              cout << "SelTrig: doing XON suppression!"<<endl;
            storefg = 0;
          }
        }

        //SV - RON suppression for low triggers
        if( RON==0 ){
          if( eq==19 ){
            if(config()->debug()>3)
              cout << "SelTrig: doing RON low triggers suppression!"<<endl;
            storefg = 0;
          }
        }

        if(storefg){ 
          //SV: K value is stored in 6 bits, so K=64->0, etc
          if(LKMax>63)
            LKMax-=64;
          int X;
          if( eq==0 )
            X=int(_Xeq[eq][0]);
          else
            X=int(_Xeq[eq][1]);
          store(eq,code,LKMax,X,Keqfloat[0],Keqfloat[1],
           Keqfloat[2],Keqfloat[3],Keqfloat[4],Keqfloat[5]);
          return 1;
        }  
      }
    } //end -B Low

    //hit in C is missing
    if(config()->debug()>3)
      cout << "SelTrig: searching low-C" << endl;
     for(i=0;i<80;i++){
      int val = val0.element(i)+val4.element(i)+val5.element(i);
      if(val==3){ //Ref. is eqAD
        code=3;
        LKMax=i;
        int storefg = 1;
 
        //SV - XON time-ind.Keq suppr.
        if(tiKes==0) {
          if(eq==0 || eq==6 || eq==10 || eq==16 || eq==26 || eq==30){  
            if(config()->debug()>3)
              cout << "SelTrig: doing XON suppression!"<<endl;
            storefg = 0;
          }
        }

        if(storefg){
          //SV: K value is stored in 6 bits, so K=64->0, etc
          if(LKMax>63)
            LKMax-=64;

          int X;
          if( eq==0 )
            X=int(_Xeq[eq][0]);
          else
            X=int(_Xeq[eq][1]);
          store(eq,code,LKMax,X,Keqfloat[0],Keqfloat[1],
           Keqfloat[2],Keqfloat[3],Keqfloat[4],Keqfloat[5]);
          return 1;
        }
      }
    } // end -C Low

    //for -A and -D low acceptance is +-1
//    if(pattType==1){
      val0.reset();
      val1.reset();
      val2.reset();
      acceptMask(&val0,mm[0],AC1); //eqAB
      acceptMask(&val1,mm[1],AC1); //eqBC
      acceptMask(&val2,mm[2],AC1); //eqCD
//    }

    //hit in A is missing
    if(config()->debug()>4)
      cout << "SelTrig: searching low-A" << endl;
    for(i=0;i<80;i++){
      int val = val1.element(i)+val2.element(i)+val4.element(i);
      if(val==3 && i==mm[4]){ //Ref. is eqBD
        code=1;
        LKMax=i;
        int storefg = 1;

        //SV - XON time-ind.Keq suppr.
        if(tiKes==0) {
          if(eq==9 || eq==19 || eq==22 || eq==25 ){  
            if(config()->debug()>3)
              cout << "SelTrig: doing low-A XON suppression!"<<endl;
          storefg = 0;
          }
        }

        if( RON==0 ){ //SV - RON suppression
          if( eq==26 ){
            if(config()->debug()>3)
              cout << "SelTrig: doing RON low triggers suppression!"<<endl;
            storefg = 0;
          }
        }

        if(storefg){
          //SV: K value is stored in 6 bits, so K=64->0, etc
          if(LKMax>63)
            LKMax-=64;

          store(eq,code,LKMax,int(_Xeq[eq][0]),Keqfloat[0],Keqfloat[1],
           Keqfloat[2],Keqfloat[3],Keqfloat[4],Keqfloat[5]);
          return 1;
        }
      }
    } //end -A Low
 

    //hit in D is missing
    if(config()->debug()>4)
      cout << "SelTrig: searching low-D" << endl;
    for(i=0;i<80;i++){
      int val = val0.element(i)+val1.element(i)+val3.element(i);
      if(val==3 && i==mm[3]){ //Ref. is eqAC
        code=4;
        LKMax=i;
        int storefg = 1;

        //SV - XON time-ind.Keq suppr.
        if(tiKes==0){
          if(eq==4 || eq==8 || eq==12 || eq==23){  
            if(config()->debug()>3)
              cout << "SelTrig: doing XON suppression!"<<endl;
            storefg = 0;
          }
        }

        //SV - RON suppression for low triggers
        if( RON==0 ){
          if(eq==1 || eq==2 || eq==3 || eq==24 || eq==25){
            if(config()->debug()>3)
              cout << "SelTrig: doing RON low triggers suppression!"<<endl;
            storefg = 0;
          }
        }

        if(storefg){ // && _Xeq[eq][1] >=0){
          //SV: K value is stored in 6 bits, so K=64->0, etc
          if(LKMax>63)
            LKMax-=64;

          store(eq,code,LKMax,int(_Xeq[eq][0]),Keqfloat[0],Keqfloat[1],
           Keqfloat[2],Keqfloat[3],Keqfloat[4],Keqfloat[5]);
          return 1;
        }
      }
    } //end -D Low
 
    return 0; 
  } //end Low
  return 0;
}

LocalPoint DTBtiChip::localPosition ( ) const

inline

Position in chamber frame (x is the one of first traco in slave plane)

Definition at line 95 of file DTBtiChip.h.

References _geom, _id, and DTTrigGeom::localPosition().

{ return _geom->localPosition(_id); }

bool DTBtiChip::matchEq ( float  eqA, float  eqB, int  AC  )

private

int DTBtiChip::nCellHit

(

)

const

Number of cells with hits.

Definition at line 390 of file DTBtiChip.cc.

References _digis, _id, begin, DTBtiId::bti(), config(), gather_cfg::cout, debug, i, n, AlCaHLTBitMon_ParallelJobs::p, and findQualityFiles::size.

Referenced by run().

                          {
  int n=0;
  int i=0;
  for(i=0;i<9;i++) {
    if( _digis[i].size() >0 ) n++;
  }
  if(config()->debug()>2) {
    cout << n << " cells with hits found:" << endl;
  }
  if(config()->debug()>2) {
    for(i=0;i<9;i++) {
      vector<const DTDigi*>::const_iterator p;
      for(p=_digis[i].begin();p<_digis[i].end();p++) {
        cout << "DTBtiChip # " << 
        _id.bti() << 
        " cell " << i+1;
        cout << " --> drift time (tdc units): " << (*p)->countsTDC() << endl;
        (*p)->print();
      }
    }
  }
  return n;
/*
 //SV 2/IV/03 counting hits from _hits
  int n=0;
  int i=0;
  for(i=0;i<9;i++) {
    if( _hits[i].size() >0 ) n++;
  }
  if(config()->debug()>2) {
    cout << n << " cells with hits found:" << endl;
  }
  if(config()->debug()>2) {
    for(i=0;i<9;i++) {
      vector<const DTBtiHit*>::const_iterator p;
      for(p=_hits[i].begin();p<_hits[i].end();p++) {
        cout << "DTBtiChip # " << 
        number() << 
        " cell " << i+1;
        if((*p)->curTime()!=4000)
          cout << " --> drift time: " << (*p)->curTime() << endl;
        else
          cout << " --> clock time: " << (*p)->clockTime() << endl;
      }
    }
  }
  return n;
*/
}

int DTBtiChip::nTrig

(

int

step

)

const

Number of triggers found.

Definition at line 455 of file DTBtiChip.cc.

References _trigs, gather_cfg::cout, DTConfig::NSTEPF, DTConfig::NSTEPL, and relval_parameters_module::step.

Referenced by doLTS(), and DTBtiCard::runBTI().

                               {
  if(step<DTConfig::NSTEPF||step>DTConfig::NSTEPL){
    cout << "DTBtiChip::nTrig: step out of range: " << step ;
    cout << " 0 returned" << endl;
    return 0;
  }
  return _trigs[step-DTConfig::NSTEPF].size(); 
}

int DTBtiChip::number ( ) const

inline

Return DTBtiChip number.

Definition at line 89 of file DTBtiChip.h.

References _id, and DTBtiId::bti().

Referenced by add_digi_clock(), and DTBtiCard::runBTI().

{ return _id.bti(); }

DTBtiChip & DTBtiChip::operator=

(

const DTBtiChip &

bti

)

Assignment operator.

Definition at line 275 of file DTBtiChip.cc.

References _curStep, _difs, _digis, _geom, _hits, _id, _JTR, _Keq, _KTR, _MaxKAcc, _MinKAcc, _nStepUsedHits, _sums, _thisStepUsedHit, _thisStepUsedTimes, _trigs, _Xeq, i, j, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, and p1.

                                         {
  if(this != &bti){
    _geom = bti._geom;
    _id = bti._id;
    _curStep = bti._curStep;
    _nStepUsedHits = bti._nStepUsedHits;
    int i=0;
    for(i=0;i<DTConfig::NSTEPL - DTConfig::NSTEPF;i++) {
      _trigs[i].reserve(2);
      vector<DTBtiTrig*>::const_iterator p;
      for(p=bti._trigs[i].begin();p<bti._trigs[i].end();p++){
    _trigs[i].push_back(*p);
      }
    }
    for(i=0;i<9;i++) {
      _digis[i].reserve(10);
      vector<const DTDigi*>::const_iterator p;
      for(p=bti._digis[i].begin();p<bti._digis[i].end();p++){
    _digis[i].push_back(*p);
      }
      _hits[i].reserve(10);
      vector<DTBtiHit*>::const_iterator p1;
      for(p1=bti._hits[i].begin();p1<bti._hits[i].end();p1++){
    _hits[i].push_back(*p1);
      }
      _thisStepUsedTimes[i] = bti._thisStepUsedTimes[i];
      _thisStepUsedHit[i] = _thisStepUsedHit[i];
    }
    for(i=0;i<25;i++){
      _sums[i] = bti._sums[i];
      _difs[i] = bti._difs[i];
    }
    for(i=0;i<26;i++){
      int j = 0;
      for(j=0;j<6;j++){
    _Keq[i][j] = bti._Keq[i][j];
      }
      for(j=0;j<3;j++){
    _JTR[i][j] = bti._JTR[i][j];
      }
      for(j=0;j<2;j++){
    _Xeq[i][j] = bti._Xeq[i][j];
    _KTR[i][j] = bti._KTR[i][j];
      }
    }
    _MinKAcc = bti._MinKAcc;
    _MaxKAcc = bti._MaxKAcc;
  }
  return *this;
}

int DTBtiChip::reSum ( int  a, int  b  )

inlineprivate

Definition at line 163 of file DTBtiChip.h.

References a, and reSumAr.

{ return reSumAr[a][b+2];}

int DTBtiChip::reSum23 ( int  a, int  b  )

inlineprivate

Definition at line 164 of file DTBtiChip.h.

References a, and reSumAr23.

{ return reSumAr23[a][b+2];}

void DTBtiChip::reSumSet ( )

private

Definition at line 103 of file DTBtiChipEquations.cc.

Referenced by DTBtiChip().

                   {

  reSumAr[2][ 2 +2]=0;
  reSumAr[2][ 1 +2]=0;
  reSumAr[2][ 0 +2]=0;
  reSumAr[2][-1 +2]=-1;
  reSumAr[2][-2 +2]=-1;

  reSumAr[1][ 2 +2]=1;
  reSumAr[1][ 1 +2]=1;
  reSumAr[1][ 0 +2]=0;
  reSumAr[1][-1 +2]=0;
  reSumAr[1][-2 +2]=0;

  reSumAr[0][ 2 +2]=1;
  reSumAr[0][ 1 +2]=0;
  reSumAr[0][ 0 +2]=0;
  reSumAr[0][-1 +2]=0;
  reSumAr[0][-2 +2]=-1;

  reSumAr23[2][ 2 +2]=1;
  reSumAr23[2][ 1 +2]=1;
  reSumAr23[2][ 0 +2]=1;
  reSumAr23[2][-1 +2]=0;
  reSumAr23[2][-2 +2]=0;

  reSumAr23[1][ 2 +2]=1;
  reSumAr23[1][ 1 +2]=1;
  reSumAr23[1][ 0 +2]=0;
  reSumAr23[1][-1 +2]=0;
  reSumAr23[1][-2 +2]=0;

  reSumAr23[0][ 2 +2]=1;
  reSumAr23[0][ 1 +2]=0;
  reSumAr23[0][ 0 +2]=0;
  reSumAr23[0][-1 +2]=0;
  reSumAr23[0][-2 +2]=-1;
 }

void DTBtiChip::run

(

void

)

Run DTBtiChip algorithm.

Definition at line 630 of file DTBtiChip.cc.

References _geom, _id, _nStepUsedHits, DTBtiId::bti(), computeEqs(), computeSums(), config(), gather_cfg::cout, currentIntStep(), currentStep(), debug, DTConfigBti::debug(), doLTS(), findTrig(), init(), init_done, DTTrigGeom::nCell(), nCellHit(), DTConfig::NSTEPF, DTConfig::NSTEPL, DTBtiId::superlayer(), superlayer(), and tick().

Referenced by DTBtiCard::runBTI().

               {

  // Debugging...
  if(config()->debug()>2){
    cout << "DTBtiChip::run: Processing BTI " << _id.bti() << endl;
    cout << " in SL " << _id.superlayer() << endl;
  }
  // End debugging

  if(_id.bti()<1 || _id.bti() >_geom->nCell(superlayer())) {
    if(config()->debug()>1)
      cout << "DTBtiChip::run : wrong BTI number: " << _id.bti() << endl;
    return;
  }

  // run algorithm
  if(!init_done)
    init();
  if( nCellHit()<3 ) return;   // check that at least 3 cell have hits

  for(int ints=0; ints<2*DTConfig::NSTEPL; ints++) { // 80 MHz 
    tick(); // Do a 12.5 ns step

    // In electronics equations are computed every 12.5 ns
    // but since triggers are searched for only every 25 ns, skip
    // also equation's computing at odd values of internal step
    if((currentIntStep()/2)*2!=currentIntStep())continue;
    //if((currentIntStep()/2)*2==currentIntStep())continue; 


    if(config()->debug()>2){
      cout << "DTBtiChip::run : internal step " << currentIntStep();
      cout << " number of JTRIG hits is " << _nStepUsedHits << endl;
    }
    if(currentStep()>=DTConfig::NSTEPF && _nStepUsedHits>2) { 
      // at least 3 good hits in this step -> run algorithm
      computeSums();
      computeEqs();
      findTrig();
    }
  }
  if( config()->LTS()>0 ) doLTS(); // low trigger suppression
}

int DTBtiChip::sector ( void  ) const

inline

Return sector number.

Definition at line 135 of file DTBtiChip.h.

References _geom, and DTTrigGeom::sector().

Referenced by geometryXMLparser.DTAlignable::index().

{ return _geom->sector(); }

void DTBtiChip::setSnap ( )

private

Definition at line 852 of file DTBtiChip.cc.

References config(), gather_cfg::cout, debug, RE23, DTConfigBti::RE43(), RE43, ST, ST2, ST23, ST3, ST4, DTConfigBti::ST43(), ST43, ST5, and ST7.

Referenced by DTBtiChip().

                  {

 //set the internally calculated drift velocity parameters
  ST43 = config()->ST43();
  RE43 = config()->RE43();
  ST23 = int(double(ST43)/2.);
  RE23 = (RE43==1) ? 2 : int(double(RE43)/2.);


  ST =  int(  double(ST43) * 3./4. + double(RE43) * 1./4.     );
  ST2 = int( (double(ST43) * 3./4. + double(RE43) * 1./4.)*2. );
  ST3 = int( (double(ST43) * 3./4. + double(RE43) * 1./4.)*3. );
  ST4 = int( (double(ST43) * 3./4. + double(RE43) * 1./4.)*4. );
  ST5 = int( (double(ST43) * 3./4. + double(RE43) * 1./4.)*5. );
  ST7 = int( (double(ST43) * 3./4. + double(RE43) * 1./4.)*7. );

  if(config()->debug()>3){
    cout << "Snap register dump: " << endl;
    cout << "ST43 = " << ST43 << endl;
    cout << "RE43 = " << RE43 << endl;
    cout << "ST23 = " << ST23 << endl;
    cout << "RE23 = " << RE23 << endl;
    cout << "ST = " << ST << endl;
    cout << "ST2 = " << ST2 << endl;
    cout << "ST3 = " << ST3 << endl;
    cout << "ST4 = " << ST4 << endl;
    cout << "ST5 = " << ST5 << endl;
    cout << "ST7 = " << ST7 << endl;
  }
}

int DTBtiChip::station ( ) const

inline

Return station number.

Definition at line 132 of file DTBtiChip.h.

References _geom, and DTTrigGeom::station().

Referenced by geometryXMLparser.DTAlignable::index(), and geometryXMLparser.CSCAlignable::index().

{ return _geom->station(); }

int DTBtiChip::store ( const int  eq, const int  code, const int  K, const int  X, float  KeqAB = 0., float  KeqBC = 0., float  KeqCD = 0., float  KeqAC = 0., float  KeqBD = 0., float  KeqAD = 0.  )

private

Definition at line 785 of file DTBtiChip.cc.

References _MaxKAcc, _MinKAcc, _thisStepUsedHit, DTBtiTrig::addDigi(), addTrig(), trackerHits::c, config(), gather_cfg::cout, currentStep(), debug, DTBtiHit::hitDigi(), DTBtiTrig::print(), and generateDummyData::trg().

                                                            {

  // remove negative position triggers
  if(X<0)return 0;

  
  // accept in range triggers (acceptances defined in constructor)
  if(K>=_MinKAcc && K<=_MaxKAcc) 
  {
    int trig_step = currentStep();

/*
    //SV test 27/I/2003 1-clock delay for critical patterns in default ACx configuration 
    int AC1 = config()->AccPattAC1(); //default 0
    int AC2 = config()->AccPattAC2(); //default 3
    int ACH = config()->AccPattACH(); //default 1
    int ACL = config()->AccPattACL(); //default 2     

    if(AC1==0 && AC2==3 && ACH==1 && ACL==2){
      if(eq==1 || eq==4 || eq==7 || eq==8 || eq==9 || eq==12 || eq==15
    || eq==19 || eq==22 || eq==23 || eq==24 || eq==25 || eq==27 )
       trig_step = currentStep()+1;
    }
*/     
    //store strobe: SV no strobe defined for this setup SV 15/I/2007
    int strobe=-1;

    // create a new trigger
    float Keq[6] = {KeqAB,KeqBC,KeqCD,KeqAC,KeqBD,KeqAD};
    //DTBtiTrig* trg = new DTBtiTrig(this,code,K,X,currentStep(),eq);
    DTBtiTrig* trg = new DTBtiTrig(this,code,K,X,trig_step,eq,strobe,Keq);

    // store also the digis
    for(int c=0; c<9; c++) {
      if(_thisStepUsedHit[c]) {
        const DTDigi* digi = _thisStepUsedHit[c]->hitDigi();
        if(digi)
          trg->addDigi(digi);
      }
    }

    //addTrig(currentStep(),trg);
    addTrig(trig_step,trg);

    // Debugging...
    if(config()->debug()>1)
      trg->print();
    // end debugging
    
    return 1;
  }
  else{
    // remove out of range triggers (acceptances defined in constructor)
    if(config()->debug()>2){
      cout << "DTBtiChip::store, REJECTED TRIGGER at step "<< currentStep();
      cout << " allowed K range in theta view is: [";
      cout << _MinKAcc << ","<< _MaxKAcc << "]";
      cout << "K value is " << K << endl; 
    }
    return 0;
  }//end else
}//end store

void DTBtiChip::sum ( const int  s, const int  a, const int  b  )

private

Definition at line 82 of file DTBtiChipEquations.cc.

                                                    {
  //
  // fill the sums and difs arrays
  // a and b are the cell numbers (in the range [1,9])
  // s is the sum number (in the range [1,25])
  //

  if( _thisStepUsedHit[a-1]!=0 && _thisStepUsedHit[b-1]!=0 ){
    _sums[s-1] = (float)(_thisStepUsedHit[a-1]->jtrig() +
                         _thisStepUsedHit[b-1]->jtrig()  );
    _difs[s-1] = (float)(_thisStepUsedHit[a-1]->jtrig() -
                         _thisStepUsedHit[b-1]->jtrig()  );
  } else {
    _sums[s-1] = 1000;
    _difs[s-1] = 1000;
  }

}

int DTBtiChip::superlayer ( ) const

inline

Return superlayer.

Definition at line 92 of file DTBtiChip.h.

References _id, and DTBtiId::superlayer().

Referenced by add_digi(), geometryXMLparser.DTAlignable::index(), and run().

{ return _id.superlayer(); }

void DTBtiChip::tick ( )

private

Definition at line 675 of file DTBtiChip.cc.

References _curStep, _hits, _nStepUsedHits, _thisStepUsedHit, begin, config(), gather_cfg::cout, currentIntStep(), DTBtiHit::curTime(), debug, DTConfigBti::debug(), DTBtiHit::jtrig(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by run().

                {
  //
  // fills the DTBtiChip registers ( _thisStepUsedHit[cell] )
  // for a given clock (Syncronizer and Shaper functionalities)
  //

  _curStep++; // increase internal step (12.5 ns --> 80 MHz)

  // debugging
  if(config()->debug()>2){
    cout << "DTBtiChip::tick: internal step is now " << currentIntStep()<< endl; 
  }
  // end debugging

  // Loop on cells
  _nStepUsedHits=0;
  for(int cell=0;cell<9;cell++) {

    // decrease drift time by 12.5 ns for each hit
    vector<DTBtiHit*>::const_iterator p;
    for(p=_hits[cell].begin();p<_hits[cell].end();p++){
      (*p)->stepDownTime();
    }

    // loop on hits
    _thisStepUsedHit[cell]=0;
    for(p=_hits[cell].begin();p<_hits[cell].end();p++){
      if       ( (*p)->isDrifting() ) { // hit is drifting
        break;                          //   --> don't consider others
      } else if( (*p)->isInsideReg() ) {  // hit is already in registers
    _thisStepUsedHit[cell]=(*p);
    _nStepUsedHits++;
    // debugging
    if(config()->debug()>2){
          if((*p)->curTime() != 4000)
            cout << "DTBtiChip::tick: hit in register: time=" << (*p)->curTime();
          else
            cout << "DTBtiChip::tick: hit in register! " << endl;
          cout <<                           " jtrig=" << (*p)->jtrig() << endl;

    }
    // end debugging
        break;                          //   --> don't consider other triggers
      }
      // hit is not drifting and not in registers: it is gone out of register, but
      // jtrig value is still=ST ; save in array and consider next one if exists
    } // end loop on cell hits

    // debugging...
    if(config()->debug()>2){
      if(_thisStepUsedHit[cell]!=0){
    cout << "int. step=" << currentIntStep() << " cell=" << cell+1;
    cout << " jtrig=" << _thisStepUsedHit[cell]->jtrig();
        if( _thisStepUsedHit[cell]->curTime() != 4000 )  
      cout << " (time=" << _thisStepUsedHit[cell]->curTime() << ")" << endl;
        else 
          cout << endl;
      }
    } 
    // end debugging

  } // end loop on cells

}

DTBtiTrig * DTBtiChip::trigger	(	int	step,
		unsigned	n
	)		const

Return the requested trigger.

Definition at line 475 of file DTBtiChip.cc.

References _trigs, gather_cfg::cout, n, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, findQualityFiles::size, and relval_parameters_module::step.

Referenced by doLTS().

                                             {
  if(step<DTConfig::NSTEPF||step>DTConfig::NSTEPL){
    cout << "DTBtiChip::trigger: step out of range: " << step ;
    cout << " empty pointer returned" << endl;
    return 0;
  } 
  if(n<1 || n>_trigs[step-DTConfig::NSTEPF].size()) {
    cout << "DTBtiChip::trigger: requested trigger does not exist: " << n;
    cout << " empty pointer returned!" << endl;
    return 0;
  }
  vector<DTBtiTrig*>::const_iterator p = _trigs[step-DTConfig::NSTEPF].begin();
  return (*(p+n-1));
}

DTBtiTrigData DTBtiChip::triggerData	(	int	step,
		unsigned	n
	)		const

Return the data part of the requested trigger.

Definition at line 491 of file DTBtiChip.cc.

References _trigs, gather_cfg::cout, n, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, findQualityFiles::size, and relval_parameters_module::step.

Referenced by DTBtiCard::runBTI().

                                                 {
  if(step<DTConfig::NSTEPF||step>DTConfig::NSTEPL){
    cout << "DTBtiChip::triggerData: step out of range: " << step ;
    cout << " dummy trigger returned" << endl;
    return DTBtiTrigData();
  } 
  if(n<1 || n>_trigs[step-DTConfig::NSTEPF].size()) {
    cout << "DTBtiChip::triggerData: requested trig. doesn't exist: " << n;
    cout << " dummy trigger returned!" << endl;
    return DTBtiTrigData();
  }
  vector<DTBtiTrig*>::const_iterator p = _trigs[step-DTConfig::NSTEPF].begin();
  return (*(p+n-1))->data();
}

vector< DTBtiTrig * > DTBtiChip::trigList

(

int

step

)

const

Definition at line 465 of file DTBtiChip.cc.

References _trigs, gather_cfg::cout, DTConfig::NSTEPF, DTConfig::NSTEPL, and relval_parameters_module::step.

                                  {
  if(step<DTConfig::NSTEPF||step>DTConfig::NSTEPL){
    cout << "DTBtiChip::trigList: step out of range: " << step ;
    cout << " empty pointer returned" << endl;
    //return 0;
  } 
  return _trigs[step-DTConfig::NSTEPF]; 
}

int DTBtiChip::wheel ( ) const

inline

Return wheel number.

Definition at line 129 of file DTBtiChip.h.

References _geom, and DTTrigGeom::wheel().

Referenced by geometryXMLparser.DTAlignable::index().

{ return _geom->wheel(); }

Member Data Documentation

int DTBtiChip::_busyStart_clock[9]

private

Definition at line 208 of file DTBtiChip.h.

Referenced by add_digi(), add_digi_clock(), and DTBtiChip().

DTBtiCard* DTBtiChip::_card

private

Definition at line 169 of file DTBtiChip.h.

Referenced by DTBtiChip().

DTConfigBti* DTBtiChip::_config

private

Definition at line 172 of file DTBtiChip.h.

Referenced by config().

int DTBtiChip::_curStep

private

Definition at line 185 of file DTBtiChip.h.

Referenced by currentIntStep(), currentStep(), init(), init_clock(), operator=(), and tick().

float DTBtiChip::_difs[25]

private

Definition at line 190 of file DTBtiChip.h.

Referenced by DTBtiChip(), init(), init_clock(), and operator=().

std::vector<const DTDigi*> DTBtiChip::_digis[9]

private

Definition at line 177 of file DTBtiChip.h.

Referenced by add_digi(), clear(), DTBtiChip(), get_CellDigis(), init(), nCellHit(), and operator=().

std::vector<int > DTBtiChip::_digis_clock[9]

private

Definition at line 179 of file DTBtiChip.h.

Referenced by add_digi_clock(), clear(), and init_clock().

DTTrigGeom* DTBtiChip::_geom

private

Definition at line 171 of file DTBtiChip.h.

Referenced by add_digi(), CMSPosition(), DTBtiChip(), geom(), localPosition(), operator=(), run(), sector(), station(), and wheel().

std::vector<DTBtiHit*> DTBtiChip::_hits[9]

private

Definition at line 186 of file DTBtiChip.h.

Referenced by clear(), DTBtiChip(), init(), init_clock(), operator=(), and tick().

DTBtiId DTBtiChip::_id

private

Definition at line 174 of file DTBtiChip.h.

Referenced by add_digi(), CMSPosition(), DTBtiChip(), id(), localPosition(), nCellHit(), number(), operator=(), run(), and superlayer().

float DTBtiChip::_JTR[32][3]

private

Definition at line 206 of file DTBtiChip.h.

Referenced by DTBtiChip(), and operator=().

float DTBtiChip::_Keq[32][6]

private

Definition at line 191 of file DTBtiChip.h.

Referenced by DTBtiChip(), and operator=().

float DTBtiChip::_KTR[32][2]

private

Definition at line 205 of file DTBtiChip.h.

Referenced by DTBtiChip(), and operator=().

int DTBtiChip::_MaxKAcc

private

Definition at line 194 of file DTBtiChip.h.

Referenced by DTBtiChip(), operator=(), and store().

int DTBtiChip::_MaxKcenterTraco

private

Definition at line 198 of file DTBtiChip.h.

int DTBtiChip::_MaxKleftTraco

private

Definition at line 196 of file DTBtiChip.h.

int DTBtiChip::_MaxKrightTraco

private

Definition at line 200 of file DTBtiChip.h.

int DTBtiChip::_MinKAcc

private

Definition at line 193 of file DTBtiChip.h.

Referenced by DTBtiChip(), operator=(), and store().

int DTBtiChip::_MinKcenterTraco

private

Definition at line 197 of file DTBtiChip.h.

int DTBtiChip::_MinKleftTraco

private

Definition at line 195 of file DTBtiChip.h.

int DTBtiChip::_MinKrightTraco

private

Definition at line 199 of file DTBtiChip.h.

int DTBtiChip::_nStepUsedHits

private

Definition at line 189 of file DTBtiChip.h.

Referenced by operator=(), run(), and tick().

float DTBtiChip::_sums[25]

private

Definition at line 190 of file DTBtiChip.h.

Referenced by DTBtiChip(), init(), init_clock(), and operator=().

DTBtiHit* DTBtiChip::_thisStepUsedHit[9]

private

Definition at line 188 of file DTBtiChip.h.

Referenced by DTBtiChip(), init(), init_clock(), operator=(), store(), and tick().

int DTBtiChip::_thisStepUsedTimes[9]

private

Definition at line 187 of file DTBtiChip.h.

Referenced by DTBtiChip(), and operator=().

std::vector<DTBtiTrig*> DTBtiChip::_trigs[DTConfig::NSTEPL-DTConfig::NSTEPF+1]

private

Definition at line 182 of file DTBtiChip.h.

Referenced by addTrig(), clear(), DTBtiChip(), eraseTrigger(), nTrig(), operator=(), trigger(), triggerData(), and trigList().

float DTBtiChip::_Xeq[32][2]

private

Definition at line 192 of file DTBtiChip.h.

Referenced by DTBtiChip(), and operator=().

float DTBtiChip::_XeqAB_patt0

private

Definition at line 203 of file DTBtiChip.h.

float DTBtiChip::_XeqAC_patt0

private

Definition at line 202 of file DTBtiChip.h.

float DTBtiChip::_XeqBD_patt0

private

Definition at line 202 of file DTBtiChip.h.

float DTBtiChip::_XeqCD_patt0

private

Definition at line 203 of file DTBtiChip.h.

int DTBtiChip::init_done

private

Definition at line 207 of file DTBtiChip.h.

Referenced by DTBtiChip(), init_clock(), and run().

int DTBtiChip::RE23

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::RE43

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::reSumAr[3][5]

private

Definition at line 213 of file DTBtiChip.h.

Referenced by reSum().

int DTBtiChip::reSumAr23[3][5]

private

Definition at line 214 of file DTBtiChip.h.

Referenced by reSum23().

int DTBtiChip::ST

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::ST2

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::ST23

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::ST3

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::ST4

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::ST43

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::ST5

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().

int DTBtiChip::ST7

private

Definition at line 211 of file DTBtiChip.h.

Referenced by setSnap().