|
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
|
///////////////////////////////////////////////////////////////////////////////
// File: DDTECCoolAlgo.cc
// Description: Placing cooling pieces in the petal material of a TEC petal
// * in each call all objects are placed at the same radial position.
// * Inserts are placed into the parent object
// * for all i: CoolInsert[i] goes to PhiPosition[i]
///////////////////////////////////////////////////////////////////////////////
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DetectorDescription/Core/interface/DDutils.h"
#include "DetectorDescription/Core/interface/DDCurrentNamespace.h"
#include "DetectorDescription/Core/interface/DDSplit.h"
#include "DetectorDescription/Core/interface/DDTypes.h"
#include "DetectorDescription/Core/interface/DDAlgorithm.h"
#include "DetectorDescription/Core/interface/DDAlgorithmFactory.h"
#include <CLHEP/Units/GlobalPhysicalConstants.h>
#include <CLHEP/Units/SystemOfUnits.h>
#include <cmath>
#include <algorithm>
#include <map>
#include <string>
#include <vector>
using namespace std;
class DDTECCoolAlgo : public DDAlgorithm {
public:
//Constructor and Destructor
DDTECCoolAlgo();
~DDTECCoolAlgo() override;
void initialize(const DDNumericArguments& nArgs,
const DDVectorArguments& vArgs,
const DDMapArguments& mArgs,
const DDStringArguments& sArgs,
const DDStringVectorArguments& vsArgs) override;
void execute(DDCompactView& cpv) override;
private:
string idNameSpace; //Namespace of this and ALL parts
int startCopyNo; //Start copy number
double rPosition; // Position of the Inserts in R
vector<double> phiPosition; // Position of the Inserts in Phi
vector<string> coolInsert; //Name of cooling pieces
};
DDTECCoolAlgo::DDTECCoolAlgo() : phiPosition(0), coolInsert(0) {
LogDebug("TECGeom") << "DDTECCoolAlgo info: Creating an instance";
}
DDTECCoolAlgo::~DDTECCoolAlgo() {}
void DDTECCoolAlgo::initialize(const DDNumericArguments& nArgs,
const DDVectorArguments& vArgs,
const DDMapArguments&,
const DDStringArguments& sArgs,
const DDStringVectorArguments& vsArgs) {
idNameSpace = DDCurrentNamespace::ns();
startCopyNo = int(nArgs["StartCopyNo"]);
DDName parentName = parent().name();
rPosition = nArgs["RPosition"];
LogDebug("TECGeom") << "DDTECCoolAlgo debug: Parent " << parentName << " NameSpace " << idNameSpace
<< " at radial Position " << rPosition;
phiPosition = vArgs["PhiPosition"];
coolInsert = vsArgs["CoolInsert"];
if (phiPosition.size() == coolInsert.size()) {
for (int i = 0; i < (int)(phiPosition.size()); i++)
LogDebug("TECGeom") << "DDTECCoolAlgo debug: Insert[" << i << "]: " << coolInsert.at(i) << " at Phi "
<< phiPosition.at(i) / CLHEP::deg;
} else {
LogDebug("TECGeom") << "ERROR: Number of inserts does not match the numer of PhiPositions!";
}
LogDebug("TECGeom") << " Done creating instance of DDTECCoolAlgo ";
}
void DDTECCoolAlgo::execute(DDCompactView& cpv) {
LogDebug("TECGeom") << "==>> Constructing DDTECCoolAlgo...";
int copyNo = startCopyNo;
// loop over the inserts to be placed
for (int i = 0; i < (int)(coolInsert.size()); i++) {
// get objects
DDName child = DDName(DDSplit(coolInsert.at(i)).first, DDSplit(coolInsert.at(i)).second);
DDName mother = parent().name();
// get positions
double xpos = rPosition * cos(phiPosition.at(i));
double ypos = -rPosition * sin(phiPosition.at(i));
// place inserts
DDTranslation tran(xpos, ypos, 0.0);
DDRotation rotation;
cpv.position(child, mother, copyNo, tran, rotation);
LogDebug("TECGeom") << "DDTECCoolAlgo test " << child << "[" << copyNo << "] positioned in " << mother << " at "
<< tran << " with " << rotation << " phi " << phiPosition.at(i) / CLHEP::deg << " r "
<< rPosition;
copyNo++;
}
LogDebug("TECGeom") << "<<== End of DDTECCoolAlgo construction ...";
}
DEFINE_EDM_PLUGIN(DDAlgorithmFactory, DDTECCoolAlgo, "track:DDTECCoolAlgo");
|