|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
/**
* \file
* A test for EEDetId::hashedIndex()
*
*/
#include <iostream>
#include <string>
#include <stdexcept>
#include <assert.h>
#include <cstdio>
#include "DataFormats/EcalDetId/interface/EBDetId.h"
int main(int argc, char *argv[]) {
FILE *ofile = fopen("eb_next_to_boundary.dat", "w");
int hi = -1;
try {
for (int ieta = EBDetId::MIN_IETA; ieta <= EBDetId::MAX_IETA; ieta++) {
for (int iphi = EBDetId::MIN_IPHI; iphi <= EBDetId::MAX_IPHI; iphi++) {
if (EBDetId::validDetId(ieta, iphi)) {
EBDetId id = EBDetId(ieta, iphi);
if (EBDetId::isNextToBoundary(id)) {
fprintf(ofile, "%d %d %d %d\n", id.ieta(), id.iphi(), hi, 1);
} else {
fprintf(ofile, "%d %d %d %d\n", id.ieta(), id.iphi(), hi, 0);
}
}
if (EBDetId::validDetId(-ieta, iphi)) {
EBDetId id = EBDetId(-ieta, iphi);
if (EBDetId::isNextToBoundary(id)) {
fprintf(ofile, "%d %d %d %d\n", id.ieta(), id.iphi(), hi, 1);
} else {
fprintf(ofile, "%d %d %d %d\n", id.ieta(), id.iphi(), hi, 0);
}
}
}
}
} catch (std::exception &e) {
std::cerr << e.what();
}
fclose(ofile);
}
// to plot the output file:
//gnuplot> set terminal postscript enhanced eps color colourtext dl 1.0 lw 1.5 "Helvetica" 21
//gnuplot> set out 'ee_numbering.eps'
//gnuplot> set size square
//gnuplot> set xlabel 'ix'
//gnuplot> set ylabel 'iy'
//gnuplot> p [0:101][0:101] 'ee_numbering.dat' u 1:($5>0 ? $2 : 1/0) not, '' u 1:($5==0 ? $2 : 1/0) not
|