PrimitiveConverter_Neighbor.h

Go to the documentation of this file.

#ifndef ADD_PrimitiveConverter_Neighbor
#define ADD_PrimitiveConverter_Neighbor

#include "DataFormats/MuonDetId/interface/CSCDetId.h"

#include "L1Trigger/L1TMuonEndCap/interface/EmulatorClasses.h"
#include "L1Trigger/L1TMuonEndCap/interface/PhThLUT_Neighbor.h"
#include "L1Trigger/L1TMuon/interface/deprecate/MuonTriggerPrimitive.h"

bool neighbor(int endcap,int sector,int SectIndex,int id,int sub,int station){
  
  bool neighbor = false;
  
  int CompIndex = (endcap - 1)*6 + sector - 1;
  bool AdSector = false;
  if((SectIndex == 0 && CompIndex == 5) || 
     (SectIndex == 1 && CompIndex == 0) ||
     (SectIndex == 2 && CompIndex == 1) ||
     (SectIndex == 3 && CompIndex == 2) ||
     (SectIndex == 4 && CompIndex == 3) ||
     (SectIndex == 5 && CompIndex == 4) ||
     (SectIndex == 6 && CompIndex == 11) || 
     (SectIndex == 7 && CompIndex == 6) ||
     (SectIndex == 8 && CompIndex == 7) ||
     (SectIndex == 9 && CompIndex == 8) ||
     (SectIndex == 10 && CompIndex == 9) ||
     (SectIndex == 11 && CompIndex == 10) ){AdSector = true;}
  
  if(AdSector && sub == 2 && station == 1 && (id == 3 || id == 6 || id == 9) )
    neighbor = true;
  
  if(AdSector && station > 1 && (id == 3 || id == 9) )
    neighbor = true;
  
  return neighbor;
  
}

//neighboring chambers//

int ph_offsets_neighbor[5][10][3] = 
  {{{39,39,-99}, {57,57,-99}, {76,76,-99}, {39,-99,-99}, {58,-99,-99}, 
    {76,-99,-99}, {41,-99,-99}, {60,-99,-99}, {79,-99,-99}, {21,21,23}},//not sure if 23 is done right yet
   {{95,95,-99}, {114,114,-99}, {132,132,-99}, {95,-99,-99}, {114,-99,-99}, 
    {133,-99,-99}, {98,-99,-99}, {116,-99,-99}, {135,-99,-99}, {21,21,23}},//not sure if 23 is done right yet
   {{38,38,-99}, {76,76,-99}, {113,113,-99}, {39,39,-99}, {58,58,-99}, 
    {76,76,-99}, {95,95,-99}, {114,114,-99}, {132,132,-99}, {1,21,21}},
   {{38,-99,-99}, {76,-99,-99}, {113,-99,-99}, {39,39,39}, {58,58,58}, 
    {76,76,76}, {95,95,95}, {114,114,114}, {132,132,132}, {1,21,21}},
   {{38,-99,-99}, {76,-99,-99}, {113,-99,-99}, {38,38,-99}, {57,57,-99}, 
    {76,76,-99}, {95,95,-99}, {113,113,-99}, {132,132,-99}, {1,20,20}}};//[station][id][phzvl look up #(-99 indicates invaled entry)]

std::vector<ConvertedHit> PrimConv(std::vector<L1TMuon::TriggerPrimitive> TriggPrim, int SectIndex){
  
  //bool verbose = false;
  
  std::vector<ConvertedHit> ConvHits;
  for(std::vector<L1TMuon::TriggerPrimitive>::iterator C1 = TriggPrim.begin();C1 != TriggPrim.end();C1++){
    
    L1TMuon::TriggerPrimitive C3 = *C1;
    CSCDetId Det = C3.detId<CSCDetId>();
    
    int station = Det.station(), chamber = Det.chamber(), ring = Det.ring(), wire = C3.getCSCData().keywire;
    int sector = Det.triggerSector(), strip = C3.getCSCData().strip, pattern = C3.getPattern(), Id = C3.Id();
    int quality = C3.getCSCData().quality, BX = C3.getCSCData().bx, endcap = Det.endcap();
    
    if(station == 1 && ring == 1 && strip > 127){
      ring = 4;
    }
    
    int sub = 0;
    
    
    if(station == 1)
      {
    if(chamber%6 > 2)
      sub = 1;
    else
      sub = 2;
      }
    
    bool IsNeighbor = neighbor(endcap,sector,SectIndex,Id,sub,station);
    
    if(ring == 4){Id += 9;}
    
    //if(endcap == 1 && sector == 1)//
    if( (SectIndex ==  (endcap - 1)*6 + sector - 1 )  || IsNeighbor )
      {
    
    
    //if(verbose){
    //  std::cout<<"\n\nSECTOR "<<SectIndex<<"\n\n";
    //  std::cout<<"\n\nRING = "<<ring<<"\n\n";
    //}
    
    
    int ph_tmp = -999, th_tmp = -999;
    int clctpatcor = -999, clctpatsign = -999;
    int eightstrip = -999;
    
    int factor =
      (station == 1 && ring == 4) ? 1707://ME1/1a1707?
      (station == 1 && ring == 3) ? 947: //ME1/3//changed Id > 6 to ring == 3.
      (station == 1 && ring == 1) ? 1301://ME1/1b
      1024;//all other chambers
    
    //std::cout<<"factor = "<<factor<<std::endl;
    
    bool ph_reverse = 
      (endcap == 1 && station >= 3) ? 1:
      (endcap == 2 && station < 3) ? 1: 0;
    
    int ph_coverage  = 
      (station <= 1 && Id > 6 && Id < 10) ? 15 : //30 :
      (station >= 2 && Id <= 3) ? 40 : 20; //80 : 40;
    
    int th_coverage =  
      (station == 1 && Id <= 3) ? 45 :
      (station == 1 && Id > 6 && Id < 10) ? 27 :
      (station == 1 && Id > 3) ? 39 :
      (station == 2 && Id <= 3) ? 43 :
      (station == 2 && Id > 3) ? 56 :
      (station == 3 && Id <= 3) ? 34 :
      (station == 3 && Id > 3) ? 52 :
      (station == 4 && Id <= 3) ? 28 :
      (station == 4 && Id > 3) ? 50 : 0;
    
    int ph_zone_bnd1 = 
      (station <= 1 && (Id <= 3 || Id > 9)) ? 41 :
      (station == 2 && Id <= 3) ? 41 :
      (station == 2 && Id >  3) ? 87 :
      (station == 3 && Id >  3) ? 49 :
      (station == 4 && Id >  3) ? 49 : 127;
                                                   
    int ph_zone_bnd2 = (station == 3 && Id >  3) ? 87 : 127;
                        
    int zone_overlap = 2;
        
    int fph = -999, th = -999, ph_hit = -999, phzvl = -999;// th_hit = -999, ////
    
    //if(sub && verbose)
    //  std::cout<<"\nsub = "<<sub<<"\n";

    
    int LUTi = -999;
    int nId = Id;
    if(IsNeighbor){
      
      if(station < 2){
        
        nId = 12 + Id/3;
        if(ring == 4)
          nId ++;
        
      }
      else{
        int xx = Id;
        if(xx > 6)
          xx = 6;
        
        nId =  9 + xx/3;
        
      }
      
    }
    switch(station)
      {
      case 1: 
        switch(sub)
          {
          case 1: LUTi = nId - 1;break;
          case 2: LUTi = 15 + nId;break;
          default:;//std::cout<<"Sub is out of range"<<std::endl;
          }
        break;
      case 2: LUTi = 27 + nId;break;
      case 3: LUTi = 38 + nId;break;
      case 4: LUTi = 49 + nId;break;
      default:;//std::cout<<"station is out of range"<<std::endl;
      }
    
    if(IsNeighbor && station == 1){
      
      switch(sub)
        {
        case 1: LUTi = 15 + nId;break;
        case 2: LUTi = nId - 1;break;
        default:;//std::cout<<"Sub is out of range"<<std::endl;
        }
      
    }
    
    
    clctpatcor = 0;
    clctpatsign = 0;
    
    if(pattern > 0 && pattern < 11){
      
      clctpatsign = ((pattern%2) == 0);
      
      if(pattern >= 2) {clctpatcor = 5;}
      if(pattern >= 6) {clctpatcor = 2;}
      if(pattern == 10) {clctpatcor = 0;clctpatsign = 0;}
    }
    
    
    eightstrip = 8*strip;
    int patcor = clctpatcor;
    
    //std::cout<<"8strip = "<<eightstrip<<"\n";
    
    if(station == 1 || Id > 3){//10 Degree Chambers
      
      eightstrip = (eightstrip>>1);
      patcor = (patcor>>1);
      if(ring == 4 && strip > 127) eightstrip -= 512;
    }
    
    //std::cout<<"8strip = "<<eightstrip<<"\n"; 
    
    if(clctpatsign) patcor = -patcor;
    eightstrip += patcor;
    
    
    ph_tmp = ((eightstrip*factor)>>10);
    int phShift = (ph_tmp>>5);
    int phLow = 0;
    
    if(ph_reverse){
      
      ph_tmp = -ph_tmp;
      phShift = -phShift;
      phLow = ph_coverage;
    }
    
    int phInitIndex = Id;
    //std::cout<<"phInitIndex = "<<phInitIndex<<" and ph_tmp = "<<ph_tmp<<"\n";
    
    if(station == 1){
      
      int neighborId = C3.Id()/3;
      int subId = sub;
      if(IsNeighbor ){
        
        subId = 1;
        phInitIndex = 12 + neighborId;
        if(ring == 4)
          phInitIndex = 16;//phInitIndex++;
        
      }
      
      fph = Ph_Init_Neighbor[SectIndex][subId-1][phInitIndex - 1] + ph_tmp;
      //std::cout<<"ph init = "<<Ph_Init_Neighbor[SectIndex][subId-1][phInitIndex - 1]<<", index = "<<phInitIndex<<", neighborId = "<<neighborId<<", Id = "<<Id<<"\n";
    }
    else{
      
      int neighborId = Id/3;
      if(neighborId > 2)
        neighborId = 2;
      
      if(IsNeighbor)
        phInitIndex = 9 + neighborId;
      
      fph = Ph_Init_Neighbor[SectIndex][station][phInitIndex - 1] + ph_tmp;
    }
    
    //std::cout<<"pl = "<<phLow<<", ps = "<<phShift<<", ph disp = "<<Ph_Disp_Neighbor[SectIndex][LUTi]<<", >>1 = "<<(Ph_Disp_Neighbor[SectIndex][LUTi]>>1)<<", LUTi = "<<LUTi<<"\n";
    ph_hit = phLow + phShift + (Ph_Disp_Neighbor[SectIndex][LUTi]>>1);
    
    
    int index = -999;
    int th_corr = -999; 
    
    int idl = Id;
    
    if(station == 1){
      
      int neighborId = C3.Id()/3;
      int subId = sub;
      if(IsNeighbor){
        
        subId = 1;
        idl = 12 + neighborId;
        if(ring == 4)
          idl = 16;
      }
      
      //std::cout<<"nid = "<<neighborId<<", idl = "<<idl<<"\n";
      
      th_tmp = Th_LUT_St1_Neighbor[subId-1][SectIndex][idl -1][wire];
      //std::cout<<"th_tmpr = "<<th_tmp<<"\n";
    }
    else{
      
      int neighborId = Id/3;
      if(neighborId > 2)
        neighborId = 2;
      
      if(IsNeighbor)
        idl = 9 + neighborId;
      
      th_tmp = Th_LUT_St234_Neighbor[station-2][SectIndex][idl-1][wire];
      //if(verbose) std::cout<<"\n\nth_tmpr = "<<th_tmp<<"\n\n";
    }
    
    th = th_tmp + Th_Init_Neighbor[SectIndex][LUTi];
    int rth = th;
    //std::cout<<"Th_Init_Neighbor["<<SectIndex<<"]["<<LUTi<<"] = "<<Th_Init_Neighbor[SectIndex][LUTi]<<"\n";
    
    if(station == 1 && (ring == 1 || ring == 4) /*&& endcap == 1*/){
      
      index = (wire>>4)*32 + (eightstrip>>4);
      
      int corrIndex = Id;
      int subId = sub;
      //std::cout<<"corrIndex = "<<corrIndex<<"\n";
      if(corrIndex > 3)
        corrIndex -= 9;
      
      //std::cout<<"corrIndex = "<<corrIndex<<"\n";
      
      if(IsNeighbor && ring == 4){
        subId = 1;
        corrIndex++;
      }
      
      //std::cout<<"corrIndex = "<<corrIndex<<"\n";
      
      th_corr = Th_Corr_Neighbor[subId-1][SectIndex][corrIndex-1][index];
      //std::cout<<"th_corr["<<subId-1<<"]["<<SectIndex<<"]["<<corrIndex-1<<"]["<<index<<"] = "<<th_corr<<"\n";
      
      if(ph_reverse) th_corr = -th_corr;
      //std::cout<<"th_tmp = "<<th_tmp<<"\n";
      
      th_tmp += th_corr;                  //add correction to th_tmp
      //std::cout<<"th_tmp = "<<th_tmp<<"\n";
      if(th_tmp < 0 || wire == 0)
        th_tmp = 0;
      
      if(th_tmp > th_coverage)//this is one change that I'm not sure if it does anything good or not
        th_tmp = th_coverage;
      
      th_tmp &= 0x3f;                     //keep only lowest 6 bits
      //std::cout<<"th_tmp = "<<th_tmp<<"\n";
      //std::cout<<"coverage = "<<th_coverage<<"\n";
      
      if(th_tmp <= th_coverage){
        
        //if(ring == 1){LUTi += 9;}  //change because new Verilog3 sp_tf treats ME11b with LUT's of ME11a
        
        th = th_tmp + Th_Init_Neighbor[SectIndex][LUTi];
        //std::cout<<"th_init["<<SectIndex<<"]["<<LUTi<<"] = "<<Th_Init_Neighbor[SectIndex][LUTi]<<"\n";
      }
      else{th = rth;}//was -999
      
    }
    
    
    if(th != -999){
      
      phzvl = 0;
      
      if(th <= (ph_zone_bnd1 + zone_overlap)) 
        {phzvl |= 1;}
      
      if(th > (ph_zone_bnd2 - zone_overlap)) 
        {phzvl |= 4;}
      
      if((th > (ph_zone_bnd1 - zone_overlap)) && (th <= (ph_zone_bnd2 + zone_overlap)))
        {phzvl |= 2;}
    }
    
    
    int zhit = -99, pz = -99;
    //Each hit could go in more than one zone so we make a vector which stores all the zones for which this hit will contribute
    std::vector<int> zonecontribution; 
    
    if(ring == 4){
      Id -= 9;
      
      if(strip < 128)
        strip += 128;
    }
    
    //determination of zone contribution
    if((phzvl & 1) && (Id < 4 || Id > 9)){pz=0;zonecontribution.push_back(0);}
    if((phzvl & 2) && (Id < 4)){pz=1;zonecontribution.push_back(1);}
    if((phzvl & 1) && (Id > 3) && (station > 2)){pz=0;zonecontribution.push_back(1);}
    if((phzvl & 1) && (Id > 3) && (Id < 7) && (station == 1)){pz=0;zonecontribution.push_back(2);}
    if((phzvl & 1) && (Id > 3) && (station == 2)){pz=0;zonecontribution.push_back(2);}
    if((phzvl & 2) && (Id > 3) && (station > 2)){pz=1;zonecontribution.push_back(2);}
    if((phzvl & 1) && (Id > 4) && (station < 2)){pz=0;zonecontribution.push_back(3);}
    if(phzvl & 4){pz=2;zonecontribution.push_back(3);}
    if((phzvl & 2) && (Id > 3) && (station < 3)){pz=1;zonecontribution.push_back(3);}
    
    int phOffIndex = Id;
    if(IsNeighbor)
      phOffIndex = 10;
    
    //applying ph_offsets
    if(sub == 1){
      zhit = ph_hit + ph_offsets_neighbor[station-1][phOffIndex-1][pz];
      //std::cout<<"\nph_hit = "<<ph_hit<<" and ph_offsets_neighbor["<<station-1<<"]["<<phOffIndex-1<<"]["<<pz<<"] = "<<ph_offsets_neighbor[station-1][phOffIndex-1][pz]<<"\n";
    }
    else{
      
      zhit = ph_hit + ph_offsets_neighbor[station][phOffIndex-1][pz];
      //std::cout<<"ph_hit = "<<ph_hit<<" and ph_offsets_neighbor["<<station<<"]["<<phOffIndex-1<<"]["<<pz<<"] = "<<ph_offsets_neighbor[station][phOffIndex-1][pz]<<"\n";
    }
        
    
    
    //if(SectIndex == 8){
    //std::cout<<"phi = "<<fph<<", theta = "<<th<<", ph_hit = "<<ph_hit<<",zhit = "<<zhit<<", station = "<<station<<", ring = "<<ring<<", id = "<<Id<<", sector "<<SectIndex<<",sub = "<<sub<<", strip = "<<strip<<", wire = "<<wire<<", IsNeighbor = "<<IsNeighbor<<"\n";
    
    //std::cout<<BX-3<<" "<<endcap<<" "<<sector<<" "<<sub<<" "<<station<<" 1 "<<quality<<" "<<pattern<<" "<<wire<<" "<<C3.Id()<<" 0 "<<strip<<"\n";
    //}
    
    /* if(station != 1) */
    /*  sub = 1; */
    /* std::cout<<"proper FR[0] = "<<FRLUT[endcap-1][sector-1][station-1][sub-1][Id-1]<<"\n"; */
    /* if(station != 1) */
    /*  sub = 0; */
    
    ConvertedHit Hit;
    
    int in = 0;
    if(IsNeighbor)
      in = 1;
    
    Hit.SetValues(fph,th,ph_hit,phzvl,station,sub,Id,quality,pattern,wire,strip,BX);
    Hit.SetTP(C3);
    Hit.SetZhit(zhit);
    //Hit.SetZoneContribution(zonecontribution); Function doesn't exist in this class...
    Hit.SetSectorIndex(SectIndex);
    Hit.SetNeighbor(in);
    
    if(Hit.Theta() != -999 && Hit.Phi() > 0){//if theta is valid
      ConvHits.push_back(Hit);
      /*if(verbose){    
        std::cout<<"Phzvl() = "<<Hit.Phzvl()<<", ph_hit = "<<Hit.Ph_hit()<<", station = "<<Hit.Station()<<" and id = "<<Hit.Id()<<std::endl;
        std::cout<<"strip = "<<strip<<", wire = "<<wire<<" and zhit = "<<zhit<<std::endl;
        std::cout<<"\n\nIn Zones: ";
        for(std::vector<int>::iterator in = zonecontribution.begin();in!=zonecontribution.end();in++){
        std::cout<<" "<<*in<<" ";
        }
        }*/
    }   
    
      }//if sector == sectIndex
    
  }
  return ConvHits;
  
}

#endif