/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/SimCalorimetry/EcalTrigPrimAlgos/src/EcalFenixStripFormatEE.cc

Go to the documentation of this file.

#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixStripFormatEE.h>
#include <CondFormats/EcalObjects/interface/EcalTPGSlidingWindow.h>
#include <CondFormats/EcalObjects/interface/EcalTPGStripStatus.h>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>


//-----------------------------------------------------------------------------------------
EcalFenixStripFormatEE::EcalFenixStripFormatEE() 
  : shift_(0)
{
}
//------------------------------------------------------------------------------------------

EcalFenixStripFormatEE::~EcalFenixStripFormatEE() {
}

//-----------------------------------------------------------------------------------------


int EcalFenixStripFormatEE::setInput(int input, int inputPeak, int fgvb ) {
  inputPeak_=inputPeak;
  input_=input;
  fgvb_=fgvb;
  return 0;
}  
//-----------------------------------------------------------------------------------------
  
int EcalFenixStripFormatEE::process(){
  // Bad strip - zero everything
  if(stripStatus_ != 0) return 0;

  // Peak not found - only return fgvb
  if(inputPeak_==0) return ((fgvb_ & 0x1) << 12);
    
  int output=input_>>shift_;

  //barrel saturates at 12 bits, endcap at 10!
  // Pascal: finally no,endcap has 12 bits as in EB (bug in FENIX!!!!)
  if(output>0XFFF) output=0XFFF; 
  output=output|((fgvb_ & 0x1) << 12); //Pascal (was 10)

  return output;    
} 
//------------------------------------------------------------------------------------------

void EcalFenixStripFormatEE::process(std::vector<int> &fgvbout,std::vector<int> &peakout,
                                     std::vector<int> &filtout, std::vector<int> &output){
  if  (peakout.size()!=filtout.size()){
    edm::LogWarning("EcalTPG")<<" problem in EcalFenixStripFormatEE: peak_out and filt_out don't have the same size";
    std::cout<<" Size peak_out"<< peakout.size()<<", size filt_out:"<<filtout.size()<<std::flush<<std::endl;
  }
  for  (unsigned int i =0;i<filtout.size();i++){
    setInput(filtout[i],peakout[i],fgvbout[i]);
    output[i]=process();
  }
  return;
}
//-----------------------------------------------------------------------------------------

void EcalFenixStripFormatEE::setParameters(uint32_t id,const EcalTPGSlidingWindow*& slWin,const EcalTPGStripStatus *stripStatus){

  const EcalTPGSlidingWindowMap &slwinmap = slWin -> getMap();
  EcalTPGSlidingWindowMapIterator it=slwinmap.find(id);
  if (it!=slwinmap.end()) shift_=(*it).second;
  else edm::LogWarning("EcalTPG")<<" could not find EcalTPGSlidingWindowMap entry for "<<id;

  const EcalTPGStripStatusMap &statusMap = stripStatus->getMap();
  EcalTPGStripStatusMapIterator sit = statusMap.find(id);
  if(sit != statusMap.end())
  {
    stripStatus_ = (*sit).second;
  }
  else
  {
    stripStatus_ = 0; // Assume strip OK
  }
}
//-----------------------------------------------------------------------------------------