Resolve "Implement ci_GEM"

Closes #20 (closed)

src/ci_GEM.cpp

@@ -7,6 +7,7 @@ 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();
@@ -20,15 +21,17 @@ static Ref_t createDetector(Detector& desc, xml_h handle, SensitiveDetector sens
   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, mat);
+  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
@@ -44,6 +47,56 @@ static Ref_t createDetector(Detector& desc, xml_h handle, SensitiveDetector sens
 
   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);
+
+  //detVol.setVisAttributes(desc.visAttributes(detElem.visStr()));
+  
+ 
+  
+  
+  //Adding layers to placed volume
+  int layerNum = 0;
+  for (cml_coll_t c(detElem, _U(layer)); c; c++)
+  {
+  	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, lz ), mat);
+  	Position layer_pos(0, 0, lz);
+  	DetElement layer(det, layer_name, detID);
+  	//layer.setAttributes(desc, layer_vol, , , );//(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