PatternRecognition.cc File Reference

#include "L1Trigger/L1TMuonEndCap/interface/PatternRecognition.h"

Go to the source code of this file.

Macros
#define	PATTERN_SIZE   9

Functions
PatternOutput	DetectPatterns (ZonesOutput Eout)
std::vector< PatternOutput >	Patterns (std::vector< ZonesOutput > Zones)
void	PrintQuality (QualityOutput out)

Variables
PhiMemoryImage	pattern1 (0x8000, 0, 0, 0x8000, 0, 0, 0x8000, 0, 0, 0x8000, 0, 0)
PhiMemoryImage	pattern2 (0x10000, 0, 0, 0x8000, 0, 0, 0xc000, 0, 0, 0xc000, 0, 0)
PhiMemoryImage	pattern3 (0x4000, 0, 0, 0x8000, 0, 0, 0x18000, 0, 0, 0x18000, 0, 0)
PhiMemoryImage	pattern4 (0x60000, 0, 0, 0x8000, 0, 0, 0xe000, 0, 0, 0xe000, 0, 0)
PhiMemoryImage	pattern5 (0x3000, 0, 0, 0x8000, 0, 0, 0x78000, 0, 0, 0x78000, 0, 0)
PhiMemoryImage	pattern6 (0x780000, 0, 0, 0x8000, 0, 0, 0xff00, 0, 0, 0xff00, 0, 0)
PhiMemoryImage	pattern7 (0xf00, 0, 0, 0x8000, 0, 0, 0x7f8000, 0, 0, 0x7f8000, 0, 0)
PhiMemoryImage	pattern8 (0x7f800000, 0, 0, 0x8000, 0, 0, 0xff00, 0, 0, 0xff00, 0, 0)
PhiMemoryImage	pattern9 (0xff, 0, 0, 0x8000, 0, 0, 0x7f8000, 0, 0, 0x7f8000, 0, 0)
PhiMemoryImage	patterns [9] = {pattern8, pattern9, pattern6, pattern7, pattern4, pattern5, pattern2, pattern3, pattern1}

Macro Definition Documentation

#define PATTERN_SIZE   9

Definition at line 12 of file PatternRecognition.cc.

Referenced by DetectPatterns().

Function Documentation

PatternOutput DetectPatterns

(

ZonesOutput

Eout

)

variable declarations///

Definition at line 24 of file PatternRecognition.cc.

References a, b, PhiMemoryImage::BitShift(), ZonesOutput::convertedhits, PatternOutput::detected, PatternOutput::hits, i, relval_2017::k, QualityOutput::layer, convertSQLitetoXML_cfg::output, patt, PATTERN_SIZE, QualityOutput::rank, QualityOutput::straightness, tmp, detailsBasic3DVector::y, and ZonesOutput::zone.

Referenced by Patterns().

                                              {
  
  //bool verbose = false;
  std::vector<int> tmp (192, 0);//was 128
  std::vector<std::vector<int>> lya (4, tmp), stra (4, tmp), ranka_t (4, tmp), ranka (4, tmp);
  std::vector<PhiMemoryImage> Merged = Eout.zone;
  
  
  for(int zone=0;zone<4;zone++){
  
    for(int b=0;b<192;b++){//loop over stips of detector zones//was 128 now 192 to accomodate 
      //larger phi scale used in neighboring sectors algorithm
      int ly[PATTERN_SIZE] = {0}, srt[PATTERN_SIZE] = {0}, qu[PATTERN_SIZE] = {0};
      PhiMemoryImage patt[PATTERN_SIZE];
      for (int i=0; i<PATTERN_SIZE; i++) patt[i] = patterns[i];
      for(int y=0;y != PATTERN_SIZE;y++){//loop over patterns
    
    bool zona[12] = {false}; //Clear out station presence
    
    if((b-15) < 63){  
      patt[y].BitShift(b-15);  
    }   
    else if((b-15) < 127){  
      patt[y].BitShift(63);patt[y].BitShift(1);patt[y].BitShift(b-79);
    }   
    else{   
      patt[y].BitShift(63);patt[y].BitShift(1);patt[y].BitShift(63);patt[y].BitShift(1);patt[y].BitShift(b-143);  
    }           
    
    for(int yy=0;yy != 12;yy++){//loop over 8 long integers of each pattern
      zona[yy] = patt[y][yy] & Merged[zone][yy];
    }
    
    if(zona[0] || zona[1] || zona[2]){ly[y] += 4;}//if station 1 is present
    if(zona[3] || zona[4] || zona[5]){ly[y] += 2;}//if station 2 is present
    if(zona[6] || zona[7] || zona[8] || zona[9]  || zona[10] || zona[11]){ly[y] += 1;}//if station 3 or 4 is present
    srt[y] = y/2;//straightness code::pattern specific
    
    if(  (ly[y] != 0) && (ly[y] != 1) && (ly[y] != 2) && (ly[y] != 4)  ){//Removes single and ME34 hit combinations//
      
      //creates quality code by interleving straightness and layer bit words///
      if((srt[y] & 1) != 0){qu[y] += 2;}
      if((srt[y] & 2) != 0){qu[y] += 8;}
      if((srt[y] & 4) != 0){qu[y] += 32;}
      if((ly[y] & 1) != 0){qu[y] += 1;}
      if((ly[y] & 2) != 0){qu[y] += 4;}
      if((ly[y] & 4) != 0){qu[y] += 16;}
      
    }
    
      }//pattern loop
      
      for(int a=0;a<PATTERN_SIZE;a++)
    {
      if(qu[a] > ranka_t[zone][b]){//take highest quality from all patterns at given key strip
        
        ranka_t[zone][b] = qu[a];
        lya[zone][b] = ly[a];
        stra[zone][b] = srt[a];
      }
    }//pattern loop to assign quality
      
    }//strip loop
  }//zone loop
  
  
  for(int zone=0;zone<4;zone++){
    
    for(int k=0;k<192;k++){//was 128
      
      int qr = 0, ql = 0, qc = ranka_t[zone][k];
      
      if(k>0){qr=ranka_t[zone][k-1];}
      if(k<191){ql=ranka_t[zone][k+1];}//was 127

      
      if((qc <= ql) || (qc < qr)){qc = 0;}
      
      ranka[zone][k] = qc;
    }
  }
  
  QualityOutput qout;
  qout.rank = ranka;
  qout.layer = lya;
  qout.straightness = stra;
  
  PatternOutput output;
  output.detected = qout;
  output.hits = Eout.convertedhits;
  
  return output;
  
}

std::vector<PatternOutput> Patterns

(

std::vector< ZonesOutput >

Zones

)

Definition at line 125 of file PatternRecognition.cc.

References DetectPatterns(), i, convertSQLitetoXML_cfg::output, and tmp.

Referenced by L1TMuonEndCapTrackProducer::produce().

                                                               {
 
  PatternOutput tmp;
  std::vector<PatternOutput> output (3,tmp);
  
  for(int i=0;i<3;i++)
    output[i] = DetectPatterns(Zones[i]);
  
  return output;
  
}

void PrintQuality

(

QualityOutput

out

)

Print Quality for comparison ///

Definition at line 138 of file PatternRecognition.cc.

References gather_cfg::cout, QualityOutput::layer, QualityOutput::rank, QualityOutput::straightness, and detailsBasic3DVector::z.

                                     {
  
  for(int zz=0;zz<4;zz++){
    for(int z = 0;z<192;z++){//was 128
      
      if((out.rank)[zz][z] && ((out.layer)[zz][z] || (out.straightness)[zz][z])){std::cout<<"Found Pattern Zone: "<<zz<<"::new "<<(z+1)<<": "<<(out.layer)[zz][z]<<", "<<(out.straightness)[zz][z]<<", "<<(out.rank)[zz][z]<<" ";std::cout<<"\n\n";}//changing zones to Merged
      
    }
  }
}

Variable Documentation

PhiMemoryImage pattern1(0x8000, 0, 0, 0x8000, 0, 0, 0x8000, 0, 0, 0x8000, 0, 0)

Referenced by CSCCLCTDigi::operator>().

PhiMemoryImage pattern2(0x10000, 0, 0, 0x8000, 0, 0, 0xc000, 0, 0, 0xc000, 0, 0)

Referenced by CSCCLCTDigi::operator>().

PhiMemoryImage pattern3(0x4000, 0, 0, 0x8000, 0, 0, 0x18000, 0, 0, 0x18000, 0, 0)

PhiMemoryImage pattern4(0x60000, 0, 0, 0x8000, 0, 0, 0xe000, 0, 0, 0xe000, 0, 0)

PhiMemoryImage pattern5(0x3000, 0, 0, 0x8000, 0, 0, 0x78000, 0, 0, 0x78000, 0, 0)

PhiMemoryImage pattern6(0x780000, 0, 0, 0x8000, 0, 0, 0xff00, 0, 0, 0xff00, 0, 0)

PhiMemoryImage pattern7(0xf00, 0, 0, 0x8000, 0, 0, 0x7f8000, 0, 0, 0x7f8000, 0, 0)

PhiMemoryImage pattern8(0x7f800000, 0, 0, 0x8000, 0, 0, 0xff00, 0, 0, 0xff00, 0, 0)

PhiMemoryImage pattern9(0xff, 0, 0, 0x8000, 0, 0, 0x7f8000, 0, 0, 0x7f8000, 0, 0)

PhiMemoryImage patterns[9] = {pattern8, pattern9, pattern6, pattern7, pattern4, pattern5, pattern2, pattern3, pattern1}

Definition at line 22 of file PatternRecognition.cc.

Referenced by TotemRPDQMSource::analyze(), MuonSeedOrcaPatternRecognition::endcapPatterns(), MuonSeedGenerator::produce(), and TotemRPUVPatternFinder::produce().