Project CMSSW displayed by LXR CMSSW_14_2_X_2024-10-02-2300/​L1Trigger/​DTTraco/​src/​DTTracoChip.cc	Source navigation Diff markup Identifier search General search
	Version: CMSSW_14_2_X_2024-10-02-2300 CMSSW_14_2_X_2024-10-01-2300 CMSSW_14_2_X_2024-09-30-2300 CMSSW_14_2_X_2024-09-29-2300 CMSSW_14_2_X_2024-09-27-2300 CMSSW_14_2_X_2024-09-26-2300 Revert   Change

File indexing completed on 2024-09-07 04:36:45

 //-------------------------------------------------------------
 //
 //   Class: DTTracoChip
 //
 //   Description: Implementation of TRACO
 //                trigger algorithm
 //
 //
 //   Author List:
 //   SV 4/II/2003
 //   Modifications:
 //   22/VI/04 SV : last trigger code update
 //   16/I/07  SV : new DTConfig update
 //   3/IV/07  SV : setTracoAcceptance moved from card to chip
 //   30/10/09 SV : lut parameters from DB are used in code
 //   110208 SV   : TRACO hardware bug included
 //------------------------------------------------------------
 
 //-----------------------
 // This Class's Header --
 //-----------------------
 #include "L1Trigger/DTTraco/interface/DTTracoChip.h"
 
 //-------------------------------
 // Collaborating Class Headers --
 //-------------------------------
 #include "L1Trigger/DTBti/interface/DTBtiTrigData.h"
 #include "L1Trigger/DTTraco/interface/DTTracoCand.h"
 #include "L1Trigger/DTTraco/interface/DTTracoCard.h"
 #include "L1Trigger/DTTraco/interface/DTTracoTrig.h"
 #include "L1Trigger/DTTraco/interface/DTTracoTrigData.h"
 #include "L1Trigger/DTTriggerServerTheta/interface/DTTSTheta.h"
 #include "L1TriggerConfig/DTTPGConfig/interface/BitArray.h"
 
 //---------------
 // C++ Headers --
 //---------------
 #include <algorithm>
 #include <cmath>
 #include <iostream>
 #include <string>
 
 //----------------
 // Constructors --
 //----------------
 
 DTTracoChip::DTTracoChip(DTTracoCard *card, int n, DTConfigTraco *conf) : _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 = nullptr;
   } else
   // this is always the case with new DTConfig SV 15/I/2007
   // if( config()->trigSetupGeom()==0 ){
   {
     _luts = nullptr;
     _lutsCCB = nullptr;
   }
   /*
     //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)
     : _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];
   }
 }
 
 //--------------
 // Destructor --
 //--------------
 DTTracoChip::~DTTracoChip() {
   clear();
   /*
   if(config()->trigSetupGeom()==1){
     _luts->reset();
     delete _luts;
   }
 
   if(config()->trigSetupGeom()==2)
     delete _lutsCCB;
   */
 
   if (_card->lutFromDBFlag() == 1)
     delete _lutsCCB;
 }
 
 //--------------
 // Operations --
 //--------------
 
 DTTracoChip &DTTracoChip::operator=(const DTTracoChip &traco) {
   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;
 }
 
 void DTTracoChip::clear() {
   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();
 }
 
 void DTTracoChip::run() {
   // 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].empty() && _outerCand[is - DTConfigTraco::NSTEPF].empty())
       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 == nullptr && outer == nullptr)
         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].empty() && 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
 }
 
 void DTTracoChip::raiseOverlap(int step) {
   _flag[step - DTConfigTraco::NSTEPF].set(1);      // overlap flag raised
   _flag[step - DTConfigTraco::NSTEPF - 1].set(2);  // mark II rej.
 }
 
 void DTTracoChip::setFlag(int step, int ext) {
   if (ext == 0) {
     // this is the original: flags from card
     DTTracoChip *prevTraco = _card->getTRACO(_id.traco() - 1);
     if (prevTraco != nullptr) {
       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 != nullptr) {
       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
 }
 
 DTTracoCand *DTTracoChip::bestCand(int itk, std::vector<DTTracoCand> &tclist) {
   // Return if no candidates
   if (tclist.empty())
     return nullptr;
 
   // 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 = nullptr;
   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 == nullptr)
         bestltrig = &(*p);
     }
   }
   return bestltrig;
 }
 
 void DTTracoChip::DoAdjBtiLts(DTTracoCand *candidate, std::vector<DTTracoCand> &tclist) {
   // 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::AdjBtiLTSuppressed(DTTracoCand *candidate) {
   // 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;
 }
 
 DTTracoTrig *DTTracoChip::setPV(int itk, DTTracoCand *inner, DTTracoCand *outer) {
   // 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 = nullptr;
   if (inner != nullptr && outer != nullptr) {
     //    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 == nullptr && outer != nullptr) {
     candidate = outer;
   } else if (inner != nullptr && outer == nullptr) {
     candidate = inner;
   } else {
     return nullptr;  // 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;
 }
 
 int DTTracoChip::storeCorr(DTTracoTrig *tctrig, DTTracoCand *inner, DTTracoCand *outer, int tkn) {
   // 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) {
   // 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 = nullptr;
   if (inner != nullptr && outer != nullptr) {
     //    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 == nullptr && outer != nullptr) {
     candidate = outer;
   } else if (inner != nullptr && outer == nullptr) {
     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;
 }
 
 void DTTracoChip::add_btiT(int step, int pos, const DTBtiTrigData *btitrig) {
   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) {
   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::nTrig(int step) const {
   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();
 }
 
 DTTracoTrig *DTTracoChip::trigger(int step, unsigned n) const {
   if (step < DTConfigTraco::NSTEPF || step > DTConfigTraco::NSTEPL) {
     std::cout << "DTTracoChip::trigger: step out of range: " << step;
     std::cout << " empty pointer returned!" << std::endl;
     return nullptr;
   }
   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 nullptr;
   }
   std::vector<DTTracoTrig *>::const_iterator p = _tracotrig[step - DTConfigTraco::NSTEPF].begin() + n - 1;
   return *p;
 }
 
 DTTracoTrigData DTTracoChip::triggerData(int step, unsigned n) const {
   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();
 }
 
 float DTTracoChip::psiRad(int sl) const {
   /*
     // 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;
 }
 
 int DTTracoChip::KRad() const {
   // K parameter of the radial angle of correlator center
   // float distp2 = (int)(2*_geom->cellH()*config()->ST()/_geom->cellPitch());
   // return -(int)(tan(psiRad(sl)/1000)*distp2); // sign is reverted!!!
   // return _krad;
 
   // SV V/03: for harware bug in traco....
   int KRad = 0;
   return KRad;
 }
 
 int DTTracoChip::useSecondTrack(int step) const {
   // 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::edgeBTI(int step, int io, int lr) const {
   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].empty()) {
       // 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].empty()) {
       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;
 }
 
 void DTTracoChip::calculateAngles(DTTracoTrig *tct) {
   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
 }
 
 int DTTracoChip::insideAngWindow(DTTracoTrig *tctrig) const {
   // 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;
 }
 
 void DTTracoChip::setTracoAcceptances() {
   // 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&&sector()==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
 }

This page was automatically generated by the 2.2.1 LXR engine. The LXR team