|
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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
|
#include "DD4hep/DetFactoryHelper.h"
#include "DataFormats/Math/interface/CMSUnits.h"
#include "DetectorDescription/DDCMS/interface/DDPlugins.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
using namespace std;
using namespace dd4hep;
using namespace cms;
using namespace cms_units::operators;
static long algorithm(Detector& /* description */, cms::DDParsingContext& ctxt, xml_h e) {
cms::DDNamespace ns(ctxt, e, true);
DDAlgoArguments args(ctxt, e);
vector<string> sectorNumber = args.vecStr("SectorNumber"); // Id. Number of the sectors
double sectorRin = args.dble("SectorRin"); // Inner radius of service sectors
double sectorRout = args.dble("SectorRout"); // Outer radius of service sectors
double sectorDz = args.dble("SectorDz"); // Sector half-length
double sectorDeltaPhi_B = args.dble("SectorDeltaPhi_B"); // Sector B phi width [A=C=0.5*(360/sectors)]
vector<double> sectorStartPhi = args.vecDble("SectorStartPhi"); // Starting phi for the service sectors
vector<string> sectorMaterial_A = args.vecStr("SectorMaterial_A"); // Material for the A sectors
vector<string> sectorMaterial_B = args.vecStr("SectorMaterial_B"); // Material for the B sectors
vector<string> sectorMaterial_C = args.vecStr("SectorMaterial_C"); // Material for the C sectors
for (int i = 0; i < (int)(sectorNumber.size()); i++)
edm::LogVerbatim("TOBGeom") << "DDTOBAxCableAlgo debug: sectorNumber[" << i << "] = " << sectorNumber[i];
edm::LogVerbatim("TOBGeom") << "DDTOBAxCableAlgo debug: Axial Service Sectors half-length " << sectorDz << "\tRin "
<< sectorRin << "\tRout = " << sectorRout << "\tPhi of sectors position:";
for (int i = 0; i < (int)(sectorNumber.size()); i++)
edm::LogVerbatim("TOBGeom") << "\t[" << i << "]\tPhi = " << sectorStartPhi[i];
edm::LogVerbatim("TOBGeom") << "DDTOBAxCableAlgo debug: List of materials for the sectors/3 parts";
//
edm::LogVerbatim("TOBGeom") << "DDTOBAxCableAlgo debug: Sector/3 A";
for (int i = 0; i < (int)(sectorNumber.size()); i++)
edm::LogVerbatim("TOBGeom") << "\t[" << i << "]\tsectorMaterial_A = " << sectorMaterial_A[i];
//
edm::LogVerbatim("TOBGeom") << "DDTOBAxCableAlgo debug: Sector/3 B";
for (int i = 0; i < (int)(sectorNumber.size()); i++)
edm::LogVerbatim("TOBGeom") << "\t[" << i << "]\tsectorMaterial_B = " << sectorMaterial_B[i];
//
edm::LogVerbatim("TOBGeom") << "DDTOBAxCableAlgo debug: Sector/3 C";
for (int i = 0; i < (int)(sectorNumber.size()); i++)
edm::LogVerbatim("TOBGeom") << "\t[" << i << "]\tsectorMaterial_C = " << sectorMaterial_C[i];
string tubsName = args.parentName();
Volume tubsVol = ns.volume(tubsName);
// Loop over sectors (sectorNumber vector)
for (int i = 0; i < (int)(sectorNumber.size()); i++) {
Solid solid;
string name;
double dz, rin, rout, startphi, widthphi, deltaphi;
// Axial Services
// Each sector is divided in 3 parts from phi[i] to phi[i+1]
widthphi = ((i + 1 == (int)(sectorStartPhi.size())) ? (sectorStartPhi[0] + 2_pi) - sectorStartPhi[i]
: (sectorStartPhi[i + 1] - sectorStartPhi[i]));
// First Part: A
name = "TOBAxService_" + sectorNumber[i] + "A";
dz = sectorDz;
rin = sectorRin;
rout = sectorRout;
startphi = sectorStartPhi[i];
deltaphi = 0.5 * (widthphi - sectorDeltaPhi_B);
solid = ns.addSolid(ns.prepend(name), Tube(rin, rout, dz, startphi, startphi + deltaphi));
edm::LogVerbatim("TOBGeom") << solid.name() << " Tubs made of " << sectorMaterial_A[i] << " from "
<< convertRadToDeg(startphi) << " to " << convertRadToDeg((startphi + deltaphi))
<< " with Rin " << rin << " Rout " << rout << " ZHalf " << dz;
Volume sectorLogic = ns.addVolume(Volume(solid.name(), solid, ns.material(sectorMaterial_A[i])));
tubsVol.placeVolume(sectorLogic, i + 1); // copyNr: i+1
edm::LogVerbatim("TOBGeom") << sectorLogic.name() << " number " << i + 1 << " positioned in " << tubsName
<< " with no translation and no rotation";
// Second Part: B
name = "TOBAxService_" + sectorNumber[i] + "B";
startphi += deltaphi;
deltaphi = sectorDeltaPhi_B;
solid = ns.addSolid(ns.prepend(name), Tube(rin, rout, dz, startphi, startphi + deltaphi));
edm::LogVerbatim("TOBGeom") << solid.name() << " Tubs made of " << sectorMaterial_B[i] << " from "
<< convertRadToDeg(startphi) << " to " << convertRadToDeg((startphi + deltaphi))
<< " with Rin " << rin << " Rout " << rout << " ZHalf " << dz;
sectorLogic = ns.addVolume(Volume(solid.name(), solid, ns.material(sectorMaterial_B[i])));
tubsVol.placeVolume(sectorLogic, i + 1); // copyNr: i+1
edm::LogVerbatim("TOBGeom") << sectorLogic.name() << " number " << i + 1 << " positioned in " << tubsName
<< " with no translation and no rotation";
// Third Part: C
name = "TOBAxService_" + sectorNumber[i] + "C";
startphi += deltaphi;
deltaphi = 0.5 * (widthphi - sectorDeltaPhi_B);
solid = ns.addSolid(ns.prepend(name), Tube(rin, rout, dz, startphi, startphi + deltaphi));
edm::LogVerbatim("TOBGeom") << solid.name() << " Tubs made of " << sectorMaterial_C[i] << " from "
<< convertRadToDeg(startphi) << " to " << convertRadToDeg((startphi + deltaphi))
<< " with Rin " << rin << " Rout " << rout << " ZHalf " << dz;
sectorLogic = ns.addVolume(Volume(solid.name(), solid, ns.material(sectorMaterial_C[i])));
tubsVol.placeVolume(sectorLogic, i + 1); // copyNr: i+1
edm::LogVerbatim("TOBGeom") << sectorLogic.name() << " number " << i + 1 << " positioned in " << tubsName
<< " with no translation and no rotation";
}
edm::LogVerbatim("TOBGeom") << "<<== End of DDTOBAxCableAlgo construction ...";
return 1;
}
// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_track_DDTOBAxCableAlgo, algorithm)
|