#include <L1MuBMMuonSorter.h>

Inheritance diagram for L1MuBMMuonSorter:

Public Member Functions
	L1MuBMMuonSorter (const L1MuBMTrackFinder &)
	constructor More...

int	numberOfTracks () const
	return number of found muon candidates after sorter More...

void	print () const
	print results after sorting More...

virtual void	reset ()
	reset Muon Sorter More...

virtual void	run ()
	run Muon Sorter More...

const L1MuBMTrack *	track (int id) const
	return pointer to a muon candidate More...

const std::vector< const L1MuBMTrack * > &	tracks () const
	return vector of muon candidates More...

virtual	~L1MuBMMuonSorter ()
	destructor More...

Public Member Functions inherited from L1AbstractProcessor
virtual void	run (const edm::EventSetup &c)

virtual void	run (const edm::EventSetup &c)

virtual	~L1AbstractProcessor ()
	destructor More...

virtual	~L1AbstractProcessor ()
	destructor More...

Private Member Functions
void	runCOL (std::vector< L1MuBMTrack * > &) const
	run the Cancel Out Logic of the muon sorter More...

Static Private Member Functions
static int	neighbour (const L1MuBMSecProcId &spid1, const L1MuBMSecProcId &spid2)
	find out if two Sector Processors are neighbours More...

Private Attributes
const L1MuBMTrackFinder &	m_tf

std::vector< const L1MuBMTrack * >	m_TrackCands

Detailed Description

BM Muon Sorter:

The BM Muon Sorter receives 2 muon candidates from each of the 12 Wedge Sorters and sorts out the 4 best (highest pt, highest quality) muons

N. Neumeister CERN EP

Definition at line 45 of file L1MuBMMuonSorter.h.

Constructor & Destructor Documentation

L1MuBMMuonSorter::L1MuBMMuonSorter ( const L1MuBMTrackFinder & tf )

constructor

Definition at line 48 of file L1MuBMMuonSorter.cc.

References m_TrackCands.

                                                               :
       m_tf(tf), m_TrackCands() {
 
   m_TrackCands.reserve(4);
 
 }

L1MuBMMuonSorter::~L1MuBMMuonSorter ( )

virtual

destructor

Definition at line 60 of file L1MuBMMuonSorter.cc.

                                     {
 
 }

Member Function Documentation

int L1MuBMMuonSorter::neighbour	(	const L1MuBMSecProcId &	spid1,
		const L1MuBMSecProcId &	spid2
	)

staticprivate

find out if two Sector Processors are neighbours

Definition at line 298 of file L1MuBMMuonSorter.cc.

References funct::abs(), L1MuBMSecProcId::sector(), ecaldqm::topology(), and L1MuBMSecProcId::wheel().

Referenced by runCOL(), and tracks().

                                                               {
 
   // definition of valid topologies:
 
   //              E T A
   //        -----------------
   //   +
   //        ---               ---
   //   |   | 2 |             | 2 |
   // P |   |___|             |___|
   //   |    ---    ---    ---        ---
   // H |   | 1 |  | 1 |  | 1 |      | 1 |
   //   |   |___|  |___|  |___|      |___|
   // I |           ---                   ---
   //   |          | 2 |                 | 2 |
   //              |___|                 |___|
   //   -
   // result: 1      2        3          4         5    6  otherwise : -1
 
   int topology = -1;
 
   int sector1 = spid1.sector();
   int wheel1  = spid1.wheel();
 
   int sector2 = spid2.sector();
   int wheel2  = spid2.wheel();
 
   int sectordiff = (sector2 - sector1)%12;
   if ( sectordiff >= 6 ) sectordiff -= 12;
   if ( sectordiff < -6 ) sectordiff += 12;
 
   if ( abs(sectordiff) == 1 ) {
 
     if ( wheel1 == wheel2 ) topology = (sectordiff > 0) ? 1 : 2;
     if ( wheel1 == +1 && wheel2 == -1 )  topology = (sectordiff > 0) ? 5 : 6;
     if ( ( wheel1 == -2 && wheel2 == -3 ) ||
          ( wheel1 == -1 && wheel2 == -2 ) ||
          ( wheel1 == +1 && wheel2 == +2 ) ||
          ( wheel1 == +2 && wheel2 == +3 ) ) topology = (sectordiff > 0) ? 3 : 4;
 
   }
 
   return topology;
 
 }

int L1MuBMMuonSorter::numberOfTracks ( ) const

inline

return number of found muon candidates after sorter

Definition at line 65 of file L1MuBMMuonSorter.h.

References m_TrackCands.

Referenced by print().

65 { return m_TrackCands.size(); }

L1MuBMMuonSorter::m_TrackCands

std::vector< const L1MuBMTrack * > m_TrackCands

Definition: L1MuBMMuonSorter.h:84

void L1MuBMMuonSorter::print ( void ) const

print results after sorting

Definition at line 145 of file L1MuBMMuonSorter.cc.

References gather_cfg::cout, m_TrackCands, and numberOfTracks().

                                    {
 
   cout << endl;
   cout << "Muon candidates found by the barrel MTTF : "
        << numberOfTracks() << endl;
   vector<const L1MuBMTrack*>::const_iterator iter = m_TrackCands.begin();
   while ( iter != m_TrackCands.end() ) {
     if ( *iter ) cout << *(*iter) << endl;
     iter++;
   }
   cout << endl;
 
 }

void L1MuBMMuonSorter::reset ( void )

virtual

reset Muon Sorter

Implements L1AbstractProcessor.

Definition at line 131 of file L1MuBMMuonSorter.cc.

References m_TrackCands.

Referenced by L1MuBMTrackFinder::reset().

                              {
 
   m_TrackCands.clear();
   vector<const L1MuBMTrack*>::iterator iter;
   for ( iter = m_TrackCands.begin(); iter != m_TrackCands.end(); iter++ ) {
     *iter = 0;
   }
 
 }

void L1MuBMMuonSorter::run ( void )

virtual

run Muon Sorter

Reimplemented from L1AbstractProcessor.

Definition at line 72 of file L1MuBMMuonSorter.cc.

References gather_cfg::cout, L1MuBMTFConfig::Debug(), m_tf, m_TrackCands, runCOL(), L1MuBMWedgeSorter::tracks(), and L1MuBMTrackFinder::ws().

                            {
 
   // get track candidates from Wedge Sorters
   vector<L1MuBMTrack*> mycands;
   mycands.reserve(24);
 
   for ( int wedge = 0; wedge < 12; wedge++ ) {
     vector<const L1MuBMTrack*> wscand = m_tf.ws(wedge)->tracks();
     vector<const L1MuBMTrack*>::iterator iter = wscand.begin();
     while ( iter != wscand.end() ) {
       if ( *iter && !(*iter)->empty() )
         mycands.push_back(const_cast<L1MuBMTrack*>(*iter) );
       iter++;
     }
   }
 
   // print input data
   if ( L1MuBMTFConfig::Debug(4) ) {
     cout << "BM Muon Sorter input: "
          << mycands.size() << endl;
     vector<L1MuBMTrack*>::const_iterator iter;
     for ( iter = mycands.begin(); iter != mycands.end(); iter++ ) {
       if (*iter ) (*iter)->print();
     }
   }
 
   // run Cancel Out Logic
   runCOL(mycands);
 
   // remove disabled candidates
   vector<L1MuBMTrack*>::iterator it = mycands.begin();
   while ( it != mycands.end() ) {
     if ( *it && (*it)->empty() ) {
       mycands.erase(it);
       it = mycands.begin(); continue;
     }
     it++;
   }
 
   // sort pt and quality
   stable_sort( mycands.begin(), mycands.end(), L1MuBMTrack::Rank() );
 
   // copy the best 4 candidates
   int number_of_tracks = 0;
   vector<L1MuBMTrack*>::const_iterator iter1 = mycands.begin();
   while ( iter1 != mycands.end() ) {
     if ( *iter1 && number_of_tracks < 4 ) {
       m_TrackCands.push_back(*iter1);
       number_of_tracks++;
     }
     iter1++;
   }
 
 }

void L1MuBMMuonSorter::runCOL ( std::vector< L1MuBMTrack * > & cands ) const

private

run the Cancel Out Logic of the muon sorter

Definition at line 163 of file L1MuBMMuonSorter.cc.

References funct::abs(), gather_cfg::cout, L1MuBMTFConfig::Debug(), neighbour(), hitfit::phidiff(), hiDetachedQuadStep_cff::pt1, hiDetachedQuadStep_cff::pt2, ecaldqm::topology(), and w2.

Referenced by run(), and tracks().

                                                                {
 
   // compare candidates which were found in nearby sectors and wheels
   // if 2 candidates have at least one track segment in common
   // disable the one with lower quality
   // compare addresses from stations 2, 3 and 4
 
   typedef vector<L1MuBMTrack*>::iterator TI;
   for ( TI iter1 = cands.begin(); iter1 != cands.end(); iter1++ ) {
     if ( *iter1 == 0 ) continue;
     if ( (*iter1)->empty() ) continue;
     L1MuBMSecProcId sp1 = (*iter1)->spid();
     int qual1 = (*iter1)->quality();
     for ( TI iter2 = cands.begin(); iter2 != cands.end(); iter2++ ) {
       if ( *iter2 == 0 ) continue;
       if ( *iter1 == *iter2 ) continue;
       if ( (*iter2)->empty() ) continue;
       L1MuBMSecProcId sp2 = (*iter2)->spid();
       int qual2 = (*iter2)->quality();
       if (sp1 == sp2 ) continue;
       int topology = neighbour(sp1,sp2);
       if ( topology == -1 ) continue;
       int countTS = 0;
       for ( int stat = 2; stat <= 4; stat++ ) {
         int adr1 = (*iter1)->address(stat);
         int adr2 = (*iter2)->address(stat);
         if ( adr1 == 15 || adr2 == 15 ) continue;
         switch ( topology ) {
           case 1 : {
                     if ( adr1 > 7 ) continue;
                     if ( adr2 > 3 && adr2 < 8 ) continue;
                     int adr_shift = ( adr2 > 7 ) ? -8 : 4;
                     if ( adr1 == adr2+adr_shift ) countTS++;
                     break;
                    }
           case 2 : {
                     if ( adr2 > 7 ) continue;
                     if ( adr1 > 3 && adr1 < 8 ) continue;
                     int adr_shift = ( adr2 > 3 ) ? -4 : 8;
                     if ( adr1 == adr2+adr_shift ) countTS++;
                     break;
                    }
           case 3 : {
                     if ( ( adr1 == 6 && adr2 == 0 ) ||
                          ( adr1 == 7 && adr2 == 1 ) ||
                          ( adr1 == 2 && adr2 == 8 ) ||
                          ( adr1 == 3 && adr2 == 9 ) ) countTS++;
                     break;
                    }
           case 4 : {
                     if ( ( adr1 == 2  && adr2 == 4 ) ||
                          ( adr1 == 3  && adr2 == 5 ) ||
                          ( adr1 == 10 && adr2 == 0 ) ||
                          ( adr1 == 11 && adr2 == 1 ) ) countTS++;
                     break;
                    }
           case 5 : {
                     if ( ( adr1 == 0 && adr2 == 8 ) ||
                          ( adr1 == 1 && adr2 == 9 ) ||
                          ( adr1 == 4 && adr2 == 0 ) ||
                          ( adr1 == 5 && adr2 == 1 ) ) countTS++;
                     break;
                    }
           case 6 : {
                     if ( ( adr1 == 0 && adr2 == 4 ) ||
                          ( adr1 == 1 && adr2 == 5 ) ||
                          ( adr1 == 8 && adr2 == 0 ) ||
                          ( adr1 == 9 && adr2 == 1 ) ) countTS++;
                     break;
                    }
           default : break;
         }
       }
       if (  countTS > 0 ) {
         if ( qual1 < qual2 ) {
           if ( L1MuBMTFConfig::Debug(5) ) {
             cout << "Muon Sorter cancel : "; (*iter1)->print();
           }
           (*iter1)->disable();
           break;
         }
         else {
           if ( L1MuBMTFConfig::Debug(5) ) {
             cout << "Muon Sorter cancel : "; (*iter2)->print();
           }
          (*iter2)->disable();
         }
       }
     }
   }
 
 
   // if two candidates have exactly the same phi and eta values
   // remove the one with lower rank
 
   for ( TI iter1 = cands.begin(); iter1 != cands.end(); iter1++ ) {
     if ( *iter1 == 0 ) continue;
     if ( (*iter1)->empty() ) continue;
     int phi1 = (*iter1)->phi();
     int pt1 = (*iter1)->pt();
     int qual1 = (*iter1)->quality();
     for ( TI iter2 = cands.begin(); iter2 != cands.end(); iter2++ ) {
       if ( *iter2 == 0 ) continue;
       if ( *iter1 == *iter2 ) continue;
       if ( (*iter2)->empty() ) continue;
       int phi2 = (*iter2)->phi();
       int pt2 = (*iter2)->pt();
       int qual2 = (*iter2)->quality();
       int w1 = (*iter1)->getStartTSphi().wheel();
       int w2 = (*iter2)->getStartTSphi().wheel();
       int phidiff = (phi2 - phi1)%144;
       if ( phidiff >= 72 ) phidiff -= 144;
       if ( phidiff < -72 ) phidiff += 144;
       if ( abs(phidiff) < 2 && (w1 == w2) ) {
         int rank1 = 10 * pt1 + qual1;
         int rank2 = 10 * pt2 + qual2;
         if ( L1MuBMTFConfig::Debug(5) ) {
           cout << "========================================" << endl;
           cout << " special cancellation : " << endl;
           (*iter1)->print(); if ( rank1 <  rank2 ) cout << "cancelled" << endl;
           (*iter2)->print(); if ( rank1 >= rank2 ) cout << "cancelled" << endl;
           cout << "========================================" << endl;
         }
         if ( rank1 >= rank2 ) (*iter2)->disable();
         if ( rank1 <  rank2 ) { (*iter1)->disable(); break; }
       }
     }
   }
 
 }