/data/refman/pasoursint/CMSSW_5_3_9_patch3/src/L1Trigger/DTSectorCollector/src/DTSC.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: DTSC.cpp
//
//   Description: Implementation of DTSectColl trigger algorithm
//
//
//   Author List:
//   S. Marcellini
//   Modifications: 
//   11/11/06 C. Battilana : theta cand added
//   12/12/06 C. Battilana : _stat added
//   09/01/07 C. Battilana : updated to local conf
//
//
//--------------------------------------------------


//-----------------------
// This Class's Header --
//-----------------------
#include "L1Trigger/DTSectorCollector/interface/DTSC.h"

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "L1TriggerConfig/DTTPGConfig/interface/DTConfigSectColl.h"
#include "L1Trigger/DTSectorCollector/interface/DTSectColl.h"
#include "L1Trigger/DTSectorCollector/interface/DTSectCollPhCand.h"
#include "L1Trigger/DTSectorCollector/interface/DTSectCollThCand.h"

//---------------
// C++ Headers --
//---------------
#include <iostream>
#include <algorithm>

//----------------
// Constructors --
//----------------

DTSC::DTSC(int istat) : _ignoreSecondTrack(0) ,_stat(istat){

  // reserve the appropriate amount of space for vectors
  // test _incand[0].reserve(DTConfigSectColl::NTSMSC);
  // test_incand[1].reserve(DTConfigSectColl::NTSMSC);
  // test _outcand.reserve(2);
 
}


//--------------
// Destructor --
//--------------
DTSC::~DTSC() {

  clear();

}


//--------------
// Operations --
//--------------

void
DTSC::clear() {

  _ignoreSecondTrack=0;

  for(int itk=0;itk<=1;itk++){
 
    _incand_ph[itk].clear();
  }

  _outcand_ph.clear();
  _cand_th.clear();

}


// 
void
DTSC::run() {

  if(config()->debug()){
    std::cout << "DTSC::run: Processing DTSectColl: ";
    std::cout << nFirstTPh() << " first & " << nSecondTPh() << " second Phi tracks ";
    std::cout << " - "<< nCandTh() << " Theta tracks" << std::endl;
  
  }

  if(nFirstTPh()<1)return; // skip if no first tracks
  //
  // SORT 1
  //

  // debugging
    if(config()->debug()){
     std::cout << "Vector of first Phi tracks in DTSectColl: " << std::endl;
    std::vector<DTSectCollPhCand*>::const_iterator p;
    for(p=_incand_ph[0].begin(); p!=_incand_ph[0].end(); p++) {
           (*p)->print();
    }
   }
  // end debugging
 
  DTSectCollPhCand* first=DTSectCollsort1();
  if(config()->debug()){
    std::cout << "SC: DTSC::run: first Phi track is = " << first << std::endl;
  }
  if(first!=0) {
    _outcand_ph.push_back(first); 

  }
  if(nSecondTPh()<1)return; // skip if no second tracks

  //
  // SORT 2
  //

  // debugging
  if(config()->debug()){
    std::vector<DTSectCollPhCand*>::const_iterator p;
    std::cout << "Vector of second Phi tracks in DTSectColl: " << std::endl;
    for(p=_incand_ph[1].begin(); p!=_incand_ph[1].end(); p++) {
       (*p)->print();
    }
  }
  // end debugging

  DTSectCollPhCand* second=DTSectCollsort2();
  if(second!=0) {
    _outcand_ph.push_back(second); 
  }
  
}


DTSectCollPhCand*
DTSC::DTSectCollsort1() {

  // Do a sort 1
  DTSectCollPhCand* best=0;
  DTSectCollPhCand* carry=0;
  std::vector<DTSectCollPhCand*>::iterator p;
  for(p=_incand_ph[0].begin(); p!=_incand_ph[0].end(); p++) {
    DTSectCollPhCand* curr=(*p);

    curr->setBitsSectColl();    // SM sector collector set bits in dataword to make SC sorting
    
    // NO Carry in Sector Collector sorting in default 
    if(config()->SCGetCarryFlag(_stat)) {  // get carry

      if(best==0){
        best=curr;
      } 
      else if((*curr)<(*best)){
        carry=best;
        best=curr;
      } 
      else if(carry==0){
        carry=curr;
      } 
      else if((*curr)<(*carry)){
        carry=curr;
      } 

    }
    else if(config()->SCGetCarryFlag(_stat)==0){ // no carry (default)
      if(best==0){
        best=curr;
      } 
      else if((*curr)<(*best)){
        
        best=curr;
      } 
      
    }
    
    if(carry!=0 && config()->SCGetCarryFlag(_stat)) { // reassign carry to sort 2 candidates
      carry->setSecondTrack(); // change value of 1st/2nd track bit
      _incand_ph[1].push_back(carry); // add to list of 2nd track
 
    }
  } 
  
  return best;

}


DTSectCollPhCand*
DTSC::DTSectCollsort2() {

  // Check if there are second tracks

  if(nTracksPh()<1){
    std::cout << "DTSC::DTSectCollsort2: called with no first Phi track.";
    std::cout << " empty pointer returned!" << std::endl;
    return 0;
  }
  // If a first track at the following BX is present, ignore second tracks of any kind
  if(_ignoreSecondTrack){

    for(std::vector<DTSectCollPhCand*>::iterator p=_incand_ph[1].begin(); p!=_incand_ph[1].end(); p++) {

    }
    return 0;
  }

  // If no first tracks at the following BX, do a sort 2
  //  DTSectCollCand* best=getTrack(1);  ! not needed as lons as there is no comparison with best in sort 2
  DTSectCollPhCand* second=0;
  std::vector<DTSectCollPhCand*>::iterator p;
  for(p=_incand_ph[1].begin(); p!=_incand_ph[1].end(); p++) {
    DTSectCollPhCand* curr=(*p);
    curr->setBitsSectColl();    // SM sector collector set bits in dataword to make SC sorting
    
    if(second==0){
      second=curr;
    } 
    else if((*curr)<(*second)){
      second=curr;
    } 
    
  }

  return second;

}


void
DTSC::addPhCand(DTSectCollPhCand* cand) {

  _incand_ph[(1-cand->isFirst())].push_back(cand); 

}

void
DTSC::addThCand(DTSectCollThCand* cand) {

  _cand_th.push_back(cand); 

}


unsigned
DTSC::nCandPh(int ifs) const {

  if(ifs<1||ifs>2){
    std::cout << "DTSC::nCandPh: wrong track number: " << ifs;
    std::cout << " 0 returned!" << std::endl;
    return 0;
  }
  return _incand_ph[ifs-1].size();

}

unsigned
DTSC::nCandTh() const {

  return _cand_th.size();

}


DTSectCollPhCand*
DTSC::getDTSectCollPhCand(int ifs, unsigned n) const {

  if(ifs<1||ifs>2){
    std::cout << "DTSC::getDTSectCollPhCand: wrong track number: " << ifs;
    std::cout << " empty pointer returned!" << std::endl;
    return 0;
  }
  if(n<1 || n>nCandPh(ifs)) {
    std::cout << "DTSC::getDTSectCollPhCand: requested trigger not present: " << n;
    std::cout << " empty pointer returned!" << std::endl;
    return 0;
  }

  std::vector<DTSectCollPhCand*>::const_iterator p = _incand_ph[ifs-1].begin()+n-1;
  return (*p);

}

DTSectCollThCand*
DTSC::getDTSectCollThCand(unsigned n) const {

  if(n<1 || n>nCandTh()) {
    std::cout << "DTSC::getDTSectCollThCand: requested trigger not present: " << n;
    std::cout << " empty pointer returned!" << std::endl;
    return 0;
  }

  std::vector<DTSectCollThCand*>::const_iterator p = _cand_th.begin()+n-1;
  return (*p);

}

void
DTSC::addDTSectCollPhCand(DTSectCollPhCand* cand) {

  int ifs = (cand->isFirst()) ? 0 : 1;
 
  _incand_ph[ifs].push_back(cand); 

}

DTSectCollPhCand*
DTSC::getTrackPh(int n) const {

  if(n<1 || n>nTracksPh()) {
    std::cout << "DTSC::getTrackPh: requested track not present: " << n;
    std::cout << " empty pointer returned!" << std::endl;
    return 0;
  }

  std::vector<DTSectCollPhCand*>::const_iterator p = _outcand_ph.begin()+n-1;

  return (*p);

}


DTSectCollThCand*
DTSC::getTrackTh(int n) const {

  if(n<1 || n>nTracksTh()) {
    std::cout << "DTSC::getTrackTh: requested track not present: " << n;
    std::cout << " empty pointer returned!" << std::endl;
    return 0;
  }

  std::vector<DTSectCollThCand*>::const_iterator p = _cand_th.begin()+n-1;

  return (*p);

}