/data/refman/pasoursint/CMSSW_5_3_0/src/L1Trigger/RegionalCaloTrigger/src/L1RCTORCAMap.cc

Go to the documentation of this file.

#include "L1Trigger/RegionalCaloTrigger/interface/L1RCTORCAMap.h"

#include <vector>
using std::vector;

L1RCTORCAMap::L1RCTORCAMap(){
  rawEMET = std::vector<unsigned short>(56*72);
  rawHDET = std::vector<unsigned short>(56*72);
  rawEMFG = std::vector<unsigned short>(56*72);
  rawHDFG = std::vector<unsigned short>(56*72);
  rawHFET = std::vector<unsigned short>(8*18);
  combEM = std::vector<unsigned short>(56*72);
  combHD = std::vector<unsigned short>(56*72);

  barrelData = std::vector<std::vector<std::vector<unsigned short> > >(18,std::vector<std::vector<unsigned short> >(7,
                                                           std::vector<unsigned short>(64)));
  hfData = std::vector<std::vector<unsigned short> >(18,std::vector<unsigned short>(8));
  
}
vector<std::vector<std::vector<unsigned short> > > L1RCTORCAMap::giveBarrel(){
  return barrelData;
}

vector<std::vector<unsigned short> > L1RCTORCAMap::giveHF(){
  return hfData;
}
void L1RCTORCAMap::makeHFData(){
  std::vector<unsigned short> crate(16);
  for(int i = 0; i< 9; i++){
    for(int j = 0; j<16; j++){
      crate.at(j) = rawHFET.at(16*i+j);
    }
    hfData.at(i) = crate;
  }  
}  

void L1RCTORCAMap::makeBarrelData(){
  std::vector<int> indices;
  for(int phi = 0; phi<72;phi++){
    for(int eta = 0; eta<56;eta++){
      indices = orcamap(eta,phi);
      (barrelData.at(indices.at(0))).at(indices.at(1)).at(indices.at(2)) = combEM.at(phi*56 + eta);
      (barrelData.at(indices.at(0))).at(indices.at(1)).at(indices.at(2)+32) = combHD.at(phi*56 + eta);
    }
  }
}

void L1RCTORCAMap::readData(std::vector<unsigned> emet, std::vector<unsigned>  hdet,
                            std::vector<bool> emfg, std::vector<bool> hdfg,
                            std::vector<unsigned> hfet){
  for(int i = 0; i<4032; i++){
    rawEMET.at(i) = emet.at(i);
    rawHDET.at(i) = hdet.at(i);
    rawEMFG.at(i) = emfg.at(i);
    rawHDFG.at(i) = hdfg.at(i);
  }
  for(int i = 0; i<144; i++)
    rawHFET.at(i) = hfet.at(i);

  combEM = combVec(rawEMET,rawEMFG);
  combHD = combVec(rawHDET,rawHDFG);

  makeBarrelData();
  makeHFData();
}

unsigned short L1RCTORCAMap::combine(unsigned short et, unsigned short fg){
  unsigned short newfg = fg << 8;
  return newfg + et;
}

vector<unsigned short> L1RCTORCAMap::combVec(std::vector<unsigned short> et, 
                                             std::vector<unsigned short> fg){
  std::vector<unsigned short> comb(56*72);
  for(int i = 0; i<(int)et.size(); i++)
    comb.at(i) = combine(et.at(i),fg.at(i));
  return comb;
}

vector<int> L1RCTORCAMap::orcamap(int eta, int phi){
  int crateNum(20);
  std::vector<int> cardTower(2,0);
  std::vector<int> returnVec(3,0);
  int modEta = eta%28;
  int modPhi = phi%8;
  if(phi < 8)
    crateNum = 0;  
  else if( phi < 16)
    crateNum = 1;
  else if( phi < 24)
    crateNum = 2;
  else if( phi < 32)
    crateNum = 3;
  else if( phi < 40)
    crateNum = 4;
  else if( phi < 48)
    crateNum = 5;
  else if( phi < 56)
    crateNum = 6;
  else if( phi < 64)
    crateNum = 7;
  else if( phi < 72)
    crateNum = 8;

  if(eta < 28)
    cardTower = lowEtaMap(modEta,modPhi);
  else {
    cardTower = highEtaMap(modEta,modPhi);
    crateNum = crateNum + 9;
  }

  returnVec.at(0) = crateNum;
  for(int i =0; i<2; i++){
    returnVec.at(i+1) = cardTower.at(i);
  }
  return returnVec;
}
    

vector<int> L1RCTORCAMap::lowEtaMap(int eta, int phi){
  int cardnum = 0;
  int towernum = 0;
  std::vector<int> returnVec(2);
  if(eta < 4){
    cardnum = 6;
    if(phi < 4)
      towernum = (3-eta)*4 + phi;
    else
      towernum = eta*4 + phi + 12;
  }
  else if(eta <12){
    if(phi < 4){
      cardnum = 4;
      if(eta < 8)
        towernum = (7-eta)*4+phi+16;
      else
        towernum = (11-eta)*4+phi;
    }
    else{
      cardnum = 5;
      if(eta < 8)
        towernum = (7-eta)*4+(phi-4)+16;
      else
        towernum = (11-eta)*4+(phi-4);
    }
  }
  else if(eta < 20){
    if(phi < 4){
      cardnum = 2;
      if(eta < 16)
        towernum = (15-eta)*4+phi+16;
      else
        towernum = (19-eta)*4+phi;
    }
    else{
      cardnum = 3;
      if(eta < 16)
        towernum = (15-eta)*4+(phi-4)+16;
      else
        towernum = (19-eta)*4+(phi-4);
    }
  }
  else if(eta < 28){
    if(phi < 4){
      cardnum = 0;
      if(eta < 24)
        towernum = (23-eta)*4+phi+16;
      else
        towernum = (27-eta)*4+phi;
    }
    else{
      cardnum = 1;
      if(eta < 24)
        towernum = (23-eta)*4+(phi-4)+16;
      else
        towernum = (27-eta)*4+(phi-4);
    }
  }
  
  returnVec.at(0) = cardnum;
  returnVec.at(1) = towernum;
  return returnVec;
}


vector<int> L1RCTORCAMap::highEtaMap(int eta, int phi){
  int cardnum;
  int towernum;
  std::vector<int> returnVec(2);
  if(eta < 8){
    if(phi < 4){
      cardnum = 0;
      towernum = eta*4+phi;
    }
    else{
      cardnum = 1;
      towernum = eta*4+(phi-4);
    }
  }
  else if(eta < 16){
    if(phi < 4){
      cardnum = 2;
      towernum = (eta-8)*4+phi;
    }
    else{
      cardnum = 3;
      towernum = (eta-8)*4+(phi-4);
    }
  }
  else if(eta < 24){
    if(phi < 4){
      cardnum = 4;
      towernum = (eta-16)*4+phi;
    }
    else{
      cardnum = 5;
      towernum = (eta-16)*4+(phi-4);
    }
  }
  else{
    cardnum = 6;
    if(phi < 4)
      towernum = (27-eta)*4+phi;
    else
      towernum = (27-eta)*4+phi+12;
  }

  returnVec.at(0)=cardnum;
  returnVec.at(1)=towernum;
  return returnVec;
}