MakeRegionalCand.h

Go to the documentation of this file.

#ifndef ADD_MakeRegionalCand
#define ADD_MakeRegionalCand


#include "DataFormats/L1TMuon/interface/RegionalMuonCand.h"
#include "DataFormats/L1TMuon/interface/RegionalMuonCandFwd.h"

#include <cmath>

int GMT_eta_from_theta[128] = {239, 235, 233, 230, 227, 224, 222, 219, 217, 214, 212, 210, 207, 205, 203, 201,
                   199, 197, 195, 193, 191, 189, 187, 186, 184, 182, 180, 179, 177, 176, 174, 172,
                   171, 169, 168, 166, 165, 164, 162, 161, 160, 158, 157, 156, 154, 153, 152, 151,
                   149, 148, 147, 146, 145, 143, 142, 141, 140, 139, 138, 137, 136, 135, 134, 133,
                   132, 131, 130, 129, 128, 127, 126, 125, 124, 123, 122, 121, 120, 119, 118, 117,
                   116, 116, 115, 114, 113, 112, 111, 110, 110, 109, 108, 107, 106, 106, 105, 104,
                   103, 102, 102, 101, 100,  99,  99,  98,  97,  96,  96,  95,  94,  93,  93,  92,
                   91,  91,  90,  89,  89,  88,  87,  87,  86,  85,  84,  84,  83,  83,  82,  81};

int GetPackedEta (int theta, int sector) {
  
  if (theta < 0) return 0;
  if (sector > 5) {
    if (theta > 127) return -240;
    else return -1*GMT_eta_from_theta[theta] - 1;
  }
  else {
    if (theta > 127) return  239;
    else return   GMT_eta_from_theta[theta];
  }
}

float GetGlobalEta(int theta, int sector){ // Should use calc_eta to compute - AWB 07.05.16
  
  // Underlying calculation:
  // float theta_angle = (theta*0.2851562 + 8.5)*(pi/180);
  // float eta = (-1)*log(tan(theta_angle/2));
  
  return 0.010875 * GetPackedEta(theta, sector);
}

int GetPackedPhi(int phi){
  
  // Experimentally, this factor must be between 0.979752 and 0.979758 - AWB 11.05.16 
  float phiDeg = (phi / (0.625 * 60)) * 0.979755;
  int PackedPhi = floor( phiDeg - 35 );
  
  return PackedPhi;
  
}


l1t::RegionalMuonCand MakeRegionalCand(float pt, int phi, int theta,
                       int sign, int mode,
                       int trackaddress, int sector){
  
  l1t::RegionalMuonCand Cand;
  
  int iEta = GetPackedEta(theta,sector);
  int iPhi = GetPackedPhi(phi);
  
  l1t::tftype TFtype = l1t::tftype::emtf_pos;
  if(sector > 5){
    TFtype = l1t::tftype::emtf_neg;
    sector -= 6;
  }
  
  // compressed pt = pt*2 (scale) + 1 (iPt = 0 is empty candidate)
  int iPt = pt*2 + 1;
  if(iPt > 511)
    iPt = 511;
  
  if(iPt < 0)
    iPt = 0;
  
  int iQual = -999;
  
  int LSB = mode & 3;
  
  float eta = GetGlobalEta(theta,sector);
  
  if(eta < 1.2){
    
    switch(mode){
    case(15): iQual = 8;break;
    case(14): iQual = 4;break;
    case(13): iQual = 4;break;
    case(12): iQual = 4;break;
    case(11): iQual = 4;break;
    default: iQual = 4;break;
    }
    
  }
  else{
    
    switch(mode){
    case(15): iQual = 12;break;
    case(14): iQual = 12;break;
    case(13): iQual = 12;break;
    case(12): iQual = 8;break;
    case(11): iQual = 12;break;
    case(10): iQual = 8;break;
    case(7): iQual = 8;break;
    default: iQual = 4;break;
    }
    
  }
  iQual |= LSB;
  
  Cand.setHwPt(iPt);
  Cand.setHwEta(iEta);
  Cand.setHwPhi(iPhi);
  Cand.setHwSign(sign);
  Cand.setHwSignValid(1);
  Cand.setHwQual(iQual);
  // jl: FIXME this has to be adapted to the new schema of saving track addresses
  //Cand.setTrackSubAddress(l1t::RegionalMuonCand::kME12, trackaddress&0xf);
  //Cand.setTrackSubAddress(l1t::RegionalMuonCand::kME22, trackaddress>>4);
  Cand.setTFIdentifiers(sector,TFtype);
  
  return Cand;
  
}

#endif