L1GtVhdlWriterBitManager.cc

Go to the documentation of this file.

// this class header
#include "L1TriggerConfig/L1GtConfigProducers/interface/L1GtVhdlWriterBitManager.h"

#include "L1TriggerConfig/L1GtConfigProducers/interface/L1GtVhdlTemplateFile.h"

#include "CondFormats/L1TObjects/interface/L1GtCaloTemplate.h"
#include "CondFormats/L1TObjects/interface/L1GtMuonTemplate.h"
#include "CondFormats/L1TObjects/interface/L1GtEnergySumTemplate.h"

// system include files
#include <iostream>
#include <fstream>
#include <map>
#include <string>
#include <sstream>
#include <vector>
#include <iomanip>

L1GtVhdlWriterBitManager::L1GtVhdlWriterBitManager()
{
    hex2binMap_["0"]="0000";
    hex2binMap_["1"]="0001";
    hex2binMap_["2"]="0010";
    hex2binMap_["3"]="0011";
    hex2binMap_["4"]="0100";
    hex2binMap_["5"]="0101";
    hex2binMap_["6"]="0110";
    hex2binMap_["7"]="0111";
    hex2binMap_["8"]="1000";
    hex2binMap_["9"]="1001";

    hex2binMap_["A"]="1010";
    hex2binMap_["B"]="1011";
    hex2binMap_["C"]="1100";
    hex2binMap_["D"]="1101";
    hex2binMap_["E"]="1110";
    hex2binMap_["F"]="1111";

    hex2binMap_["a"]="1010";
    hex2binMap_["b"]="1011";
    hex2binMap_["c"]="1100";
    hex2binMap_["d"]="1101";
    hex2binMap_["e"]="1110";
    hex2binMap_["f"]="1111";
}


std::string L1GtVhdlWriterBitManager::readMapInverse(const std::map<std::string,std::string>& map,std::string value)
{
    std::map<std::string,std::string>::const_iterator iter = map.begin();
    while (iter!=map.end())
    {
        if ((*iter).second == value)
            return (*iter).first;
        iter++;
    }
    return "";
}


std::string L1GtVhdlWriterBitManager::hex2bin(std::string hexString)
{

    std::string temp;
    for (unsigned int i=0; i<hexString.length(); i++)
    {
        std::string str;
        str=hexString[i];

        temp+=hex2binMap_[str];
    }

    return temp;
}


std::string L1GtVhdlWriterBitManager::bin2hex(std::string binString)
{

    std::string temp;
    for (unsigned int i=1; i<=binString.length(); i++)
    {
        //std::cout<<i%4<<std::endl;
        if(i%4==0)
        {
            //std::cout<<"I'm here!"<<std::endl;
            std::string str;
            str = binString.substr(i-4,4);
            //std::cout<<str<<std::cout<<std::endl;
            temp+=readMapInverse(hex2binMap_,str);
        }
    }

    return temp;
}


std::string L1GtVhdlWriterBitManager::mirror(unsigned int offset, std::string hexString, bool hexOutput)
{
    std::string temp,binString;

    char digit;
    bool hexInput=false;

    // check weather input hex or binary
    for(unsigned int i = 0; i<hexString.length(); i++)
    {
        if (hexString[i]!='0' || hexString[i]!='1' )
        {
            hexInput=true;
            break;
        }
    }

    if (hexInput)
        binString = hex2bin(hexString);
    else binString=hexString;

    unsigned int i=0, len=0;
    len = binString.length();

    if(offset > len)
        return binString;

    for(i = 0; i < (len - offset)/2; i++)
    {
        digit = binString[i + offset];

        binString[i + offset] = binString[len - 1 - i];
        binString[len - 1 - i] = digit;

    }

    if (hexOutput)
        return bin2hex(binString);
    else return binString;

}


std::string L1GtVhdlWriterBitManager::capitalLetters(std::string hexString)
{

    unsigned int i = 0;
    while(i<hexString.length())
    {
        if(hexString[i] == 'a') hexString[i] = 'A';
        else if(hexString[i] == 'b') hexString[i] = 'B';
        else if(hexString[i] == 'c') hexString[i] = 'C';
        else if(hexString[i] == 'd') hexString[i] = 'D';
        else if(hexString[i] == 'e') hexString[i] = 'E';
        else if(hexString[i] == 'f') hexString[i] = 'F';
        i++;
    }

    return hexString;
}


std::string L1GtVhdlWriterBitManager::shiftLeft(std::string hexString)
{
    std::string binString = hex2bin(hexString);

    binString.erase(0,1);
    binString+="0";

    return bin2hex(binString);
}


std::string L1GtVhdlWriterBitManager::buildEtaMuon(const std::vector<L1GtMuonTemplate::ObjectParameter>* op, const unsigned int &num, const unsigned int &counter)
{
    std::ostringstream ossEta;

    ossEta
        <<counter
        <<" => X\"";

    for (unsigned int i =0; i<num; i++)
    {
        std::ostringstream ossEtaRange;
        ossEtaRange << std::hex<< std::setw(16)<<std::setfill('0')<<(*op).at(i).etaRange<<std::dec;

        /*
        ossEta
            <<mirror(32,ossEtaRange.str());
        */
        
        ossEta << ossEtaRange.str();
        
        if (num>0 && i!=num-1) ossEta <<"_";

    }

    ossEta<<"\",\n";
    return ossEta.str();

}


std::string L1GtVhdlWriterBitManager::buildEtaCalo(const std::vector<L1GtCaloTemplate::ObjectParameter>* op, const unsigned int &num, const unsigned int &counter)
{
    std::ostringstream ossEta;

    ossEta
        <<counter
        <<" => X\"";

    // loop over all relevant components of object parameters

    for (unsigned int i =0; i<num; i++)
    {
        std::ostringstream ossEtaRange;

        ossEtaRange << std::hex <<std::setw(4)<<std::setfill('0')<< (*op).at(i).etaRange<<std::dec;

        /*
        std::string tempstr=hex2bin(ossEtaRange.str());
        tempstr[0]='0';
        tempstr[15]='0';
        */
        
        //ossEta<<std::setw(4)<<std::setfill('0')<<mirror(8,bin2hex(tempstr));

        ossEta << ossEtaRange.str();
        
        if (num>0 && i!=num-1) ossEta <<"_";

    }

    ossEta<<"\",\n";
    return ossEta.str();
}


std::string L1GtVhdlWriterBitManager::buildPhiCalo(const std::vector<L1GtCaloTemplate::ObjectParameter>* op, const unsigned int &num, const unsigned int &counter)
{
    std::ostringstream ossPhi;

    ossPhi
        <<counter
        <<" => X\"";

    for (unsigned int i =0; i<num; i++)
    {
        ossPhi << std::hex<< std::setw(8)<<std::setfill('0')<<(*op).at(i).phiRange<<std::dec;

        if (num>0 && i!=num-1) ossPhi <<"_";
    }

    ossPhi<<"\",\n";
    return capitalLetters(ossPhi.str());
}


std::string L1GtVhdlWriterBitManager::buildPhiEnergySum(const std::vector<L1GtEnergySumTemplate::ObjectParameter>* op, const unsigned int &num, const unsigned int &counter)
{
    std::ostringstream ossPhi,count;

    if ((*op).at(0).phiRange0Word!=0)
    {
        ossPhi << "000000000000000000"<<std::hex<<(*op).at(0).phiRange1Word<<(*op).at(0).phiRange0Word<<std::dec;

        count<<counter;

        return (count.str()+" => X\""+capitalLetters(ossPhi.str())+"\",\n");

    } else
    return "";
}


std::string L1GtVhdlWriterBitManager::buildPhiMuon(const std::vector<L1GtMuonTemplate::ObjectParameter>* op, const unsigned int &num, const unsigned int &counter,bool high)
{
    std::ostringstream ossPhi;

    ossPhi
        <<counter
        <<" => X\"";

    for (unsigned int i =0; i<num; i++)
    {
        if (high)
            ossPhi << std::hex<< std::setw(2)<<std::setfill('0')<<(*op).at(i).phiHigh<<std::dec;
        else
            ossPhi << std::hex<< std::setw(2)<<std::setfill('0')<<(*op).at(i).phiLow<<std::dec;

        if (num>0 && i!=num-1) ossPhi <<"_";
    }

    ossPhi<<"\",\n";
    return capitalLetters(ossPhi.str());

}


std::string L1GtVhdlWriterBitManager::buildDeltaEtaCalo(const L1GtCaloTemplate::CorrelationParameter* &cp,const unsigned int &counter)
{
    std::ostringstream deltaEtaRange, dEta,res;

    //std::cout<<"b:"<<.hex2bin("00383800")<<std::endl;

    deltaEtaRange << std::hex<< std::setw(4)<<std::setfill('0')<<(*cp).deltaEtaRange<<std::dec;

    // mirror deltaEtaRang and shift the mirrored value to the left by one bit;
    // add the original value to the calculated value
    dEta<<hex2bin(shiftLeft(mirror(0,deltaEtaRange.str())))<<hex2bin(deltaEtaRange.str());

    std::string result = capitalLetters(bin2hex(dEta.str()));

    res<<counter<<" => X\""<<result<<"\",\n";

    return res.str();

}


std::string L1GtVhdlWriterBitManager::buildDeltaEtaMuon(const L1GtMuonTemplate::CorrelationParameter* &cp,const unsigned int &counter)
{
    std::ostringstream deltaEtaRange, dEta;
    deltaEtaRange << std::hex<< std::setw(16)<<std::setfill('0')<<(*cp).deltaEtaRange<<std::dec;

    // mirror deltaEtaRang and shift the mirrored value to the left by one bit;
    // add the original value to the calculated value

    std::string result = capitalLetters((shiftLeft(mirror(0,deltaEtaRange.str()))+deltaEtaRange.str()));

    dEta<<counter<<" => X\""<<result<<"\",\n";

    return dEta.str();

}


std::string L1GtVhdlWriterBitManager::buildDeltaPhiCalo(const L1GtCaloTemplate::CorrelationParameter* &cp,const unsigned int &counter)
{
    std::ostringstream dPhi,deltaPhiRange, result;
    //std::cout<<.hex2bin("03E0000000000F81")<<std::endl;
    //std::cout<<.hex2bin("0080000000000200")<<std::endl;

    deltaPhiRange << std::hex<<std::setw(3)<<std::setfill('0')<<(*cp).deltaPhiRange<<std::dec;

    std::string binString = hex2bin(deltaPhiRange.str());

    //std::cout <<"========================" <<std::endl;

    std::string help2 = binString.substr(2,binString.length()-2);
    std::string help1 = help2.substr(1);
    help1 = mirror(0,bin2hex(help1),false);

    // here delta phi is built
    result<<help1;

    // might be wrong - has to be tested with more reference values!
    result<<help2.substr(0,8);
    result<<"0";

    result<<"00000000000000000000000000000";

    result<<help1;

    result<<binString.substr(2,binString.length()-2);

    dPhi<<counter<<" => X\""<<bin2hex(result.str())<<"\",\n";

    //std::cout<<result<<std::endl;

    /*
     * Code from old GTS:
      bm_dphi = bm_dphi.Get(2);
      BitMngr help1, help2 = bm_dphi;
      help2.SetAt(9, '\0');
      help1 = help2.Get(1);
      help1.Mirror();
      BitMngr nuller;
      nuller.InitBin("00 00 00 00 00 00 00 00 00 00 00 00 00 00 0");
      BitMngr result;
      result = result + help1 + help2 + nuller + help1 + bm_dphi;
    dphi = result.GetHex();
    */

    return capitalLetters(dPhi.str());
}


std::string L1GtVhdlWriterBitManager::buildDeltaPhiMuon(const L1GtMuonTemplate::CorrelationParameter* &cp,const unsigned int &counter)
{
    std::ostringstream dPhi,deltaPhiRange0,deltaPhiRange1,temp,result ;
    std::string tempstr;

    deltaPhiRange0 << std::hex<< std::setw(16)<<std::setfill('0')<<(*cp).deltaPhiRange0Word<<std::dec;
    //deltaPhiRange1  /*<< std::hex<< std::setw(3)<<std::setfill('0')*/<<(*cp).deltaPhiRange1Word<<std::dec;

    //mirror deltaPhiRange, shift the mirrored value left and convert it to binary format
    temp<<hex2bin(shiftLeft(mirror(0,deltaPhiRange0.str())));
    tempstr=temp.str();
    // set last bit of tempstr to one;
    tempstr[tempstr.length()-1]='1';

    // build delta eta as stringstreamtau
    result
        <<tempstr
    // insert 16 ones
        <<hex2bin("FFFF")
        <<hex2bin(deltaPhiRange0.str());

    dPhi<<counter<<" => X\""<<bin2hex(result.str())<<"\",\n";

    return dPhi.str();

}