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(b-78);                        
                }                                                                       
                else{                                                                   
                    patt[y].BitShift(63);patt[y].BitShift(63);patt[y].BitShift(b-141);  
                }                                                                       
            
            
                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 = ranka_t[zone][k-1], ql = ranka_t[zone][k+1], qc = ranka_t[zone][k];
            
            //if(qc && verbose)
            //  std::cout<<"\n"<<k<<":qc = "<<qc<<" straight: "<<stra[zone][k]<<"  lya: "<<lya[zone][k]<<std::endl; 
        
            if(k==0){qr=0;}
            if(k==191){ql=0;}//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 143 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 157 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
  
       }
     }
 }