|
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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
#ifndef Geometry_HcalTowerAlgo_HcalLayerDepthMap_h
#define Geometry_HcalTowerAlgo_HcalLayerDepthMap_h
/** \class HcalLayerDepthMap
*
* this class stores the map of layer to depth for special phi sections
*
* \author Sunanda Banerjee, SINP <sunanda.banerjee@cern.ch>
*
*/
#include <iostream>
#include <iomanip>
#include <map>
#include <string>
#include <vector>
class HcalLayerDepthMap {
public:
HcalLayerDepthMap();
~HcalLayerDepthMap();
void initialize(const int subdet,
const int ietaMax,
const int dep16C,
const int dep29C,
const double wtl0C,
std::vector<int> const& iphi,
std::vector<int> const& ieta,
std::vector<int> const& layer,
std::vector<int> const& depth);
int getSubdet() const { return subdet_; }
int getDepth(const int subdet, const int ieta, const int iphi, const int zside, const int layer) const;
int getDepth16(const int subdet, const int iphi, const int zside) const;
int getDepthMin(const int subdet, const int iphi, const int zside) const;
int getDepthMax(const int subdet, const int iphi, const int zside) const;
int getDepthMax(const int subdet, const int ieta, const int iphi, const int zside) const;
std::pair<int, int> getDepths(const int eta) const;
int getLayerFront(const int subdet, const int ieta, const int iphi, const int zside, const int depth) const;
int getLayerBack(const int subdet, const int ieta, const int iphi, const int zside, const int depth) const;
void getLayerDepth(
const int subdet, const int ieta, const int iphi, const int zside, std::map<int, int>& layers) const;
void getLayerDepth(const int ieta, std::map<int, int>& layers) const;
double getLayer0Wt(const int subdet, const int iphi, const int zside) const;
int getMaxDepthLastHE(const int subdet, const int iphi, const int zside) const;
const std::vector<int>& getPhis() const { return iphi_; }
bool isValid(const int det, const int phi, const int zside) const;
int validDet(std::vector<int>& phis) const;
std::pair<int, int> validEta() const { return std::pair<int, int>(ietaMin_, ietaMax_); }
private:
static const int maxLayers_ = 18;
int subdet_; // Subdet (HB=1, HE=2)
int ietaMin_; // Minimum eta value
int ietaMax_; // Maximum eta value
int depthMin_; // Minimum depth
int depthMax_; // Maximum depth
int dep16C_; // Max/Min layer # for HB/HE (ieta=16)
int dep29C_; // Max Depth of the last HE
double wtl0C_; // Layer 0 weight
std::vector<int> iphi_; // phi*zside values
std::map<std::pair<int, int>, int> layer2Depth_; // Layer to depth map
std::map<std::pair<int, int>, int> depth2LayerF_; // Depth to front layer map
std::map<std::pair<int, int>, int> depth2LayerB_; // Depth to back layer map
std::map<int, std::pair<int, int>> depthsEta_; // Depth range for each eta
};
#endif
|