#include <L1MuBMEtaProcessor.h>

Public Member Functions
int	eta (int id) const
	return eta values, index [0,11] More...

bool	fineBit (int id) const
	return fine bit, index [0,11] More...

int	id () const
	return Eta Processor identifier (0-11) More...

	L1MuBMEtaProcessor (const L1MuBMTrackFinder &, int id, edm::ConsumesCollector &&iC)
	constructor More...

void	print () const
	print muon candidates found by the Eta Processor More...

virtual void	reset ()
	reset the Eta Processor More...

virtual void	run (int bx, const edm::Event &e, const edm::EventSetup &c)
	run the Eta Processor More...

const L1MuBMTrackFinder &	tf () const
	return reference to barrel MTTF More...

virtual	~L1MuBMEtaProcessor ()
	destructor More...

Private Member Functions
void	assign ()
	assign eta and etaFineBit More...

void	receiveAddresses ()
	receive addresses (from 6 Sector Processors) More...

void	receiveData (int bx, const edm::Event &e, const edm::EventSetup &c)
	receive data (eta trigger primitives) More...

void	runEtaMatchingUnit (const edm::EventSetup &c)
	run Eta Matching Unit (EMU) More...

void	runEtaTrackFinder (const edm::EventSetup &c)
	run Eta Track Finder (ETF) More...

Static Private Member Functions
static int	quality (int id, int stat)
	get quality code; id [0,26], stat [1,3] More...

Private Attributes
edm::ESHandle< L1TMuonBarrelParams >	bmtfParamsHandle

int	m_address [12]

edm::EDGetTokenT< L1MuDTChambThContainer >	m_DTDigiToken

int	m_epid

int	m_eta [12]

bool	m_fine [12]

std::vector< int >	m_foundPattern

int	m_mask

int	m_pattern [12]

const L1MuBMTrackFinder &	m_tf

L1MuBMTrack *	m_TrackCand [12]

L1MuBMTrack *	m_TracKCand [12]

std::vector< const L1MuBMTrackSegEta * >	m_tseta

L1MuDTTFMasks	msks

L1MuBMTEtaPatternLut	theEtaPatternLUT

L1MuBMTQualPatternLut	theQualPatternLUT

Detailed Description

Eta Processor:

An Eta Processor consists of :

one receiver unit,
one Eta Track Finder (ETF) and
one Eta Matching Unit (EMU)

N. Neumeister CERN EP J. Troconiz UAM Madrid

Definition at line 62 of file L1MuBMEtaProcessor.h.

Constructor & Destructor Documentation

L1MuBMEtaProcessor::L1MuBMEtaProcessor	(	const L1MuBMTrackFinder &	tf,
		int	id,
		edm::ConsumesCollector &&	iC
	)

constructor

Definition at line 64 of file L1MuBMEtaProcessor.cc.

References m_tseta.

                                                                                                      :
       m_tf(tf), m_epid(id), m_foundPattern(0), m_tseta(15), m_DTDigiToken(iC.consumes<L1MuDTChambThContainer>(L1MuBMTFConfig::getBMThetaDigiInputTag()))
  {
 
   m_tseta.reserve(15);
 
 }

L1MuBMEtaProcessor::~L1MuBMEtaProcessor ( )

virtual

destructor

Definition at line 77 of file L1MuBMEtaProcessor.cc.

77 {}

Member Function Documentation

void L1MuBMEtaProcessor::assign ( )

private

assign eta and etaFineBit

Definition at line 434 of file L1MuBMEtaProcessor.cc.

References L1MuBMTEtaPatternLut::getPattern(), mps_fire::i, m_eta, m_fine, m_pattern, m_TrackCand, m_TracKCand, m_tseta, AlCaHLTBitMon_ParallelJobs::p, L1MuDTEtaPattern::position(), L1MuBMTrack::setEta(), L1MuBMTrack::setFineEtaBit(), L1MuBMTrack::setTSeta(), theEtaPatternLUT, and L1MuDTEtaPattern::wheel().

Referenced by fineBit(), and run().

                                 {
 
   for ( int i = 0; i < 12; i++ ) {
     if ( m_TrackCand[i] ) {
       if ( m_eta[i] != 99 ) {
         m_TrackCand[i]->setEta(m_eta[i]);
         m_TracKCand[i]->setEta(m_eta[i]);
       }
       else {
     //        if ( i/2 != 2 ) cerr << "L1MuBMEtaProcessor: assign invalid eta" << " " << m_address[i] << endl;
       }
       if ( m_fine[i] ) {
         m_TrackCand[i]->setFineEtaBit();
         m_TracKCand[i]->setFineEtaBit();
         // find all contributing track segments
         const L1MuDTEtaPattern p = theEtaPatternLUT.getPattern(m_pattern[i]);
         vector<const L1MuBMTrackSegEta*> TSeta;
         const L1MuBMTrackSegEta* ts = nullptr;
         for ( int stat = 0; stat < 3; stat++ ) {
           int wh = p.wheel(stat+1);
           int pos = p.position(stat+1);
           if ( pos == 0 ) continue;
           ts = m_tseta[wh+2 + stat*5];
           TSeta.push_back(ts);
         }
         m_TrackCand[i]->setTSeta(TSeta);
         m_TracKCand[i]->setTSeta(TSeta);
       }
     }
   }
 
 }

int L1MuBMEtaProcessor::eta ( int id ) const

inline

return eta values, index [0,11]

Definition at line 88 of file L1MuBMEtaProcessor.h.

References id(), and m_eta.

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

88 { return m_eta[id]; }

L1MuBMEtaProcessor::m_eta

int m_eta[12]

Definition: L1MuBMEtaProcessor.h:120

L1MuBMEtaProcessor::id

int id() const

return Eta Processor identifier (0-11)

Definition: L1MuBMEtaProcessor.h:73

bool L1MuBMEtaProcessor::fineBit ( int id ) const

inline

return fine bit, index [0,11]

Definition at line 91 of file L1MuBMEtaProcessor.h.

References assign(), id(), m_fine, quality(), receiveAddresses(), receiveData(), runEtaMatchingUnit(), and runEtaTrackFinder().

91 { return m_fine[id]; }

L1MuBMEtaProcessor::m_fine

bool m_fine[12]

Definition: L1MuBMEtaProcessor.h:121

L1MuBMEtaProcessor::id

int id() const

return Eta Processor identifier (0-11)

Definition: L1MuBMEtaProcessor.h:73

int L1MuBMEtaProcessor::id ( void ) const

inline

return Eta Processor identifier (0-11)

Definition at line 73 of file L1MuBMEtaProcessor.h.

References EnergyCorrector::c, MillePedeFileConverter_cfg::e, m_epid, print(), reset(), and run().

Referenced by eta(), fineBit(), and quality().

73 { return m_epid; }

L1MuBMEtaProcessor::m_epid

int m_epid

Definition: L1MuBMEtaProcessor.h:116

void L1MuBMEtaProcessor::print ( void ) const

print muon candidates found by the Eta Processor

Definition at line 137 of file L1MuBMEtaProcessor.cc.

References gather_cfg::cout, L1MuDTEtaPattern::eta(), L1MuBMTQualPatternLut::getCoarseEta(), L1MuBMTEtaPatternLut::getPattern(), mps_fire::i, L1MuDTEtaPattern::id(), m_address, m_epid, m_eta, m_fine, m_foundPattern, m_pattern, m_tseta, AlCaHLTBitMon_ParallelJobs::p, position, L1MuDTEtaPattern::quality(), quality(), relativeConstraints::station, theEtaPatternLUT, and theQualPatternLUT.

Referenced by id().

                                      {
 
   bool empty1 = true;
   for ( int i = 0; i < 15; i++ ) {
     empty1 &= ( m_tseta[i] == nullptr || m_tseta[i]->empty() );
   }
 
   bool empty2 = true;
   for ( int i = 0; i < 12; i++ ) {
     empty2 &= ( m_address[i] == 0 );
   }
 
   if ( !empty1 || !empty2 ) {
     cout << "Eta processor " << m_epid << " : " << endl;
 
    // print local pattern
    if ( !empty1 ) {
      cout << "Local pattern : " << endl;
      for ( int i = 0; i < 15; i++ ) {
        if ( (i+5)%5 == 0 ) cout << "station " << m_tseta[i]->station() << " : ";
        bitset<7> pos(m_tseta[i]->position());
        bitset<7> qua(m_tseta[i]->quality());
        for ( int j = 6; j >= 0; j-- ) {
          cout << pos[j]+qua[j];
        }
        cout << " ";
        if ( (i+1)%5 == 0 ) cout << endl;
      }
      cout << "Found patterns :" << endl;
      vector<int>::const_iterator iter;
      for ( iter = m_foundPattern.begin(); iter != m_foundPattern.end(); iter++ ) {
         const L1MuDTEtaPattern p = theEtaPatternLUT.getPattern(*iter);
         int qualitycode = p.quality();
         cout << "ID = " << setw(4) << p.id() << "  "
              << "eta = " << setw(3) << p.eta() << "  "
              << "quality = " << setw(2) << qualitycode << " ("
              << quality(qualitycode,1) << " "
              << quality(qualitycode,2) << " "
              << quality(qualitycode,3) << ")";
         for ( int i = 0; i < 12; i++ ) {
           if ( m_pattern[i] ==  p.id() ) cout << " <--";
         }
         cout << endl;
       }
     }
 
     cout << "Received addresses : " << endl;
     for ( int i = 0; i < 12; i++ ) cout << setw(3) << m_address[i] << " ";
     cout << endl;
 
     if ( !empty1 ) {
       cout << "Matched patterns : " << endl;
       for ( int i = 0; i < 12; i++ ) {
         if ( m_fine[i] ) {
           const L1MuDTEtaPattern p = theEtaPatternLUT.getPattern(m_pattern[i]);
           int fineeta = p.eta();
           int coarseeta = theQualPatternLUT.getCoarseEta(i/2+1,m_address[i]);
           cout << "Index = " << setw(2) << i << ", "
                << "address = " << setw(2) << m_address[i] << " --> "
                << "pattern = " << setw(4) << m_pattern[i] << " "
                << "eta (coarse) = " << setw(3) << coarseeta << " "
                << "eta (fine) = " << setw(3) << fineeta << " "
                << "quality = " << setw(2) << p.quality() << endl;
         }
       }
     }
 
     cout << "Eta values and fine bits : " << endl;
     for ( int i = 0; i < 12; i++ ) cout << setw(3) << m_eta[i] << " ";
     cout << endl;
     for ( int i = 0; i < 12; i++ ) cout << setw(3) << m_fine[i] << " ";
     cout << endl;
   }
 
 }

int L1MuBMEtaProcessor::quality	(	int	id,
		int	stat
	)

staticprivate

get quality code; id [0,26], stat [1,3]

Definition at line 471 of file L1MuBMEtaProcessor.cc.

References id().

Referenced by fineBit(), print(), and runEtaTrackFinder().

                                                 {
 
   // quality codes as defined in CMS Note 2001/027
   // This QualityPatterns are used to have a defined Quality-Identifier for
   // all possible found tracks.
   // Therefore three integer numbers ( Station 1, 2, 3 from left to right )
   // can have numbers like 0, 1 or 2
   //    0 ... no hit is given
   //    1 ... a LTRG is given
   //    2 ... a HTRG is given
 
   const int qualcode[27][3] = { {0,0,0},{1,0,0},{0,1,0},{0,0,1},{2,0,0},
                                 {0,2,0},{0,0,2},{1,1,0},{1,0,1},{0,1,1},
                                 {2,1,0},{1,2,0},{2,0,1},{1,0,2},{0,2,1},
                                 {0,1,2},{2,2,0},{2,0,2},{0,2,2},{1,1,1},
                                 {2,1,1},{1,2,1},{1,1,2},{2,2,1},{2,1,2},
                                 {1,2,2},{2,2,2} };
 
   return qualcode[id][stat-1];
 
 }

void L1MuBMEtaProcessor::receiveAddresses ( )

private

receive addresses (from 6 Sector Processors)

Definition at line 282 of file L1MuBMEtaProcessor.cc.

References L1MuBMTrack::address(), L1MuBMTrack::empty(), mps_fire::i, m_address, m_epid, m_tf, m_TrackCand, m_TracKCand, L1MuBMTrackFinder::sp(), L1MuBMSectorProcessor::track(), L1MuBMSectorProcessor::tracK(), L1MuBMAddressArray::trackAddressCode(), and makeMuonMisalignmentScenario::wheel.

Referenced by fineBit(), and run().

                                           {
 
   // get track address code of all track segments
   // 6*2 times 5 bits; valid range [1-22]
 
   int sector = m_epid;
 
   int i = 0;
   for ( int wheel = -3; wheel <= 3; wheel++ ) {
     if ( wheel == 0 ) continue;
     L1MuBMSecProcId tmpspid(wheel,sector);
     for ( int number = 0; number < 2; number++ ) {
       const L1MuBMTrack* cand = m_tf.sp(tmpspid)->track(number);
       const L1MuBMTrack* canD = m_tf.sp(tmpspid)->tracK(number);
       if ( cand ) {
         m_address[i] = cand->address().trackAddressCode();
         if ( !cand->empty() ) {
           m_TrackCand[i] = const_cast<L1MuBMTrack*>(cand);
           m_TracKCand[i] = const_cast<L1MuBMTrack*>(canD);
         }
         i++;
       }
     }
   }
 
 }

void L1MuBMEtaProcessor::receiveData	(	int	bx,
		const edm::Event &	e,
		const edm::EventSetup &	c
	)

private

receive data (eta trigger primitives)

Definition at line 217 of file L1MuBMEtaProcessor.cc.

References bmtfParamsHandle, L1MuDTChambThContainer::chThetaSegm(), L1TMuonBarrelParams::LUTParams::eta_lut_, edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), L1MuDTTFMasks::get_etsoc_chdis_st1(), L1MuDTTFMasks::get_etsoc_chdis_st2(), L1MuDTTFMasks::get_etsoc_chdis_st3(), edm::Event::getByToken(), mps_fire::i, L1TMuonBarrelParams::l1mudttfmasks, L1TMuonBarrelParams::lutparams_, m_DTDigiToken, m_epid, L1MuBMTEtaPatternLut::m_lut, L1MuBMTQualPatternLut::m_lut, m_mask, m_tseta, msks, L1MuDTChambThDigi::position(), edm::ESHandle< T >::product(), L1TMuonBarrelParams::LUTParams::qp_lut_, L1MuDTChambThDigi::quality(), theEtaPatternLUT, theQualPatternLUT, and makeMuonMisalignmentScenario::wheel.

Referenced by fineBit(), and run().

                                                                                       {
 
   //c.get< L1MuDTTFMasksRcd >().get( msks );
   const L1TMuonBarrelParamsRcd& bmtfParamsRcd = c.get<L1TMuonBarrelParamsRcd>();
   bmtfParamsRcd.get(bmtfParamsHandle);
   const L1TMuonBarrelParams& bmtfParams = *bmtfParamsHandle.product();
   msks =  bmtfParams.l1mudttfmasks;
   theEtaPatternLUT.m_lut  = bmtfParams.lutparams_.eta_lut_; //l1mudttfetaplut;  // ETF look-up table
   theQualPatternLUT.m_lut = bmtfParams.lutparams_.qp_lut_; //l1mudttfqualplut; // EMU look-up tables
 
   edm::Handle<L1MuDTChambThContainer> dttrig;
 //  e.getByLabel(L1MuBMTFConfig::getBMThetaDigiInputTag(),dttrig);
   e.getByToken(m_DTDigiToken,dttrig);
 
   // const int bx_offset = dttrig->correctBX();
   int bx_offset=0;
   bx = bx + bx_offset;
 
   //
   // get 5*3 eta track segments
   //
   int sector = m_epid;
   for ( int stat = 1; stat <= 3; stat++ ) {
     for ( int wheel = -2; wheel <= 2; wheel++ ) {
       L1MuDTChambThDigi const* tseta = dttrig->chThetaSegm(wheel,stat,sector,bx);
       bitset<7> pos;
       bitset<7> qual;
 
       int lwheel = wheel+1;
       if ( wheel < 0 ) lwheel = wheel-1;
 
       bool masked = false;
       if ( stat == 1 ) masked = msks.get_etsoc_chdis_st1(lwheel, sector);
       if ( stat == 2 ) masked = msks.get_etsoc_chdis_st2(lwheel, sector);
       if ( stat == 3 ) masked = msks.get_etsoc_chdis_st3(lwheel, sector);
 
       if ( !masked ) m_mask = false;
 
       if ( tseta && !masked ) {
 
         if ( wheel == -2 || wheel == -1 ||
              ( wheel == 0 && (sector == 0 || sector == 3 || sector == 4 || sector == 7 || sector == 8 || sector == 11) ) ) {
           for ( int i = 0; i < 7; i++ ) {
             if ( tseta->position(i) ) pos.set(6-i);
             if ( tseta->quality(i) ) qual.set(6-i);
           }
         } else {
           for ( int i = 0; i < 7; i++ ) {
             if ( tseta->position(i) ) pos.set(i);
             if ( tseta->quality(i) ) qual.set(i);
           }
         }
       }
 
       const L1MuBMTrackSegEta* tmpts = new L1MuBMTrackSegEta(wheel,sector,stat,pos.to_ulong(),qual.to_ulong(),bx-bx_offset);
       m_tseta.push_back(tmpts);
     }
   }
 
 }

void L1MuBMEtaProcessor::reset ( void )

virtual

reset the Eta Processor

Definition at line 105 of file L1MuBMEtaProcessor.cc.

References mps_fire::i, m_address, m_eta, m_fine, m_foundPattern, m_mask, m_pattern, m_TrackCand, m_TracKCand, and m_tseta.

Referenced by id().

                                {
 
   vector<const L1MuBMTrackSegEta*>::iterator iter = m_tseta.begin();
   while ( iter != m_tseta.end() ) {
     if ( *iter ) {
       delete *iter;
       *iter = nullptr;
     }
     iter++;
   }
 
   m_tseta.clear();
 
   for ( int i = 0; i < 12; i++ ) {
     m_eta[i] = 99;
     m_fine[i] = false;
     m_pattern[i] = 0;
     m_address[i] = 0;
     m_TrackCand[i] = nullptr;
     m_TracKCand[i] = nullptr;
   }
 
   m_foundPattern.clear();
 
   m_mask = true;
 
 }

void L1MuBMEtaProcessor::run	(	int	bx,
		const edm::Event &	e,
		const edm::EventSetup &	c
	)

virtual

run the Eta Processor

Definition at line 87 of file L1MuBMEtaProcessor.cc.

References assign(), L1MuBMTFConfig::getEtaTF(), receiveAddresses(), receiveData(), runEtaMatchingUnit(), and runEtaTrackFinder().

Referenced by id().

                                                                               {
 
   if ( L1MuBMTFConfig::getEtaTF() ) {
     receiveData(bx,e,c);
     runEtaTrackFinder(c);
   }
 
   receiveAddresses();
   runEtaMatchingUnit(c);
 
   assign();
 
 }

void L1MuBMEtaProcessor::runEtaMatchingUnit ( const edm::EventSetup & c )

private

run Eta Matching Unit (EMU)

Definition at line 364 of file L1MuBMEtaProcessor.cc.

References bmtfParamsHandle, L1MuDTEtaPattern::eta(), spr::find(), edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), L1MuBMTQualPatternLut::getCoarseEta(), L1MuBMTEtaPatternLut::getPattern(), L1MuBMTQualPatternLut::getQualifiedPatterns(), mps_fire::i, L1TMuonBarrelParams::lutparams_, m_address, m_eta, m_fine, m_foundPattern, L1MuBMTQualPatternLut::m_lut, m_mask, m_pattern, AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< T >::product(), L1TMuonBarrelParams::LUTParams::qp_lut_, theEtaPatternLUT, and theQualPatternLUT.

Referenced by fineBit(), and run().

                                                                   {
 
   //c.get< L1MuDTQualPatternLutRcd >().get( theQualPatternLUT );
   const L1TMuonBarrelParamsRcd& bmtfParamsRcd = c.get<L1TMuonBarrelParamsRcd>();
   bmtfParamsRcd.get(bmtfParamsHandle);
   const L1TMuonBarrelParams& bmtfParams = *bmtfParamsHandle.product();
   theQualPatternLUT.m_lut = bmtfParams.lutparams_.qp_lut_; //l1mudttfqualplut; // EMU look-up tables
 
   // loop over all addresses
   for ( int i = 0; i < 12; i++ ) {
 
     int adr = m_address[i];
     if ( adr == 0 ) continue;
     int sp = i/2 + 1;       //sector processor [1,6]
 
     // assign coarse eta value
     if ( !m_mask ) m_eta[i] = theQualPatternLUT.getCoarseEta(sp,adr);
 
 
     if ( m_foundPattern.empty() ) continue;
 
     // get list of qualified patterns ordered by quality
     // and compare with found patterns
     const vector<short>& qualifiedPatterns = theQualPatternLUT.getQualifiedPatterns(sp,adr);
     vector<short>::const_iterator iter;
     vector<int>::const_iterator f_iter;
     for ( iter = qualifiedPatterns.begin(); iter != qualifiedPatterns.end(); iter++ ) {
       f_iter = find(m_foundPattern.begin(),m_foundPattern.end(),(*iter));
       // found
       if ( f_iter != m_foundPattern.end() ) {
         const L1MuDTEtaPattern p = theEtaPatternLUT.getPattern(*f_iter);
         // assign fine eta value
         m_fine[i] = true;
         m_eta[i]  = p.eta();  // improved eta
 ;
         m_pattern[i] = (*f_iter);
         break;
       }
     }
 
   }
 
   // if both tracks from one sector processor deliver the same track address
   // both tracks get only a coarse eta value!
 
   // loop over sector processors
   for ( int i = 0; i < 6; i++ ) {
     int idx1 = 2*i;
     int idx2 = 2*i+1;
     int adr1 = m_address[idx1];
     int adr2 = m_address[idx2];
     if ( adr1 == 0 || adr2 == 0 ) continue;
     if ( adr1 == adr2 && !m_mask ) {
       // both tracks get coarse (default) eta value
       m_eta[idx1]  = theQualPatternLUT.getCoarseEta(i+1,adr1);
       m_pattern[idx1] = 0;
       m_fine[idx1] = false;
       m_eta[idx2]  = theQualPatternLUT.getCoarseEta(i+1,adr2);
 
       m_pattern[idx2] = 0;
       m_fine[idx2] = false;
     }
   }
 
 }

void L1MuBMEtaProcessor::runEtaTrackFinder ( const edm::EventSetup & c )

private

run Eta Track Finder (ETF)

Definition at line 313 of file L1MuBMEtaProcessor.cc.

References L1MuBMTEtaPatternLut::begin(), stringResolutionProvider_cfi::bin, bmtfParamsHandle, relativeConstraints::empty, L1MuBMTEtaPatternLut::end(), L1TMuonBarrelParams::LUTParams::eta_lut_, edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), mps_fire::i, L1MuDTEtaPattern::id(), L1TMuonBarrelParams::lutparams_, m_foundPattern, L1MuBMTEtaPatternLut::m_lut, m_tseta, listBenchmarks::pattern, L1MuDTEtaPattern::position(), edm::ESHandle< T >::product(), lumiQueryAPI::q, L1MuDTEtaPattern::quality(), quality(), relativeConstraints::station, theEtaPatternLUT, L1MuDTEtaPattern::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by fineBit(), and run().

                                                                  {
 
   //c.get< L1MuDTEtaPatternLutRcd >().get( theEtaPatternLUT );
   const L1TMuonBarrelParamsRcd& bmtfParamsRcd = c.get<L1TMuonBarrelParamsRcd>();
   bmtfParamsRcd.get(bmtfParamsHandle);
   const L1TMuonBarrelParams& bmtfParams = *bmtfParamsHandle.product();
   theEtaPatternLUT.m_lut  = bmtfParams.lutparams_.eta_lut_; //l1mudttfetaplut;  // ETF look-up table
 
   // check if there are any data
   bool empty = true;
   for ( int i = 0; i < 15; i++ ) {
     empty &= m_tseta[i]->empty();
   }
   if ( empty ) return;
 
   // Pattern comparator:
   // loop over all patterns and compare with local chamber pattern
   // result : list of valid pattern IDs ( m_foundPattern )
   L1MuBMTEtaPatternLut::ETFLut_iter it = theEtaPatternLUT.begin();
   while ( it != theEtaPatternLUT.end() ) {
 
     const L1MuDTEtaPattern pattern = (*it).second;
     int qualitycode = pattern.quality();
 
     bool good = true;
 
     for ( int station = 0; station < 3; station++) {
       int q = quality(qualitycode,station+1);
       int wheel = pattern.wheel(station+1);
       int bin = pattern.position(station+1);
       if ( bin == 0 ) continue;
       bitset<7> pos  = m_tseta[wheel+2 + station*5]->position();
       bitset<7> qual = m_tseta[wheel+2 + station*5]->quality();
       // compare position
       good &= pos.test(bin-1);
       // compare quality
       if ( q == 2 ) good &= qual.test(bin-1);
     }
 
     if ( good ) m_foundPattern.push_back(pattern.id());
 
     it++;
 
   }
 
 }