PatternRecognition.h File Reference

#include "L1Trigger/L1TMuonEndCap/interface/PhiMemoryImage.h"
#include "L1Trigger/L1TMuonEndCap/interface/EmulatorClasses.h"

Go to the source code of this file.

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

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
      
    }
  }
}