Newer
Older
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
#include "DataHelper/GeomeryHelper.h"
#include "Identifier/CEPCConf.h"
//Externals
#include "DD4hep/DD4hepUnits.h"
#include "DD4hep/Detector.h"
#include "DD4hep/DetElement.h"
#include "DDRec/ISurface.h"
#include "DDRec/SurfaceManager.h"
#include "DDRec/Vector3D.h"
#include "TMatrixF.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/Vector/Rotation.h"
#include <bitset>
// must initialize before using Get()!!!!
const GeomeryWire* GeomeryWire::Instance = nullptr;
const GeomeryWire* GeomeryWire::Get(){
if(nullptr==Instance){
Instance = new GeomeryWire();
}
return Instance;
}
GeomeryWire::GeomeryWire()
{
//m_geosvc = service<IGeomSvc>("GeomSvc");
m_geomsvc = Gaudi::svcLocator()->service("GeomSvc");
m_dd4hepDetector=m_geomsvc->lcdd();
dd4hep::Readout readout=m_dd4hepDetector->readout("DriftChamberHitsCollection");
m_gridDriftChamber=dynamic_cast<dd4hep::DDSegmentation::GridDriftChamber*>(readout.segmentation().segmentation());
if(nullptr==m_gridDriftChamber){
std::cout << "Failed to get the GridDriftChamber" << std::endl;
}
Instance = this;
}
GeomeryWire::~GeomeryWire()
{
}
void GeomeryWire::getWirePos(int layerID,int wireID,TVector3& Wstart,TVector3& Wend) const
{
assert(m_gridDriftChamber);
m_gridDriftChamber->cellposition2(0,layerID,wireID,Wstart,Wend);
}
dd4hep::DDSegmentation::GridDriftChamber* GeomeryWire::getSeg() const
{
return m_gridDriftChamber;
}