Resolve "Implement ci_GEM"

Closes #20 (closed)

src/ci_GEM.cpp

+#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)