/DetectorDescription/Parser/src/DDDividedCons.cc

https://github.com/aivanov-cern/cmssw · C++ · 323 lines · 260 code · 46 blank · 17 comment · 22 complexity · 2bf37b7ce7832c84cd52e880d5a090f2 MD5 · raw file 

    

  1. //
    2. 

  2. // ********************************************************************
    3. 

  3. // 25.04.04 - M. Case ddd-ize G4ParameterisationCons*
    4. 

  4. // ********************************************************************
    5. 

  5. 

  6. #include "DetectorDescription/Parser/src/DDDividedCons.h"
    7. 

  7. #include "DetectorDescription/Parser/src/DDXMLElement.h"
    8. 

  8. 

  9. #include "DetectorDescription/Core/interface/DDLogicalPart.h"
    10. 

  10. #include "DetectorDescription/Core/interface/DDName.h"
    11. 

  11. #include "DetectorDescription/Core/interface/DDAxes.h"
    12. 

  12. #include "DetectorDescription/Core/interface/DDSolid.h"
    13. 

  13. #include "DetectorDescription/Core/interface/DDMaterial.h"
    14. 

  14. 

  15. #include "DetectorDescription/Base/interface/DDdebug.h"
    16. 

  16. #include "DetectorDescription/Base/interface/DDRotationMatrix.h"
    17. 

  17. 

  18. #include "CLHEP/Units/GlobalSystemOfUnits.h"
    19. 

  19. 

  20. DDDividedConsRho::DDDividedConsRho( const DDDivision& div, DDCompactView* cpv )
    21. 

  21.   : DDDividedGeometryObject::DDDividedGeometryObject( div, cpv )
    22. 

  22. {
    23. 

  23.   checkParametersValidity();
    24. 

  24.   setType( "DivisionConsRho" );
    25. 

  25.   DDCons msol = (DDCons)(div_.parent().solid());
    26. 

  26. 

  27.   if( msol.rInPlusZ() == 0. )
    28. 

  28.   {
    29. 

  29.     std::cout << "WARNING - DDDividedConsRho, rOutMinusZ = 0. "
    30. 

  30. 	      << std::endl 
    31. 

  31. 	      << "          Width is calculated as that of rOutMinusZ !"
    32. 

  32. 	      << std::endl;
    33. 

  33.   } 
    34. 

  34.       
    35. 

  35.   if( divisionType_ == DivWIDTH )
    36. 

  36.   {
    37. 

  37.     compNDiv_ = calculateNDiv( msol.rOutMinusZ()
    38. 

  38. 			       - msol.rInMinusZ(), div_.width(), div_.offset() );
    39. 

  39.   }
    40. 

  40.   else if( divisionType_ == DivNDIV )
    41. 

  41.   {
    42. 

  42.     DDCons msol = (DDCons)(div_.parent().solid());
    43. 

  43.     compWidth_ = calculateWidth( msol.rOutMinusZ() - msol.rInMinusZ()
    44. 

  44. 				 , div_.nReplicas()
    45. 

  45. 				 , div_.offset() );
    46. 

  46.   }
    47. 

  47. 

  48.   DCOUT_V( 'P', " DDDividedConsRho - no divisions " << compNDiv_ << " = " << div_.nReplicas()
    49. 

  49. 	   << "\n Offset " << div_.offset()
    50. 

  50. 	   << "\n Width " << compWidth_ << " = " << div_.width()
    51. 

  51. 	   << "\n DivType " << divisionType_);  
    52. 

  52. }
    53. 

  53. 

  54. DDDividedConsRho::~DDDividedConsRho( void )
    55. 

  55. {}
    56. 

  56. 

  57. double
    58. 

  58. DDDividedConsRho::getMaxParameter( void ) const
    59. 

  59. {
    60. 

  60.   DDCons msol = (DDCons)(div_.parent().solid());
    61. 

  61.   return msol.rOutMinusZ() - msol.rInMinusZ();
    62. 

  62. 

  63. }
    64. 

  64. 

  65. DDRotation
    66. 

  66. DDDividedConsRho::makeDDRotation( const int copyNo ) const
    67. 

  67. {
    68. 

  68.   DDRotation myddrot; // sets to identity.
    69. 

  69.   DCOUT_V ('P', "DDDividedConsRho::makeDDRotation : " << myddrot);
    70. 

  70.   return myddrot;
    71. 

  71. }
    72. 

  72. 

  73. DDTranslation
    74. 

  74. DDDividedConsRho::makeDDTranslation( const int copyNo ) const
    75. 

  75. {
    76. 

  76.   //----- translation 
    77. 

  77.   DDTranslation translation;
    78. 

  78.   DCOUT_V ('P', " DDDividedConsRho " << "\n\t Position: " << translation << " - Width: " << compWidth_ << " - Axis " << DDAxesNames::name(div_.axis()));
    79. 

  79.   return translation;
    80. 

  80. }
    81. 

  81. 

  82. DDLogicalPart
    83. 

  83. DDDividedConsRho::makeDDLogicalPart( const int copyNo ) const
    84. 

  84. { 
    85. 

  85.   DDName solname(div_.parent().ddname().name() + "_DIVCHILD" 
    86. 

  86. 		 + DDXMLElement::itostr(copyNo) 
    87. 

  87. 		 , div_.parent().ddname().ns());
    88. 

  88.   DDSolid ddcons(solname);
    89. 

  89.   DDMaterial usemat(div_.parent().material());
    90. 

  90.   DDCons msol = (DDCons)(div_.parent().solid());
    91. 

  91. 

  92.   double pRMin1 = msol.rInMinusZ() + div_.offset() + compWidth_ * copyNo;
    93. 

  93.   double pRMax1 = msol.rInMinusZ() + div_.offset() + compWidth_ * (copyNo+1);
    94. 

  94.  
    95. 

  95.   //width at Z Plus
    96. 

  96.   //- double compWidth_Plus =
    97. 

  97.   //   compWidth_ * ( msol.rOutPlusZ()/ msol.rInPlusZ())
    98. 

  98.   //-         / ( msol.rOutMinusZ() - msol.rInMinusZ());
    99. 

  99.   double compWidth_Plus = calculateWidth( msol.rOutPlusZ()
    100. 

  100. 					  - msol.rInPlusZ(), compNDiv_, div_.offset() );
    101. 

  101.   double pRMin2 = msol.rInPlusZ()
    102. 

  102. 		  + div_.offset() + compWidth_Plus * copyNo;
    103. 

  103.   double pRMax2 = msol.rInPlusZ()
    104. 

  104. 		  + div_.offset() + compWidth_Plus * (copyNo+1);
    105. 

  105.   double pDz = msol.zhalf();
    106. 

  106. 

  107.   double pSPhi = msol.phiFrom();
    108. 

  108.   double pDPhi = msol.deltaPhi();
    109. 

  109. 

  110.   ddcons = DDSolidFactory::cons(DDName(solname), pDz, pRMin1, pRMax1
    111. 

  111. 				, pRMin2, pRMax2, pSPhi, pDPhi);      
    112. 

  112.   
    113. 

  113.   DDLogicalPart ddlp = DDLogicalPart(solname, usemat, ddcons);
    114. 

  114.   DCOUT_V( 'P', " DDDividedConsRho::makeDDLogicalPart() lp:" << ddlp );
    115. 

  115.   return ddlp;
    116. 

  116. }
    117. 

  117. 

  118. //--------------------------------------------------------------------------
    119. 

  119. //--------------------------------------------------------------------------
    120. 

  120. DDDividedConsPhi::DDDividedConsPhi( const DDDivision& div, DDCompactView* cpv )
    121. 

  121.   : DDDividedGeometryObject::DDDividedGeometryObject( div, cpv )
    122. 

  122. {
    123. 

  123.   checkParametersValidity();
    124. 

  124.   setType( "DivisionConsPhi" );
    125. 

  125.   DDCons msol = (DDCons)(div_.parent().solid());
    126. 

  126. 

  127.   if( divisionType_ == DivWIDTH )
    128. 

  128.   {
    129. 

  129.     DDCons msol = (DDCons)(div_.parent().solid());
    130. 

  130.     //If you divide a tube of 360 degrees the offset displaces the starting angle, but you still fill the 360 degrees
    131. 

  131.     if( msol.deltaPhi() == 360.*deg )
    132. 

  132.     {
    133. 

  133.       compNDiv_ = calculateNDiv( msol.deltaPhi(), div_.width(), 0. );
    134. 

  134.     }
    135. 

  135.     else
    136. 

  136.     {
    137. 

  137.       compNDiv_ = calculateNDiv( msol.deltaPhi(), div_.width(), div_.offset() );
    138. 

  138.     }
    139. 

  139.   }
    140. 

  140.   else if( divisionType_ == DivNDIV )
    141. 

  141.   {
    142. 

  142.     DDCons msol = (DDCons)(div_.parent().solid());
    143. 

  143.     if( msol.deltaPhi() == 360.*deg )
    144. 

  144.     {
    145. 

  145.       compWidth_ = calculateWidth( msol.deltaPhi(), div_.nReplicas(), 0. );
    146. 

  146.     }
    147. 

  147.     else
    148. 

  148.     {
    149. 

  149.       compWidth_ = calculateWidth( msol.deltaPhi(), div_.nReplicas(), div_.offset() );
    150. 

  150.     }
    151. 

  151.   }
    152. 

  152. 

  153.   DCOUT_V( 'P', " DDDividedConsPhi - no divisions " << compNDiv_ << " = " << div_.nReplicas()
    154. 

  154. 	   << "\n Offset " << div_.offset()
    155. 

  155. 	   << "\n Width " << compWidth_ << " = " << div_.width()
    156. 

  156. 	   << "\n DivType " << divisionType_); 
    157. 

  157. }
    158. 

  158. 

  159. DDDividedConsPhi::~DDDividedConsPhi( void )
    160. 

  160. {}
    161. 

  161. 

  162. double
    163. 

  163. DDDividedConsPhi::getMaxParameter( void ) const
    164. 

  164. {
    165. 

  165.   DDCons msol = (DDCons)(div_.parent().solid());
    166. 

  166.   return msol.deltaPhi();
    167. 

  167. }
    168. 

  168. 

  169. DDRotation
    170. 

  170. DDDividedConsPhi::makeDDRotation( const int copyNo ) const
    171. 

  171. {
    172. 

  172.   DDRotation myddrot; // sets to identity.
    173. 

  173.   double posi = ( copyNo - 1 ) * compWidth_;
    174. 

  174.   DDRotationMatrix* rotMat = changeRotMatrix( posi );
    175. 

  175.   // how to name the rotation??
    176. 

  176.   // i hate this crap :-)
    177. 

  177.   DDName ddrotname(div_.parent().ddname().name() + "_DIVCHILD_ROT" + DDXMLElement::itostr(copyNo)
    178. 

  178. 		   , div_.parent().ddname().ns());
    179. 

  179.   myddrot = DDrot(ddrotname, rotMat);
    180. 

  180. 

  181.   DCOUT_V ('P', "DDDividedConsPhi::makeDDRotation : " << myddrot);
    182. 

  182.   return myddrot;
    183. 

  183. }
    184. 

  184. 

  185. DDTranslation
    186. 

  186. DDDividedConsPhi::makeDDTranslation( const int copyNo ) const
    187. 

  187. {
    188. 

  188.   //----- translation 
    189. 

  189.   DDTranslation translation;  
    190. 

  190.   DCOUT_V( 'P', " DDDividedConsPhi " << "\n\t Position: " << translation
    191. 

  191. 	   << " - Width: " << compWidth_
    192. 

  192. 	   << " - Axis " << DDAxesNames::name(div_.axis()));
    193. 

  193.   return translation;
    194. 

  194. }
    195. 

  195. 

  196. DDLogicalPart
    197. 

  197. DDDividedConsPhi::makeDDLogicalPart( const int copyNo ) const
    198. 

  198. { 
    199. 

  199.   DDName solname(div_.parent().ddname().name() + "_DIVCHILD"
    200. 

  200. 		 , div_.parent().ddname().ns());
    201. 

  201.   DDSolid ddcons(solname);
    202. 

  202.   DDMaterial usemat(div_.parent().material());
    203. 

  203.   DDCons msol = (DDCons)(div_.parent().solid());
    204. 

  204. 

  205.   if (!ddcons.isDefined().second)
    206. 

  206.   {
    207. 

  207.     double pRMin1 = msol.rInMinusZ();
    208. 

  208.     double pRMax1 = msol.rOutMinusZ();
    209. 

  209.     double pRMin2 = msol.rInPlusZ();
    210. 

  210.     double pRMax2 = msol.rOutPlusZ();
    211. 

  211.     double pDz = msol.zhalf();
    212. 

  212. 	  
    213. 

  213.     //- already rotated  double pSPhi = div_.offset() + copyNo*compWidth_;
    214. 

  214.     double pSPhi = div_.offset() + msol.phiFrom();
    215. 

  215.     double pDPhi = compWidth_;
    216. 

  216.     ddcons = DDSolidFactory::cons(DDName(solname), pDz, pRMin1, pRMax1
    217. 

  217. 				  , pRMin2, pRMax2, pSPhi, pDPhi);      
    218. 

  218.   }
    219. 

  219.   
    220. 

  220.   DDLogicalPart ddlp = DDLogicalPart(solname, usemat, ddcons);
    221. 

  221.   DCOUT_V ('P', " DDDividedConsPhi::makeDDLogicalPart() lp:" << ddlp);
    222. 

  222.   return ddlp;
    223. 

  223. }
    224. 

  224. 

  225. DDDividedConsZ::DDDividedConsZ( const DDDivision& div, DDCompactView* cpv )
    226. 

  226.   :  DDDividedGeometryObject::DDDividedGeometryObject( div, cpv )
    227. 

  227. {
    228. 

  228.   checkParametersValidity();
    229. 

  229. 

  230.   DDCons msol = (DDCons)(div_.parent().solid());
    231. 

  231.   setType( "DivisionConsZ" );
    232. 

  232.       
    233. 

  233.   if( divisionType_ == DivWIDTH )
    234. 

  234.   {
    235. 

  235.     DDCons msol = (DDCons)(div_.parent().solid());
    236. 

  236.     compNDiv_ = calculateNDiv( 2*msol.zhalf()
    237. 

  237. 			       , div_.width(), div_.offset() );
    238. 

  238.   }
    239. 

  239.   else if( divisionType_ == DivNDIV )
    240. 

  240.   {
    241. 

  241.     DDCons msol = (DDCons)(div_.parent().solid());
    242. 

  242.     compWidth_ = calculateWidth( 2*msol.zhalf()
    243. 

  243. 				 , div_.nReplicas(), div_.offset() );
    244. 

  244.   }
    245. 

  245. 

  246.   DCOUT_V( 'P', " DDDividedConsZ - no divisions " << compNDiv_ << " = " << div_.nReplicas()
    247. 

  247. 	   << "\n Offset " << div_.offset()
    248. 

  248. 	   << "\n Width " << compWidth_ << " = " << div_.width()
    249. 

  249. 	   << "\n DivType " << divisionType_ );
    250. 

  250. }
    251. 

  251. 

  252. DDDividedConsZ::~DDDividedConsZ( void )
    253. 

  253. {}
    254. 

  254. 

  255. double
    256. 

  256. DDDividedConsZ::getMaxParameter( void ) const
    257. 

  257. {
    258. 

  258.   DDCons msol = (DDCons)(div_.parent().solid());
    259. 

  259.   return 2*msol.zhalf();
    260. 

  260. }
    261. 

  261. 

  262. DDRotation
    263. 

  263. DDDividedConsZ::makeDDRotation( const int copyNo ) const
    264. 

  264. {
    265. 

  265.   DDRotation myddrot; // sets to identity.
    266. 

  266.   DCOUT_V ('P', "DDDividedConsZ::makeDDRotation : " << myddrot);
    267. 

  267.   return myddrot;
    268. 

  268. }
    269. 

  269. 

  270. DDTranslation
    271. 

  271. DDDividedConsZ::makeDDTranslation( const int copyNo ) const
    272. 

  272. {
    273. 

  273.   //----- translation 
    274. 

  274.   DDTranslation translation;
    275. 

  275. 

  276.   DDCons motherCons = (DDCons)(div_.parent().solid());
    277. 

  277.   double posi = - motherCons.zhalf() + div_.offset() 
    278. 

  278. 		+ compWidth_/2 + copyNo*compWidth_;
    279. 

  279.   translation.SetZ(posi); 
    280. 

  280.   
    281. 

  281.   DCOUT_V ('P', " DDDividedConsZ " << "\n\t Position: " << translation << " - Width: " << compWidth_ << " - Axis " << DDAxesNames::name(div_.axis()));
    282. 

  282.   return translation;
    283. 

  283. }
    284. 

  284. 

  285. DDLogicalPart
    286. 

  286. DDDividedConsZ::makeDDLogicalPart( const int copyNo ) const
    287. 

  287. { 
    288. 

  288.   DDName solname(div_.parent().ddname().name() + "_DIVCHILD" + DDXMLElement::itostr(copyNo) 
    289. 

  289. 		 , div_.parent().ddname().ns());
    290. 

  290.   DDSolid ddcons(solname);
    291. 

  291.   DDMaterial usemat(div_.parent().material());
    292. 

  292.   DDCons msol = (DDCons)(div_.parent().solid());
    293. 

  293. 

  294.   double mHalfLength = msol.zhalf();
    295. 

  295.   double aRInner = (msol.rInPlusZ()
    296. 

  296. 		    - msol.rInMinusZ()) / (2*mHalfLength);
    297. 

  297.   double bRInner = (msol.rInPlusZ()
    298. 

  298. 		    + msol.rInMinusZ()) / 2;
    299. 

  299.   double aROuter = (msol.rOutPlusZ()
    300. 

  300. 		    - msol.rOutMinusZ()) / (2*mHalfLength);
    301. 

  301.   double bROuter = (msol.rOutPlusZ()
    302. 

  302. 		    + msol.rOutMinusZ()) / 2;
    303. 

  303.   double xMinusZ = -mHalfLength + div_.offset() + compWidth_*copyNo;
    304. 

  304.   double xPlusZ  = -mHalfLength + div_.offset() + compWidth_*(copyNo+1);
    305. 

  305. 

  306.   double pDz = compWidth_ / 2.;
    307. 

  307.   double pSPhi = msol.phiFrom();
    308. 

  308.   double pDPhi = msol.deltaPhi();
    309. 

  309. 

  310.   ddcons = DDSolidFactory::cons(DDName(solname)
    311. 

  311. 				, pDz
    312. 

  312. 				, aRInner * xMinusZ + bRInner
    313. 

  313. 				, aROuter * xMinusZ + bROuter
    314. 

  314. 				, aRInner * xPlusZ + bRInner
    315. 

  315. 				, aROuter * xPlusZ + bROuter
    316. 

  316. 				, pSPhi
    317. 

  317. 				, pDPhi
    318. 

  318.     );
    319. 

  319.   
    320. 

  320.   DDLogicalPart ddlp = DDLogicalPart(solname, usemat, ddcons);
    321. 

  321.   DCOUT_V( 'P', " DDDividedConsZ::makeDDLogicalPart() lp:" << ddlp );
    322. 

  322.   return ddlp;
    323. 

  323. }
    324.