CMS 3D CMS Logo

Public Member Functions | Private Member Functions | Private Attributes

L1RCT Class Reference

#include <L1RCT.h>

List of all members.

Public Member Functions

void digiInput (const EcalTrigPrimDigiCollection &ecalCollection, const HcalTrigPrimDigiCollection &hcalCollection)
unsigned short ecalCompressedET (int crate, int card, int tower)
unsigned short ecalFineGrainBit (int crate, int card, int tower)
void fileInput (const char *filename)
L1CaloEmCollection getIsolatedEGObjects (unsigned crate)
std::vector< unsigned short > getJetRegions (unsigned crate)
L1CaloEmCollection getNonisolatedEGObjects (unsigned crate)
std::vector< L1CaloRegiongetRegions (unsigned crate)
unsigned short hcalCompressedET (int crate, int card, int tower)
unsigned short hcalFineGrainBit (int crate, int card, int tower)
unsigned short hfCompressedET (int crate, int tower)
unsigned short hfFineGrainBit (int crate, int tower)
void input ()
void input (std::vector< std::vector< std::vector< unsigned short > > > barrelIn, std::vector< std::vector< unsigned short > > hfIn)
 L1RCT (const L1RCTLookupTables *rctLookupTables)
void print ()
void printCrate (int i)
void printEIC (int i, int j)
void printEICEdges (int i, int j)
void printJSC ()
void printJSC (int i)
void printRC (int i, int j)
void processEvent ()
void randomInput ()

Private Member Functions

void configureCards ()
 L1RCT ()
void makeCrates ()
void shareNeighbors ()

Private Attributes

std::vector< std::vector
< std::vector< unsigned short > > > 
barrel
std::vector< L1RCTCratecrates
L1RCTRegion empty
std::vector< std::vector
< unsigned short > > 
hf
L1RCTNeighborMap neighborMap
const L1RCTLookupTablesrctLookupTables_

Detailed Description

Definition at line 20 of file L1RCT.h.


Constructor & Destructor Documentation

L1RCT::L1RCT ( const L1RCTLookupTables rctLookupTables)

Definition at line 56 of file L1RCT.cc.

References makeCrates().

                                                     : 
  rctLookupTables_(rctLookupTables),
  empty(),
  neighborMap(),
  barrel(18,std::vector<std::vector<unsigned short> >(7,std::vector<unsigned short>(64))),
  hf(18,std::vector<unsigned short>(8))
{
  makeCrates();
}
L1RCT::L1RCT ( ) [private]

Member Function Documentation

void L1RCT::configureCards ( ) [private]
void L1RCT::digiInput ( const EcalTrigPrimDigiCollection ecalCollection,
const HcalTrigPrimDigiCollection hcalCollection 
)

Definition at line 116 of file L1RCT.cc.

References abs, barrel, L1RCTParameters::calcCard(), L1RCTParameters::calcCrate(), L1RCTParameters::calcTower(), dtNoiseDBValidation_cfg::cerr, gather_cfg::cout, relval_parameters_module::energy, hf, i, input(), createXMLFile::iphi, j, gen::k, rctLookupTables_, L1RCTLookupTables::rctParameters(), and edm::SortedCollection< T, SORT >::size().

Referenced by L1RCTSaveInput::analyze(), L1RCTInputProducer::produce(), and L1RCTProducer::produce().

{
  // fills input vectors with 0's in case ecal or hcal collection not used
  for (int i = 0; i < 18; i++)
    {
      for (int j = 0; j < 7; j++)
        {
          for (int k = 0; k < 64; k++)
            {
              barrel.at(i).at(j).at(k) = 0;
            }
        }
      for (int j = 0; j < 8; j++)
        {
          hf.at(i).at(j) = 0;
        }
    }

  int nEcalDigi = ecalCollection.size();
  if (nEcalDigi>4032) {nEcalDigi=4032;}
  for (int i = 0; i < nEcalDigi; i++){
    short ieta = (short) ecalCollection[i].id().ieta(); 
    // Note absIeta counts from 1-28 (not 0-27)
    unsigned short absIeta = (unsigned short) abs(ieta);
    unsigned short cal_iphi = (unsigned short) ecalCollection[i].id().iphi(); 
    unsigned short iphi = (72 + 18 - cal_iphi) % 72; // transform TOWERS (not regions) into local rct (intuitive) phi bins

    //map digis to crates, cards, and towers
    unsigned short crate = 999, card = 999, tower = 999;
    crate = rctLookupTables_->rctParameters()->calcCrate(iphi, ieta);
    card = rctLookupTables_->rctParameters()->calcCard(iphi, absIeta);
    tower = rctLookupTables_->rctParameters()->calcTower(iphi, absIeta);

    unsigned short energy = ecalCollection[i].compressedEt();
    unsigned short fineGrain = (unsigned short) ecalCollection[i].fineGrain();  // 0 or 1
    unsigned short ecalInput = energy*2 + fineGrain;

    // put input into correct crate/card/tower of barrel
    if ((crate<18) && (card<7) && (tower<32)) {             // changed 64 to 32 Sept. 19 J. Leonard, rm -1 7Nov07
      barrel.at(crate).at(card).at(tower) = ecalInput;        // rm -1
    }
    else { std::cerr << "L1RCT: ecal out of range! tower = " << tower << " iphi is " << iphi << " absieta is " << absIeta << std::endl; }
  }

//same for hcal, once we get the hcal digis, just need to add 32 to towers:
// just copied and pasted and changed names where necessary
  int nHcalDigi = hcalCollection.size();
  //if (nHcalDigi != 4176){ std::cout << "L1RCT: Warning: There are " << nHcalDigi << "hcal digis instead of 4176!" << std::endl;}
  // incl HF 4032 + 144 = 4176
  for (int i = 0; i < nHcalDigi; i++){
    short ieta = (short) hcalCollection[i].id().ieta(); 
    unsigned short absIeta = (unsigned short) abs(ieta);
    unsigned short cal_iphi = (unsigned short) hcalCollection[i].id().iphi();
    // All Hcal primitives (including HF) are reported
    // with phi bin numbering in the range 0-72.
    unsigned short iphi = (72 + 18 - cal_iphi) % 72;
    // transform Hcal TOWERS (1-72)into local rct (intuitive) phi bins (72 bins) 0-71
    // Use local iphi to work out the region and crate (for HB/HE and HF)
    // HF regions need to have local iphi 0-17
    if (absIeta >= 29) {
      iphi = iphi/4;
    }

    //map digis to crates, cards, and towers
    unsigned short crate = 999, card = 999, tower = 999;
    crate = rctLookupTables_->rctParameters()->calcCrate(iphi, ieta);
    if (absIeta < 29){
      card = rctLookupTables_->rctParameters()->calcCard(iphi, absIeta);
    }
    tower = rctLookupTables_->rctParameters()->calcTower(iphi, absIeta);


    unsigned short energy = hcalCollection[i].SOI_compressedEt();     // access only sample of interest
    unsigned short fineGrain = (unsigned short) hcalCollection[i].SOI_fineGrain();
    unsigned short hcalInput = energy*2 + fineGrain;
    if (absIeta <= 28){
      // put input into correct crate/card/tower of barrel
      if ((crate<18) && (card<7) && (tower<32)) {               // changed 64 to 32 Sept. 19 J. Leonard, rm -1 7Nov07
        barrel.at(crate).at(card).at(tower + 32) = hcalInput;  // hcal towers are ecal + 32 see RC.cc; rm -1 7Nov07
      }
      else { std::cout << "L1RCT: hcal out of range!  tower = " << tower << std::endl; }
    }
    else if ((absIeta >= 29) && (absIeta <= 32)){
      // put input into correct crate/region of HF
      if ((crate<18) && (tower<8)) {
        hf.at(crate).at(tower) = hcalInput;
      }
      else { std::cout << "L1RCT: hf out of range!  region = " << tower << std::endl; }
    }

  }
  
  input();

  return;

}
unsigned short L1RCT::ecalCompressedET ( int  crate,
int  card,
int  tower 
) [inline]

Definition at line 78 of file L1RCT.h.

References barrel.

Referenced by L1RCTSaveInput::analyze(), and L1RCTInputProducer::produce().

    {return barrel.at(crate).at(card).at(tower) / 2;}
unsigned short L1RCT::ecalFineGrainBit ( int  crate,
int  card,
int  tower 
) [inline]

Definition at line 80 of file L1RCT.h.

References barrel.

Referenced by L1RCTSaveInput::analyze(), and L1RCTInputProducer::produce().

    {return barrel.at(crate).at(card).at(tower) & 1;}
void L1RCT::fileInput ( const char *  filename)

Definition at line 85 of file L1RCT.cc.

References barrel, relval_parameters_module::energy, hf, i, input(), j, gen::k, and x.

                                         {            // added "const" also in .h
  unsigned short x;
  std::ifstream instream(filename);
  if(instream){
    for(int i = 0; i<18;i++){
      for(int j = 0; j<7; j++){
        for(int k = 0; k<64; k++){
          if(instream >> x){
            unsigned short bit = x/256;             // added J.Leonard Aug. 16 06
            unsigned short energy = x&255;          //
            unsigned short input = energy*2 + bit;  //
            barrel.at(i).at(j).at(k) = input;
          }
          else
            break;
        }
      }
      for(int j = 0; j<8; j++){
        if(instream >> x){
          hf.at(i).at(j) = x;
        }
        else
          break;
      }
    }
  }
  input();
}
L1CaloEmCollection L1RCT::getIsolatedEGObjects ( unsigned  crate)

Definition at line 332 of file L1RCT.cc.

References crates, L1RCTLookupTables::emRank(), relval_parameters_module::energy, i, and rctLookupTables_.

Referenced by L1RCTProducer::produce().

                                                            {
  std::vector<unsigned short> isoEmObjects = crates.at(crate).getIsolatedEGObjects();
  L1CaloEmCollection isoEmCands;
  for (uint16_t i = 0; i < 4; i++){
    unsigned rgn = ((isoEmObjects.at(i)) & 1);
    unsigned crd = (((isoEmObjects.at(i))/2) & 7);
    unsigned energy = ((isoEmObjects.at(i))/16);
    unsigned rank = rctLookupTables_->emRank(energy);
    L1CaloEmCand isoCand(rank, rgn, crd, crate, 1, i, 0);  // includes emcand index
    isoEmCands.push_back(isoCand);
  }
  return isoEmCands;
}
std::vector<unsigned short> L1RCT::getJetRegions ( unsigned  crate) [inline]

Definition at line 72 of file L1RCT.h.

References crates.

                                                         {
    return crates.at(crate).getJetRegions();
  }
L1CaloEmCollection L1RCT::getNonisolatedEGObjects ( unsigned  crate)

Definition at line 349 of file L1RCT.cc.

References crates, L1RCTLookupTables::emRank(), relval_parameters_module::energy, i, and rctLookupTables_.

Referenced by L1RCTProducer::produce().

                                                               {
  std::vector<unsigned short> nonIsoEmObjects = crates.at(crate).getNonisolatedEGObjects();
  L1CaloEmCollection nonIsoEmCands;
  for (uint16_t i = 0; i < 4; i++){
    unsigned rgn = ((nonIsoEmObjects.at(i)) & 1);
    unsigned crd = (((nonIsoEmObjects.at(i))/2) & 7);
    unsigned energy = ((nonIsoEmObjects.at(i))/16);
    unsigned rank = rctLookupTables_->emRank(energy);
    L1CaloEmCand nonIsoCand(rank, rgn, crd, crate, 0, i, 0);  // includes emcand index
    nonIsoEmCands.push_back(nonIsoCand);
  }
  return nonIsoEmCands;
}
vector< L1CaloRegion > L1RCT::getRegions ( unsigned  crate)

Definition at line 363 of file L1RCT.cc.

References crates, relval_parameters_module::energy, and metsig::tau.

Referenced by L1RCTProducer::produce().

                                                    {
  // barrel regions
  std::bitset<14> taus( (long) crates.at(crate).getTauBits());
  std::bitset<14> mips( (long) crates.at(crate).getMIPBits());
  std::bitset<14> quiets( (long) crates.at(crate).getQuietBits());
  std::bitset<14> overflows( (long) crates.at(crate).getOverFlowBits());
  std::vector<unsigned short> barrelEnergies = crates.at(crate).getBarrelRegions();
  std::vector<L1CaloRegion> regionCollection;
  for (unsigned card = 0; card < 7; card++){
    for (unsigned rgn = 0; rgn < 2; rgn++){
      bool tau = taus[card*2+rgn];
      bool mip = mips[card*2+rgn];
      bool quiet = quiets[card*2+rgn];
      bool overflow = overflows[card*2+rgn];
      unsigned barrelEnergy = barrelEnergies.at(card*2+rgn);
      L1CaloRegion region(barrelEnergy, overflow, tau, mip, quiet, crate, card, rgn);
      regionCollection.push_back(region);
    }
  }

  // hf regions
  std::vector<unsigned short> hfEnergies = crates.at(crate).getHFRegions();
  // fine grain bits -- still have to work out digi input
  std::vector<unsigned short> hfFineGrainBits = crates.at(crate).getHFFineGrainBits();
  for (unsigned hfRgn = 0; hfRgn<8; hfRgn++){  // region number, see diagram on paper.  make sure know how hf regions come in. 
    unsigned energy = hfEnergies.at(hfRgn);
    bool fineGrain = hfFineGrainBits.at(hfRgn);
    L1CaloRegion hfRegion(energy, fineGrain, crate, hfRgn);
    regionCollection.push_back(hfRegion);
  }
  return regionCollection;
}
unsigned short L1RCT::hcalCompressedET ( int  crate,
int  card,
int  tower 
) [inline]

Definition at line 82 of file L1RCT.h.

References barrel.

Referenced by L1RCTSaveInput::analyze(), and L1RCTInputProducer::produce().

    {return barrel.at(crate).at(card).at(tower+32) / 2;}
unsigned short L1RCT::hcalFineGrainBit ( int  crate,
int  card,
int  tower 
) [inline]

Definition at line 84 of file L1RCT.h.

References barrel.

Referenced by L1RCTSaveInput::analyze().

    {return barrel.at(crate).at(card).at(tower+32) & 1;}
unsigned short L1RCT::hfCompressedET ( int  crate,
int  tower 
) [inline]

Definition at line 86 of file L1RCT.h.

References hf.

Referenced by L1RCTInputProducer::produce().

    {return hf.at(crate).at(tower) / 2;}
unsigned short L1RCT::hfFineGrainBit ( int  crate,
int  tower 
) [inline]

Definition at line 88 of file L1RCT.h.

References hf.

    {return hf.at(crate).at(tower) & 1;}
void L1RCT::input ( )

Definition at line 66 of file L1RCT.cc.

References barrel, crates, hf, and i.

Referenced by digiInput(), fileInput(), and randomInput().

{
  for(int i = 0; i<18; i++){
    crates.at(i).input(barrel.at(i),hf.at(i));
  }
}
void L1RCT::input ( std::vector< std::vector< std::vector< unsigned short > > >  barrelIn,
std::vector< std::vector< unsigned short > >  hfIn 
)

Definition at line 73 of file L1RCT.cc.

References crates, and i.

{
  for(int i = 0; i<18; i++){
    crates.at(i).input(barrelIn.at(i),hfIn.at(i));
  }
}
void L1RCT::makeCrates ( ) [private]

Definition at line 48 of file L1RCT.cc.

References trackerHits::c, crates, i, and rctLookupTables_.

Referenced by L1RCT().

{
  for(int i = 0; i<18; i++){
    L1RCTCrate c(i, rctLookupTables_);
    crates.push_back(c);
  }
}
void L1RCT::print ( void  )

Definition at line 323 of file L1RCT.cc.

References gather_cfg::cout, crates, and i.

                 {
  for(int i = 0; i<18; i++){
    std::cout << "Crate " << i << std::endl;
    crates.at(i).print();
  } 
}
void L1RCT::printCrate ( int  i) [inline]

Definition at line 46 of file L1RCT.h.

References crates.

                        {
    crates.at(i).print();
  }
void L1RCT::printEIC ( int  i,
int  j 
) [inline]

Definition at line 61 of file L1RCT.h.

References crates.

                             {
    crates.at(i).printEIC(j);
  }
void L1RCT::printEICEdges ( int  i,
int  j 
) [inline]

Definition at line 64 of file L1RCT.h.

References crates.

                                  {
    crates.at(i).printEICEdges(j);
  }
void L1RCT::printJSC ( int  i) [inline]

Definition at line 49 of file L1RCT.h.

References crates.

                      {
    crates.at(i).printJSC();
  }
void L1RCT::printJSC ( ) [inline]

Definition at line 52 of file L1RCT.h.

References gather_cfg::cout, crates, and i.

                 {
    for(int i=0;i<18;i++){
      std::cout << "JSC for Crate " << i << std::endl;
      crates.at(i).printJSC();
    }
  }
void L1RCT::printRC ( int  i,
int  j 
) [inline]

Definition at line 58 of file L1RCT.h.

References crates.

                            {
    crates.at(i).printRC(j);
  }
void L1RCT::processEvent ( )

Definition at line 36 of file L1RCT.cc.

References crates, i, and shareNeighbors().

Referenced by L1RCTProducer::produce().

                        {
  for(int i=0; i<18;i++)
    crates.at(i).processReceiverCards();
  shareNeighbors();  
  for(int i=0; i<18;i++){
    crates.at(i).fillElectronIsolationCards();
    crates.at(i).processElectronIsolationCards();
    crates.at(i).fillJetSummaryCard();
    crates.at(i).processJetSummaryCard();
  }
}
void L1RCT::randomInput ( )

Definition at line 217 of file L1RCT.cc.

References barrel, hf, i, input(), j, and gen::k.

{
  for(int i = 0; i<18;i++){
    for(int j = 0; j<7;j++){
      for(int k = 0; k<64; k++){
        barrel.at(i).at(j).at(k) = rand()%511;
      }
    }
    for(int j = 0; j<8;j++){
      hf.at(i).at(j) = rand()%255;  // changed from 1023 (10 bits)
    }
  }
  input();
  return;
}
void L1RCT::shareNeighbors ( ) [private]

Definition at line 240 of file L1RCT.cc.

References crates, east, L1RCTNeighborMap::east(), empty, L1RCTRegion::giveEastEt(), L1RCTRegion::giveEastHE_FG(), L1RCTRegion::giveNEEt(), L1RCTRegion::giveNEHE_FG(), L1RCTRegion::giveNorthEt(), L1RCTRegion::giveNorthHE_FG(), L1RCTRegion::giveNWEt(), L1RCTRegion::giveNWHE_FG(), L1RCTRegion::giveSEEt(), L1RCTRegion::giveSEHE_FG(), L1RCTRegion::giveSouthEt(), L1RCTRegion::giveSouthHE_FG(), L1RCTRegion::giveSWEt(), L1RCTRegion::giveSWHE_FG(), L1RCTRegion::giveWestEt(), L1RCTRegion::giveWestHE_FG(), i, j, gen::k, Vispa::Plugins::EdmBrowser::EdmDataAccessor::ne(), L1RCTNeighborMap::ne(), neighborMap, north, L1RCTNeighborMap::north(), L1RCTNeighborMap::nw(), L1RCTNeighborMap::se(), L1RCTRegion::setEastEt(), L1RCTRegion::setEastHE_FG(), L1RCTRegion::setNEEt(), L1RCTRegion::setNEHE_FG(), L1RCTRegion::setNorthEt(), L1RCTRegion::setNorthHE_FG(), L1RCTRegion::setNWEt(), L1RCTRegion::setNWHE_FG(), L1RCTRegion::setSEEt(), L1RCTRegion::setSEHE_FG(), L1RCTRegion::setSouthEt(), L1RCTRegion::setSouthHE_FG(), L1RCTRegion::setSWEt(), L1RCTRegion::setSWHE_FG(), L1RCTRegion::setWestEt(), L1RCTRegion::setWestHE_FG(), south, L1RCTNeighborMap::south(), L1RCTNeighborMap::sw(), L1RCTNeighborMap::west(), and west.

Referenced by processEvent().

                          {
  L1RCTRegion *north;
  L1RCTRegion *south;
  L1RCTRegion *west;
  L1RCTRegion *east;
  L1RCTRegion *se;
  L1RCTRegion *sw;
  L1RCTRegion *nw;
  L1RCTRegion *ne;
  L1RCTRegion *primary;


  for(int i = 0; i < 18; i++){
    for(int j = 0; j < 7; j++){
      for(int k = 0; k < 2; k++){

        primary = crates.at(i).getReceiverCard(j)->getRegion(k);

        vector<int> northIndices = neighborMap.north(i,j,k);
        if(northIndices.at(0) != -1)
          north = crates.at(northIndices.at(0)).getReceiverCard(northIndices.at(1))->getRegion(northIndices.at(2));
        else north = &empty;

        vector<int> southIndices = neighborMap.south(i,j,k);
        if(southIndices.at(0) != -1)
          south = crates.at(southIndices.at(0)).getReceiverCard(southIndices.at(1))->getRegion(southIndices.at(2));
        else south = &empty;

        vector<int> westIndices = neighborMap.west(i,j,k);
        if(westIndices.at(0) != -1)
          west = crates.at(westIndices.at(0)).getReceiverCard(westIndices.at(1))->getRegion(westIndices.at(2));
        else west = &empty;

        vector<int> eastIndices = neighborMap.east(i,j,k);
        if(eastIndices.at(0) != -1)
          east = crates.at(eastIndices.at(0)).getReceiverCard(eastIndices.at(1))->getRegion(eastIndices.at(2));
        else east = &empty;

        vector<int> seIndices = neighborMap.se(i,j,k);
        if(seIndices.at(0) != -1)
          se = crates.at(seIndices.at(0)).getReceiverCard(seIndices.at(1))->getRegion(seIndices.at(2));
        else se = &empty;

        vector<int> swIndices = neighborMap.sw(i,j,k);
        if(swIndices.at(0) != -1)
          sw= crates.at(swIndices.at(0)).getReceiverCard(swIndices.at(1))->getRegion(swIndices.at(2));
        else sw = &empty;

        vector<int> neIndices = neighborMap.ne(i,j,k);
        if(neIndices.at(0) != -1)
          ne= crates.at(neIndices.at(0)).getReceiverCard(neIndices.at(1))->getRegion(neIndices.at(2));
        else ne = &empty;

        vector<int> nwIndices = neighborMap.nw(i,j,k);
        if(nwIndices.at(0) != -1)
          nw= crates.at(nwIndices.at(0)).getReceiverCard(nwIndices.at(1))->getRegion(nwIndices.at(2));
        else nw = &empty;

        primary->setNorthEt(north->giveNorthEt());
        primary->setNorthHE_FG(north->giveNorthHE_FG());
        primary->setSouthEt(south->giveSouthEt());
        primary->setSouthHE_FG(south->giveSouthHE_FG());
        primary->setEastEt(east->giveEastEt());
        primary->setEastHE_FG(east->giveEastHE_FG());
        primary->setWestEt(west->giveWestEt());
        primary->setWestHE_FG(west->giveWestHE_FG());
        primary->setSEEt(se->giveSEEt());
        primary->setSEHE_FG(se->giveSEHE_FG());
        primary->setSWEt(sw->giveSWEt());
        primary->setSWHE_FG(sw->giveSWHE_FG());
        primary->setNWEt(nw->giveNWEt());
        primary->setNWHE_FG(nw->giveNWHE_FG());
        primary->setNEEt(ne->giveNEEt());
        primary->setNEHE_FG(ne->giveNEHE_FG());



      }
    }
  }

}

Member Data Documentation

std::vector<std::vector<std::vector<unsigned short> > > L1RCT::barrel [private]
std::vector<L1RCTCrate> L1RCT::crates [private]

Definition at line 104 of file L1RCT.h.

Referenced by shareNeighbors().

std::vector<std::vector<unsigned short> > L1RCT::hf [private]

Definition at line 129 of file L1RCT.h.

Referenced by digiInput(), fileInput(), hfCompressedET(), hfFineGrainBit(), input(), and randomInput().

Definition at line 111 of file L1RCT.h.

Referenced by shareNeighbors().

Definition at line 95 of file L1RCT.h.

Referenced by digiInput(), getIsolatedEGObjects(), getNonisolatedEGObjects(), and makeCrates().