Calibration/src/big.cc

0001 #include "CondFormats/Calibration/interface/big.h"
0002 #include <iostream>
0003 //fill big
0004 void big::fill(size_t tVectorSize, size_t thVectorSize, size_t sVectorSize, const std::string& atitle) {
0005   for (size_t i = 0; i < tVectorSize; ++i) {
0006     big::bigEntry b;
0007     b.fill(i, 1.0);
0008     tVector_.push_back(b);
0009   }
0010   for (size_t i = 0; i < thVectorSize; ++i) {
0011     big::bigHeader h;
0012     h.fill(atitle);
0013     thVector_.push_back(h);
0014   }
0015   for (size_t i = 0; i < sVectorSize; ++i) {
0016     big::bigStore s;
0017     s.fill(atitle);
0018     sVector_.push_back(s);
0019   }
0020 }
0021
0022 //fill bigEntry
0023 void big::bigEntry::fill(int r, float seed) {
0024   runnum = r;
0025   alpha = seed;
0026   cotalpha = seed;
0027   beta = seed;
0028   cotbeta = seed;
0029   costrk[0] = seed * 0.1;
0030   costrk[1] = seed * 0.2;
0031   costrk[2] = seed * 0.3;
0032   qavg = seed;
0033   symax = seed;
0034   dyone = seed;
0035   syone = seed;
0036   sxmax = seed;
0037   dxone = seed;
0038   sxone = seed;
0039   dytwo = seed;
0040   sytwo = seed;
0041   dxtwo = seed;
0042   sxtwo = seed;
0043   qmin = seed;
0044   for (int i = 0; i < parIDX::LEN1; ++i) {
0045     for (int j = 0; j < parIDX::LEN2; ++j) {
0046       for (int k = 0; k < parIDX::LEN3; ++k) {
0047         par[parIDX::indexOf(i, j, k)] = seed;
0048       }
0049     }
0050   }
0051   for (int i = 0; i < ytempIDX::LEN1; ++i) {
0052     for (int j = 0; j < ytempIDX::LEN2; ++j) {
0053       ytemp[ytempIDX::indexOf(i, j)] = seed;
0054     }
0055   }
0056   for (int i = 0; i < xtempIDX::LEN1; ++i) {
0057     for (int j = 0; j < xtempIDX::LEN2; ++j) {
0058       xtemp[xtempIDX::indexOf(i, j)] = seed;
0059     }
0060   }
0061   for (int i = 0; i < avgIDX::LEN1; ++i) {
0062     for (int j = 0; j < avgIDX::LEN2; ++j) {
0063       for (int k = 0; k < avgIDX::LEN3; ++k) {
0064         avg[avgIDX::indexOf(i, j, k)] = seed;
0065       }
0066     }
0067   }
0068   for (int i = 0; i < aqflIDX::LEN1; ++i) {
0069     for (int j = 0; j < aqflIDX::LEN2; ++j) {
0070       for (int k = 0; k < aqflIDX::LEN3; ++k) {
0071         aqfl[aqflIDX::indexOf(i, j, k)] = seed;
0072       }
0073     }
0074   }
0075   for (int i = 0; i < chi2IDX::LEN1; ++i) {
0076     for (int j = 0; j < chi2IDX::LEN2; ++j) {
0077       for (int k = 0; k < chi2IDX::LEN3; ++k) {
0078         chi2[chi2IDX::indexOf(i, j, k)] = seed;
0079       }
0080     }
0081   }
0082   for (int i = 0; i < spareIDX::LEN1; ++i) {
0083     for (int j = 0; j < spareIDX::LEN2; ++j) {
0084       spare[spareIDX::indexOf(i, j)] = seed;
0085     }
0086   }
0087 }
0088
0089 //fill bigHeader
0090 void big::bigHeader::fill(const std::string& atitle) {
0091   title = std::string("atitle");
0092   ID = 0;
0093   NBy = 1;
0094   NByx = 2;
0095   NBxx = 3;
0096   NFy = 4;
0097   NFyx = 5;
0098   NFxx = 6;
0099   vbias = 0.1;
0100   temperature = 0.2;
0101   fluence = 0.3;
0102   qscale = 0.4;
0103   s50 = 0.5;
0104   templ_version = 1;
0105 }
0106 //fill bigStore
0107 void big::bigStore::fill(const std::string& atitle) {
0108   head.fill(atitle);
0109   for (int i = 0; i < entbyIDX::LEN1; ++i) {
0110     bigEntry b;
0111     b.fill(i, 0.5 * i);
0112     entby[entbyIDX::indexOf(i)] = b;  //or use push_back as prefer
0113   }
0114   std::cout << "length of entbx 1 " << entbxIDX::LEN1 << std::endl;
0115   std::cout << "length of entbx 2 " << entbxIDX::LEN2 << std::endl;
0116   std::cout << "total size of entbx " << entbxIDX::SIZE << std::endl;
0117   for (int i = 0; i < entbxIDX::LEN1; ++i) {
0118     for (int j = 0; j < entbxIDX::LEN2; ++j) {
0119       bigEntry c;
0120       c.fill(i * j, 0.3 * j);
0121       entbx[entbxIDX::indexOf(i, j)] = c;  //or use push_back as prefer
0122     }
0123   }
0124   for (int i = 0; i < entfyIDX::LEN1; ++i) {
0125     bigEntry f;
0126     f.fill(i, 0.4 * i);
0127     entfy[entfyIDX::indexOf(i)] = f;  //or use push_back as prefer
0128   }
0129   for (int i = 0; i < entfxIDX::LEN1; ++i) {
0130     for (int j = 0; j < entfxIDX::LEN2; ++j) {
0131       bigEntry f;
0132       f.fill(i * j, 0.25 * j);
0133       entfx[entfxIDX::indexOf(i, j)] = f;  //or use push_back as prefer
0134     }
0135   }
0136 }