Linearisation for Fenix strip input: 16 bits corresponding to input EBDataFrame output: 18 bits. More...

#include <EcalFenixLinearizer.h>

Public Member Functions
	EcalFenixLinearizer (bool famos)

	EcalFenixLinearizer (bool famos)

template<class T >
void	process (const T &, std::vector< int > &)

template<class T >
void	process (const T &, std::vector< int > &)

void	setParameters (uint32_t raw, const EcalTPGPedestals ecaltpPed, const EcalTPGLinearizationConst ecaltpLin, const EcalTPGCrystalStatus *ecaltpBadX)

void	setParameters (uint32_t raw, const EcalTPGPedestals ecaltpPed, const EcalTPGLinearizationConst ecaltpLin, const EcalTPGCrystalStatus *ecaltpBadX)

virtual	~EcalFenixLinearizer ()

virtual	~EcalFenixLinearizer ()

Private Member Functions
int	process ()

int	process ()

int	setInput (const EcalMGPASample &RawSam)

int	setInput (const EcalMGPASample &RawSam)

Private Attributes
const EcalTPGCrystalStatusCode *	badXStatus_

int	base_

bool	famos_

int	gainID_

bool	init_

const EcalTPGLinearizationConstant *	linConsts_

int	mult_

const EcalTPGPedestal *	peds_

int	shift_

int	strip_

int	uncorrectedSample_

std::vector< const EcalTPGCrystalStatusCode * >	vectorbadXStatus_

Detailed Description

Linearisation for Fenix strip input: 16 bits corresponding to input EBDataFrame output: 18 bits.

Definition at line 19 of file EcalFenixLinearizer.h.

Constructor & Destructor Documentation

◆ EcalFenixLinearizer() [1/2]

EcalFenixLinearizer::EcalFenixLinearizer ( bool famos )

Definition at line 9 of file EcalFenixLinearizer.cc.

10 : famos_(famos), init_(false), linConsts_(nullptr), peds_(nullptr), badXStatus_(nullptr) {}

◆ ~EcalFenixLinearizer() [1/2]

EcalFenixLinearizer::~EcalFenixLinearizer ( )

virtual

Definition at line 12 of file EcalFenixLinearizer.cc.

                                           {
   if (init_) {
     for (int i = 0; i < (int)vectorbadXStatus_.size(); i++) {
       delete vectorbadXStatus_[i];
     }
   }
 }

References mps_fire::i, init_, createfilelist::int, and vectorbadXStatus_.

◆ EcalFenixLinearizer() [2/2]

EcalFenixLinearizer::EcalFenixLinearizer ( bool famos )

◆ ~EcalFenixLinearizer() [2/2]

virtual EcalFenixLinearizer::~EcalFenixLinearizer ( )

virtual

Member Function Documentation

◆ process() [1/4]

int EcalFenixLinearizer::process ( )

private

Definition at line 51 of file EcalFenixLinearizer.cc.

                                  {
   int output = (uncorrectedSample_ - base_);  //Substract base
   //std::cout << " EcalFenixLinearizer::process() output non bit shifted " << output << std::endl;
   if (famos_ || output < 0)
     return 0;
  
   if (output < 0)
     return shift_ << 12;                      // FENIX bug(!)
   output = (output * mult_) >> (shift_ + 2);  //Apply multiplicative factor
   //std::cout << " EcalFenixLinearizer::process() output 2nd step " << output << std::endl;
   if (output > 0X3FFFF)
     output = 0X3FFFF;  //Saturation if too high
   //std::cout << " EcalFenixLinearizer::process() output 3rd step " << output << std::endl;
   return output;
 }

References base_, famos_, mult_, convertSQLitetoXML_cfg::output, shift_, and uncorrectedSample_.

Referenced by process(), EcalFenixStrip::process_part1(), and EcalEBTrigPrimTestAlgo::run().

◆ process() [2/4]

int EcalFenixLinearizer::process ( )

private

◆ process() [3/4]

template<class T >

void EcalFenixLinearizer::process	(	const T &	,
		std::vector< int > &
	)

◆ process() [4/4]

template<class T >

void EcalFenixLinearizer::process	(	const T &	df,
		std::vector< int > &	output_percry
	)

Definition at line 52 of file EcalFenixLinearizer.h.

                                                                             {
   //We know a tower numbering is:
   // S1 S2 S3 S4 S5
   //
   // 4  5  14 15 24
   // 3  6  13 16 23
   // 2  7  12 17 22
   // 1  8  11 18 21
   // 0  9  10 19 20
  
   //std::cout << " EcalFenixLinearizer::process(const  .. DataFrame " << std::endl;
   for (int i = 0; i < df.size(); i++) {
     //std::cout <<  df[i] << " ";
   }
  
   for (int i = 0; i < df.size(); i++) {
     setInput(df[i]);
     output_percry[i] = process();
   }
  
   //std::cout << " EcalFenixLinearizer::process(const  .. Final output " << std::endl;
   for (int i = 0; i < df.size(); i++) {
     //std::cout << " output_percry " << output_percry[i]<< " ";
   }
  
   return;
 }

References mps_fire::i, process(), and setInput().

◆ setInput() [1/2]

int EcalFenixLinearizer::setInput ( const EcalMGPASample & RawSam )

private

◆ setInput() [2/2]

int EcalFenixLinearizer::setInput ( const EcalMGPASample & RawSam )

private

Definition at line 67 of file EcalFenixLinearizer.cc.

                                                               {
   //std::cout << "  EcalFenixLinearizer::setInput RawSam.raw() " << RawSam.raw() << std::endl;
   if (RawSam.raw() > 0X3FFF) {
     LogDebug("EcalTPG") << "ERROR IN INPUT SAMPLE OF FENIX LINEARIZER";
     return -1;
   }
  
   uncorrectedSample_ = RawSam.adc();  //uncorrectedSample_ is coded in the 12 LSB
   gainID_ = RawSam.gainId();          //uncorrectedSample_ is coded in the 2 next bits!
  
   if (gainID_ == 0) {
     base_ = 0;
     shift_ = 0;
     mult_ = 0xFF;
     if ((linConsts_->mult_x12 == 0) && (linConsts_->mult_x6 == 0) && (linConsts_->mult_x1 == 0)) {
       mult_ = 0;  // Implemented in CCSSupervisor to
                   // reject overflow cases in rejected channels
     }
   } else if (gainID_ == 1) {
     base_ = peds_->mean_x12;
     shift_ = linConsts_->shift_x12;
  
     // take into account the badX
     // badXStatus_ == 0 if the crystal works
     // badXStatus_ !=0 some problem with the crystal
     if (badXStatus_->getStatusCode() != 0) {
       mult_ = 0;
     } else {
       mult_ = linConsts_->mult_x12;
     }
   } else if (gainID_ == 2) {
     base_ = peds_->mean_x6;
     shift_ = linConsts_->shift_x6;
  
     // take into account the badX
     // check if the badX has a status code=0 or 1
     if (badXStatus_->getStatusCode() != 0) {
       mult_ = 0;
     } else {
       mult_ = linConsts_->mult_x6;
     }
   } else if (gainID_ == 3) {
     base_ = peds_->mean_x1;
     shift_ = linConsts_->shift_x1;
  
     // take into account the badX
     // check if the badX has a status code=0 or 1
     if (badXStatus_->getStatusCode() != 0) {
       mult_ = 0;
     } else {
       mult_ = linConsts_->mult_x1;
     }
   }
  
   //std::cout << "  EcalFenixLinearizer::setInput   uncorrectedSample_ " << RawSam.adc() << " gainID " << gainID_ << " baseline " << base_ << std::endl;
  
   if (famos_)
     base_ = 200;  //FIXME by preparing a correct TPG.txt for Famos
  
   return 1;
 }

References EcalMGPASample::adc(), badXStatus_, base_, famos_, EcalMGPASample::gainId(), gainID_, EcalTPGCrystalStatusCode::getStatusCode(), linConsts_, LogDebug, EcalTPGPedestal::mean_x1, EcalTPGPedestal::mean_x12, EcalTPGPedestal::mean_x6, mult_, EcalTPGLinearizationConstant::mult_x1, EcalTPGLinearizationConstant::mult_x12, EcalTPGLinearizationConstant::mult_x6, peds_, EcalMGPASample::raw(), shift_, EcalTPGLinearizationConstant::shift_x1, EcalTPGLinearizationConstant::shift_x12, EcalTPGLinearizationConstant::shift_x6, and uncorrectedSample_.

Referenced by process().

◆ setParameters() [1/2]

void EcalFenixLinearizer::setParameters	(	uint32_t	raw,
		const EcalTPGPedestals *	ecaltpPed,
		const EcalTPGLinearizationConst *	ecaltpLin,
		const EcalTPGCrystalStatus *	ecaltpBadX
	)

◆ setParameters() [2/2]

void EcalFenixLinearizer::setParameters	(	uint32_t	raw,
		const EcalTPGPedestals *	ecaltpPed,
		const EcalTPGLinearizationConst *	ecaltpLin,
		const EcalTPGCrystalStatus *	ecaltpBadX
	)

Definition at line 20 of file EcalFenixLinearizer.cc.

                                                                                 {
   const EcalTPGLinearizationConstMap &linMap = ecaltpLin->getMap();
   EcalTPGLinearizationConstMapIterator it = linMap.find(raw);
  
   if (it != linMap.end()) {
     linConsts_ = &(*it);
   } else
     std::cout << " could not find EcalTPGLinearizationConstMap entry for " << raw << std::endl;
  
   EcalTPGPedestalsMapIterator itped = ecaltpPed->find(raw);
   if (itped != ecaltpPed->end())
     peds_ = &(*itped);
   else
     std::cout << " could not find EcalTPGPedestalsMap entry for " << raw << std::endl;
  
   const EcalTPGCrystalStatusMap &badXMap = ecaltpBadX->getMap();
   EcalTPGCrystalStatusMapIterator itbadX = badXMap.find(raw);
  
   if (itbadX != badXMap.end()) {
     badXStatus_ = &(*itbadX);
   } else {
     edm::LogWarning("EcalTPG") << " could not find EcalTPGCrystalStatusMap entry for " << raw;
     badXStatus_ = new EcalTPGCrystalStatusCode();
     vectorbadXStatus_.push_back(&(*badXStatus_));
     init_ = true;
   }
 }

References badXStatus_, gather_cfg::cout, EcalCondObjectContainer< T >::end(), EcalCondObjectContainer< T >::find(), EcalCondObjectContainer< T >::getMap(), init_, linConsts_, peds_, and vectorbadXStatus_.

Referenced by EcalFenixStrip::process_part1(), and EcalEBTrigPrimTestAlgo::run().

Member Data Documentation

◆ badXStatus_

const EcalTPGCrystalStatusCode * EcalFenixLinearizer::badXStatus_

private

Definition at line 32 of file EcalFenixLinearizer.h.

Referenced by setInput(), and setParameters().

◆ base_

int EcalFenixLinearizer::base_

private

Definition at line 24 of file EcalFenixLinearizer.h.

Referenced by process(), and setInput().

◆ famos_

bool EcalFenixLinearizer::famos_

private

Definition at line 21 of file EcalFenixLinearizer.h.

Referenced by process(), and setInput().

◆ gainID_

int EcalFenixLinearizer::gainID_

private

Definition at line 23 of file EcalFenixLinearizer.h.

Referenced by setInput().

◆ init_

bool EcalFenixLinearizer::init_

private

Definition at line 28 of file EcalFenixLinearizer.h.

Referenced by setParameters(), and ~EcalFenixLinearizer().

◆ linConsts_

const EcalTPGLinearizationConstant * EcalFenixLinearizer::linConsts_

private

Definition at line 30 of file EcalFenixLinearizer.h.

Referenced by setInput(), and setParameters().

◆ mult_

int EcalFenixLinearizer::mult_

private

Definition at line 25 of file EcalFenixLinearizer.h.

Referenced by process(), and setInput().

◆ peds_

const EcalTPGPedestal * EcalFenixLinearizer::peds_

private

Definition at line 31 of file EcalFenixLinearizer.h.

Referenced by setInput(), and setParameters().

◆ shift_

int EcalFenixLinearizer::shift_

private

Definition at line 26 of file EcalFenixLinearizer.h.

Referenced by process(), and setInput().

◆ strip_

int EcalFenixLinearizer::strip_

private

Definition at line 27 of file EcalFenixLinearizer.h.

◆ uncorrectedSample_

int EcalFenixLinearizer::uncorrectedSample_

private

Definition at line 22 of file EcalFenixLinearizer.h.

Referenced by process(), and setInput().

◆ vectorbadXStatus_

std::vector< const EcalTPGCrystalStatusCode * > EcalFenixLinearizer::vectorbadXStatus_

private

Definition at line 34 of file EcalFenixLinearizer.h.

Referenced by setParameters(), and ~EcalFenixLinearizer().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ EcalFenixLinearizer() [1/2]

◆ ~EcalFenixLinearizer() [1/2]

◆ EcalFenixLinearizer() [2/2]

◆ ~EcalFenixLinearizer() [2/2]

Member Function Documentation

◆ process() [1/4]

◆ process() [2/4]

◆ process() [3/4]

◆ process() [4/4]

◆ setInput() [1/2]

◆ setInput() [2/2]

◆ setParameters() [1/2]

◆ setParameters() [2/2]

Member Data Documentation

◆ badXStatus_

◆ base_

◆ famos_

◆ gainID_

◆ init_

◆ linConsts_

◆ mult_

◆ peds_

◆ shift_

◆ strip_

◆ uncorrectedSample_

◆ vectorbadXStatus_