|
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
68
69
70
71
72
73
74
75
76
|
#ifndef FastSimulation_CaloGeometryTools_CaloGeometryHelper
#define FastSimulation_CaloGeometryTools_CaloGeometryHelper
#include "FastSimulation/CalorimeterProperties/interface/Calorimeter.h"
#include "Geometry/CaloTopology/interface/CaloDirection.h"
#include "FastSimulation/CaloGeometryTools/interface/BaseCrystal.h"
#include <vector>
#include <array>
class DetId;
class Crystal;
namespace edm {
class ParameterSet;
}
class CaloGeometryHelper : public Calorimeter {
public:
typedef math::XYZVector XYZVector;
typedef math::XYZVector XYZPoint;
CaloGeometryHelper();
CaloGeometryHelper(const edm::ParameterSet& fastCalo);
~CaloGeometryHelper();
// more user friendly getClosestCell
DetId getClosestCell(const XYZPoint& point, bool ecal, bool central) const;
// more user friendly getWindow
void getWindow(const DetId& pivot, int s1, int s2, std::vector<DetId>&) const;
double preshowerZPosition(int layer) const { return (layer == 1) ? psLayer1Z_ : psLayer2Z_; }
// the Crystal constructor
void buildCrystal(const DetId& id, Crystal&) const;
void initialize(double bField);
// get the <=8 neighbours
typedef std::array<DetId, 8> NeiVect;
const NeiVect& getNeighbours(const DetId& det) const;
inline double magneticField() const { return bfield_; }
// temporary. Requires a missing geometry tool
bool borderCrossing(const DetId&, const DetId&) const;
bool move(DetId& cell, const CaloDirection& dir, bool fast = true) const;
inline bool preshowerPresent() const { return preshowerPresent_; };
private:
void buildNeighbourArray();
void buildCrystalArray();
bool simplemove(DetId& cell, const CaloDirection& dir) const;
bool diagonalmove(DetId& cell, const CaloDirection& dir) const;
private:
// Preshower layer positions
double psLayer1Z_, psLayer2Z_;
// array of neighbours the hashed index is used for the first vector
std::vector<NeiVect> barrelNeighbours_;
std::vector<NeiVect> endcapNeighbours_;
std::vector<BaseCrystal> barrelCrystals_;
std::vector<BaseCrystal> endcapCrystals_;
bool neighbourmapcalculated_;
//mag field at 0,0,0
double bfield_;
bool preshowerPresent_;
};
#endif
|