CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiStripNoises.h
Go to the documentation of this file.
1 #ifndef SiStripNoises_h
2 #define SiStripNoises_h
3 
5 
6 #include<vector>
7 #include<utility>
8 #include<iostream>
9 #include<boost/cstdint.hpp>
10 
11 #include<cassert>
12 #include<cstring>
13 
24 {
25  public:
26 
27  struct ratioData{
28  uint32_t detid;
29  std::vector<float> values;
30  };
31 
32  struct DetRegistry{
33  uint32_t detid;
34  uint32_t ibegin;
35  uint32_t iend;
36 
38 };
39 
41  {
42  public:
43  bool operator() (const DetRegistry& p,const uint32_t& i) const {return p.detid < i;}
44  };
45 
46  typedef std::vector<unsigned char> Container;
47  typedef std::vector<unsigned char>::const_iterator ContainerIterator;
48  typedef std::pair<ContainerIterator, ContainerIterator> Range;
49  typedef std::vector<DetRegistry> Registry;
50  typedef Registry::const_iterator RegistryIterator;
51  typedef std::vector<uint16_t> InputVector;
52 
56 
57  bool put(const uint32_t& detID,const InputVector &input);
58  const Range getRange(const uint32_t detID) const;
59  Range getRangeByPos(unsigned short pos) const;
60  void getDetIds(std::vector<uint32_t>& DetIds_) const;
61 
62  ContainerIterator getDataVectorBegin() const {return v_noises.begin();}
66 
67  static inline float getNoiseFast(const uint16_t& strip, const Range& range) {
68  return 0.1f*float(decode(strip,range));
69  }
70 
71 #ifdef EDM_ML_DEBUG
72  static float getNoise(uint16_t strip, const Range& range);
73 #else
74  static float getNoise(uint16_t strip, const Range& range) { return getNoiseFast(strip,range);}
75 #endif
76 
77 
78  void allNoises (std::vector<float> & noises, const Range& range) const;
79  void setData(float noise_, InputVector& vped);
80 
81  void printDebug(std::stringstream& ss) const;
82  void printSummary(std::stringstream& ss) const;
83 
84  std::vector<ratioData> operator / (const SiStripNoises& d) ;
85 
86  private:
87  static void encode(const InputVector& Vi, std::vector<unsigned char>& Vo_CHAR);
88 
89  static inline uint16_t decode (uint16_t strip, const Range& range);
90 
93  static inline uint16_t get9bits(const uint8_t * &ptr, int8_t skip);
94 
97 
98 
99  /*
100  const std::string print_as_binary(const uint8_t ch) const;
101  std::string print_char_as_binary(const unsigned char ch) const;
102  std::string print_short_as_binary(const short ch) const;
103  */
104 
106 };
107 
110 inline uint16_t SiStripNoises::get9bits(const uint8_t * &ptr, int8_t skip) {
111  uint8_t maskThis = (0xFF << skip);
112  uint8_t maskThat = ((2 << skip) - 1);
113  uint16_t ret = ( ((*ptr) & maskThis) >> skip );
114  --ptr;
115  return ret | ( ((*ptr) & maskThat) << (8 - skip) );
116 }
117 
118 inline uint16_t SiStripNoises::decode (uint16_t strip, const Range& range) {
119  const unsigned char *data = &*(range.second -1); // pointer to the last byte of data
120  static const uint16_t BITS_PER_STRIP = 9;
121 
122  uint32_t lowBit = strip * BITS_PER_STRIP;
123  uint8_t firstByteBit = (lowBit & 7);//module 8
124 
125  uint16_t vin = uint16_t(*(data-lowBit/8)) | (uint16_t(*(data-lowBit/8-1))<<8);
126  vin = vin >> firstByteBit; vin &= 0x1FF;
127  return vin;
128 
129  /*
130  uint8_t firstByteNBits = 8 - firstByteBit;
131  uint8_t firstByteMask = 0xffu << firstByteBit;
132  uint8_t secondByteMask = ~(0xffu << (BITS_PER_STRIP - firstByteNBits));
133  uint16_t value = ((uint16_t(*(data-lowBit/8 )) & firstByteMask) >> firstByteBit) | ((uint16_t(*(data-lowBit/8-1)) & secondByteMask) << firstByteNBits);
134 
135  if(vin!=value) std::cout << vin << ',' <<value << std::endl;
136  */
137  /*
138  if(strip < 25){
139  std::cout << "***************DECODE*********************"<<"\n"
140  << "strip "<<strip << " "
141  << value
142  <<"\t :"<<print_as_binary(value)
143  <<"\t :"<<print_as_binary( ((uint16_t(*(data-lowBit/8 )) & firstByteMask) >> firstByteBit) )
144  << "-"<<print_as_binary( ((uint16_t(*(data-lowBit/8-1)) & secondByteMask) <<firstByteNBits) )
145  << "\t *(data-lowBit/8) " << print_as_binary( *(data-lowBit/8 ))
146  << "\t *(data-lowBit/8-1) " << print_as_binary( *(data-lowBit/8 -1 ))
147  << "\tlowBit:"<< lowBit
148  << "\tfirstByteMask :"<<print_as_binary(firstByteMask)
149  << "\tsecondByteMask:"<<print_as_binary(secondByteMask)
150  << "\tfirstByteBit:"<<print_as_binary(firstByteBit)
151  << std::endl;
152  }
153  */
154  //return value;
155 }
156 
157 
158 #endif
static const char noise_[]
int i
Definition: DBlmapReader.cc:9
tuple ret
prodAgent to be discontinued
ContainerIterator getDataVectorEnd() const
Definition: SiStripNoises.h:63
Registry indexes
Definition: SiStripNoises.h:96
static uint16_t get9bits(const uint8_t *&ptr, int8_t skip)
std::vector< float > values
Definition: SiStripNoises.h:29
std::vector< uint16_t > InputVector
Definition: SiStripNoises.h:51
Container v_noises
Definition: SiStripNoises.h:95
std::vector< ratioData > operator/(const SiStripNoises &d)
std::vector< DetRegistry > Registry
Definition: SiStripNoises.h:49
static void encode(const InputVector &Vi, std::vector< unsigned char > &Vo_CHAR)
static std::string const input
Definition: EdmProvDump.cc:44
tuple d
Definition: ztail.py:151
static float getNoiseFast(const uint16_t &strip, const Range &range)
Definition: SiStripNoises.h:67
static float getNoise(uint16_t strip, const Range &range)
Definition: SiStripNoises.h:74
std::vector< unsigned char >::const_iterator ContainerIterator
Definition: SiStripNoises.h:47
static uint16_t decode(uint16_t strip, const Range &range)
Registry::const_iterator RegistryIterator
Definition: SiStripNoises.h:50
std::vector< unsigned char > Container
Definition: SiStripNoises.h:46
bool put(const uint32_t &detID, const InputVector &input)
RegistryIterator getRegistryVectorEnd() const
Definition: SiStripNoises.h:65
void getDetIds(std::vector< uint32_t > &DetIds_) const
RegistryIterator getRegistryVectorBegin() const
Definition: SiStripNoises.h:64
#define COND_SERIALIZABLE
Definition: Serializable.h:37
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
Range getRangeByPos(unsigned short pos) const
const Range getRange(const uint32_t detID) const
void printDebug(std::stringstream &ss) const
void printSummary(std::stringstream &ss) const
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripNoises.h:48
void allNoises(std::vector< float > &noises, const Range &range) const
ContainerIterator getDataVectorBegin() const
Definition: SiStripNoises.h:62
bool operator()(const DetRegistry &p, const uint32_t &i) const
Definition: SiStripNoises.h:43
void setData(float noise_, InputVector &vped)