6 #include "CLHEP/Units/GlobalSystemOfUnits.h" 7 #include "CLHEP/Units/SystemOfUnits.h" 22 #include "Math/GenVector/AxisAngle.h" 23 #include "Math/GenVector/Cartesian3D.h" 24 #include "Math/GenVector/DisplacementVector3D.h" 25 #include "Math/GenVector/Rotation3D.h" 26 #include "Math/GenVector/RotationZ.h" 47 eval.
set(ns,
"length",
"20.*m");
50 eval.
set(ns,
"corner",
"[length]/4.");
53 DDName worldName(
"world",ns);
55 DDName nitrogenName(
"Nitrogen",
"elements");
56 DDName oxygenName(
"Oxygen",
"elements");
59 eval.
eval(ns,
"[length]/2."),
60 eval.
eval(ns,
"[length]/2."));
75 eval.
eval(ns,
"14.007*g/mole"),
76 eval.
eval(ns,
"0.808*g/cm3") );
80 eval.
eval(ns,
"15.999*g/mole"),
81 eval.
eval(ns,
"1.43*g/cm3") );
118 cout <<
"main::initialize DDL parser" << endl;
121 cout <<
"main::about to set configuration" << endl;
126 eval.
eval(ns,
"[setup:corner]/4."),
127 eval.
eval(ns,
"[setup:corner]/8."),
128 eval.
eval(ns,
"[setup:corner]/4.")
131 eval.
eval(ns,
"[setup:corner]/16."),
132 eval.
eval(ns,
"[setup:corner]/16."),
133 eval.
eval(ns,
"[setup:corner]/16.")
137 DDName(
"Oxygen",
"elements"),
141 DDName(
"Nitrogen",
"elements"),
155 eval.
eval(ns,
"[setup:corner]/16."),
156 eval.
eval(ns,
"[setup:corner]/8.")
159 eval.
eval(ns,
"[setup:corner]*1.25*sin(0.)"),
162 eval.
eval(ns,
"[setup:corner]*1.25*sin(30.*deg)"),
165 eval.
eval(ns,
"[setup:corner]*1.25*sin(60.*deg)"),
168 eval.
eval(ns,
"[setup:corner]*1.25*sin(90.*deg)"),
188 os <<
"Starting Regressiontest Output" << endl;
191 cout <<
"main::initialize DDL parser" << endl;
194 cout <<
"main::about to set configuration" << endl;
196 cf.
readConfig(
"DetectorDescription/RegressionTest/test/configuration.xml");
198 cout <<
"main::about to start parsing" << endl;
202 cout <<
"main::completed Parser" << endl;
205 vector<DDTranslation> tvec;
207 std::cout <<
"Before the loop..." << std::endl;
209 ROOT::Math::AxisAngle ra(exv.
rotation());
214 os <<
" " << ra.Axis() << ra.Angle()/deg << endl;
219 vector<DDTranslation>::iterator it = tvec.begin();
220 os << endl <<
"center points of all solids" << endl;
221 for (; it != tvec.end(); ++it) {
222 os << (*it).x() <<
" " << (*it).y() <<
" " << (*it).z() << endl;
229 cout <<
"main:: initialize" << endl;
231 cout <<
"main::initialize DDL parser" << endl;
234 cout <<
"main::about to set configuration" << endl;
237 cf.
readConfig(
"DetectorDescription/RegressionTest/test/configuration.xml");
239 cout <<
"main::about to start parsing" << endl;
243 cout <<
"main::completed Parser" << endl;
245 cout << endl << endl <<
"main::Start checking!" << endl << endl;
250 cout <<
"main::PROBLEM:" << endl
256 std::cout <<
"rot asis\n" << rot << std::endl;
258 rot.GetComponents(x,y,z);
263 cout <<
"phiX=" << x.phi() <<
" or in degrees = " 264 << x.phi()/deg << endl;
265 cout <<
"thetaX=" << x.theta() <<
" or in degrees = " 266 << x.theta()/deg << endl;
267 cout <<
"phiY=" << y.phi() <<
" or in degrees = " 268 << y.phi()/deg << endl;
269 cout <<
"thetaY=" << y.theta() <<
" or in degrees = " 270 << y.theta()/deg << endl;
271 cout <<
"phiZ=" << z.phi() <<
" or in degrees = " 272 << z.phi()/deg << endl;
273 cout <<
"thetaZ=" << z.theta() <<
" or in degrees = " 274 << z.theta()/deg << endl;
276 cout <<
"some factor/equations..." << endl;
277 cout <<
" sin(thetaX()) * cos(phiX()) = " 278 <<
sin(x.theta()) *
cos(x.phi()) << endl;
285 ROOT::Math::AxisAngle aa(
DD3Vector(1.,1.,1.), 20.*deg);
287 cout <<
"DD3Vector was " <<
DD3Vector(1.,1.,1.) <<
" and the rotation was 20*deg around that axis." << endl;
292 cout <<
"(1,0,0, 0,-1,0, 0,0,1)" << endl;
bool next()
set current node to the next node in the expanded tree
const DDRotationMatrix & rotation() const
The absolute rotation of the current node.
int parse(const DDLDocumentProvider &dp)
Parse all files. Return is meaningless.
DDMaterial is used to define and access material information.
Sin< T >::type sin(const T &t)
DDName is used to identify DDD entities uniquely.
const DDSolid & solid(void) const
Returns a reference object of the solid being the shape of this LogicalPart.
int readConfig(const std::string &filename) override
Read in the configuration file.
type of data representation of DDCompactView
void regressionTest_first()
A DDSolid represents the shape of a part.
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
Represents a uniquely identifyable rotation matrix.
void printRot(const DDRotationMatrix &rot)
int addMaterial(const DDMaterial &m, double fm)
adds a material to the mixture proportional to its fraction-mass fm.
T x() const
Cartesian x coordinate.
static value_type & instance()
void regressionTest_setup()
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DD3Vector
A DD Translation is currently implemented with Root Vector3D.
auto const T2 &decltype(t1.eta()) t2
Cos< T >::type cos(const T &t)
A DDLogicalPart aggregates information concerning material, solid and sensitveness ...
static DDSolid box(const DDName &name, double xHalf, double yHalf, double zHalf)
Creates a box with side length 2*xHalf, 2*yHalf, 2*zHalf.
const DDTranslation & translation() const
The absolute translation of the current node.
DDLParser is the main class of Detector Description Language Parser.
FIPConfiguration reads in the configuration file for the DDParser.
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string ©no, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=0)
void output(string filename)
static DDSolid shapeless(const DDName &name)
void set(const std::string &ns, const std::string &name, const std::string &exprValue)
double eval(const std::string &ns, const std::string &expr)
const DDLogicalPart & logicalPart() const
The logical-part of the current node in the expanded-view.
Provides an exploded view of the detector (tree-view)
ROOT::Math::Rotation3D DDRotationMatrix
A DDRotationMatrix is currently implemented with a ROOT Rotation3D.
const DDMaterial & material(void) const
Returns a reference object of the material this LogicalPart is made of.