|
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
|
#ifndef EcalTestBeam_EcalTBCrystalMap_h
#define EcalTestBeam_EcalTBCrystalMap_h
/*
*
*
*/
#include <map>
#include <iostream>
#include <fstream>
#include <string>
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
class EcalTBCrystalMap {
public:
typedef std::map<std::pair<double, double>, int> CrystalTBIndexMap;
EcalTBCrystalMap(std::string const& MapFileName);
~EcalTBCrystalMap();
int CrystalIndex(double thisEta, double thisPhi);
void findCrystalAngles(const int thisCrysIndex, double& thisEta, double& thisPhi);
static const int NCRYSTAL = 1700;
private:
double crysEta, crysPhi;
int crysIndex;
CrystalTBIndexMap map_;
};
#endif
|