Deltas.h

Go to the documentation of this file.

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



DeltaOutput Deltas(MatchingOutput Mout, int zone, int winner){

    //bool verbose = false;
    
    PhOutput phmatch = Mout.PhiMatch();
    ThOutput thmatch = Mout.ThetaMatch();
    
    /*for comparison only
    for(int xx=0;xx<4;xx++){
    
        if(phmatch[zone][winner][xx].Phi() != -999 && verbose)
            std::cout<<"phmatch["<<zone<<"]["<<winner<<"]["<<xx<<"] = "<<phmatch[zone][winner][xx].Phi()<<"\n\n";
    }
    */
    
    int dphi[6] = {-999,-999,-999,-999,-999,-999};
    int dtmp[6] = {-999,-999,-999,-999,-999,-999};
    int dtmpi[6] = {-999,-999,-999,-999,-999,-999};
    int dth[6][4] = {{-999,-999,-999,-999},{-999,-999,-999,-999},{-999,-999,-999,-999},{-999,-999,-999,-999},{-999,-999,-999,-999},{-999,-999,-999,-999}};
    

    for(int s1=0;s1<3;s1++){
    
        for(int s2=s1+1;s2<4;s2++){
        
                    
            if((s1 == 0) && (phmatch[zone][winner][s1].Phi() != -999) && (phmatch[zone][winner][s2].Phi() != -999)){ //if using station one and both hits are valid
                    
                dphi[s2-1] = phmatch[zone][winner][s1].Phi() - phmatch[zone][winner][s2].Phi();
            }
            else if((s1 != 0) && (phmatch[zone][winner][s1].Phi() != -999) && (phmatch[zone][winner][s2].Phi() != -999)){//if not using station one and both hits are valid
                    
                dphi[s1+s2] = phmatch[zone][winner][s1].Phi() - phmatch[zone][winner][s2].Phi();    
            }
            
            
                    
            for(int i=0;i<2;i++){
                        
                    
                for(int j=0;j<2;j++){
                        
                    int thi = thmatch[zone][winner][s1][i].Theta();
                    int thj = thmatch[zone][winner][s2][j].Theta();
                    int deltath = thi - thj;
                    
                    
                    if((s1 == 0) && (thi != -999) && (thj != -999)){
                                
                        if(!i){dth[s2-1][i+j] = deltath;}
                        else{dth[s2-1][i+j+1] = deltath;}
                                
                                
                    }
                    else if((s1 != 0) && (thi != -999) && (thj != -999)){
                            
                        if(!i){dth[s1+s2][i+j] = deltath;}
                        else{dth[s1+s2][i+j+1] = deltath;}
                    }
                            
                }
            }
        
        }
    } 
    
    
    int vmask[3] = {0,0,0};
    const unsigned int mask[6] = {0x3,0x5,0x9,0x6,0xa,0xc};
    const unsigned int mindex[6] = {0,0,0,1,1,2};//vmasindex-->says that the first three entries of mask belong to vmask[0] etc...
    
    for(int p=0;p<6;p++){
    
        //if(dphi[p] != -999 && verbose)
        //  std::cout<<"dphi["<<p<<"] = "<<dphi[p]<<"\n\n";
    
        for(int l=0;l<4;l++){
        
            //if(dth[p][l] != -999 && verbose){std::cout<<"dth["<<p<<"]["<<l<<"] = "<<dth[p][l]<<"\n\n";}
        
            if(abs(dth[p][l]) < fabs(dtmp[p])){//get best dtheta(i.e. the smallest)
            
                //if(verbose) std::cout<<"chose second hit theta  index-"<<p<<std::endl;
            
                dtmp[p] = dth[p][l];
                
                dtmpi[p] = l;//says which combination of dth is the one to choose (because there are 4 possible as stated above^)
            }
        }
        
        
        if((abs(dtmp[p]) <= 4) && (dtmp[p] != -999)){//if dtheta is small enought and valid
        
            //std::cout<<"valid "<<p<<std::endl;
        
            vmask[mindex[p]] |= mask[p];    
        } 
    }
    
    //for(int q=0;q<3;q++){
    //  if(vmask[q] && verbose)
    //      std::cout<<"vmask["<<q<<"] = "<<vmask[q]<<std::endl;
    //}
    
    
    unsigned int vstat = vmask[0];
    
    //if(vstat && verbose){std::cout<<"vstat = "<<vstat<<std::endl;}
    
    if( !vstat || (vstat & vmask[1])){vstat |= vmask[1];}
    
    
    if( !vstat || (vstat & vmask[2])){vstat |= vmask[2];}
    
    //if(vstat && verbose){std::cout<<"vstat = "<<vstat<<std::endl;}
    
    const unsigned int vstatindex[11] = {0xc,0xa,0x6,0xe,0x9,0x5,0xd,0x3,0xb,0x7,0xf};
    const unsigned int viindex[2][11] = {{5,4,3,3,2,1,1,0,0,0,0},{5,4,3,5,2,1,5,0,4,3,3}};
    
    std::vector<int> d (2,-999);
    std::vector<std::vector<int>> deltas (2,d);//deltas[0]->dPhi & deltas[1]->dTheta
    

    
    for(int c=0;c<11;c++){
    
        if(vstat == vstatindex[c]){
        
            deltas[0][0] = dphi[viindex[0][c]];
            deltas[0][1] = dphi[viindex[1][c]];
            deltas[1][0] = dtmp[viindex[0][c]];
            deltas[1][1] = dtmp[viindex[1][c]];
        }
    }
    
    int phi = 0, theta = 0, id = 0;
    if(vstat & 2){//ME2 Present
        
        //phi is simple, we have only one per station to choose from
        phi = phmatch[zone][winner][1].Phi();
        
        //for theta, select delta to best station, use dtmpi as index
        if(dtmp[0] != -999){
        
            if(dtmpi[0] < 2)
                id = 1;
            
            
            theta = thmatch[zone][winner][1][id].Theta();
        }
        else if(dtmp[3] != -999){
            
            if(dtmpi[3] > 1)
                id = 1;
            
            
            theta = thmatch[zone][winner][1][id].Theta();
        }
        else if(dtmp[4] != -999){
            
            if(dtmpi[4] > 1)
                id = 1;
            
            
            theta = thmatch[zone][winner][1][id].Theta();
        }
    
    }
    else if(vstat & 4){//ME3 Present, but not ME2
    
        phi = phmatch[zone][winner][2].Phi();
        
        if(dtmp[1] != -999){
        
            if(dtmpi[1] < 2)
                id = 1;
            
            
            theta = thmatch[zone][winner][2][id].Theta();
        }
        else if(dtmp[5] != -999){
            
            if(dtmpi[5] > 1)
                id = 1;
            
            
            theta = thmatch[zone][winner][2][id].Theta();
        }
    }
    else if(vstat & 8){//ME4 Present but not ME2 or ME3
    
        phi = phmatch[zone][winner][3].Phi();
        if(dtmp[2] != -999){
        
            if(dtmpi[2] < 2)
                id = 1;
            
            
            theta = thmatch[zone][winner][3][id].Theta();
        }
    }
    
    
    
    int rank = (Mout.Winners()[zone][winner].Rank()<<1);
    
    
    if(vstat & 8){rank |= 1;}else{rank &= 0x7e;}//if station 4 add to rank, else keep everythign and zero the bit which indicates station 4 is present
    if(vstat & 4){rank |= 2;}else{rank &= 0x7d;}//if station 3 add to rank, else keep everythign and zero the bit which indicates station 3 is present
    if(vstat & 2){rank |= 8;}else{rank &= 0x77;}//if station 2 add to rank, else keep everythign and zero the bit which indicates station 2 is present
    if(vstat & 1){rank |= 32;}else{rank &= 0x5f;}//if station 1 add to rank, else keep everythign and zero the bit which indicates station 1 is present
    
    if(vstat == 0 || vstat == 1 || vstat == 2 || vstat == 4 || vstat == 8){rank = 0;}//removes single, and ME3--ME4 hit combinations
    
    DeltaOutput output;
    Mout.SetPhOut(phmatch);
    Winner win = Mout.Winners()[zone][winner];
    win.SetRank(rank);
    
    output.SetValues(Mout,deltas,phi,theta,win);
    
    return output;


}

std::vector<std::vector<DeltaOutput>> CalcDeltas(MatchingOutput Mout){

    DeltaOutput output;output.SetNull();
    
    std::vector<DeltaOutput> o (3,output);
    std::vector<std::vector<DeltaOutput>> out (4,o);
    
    for(int zone=0;zone<4;zone++){
    
        for(int winner=0;winner<3;winner++){
            
            
            out[zone][winner] = Deltas(Mout, zone, winner);
        }
    }   
    
    
    
    return out;
}