DTConfigBti.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: DTConfigBti
//
//   Description: Configurable parameters and constants for Level1 Mu DT Trigger - BTI chip
//
//
//   Author List:
//   S.Vanini
//
//   Modifications:
//   April,10th : set BTI parameters from string
//-----------------------------------------------------------------------

//-----------------------
// This Class's Header --
//-----------------------
#include "L1TriggerConfig/DTTPGConfig/interface/DTConfigBti.h"
#include "FWCore/Utilities/interface/Exception.h"

//---------------
// C++ Headers --
//---------------
#include <cstdlib>
#include <string>  
#include <iostream>
#include <fstream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <iomanip>
               
//-------------------------------
// Collaborating Class Headers --
//-------------------------------

//----------------
// Constructors --
//----------------
DTConfigBti::DTConfigBti(const edm::ParameterSet& ps) { 

  setDefaults(ps);
  
}

DTConfigBti::DTConfigBti(int debugBTI, unsigned short int * buffer) {

        m_debug = debugBTI;

    // check if this is a BTI configuration string
    if (buffer[2]!=0x54){
        throw cms::Exception("DTTPG") << "===> ConfigBti constructor : not a BTI string!" << std::endl;
    }

    // decode
    unsigned short int memory_bti[31];
    for(int i=0;i<31;i++){
                memory_bti[i] = buffer[i+5];
                //std::cout << hex << memory_bti[i] << "  ";
        }
    int wmask[9];
        int st43 = memory_bti[0] & 0x3f;
    wmask[6] = memory_bti[1] & 0x01;
        wmask[7] = (memory_bti[1] >> 1 )& 0x01;
        wmask[8] = (memory_bti[1] >> 2 )& 0x01;
    int re43 = (memory_bti[1] >> 4 )& 0x03;
        wmask[0] = memory_bti[2] & 0x01;
        wmask[1] = (memory_bti[2] >> 1) & 0x01;
    wmask[2] = (memory_bti[2] >> 2 )& 0x01;
    wmask[3] = (memory_bti[2] >> 3 )& 0x01;
    wmask[4] = (memory_bti[2] >> 4 )& 0x01;
    wmask[5] = (memory_bti[2] >> 5 )& 0x01;
        int dead = memory_bti[3] & 0x3F;
        int LH =  memory_bti[4] & 0x3F;
        int LL =  memory_bti[5] & 0x3F;
        int CH =  memory_bti[6] & 0x3F;
        int CL =  memory_bti[7] & 0x3F;
        int RH =  memory_bti[8] & 0x3F;
        int RL =  memory_bti[9] & 0x3F;
        int tston = ( memory_bti[10] & 0x20 ) != 0 ;
        int test = ( memory_bti[10] & 0x10 ) != 0 ;
        int ten = ( memory_bti[10] & 0x8 ) != 0 ;
        int xon = ( memory_bti[10] & 0x2 ) != 0 ;
        int ron = ( memory_bti[10] & 0x1 ) != 0 ;
        int set = ( memory_bti[11] & 0x38 ) >> 3 ;
        int lts = ( memory_bti[11] & 0x6 ) >> 1 ;
        int ac1 = ( memory_bti[12] & 0x30 ) >> 4 ;
        int ac2 = ( memory_bti[12] & 0xc ) >> 2 ;
        int acl = ( memory_bti[12] & 0x3 ) ;
        int ach = ( memory_bti[13] & 0x30 ) >> 4 ;

    int pmask[32];
    for(int ir=0; ir<6; ir++)
    {
            pmask[0+6*ir] = memory_bti[19-ir] & 0x01;
            pmask[1+6*ir] = (memory_bti[19-ir] >> 1 )& 0x01;    
        if(ir!=5)
        {
            pmask[2+6*ir] = (memory_bti[19-ir] >> 2 )& 0x01;        
            pmask[3+6*ir] = (memory_bti[19-ir] >> 3 )& 0x01;
            pmask[4+6*ir] = (memory_bti[19-ir] >> 4 )& 0x01;    
            pmask[5+6*ir] = (memory_bti[19-ir] >> 5 )& 0x01;
        }
    }   
    
    // dump
        if(debug()==1){
               std::cout << std::dec << "st43=" << st43
                << " re43=" << re43
                << " dead=" << dead
                << " LH=" << LH
                << " LL=" << LL
                << " CH=" << CH
                << " CL=" << CL
                << " RH=" << RH
                << " RL=" << RL
                << " tston=" << tston
                << " test=" << test
                << " ten=" << ten
                << " xon=" << xon
                << " ron=" << ron
                << " set=" << set
                << " lts=" << lts
                << " ac1=" << ac1
                << " ac2=" << ac2
                << " acl=" << acl
                << " ach=" << ach
                << std::endl;
        std::cout << std::dec << " wire masks= ";
        for(int iw=0; iw<9; iw++)
            std::cout << wmask[iw] << " ";
        std::cout << std::dec << "\n pattern masks= ";
        for(int ip=0; ip<32; ip++)
            std::cout << pmask[ip] << " ";
        std::cout << std::endl;
    }
    
    // set parameters
    // default for KCut and KAccTheta
    setKCut(64);
    setKAccTheta(1);

    for(int ip=0; ip<32; ip++)
        setPTMSflag(pmask[ip],ip);

    for(int iw=0; iw<9; iw++)
        setWENflag(wmask[iw],iw+1);

    setST43(st43);
    setRE43(re43);
    setDEADpar(dead);
    setLL(LL);
    setLH(LH);
    setCL(CL);
    setCH(CH);
    setRL(RL);
    setRH(RH);
    setXON(xon);
    setRONflag(ron);
    setSET(set);
    setLTS(lts);
    setAccPattAC1(ac1);
    setAccPattAC2(ac2);
    setAccPattACH(ach);
    setAccPattACL(acl);
} 

//--------------
// Destructor --
//--------------
DTConfigBti::~DTConfigBti() {}

//--------------
// Operations --
//--------------

void
DTConfigBti::setDefaults(const edm::ParameterSet& ps) {

  // Debug flag 
  m_debug  = ps.getUntrackedParameter<int>("Debug");

  // Max K param accepted  
  m_kcut = ps.getParameter<int>("KMAX");
 
  // BTI angular acceptance in theta view
  m_kacctheta = ps.getParameter<int>("KACCTHETA");

  // Time indep. K equation suppression (XON) 
  m_xon = ps.getParameter<bool>("XON");
  // LTS and SET for low trigger suppression 
  m_lts = ps.getParameter<int>("LTS");
  m_set = ps.getParameter<int>("SET");
  // pattern acceptance AC1, AC2, ACH, ACL
  m_ac1 = ps.getParameter<int>("AC1");
  m_ac2 = ps.getParameter<int>("AC2");
  m_ach = ps.getParameter<int>("ACH");
  m_acl = ps.getParameter<int>("ACL");
  // redundant patterns flag RON
  m_ron = ps.getParameter<bool>("RON");

  // pattern masks
  for(int p=0; p<32; p++)
  {
        std::string label = "PTMS";
    char patt0 = (p/10)+'0';
    char patt1 = (p%10)+'0';
    if ( patt0 != '0' )
            label = label + patt0;
    label = label + patt1;
    
    m_pattmask.set(p,ps.getParameter<int>(label));
  }

  // wire masks
  for(int w=0; w<9; w++)
  {
        std::string label = "WEN";
    char wname = w+'0';
    label = label + wname;
    
    m_wiremask.set(w,ps.getParameter<int>(label));
  }

  // angular window limits for traco
  m_ll = ps.getParameter<int>("LL");
  m_lh = ps.getParameter<int>("LH");
  m_cl = ps.getParameter<int>("CL");
  m_ch = ps.getParameter<int>("CH");
  m_rl = ps.getParameter<int>("RL");
  m_rh = ps.getParameter<int>("RH");
  // drift velocity parameter 4ST3
  m_4st3 = ps.getParameter<int>("ST43");
  // drift velocity parameter 4RE3
  m_4re3 = ps.getParameter<int>("RE43");
  // DEAD parameter
  m_dead = ps.getParameter<int>("DEAD");
}

void 
DTConfigBti::print() const {
  std::cout << "******************************************************************************" << std::endl;
  std::cout << "*              DTTrigger configuration : BTI chips                                 *" << std::endl;
  std::cout << "******************************************************************************" << std::endl;
  std::cout << "*                                                                            *" << std::endl;
  std::cout << "Debug flag : " <<  debug() << std::endl;
  std::cout << "Max K param accepted : " << KCut() << std::endl; 
  std::cout << "BTI angular acceptance in theta view : " << KAccTheta() << std::endl;
  std::cout << "Time indep. K equation suppression (XON) : " << XON() << std::endl;
  std::cout << "LTS for low trigger suppression : " << LTS() << std::endl;
  std::cout << "SET for low trigger suppression : " << SET() << std::endl;
  std::cout << "pattern acceptance AC1, AC2, ACH, ACL : " << 
    AccPattAC1() << ", " << AccPattAC2() << " , " << AccPattACH() << ", " << AccPattACL() << std::endl;
  std::cout << "redundant patterns flag RON : " << RONflag() << std::endl;
  std::cout << "pattern masks : "; 
  for(int p=0; p<32; p++)
    std::cout << PTMSflag(p) << " ";
  std::cout << std::endl;
 
  std::cout << "wire masks : "; 
  for(int w=1; w<=9; w++)
    std::cout << WENflag(w) << " ";
  std::cout << std::endl;

  std::cout << "angular window limits for traco : " << LL() << ", " << LH() << ", " 
    << CL() << ", " << CH() << ", " << RL() << ", " << RH() << std::endl;
  std::cout << "drift velocity parameter 4ST3 : " << ST43() << std::endl;
  std::cout << "drift velocity parameter 4RE3 : " << RE43() << std::endl;  
  std::cout << "DEAD parameter : " << DEADpar() << std::endl;

  std::cout << "******************************************************************************" << std::endl;

}