#include <DTTracoChip.h>
Public Member Functions | |
void | add_btiT (int step, int pos, const DTBtiTrigData *btitrig) |
Add a BTI trigger to the TRACO. | |
void | addTrig (int step, DTTracoTrig *) |
Add a TRACO trigger. | |
int | BTIC () const |
BTIC parameter. | |
void | calculateAngles (DTTracoTrig *) |
Calculate trigger angles. | |
void | clear () |
Clear. | |
GlobalPoint | CMSPosition () const |
Position in CMS frame. | |
DTConfigTraco * | config () const |
Configuration set. | |
int | DD () const |
DD parameter. | |
DTTracoChip (DTTracoCard *card, int n, DTConfigTraco *config) | |
Constructor. | |
DTTracoChip (const DTTracoChip &traco) | |
Copy constructor. | |
int | edgeBTI (int step, int io, int lr) const |
flags for HTRIG in edge BTI | |
DTTrigGeom * | geom () const |
Return trigger geometry. | |
int | IBTIOFF () const |
IBTIOFF parameter. | |
DTTracoId | id () const |
Return TRACO id. | |
int | KRad () const |
K par of the radial angle of corr center referred to plane sl. | |
LocalPoint | localPosition () const |
Position in chamber frame. | |
int | nTrig (int step) const |
Return the number of trigger candidates. | |
int | number () const |
Return TRACO number. | |
DTTracoChip & | operator= (const DTTracoChip &traco) |
Assignment operator. | |
int | ovlFlag (int step) |
return overlap flag | |
float | psiRad (int sl=0) const |
Radial angle of correlator center in mrad referred to plane sl. | |
void | raiseOverlap (int step) |
raise overlap flag | |
void | run () |
Run TRACO algorithm. | |
int | sector () const |
Return sector number. | |
void | setFlag (int step, int ext=0) |
Set flags for multiple trigger detection between cons. TRACO's. | |
void | setPV (int step, int ntrk, int code, int K) |
Set the preview values for a TRACO trigger. | |
int | station () const |
Return station number. | |
DTTracoTrig * | trigger (int step, unsigned n) const |
Return the requested trigger. | |
DTTracoTrigData | triggerData (int step, unsigned n) const |
Return the data part of the requested trigger. | |
int | useSecondTrack (int step) const |
a flag for a usable second track | |
int | wheel () const |
Return wheel number. | |
~DTTracoChip () | |
Destructor. | |
Private Member Functions | |
int | AdjBtiLTSuppressed (DTTracoCand *candidate) |
Do suppression of LTRIG on adjacent TRACO. | |
DTTracoCand * | bestCand (int itk, std::vector< DTTracoCand > &tclist) |
Get the best inner/outer candidate. | |
void | DoAdjBtiLts (DTTracoCand *candidate, std::vector< DTTracoCand > &tclist) |
Do suppression of LTRIG on BTI close to selected HTRIG. | |
int | insideAngWindow (DTTracoTrig *) const |
Check if a trigger is inside the angular acceptance window. | |
DTTracoTrig * | setPV (int itk, DTTracoCand *inner, DTTracoCand *outer) |
Set the preview for a trigger. | |
void | setTracoAcceptances () |
Compute traco chip acceptances. | |
int | storeCorr (DTTracoTrig *tctrig, DTTracoCand *inner, DTTracoCand *outer, int tkn) |
Check correlation and store correlated trigger. | |
int | storeUncorr (DTTracoTrig *tctrig, DTTracoCand *inner, DTTracoCand *outer, int tkn) |
Store uncorrelated trigger. | |
Private Attributes | |
int | _btic |
BitArray< DTConfigTraco::NSTEPL+2 > | _bxlts |
DTTracoCard * | _card |
DTConfigTraco * | _config |
int | _dd |
BitArray< 32 > | _flag [DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1] |
DTTrigGeom * | _geom |
int | _ibtioff |
DTTracoId | _id |
std::vector< DTTracoCand > | _innerCand [DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1] |
int | _krad |
DTTracoLUTs * | _luts |
Lut * | _lutsCCB |
std::vector< DTTracoCand > | _outerCand [DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1] |
int | _PSIMAX [4 *DTConfig::NBTITC] |
int | _PSIMIN [4 *DTConfig::NBTITC] |
std::vector< DTTracoTrig * > | _tracotrig [DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1] |
Implementation of TRACO trigger algorithm. Internally uses DTTracoCand to store BTI triggers
Definition at line 55 of file DTTracoChip.h.
DTTracoChip::DTTracoChip | ( | DTTracoCard * | card, |
int | n, | ||
DTConfigTraco * | config | ||
) |
Constructor.
Constructor for passing configuration
Definition at line 47 of file DTTracoChip.cc.
References _btic, _bxlts, _card, _flag, _geom, _ibtioff, _id, _innerCand, _krad, _luts, _lutsCCB, _outerCand, _tracotrig, DTConfigTraco::BTIC(), DTTrigGeom::cellPitch(), CMSPosition(), config(), DTTracoCard::config_luts(), gather_cfg::cout, debug, DTSuperLayerId, DTGeomSupplier::geom(), i, DTConfigTraco::KRAD(), KRad(), DTTrigGeom::localPosition(), DTTracoCard::lutFromDBFlag(), n, DTConfig::NBTITC, DTConfig::NSTEPF, DTConfig::NSTEPL, DTTrigGeom::phiSLOffset(), psiRad(), DTChamberId::sector(), setTracoAcceptances(), evf::utils::sid, DTTrigGeom::statId(), DTChamberId::station(), DTChamberId::wheel(), x, and BitArray< N >::zero().
: _card(card), _config(conf) { _geom = _card->geom(); // n=traco number 1,2,... if(config()->debug()==4){ std::cout << "DTTracoChip constructor called for TRACO number " << n << std::endl; } // set acceptances from CMSSW geometry setTracoAcceptances(); // reserve the appropriate amount of space for vectors int i=0; for(i=0;i<DTConfigTraco::NSTEPL - DTConfigTraco::NSTEPF;i++) { _innerCand[i].reserve(DTConfigTraco::NBTITC); _outerCand[i].reserve(3*DTConfigTraco::NBTITC); _tracotrig[i].reserve(2); } // the identifier DTChamberId sid = _geom->statId(); _id = DTTracoId(sid,n); // Flags for LTS _bxlts.zero(); for(int is=0;is<DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1;is++){ _flag[is].zero(); } // debugging if(config()->debug()==4){ std::cout << "CMS position:" << CMSPosition() << std::endl; std::cout << " psiRad=" << psiRad() << " KRad=" << KRad() << std::endl; } //init traco parameters from traco config file _krad = config()->KRAD(); _btic=config()->BTIC(); //offset from geometry (x1-x3 FE view): converted from cm to ST units (0.9999 for rounding) _ibtioff=static_cast<int>(config()->BTIC()/(_geom->cellPitch())*(_geom->phiSLOffset()/0.9999)); // 091030 SV lut parameters from DB // SV 08/12/12 : added flag for computing luts from DB parameters if( _card->lutFromDBFlag()==1 ) { //int board = int( (n-1)/4 ); //int traco = int(fmod( double(n-1),4.)); // 110208 SV for TRACO hardware bug included SL_shift // SL shift float xBTI1_3 = _geom->localPosition( DTBtiId(DTSuperLayerId(sid.wheel(),sid.station(),sid.sector(),3),1) ).x(); float xBTI1_1 = _geom->localPosition( DTBtiId(DTSuperLayerId(sid.wheel(),sid.station(),sid.sector(),1),1) ).x(); float SL_shift = xBTI1_3 - xBTI1_1; _lutsCCB = new Lut(_card->config_luts(),n,SL_shift); _luts = 0; } else //this is always the case with new DTConfig SV 15/I/2007 //if( config()->trigSetupGeom()==0 ){ { _luts = 0; _lutsCCB = 0; } /* //SV 21/V/03 for testbeam purpose: parameters from hardware setup if(config()->trigSetupGeom()==1){ //init traco parameters _dd=config()->DD(); _krad=config()->KRAD(); _btic=config()->BTIC(); _ibtioff=config()->IBTIOFF(); //const char* testfile = "traco"; //FIXTB std::string testfile = "traco_"; int it = number() - 4 - 1; if(it<10) testfile += it+'0'; if(it>9){ testfile += int(float(it)/10.0) + '0'; //add decimal char testfile += int(fmod(float(it),float(10))) + '0'; //add unit char } //const char* name = testfile; if(config()->debug()==4) std::cout << "Loading " << testfile << " luts for traco " << number() << std::endl; _luts = new DTTracoLUTs(testfile); _luts->reset(); _luts->load(); if(config()->debug()==4) _luts->print(); }//end hardware setup //TB 2004 setup : luts from minicrate ccb equations if(config()->trigSetupGeom()==2){ int board = int( (n-1)/4 ); int traco = fmod( double(n-1),4. ); _lutsCCB = new Lut(sid.station(),board,traco); // 091030 SV this constructur is obsolete now, use setForTestBeam instead }//end TB2004 */ }
DTTracoChip::DTTracoChip | ( | const DTTracoChip & | traco | ) |
Copy constructor.
Definition at line 153 of file DTTracoChip.cc.
References _bxlts, _flag, _innerCand, _outerCand, _tracotrig, i, DTConfig::NBTITC, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, and p1.
: _geom(traco._geom), _id(traco._id), _card(traco._card), _luts(traco._luts) { int i=0; for(i=0;i<DTConfigTraco::NSTEPL - DTConfigTraco::NSTEPF;i++) { _innerCand[i].reserve(DTConfigTraco::NBTITC); std::vector<DTTracoCand>::const_iterator p; for(p=traco._innerCand[i].begin();p<traco._innerCand[i].end();p++) { _innerCand[i].push_back(*p); } _outerCand[i].reserve(3*DTConfigTraco::NBTITC); for(p=traco._outerCand[i].begin();p<traco._outerCand[i].end();p++) { _outerCand[i].push_back(*p); } _tracotrig[i].reserve(2); std::vector<DTTracoTrig*>::const_iterator p1; for(p1=traco._tracotrig[i].begin();p1<traco._tracotrig[i].end();p1++) { _tracotrig[i].push_back(*p1); } } _bxlts = traco._bxlts; for(int is=0;is<DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1;is++){ _flag[is] = traco._flag[is]; } }
DTTracoChip::~DTTracoChip | ( | ) |
Destructor.
Definition at line 183 of file DTTracoChip.cc.
References _card, _lutsCCB, clear(), and DTTracoCard::lutFromDBFlag().
{ clear(); /* if(config()->trigSetupGeom()==1){ _luts->reset(); delete _luts; } if(config()->trigSetupGeom()==2) delete _lutsCCB; */ if( _card->lutFromDBFlag()==1 ) delete _lutsCCB; }
void DTTracoChip::add_btiT | ( | int | step, |
int | pos, | ||
const DTBtiTrigData * | btitrig | ||
) |
Add a BTI trigger to the TRACO.
Definition at line 990 of file DTTracoChip.cc.
References _bxlts, _card, _id, _innerCand, _outerCand, _PSIMAX, _PSIMIN, DTBtiTrigData::code(), config(), gather_cfg::cout, debug, DTConfigTraco::debug(), BitArray< N >::element(), DTBtiTrigData::K(), DTConfig::NBTITC, DTConfig::NSTEPF, DTConfig::NSTEPL, number(), pos, DTBtiTrigData::print(), BitArray< N >::set(), launcher::step, DTTracoId::traco(), and DTTracoCard::useAcceptParamFlag().
Referenced by DTTracoCard::loadTRACO().
{ if(pos<1 || pos>4*DTConfigTraco::NBTITC) { std::cout << "DTTracoChip::add_btiT: wrong position: " << pos; std::cout << "trigger not added!" << std::endl; return; } if(step<DTConfigTraco::NSTEPF||step>DTConfigTraco::NSTEPL){ std::cout << "DTTracoChip::add_btiT: step out of range: " << step; std::cout << "trigger not added!" << std::endl; return; } if(!config()->usedBti(pos)) { if(config()->debug()==4){ std::cout << "DTTracoChip::add_btiT: position: " << pos; std::cout << "has disconnected bti" << std::endl; } return; } // 091103 SV: acceptances are taken from geometry if useAcceptParam()=false // otherwise cuts based on LL,LH,CL,CH,RL,RH taken from configuration are applied in TracoCard::loadTraco if(_card->useAcceptParamFlag()==false) { // check K inside acceptance if(btitrig->K()<_PSIMIN[pos-1] || btitrig->K()>_PSIMAX[pos-1] ) { if(config()->debug()>1){ std::cout << "In TRACO num. " << number() << " BTI trig. in pos " << pos << " outside K acceptance ("; std::cout << _PSIMIN[pos-1] << "-->"; std::cout << _PSIMAX[pos-1] << ") - Not added" << std::endl; } return; } } // Store trigger candidate if(pos<=DTConfigTraco::NBTITC){ _innerCand[step-DTConfigTraco::NSTEPF].push_back( DTTracoCand(this,btitrig,pos,step)); } else { _outerCand[step-DTConfigTraco::NSTEPF].push_back( DTTracoCand(this,btitrig,pos,step)); } // Fill array for BX LTS if(btitrig->code()==8){ for(int is=step-4;is<step;is++){ // set flag for 4 previous BX if(is>0&&is<=DTConfigTraco::NSTEPL) _bxlts.set(is); } //SV 1/IV/04 BUG FIX _bxlts.set(step+1); // Debugging if(config()->debug()==4) for(int is=0;is<DTConfigTraco::NSTEPL;is++) std::cout<<"_bxlts["<<is<<"]="<<_bxlts.element(is)<<std::endl; } // Debugging if(config()->debug()>1){ std::cout << "BTI Trigger added at step " << step; std::cout << " to TRACO " << _id.traco() << " at position " << pos << std::endl; btitrig->print(); } // End debugging }
void DTTracoChip::addTrig | ( | int | step, |
DTTracoTrig * | tctrig | ||
) |
Add a TRACO trigger.
Definition at line 1059 of file DTTracoChip.cc.
References _tracotrig, config(), gather_cfg::cout, debug, DTConfig::NSTEPF, DTConfig::NSTEPL, DTTracoTrig::print(), and launcher::step.
Referenced by run(), and DTTracoCard::storeTrigger().
{ if(step<DTConfigTraco::NSTEPF||step>DTConfigTraco::NSTEPL){ std::cout << "DTTracoChip::addTrig: step out of range: " << step; std::cout << " trigger not added!" << std::endl; return; } _tracotrig[step-DTConfigTraco::NSTEPF].push_back(tctrig); // Debugging... if(config()->debug()==4){ std::cout << "DTTracoChip::addTrig: adding trigger:"<< std::endl; tctrig->print(); } // End debugging }
int DTTracoChip::AdjBtiLTSuppressed | ( | DTTracoCand * | candidate | ) | [private] |
Do suppression of LTRIG on adjacent TRACO.
Definition at line 554 of file DTTracoChip.cc.
References _flag, DTTracoCand::BtiTrig(), DTBtiTrigData::code(), BitArray< N >::element(), DTConfig::NBTITC, DTConfig::NSTEPF, DTTracoCand::position(), and DTTracoCand::step().
Referenced by bestCand().
{ // If requested, do suppression of LTRIG on adjacent BTI -> obsolete! //if(!(config()->adjBtiLts()) && candidate->BtiTrig()->code()<8) { //SV: Ltrig always suppressed in hardware if Htrig in adj traco! if(candidate->BtiTrig()->code()<8) { if( _flag[candidate->step()-DTConfigTraco::NSTEPF].element(3) && candidate->position()==1 ) return 1; if( _flag[candidate->step()-DTConfigTraco::NSTEPF].element(4) && candidate->position()==DTConfigTraco::NBTITC ) return 1; if( _flag[candidate->step()-DTConfigTraco::NSTEPF].element(5) && candidate->position()== DTConfigTraco::NBTITC+1) return 1; if( _flag[candidate->step()-DTConfigTraco::NSTEPF].element(6) && candidate->position()==DTConfigTraco::NBTITC*4 ) return 1; } //} return 0; }
DTTracoCand * DTTracoChip::bestCand | ( | int | itk, |
std::vector< DTTracoCand > & | tclist | ||
) | [private] |
Get the best inner/outer candidate.
Definition at line 471 of file DTTracoChip.cc.
References AdjBtiLTSuppressed(), config(), gather_cfg::cout, DTConfigTraco::debug(), i, AlCaHLTBitMon_ParallelJobs::p, DTConfigTraco::prefHtrig(), and DTConfigTraco::sortKascend().
Referenced by run().
{ // Return if no candidates if(tclist.size()<1) return 0; // stl function: sort in Ktc ascending or descending order according // to user request comparing by default with user-defined < // NB don't reverse if candidates are two with same K stable_sort( tclist.begin(),tclist.end() ); //0=K ascending, 1=K descending if(config()->sortKascend(itk) && !(tclist.size()==2 && tclist[0].K()==tclist[1].K()) ) { reverse( tclist.begin(),tclist.end() ); if(config()->debug()==4) std::cout << "Reversing order of sorted candidate list..." << std::endl; } /* if(!config()->sortKascend(itk)){ stable_sort( tclist.begin(),tclist.end(),DTTracoCand::closer ); } else { stable_sort( tclist.begin(),tclist.end(),DTTracoCand::wider ); } */ // debugging... if(config()->debug()==4){ std::cout << "DTTracoChip::findBest - Looking for track number " << itk+1 << std::endl ; std::cout << "Sorted std::vector of usable track candidates is:" << std::endl; int i = 1; for(std::vector<DTTracoCand>::iterator p=tclist.begin();p<tclist.end();p++){ if((*p).usable()){ std::cout << " DTTracoChip Candidate # " << i++; (*p).print(); } } std::cout << "--------------------------------------------------" << std::endl; } // end debugging // return the best candidate int i=0; DTTracoCand* bestltrig = 0; std::vector<DTTracoCand>::iterator p; for ( p = tclist.begin(); p < tclist.end(); ++p ) { i++; // candidate must be usable and not suppressed by LTS if(AdjBtiLTSuppressed(&(*p))) if(config()->debug()==4) std::cout << "Candidate # " << i << " supp. because next to H in adiacent Tracos" << std::endl; if ( (*p).usable() && !AdjBtiLTSuppressed(&(*p)) ) { // check if preference to HTRIG is set and return first trigger if( !config()->prefHtrig(itk) ) return &(*p); if( (*p).BtiTrig()->code()==8 ) return &(*p); if( bestltrig==0 ) bestltrig=&(*p); } } return bestltrig; }
int DTTracoChip::BTIC | ( | ) | const [inline] |
BTIC parameter.
Definition at line 127 of file DTTracoChip.h.
References _btic.
Referenced by DTTracoCand::DTTracoCand(), storeCorr(), and storeUncorr().
{ return _btic; }
void DTTracoChip::calculateAngles | ( | DTTracoTrig * | tct | ) |
Calculate trigger angles.
Definition at line 1210 of file DTTracoChip.cc.
References _card, _geom, _lutsCCB, DTTracoTrig::ChamberId(), DTGeomSupplier::CMSPosition(), config(), gather_cfg::cout, DTTracoTrig::data(), debug, Lut::get_k(), Lut::get_x(), DTTracoTrig::K(), DTTracoCard::localDirection(), DTTracoCard::lutFromDBFlag(), M_PI, PV3DBase< T, PVType, FrameType >::phi(), DTTrigGeom::phiCh(), DTTracoTrig::print(), DTTracoTrigData::qdec(), DTConfigTraco::RESOLPSI, DTConfigTraco::RESOLPSIR, DTChamberId::sector(), DTTracoTrig::setAngles(), DTChamberId::wheel(), PV3DBase< T, PVType, FrameType >::x(), DTTracoTrig::X(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by storeCorr(), and storeUncorr().
{ int ipsi=0; int iphir=0; int idpsir=0; /* obsolete //TB 2004 luts formula from minicrate CCB if( config()->trigSetupGeom()==2 ){ ipsi = _lutsCCB->get_k( (tct->K()+511) ); int flag = 0; int qual=tct->data().qdec(); if(qual==3 || qual==1) //case 0:outer flag=0; if(qual==2 || qual==0) //case 1:inner flag=1; if(qual==6 || qual==5 || qual==4) //case 2:correlated flag=2; iphir = _lutsCCB->get_x( (tct->X()+512*flag) ); idpsir = ipsi - iphir/8; } //TB 2003 luts data format if( config()->trigSetupGeom()==1 ){ //compute bending angles of traco output with lookup tables //SV TB2003: first trigger board isn't connected; ipsi = _luts->getPsi(tct->K()); int flag = 0; int qual=tct->data().qdec(); if(qual==3 || qual==1) //case 0:outer flag=0; if(qual==2 || qual==0) //case 1:inner flag=1; if(qual==6 || qual==5 || qual==4) //case 2:correlated flag=2; iphir = _luts->getPhiRad( tct->X(), flag); idpsir = _luts->getBendAng( tct->X(), tct->K(), flag); } */ // 091030 SV angles computed from DB lut parameters if( _card->lutFromDBFlag()==1 ) { ipsi = _lutsCCB->get_k( (tct->K()+512) ); int flag = 0; int qual=tct->data().qdec(); if(qual==3 || qual==1) //case 0:outer flag=0; if(qual==2 || qual==0) //case 1:inner flag=1; if(qual==6 || qual==5 || qual==4) //case 2:correlated flag=2; iphir = _lutsCCB->get_x( (tct->X()+512*flag) ); idpsir = ipsi - iphir/8; } else // compute angles from CMSSW geometry //if( config()->trigSetupGeom()==0 ) { DTTracoTrigData td = tct->data(); // psi // float fpsi = atan( (float)(tct->K()) * _geom->cellPitch() / // (_geom->distSL() * config()->ST()) ); float fpsi = atan( _card->localDirection(&td).x() / // e.g. x>0 and _card->localDirection(&td).z() ); // z<0 => fpsi<0 // Change sign in case of wheel<0 or wheel==0 and station == 1,4,5,8,9,12 int mywh = tct->ChamberId().wheel(); if (mywh<0 || (mywh==0 && (tct->ChamberId().sector()%4)<2)) fpsi = -fpsi; fpsi*=DTConfigTraco::RESOLPSI; if(fpsi<=0) fpsi-=1.0; ipsi = (int)fpsi; // if outside range set to lower edge if( ipsi>= DTConfigTraco::RESOLPSI || ipsi< -DTConfigTraco::RESOLPSI ) ipsi=-DTConfigTraco::RESOLPSI; // psi_r float fpsir = _card->CMSPosition(&td).phi()-_geom->phiCh(); if(fpsir<-M_PI) fpsir+=M_PI*2; if(fpsir>M_PI) fpsir-=M_PI*2; fpsir*=DTConfigTraco::RESOLPSIR; if(fpsir<=0) fpsir-=1.0; iphir = (int)fpsir; // if outside range set to lower edge if( iphir>= DTConfigTraco::RESOLPSIR/2 || iphir <-DTConfigTraco::RESOLPSIR/2 ) iphir=-DTConfigTraco::RESOLPSIR/2; // Delta(psi_r) int dpsir = (iphir*DTConfigTraco::RESOLPSI) / DTConfigTraco::RESOLPSIR; idpsir = ipsi-dpsir; // if outside range set to lower edge if(idpsir>= DTConfigTraco::RESOLPSI || idpsir <-DTConfigTraco::RESOLPSI ) idpsir=-DTConfigTraco::RESOLPSI; } tct->setAngles(ipsi,iphir,idpsir); // debugging if(config()->debug()==4){ std::cout << "DTTracoChip::calculateAngles :" << std::endl; tct->print(); std::cout << std::dec << "K = " << tct->K() << " X = " << tct->X(); std::cout << " ipsi = " << ipsi << " iphir = " << iphir; std::cout << " idpsir = " << idpsir << std::endl; if( _card->lutFromDBFlag()==1 ) std::cout << "Angles are calculated from LUT parameters from DB!" << std::endl; }// end debugging }
void DTTracoChip::clear | ( | void | ) |
Clear.
Definition at line 237 of file DTTracoChip.cc.
References _bxlts, _flag, _innerCand, _outerCand, _tracotrig, begin, DTConfig::NSTEPF, DTConfig::NSTEPL, p1, and BitArray< N >::zero().
Referenced by ~DTTracoChip().
{ std::vector<DTTracoTrig*>::iterator p1; for(int is=0;is<DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1;is++){ for(p1=_tracotrig[is].begin();p1<_tracotrig[is].end();p1++){ delete (*p1); } _tracotrig[is].clear(); _innerCand[is].clear(); _outerCand[is].clear(); _flag[is].zero(); } _bxlts.zero(); }
GlobalPoint DTTracoChip::CMSPosition | ( | ) | const [inline] |
Position in CMS frame.
Definition at line 154 of file DTTracoChip.h.
References _geom, _id, and DTTrigGeom::CMSPosition().
Referenced by DTTracoChip().
{ return _geom->CMSPosition(_id); }
DTConfigTraco* DTTracoChip::config | ( | void | ) | const [inline] |
Configuration set.
New Configuration set
Definition at line 117 of file DTTracoChip.h.
References _config.
Referenced by add_btiT(), addTrig(), bestCand(), calculateAngles(), DoAdjBtiLts(), DTTracoChip(), insideAngWindow(), run(), setFlag(), setPV(), setTracoAcceptances(), storeCorr(), and storeUncorr().
{ return _config; }
int DTTracoChip::DD | ( | ) | const [inline] |
void DTTracoChip::DoAdjBtiLts | ( | DTTracoCand * | candidate, |
std::vector< DTTracoCand > & | tclist | ||
) | [private] |
Do suppression of LTRIG on BTI close to selected HTRIG.
Definition at line 532 of file DTTracoChip.cc.
References abs, DTBtiTrigData::btiNumber(), DTTracoCand::BtiTrig(), DTBtiTrigData::code(), config(), gather_cfg::cout, debug, and AlCaHLTBitMon_ParallelJobs::p.
Referenced by run().
{ // If requested, do suppression of LTRIG on BTI close to selected HTRIG in same traco //if(!(config()->adjBtiLts()) && candidate->BtiTrig()->code()==8) { // SV this is done always, not parametrized !! if(candidate->BtiTrig()->code()==8) { std::vector<DTTracoCand>::iterator p; for(p=tclist.begin();p<tclist.end();p++){ if( (*p).BtiTrig()->code()<8 && abs( (*p).BtiTrig()->btiNumber() - candidate->BtiTrig()->btiNumber() ) < 2 ) { (*p).setUsed(); if(config()->debug()==4){ std::cout << "Candidate :"; (*p).print(); std::cout << "Suppressed because adiacent to H trig" <<std::endl; } // end debug } // end if } // end candidate loop } // end if H }
int DTTracoChip::edgeBTI | ( | int | step, |
int | io, | ||
int | lr | ||
) | const |
flags for HTRIG in edge BTI
Definition at line 1169 of file DTTracoChip.cc.
References _innerCand, _outerCand, begin, gather_cfg::cout, DTConfig::NBTITC, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, findQualityFiles::size, and launcher::step.
Referenced by setFlag().
{ if(step<DTConfigTraco::NSTEPF||step>DTConfigTraco::NSTEPL){ std::cout << "DTTracoChip::edgeBTI: step out of range: " << step; std::cout << " 0 returned!" << std::endl; return 0; } // // inner supl ==> io=1; outer supl ==> io=2 |21 | | 22| // right edge ==> rl=1; left edge ==> rl=2 |11 12| // std::vector<DTTracoCand>::const_iterator p; if(io==1){ if(_innerCand[step-DTConfigTraco::NSTEPF].size()>0) { // SV 24/IX/03 fix: only HTRIG accepted for(p=_innerCand[step-DTConfigTraco::NSTEPF].begin(); p<_innerCand[step-DTConfigTraco::NSTEPF].end(); p++){ if(lr==1 && (*p).position()==1 && (*p).BtiTrig()->code()==8 ) return 1; if(lr==2 && (*p).position()==DTConfigTraco::NBTITC && (*p).BtiTrig()->code()==8 ) return 1; } } } else { if(_outerCand[step-DTConfigTraco::NSTEPF].size()>0) { for(p=_outerCand[step-DTConfigTraco::NSTEPF].begin(); p<_outerCand[step-DTConfigTraco::NSTEPF].end(); p++){ //SV: is the following correct???FIX if using _card to set _flag //if(lr==1 && (*p).position()==DTConfigTraco::NBTITC+1)return 1; //or pos=8?? //if(lr==2 && (*p).position()==DTConfigTraco::NBTITC*4)return 1; //or pos=13?? //SV 24/IX/03 fix if(lr==1 && (*p).position()==DTConfigTraco::NBTITC*3+1 && (*p).BtiTrig()->code()==8 ) return 1; if(lr==2 && (*p).position()==DTConfigTraco::NBTITC*2 && (*p).BtiTrig()->code()==8 ) return 1; } } } return 0; }
DTTrigGeom* DTTracoChip::geom | ( | ) | const [inline] |
Return trigger geometry.
Definition at line 99 of file DTTracoChip.h.
References _geom.
Referenced by run().
{ return _geom; }
int DTTracoChip::IBTIOFF | ( | ) | const [inline] |
IBTIOFF parameter.
Definition at line 130 of file DTTracoChip.h.
References _ibtioff.
Referenced by DTTracoCand::DTTracoCand().
{ return _ibtioff; }
DTTracoId DTTracoChip::id | ( | void | ) | const [inline] |
Return TRACO id.
Definition at line 102 of file DTTracoChip.h.
References _id.
Referenced by DTTracoTrig::DTTracoTrig().
{ return _id; }
int DTTracoChip::insideAngWindow | ( | DTTracoTrig * | tctrig | ) | const [private] |
Check if a trigger is inside the angular acceptance window.
Definition at line 1335 of file DTTracoChip.cc.
References BitArray< N >::assign(), DTConfigTraco::BendingAngleCut(), config(), gather_cfg::cout, debug, DTTracoTrig::DeltaPsiR(), BitArray< N >::element(), BitArray< N >::read(), and BitArray< N >::twoComplement().
Referenced by storeCorr(), and storeUncorr().
{ // check angular window for every station type // return 1 for accept, 0 for reject // abs value of bending angle is 9 bits BitArray<10> bendAng; bendAng.assign(0,10,tctrig->DeltaPsiR()); //bendAng.print(); if(bendAng.element(9)) //negativo! bendAng.twoComplement(); int bendAngINT = bendAng.read(0,9); //std::cout<<"abs bend angle int ="<< bendAngINT <<std::endl; if( config()->BendingAngleCut()!= -1 && bendAngINT > 2*(config()->BendingAngleCut())) { int absBendAng = tctrig->DeltaPsiR() & 0x1FF; if(config()->debug()==4) std::cout << "Attention: abs(bendAng)=" << absBendAng << " > " << 2*config()->BendingAngleCut() <<"!! reject trigger"<<std::endl; return 0; } return 1; }
int DTTracoChip::KRad | ( | ) | const |
K par of the radial angle of corr center referred to plane sl.
Definition at line 1144 of file DTTracoChip.cc.
Referenced by DTTracoCand::DTTracoCand(), DTTracoChip(), and DTTracoCand::print().
LocalPoint DTTracoChip::localPosition | ( | ) | const [inline] |
Position in chamber frame.
Definition at line 151 of file DTTracoChip.h.
References _geom, _id, and DTTrigGeom::localPosition().
{ return _geom->localPosition(_id); }
int DTTracoChip::nTrig | ( | int | step | ) | const |
Return the number of trigger candidates.
Definition at line 1078 of file DTTracoChip.cc.
References _tracotrig, gather_cfg::cout, DTConfig::NSTEPF, DTConfig::NSTEPL, and launcher::step.
Referenced by run(), and DTTracoCard::runTRACO().
{ if(step<DTConfigTraco::NSTEPF||step>DTConfigTraco::NSTEPL){ std::cout << "DTTracoChip::nTrig: step out of range: " << step; std::cout << " 0 returned!" << std::endl; return 0; } return _tracotrig[step-DTConfigTraco::NSTEPF].size(); }
int DTTracoChip::number | ( | ) | const [inline] |
Return TRACO number.
Definition at line 93 of file DTTracoChip.h.
References _id, and DTTracoId::traco().
Referenced by add_btiT().
DTTracoChip & DTTracoChip::operator= | ( | const DTTracoChip & | traco | ) |
Assignment operator.
Definition at line 205 of file DTTracoChip.cc.
References _bxlts, _card, _flag, _geom, _id, _innerCand, _outerCand, _tracotrig, i, DTConfig::NBTITC, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, and p1.
{ if(this != &traco){ _geom = traco._geom; _id = traco._id; _card = traco._card; int i=0; for(i=0;i<DTConfigTraco::NSTEPL - DTConfigTraco::NSTEPF;i++) { _innerCand[i].reserve(DTConfigTraco::NBTITC); std::vector<DTTracoCand>::const_iterator p; for(p=traco._innerCand[i].begin();p<traco._innerCand[i].end();p++) { _innerCand[i].push_back(*p); } _outerCand[i].reserve(3*DTConfigTraco::NBTITC); for(p=traco._outerCand[i].begin();p<traco._outerCand[i].end();p++) { _outerCand[i].push_back(*p); } _tracotrig[i].reserve(2); std::vector<DTTracoTrig*>::const_iterator p1; for(p1=traco._tracotrig[i].begin();p1<traco._tracotrig[i].end();p1++) { _tracotrig[i].push_back(*p1); } } _bxlts = traco._bxlts; for(int is=0;is<DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1;is++){ _flag[is] = traco._flag[is]; } } return *this; }
int DTTracoChip::ovlFlag | ( | int | step | ) | [inline] |
return overlap flag
Definition at line 160 of file DTTracoChip.h.
References _flag, BitArray< N >::element(), and DTConfig::NSTEPF.
{ return _flag[step-DTConfigTraco::NSTEPF].element(1);}
float DTTracoChip::psiRad | ( | int | sl = 0 | ) | const |
Radial angle of correlator center in mrad referred to plane sl.
Definition at line 1122 of file DTTracoChip.cc.
Referenced by DTTracoChip().
{ /* // Radial angle of correlator center in mrad in CMS frame LocalPoint p = localPosition(); float x = p.x(); float y = p.y(); float z = p.z(); if (sl==1) { z -= 0.5 * _geom->distSL(); } else if (sl==3) { z += 0.5 * _geom->distSL(); } float fpsir = _geom->stat()->toGlobal(LocalPoint(x,y,z)).phi()- _geom->phiCh(); if(fpsir<-M_PI)fpsir+=M_PI*2; if(fpsir>M_PI)fpsir-=M_PI*2; return fpsir*1000; */ return 0.0; }
void DTTracoChip::raiseOverlap | ( | int | step | ) |
raise overlap flag
Definition at line 419 of file DTTracoChip.cc.
References _flag, DTConfig::NSTEPF, and BitArray< N >::set().
Referenced by run().
{ _flag[step-DTConfigTraco::NSTEPF].set(1); //overlap flag raised _flag[step-DTConfigTraco::NSTEPF-1].set(2); //mark II rej. }
void DTTracoChip::run | ( | void | ) |
Run TRACO algorithm.
Definition at line 255 of file DTTracoChip.cc.
References _flag, _id, _innerCand, _outerCand, _tracotrig, addTrig(), bestCand(), config(), gather_cfg::cout, debug, DTConfigTraco::debug(), DoAdjBtiLts(), alignCSCRings::e, BitArray< N >::element(), geom(), SurfaceOrientation::inner, DTConfig::NBTITC, DTConfig::NSTEPF, DTConfig::NSTEPL, nTrig(), SurfaceOrientation::outer, DTTracoCand::print(), DTTracoTrig::print(), raiseOverlap(), BitArray< N >::set(), setFlag(), setPV(), DTTracoCand::setUsed(), findQualityFiles::size, storeCorr(), storeUncorr(), DTTracoId::traco(), and trigger().
Referenced by DTTracoCard::runTRACO().
{ // Debugging... if(config()->debug()>1){ std::cout << "DTTracoChip::run: Processing TRACO " << _id.traco() << std::endl; } // End debugging int maxtc = static_cast<int>(ceil( float(geom()->nCell(1)) / float(DTConfigTraco::NBTITC) )); if( _id.traco()<1 || _id.traco()>maxtc ) { if(config()->debug()==4) std::cout << "DTTracoChip::run: wrong TRACO number " << _id.traco() << std::endl; return; } // Loop on step for(int is=DTConfigTraco::NSTEPF; is<=DTConfigTraco::NSTEPL;is++) { if(config()->debug()>1){ std::cout << "\n STEP: " << is << std::endl; std::cout << " ================" << std::endl; } // skip if no cand. at this step if(_innerCand[is-DTConfigTraco::NSTEPF].size()<1 && _outerCand[is-DTConfigTraco::NSTEPF].size()<1 ) continue; // Debugging... if(config()->debug()==4){ std::cout << " --> " << _innerCand[is-DTConfigTraco::NSTEPF].size()+ _outerCand[is-DTConfigTraco::NSTEPF].size(); std::cout << " candidates " << std::endl; } // End debugging // Multiple trigger detection between consecutive TRACO's setFlag(is); //check if there is a H in bx for LVALIDIFH flag if(config()->LVALIDIFH()){ for(unsigned int e=0; e<_innerCand[is-DTConfigTraco::NSTEPF].size(); e++) { if(_innerCand[is-DTConfigTraco::NSTEPF][e].BtiTrig()->code()==8){ _flag[is-DTConfigTraco::NSTEPF].set(9); break; } } for(unsigned int e=0; e<_outerCand[is-DTConfigTraco::NSTEPF].size(); e++) { if(_outerCand[is-DTConfigTraco::NSTEPF][e].BtiTrig()->code()==8){ _flag[is-DTConfigTraco::NSTEPF].set(9); break; } } } // Loop over first/second tracks //for(int itk=0; itk < DTConfigTraco::NMAXCAND; itk++){ // FIX this hardcoded 2!! for(int itk=0; itk < 2; itk++){ // Get the best inner and outer segments if(config()->debug()==4) std::cout << "Inner:" << std::endl; DTTracoCand* inner = bestCand(itk,_innerCand[is-DTConfigTraco::NSTEPF]); if(config()->debug()==4) std::cout << "Outer:" << std::endl; DTTracoCand* outer = bestCand(itk,_outerCand[is-DTConfigTraco::NSTEPF]); //debug if(config()->debug()>1){ if(inner || outer) std::cout<<"Best candidates for track " << itk+1 << " are:"<<std::endl; if(inner){ std::cout<<"inner->"; inner->print(); } if(outer){ std::cout<<"outer->"; outer->print(); } } // Skip to next step if no suitable candidates found if(inner==0&&outer==0) break; // suppression of LTRIG on BTI close to selected HTRIG // SV 24/IX/03 : AND suppression of LTRIG close to H in adiacent Traco // SV 31/III/03 : test : only if htprf is off--> NO, it's worse // if( config()->prefHtrig(0) && config()->prefHtrig(1) ){ if(inner){ DoAdjBtiLts( inner, _innerCand[is-DTConfigTraco::NSTEPF] ); } if(outer){ DoAdjBtiLts( outer, _outerCand[is-DTConfigTraco::NSTEPF] ); } //} // set candidates unusable by further steps if(inner)inner->setUsed(); if(outer)outer->setUsed(); // Create a new TRACO trigger with preview for TS DTTracoTrig* tct = setPV(itk,inner,outer); // skip if no TRACO trigger has been created with this trigger if(!tct) break; // something nasty happened. Go to next step // try a correlation between segments int stored = 0; if(inner && outer) { stored = storeCorr(tct,inner,outer,itk); } if (!stored) { // No suitable candidate found or no correlation possible // Fill single trigger stored = storeUncorr(tct,inner,outer,itk); } // if trigger has been filled store it in TRACO, otherway delete it if (stored) { addTrig(is,tct); } else { delete tct; //break;-> II track is computed even if no I track found... } } // end loop on first/second track // Inhibit second track at previous bunch crossing if(config()->debug()==4) std::cout<<"Checking overlap I-II track..." <<std::endl; if(_tracotrig[is-DTConfigTraco::NSTEPF].size()>0 && is>DTConfigTraco::NSTEPF && (_tracotrig[is-DTConfigTraco::NSTEPF])[0]->isFirst() ) { //I track at bx if(nTrig(is-1)>0) { //there is a track at bx-1 if( !(trigger(is-1,1)->isFirst()) || //trig 1 is II track ( nTrig(is-1)==2 && !(trigger(is-1,2)->isFirst()) )) { //trig 2 is II track raiseOverlap(is); if(config()->debug()==4){ std::cout << "II track at step " << std::hex << is-1 <<std::dec << "marked rej."<< std::endl; std::cout << "I track overlap flag at step " << std::hex << is << std::dec << " setted" << std::endl; } } } } //debug... for(int isd=0;isd<=DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1;isd++) if(config()->debug()==4){ std::cout << "overlap flag step = " << isd+DTConfigTraco::NSTEPF << " " << _flag[isd].element(1) << std::endl; } // debugging... if(config()->debug()>0){ if(nTrig(is)>0) { for(int cc=1;cc<=nTrig(is);cc++){ trigger(is,cc)->print(); } } }// end debugging }// end loop on step }
int DTTracoChip::sector | ( | void | ) | const [inline] |
Return sector number.
Definition at line 111 of file DTTracoChip.h.
References _geom, and DTTrigGeom::sector().
void DTTracoChip::setFlag | ( | int | step, |
int | ext = 0 |
||
) |
Set flags for multiple trigger detection between cons. TRACO's.
Definition at line 426 of file DTTracoChip.cc.
References _card, _flag, _id, config(), gather_cfg::cout, debug, edgeBTI(), BitArray< N >::element(), DTTracoCard::getTRACO(), i, DTConfig::NSTEPF, BitArray< N >::set(), and DTTracoId::traco().
Referenced by run().
{ if(ext==0){ //this is the original: flags from card DTTracoChip* prevTraco = _card->getTRACO(_id.traco()-1); if(prevTraco!=0){ if(prevTraco->edgeBTI(step,1,2)) _flag[step-DTConfigTraco::NSTEPF].set(3); if(prevTraco->edgeBTI(step,2,2)) _flag[step-DTConfigTraco::NSTEPF].set(5); } DTTracoChip* nextTraco = _card->getTRACO(_id.traco()+1); if(nextTraco!=0){ if(nextTraco->edgeBTI(step,1,1)) _flag[step-DTConfigTraco::NSTEPF].set(4); if(nextTraco->edgeBTI(step,2,1)) _flag[step-DTConfigTraco::NSTEPF].set(6); } } else{ //SV III/03: flags from input EXT: only for testing purpose for(int i=0;i<6;i++){ int ibit = ext >> i; if(ibit & 0x01) // bit i+1 -> flag 3,4,5,6 : IL,IR,OL,OR _flag[step-DTConfigTraco::NSTEPF].set(i+1 + 2); } } //debug: if(config()->debug()==4){ std::cout << "Flags set for bx=" << step << std::endl; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(1)<< " "; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(2)<< " "; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(3)<< " "; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(4)<< " "; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(5)<< " "; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(6)<< " "; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(7)<< " "; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(8)<< " "; std::cout << _flag[step-DTConfigTraco::NSTEPF].element(9)<< " "<<std::endl; } //end debugging }
DTTracoTrig * DTTracoChip::setPV | ( | int | itk, |
DTTracoCand * | inner, | ||
DTTracoCand * | outer | ||
) | [private] |
Set the preview for a trigger.
Definition at line 573 of file DTTracoChip.cc.
References DTBtiTrigData::btiSL(), DTTracoCand::BtiTrig(), DTBtiTrigData::code(), config(), gather_cfg::cout, debug, first, SurfaceOrientation::inner, DTTracoCand::K(), SurfaceOrientation::outer, DTTracoCand::position(), DTConfigTraco::prefInner(), DTTracoCand::print(), DTTracoTrig::setPV(), and DTTracoCand::step().
{ // debugging... if(config()->debug()==4){ std::cout << "DTTracoChip::setPV called for candidates : " << std::endl; if(inner)inner->print(); if(outer)outer->print(); std::cout << "--------------------------------------------------" << std::endl; } //end debugging // first or second track. This is tricky: // itk=0 menas first track ==> first= true=1 // itk=1 menas second track ==> first=false=0 int first = (itk==0) ? 1 : 0; //preview selector: the same as priority selector I !! // select which of the inner/outer segments should be used DTTracoCand* candidate=0; if(inner!=0&&outer!=0) { // if(config()->prefHtrig(itk)){ // ---> BUG! selection is ALWAYS for H trigs // ---> fixed by Sara Vanini 25/III/03 if(inner->BtiTrig()->code()==8 && outer->BtiTrig()->code()<8 ){ candidate=inner; } else if(inner->BtiTrig()->code()<8 && outer->BtiTrig()->code()==8){ candidate=outer; } else { //for same quality tracks, pref. to in/out selection if(!config()->prefInner(itk)) { candidate=inner; } else { candidate=outer; } } // } //end if(config()->prefHtrig... /* else { if(!config()->prefInner(itk)) { candidate=inner; } else { candidate=outer; } } */ } else if(inner==0&&outer!=0) { candidate=outer; } else if(inner!=0&&outer==0) { candidate=inner; } else { return 0; // no candidates } // create new trigger with this candidate DTTracoTrig* tct = new DTTracoTrig(this, candidate->step()); // store preview for TS int cod = candidate->BtiTrig()->code(); if(candidate->BtiTrig()->btiSL()==1) cod *= 10; // k output is 5 bits!!! SV int K=candidate->K(); if(K>31) K-=32; int ioflag = 0; if(candidate->position()>4) ioflag = 1; tct->setPV(first, cod, K, ioflag); // this is already BTI_K-KRAD if(config()->debug()==4){ std::cout << "Selected candidate stored for preview is: "; candidate->print(); } return tct; }
void DTTracoChip::setPV | ( | int | step, |
int | ntrk, | ||
int | code, | ||
int | K | ||
) |
Set the preview values for a TRACO trigger.
Referenced by run().
void DTTracoChip::setTracoAcceptances | ( | ) | [private] |
Compute traco chip acceptances.
Definition at line 1361 of file DTTracoChip.cc.
References _geom, _PSIMAX, _PSIMIN, DTConfigTraco::BTIC(), DTTrigGeom::cellH(), DTTrigGeom::cellPitch(), config(), gather_cfg::cout, debug, DTTrigGeom::distSL(), h, i, DTConfig::NBTITC, and DTTrigGeom::phiSLOffset().
Referenced by DTTracoChip().
{ // Set K acceptances of DTTracoChip MT ports: Ktraco = Xinner - Xouter float h = _geom->cellH(); float pitch = _geom->cellPitch(); float distsl = _geom->distSL(); float K0 = config()->BTIC(); float shiftSL = _geom->phiSLOffset() / pitch * K0; // mt ports from orca geometry: this is always the case with new DTConfig //if(config_traco(tracoid)->trigSetupGeom() != 1){ { // Master Plane int i = 0; for(i=0;i<DTConfig::NBTITC;i++){ float Xin_min = (i + DTConfig::NBTITC) * K0 + shiftSL; float Xin_max = Xin_min + K0; float Xout_min = 0; float Xout_max = 3 * DTConfig::NBTITC * K0; _PSIMAX[i] = int( 2.*h/distsl * (Xin_max - Xout_min) + K0 + 1.01 ); _PSIMIN[i] = int( 2.*h/distsl * (Xin_min - Xout_max) + K0 ); } // Slave Plane for(i=0;i<3*DTConfig::NBTITC;i++){ float Xin_min = (DTConfig::NBTITC) * K0 + shiftSL; float Xin_max = 2. * DTConfig::NBTITC * K0 + shiftSL; float Xout_min = i * K0; float Xout_max = Xout_min + K0; _PSIMAX[DTConfig::NBTITC+i] = int( 2.*h/distsl * (Xin_max - Xout_min) + K0 + 1.01 ); _PSIMIN[DTConfig::NBTITC+i] = int( 2.*h/distsl * (Xin_min - Xout_max) + K0 ); } } // debugging if(config()->debug()==4){ //if(wheel()==2&&station()==3&§or()==1){ // only 1 chamber std::cout << "Acceptance of mt ports for offset (cell unit) " << _geom->phiSLOffset() / pitch << std::endl; for(int i=0;i<4*DTConfig::NBTITC;i++){ std::cout << "Port " << i+1 << " : "; std::cout << _PSIMIN[i] << " --> " << _PSIMAX[i] << std::endl; } //} }// end debugging }
int DTTracoChip::station | ( | ) | const [inline] |
Return station number.
Definition at line 108 of file DTTracoChip.h.
References _geom, and DTTrigGeom::station().
int DTTracoChip::storeCorr | ( | DTTracoTrig * | tctrig, |
DTTracoCand * | inner, | ||
DTTracoCand * | outer, | ||
int | tkn | ||
) | [private] |
Check correlation and store correlated trigger.
Definition at line 648 of file DTTracoChip.cc.
References _bxlts, _geom, abs, DTTracoTrig::addDTBtiTrig(), BTIC(), DTBtiTrigData::btiNumber(), DTTracoCand::BtiTrig(), calculateAngles(), DTTrigGeom::cellH(), DTBtiTrigData::code(), DTTracoTrig::code(), config(), gather_cfg::cout, debug, DTConfigTraco::debug(), DTTrigGeom::distSL(), BitArray< N >::element(), DTBtiTrigData::eq(), insideAngWindow(), DTTracoTrig::K(), DTBtiTrigData::K(), DTTracoCand::print(), DTTracoTrig::qdec(), DTTracoTrig::resetVar(), DTTracoTrig::setCodeIn(), DTTracoTrig::setCodeOut(), DTTracoTrig::setEqIn(), DTTracoTrig::setEqOut(), DTTracoTrig::setK(), DTTracoTrig::setPosIn(), DTTracoTrig::setPosOut(), DTTracoTrig::setPVCorr(), DTTracoTrig::setX(), edm::shift, DTTracoTrig::step(), DTConfigTraco::TcBxLts(), DTConfigTraco::TcKToll(), DTTracoCand::X(), and DTTracoTrig::X().
Referenced by run().
{ // Bunch crossing int is = tctrig->step(); // Debugging... if(config()->debug()==4){ std::cout << "DTTracoChip::storeCorr called with candidates: " << std::endl; if(inner)inner->print(); if(outer)outer->print(); std::cout << "--------------------------------------------------" << std::endl; } // End debugging //old orca shift definition float shift = 0.; //if( config()->trigSetupGeom()!=1 ) shift = (int)( _geom->distSL()/_geom->cellH() + 0.5 ); //else //shift = DD(); int kcor = 9999; int xcor = 0; int icor = 0; // Check correlation only if --> this cuts LL follow by H in next 4 bx // SV 1/IV/04 BUG FIX: this cuts LL preview also, traco outputs preview when LTS cut!!! //if( !config()->TcBxLts() || // BX LTS is not enabled or // !_bxlts.element(is) || // no HTRIG in next 4 BX or // inner->BtiTrig()->code()==8 || // inner track is HTRIG or // outer->BtiTrig()->code()==8 ){ // outer track is HTRIG //otherwise in and out trig are L, and necessary one is suppressed for LTS int xq1 = inner->X(); int xq2 = outer->X(); xcor = (xq2+xq1)/2; kcor = (xq1-xq2)+512; int kq1 = int(shift/2.) * (inner->BtiTrig()->K()-BTIC()) + 512; int kq2 = int(shift/2.) * (outer->BtiTrig()->K()-BTIC()) + 512; //int kd1 = abs(kcor/16-kq1/16); //int kd2 = abs(kcor/16-kq2/16); int kd1 = abs(kcor/16-kq1/16); int kd2 = abs(kcor/16-kq2/16); icor = kd1<=config()->TcKToll(tkn) && kd2<=config()->TcKToll(tkn) && xcor>0; // Debugging... if(config()->debug()==4){ std::cout << "*************************************************************"; std::cout << std::endl; std::cout << " shift = " << shift; std::cout << " xq1 = " << xq1; std::cout << " xq2 = " << xq2; std::cout << " xcor = " << xcor; std::cout << " kcor = " << kcor; std::cout << " kq1 = " << kq1; std::cout << " kq2 = " << kq2; std::cout << " kd1 = " << kd1; std::cout << " kd2 = " << kd2; std::cout << " icor = " << icor; std::cout << std::endl; std::cout << "*************************************************************"; std::cout << std::endl; }// End debugging //}//end if TcBxLts.... if(icor){ // correlation was successfull // set the preview correlation bit. It isn't reset if outside ang. window tctrig->setPVCorr(1); // set K and X tctrig->setK(kcor - 512); //std::cout<<"Set K " << kcor << " instead of " << kcor-512 << std::endl; //tctrig->setK(kcor); tctrig->setX(xcor); // set codes tctrig->setCodeIn( inner->BtiTrig()->code()); tctrig->setCodeOut(outer->BtiTrig()->code()); // set position mask //tctrig->setPosIn(inner->position()); //tctrig->setPosOut(outer->position()); //SV number of bti instead of position... tctrig->setPosIn( inner->BtiTrig()->btiNumber() ); tctrig->setPosOut( outer->BtiTrig()->btiNumber() ); //SV store also equation: pattern numbers are 1-32 tctrig->setEqIn( inner->BtiTrig()->eq() + 1 ); tctrig->setEqOut( outer->BtiTrig()->eq() + 1 ); // calculate psi, psi_r and Delta(psi_r) calculateAngles(tctrig); // check angular window for LL --> fixed by SV 27/III/03 --> NO, for all! //if( (tctrig->qdec()==4) && !insideAngWindow(tctrig)) { if( !insideAngWindow(tctrig) ) { // reset codes, K, X and angles tctrig->resetVar() ; } //SV 1/IV/04 BUG FIX: check LTS after angle cut... else if( tctrig->qdec()==4 && // cut only LL config()->TcBxLts() ){ // BX LTS is enabled or // reset codes, K, X and angles if(tkn==0 && _bxlts.element(is) ) // I track : there is H -4,+1 tctrig->resetVar() ; if(tkn==1 && _bxlts.element(is+1) ) // II track : there is H -4,+1 1 bx later tctrig->resetVar() ; } else { // set links to BTI triggers tctrig->addDTBtiTrig(inner->BtiTrig()); tctrig->addDTBtiTrig(outer->BtiTrig()); } // Debugging... if(config()->debug()>1){ std::cout << "*************************************************************"; std::cout << std::endl; std::cout << " Correlation was successfull: "; std::cout << std::endl; std::cout << " Code = " << tctrig->code(); std::cout << " K = " << tctrig->K(); std::cout << " X = " << tctrig->X(); std::cout << std::endl; std::cout << "*************************************************************"; std::cout << std::endl; } // End debugging } else { // Debugging... if(config()->debug()>1){ std::cout << "*************************************************************"; std::cout << std::endl; std::cout << " No correlation possible "; std::cout << std::endl; std::cout << "*************************************************************"; std::cout << std::endl; } // End debugging } return icor; }
int DTTracoChip::storeUncorr | ( | DTTracoTrig * | tctrig, |
DTTracoCand * | inner, | ||
DTTracoCand * | outer, | ||
int | tkn | ||
) | [private] |
Store uncorrelated trigger.
Definition at line 796 of file DTTracoChip.cc.
References _bxlts, _card, _flag, _geom, DTTracoTrig::addDTBtiTrig(), BTIC(), DTBtiTrigData::btiNumber(), DTBtiTrigData::btiSL(), DTTracoCand::BtiTrig(), calculateAngles(), DTTrigGeom::cellH(), DTBtiTrigData::code(), DTTracoTrig::code(), config(), gather_cfg::cout, debug, DTConfigTraco::debug(), DTTrigGeom::distSL(), BitArray< N >::element(), DTBtiTrigData::eq(), SurfaceOrientation::inner, insideAngWindow(), DTTracoTrig::K(), DTBtiTrigData::K(), DTConfigTraco::LVALIDIFH(), DTTSTheta::nHTrig(), DTConfig::NSTEPF, DTTSTheta::nTrig(), SurfaceOrientation::outer, DTConfigTraco::prefInner(), DTTracoCand::print(), DTTracoTrig::resetVar(), DTTracoTrig::setCodeIn(), DTTracoTrig::setCodeOut(), DTTracoTrig::setEqIn(), DTTracoTrig::setEqOut(), DTTracoTrig::setK(), DTTracoTrig::setPosIn(), DTTracoTrig::setPosOut(), DTTracoTrig::setPVCorr(), DTTracoCand::setUnused(), DTTracoTrig::setX(), edm::shift, DTConfigTraco::singleLenab(), DTConfigTraco::singleLflag(), DTTracoTrig::step(), DTConfigTraco::TcReuse(), DTTracoCard::TSTh(), DTTracoCand::X(), and DTTracoTrig::X().
Referenced by run().
{ // Bunch crossing int is = tctrig->step(); // Debugging... if(config()->debug()==4){ std::cout << "DTTracoChip::storeUncorr called with candidates: " << std::endl; if(inner)inner->print(); if(outer)outer->print(); std::cout << "--------------------------------------------------" << std::endl; } // End debugging // priority selector // select which of the inner/outer segments should be used // allow re-use of other segment according to configuration DTTracoCand* candidate=0; if(inner!=0&&outer!=0) { // if(config()->prefHtrig(tkn)){ // --> BUG: selector I preference is ALWAYS for H trig // fixed by Sara Vanini - 25/III/03 if(inner->BtiTrig()->code()==8 && outer->BtiTrig()->code()<8 ){ candidate=inner; //if(config()->TcReuse(1)) outer->setUnused(); // reusable } else if(inner->BtiTrig()->code()<8 && outer->BtiTrig()->code()==8){ candidate=outer; //if(config()->TcReuse(0)) inner->setUnused(); // reusable } else { //for the same quality triggers: if(!config()->prefInner(tkn)) { candidate=inner; //if(config()->TcReuse(1)) outer->setUnused(); // reusable } else { candidate=outer; //if(config()->TcReuse(0)) inner->setUnused(); // reusable } }//end else /* } else {//no Htrig preference if(!config()->prefInner(tkn)) { candidate=inner; if(config()->TcReuse(1)) outer->setUnused(); // reusable } else { candidate=outer; if(config()->TcReuse(0)) inner->setUnused(); // reusable } } */ } else if(inner==0&&outer!=0) { candidate=outer; } else if(inner!=0&&outer==0) { candidate=inner; } else { return 0; // no candidates } //SV *** FOR TESTBEAM OR TEST BENCH PURPOSE *** //theta trigger bin present(1) or absent(0) //int thTr = (_flag[is-DTConfigTraco::NSTEPF].element(7)) ? // _flag[is-DTConfigTraco::NSTEPF].element(7): // _flag[is-DTConfigTraco::NSTEPF].element(8); // priority selector II: accept or discard candidate according to masks: // ** LOW TRIGGERS if( candidate->BtiTrig()->code()<8 ) { // first check with LVALIDIFH: if 1, accept low if there is a H in traco at bx if(config()->LVALIDIFH() && _flag[is-DTConfigTraco::NSTEPF].element(9)){ if(config()->debug()>1) std::cout << "Low accepted because LVALIDIFH on...." << std::endl; } else {//LVALIDIFH==0 or there isn't H in traco in bx : check theta! //theta check if( !config()->singleLenab(tkn) ) { // LTF: single LTRIG not always en. Check cond.: if( config()->singleLflag(tkn)==1 || //always discarded ( config()->singleLflag(tkn)==2 && !(_card->TSTh()->nHTrig(is)) ) || ( config()->singleLflag(tkn)==0 && !(_card->TSTh()->nTrig(is)) ) ){ // SV --> for TESTS version // config()->singleLflag(tkn)==0 && thTr==0 || //only with theta trig. // config()->singleLflag(tkn)==2 && thTr==0 ){ //only with theta H trig (not hw) if(config()->debug()>1) std::cout << "Single low trigger discarded by preview and " << "priority selector for ltmsk!" <<std::endl; return 0;} // ^-------- trigger is suppressed and will not be stored }//end theta } //end else //REUSE : mark candidates reusable HERE! SV BUG FIX 6IV04 if(candidate==inner && config()->TcReuse(1) && outer) outer->setUnused(); if(candidate==outer && config()->TcReuse(0) && inner) inner->setUnused(); // LTS suppression if(config()->TcBxLts()){ if( (tkn==0 && _bxlts.element(is)) // I track : there is H -4,+1 || (tkn==1 && _bxlts.element(is+1)) ){ // II track : there is H -4,+1 1 bx later tctrig->resetVar() ; if(config()->debug()>1) std::cout << "Low trigger suppressed because H in next 4 bx "<< " and LTS flag on...." << std::endl; return 1; // trigger is suppressed but preview will be stored } }//end lts // } //end else } //Low trigs // Preview Htmsk not implemented: awaiting decision // --> implemented in priority selector by SV else { // HTRIG //if(config()->singleHflag(tkn)==1 && thTr==0 ) //this is for testing if( config()->singleHflag(tkn)==1 && !(_card->TSTh()->nTrig(is) ) ) return 0; // ^-----trigger is suppressed and will not be stored //mark candidates reusable HERE! SV BUG FIX 6IV04 if(candidate==inner && config()->TcReuse(1) && outer) outer->setUnused(); if(candidate==outer && config()->TcReuse(0) && inner) inner->setUnused(); } // set code, position, K and X float shift; //if(config()->trigSetupGeom()!=1 ) shift = (int)( _geom->distSL()/_geom->cellH() + 0.5 ); //else //shift = DD(); //SV 19/III/03 int kucor = (int)( 0.5*shift * (candidate->BtiTrig()->K()-BTIC()) ); tctrig->setK(kucor); tctrig->setX( candidate->X() ); // correlation wasn't successfull // set the preview correlation bit. tctrig->setPVCorr(0); if(candidate->BtiTrig()->btiSL()==1){ // inner track tctrig->setCodeIn(candidate->BtiTrig()->code()); tctrig->setCodeOut(0); //tctrig->setPosIn(candidate->position()); //SV number of bti instead of position... tctrig->setPosIn(candidate->BtiTrig()->btiNumber() ); tctrig->setPosOut(0); //SV store also equation tctrig->setEqIn( candidate->BtiTrig()->eq() + 1 ); tctrig->setEqOut( 0 ); } else { // outer track tctrig->setCodeIn(0); tctrig->setCodeOut(candidate->BtiTrig()->code()); tctrig->setPosIn(0); //SV number of bti instead of position... tctrig->setPosOut(candidate->BtiTrig()->btiNumber() ); //tctrig->setPosOut(candidate->position()); //SV store also equation tctrig->setEqIn( 0 ); tctrig->setEqOut( candidate->BtiTrig()->eq() + 1); } // coordinate converter LUT // calculate psi, psi_r and Delta(psi_r) calculateAngles(tctrig); // check angular window only for Low!! --> fixed SV 27/III/03--> NO, for all! //if( candidate->BtiTrig()->code() < 8 && !insideAngWindow(tctrig) ){ if( !insideAngWindow(tctrig) ){ // reset codes, K, X and angles tctrig->resetVar() ; if(config()->debug()>1) std::cout << "L rejected because outside angular window!" << std::endl; } else { // set links to BTI trigger tctrig->addDTBtiTrig(candidate->BtiTrig()); } // Debugging... if(config()->debug()>1){ std::cout << "*************************************************************"; std::cout << std::endl; std::cout << " Single trigger stored: "; std::cout << std::endl; std::cout << " Code = " << tctrig->code(); std::cout << " K = " << tctrig->K(); std::cout << " X = " << tctrig->X(); std::cout << std::endl; std::cout << "*************************************************************"; std::cout << std::endl; } // End debugging return 1; }
DTTracoTrig * DTTracoChip::trigger | ( | int | step, |
unsigned | n | ||
) | const |
Return the requested trigger.
Definition at line 1088 of file DTTracoChip.cc.
References _tracotrig, gather_cfg::cout, n, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, findQualityFiles::size, and launcher::step.
Referenced by run().
{ if(step<DTConfigTraco::NSTEPF||step>DTConfigTraco::NSTEPL){ std::cout << "DTTracoChip::trigger: step out of range: " << step; std::cout << " empty pointer returned!" << std::endl; return 0; } if(n<1 || n>_tracotrig[step-DTConfigTraco::NSTEPF].size()) { std::cout << "DTTracoChip::trigger: requested trigger doesn't exist: " << n; std::cout << " empty pointer returned!" << std::endl; return 0; } std::vector<DTTracoTrig*>::const_iterator p = _tracotrig[step-DTConfigTraco::NSTEPF].begin()+n-1; return *p; }
DTTracoTrigData DTTracoChip::triggerData | ( | int | step, |
unsigned | n | ||
) | const |
Return the data part of the requested trigger.
Definition at line 1105 of file DTTracoChip.cc.
References _tracotrig, gather_cfg::cout, n, DTConfig::NSTEPF, DTConfig::NSTEPL, AlCaHLTBitMon_ParallelJobs::p, findQualityFiles::size, and launcher::step.
Referenced by DTTracoCard::runTRACO().
{ if(step<DTConfigTraco::NSTEPF||step>DTConfigTraco::NSTEPL){ std::cout << "DTTracoChip::triggerData: step out of range: " << step; std::cout << " dummy trigger returned!" << std::endl; return DTTracoTrigData(); } if(n<1 || n>_tracotrig[step-DTConfigTraco::NSTEPF].size()) { std::cout << "DTTracoChip::trigger: requested trigger doesn't exist: " << n; std::cout << " dummy trigger returned!" << std::endl; return DTTracoTrigData(); } std::vector<DTTracoTrig*>::const_iterator p = _tracotrig[step-DTConfigTraco::NSTEPF].begin()+n-1; return (*p)->data(); }
int DTTracoChip::useSecondTrack | ( | int | step | ) | const |
a flag for a usable second track
Definition at line 1157 of file DTTracoChip.cc.
References _flag, gather_cfg::cout, BitArray< N >::element(), DTConfig::NSTEPF, DTConfig::NSTEPL, and launcher::step.
Referenced by DTTracoCard::runTRACO().
{ // return 1 if II track use is allow // return 0 if II track has been rejected if(step<DTConfigTraco::NSTEPF||step>DTConfigTraco::NSTEPL){ std::cout << "DTTracoChip::useSecondTrack: step out of range: " << step; std::cout << " 0 returned!" << std::endl; return 0; } return !(_flag[step-DTConfigTraco::NSTEPF].element(2)); }
int DTTracoChip::wheel | ( | ) | const [inline] |
Return wheel number.
Definition at line 105 of file DTTracoChip.h.
References _geom, and DTTrigGeom::wheel().
int DTTracoChip::_btic [private] |
Definition at line 199 of file DTTracoChip.h.
Referenced by BTIC(), and DTTracoChip().
BitArray<DTConfigTraco::NSTEPL+2> DTTracoChip::_bxlts [private] |
Definition at line 211 of file DTTracoChip.h.
Referenced by add_btiT(), clear(), DTTracoChip(), operator=(), storeCorr(), and storeUncorr().
DTTracoCard* DTTracoChip::_card [private] |
Definition at line 194 of file DTTracoChip.h.
Referenced by add_btiT(), calculateAngles(), DTTracoChip(), operator=(), setFlag(), storeUncorr(), and ~DTTracoChip().
DTConfigTraco* DTTracoChip::_config [private] |
Definition at line 196 of file DTTracoChip.h.
Referenced by config().
int DTTracoChip::_dd [private] |
Definition at line 201 of file DTTracoChip.h.
Referenced by DD().
BitArray<32> DTTracoChip::_flag[DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1] [private] |
Definition at line 217 of file DTTracoChip.h.
Referenced by AdjBtiLTSuppressed(), clear(), DTTracoChip(), operator=(), ovlFlag(), raiseOverlap(), run(), setFlag(), storeUncorr(), and useSecondTrack().
DTTrigGeom* DTTracoChip::_geom [private] |
Definition at line 191 of file DTTracoChip.h.
Referenced by calculateAngles(), CMSPosition(), DTTracoChip(), geom(), localPosition(), operator=(), sector(), setTracoAcceptances(), station(), storeCorr(), storeUncorr(), and wheel().
int DTTracoChip::_ibtioff [private] |
Definition at line 200 of file DTTracoChip.h.
Referenced by DTTracoChip(), and IBTIOFF().
DTTracoId DTTracoChip::_id [private] |
Definition at line 192 of file DTTracoChip.h.
Referenced by add_btiT(), CMSPosition(), DTTracoChip(), id(), localPosition(), number(), operator=(), run(), and setFlag().
std::vector<DTTracoCand> DTTracoChip::_innerCand[DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1] [private] |
Definition at line 204 of file DTTracoChip.h.
Referenced by add_btiT(), clear(), DTTracoChip(), edgeBTI(), operator=(), and run().
int DTTracoChip::_krad [private] |
Definition at line 198 of file DTTracoChip.h.
Referenced by DTTracoChip().
DTTracoLUTs* DTTracoChip::_luts [private] |
Definition at line 224 of file DTTracoChip.h.
Referenced by DTTracoChip().
Lut* DTTracoChip::_lutsCCB [private] |
Definition at line 226 of file DTTracoChip.h.
Referenced by calculateAngles(), DTTracoChip(), and ~DTTracoChip().
std::vector<DTTracoCand> DTTracoChip::_outerCand[DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1] [private] |
Definition at line 205 of file DTTracoChip.h.
Referenced by add_btiT(), clear(), DTTracoChip(), edgeBTI(), operator=(), and run().
int DTTracoChip::_PSIMAX[4 *DTConfig::NBTITC] [private] |
Definition at line 221 of file DTTracoChip.h.
Referenced by add_btiT(), and setTracoAcceptances().
int DTTracoChip::_PSIMIN[4 *DTConfig::NBTITC] [private] |
Definition at line 220 of file DTTracoChip.h.
Referenced by add_btiT(), and setTracoAcceptances().
std::vector<DTTracoTrig*> DTTracoChip::_tracotrig[DTConfigTraco::NSTEPL-DTConfigTraco::NSTEPF+1] [private] |
Definition at line 208 of file DTTracoChip.h.
Referenced by addTrig(), clear(), DTTracoChip(), nTrig(), operator=(), run(), trigger(), and triggerData().