ATLAS Offline Software
Loading...
Searching...
No Matches
Functions | Variables
CrackDMConstruction.cxx File Reference
#include "CrackRegionGeoModel/CrackDMConstruction.h"
#include "GeoModelKernel/GeoElement.h"
#include "GeoModelKernel/GeoMaterial.h"
#include "GeoModelKernel/GeoFullPhysVol.h"
#include "GeoModelKernel/GeoPhysVol.h"
#include "GeoModelKernel/GeoVPhysVol.h"
#include "GeoModelKernel/GeoLogVol.h"
#include "GeoModelKernel/GeoTransform.h"
#include "GeoModelKernel/GeoAlignableTransform.h"
#include "GeoModelKernel/GeoIdentifierTag.h"
#include "GeoModelKernel/GeoNameTag.h"
#include "GeoModelKernel/GeoSerialTransformer.h"
#include "GeoModelKernel/GeoXF.h"
#include "GeoModelKernel/GeoTube.h"
#include "GeoModelKernel/GeoPcon.h"
#include "GeoModelKernel/GeoTubs.h"
#include "GeoModelKernel/GeoCons.h"
#include "GeoModelKernel/GeoBox.h"
#include "GeoModelKernel/GeoTrap.h"
#include "GeoModelKernel/GeoTrd.h"
#include "GeoModelKernel/GeoShape.h"
#include "GeoModelKernel/GeoShapeUnion.h"
#include "GeoModelKernel/GeoShapeShift.h"
#include "GeoModelKernel/GeoShapeSubtraction.h"
#include "GeoModelKernel/Units.h"
#include "StoreGate/StoreGateSvc.h"
#include "GeoModelInterfaces/StoredMaterialManager.h"
#include "GeoModelInterfaces/IGeoModelSvc.h"
#include "GeoModelUtilities/DecodeVersionKey.h"
#include "RDBAccessSvc/IRDBAccessSvc.h"
#include "RDBAccessSvc/IRDBRecordset.h"
#include "RDBAccessSvc/IRDBRecord.h"
#include "GaudiKernel/PhysicalConstants.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/Bootstrap.h"
#include "GeoGenericFunctions/Variable.h"

Go to the source code of this file.

Functions

void createSectorEnvelopes2FromDB (GeoFullPhysVol *envelope, StoredMaterialManager *materialManager, std::map< std::string, unsigned int > &trdMap, IRDBRecordset &BarrelDMTrds, std::map< std::string, unsigned int > &trapMap, IRDBRecordset &BarrelDMTraps, std::map< std::string, unsigned int > &boxMap, IRDBRecordset &BarrelDMBoxes, GENFUNCTION &f, GeoBox *Box)
void createBridgeEnvelopesFromDB (GeoFullPhysVol *envelope, std::map< std::string, unsigned int > &trapMap, IRDBRecordset &BarrelDMTraps, const GeoMaterial *matLArServices8, GENFUNCTION &f)
void createBaseEnvelopesFromDB (GeoFullPhysVol *envelope, std::map< std::string, unsigned int > &trdMap, IRDBRecordset &BarrelDMTrds, const GeoMaterial *matLArServices8, GENFUNCTION &f)
void createFromDB (GeoFullPhysVol *envelope, IRDBAccessSvc *rdbAccess, IGeoModelSvc *geoModel, StoredMaterialManager *materialManager)

Variables

static const unsigned int NCrates =16
static const double Alfa =360*Gaudi::Units::deg/NCrates
static const double Enda =1155
static const double Endb =1695.2
static const double Endc =2771.6
static const double DYa =1155*tan(Alfa/2)
static const double DYb =1695.2*tan(Alfa/2)
static const double DYc =2771.6*tan(Alfa/2)

Function Documentation

◆ createBaseEnvelopesFromDB()

void createBaseEnvelopesFromDB ( GeoFullPhysVol * envelope,
std::map< std::string, unsigned int > & trdMap,
IRDBRecordset & BarrelDMTrds,
const GeoMaterial * matLArServices8,
GENFUNCTION & f )

Definition at line 306 of file CrackDMConstruction.cxx.

311{
312 unsigned int recordIndex = trdMap["BaseE"];
313 const IRDBRecord* r = BarrelDMTrds[recordIndex];
314 double BaseExhlen1 = r->getDouble("XHLEN1");
315 double BaseExhlen2 = r->getDouble("XHLEN2");
316 double BaseEyhlen1 = r->getDouble("YHLEN1");
317 double BaseEyhlen2 = r->getDouble("YHLEN2");
318 double BaseEzhlen = r->getDouble("ZHLEN");//
319 double BaseExtr = r->getDouble("XTR");
320 double BaseEztr = r->getDouble("ZTR");
321
322 GeoTrd *Trd1air = new GeoTrd(BaseExhlen1, BaseExhlen2, BaseEyhlen1, BaseEyhlen2, BaseEzhlen);
323 GeoLogVol *lvbe = new GeoLogVol("LAr::DM::BaseEnvelopes",Trd1air,matLArServices8); //In the end Density at least >= than SE1 because of Cryo Pipes
324 GeoPhysVol *baseenvelopes = new GeoPhysVol(lvbe);
325
326 TRANSFUNCTION beA = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(BaseExtr)*GeoTrf::TranslateZ3D(BaseEztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
327 TRANSFUNCTION beC = Pow(GeoTrf::RotateZ3D(1.0),f+(Alfa/2))*GeoTrf::TranslateX3D(BaseExtr)*GeoTrf::TranslateZ3D(-BaseEztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
328 GeoSerialTransformer *betA = new GeoSerialTransformer(baseenvelopes,&beA, NCrates);
329 GeoSerialTransformer *betC = new GeoSerialTransformer(baseenvelopes,&beC, NCrates);
330 envelope->add(betA);
331 envelope->add(betC);
332}
NCrates
static const unsigned int NCrates
Definition CrackDMConstruction.cxx:51
Alfa
static const double Alfa
Definition CrackDMConstruction.cxx:52
IRDBRecord
IRDBRecord is one record in the IRDBRecordset object.
Definition IRDBRecord.h:27
r
int r
Definition globals.cxx:22

◆ createBridgeEnvelopesFromDB()

void createBridgeEnvelopesFromDB ( GeoFullPhysVol * envelope,
std::map< std::string, unsigned int > & trapMap,
IRDBRecordset & BarrelDMTraps,
const GeoMaterial * matLArServices8,
GENFUNCTION & f )

Definition at line 270 of file CrackDMConstruction.cxx.

275{
276 unsigned int recordIndex = trapMap["BridgeE"];
277 const IRDBRecord* r = BarrelDMTraps[recordIndex];
278 double BridgeEzhlen = r->getDouble("ZHLEN");
279 double BridgeEtheta = r->getDouble("THETA");
280 double BridgeEphi = r->getDouble("PHI");
281 double BridgeEyzn = r->getDouble("YZN");
282 double BridgeExynzn = r->getDouble("XYNZN");
283 double BridgeExypzn = r->getDouble("XYPZN");
284 double BridgeEangn = r->getDouble("ANGN");
285 double BridgeEyzp = r->getDouble("YZP");
286 double BridgeExynzp = r->getDouble("XYNZP");
287 double BridgeExypzp = r->getDouble("XYPZP");
288 double BridgeEangp = r->getDouble("ANGP");
289 double BridgeExtr = r->getDouble("XTR");
290 double BridgeEztr = r->getDouble("ZTR");
291
292 GeoTrap *Trapair = new GeoTrap(BridgeEzhlen, BridgeEtheta*Gaudi::Units::deg, BridgeEphi, BridgeEyzn, BridgeExynzn, BridgeExypzn, BridgeEangn, BridgeEyzp, BridgeExynzp, BridgeExypzp, BridgeEangp);
293 GeoLogVol *lvbre = new GeoLogVol("LAr::DM::BridgeEnvelopes",Trapair,matLArServices8);//In the end Density at least >= than SE1 because of Cryo Pipes
294 GeoPhysVol *bridgeenvelopes = new GeoPhysVol(lvbre);
295
296 TRANSFUNCTION breA = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(BridgeExtr)*GeoTrf::TranslateZ3D(BridgeEztr)*GeoTrf::RotateZ3D(90*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)*GeoTrf::RotateX3D(90*Gaudi::Units::deg);
297 TRANSFUNCTION breC = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(BridgeExtr)*GeoTrf::TranslateZ3D(-BridgeEztr)*GeoTrf::RotateZ3D(-90*Gaudi::Units::deg)*GeoTrf::RotateY3D(-90*Gaudi::Units::deg)*GeoTrf::RotateX3D(-90*Gaudi::Units::deg);
298 GeoSerialTransformer *bretA = new GeoSerialTransformer(bridgeenvelopes,&breA, NCrates);
299 GeoSerialTransformer *bretC = new GeoSerialTransformer(bridgeenvelopes,&breC, NCrates);
300 envelope->add(bretA);
301 envelope->add(bretC);
302}

◆ createFromDB()

void createFromDB ( GeoFullPhysVol * envelope,
IRDBAccessSvc * rdbAccess,
IGeoModelSvc * geoModel,
StoredMaterialManager * materialManager )

Definition at line 334 of file CrackDMConstruction.cxx.

338{
339 // Use Geometry Database
340 DecodeVersionKey keyLAr(geoModel,"LAr");
341 IRDBRecordset_ptr BarrelDMTraps = rdbAccess->getRecordsetPtr("BarrelDMTraps",keyLAr.tag(),keyLAr.node());
342 IRDBRecordset_ptr BarrelDMTrds = rdbAccess->getRecordsetPtr("BarrelDMTrds",keyLAr.tag(),keyLAr.node());
343 IRDBRecordset_ptr BarrelDMTubes = rdbAccess->getRecordsetPtr("BarrelDMTubes",keyLAr.tag(),keyLAr.node());
344 IRDBRecordset_ptr BarrelDMBoxes = rdbAccess->getRecordsetPtr("BarrelDMBoxes",keyLAr.tag(),keyLAr.node());
345
346 std::map<std::string, unsigned int> tubeMap;
347 for (unsigned int i=0; i<BarrelDMTubes->size(); i++)
348 {
349 const std::string& key = (*BarrelDMTubes)[i]->getString("TUBENAME");
350 tubeMap[key] = i;
351 }
352 std::map<std::string, unsigned int> boxMap;
353 for (unsigned int j=0; j<BarrelDMBoxes->size(); j++)
354 {
355 const std::string& key = (*BarrelDMBoxes)[j]->getString("BOXNAME");
356 boxMap[key] = j;
357 }
358 std::map<std::string, unsigned int> trdMap;
359 for (unsigned int k=0; k<BarrelDMTrds->size(); k++)
360 {
361 const std::string& key = (*BarrelDMTrds)[k]->getString("TRDNAME");
362 trdMap[key] = k;
363 }
364 std::map<std::string, unsigned int> trapMap;
365 for (unsigned int l=0; l<BarrelDMTraps->size(); l++)
366 {
367 const std::string& key = (*BarrelDMTraps)[l]->getString("TRAPNAME");
368 trapMap[key] = l;
369 }
370
371 unsigned int recordIndex;
372
373 // Get materials
374 const GeoMaterial *alu = materialManager->getMaterial("std::Aluminium"); //2.7 g/cm3
375 const GeoMaterial* matBoardsEnvelope = materialManager->getMaterial("LAr::BoardsEnvelope");// 0.932*gram/cm3);
376 const GeoMaterial* matLArServices1 = materialManager->getMaterial("LAr::LArServices1");// 1.020*gram/cm3
377 const GeoMaterial* matLArServices2 = materialManager->getMaterial("LAr::LArServices2");// 0.955*gram/cm3
378 const GeoMaterial* matLArServices3 = materialManager->getMaterial("LAr::LArServices3");// 1.005*gram/cm3
379 const GeoMaterial* matLArServices4 = materialManager->getMaterial("LAr::LArServices4");// 0.460*gram/cm3
380 const GeoMaterial* matLArServices5 = materialManager->getMaterial("LAr::LArServices5");// 0.480*gram/cm3
381 const GeoMaterial* matLArServices6 = materialManager->getMaterial("LAr::LArServices6");// 1.000*gram/cm3
382 const GeoMaterial* matLArServices7 = materialManager->getMaterial("LAr::LArServices7");// 0.935*gram/cm3
383 const GeoMaterial* matLArServices8 = materialManager->getMaterial("LAr::LArServices8");// 1.070*gram/cm3
384 const GeoMaterial* matLArServices9 = materialManager->getMaterial("LAr::LArServices9");// 1.020*gram/cm3
385 const GeoMaterial* matLArServices10 = materialManager->getMaterial("LAr::LArServices10");// 0.995*gram/cm3
386 const GeoMaterial* matLArServices11 = materialManager->getMaterial("LAr::LArServices11");// 0.835*gram/cm3
387 const GeoMaterial* matLArServices12 = materialManager->getMaterial("LAr::LArServices12");// 0.640*gram/cm3
388 const GeoMaterial* matLArServices13 = materialManager->getMaterial("LAr::LArServices13");// 0.690*gram/cm3
389 const GeoMaterial* matLArServices14 = materialManager->getMaterial("LAr::LArServices14");// 0.825*gram/cm3
390 const GeoMaterial* matLArServices15 = materialManager->getMaterial("LAr::LArServices15");// 0.875*gram/cm3
391 const GeoMaterial* matLArServices16 = materialManager->getMaterial("LAr::LArServices16");// 1.035*gram/cm3
392
393 const double inv_Endab = 1. / (Endb - Enda);
394 Variable i;
395 GENFUNCTION f = Alfa*i;
396
398
399 recordIndex = tubeMap["Ped2"];
400 double ped2zhlen = (*BarrelDMTubes)[recordIndex]->getDouble("ZHLEN");
401 double ped2minr = (*BarrelDMTubes)[recordIndex]->getDouble("MINR");
402 double ped2maxr = (*BarrelDMTubes)[recordIndex]->getDouble("MAXR");
403 double ped2ytr = (*BarrelDMTubes)[recordIndex]->getDouble("YTR");
404 recordIndex = tubeMap["Ped3"];
405 double ped3zhlen = (*BarrelDMTubes)[recordIndex]->getDouble("ZHLEN");
406 double ped3minr = (*BarrelDMTubes)[recordIndex]->getDouble("MINR");
407 double ped3maxr = (*BarrelDMTubes)[recordIndex]->getDouble("MAXR");
408 double ped3xtr = (*BarrelDMTubes)[recordIndex]->getDouble("XTR");
409
410 recordIndex = boxMap["Pedest"];
411 double pedesthlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
412 double pedesthwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
413 double pedesthhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
414 double pedestxtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
415 double pedestztr = (*BarrelDMBoxes)[recordIndex]->getDouble("ZTR");
416 recordIndex = boxMap["Ped1"];
417 double ped1hlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
418 double ped1hwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
419 double ped1hhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
420 recordIndex = boxMap["Crate1"];
421 double crate1hlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
422 double crate1hwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
423 double crate1hhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
424 double crate1xtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
425 double crate1ztr = (*BarrelDMBoxes)[recordIndex]->getDouble("ZTR");
426 recordIndex = boxMap["Crate2"];
427 double crate2hlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
428 double crate2hwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
429 double crate2hhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
430 recordIndex = boxMap["Crate3"];
431 double crate3hlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
432 double crate3hwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
433 double crate3hhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
434 double crate3xtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
435 recordIndex = boxMap["BoardE"];
436 double BoardEhlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
437 double BoardEhwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
438 double BoardEhhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
439 double BoardExtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
440 double BoardEytr = (*BarrelDMBoxes)[recordIndex]->getDouble("YTR");
441 double BoardEztr = (*BarrelDMBoxes)[recordIndex]->getDouble("ZTR");
442 recordIndex = boxMap["Box"];
443 double Boxhlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
444 double Boxhwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
445 double Boxhhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
446 double Boxxtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
447 double Boxytr = (*BarrelDMBoxes)[recordIndex]->getDouble("YTR");
448 double Boxxrot = (*BarrelDMBoxes)[recordIndex]->getDouble("XROT");
449
450 // ----- build pedestal -----
451 GeoBox *Pedestal = new GeoBox(pedesthlen, pedesthwdt, pedesthhgt);
452 GeoBox *Ped1 = new GeoBox(ped1hlen, ped1hwdt, ped1hhgt);
453 GeoTube *Ped2 = new GeoTube(ped2minr, ped2maxr, ped2zhlen);
454 GeoTube *Ped3 = new GeoTube(ped3minr,ped3maxr , ped3zhlen);
455 const GeoShape & CratePed=((*Pedestal).subtract(*Ped1).
456 subtract((*Ped2) <<GeoTrf::TranslateY3D(-ped2ytr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)).
457 subtract((*Ped3) <<GeoTrf::TranslateX3D(-ped3xtr)).
458 subtract((*Ped2) <<GeoTrf::TranslateY3D(ped2ytr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)));
459
460 GeoLogVol *lvped = new GeoLogVol("LAr::DM::Ped",&CratePed,alu);
461 GeoPhysVol *pedestal = new GeoPhysVol(lvped);
462
463 // ----- build crates -----
464 GeoBox *Crate1 = new GeoBox(crate1hlen, crate1hwdt, crate1hhgt);
465 GeoBox *Crate2 = new GeoBox(crate2hlen, crate2hwdt, crate2hhgt);
466 GeoBox *Crate3 = new GeoBox(crate3hlen, crate3hwdt, crate3hhgt);
467 const GeoShape & FEBCrate=(*Crate1).subtract(*Crate2).add((*Crate3) <<GeoTrf::TranslateX3D(-crate3xtr));
468
469 GeoLogVol *lvcrate = new GeoLogVol("LAr::DM::Crate",&FEBCrate,alu);
470 GeoPhysVol *crate = new GeoPhysVol(lvcrate);
471
472 // ----- build boardenvelopes -----
473 GeoBox *BoardEnvelope = new GeoBox(BoardEhlen, BoardEhwdt, BoardEhhgt);
474 GeoLogVol *lvbenv = new GeoLogVol("LAr::DM::FEBoard",BoardEnvelope,matBoardsEnvelope);
475 GeoPhysVol *boardenvelope = new GeoPhysVol(lvbenv);
476
477 //-------------- Place volumes in envelope ----------------------------
478
479 //Crates
480 TRANSFUNCTION crA = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateX3D(crate1xtr)*GeoTrf::TranslateZ3D(crate1ztr);
481 TRANSFUNCTION crC = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateX3D(crate1xtr)*GeoTrf::TranslateZ3D(-crate1ztr);
482 GeoSerialTransformer *crtA = new GeoSerialTransformer(crate,&crA, NCrates);
483 GeoSerialTransformer *crtC = new GeoSerialTransformer(crate,&crC, NCrates);
484 envelope->add(crtA);
485 envelope->add(crtC);
486
487 //Pedestals
488 TRANSFUNCTION pedA = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateX3D(pedestxtr)*GeoTrf::TranslateZ3D(pedestztr);
489 TRANSFUNCTION pedC = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateX3D(pedestxtr)*GeoTrf::TranslateZ3D(-pedestztr);
490 GeoSerialTransformer *pedtA = new GeoSerialTransformer(pedestal,&pedA, NCrates);
491 GeoSerialTransformer *pedtC = new GeoSerialTransformer(pedestal,&pedC, NCrates);
492 envelope->add(pedtA);
493 envelope->add(pedtC);
494
495 //FEBoards
496 TRANSFUNCTION feb1A = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateY3D(BoardEytr)*GeoTrf::TranslateX3D(BoardExtr)*GeoTrf::TranslateZ3D(BoardEztr);
497 TRANSFUNCTION feb2A = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateY3D(-BoardEytr)*GeoTrf::TranslateX3D(BoardExtr)*GeoTrf::TranslateZ3D(BoardEztr);
498 TRANSFUNCTION feb1C = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateY3D(BoardEytr)*GeoTrf::TranslateX3D(BoardExtr)*GeoTrf::TranslateZ3D(-BoardEztr);
499 TRANSFUNCTION feb2C = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateY3D(-BoardEytr)*GeoTrf::TranslateX3D(BoardExtr)*GeoTrf::TranslateZ3D(-BoardEztr);
500 GeoSerialTransformer *febt1A = new GeoSerialTransformer(boardenvelope,&feb1A, NCrates);
501 GeoSerialTransformer *febt1C = new GeoSerialTransformer(boardenvelope,&feb1C, NCrates);
502 GeoSerialTransformer *febt2A = new GeoSerialTransformer(boardenvelope,&feb2A, NCrates);
503 GeoSerialTransformer *febt2C = new GeoSerialTransformer(boardenvelope,&feb2C, NCrates);
504 envelope->add(febt1A);
505 envelope->add(febt1C);
506 envelope->add(febt2A);
507 envelope->add(febt2C);
508
510 recordIndex = trdMap["SecP"];
511 double SecPxhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN1");
512 double SecPxhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN2");
513 double SecPyhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("YHLEN1");//
514 double SecPyhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("YHLEN2");//
515 double SecPzhlen = (*BarrelDMTrds)[recordIndex]->getDouble("ZHLEN");//
516 double SecPxtr = (*BarrelDMTrds)[recordIndex]->getDouble("XTR");
517 double SecPztr = (*BarrelDMTrds)[recordIndex]->getDouble("ZTR");
518 recordIndex = trdMap["BaseP"];
519 double BasePxhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN1");
520 double BasePxhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN2");
521 double BasePyhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("YHLEN1");
522 double BasePyhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("YHLEN2");
523 double BasePzhlen = (*BarrelDMTrds)[recordIndex]->getDouble("ZHLEN");
524 double BasePxtr = (*BarrelDMTrds)[recordIndex]->getDouble("XTR");
525 double BasePztr = (*BarrelDMTrds)[recordIndex]->getDouble("ZTR");
526 recordIndex = trdMap["SecE1"];
527 double SecE1xhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN1");
528 double SecE1xhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN2");
529 double SecE1ztr = (*BarrelDMTrds)[recordIndex]->getDouble("ZTR");
530 recordIndex = trapMap["BridgeP"];
531 double BridgePzhlen = (*BarrelDMTraps)[recordIndex]->getDouble("ZHLEN");
532 double BridgePtheta = (*BarrelDMTraps)[recordIndex]->getDouble("THETA");
533 double BridgePphi = (*BarrelDMTraps)[recordIndex]->getDouble("PHI");
534 double BridgePyzn = (*BarrelDMTraps)[recordIndex]->getDouble("YZN");
535 double BridgePxynzn = (*BarrelDMTraps)[recordIndex]->getDouble("XYNZN");
536 double BridgePxypzn = (*BarrelDMTraps)[recordIndex]->getDouble("XYPZN");
537 double BridgePangn = (*BarrelDMTraps)[recordIndex]->getDouble("ANGN");
538 double BridgePyzp = (*BarrelDMTraps)[recordIndex]->getDouble("YZP");
539 double BridgePxynzp = (*BarrelDMTraps)[recordIndex]->getDouble("XYNZP");
540 double BridgePxypzp = (*BarrelDMTraps)[recordIndex]->getDouble("XYPZP");
541 double BridgePangp = (*BarrelDMTraps)[recordIndex]->getDouble("ANGP");
542 double BridgePxtr = (*BarrelDMTraps)[recordIndex]->getDouble("XTR");
543 double BridgePztr = (*BarrelDMTraps)[recordIndex]->getDouble("ZTR");
544
545 // transforms
546 GeoBox *Box = new GeoBox(Boxhlen, Boxhwdt, Boxhhgt);
547
548 GeoTrf::Transform3D Cut3Boxp = GeoTrf::Translate3D(Boxxtr, Boxytr, Boxxrot)*GeoTrf::RotateX3D(-20*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
549 GeoTrf::Transform3D Cut4Boxp = GeoTrf::Translate3D(Boxxtr, -Boxytr,Boxxrot)*GeoTrf::RotateX3D(20*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
550
551 // ----- build sector envelopes -----
552 // build 16 instances of SectorEnvelopes1 each with its own material!
553 GeoTrd *Trdair1 = new GeoTrd(SecE1xhlen1, SecE1xhlen2, DYa, DYb, (Endb-Enda)/2);
554
555 GeoLogVol *lvse1g1 = new GeoLogVol("LAr::DM::SectorEnvelopes1g1",Trdair1,matLArServices1);
556 GeoPhysVol *sectorenvelopes1g1 = new GeoPhysVol(lvse1g1);
557
558 GeoLogVol *lvse1g2 = new GeoLogVol("LAr::DM::SectorEnvelopes1g3",Trdair1,matLArServices2);
559 GeoPhysVol *sectorenvelopes1g2 = new GeoPhysVol(lvse1g2);
560
561 GeoLogVol *lvse1g3 = new GeoLogVol("LAr::DM::SectorEnvelopes1g3",Trdair1,matLArServices3);
562 GeoPhysVol *sectorenvelopes1g3 = new GeoPhysVol(lvse1g3);
563
564 GeoLogVol *lvse1g4 = new GeoLogVol("LAr::DM::SectorEnvelopes1g4",Trdair1,matLArServices4);
565 GeoPhysVol *sectorenvelopes1g4 = new GeoPhysVol(lvse1g4);
566
567 GeoLogVol *lvse1g5 = new GeoLogVol("LAr::DM::SectorEnvelopes1g5",Trdair1,matLArServices5);
568 GeoPhysVol *sectorenvelopes1g5 = new GeoPhysVol(lvse1g5);
569
570 GeoLogVol *lvse1g6 = new GeoLogVol("LAr::DM::SectorEnvelopes1g6",Trdair1,matLArServices6);
571 GeoPhysVol *sectorenvelopes1g6 = new GeoPhysVol(lvse1g6);
572
573 GeoLogVol *lvse1g7 = new GeoLogVol("LAr::DM::SectorEnvelopes1g7",Trdair1,matLArServices7);
574 GeoPhysVol *sectorenvelopes1g7 = new GeoPhysVol(lvse1g7);
575
576 GeoLogVol *lvse1g8 = new GeoLogVol("LAr::DM::SectorEnvelopes1g8",Trdair1,matLArServices8);
577 GeoPhysVol *sectorenvelopes1g8 = new GeoPhysVol(lvse1g8);
578
579 GeoLogVol *lvse1g9 = new GeoLogVol("LAr::DM::SectorEnvelopes1g9",Trdair1,matLArServices9);
580 GeoPhysVol *sectorenvelopes1g9 = new GeoPhysVol(lvse1g9);
581
582 GeoLogVol *lvse1g10 = new GeoLogVol("LAr::DM::SectorEnvelopes1g10",Trdair1,matLArServices10);
583 GeoPhysVol *sectorenvelopes1g10 = new GeoPhysVol(lvse1g10);
584
585 GeoLogVol *lvse1g11 = new GeoLogVol("LAr::DM::SectorEnvelopes1g11",Trdair1,matLArServices11);
586 GeoPhysVol *sectorenvelopes1g11 = new GeoPhysVol(lvse1g11);
587
588 GeoLogVol *lvse1g12 = new GeoLogVol("LAr::DM::SectorEnvelopes1g12",Trdair1,matLArServices12);
589 GeoPhysVol *sectorenvelopes1g12 = new GeoPhysVol(lvse1g12);
590
591 GeoLogVol *lvse1g13 = new GeoLogVol("LAr::DM::SectorEnvelopes1g13",Trdair1,matLArServices13);
592 GeoPhysVol *sectorenvelopes1g13 = new GeoPhysVol(lvse1g13);
593
594 GeoLogVol *lvse1g14 = new GeoLogVol("LAr::DM::SectorEnvelopes1g14",Trdair1,matLArServices14);
595 GeoPhysVol *sectorenvelopes1g14 = new GeoPhysVol(lvse1g14);
596
597 GeoLogVol *lvse1g15 = new GeoLogVol("LAr::DM::SectorEnvelopes1g15",Trdair1,matLArServices15);
598 GeoPhysVol *sectorenvelopes1g15 = new GeoPhysVol(lvse1g15);
599
600 GeoLogVol *lvse1g16 = new GeoLogVol("LAr::DM::SectorEnvelopes1g16",Trdair1,matLArServices16);
601 GeoPhysVol *sectorenvelopes1g16 = new GeoPhysVol(lvse1g16);
602
603
604 std::vector<GeoPhysVol*> se1List;
605 se1List.push_back(sectorenvelopes1g1);
606 se1List.push_back(sectorenvelopes1g2);
607 se1List.push_back(sectorenvelopes1g3);
608 se1List.push_back(sectorenvelopes1g4);
609 se1List.push_back(sectorenvelopes1g5);
610 se1List.push_back(sectorenvelopes1g6);
611 se1List.push_back(sectorenvelopes1g7);
612 se1List.push_back(sectorenvelopes1g8);
613 se1List.push_back(sectorenvelopes1g9);
614 se1List.push_back(sectorenvelopes1g10);
615 se1List.push_back(sectorenvelopes1g11);
616 se1List.push_back(sectorenvelopes1g12);
617 se1List.push_back(sectorenvelopes1g13);
618 se1List.push_back(sectorenvelopes1g14);
619 se1List.push_back(sectorenvelopes1g15);
620 se1List.push_back(sectorenvelopes1g16);
621
622 // dedicated volumes to add in SectorEnvelopes1
623 // (a) volume in Trd shape to add slice of extra material in a given eta bin
624 IRDBRecordset_ptr BarrelDMRing = rdbAccess->getRecordsetPtr("LArDMEnv1Ring",keyLAr.tag(),keyLAr.node());
625 for (unsigned int i=0;i<BarrelDMRing->size();i++) {
626 double etaMin=(*BarrelDMRing)[i]->getDouble("ETAMIN");
627 double etaMax=(*BarrelDMRing)[i]->getDouble("ETAMAX");
628 double thicknessExtra=(*BarrelDMRing)[i]->getDouble("THICKNESS");
629 std::string ringName = "LAr::DM::SectorEnvelopes1::"+(*BarrelDMRing)[i]->getString("RINGNAME");
630
631 double radiusMin=SecE1ztr/sinh(etaMax);
632 double radiusMax=SecE1ztr/sinh(etaMin);
633 double dy1 = DYa + (DYb-DYa)*(radiusMin-Enda)*inv_Endab;
634 double dy2 = DYa + (DYb-DYa)*(radiusMax-Enda)*inv_Endab;
635 double zpos=0.5*(radiusMax+radiusMin - (Endb+Enda));
636 const GeoMaterial *matExtraTdr = materialManager->getMaterial((*BarrelDMRing)[i]->getString("MATERIAL"));
637
638 GeoTrd *extraMatTdr = new GeoTrd(thicknessExtra/2., thicknessExtra/2., dy1, dy2, (radiusMax-radiusMin)/2);
639 GeoLogVol *extraMatLog = new GeoLogVol(ringName,extraMatTdr,matExtraTdr);
640 GeoPhysVol *extraMatPhys = new GeoPhysVol(extraMatLog);
641 for (unsigned int isect=0;isect<se1List.size();isect++) {
642 se1List[isect]->add(new GeoTransform(GeoTrf::TranslateZ3D(zpos)));
643 se1List[isect]->add(extraMatPhys);
644 }
645 }
646
647 // (b) extra material at fixed phi locations in the PPF1 area
648 IRDBRecordset_ptr BarrelDMPhiBox = rdbAccess->getRecordsetPtr("LArDMEnv1PhiBox",keyLAr.tag(),keyLAr.node());
649 for (unsigned int i=0;i<BarrelDMPhiBox->size();i++) {
650 double eta=(*BarrelDMPhiBox)[i]->getDouble("ETA");
651 double phi0=(*BarrelDMPhiBox)[i]->getDouble("PHI0");
652 double deltaR=(*BarrelDMPhiBox)[i]->getDouble("DELTAR");
653 double deltaRphi=(*BarrelDMPhiBox)[i]->getDouble("DELTARPHI");
654 double thickness=(*BarrelDMPhiBox)[i]->getDouble("THICKNESS");
655 int nphi=(*BarrelDMPhiBox)[i]->getInt("NPHI");
656 int noHorizontal = (*BarrelDMPhiBox)[i]->getInt("NOHORIZ");
657
658 const GeoMaterial* matExtraPPF1 = materialManager->getMaterial((*BarrelDMPhiBox)[i]->getString("MATERIAL"));
659 std::string boxName = "LAr::DM::SectorEnvelopes1::"+(*BarrelDMPhiBox)[i]->getString("BOXNAME");
660
661 GeoBox *ppf1Box = new GeoBox(thickness/2.,deltaRphi/2.,deltaR/2.);
662 GeoLogVol *ppf1Log = new GeoLogVol(boxName,ppf1Box,matExtraPPF1);
663 GeoPhysVol *ppf1Phys = new GeoPhysVol(ppf1Log);
664
665 int nPerEnv1 = nphi/NCrates;
666 double dphi=2.*M_PI/((float)(nphi));
667 double radius=SecE1ztr/sinh(eta);
668 for (int iphi=0;iphi<nPerEnv1;iphi++) {
669 double xpos=0.;
670 double ypos=radius*sin(phi0+((float)(iphi))*dphi);
671 double zpos=radius*cos(phi0+((float)(iphi))*dphi) - (Endb+Enda)/2.;
672 for (unsigned int isect=0;isect<se1List.size();isect++) {
673 // no PPF1 box around phi=0 and phi=pi
674 if (noHorizontal>0 && ((isect==7 && iphi==1) || (isect==8 && iphi==0) || (isect==15 && iphi==1) || (isect==0 && iphi==0) ) ) continue;
675 se1List[isect]->add(new GeoTransform(GeoTrf::Translate3D(xpos,ypos,zpos)));
676 se1List[isect]->add(ppf1Phys);
677 }
678 }
679 }
680
681 // ----- build base plates -----
682 GeoTrd *Trd1alu = new GeoTrd(BasePxhlen1, BasePxhlen2, BasePyhlen1, BasePyhlen2, BasePzhlen);
683 GeoLogVol *lvbp = new GeoLogVol("LAr::DM::BasePlates",Trd1alu,alu);
684 GeoPhysVol *baseplates = new GeoPhysVol(lvbp);
685
686 // ----- build bridge plates -----
687 GeoTrap *Trapalu = new GeoTrap(BridgePzhlen, BridgePtheta*Gaudi::Units::deg, BridgePphi, BridgePyzn, BridgePxynzn, BridgePxypzn, BridgePangn, BridgePyzp, BridgePxynzp, BridgePxypzp, BridgePangp);
688 GeoLogVol *lvbrp = new GeoLogVol("LAr::DM::BridgePlates",Trapalu,alu);
689 GeoPhysVol *bridgeplates = new GeoPhysVol(lvbrp);
690
691
692 // ----- build sector plates -----
693 GeoTrd *Trd2alu = new GeoTrd(SecPxhlen1, SecPxhlen2, SecPyhlen1, SecPyhlen2, SecPzhlen );
694 const GeoShape & SectorPlates= ((*Trd2alu).
695 subtract((*Box) <<GeoTrf::Transform3D(Cut3Boxp)).
696 subtract((*Box) <<GeoTrf::Transform3D(Cut4Boxp)));
697 GeoLogVol *lvsp = new GeoLogVol("LAr::DM::SectorPlates",&SectorPlates,alu);
698 GeoPhysVol *sectorplates = new GeoPhysVol(lvsp);
699
700
701 //-------------- Place volumes in LAr Envelope -------------------
702
703 //sectorPlates
704 TRANSFUNCTION spA = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(SecPxtr)*GeoTrf::TranslateZ3D(SecPztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
705 TRANSFUNCTION spC = Pow(GeoTrf::RotateZ3D(1.0),f+(Alfa/2))*GeoTrf::TranslateX3D(SecPxtr)*GeoTrf::TranslateZ3D(-SecPztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
706 GeoSerialTransformer *sptA = new GeoSerialTransformer(sectorplates,&spA, NCrates);
707 GeoSerialTransformer *sptC = new GeoSerialTransformer(sectorplates,&spC, NCrates);
708 envelope->add(sptA);
709 envelope->add(sptC);
710
711 //bridgePlates
712 TRANSFUNCTION brpA1 = Pow(GeoTrf::RotateZ3D(1.0),f-(5*Alfa/2))*GeoTrf::TranslateX3D(BridgePxtr)*GeoTrf::TranslateZ3D(BridgePztr)*GeoTrf::RotateZ3D(90*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)*GeoTrf::RotateX3D(90*Gaudi::Units::deg);
713 TRANSFUNCTION brpA2 = Pow(GeoTrf::RotateZ3D(1.0),f+(13*Alfa/2))*GeoTrf::TranslateX3D(BridgePxtr)*GeoTrf::TranslateZ3D(BridgePztr)*GeoTrf::RotateZ3D(90*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)*GeoTrf::RotateX3D(90*Gaudi::Units::deg);
714 TRANSFUNCTION brpC1 = Pow(GeoTrf::RotateZ3D(1.0),f-(5*Alfa/2))*GeoTrf::TranslateX3D(BridgePxtr)*GeoTrf::TranslateZ3D(-BridgePztr)*GeoTrf::RotateZ3D(-90*Gaudi::Units::deg)*GeoTrf::RotateY3D(-90*Gaudi::Units::deg)*GeoTrf::RotateX3D(-90*Gaudi::Units::deg);
715 TRANSFUNCTION brpC2 = Pow(GeoTrf::RotateZ3D(1.0),f+(13*Alfa/2))*GeoTrf::TranslateX3D(BridgePxtr)*GeoTrf::TranslateZ3D(-BridgePztr)*GeoTrf::RotateZ3D(-90*Gaudi::Units::deg)*GeoTrf::RotateY3D(-90*Gaudi::Units::deg)*GeoTrf::RotateX3D(-90*Gaudi::Units::deg); GeoSerialTransformer *brptA1 = new GeoSerialTransformer(bridgeplates,&brpA1, 5);
716 GeoSerialTransformer *brptA2 = new GeoSerialTransformer(bridgeplates,&brpA2, 5);
717 GeoSerialTransformer *brptC1 = new GeoSerialTransformer(bridgeplates,&brpC1, 5);
718 GeoSerialTransformer *brptC2 = new GeoSerialTransformer(bridgeplates,&brpC2, 5);
719 envelope->add(brptA1);
720 envelope->add(brptA2);
721 envelope->add(brptC1);
722 envelope->add(brptC2);
723
724 //basePlates
725 TRANSFUNCTION bpA = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(BasePxtr)*GeoTrf::TranslateZ3D(BasePztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
726 TRANSFUNCTION bpC = Pow(GeoTrf::RotateZ3D(1.0),f+(Alfa/2))*GeoTrf::TranslateX3D(BasePxtr)*GeoTrf::TranslateZ3D(-BasePztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
727 GeoSerialTransformer *bptA = new GeoSerialTransformer(baseplates,&bpA, NCrates);
728 GeoSerialTransformer *bptC = new GeoSerialTransformer(baseplates,&bpC, NCrates);
729 envelope->add(bptA);
730 envelope->add(bptC);
731
732 //sectorEnvelopes1
733 //counter-clockwise from top if taking sideA for reference (clockwise for sideC)
734 TRANSFUNCTION seA1G5 = Pow(GeoTrf::RotateZ3D(1.0),f+(9*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
735 TRANSFUNCTION seC1G5 = Pow(GeoTrf::RotateZ3D(1.0),f+(9*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
736 TRANSFUNCTION seA1G6 = Pow(GeoTrf::RotateZ3D(1.0),f+(11*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
737 TRANSFUNCTION seC1G6 = Pow(GeoTrf::RotateZ3D(1.0),f+(11*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
738 TRANSFUNCTION seA1G7 = Pow(GeoTrf::RotateZ3D(1.0),f+(13*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
739 TRANSFUNCTION seC1G7 = Pow(GeoTrf::RotateZ3D(1.0),f+(13*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
740 TRANSFUNCTION seA1G8 = Pow(GeoTrf::RotateZ3D(1.0),f+(15*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
741 TRANSFUNCTION seC1G8 = Pow(GeoTrf::RotateZ3D(1.0),f+(15*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
742 TRANSFUNCTION seA1G9 = Pow(GeoTrf::RotateZ3D(1.0),f+(17*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
743 TRANSFUNCTION seC1G9 = Pow(GeoTrf::RotateZ3D(1.0),f+(17*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
744 TRANSFUNCTION seA1G10 = Pow(GeoTrf::RotateZ3D(1.0),f+(19*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
745 TRANSFUNCTION seC1G10 = Pow(GeoTrf::RotateZ3D(1.0),f+(19*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
746 TRANSFUNCTION seA1G11 = Pow(GeoTrf::RotateZ3D(1.0),f+(21*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
747 TRANSFUNCTION seC1G11 = Pow(GeoTrf::RotateZ3D(1.0),f+(21*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
748 TRANSFUNCTION seA1G12 = Pow(GeoTrf::RotateZ3D(1.0),f+(23*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
749 TRANSFUNCTION seC1G12 = Pow(GeoTrf::RotateZ3D(1.0),f+(23*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
750 //clockwise from top if taking sideA for reference (counter-clockwise for sideC)
751 TRANSFUNCTION seA1G4 = Pow(GeoTrf::RotateZ3D(1.0),f+(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
752 TRANSFUNCTION seC1G4 = Pow(GeoTrf::RotateZ3D(1.0),f+(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
753 TRANSFUNCTION seA1G3 = Pow(GeoTrf::RotateZ3D(1.0),f+(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
754 TRANSFUNCTION seC1G3 = Pow(GeoTrf::RotateZ3D(1.0),f+(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
755 TRANSFUNCTION seA1G2 = Pow(GeoTrf::RotateZ3D(1.0),f+(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
756 TRANSFUNCTION seC1G2 = Pow(GeoTrf::RotateZ3D(1.0),f+(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
757 TRANSFUNCTION seA1G1 = Pow(GeoTrf::RotateZ3D(1.0),f+(1*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
758 TRANSFUNCTION seC1G1 = Pow(GeoTrf::RotateZ3D(1.0),f+(1*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
759 TRANSFUNCTION seA1G16 = Pow(GeoTrf::RotateZ3D(1.0),f-(1*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
760 TRANSFUNCTION seC1G16 = Pow(GeoTrf::RotateZ3D(1.0),f-(1*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
761 TRANSFUNCTION seA1G15 = Pow(GeoTrf::RotateZ3D(1.0),f-(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
762 TRANSFUNCTION seC1G15 = Pow(GeoTrf::RotateZ3D(1.0),f-(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
763 TRANSFUNCTION seA1G14 = Pow(GeoTrf::RotateZ3D(1.0),f-(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
764 TRANSFUNCTION seC1G14 = Pow(GeoTrf::RotateZ3D(1.0),f-(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
765 TRANSFUNCTION seA1G13 = Pow(GeoTrf::RotateZ3D(1.0),f-(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
766 TRANSFUNCTION seC1G13 = Pow(GeoTrf::RotateZ3D(1.0),f-(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
767
768 GeoSerialTransformer *setA1G5 = new GeoSerialTransformer(sectorenvelopes1g5,&seA1G5, 1);
769 GeoSerialTransformer *setC1G5 = new GeoSerialTransformer(sectorenvelopes1g5,&seC1G5, 1);
770 GeoSerialTransformer *setA1G6 = new GeoSerialTransformer(sectorenvelopes1g6,&seA1G6, 1);
771 GeoSerialTransformer *setC1G6 = new GeoSerialTransformer(sectorenvelopes1g6,&seC1G6, 1);
772 GeoSerialTransformer *setA1G7 = new GeoSerialTransformer(sectorenvelopes1g7,&seA1G7, 1);
773 GeoSerialTransformer *setC1G7 = new GeoSerialTransformer(sectorenvelopes1g7,&seC1G7, 1);
774 GeoSerialTransformer *setA1G8 = new GeoSerialTransformer(sectorenvelopes1g8,&seA1G8, 1);
775 GeoSerialTransformer *setC1G8 = new GeoSerialTransformer(sectorenvelopes1g8,&seC1G8, 1);
776 GeoSerialTransformer *setA1G9 = new GeoSerialTransformer(sectorenvelopes1g9,&seA1G9, 1);
777 GeoSerialTransformer *setC1G9 = new GeoSerialTransformer(sectorenvelopes1g9,&seC1G9, 1);
778 GeoSerialTransformer *setA1G10 = new GeoSerialTransformer(sectorenvelopes1g10,&seA1G10, 1);
779 GeoSerialTransformer *setC1G10 = new GeoSerialTransformer(sectorenvelopes1g10,&seC1G10, 1);
780 GeoSerialTransformer *setA1G11 = new GeoSerialTransformer(sectorenvelopes1g11,&seA1G11, 1);
781 GeoSerialTransformer *setC1G11 = new GeoSerialTransformer(sectorenvelopes1g11,&seC1G11, 1);
782 GeoSerialTransformer *setA1G12 = new GeoSerialTransformer(sectorenvelopes1g12,&seA1G12, 1);
783 GeoSerialTransformer *setC1G12 = new GeoSerialTransformer(sectorenvelopes1g12,&seC1G12, 1);
784
785 GeoSerialTransformer *setA1G4 = new GeoSerialTransformer(sectorenvelopes1g4,&seA1G4, 1);
786 GeoSerialTransformer *setC1G4 = new GeoSerialTransformer(sectorenvelopes1g4,&seC1G4, 1);
787 GeoSerialTransformer *setA1G3 = new GeoSerialTransformer(sectorenvelopes1g3,&seA1G3, 1);
788 GeoSerialTransformer *setC1G3 = new GeoSerialTransformer(sectorenvelopes1g3,&seC1G3, 1);
789 GeoSerialTransformer *setA1G2 = new GeoSerialTransformer(sectorenvelopes1g2,&seA1G2, 1);
790 GeoSerialTransformer *setC1G2 = new GeoSerialTransformer(sectorenvelopes1g2,&seC1G2, 1);
791 GeoSerialTransformer *setA1G1 = new GeoSerialTransformer(sectorenvelopes1g1,&seA1G1, 1);
792 GeoSerialTransformer *setC1G1 = new GeoSerialTransformer(sectorenvelopes1g1,&seC1G1, 1);
793 GeoSerialTransformer *setA1G16 = new GeoSerialTransformer(sectorenvelopes1g16,&seA1G16, 1);
794 GeoSerialTransformer *setC1G16 = new GeoSerialTransformer(sectorenvelopes1g16,&seC1G16, 1);
795 GeoSerialTransformer *setA1G15 = new GeoSerialTransformer(sectorenvelopes1g15,&seA1G15, 1);
796 GeoSerialTransformer *setC1G15 = new GeoSerialTransformer(sectorenvelopes1g15,&seC1G15, 1);
797 GeoSerialTransformer *setA1G14 = new GeoSerialTransformer(sectorenvelopes1g14,&seA1G14, 1);
798 GeoSerialTransformer *setC1G14 = new GeoSerialTransformer(sectorenvelopes1g14,&seC1G14, 1);
799 GeoSerialTransformer *setA1G13 = new GeoSerialTransformer(sectorenvelopes1g13,&seA1G13, 1);
800 GeoSerialTransformer *setC1G13 = new GeoSerialTransformer(sectorenvelopes1g13,&seC1G13, 1);
801
802 envelope->add(setA1G5);
803 envelope->add(setC1G5);
804 envelope->add(setA1G6);
805 envelope->add(setC1G6);
806 envelope->add(setA1G7);
807 envelope->add(setC1G7);
808 envelope->add(setA1G8);
809 envelope->add(setC1G8);
810 envelope->add(setA1G9);
811 envelope->add(setC1G9);
812 envelope->add(setA1G10);
813 envelope->add(setC1G10);
814 envelope->add(setA1G11);
815 envelope->add(setC1G11);
816 envelope->add(setA1G12);
817 envelope->add(setC1G12);
818
819 envelope->add(setA1G4);
820 envelope->add(setC1G4);
821 envelope->add(setA1G3);
822 envelope->add(setC1G3);
823 envelope->add(setA1G2);
824 envelope->add(setC1G2);
825 envelope->add(setA1G1);
826 envelope->add(setC1G1);
827 envelope->add(setA1G16);
828 envelope->add(setC1G16);
829 envelope->add(setA1G15);
830 envelope->add(setC1G15);
831 envelope->add(setA1G14);
832 envelope->add(setC1G14);
833 envelope->add(setA1G13);
834 envelope->add(setC1G13);
835
836 createSectorEnvelopes2FromDB (envelope, materialManager,
837 trdMap, *BarrelDMTrds,
838 trapMap, *BarrelDMTraps,
839 boxMap, *BarrelDMBoxes,
840 f, Box);
841 createBridgeEnvelopesFromDB (envelope, trapMap, *BarrelDMTraps,
842 matLArServices8, f);
843 createBaseEnvelopesFromDB (envelope, trdMap, *BarrelDMTrds,
844 matLArServices8, f);
845}
M_PI
#define M_PI
Definition ActiveFraction.h:11
Box
Acts::Volume::BoundingBox Box
Definition ActsCaloTrackingVolumeBuilder.cxx:28
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
deltaR
Scalar deltaR(const MatrixBase< Derived > &vec) const
Definition AmgMatrixBasePlugin.h:105
Enda
static const double Enda
Definition CrackDMConstruction.cxx:53
DYa
static const double DYa
Definition CrackDMConstruction.cxx:56
createBridgeEnvelopesFromDB
void createBridgeEnvelopesFromDB(GeoFullPhysVol *envelope, std::map< std::string, unsigned int > &trapMap, IRDBRecordset &BarrelDMTraps, const GeoMaterial *matLArServices8, GENFUNCTION &f)
Definition CrackDMConstruction.cxx:270
createBaseEnvelopesFromDB
void createBaseEnvelopesFromDB(GeoFullPhysVol *envelope, std::map< std::string, unsigned int > &trdMap, IRDBRecordset &BarrelDMTrds, const GeoMaterial *matLArServices8, GENFUNCTION &f)
Definition CrackDMConstruction.cxx:306
Endb
static const double Endb
Definition CrackDMConstruction.cxx:54
DYb
static const double DYb
Definition CrackDMConstruction.cxx:57
createSectorEnvelopes2FromDB
void createSectorEnvelopes2FromDB(GeoFullPhysVol *envelope, StoredMaterialManager *materialManager, std::map< std::string, unsigned int > &trdMap, IRDBRecordset &BarrelDMTrds, std::map< std::string, unsigned int > &trapMap, IRDBRecordset &BarrelDMTraps, std::map< std::string, unsigned int > &boxMap, IRDBRecordset &BarrelDMBoxes, GENFUNCTION &f, GeoBox *Box)
Definition CrackDMConstruction.cxx:62
IRDBRecordset_ptr
std::shared_ptr< IRDBRecordset > IRDBRecordset_ptr
Definition IRDBAccessSvc.h:27
DecodeVersionKey
This is a helper class to query the version tags from GeoModelSvc and determine the appropriate tag a...
Definition DecodeVersionKey.h:18
IRDBAccessSvc::getRecordsetPtr
virtual IRDBRecordset_ptr getRecordsetPtr(const std::string &node, const std::string &tag, const std::string &tag2node="", const std::string &connName="ATLASDD")=0
Provides access to the Recordset object containing HVS-tagged data.
IRDBRecordset::size
virtual unsigned int size() const =0
StoredMaterialManager::getMaterial
virtual const GeoMaterial * getMaterial(const std::string &name)=0
Variable
Wrapper around a histogram which allows for some additional filling patterns and data manipulation.
Definition Trigger/TrigCost/TrigCostAnalysis/src/Variable.h:39
ParticleGun_SamplingFraction.radius
radius
Definition ParticleGun_SamplingFraction.py:96
PrintTrkAnaSummary.l
l
Printing final latex table to .tex output file.
Definition PrintTrkAnaSummary.py:370
hist_file_dump.f
f
Definition hist_file_dump.py:140

◆ createSectorEnvelopes2FromDB()

void createSectorEnvelopes2FromDB ( GeoFullPhysVol * envelope,
StoredMaterialManager * materialManager,
std::map< std::string, unsigned int > & trdMap,
IRDBRecordset & BarrelDMTrds,
std::map< std::string, unsigned int > & trapMap,
IRDBRecordset & BarrelDMTraps,
std::map< std::string, unsigned int > & boxMap,
IRDBRecordset & BarrelDMBoxes,
GENFUNCTION & f,
GeoBox * Box )

Definition at line 62 of file CrackDMConstruction.cxx.

72{
73 unsigned int recordIndex = trdMap["SecE2"];
74 double SecE2xhlen1 = BarrelDMTrds[recordIndex]->getDouble("XHLEN1");
75 double SecE2xhlen2 = BarrelDMTrds[recordIndex]->getDouble("XHLEN2");
76 double SecE2ztr = BarrelDMTrds[recordIndex]->getDouble("ZTR");
77
78 recordIndex = trdMap["Spb0"];
79 double Spb0xhlen1 = BarrelDMTrds[recordIndex]->getDouble("XHLEN1");
80 double Spb0xhlen2 = BarrelDMTrds[recordIndex]->getDouble("XHLEN2");
81 double Spb0yhlen1 = BarrelDMTrds[recordIndex]->getDouble("YHLEN1");
82 double Spb0yhlen2 = BarrelDMTrds[recordIndex]->getDouble("YHLEN2");
83 double Spb0zhlen = BarrelDMTrds[recordIndex]->getDouble("ZHLEN");
84 double Spb0ytr = BarrelDMTrds[recordIndex]->getDouble("YTR");
85 double Spb0xrot = BarrelDMTrds[recordIndex]->getDouble("XROT");
86
87 recordIndex = boxMap["Box"];
88 //double Boxhlen = BarrelDMBoxes[recordIndex]->getDouble("HLEN");
89 //double Boxhwdt = BarrelDMBoxes[recordIndex]->getDouble("HWDT");
90 //double Boxhhgt = BarrelDMBoxes[recordIndex]->getDouble("HHGT");
91 double Boxxtr = BarrelDMBoxes[recordIndex]->getDouble("XTR");
92 double Boxytr = BarrelDMBoxes[recordIndex]->getDouble("YTR");
93 double Boxztr = BarrelDMBoxes[recordIndex]->getDouble("ZTR");
94 //double Boxxrot = BarrelDMBoxes[recordIndex]->getDouble("XROT");
95
96 recordIndex = boxMap["SplBox"];
97 double SplBoxhlen = BarrelDMBoxes[recordIndex]->getDouble("HLEN");
98 double SplBoxhwdt = BarrelDMBoxes[recordIndex]->getDouble("HWDT");
99 double SplBoxhhgt = BarrelDMBoxes[recordIndex]->getDouble("HHGT");
100 double SplBoxytr = BarrelDMBoxes[recordIndex]->getDouble("YTR");
101 double SplBoxztr = BarrelDMBoxes[recordIndex]->getDouble("ZTR");
102 double SplBoxxrot = BarrelDMBoxes[recordIndex]->getDouble("XROT");
103
104 recordIndex = trapMap["Spb1"];
105 double Spb1zhlen = BarrelDMTraps[recordIndex]->getDouble("ZHLEN");
106 double Spb1theta = BarrelDMTraps[recordIndex]->getDouble("THETA");
107 double Spb1phi = BarrelDMTraps[recordIndex]->getDouble("PHI");
108 double Spb1yzn = BarrelDMTraps[recordIndex]->getDouble("YZN");
109 double Spb1xynzn = BarrelDMTraps[recordIndex]->getDouble("XYNZN");
110 double Spb1xypzn = BarrelDMTraps[recordIndex]->getDouble("XYPZN");
111 double Spb1angn = BarrelDMTraps[recordIndex]->getDouble("ANGN");
112 double Spb1yzp = BarrelDMTraps[recordIndex]->getDouble("YZP");
113 double Spb1xynzp = BarrelDMTraps[recordIndex]->getDouble("XYNZP");
114 double Spb1xypzp = BarrelDMTraps[recordIndex]->getDouble("XYPZP");
115 double Spb1angp = BarrelDMTraps[recordIndex]->getDouble("ANGP");
116 double Spb1xtr = BarrelDMTraps[recordIndex]->getDouble("XTR");
117 double Spb1ytr = BarrelDMTraps[recordIndex]->getDouble("YTR");
118 double Spb1ztr = BarrelDMTraps[recordIndex]->getDouble("ZTR");
119 double Spb1xrot = BarrelDMTraps[recordIndex]->getDouble("XROT");
120
121 recordIndex = trapMap["Spb3"];
122 double Spb3zhlen = BarrelDMTraps[recordIndex]->getDouble("ZHLEN");
123 double Spb3theta = BarrelDMTraps[recordIndex]->getDouble("THETA");
124 double Spb3phi = BarrelDMTraps[recordIndex]->getDouble("PHI");
125 double Spb3yzn = BarrelDMTraps[recordIndex]->getDouble("YZN");
126 double Spb3xynzn = BarrelDMTraps[recordIndex]->getDouble("XYNZN");
127 double Spb3xypzn = BarrelDMTraps[recordIndex]->getDouble("XYPZN");
128 double Spb3angn = BarrelDMTraps[recordIndex]->getDouble("ANGN");
129 double Spb3yzp = BarrelDMTraps[recordIndex]->getDouble("YZP");
130 double Spb3xynzp = BarrelDMTraps[recordIndex]->getDouble("XYNZP");
131 double Spb3xypzp = BarrelDMTraps[recordIndex]->getDouble("XYPZP");
132 double Spb3angp = BarrelDMTraps[recordIndex]->getDouble("ANGP");
133 double Spb3ztr = BarrelDMTraps[recordIndex]->getDouble("ZTR");
134
135 recordIndex = trapMap["Spb2"];
136 double Spb2zhlen = BarrelDMTraps[recordIndex]->getDouble("ZHLEN");
137 double Spb2theta = BarrelDMTraps[recordIndex]->getDouble("THETA");
138 double Spb2phi = BarrelDMTraps[recordIndex]->getDouble("PHI");
139 double Spb2yzn = BarrelDMTraps[recordIndex]->getDouble("YZN");
140 double Spb2xynzn = BarrelDMTraps[recordIndex]->getDouble("XYNZN");
141 double Spb2xypzn = BarrelDMTraps[recordIndex]->getDouble("XYPZN");
142 double Spb2angn = BarrelDMTraps[recordIndex]->getDouble("ANGN");
143 double Spb2yzp = BarrelDMTraps[recordIndex]->getDouble("YZP");
144 double Spb2xynzp = BarrelDMTraps[recordIndex]->getDouble("XYNZP");
145 double Spb2xypzp = BarrelDMTraps[recordIndex]->getDouble("XYPZP");
146 double Spb2angp = BarrelDMTraps[recordIndex]->getDouble("ANGP");
147 double Spb2ytr = BarrelDMTraps[recordIndex]->getDouble("YTR");
148 double Spb2ztr = BarrelDMTraps[recordIndex]->getDouble("ZTR");
149
150 const GeoMaterial* matLArServices17 = materialManager->getMaterial("LAr::LArServices17");// 0.035*gram/cm3
151 const GeoMaterial* matLArServices18 = materialManager->getMaterial("LAr::LArServices18");// 0.240*gram/cm3
152 const GeoMaterial* matLArServices19 = materialManager->getMaterial("LAr::LArServices19");// 0.469*gram/cm3
153 const GeoMaterial* matLArServices20 = materialManager->getMaterial("LAr::LArServices20");// 0.353*gram/cm3
154 const GeoMaterial *alu = materialManager->getMaterial("std::Aluminium"); //2.7 g/cm3
155 const GeoMaterial *air = materialManager->getMaterial("std::Air"); //0.001214 g/cm3
156
157 GeoTrf::Transform3D Cut3Boxe = GeoTrf::Translate3D(Boxxtr, Boxytr, Boxztr)*GeoTrf::RotateX3D(-20*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
158 GeoTrf::Transform3D Cut4Boxe = GeoTrf::Translate3D(Boxxtr, -Boxytr,Boxztr)*GeoTrf::RotateX3D(20*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
159
160 // build 5 instances of SectorEnvelopes1 with 3 different materials!
161 GeoTrd *Trdair2 = new GeoTrd(SecE2xhlen1, SecE2xhlen2, DYb, DYc, (Endc-Endb)/2);
162 const GeoShape & SectorEnvelope= ((*Trdair2).
163 subtract((*Box) <<GeoTrf::Transform3D(Cut3Boxe)).
164 subtract((*Box) <<GeoTrf::Transform3D(Cut4Boxe)));
165
166 const GeoShape & SectorEnvelopes= ((SectorEnvelope).
167 add(SectorEnvelope << GeoTrf::TranslateY3D(-(DYb+DYc)*cos(Alfa/2)*cos(Alfa/2))*GeoTrf::TranslateZ3D(-(DYb+DYc)*0.5*sin(Alfa))*GeoTrf::RotateX3D(Alfa)));
168
169 GeoLogVol *lvse2r = new GeoLogVol("LAr::DM::SectorEnvelopes2r",&SectorEnvelopes,matLArServices20);
170 GeoPhysVol *sectorenvelopes2r = new GeoPhysVol(lvse2r); // for right-handed splice boxes
171
172 GeoLogVol *lvse2l = new GeoLogVol("LAr::DM::SectorEnvelopes2l",&SectorEnvelopes,matLArServices20);
173 GeoPhysVol *sectorenvelopes2l = new GeoPhysVol(lvse2l); // for left-handed splice boxes
174
175 GeoLogVol *lvse2h = new GeoLogVol("LAr::DM::SectorEnvelopes2h",&SectorEnvelopes,matLArServices19);
176 GeoPhysVol *sectorenvelopes2h = new GeoPhysVol(lvse2h); // no splice boxes horizontal at 0 & 180 Gaudi::Units::deg.
177
178 GeoLogVol *lvse2vup = new GeoLogVol("LAr::DM::SectorEnvelopes2vup",&SectorEnvelopes,matLArServices17);
179 GeoPhysVol *sectorenvelopes2vup = new GeoPhysVol(lvse2vup); // no splice boxes vertical up at 90 Gaudi::Units::deg
180
181 GeoLogVol *lvse2vd = new GeoLogVol("LAr::DM::SectorEnvelopes2Vd",&SectorEnvelopes,matLArServices18);
182 GeoPhysVol *sectorenvelopes2vd = new GeoPhysVol(lvse2vd); // no splice boxes vertical down at 270 Gaudi::Units::deg
183
184 //---------- Build Splice boxes for InDet optical fibers--------
185
186 GeoTrap *GeoTrap1 = new GeoTrap(Spb1zhlen, Spb1theta, Spb1phi, Spb1yzn, Spb1xynzn, Spb1xypzn, Spb1angn, Spb1yzp, Spb1xynzp, Spb1xypzp, Spb1angp);
187 GeoBox *Box1 = new GeoBox(SplBoxhlen, SplBoxhwdt, SplBoxhhgt);
188 const GeoShape & SpliceBox = ((*GeoTrap1).
189 subtract(*Box1 << GeoTrf::TranslateZ3D(SplBoxztr)*GeoTrf::TranslateY3D(-SplBoxytr)*GeoTrf::RotateX3D(SplBoxxrot*Gaudi::Units::deg)));
190
191 GeoTransform *xtr = new GeoTransform (GeoTrf::TranslateZ3D(Spb1ztr)*GeoTrf::TranslateY3D(-Spb1ytr)*GeoTrf::TranslateX3D(Spb1xtr)*GeoTrf::RotateX3D(Spb1xrot*Gaudi::Units::deg));
192 sectorenvelopes2r->add(xtr);
193 GeoLogVol *lvspbr = new GeoLogVol("LAr::DM::SPliceBoxr",&SpliceBox,alu);
194 GeoPhysVol *spliceboxr = new GeoPhysVol(lvspbr);
195 sectorenvelopes2r->add(spliceboxr);
196
197 GeoTransform *xtl = new GeoTransform (GeoTrf::TranslateZ3D(Spb1ztr)*GeoTrf::TranslateY3D(-Spb1ytr)*GeoTrf::TranslateX3D(Spb1xtr)*GeoTrf::RotateY3D(-180*Gaudi::Units::deg)*GeoTrf::RotateX3D(-(Alfa/2)));
198 sectorenvelopes2l->add(xtl);
199 GeoLogVol *lvspbl = new GeoLogVol("LAr::DM::SpliceBoxl",&SpliceBox,alu);
200 GeoPhysVol *spliceboxl = new GeoPhysVol(lvspbl);
201 sectorenvelopes2l->add(spliceboxl);
202
204 GeoTrd *Trd1 = new GeoTrd(Spb0xhlen1, Spb0xhlen2, Spb0yhlen1, Spb0yhlen2, Spb0zhlen);
205 GeoTrap *GeoTrap2 = new GeoTrap(Spb2zhlen, Spb2theta, Spb2phi, Spb2yzn, Spb2xynzn, Spb2xypzn, Spb2angn, Spb2yzp, Spb2xynzp, Spb2xypzp, Spb2angp);
206 GeoTrap *GeoTrap3 = new GeoTrap(Spb3zhlen, Spb3theta, Spb3phi, Spb3yzn, Spb3xynzn, Spb3xypzn, Spb3angn, Spb3yzp, Spb3xynzp, Spb3xypzp, Spb3angp);
207
208 GeoTransform *xt1 = new GeoTransform (GeoTrf::TranslateY3D(-Spb0ytr)*GeoTrf::RotateX3D(Spb0xrot*Gaudi::Units::deg));
209 spliceboxr->add(xt1);
210 spliceboxl->add(xt1);
211 GeoLogVol *lt1 = new GeoLogVol("LAr::DM::TBox1",Trd1,air);
212 GeoPhysVol *tbox1 = new GeoPhysVol(lt1);
213 spliceboxr->add(tbox1);
214 spliceboxl->add(tbox1);
215
216 GeoTransform *xt2 = new GeoTransform (GeoTrf::TranslateZ3D(Spb2ztr)*GeoTrf::TranslateY3D(Spb2ytr));
217 spliceboxr->add(xt2);
218 spliceboxl->add(xt2);
219 GeoLogVol *lt2 = new GeoLogVol("LAr::DM::TBox2",GeoTrap2,air);
220 GeoPhysVol *tbox2 = new GeoPhysVol(lt2);
221 spliceboxr->add(tbox2);
222 spliceboxl->add(tbox2);
223
224 GeoTransform *xt3 = new GeoTransform (GeoTrf::TranslateZ3D(-Spb3ztr));
225 spliceboxr->add(xt3);
226 spliceboxl->add(xt3);
227 GeoLogVol *lt3 = new GeoLogVol("LAr::DM::TBox3",GeoTrap3,air);
228 GeoPhysVol *tbox3 = new GeoPhysVol(lt3);
229 spliceboxr->add(tbox3);
230 spliceboxl->add(tbox3);
231
232 //-------------- Place volumes in LAr Envelope -------------------
233
234 TRANSFUNCTION seA2r = Pow(GeoTrf::RotateZ3D(1.0),8*f-(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
235 TRANSFUNCTION seA2l = Pow(GeoTrf::RotateZ3D(1.0),8*f+(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
236 TRANSFUNCTION seC2r = Pow(GeoTrf::RotateZ3D(1.0),8*f-(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
237 TRANSFUNCTION seC2l = Pow(GeoTrf::RotateZ3D(1.0),8*f+(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
238 TRANSFUNCTION seA2Vup = Pow(GeoTrf::RotateZ3D(1.0),f+(9*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
239 TRANSFUNCTION seA2Vd = Pow(GeoTrf::RotateZ3D(1.0),f-(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
240 TRANSFUNCTION seA2H = Pow(GeoTrf::RotateZ3D(1.0),8*f+(Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
241 TRANSFUNCTION seC2Vup = Pow(GeoTrf::RotateZ3D(1.0),f+(9*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
242 TRANSFUNCTION seC2Vd = Pow(GeoTrf::RotateZ3D(1.0),f-(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
243 TRANSFUNCTION seC2H = Pow(GeoTrf::RotateZ3D(1.0),8*f+(Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
244
245 GeoSerialTransformer *setA2r = new GeoSerialTransformer(sectorenvelopes2r,&seA2r, 2);
246 GeoSerialTransformer *setA2l = new GeoSerialTransformer(sectorenvelopes2l,&seA2l, 2);
247 GeoSerialTransformer *setC2r = new GeoSerialTransformer(sectorenvelopes2r,&seC2r, 2);
248 GeoSerialTransformer *setC2l = new GeoSerialTransformer(sectorenvelopes2l,&seC2l, 2);
249 GeoSerialTransformer *setA2Vup = new GeoSerialTransformer(sectorenvelopes2vup,&seA2Vup, 1);
250 GeoSerialTransformer *setA2Vd = new GeoSerialTransformer(sectorenvelopes2vd,&seA2Vd, 1);
251 GeoSerialTransformer *setA2H = new GeoSerialTransformer(sectorenvelopes2h,&seA2H, 2);
252 GeoSerialTransformer *setC2Vup = new GeoSerialTransformer(sectorenvelopes2vup,&seC2Vup, 1);
253 GeoSerialTransformer *setC2Vd = new GeoSerialTransformer(sectorenvelopes2vd,&seC2Vd, 1);
254 GeoSerialTransformer *setC2H = new GeoSerialTransformer(sectorenvelopes2h,&seC2H, 2);
255
256 envelope->add(setA2r);
257 envelope->add(setA2l);
258 envelope->add(setC2r);
259 envelope->add(setC2l);
260 envelope->add(setA2Vup);
261 envelope->add(setA2Vd);
262 envelope->add(setA2H);
263 envelope->add(setC2Vup);
264 envelope->add(setC2Vd);
265 envelope->add(setC2H);
266}
Endc
static const double Endc
Definition CrackDMConstruction.cxx:55
DYc
static const double DYc
Definition CrackDMConstruction.cxx:58
add
bool add(const std::string &hname, TKey *tobj)
Definition fastadd.cxx:55

Variable Documentation

◆ Alfa

const double Alfa =360*Gaudi::Units::deg/NCrates
static

Definition at line 52 of file CrackDMConstruction.cxx.

◆ DYa

const double DYa =1155*tan(Alfa/2)
static

Definition at line 56 of file CrackDMConstruction.cxx.

◆ DYb

const double DYb =1695.2*tan(Alfa/2)
static

Definition at line 57 of file CrackDMConstruction.cxx.

◆ DYc

const double DYc =2771.6*tan(Alfa/2)
static

Definition at line 58 of file CrackDMConstruction.cxx.

◆ Enda

const double Enda =1155
static

Definition at line 53 of file CrackDMConstruction.cxx.

◆ Endb

const double Endb =1695.2
static

Definition at line 54 of file CrackDMConstruction.cxx.

◆ Endc

const double Endc =2771.6
static

Definition at line 55 of file CrackDMConstruction.cxx.

◆ NCrates

const unsigned int NCrates =16
static

Definition at line 51 of file CrackDMConstruction.cxx.

Generated on for ATLAS Offline Software by doxygen 1.14.0