Skip to content
GitLab
Explore
Sign in
Primary navigation
Search or go to…
Project
R
Reference Detector
Manage
Activity
Members
Labels
Plan
Issues
Issue boards
Milestones
Wiki
Code
Merge requests
Repository
Branches
Commits
Tags
Repository graph
Compare revisions
Snippets
Build
Pipelines
Jobs
Pipeline schedules
Artifacts
Deploy
Releases
Package registry
Container registry
Model registry
Operate
Terraform modules
Analyze
Value stream analytics
Contributor analytics
CI/CD analytics
Repository analytics
Model experiments
Help
Help
Support
GitLab documentation
Compare GitLab plans
Community forum
Contribute to GitLab
Provide feedback
Keyboard shortcuts
?
Snippets
Groups
Projects
This is an archived project. Repository and other project resources are read-only.
Show more breadcrumbs
EIC
detectors
Reference Detector
Merge requests
!26
Fix material
Code
Review changes
Check out branch
Download
Patches
Plain diff
Closed
Fix material
fix_material
into
master
Overview
1
Commits
7
Pipelines
2
Changes
3
Closed
Marshall Scott
requested to merge
fix_material
into
master
4 years ago
Overview
1
Commits
7
Pipelines
2
Changes
3
Expand
Added Ar10CO2 to material list
0
0
Merge request reports
Compare
master
version 1
614867cf
4 years ago
master (base)
and
latest version
latest version
32fc5584
7 commits,
4 years ago
version 1
614867cf
6 commits,
4 years ago
3 files
+
175
−
1
Inline
Compare changes
Side-by-side
Inline
Show whitespace changes
Show one file at a time
Files
3
Search (e.g. *.vue) (Ctrl+P)
src/ci_GEM.cpp
0 → 100644
+
98
−
0
Options
#include
<XML/Helper.h>
///////////////////////////
// Central Ion GEM
///////////////////////////
using
namespace
dd4hep
;
static
Ref_t
createDetector
(
Detector
&
desc
,
xml_h
handle
,
SensitiveDetector
sens
)
{
/*
xml::DetElement detElem = handle;
std::string detName = detElem.nameStr();
int detID = detElem.id();
xml::Component dims = detElem.dimensions();
double SizeZ = dims.sizez(); // Size in Z direction
double ROut = dims.rmax(); // Outer radius
double RIn = dims.rmin(); // Inner radius
double ShiftZ = dims.delta();
double X = dims.x();
double Z = dims.z();
int Nlayers = dims.nlayers();
Material mat = desc.material(detElem.materialStr());
Material air = desc.material("Air");
Tube ci_GEM_GVol_Solid(RIn, ROut, SizeZ / 2., 0., 360 * deg);
Volume detVol("ci_GEM_GVol_Logic", ci_GEM_GVol_Solid, air);
detVol.setVisAttributes(desc.visAttributes(detElem.visStr()));
DetElement det(detName, detID);
Volume motherVol = desc.pickMotherVolume(det);
Transform3D tr(RotationZYX(0., 0., 0.), Position(0., 0., ShiftZ));
PlacedVolume detPV = motherVol.placeVolume(detVol, tr);
det.setPlacement(detPV);
//Adding layers to placed volume
for (cml_coll_t c(detElem, _U(layer)); c; c++)
{
xml_comp_t x_layer = c;
string layer_name = detName + _toString(c.id(), "_layer%d");
Volume layer_vol(layer_name, Tube(c.rmin(), c.rmax(), c.z() ), mat );
layer_vol.
Position layer_pos(0, 0, 0);
detPV.placeVolume(layer_vol, layer_pos);
}
det.setPlacement(detPV);
return det;
*/
xml
::
DetElement
detElem
=
handle
;
std
::
string
detName
=
"ci_GEM"
;
int
detID
=
1
;
xml
::
Component
dims
=
detElem
.
dimensions
();
double
SizeZ
=
30.
*
dd4hep
::
cm
;
// Size in Z direction
double
ROut
=
95.
*
dd4hep
::
cm
;
// Outer radius
double
RIn
=
10.
*
dd4hep
::
cm
;
// Inner radius
double
ShiftZ
=
0.
*
dd4hep
::
cm
;
double
X
=
0.
*
dd4hep
::
cm
;
double
Z
=
0.
*
dd4hep
::
cm
;
int
Nlayers
=
8
;
Material
mat
=
desc
.
material
(
detElem
.
materialStr
());
Material
air
=
desc
.
material
(
"Air"
);
DetElement
det
(
detName
,
detID
);
Volume
motherVol
=
desc
.
pickMotherVolume
(
det
);
Tube
ci_GEM_GVol_Solid
(
RIn
,
ROut
,
SizeZ
/
2.
,
0.
,
360
*
deg
);
Volume
detVol
(
"ci_GEM_GVol_Logic"
,
ci_GEM_GVol_Solid
,
air
);
PlacedVolume
detPV
=
motherVol
.
placeVolume
(
detVol
,
RotationZYX
(
0
,
0
,
0
));
detPV
.
addPhysVolID
(
"system"
,
detID
);
det
.
setPlacement
(
detPV
);
//Adding layers to placed volume
for
(
int
layerNum
=
0
;
layerNum
<
Nlayers
;
layerNum
++
){
xml_comp_t
x_layer
=
c
;
string
layer_name
=
detName
+
_toString
(
layerNum
,
"_layer%d"
);
double
lrmin
,
lrmax
,
lz
;
lrmin
=
RIn
+
(
1.0
+
0.5
*
layerNum
)
*
dd4hep
::
cm
;
lrmax
=
ROut
+
(
-
25.0
+
2.0
*
layerNum
)
*
dd4hep
::
cm
;
lz
=
sizez
*
0.5
+
(
5.0
+
3.
*
layerNum
)
*
dd4hep
::
cm
;
Volume
layer_vol
(
layer_name
,
Tube
(
lrmin
,
lrmax
,
1.0
*
dd4hep
::
cm
),
mat
);
Position
layer_pos
(
0
,
0
,
lz
);
DetElement
layer
(
det
,
layer_name
,
detID
);
layer
.
setVisAttributes
(
desc
.
visAttributes
(
det
.
visStr
()));
//layer.setAttributes(desc, layer_vol, desc.region(), desc.limits(), );//(region, limits, visable).
PlacedVolume
layer_phv
=
detVol
.
placeVolume
(
layer_vol
,
layer_pos
);
layer_phv
.
addPhysVolID
(
"layer"
,
layerNum
);
layer
.
setPlacement
(
layer_phv
);
layerNum
++
;
}
return
det
;
}
// clang-format off
DECLARE_DETELEMENT
(
ci_GEM
,
createDetector
)
Loading