/data/refman/pasoursint/CMSSW_5_3_9_patch3/src/CalibFormats/SiPixelObjects/src/PixelTrimAllPixels.cc

Go to the documentation of this file.

//
// This class provide a base class for the
// pixel trim data for the pixel FEC configuration
// This is a pure interface (abstract class) that
// needs to have an implementation.
//
// Need to figure out what is 'VMEcommand' below! 
//
// All applications should just use this 
// interface and not care about the specific
// implementation.
//

#include <sstream>
#include <iostream>
#include <ios>
#include <assert.h>
#include <stdexcept>
#include "CalibFormats/SiPixelObjects/interface/PixelTrimAllPixels.h"
#include "CalibFormats/SiPixelObjects/interface/PixelTimeFormatter.h"
#include "CalibFormats/SiPixelObjects/interface/PixelBase64.h"

using namespace pos;

PixelTrimAllPixels::PixelTrimAllPixels( std::vector <std::vector<std::string> >& tableMat):
  PixelTrimBase("","","")
{
    std::string mthn = "]\t[PixelTrimAllPixels::PixelTrimAllPixels()]\t\t    " ; 
    std::stringstream currentRocName;
    std::map<std::string , int > colM;
    std::vector<std::string > colNames;
    colNames.push_back("CONFIG_KEY"  );
    colNames.push_back("KEY_TYPE"    );
    colNames.push_back("KEY_ALIAS"   );
    colNames.push_back("VERSION"     );
    colNames.push_back("KIND_OF_COND");
    colNames.push_back("ROC_NAME"    );
    colNames.push_back("TRIM_BITS"   );
 

    for(unsigned int c = 0 ; c < tableMat[0].size() ; c++)
      {
        for(unsigned int n=0; n<colNames.size(); n++)
          {
            if(tableMat[0][c] == colNames[n])
              {
                colM[colNames[n]] = c;
                break;
              }
          }
      }//end for
    for(unsigned int n=0; n<colNames.size(); n++)
      {
      if(colM.find(colNames[n]) == colM.end())
        {
          std::cerr << __LINE__ << mthn << "Couldn't find in the database the column with name " << colNames[n] << std::endl;
          assert(0);
        }
      }
    
    //unsigned char *bits ;        /// supose to be " unsigned  char bits[tableMat[1][colM["TRIM_BLOB"]].size()] ;  "
    //char c[2080];
    std::string bits;
        trimbits_.clear() ;
    for(unsigned int r = 1 ; r < tableMat.size() ; r++)    //Goes to every row of the Matrix
      {
        PixelROCName rocid( tableMat[r][colM["ROC_NAME"]] );
        PixelROCTrimBits tmp;     
        tmp.read(rocid, base64_decode(tableMat[r][colM["TRIM_BITS"]])) ;
        trimbits_.push_back(tmp);
      }//end for r 
    //std::cout<<trimbits_.size()<<std::endl;
} //end contructor with databasa table

 
PixelTrimAllPixels::PixelTrimAllPixels(std::string filename):
  PixelTrimBase("","",""){

    if (filename[filename.size()-1]=='t'){

      std::ifstream in(filename.c_str());
        
      if (!in.good())  throw std::runtime_error("Failed to open file "+filename);
      //        std::cout << "filename =" << filename << std::endl;

      std::string s1;
      in >> s1; 

      trimbits_.clear();


      while (in.good()){
          
        std::string s2;
        in>>s2;

        //          std::cout << "PixelTrimAllPixels::PixelTrimAllPixels read s1:"<<s1<< " s2:" << s2 << std::endl;

        assert( s1 == "ROC:" );

        PixelROCName rocid(s2);

        //std::cout << "PixelTrimAllPixels::PixelTrimAllPixels read rocid:"<<rocid<<std::endl;
            
        PixelROCTrimBits tmp;
      
        tmp.read(rocid, in);

        trimbits_.push_back(tmp);

        in >> s1;

      }

      in.close();

    }
    else{

      std::ifstream in(filename.c_str(),std::ios::binary);
      if (!in.good())  throw std::runtime_error("Failed to open file "+filename);

      char nchar;

      in.read(&nchar,1);

      std::string s1;

      //wrote these lines of code without ref. needs to be fixed
      for(int i=0;i< nchar; i++){
        char c;
        in >>c;
        s1.push_back(c);
      }

      //std::cout << "READ ROC name:"<<s1<<std::endl;

      trimbits_.clear();


      while (!in.eof()){

        //std::cout << "PixelTrimAllPixels::PixelTrimAllPixels read s1:"<<s1<<std::endl;

        PixelROCName rocid(s1);

        //std::cout << "PixelTrimAllPixels::PixelTrimAllPixels read rocid:"<<rocid<<std::endl;
            
        PixelROCTrimBits tmp;
      
        tmp.readBinary(rocid, in);

        trimbits_.push_back(tmp);


        in.read(&nchar,1);

        s1.clear();

        if (in.eof()) continue;

        //wrote these lines of code without ref. needs to be fixed
        for(int i=0;i< nchar; i++){
          char c;
          in >>c;
          s1.push_back(c);
        }


      }

      in.close();



    }

    //std::cout << "Read trimbits for "<<trimbits_.size()<<" ROCs"<<std::endl;

  }


//std::string PixelTrimAllPixels::getConfigCommand(PixelMaskBase& pixelMask){
//
//  std::string s;
//  return s;
//
//}

PixelROCTrimBits PixelTrimAllPixels::getTrimBits(int ROCId) const {

  return trimbits_[ROCId];

}

PixelROCTrimBits* PixelTrimAllPixels::getTrimBits(PixelROCName name){

  for(unsigned int i=0;i<trimbits_.size();i++){
    if (trimbits_[i].name()==name) return &(trimbits_[i]);
  }

  return 0;

}



void PixelTrimAllPixels::generateConfiguration(PixelFECConfigInterface* pixelFEC,
                                               PixelNameTranslation* trans,
                                               const PixelMaskBase& pixelMask) const{

  for(unsigned int i=0;i<trimbits_.size();i++){

    std::vector<unsigned char> trimAndMasks(4160);

    const PixelROCMaskBits& maskbits=pixelMask.getMaskBits(i);

    for (unsigned int col=0;col<52;col++){
      for (unsigned int row=0;row<80;row++){
        unsigned char tmp=trimbits_[i].trim(col,row);
        if (maskbits.mask(col,row)!=0) tmp|=0x80;
        trimAndMasks[col*80+row]=tmp;
      }
    }

    // the slow way, one pixel at a time
    //pixelFEC->setMaskAndTrimAll(*(trans->getHdwAddress(trimbits_[i].name())),trimAndMasks);
    // the fast way, a full roc in column mode (& block xfer)
    const PixelHdwAddress* theROC = trans->getHdwAddress(trimbits_[i].name());
    pixelFEC->roctrimload(theROC->mfec(),
                          theROC->mfecchannel(),
                          theROC->hubaddress(),
                          theROC->portaddress(),
                          theROC->rocid(),
                          trimAndMasks);
  }
}

void PixelTrimAllPixels::writeBinary(std::string filename) const{

  
  std::ofstream out(filename.c_str(),std::ios::binary);

  for(unsigned int i=0;i<trimbits_.size();i++){
    trimbits_[i].writeBinary(out);
  }


}


void PixelTrimAllPixels::writeASCII(std::string dir) const{

  if (dir!="") dir+="/";
  PixelModuleName module(trimbits_[0].name().rocname());
  std::string filename=dir+"ROC_Trims_module_"+module.modulename()+".dat";
  

  std::ofstream out(filename.c_str());

  for(unsigned int i=0;i<trimbits_.size();i++){
    trimbits_[i].writeASCII(out);
  }


}
//=============================================================================================
void PixelTrimAllPixels::writeXMLHeader(pos::PixelConfigKey key, 
                                        int version, 
                                        std::string path, 
                                        std::ofstream *outstream,
                                        std::ofstream *out1stream,
                                        std::ofstream *out2stream) const
{
  std::string mthn = "[PixelTrimAllPixels::writeXMLHeader()]\t\t\t    " ;
  std::stringstream maskFullPath ;

  maskFullPath << path << "/Pixel_RocTrims_" << PixelTimeFormatter::getmSecTime() << ".xml";
  std::cout << mthn << "Writing to: " << maskFullPath.str() << std::endl ;

  outstream->open(maskFullPath.str().c_str()) ;
  
  *outstream << "<?xml version='1.0' encoding='UTF-8' standalone='yes'?>"                                 << std::endl ;
  *outstream << "<ROOT xmlns:xsi='https://www.w3.org/2001/XMLSchema-instance'>"                            << std::endl ;
  *outstream << ""                                                                                        << std::endl ; 
  *outstream << " <HEADER>"                                                                               << std::endl ; 
  *outstream << "  <TYPE>"                                                                                << std::endl ; 
  *outstream << "   <EXTENSION_TABLE_NAME>ROC_TRIMS</EXTENSION_TABLE_NAME>"                               << std::endl ; 
  *outstream << "   <NAME>ROC Trim Bits</NAME>"                                                           << std::endl ; 
  *outstream << "  </TYPE>"                                                                               << std::endl ; 
  *outstream << "  <RUN>"                                                                                 << std::endl ; 
  *outstream << "   <RUN_TYPE>ROC Trim Bits</RUN_TYPE>"                                                   << std::endl ; 
  *outstream << "   <RUN_NUMBER>1</RUN_NUMBER>"                                                           << std::endl ; 
  *outstream << "   <RUN_BEGIN_TIMESTAMP>" << PixelTimeFormatter::getTime() << "</RUN_BEGIN_TIMESTAMP>"   << std::endl ; 
  *outstream << "   <LOCATION>CERN P5</LOCATION>"                                                         << std::endl ; 
  *outstream << "  </RUN>"                                                                                << std::endl ; 
  *outstream << " </HEADER>"                                                                              << std::endl ; 
  *outstream << ""                                                                                        << std::endl ; 
  *outstream << " <DATA_SET>"                                                                             << std::endl ;
  *outstream << ""                                                                                        << std::endl ;
  *outstream << "  <VERSION>"             << version      << "</VERSION>"                                 << std::endl ;
  *outstream << "  <COMMENT_DESCRIPTION>" << getComment() << "</COMMENT_DESCRIPTION>"                     << std::endl ;
  *outstream << "  <CREATED_BY_USER>"   << getAuthor()  << "</CREATED_BY_USER>"                   << std::endl ;
  *outstream << ""                                                                                        << std::endl ;
  *outstream << "  <PART>"                                                                                << std::endl ;
  *outstream << "   <NAME_LABEL>CMS-PIXEL-ROOT</NAME_LABEL>"                                              << std::endl ;      
  *outstream << "   <KIND_OF_PART>Detector ROOT</KIND_OF_PART>"                                           << std::endl ;         
  *outstream << "  </PART>"                                                                               << std::endl ;
  *outstream << ""                                                                                        << std::endl ;

}

//=============================================================================================
void PixelTrimAllPixels::writeXML( std::ofstream *outstream,
                                   std::ofstream *out1stream,
                                   std::ofstream *out2stream) const 
{
  std::string mthn = "[PixelTrimAllPixels::writeXML()]\t\t\t    " ;

  for(unsigned int i=0;i<trimbits_.size();i++){
      trimbits_[i].writeXML(outstream);
  }
}

//=============================================================================================
void PixelTrimAllPixels::writeXMLTrailer(std::ofstream *outstream,
                                         std::ofstream *out1stream,
                                         std::ofstream *out2stream ) const 
{
  std::string mthn = "[PixelTrimAllPixels::writeXMLTrailer()]\t\t\t    " ;
  
  *outstream << " </DATA_SET>"                                                                            << std::endl ;
  *outstream << "</ROOT>"                                                                                 << std::endl ;
  
  outstream->close() ;
  std::cout << mthn << "Data written "                                                                    << std::endl ;

}