+ // Upper part of the readout chambers (readout plane + fee)
+ //
+ // UDxx Wacosit frames of amp. region (C)
+ // UExx Inner volumes of the frame (Air)
+ // UFxx Aluminum frame of back panel (Al)
+ // UGxx Inner volumes of the back panel (Air)
+ //
+ // Inner material layers
+ //
+ // UHxx Radiator (Rohacell)
+ // UJxx Drift volume (Xe/CO2)
+ // UKxx Amplification volume (Xe/CO2)
+ // UWxx Wire plane (Cu)
+ // ULxx Pad plane (Cu)
+ // UYxx Glue layer (Epoxy)
+ // UMxx Support structure (Rohacell)
+ // UNxx ROB base material (C)
+ // UOxx ROB copper (Cu)
+ // UVxx ROB other materials (Cu)
+ //
+
+ const Int_t kNparTrd = 4;
+ const Int_t kNparCha = 3;
+
+ Float_t xpos;
+ Float_t ypos;
+ Float_t zpos;
+
+ Float_t parTrd[kNparTrd];
+ Float_t parCha[kNparCha];
+
+ Char_t cTagV[6];
+ Char_t cTagM[5];
+
+ // The TRD mother volume for one sector (Air), full length in z-direction
+ // Provides material for side plates of super module
+ parTrd[0] = fgkSwidth1/2.0;
+ parTrd[1] = fgkSwidth2/2.0;
+ parTrd[2] = fgkSlength/2.0;
+ parTrd[3] = fgkSheight/2.0;
+ gMC->Gsvolu("UTR1","TRD1",idtmed[1302-1],parTrd,kNparTrd);
+
+ // The outer aluminum plates of the super module (Al)
+ parTrd[0] = fgkSwidth1/2.0;
+ parTrd[1] = fgkSwidth2/2.0;
+ parTrd[2] = fgkSlength/2.0;
+ parTrd[3] = fgkSheight/2.0;
+ gMC->Gsvolu("UTS1","TRD1",idtmed[1301-1],parTrd,kNparTrd);
+
+ // The inner part of the TRD mother volume for one sector (Air),
+ // full length in z-direction
+ parTrd[0] = fgkSwidth1/2.0 - fgkSMpltT;
+ parTrd[1] = fgkSwidth2/2.0 - fgkSMpltT;
+ parTrd[2] = fgkSlength/2.0;
+ parTrd[3] = fgkSheight/2.0 - fgkSMpltT;
+ gMC->Gsvolu("UTI1","TRD1",idtmed[1302-1],parTrd,kNparTrd);
+
+ for (Int_t icham = 0; icham < kNcham; icham++) {
+ for (Int_t iplan = 0; iplan < kNplan; iplan++) {
+
+ Int_t iDet = GetDetectorSec(iplan,icham);
+
+ // The lower part of the readout chambers (drift volume + radiator)
+ // The aluminum frames
+ sprintf(cTagV,"UA%02d",iDet);
+ parCha[0] = fCwidth[iplan]/2.0;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0;
+ parCha[2] = fgkCraH/2.0 + fgkCdrH/2.0;
+ fChamberUAboxd[iDet][0] = parCha[0];
+ fChamberUAboxd[iDet][1] = parCha[1];
+ fChamberUAboxd[iDet][2] = parCha[2];
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
+ // The additional aluminum on the frames
+ // This part has not the correct postion but is just supposed to
+ // represent the missing material. The correct from of the L-shaped
+ // profile would not fit into the alignable volume.
+ sprintf(cTagV,"UZ%02d",iDet);
+ parCha[0] = fgkCroW/2.0;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0;
+ parCha[2] = fgkCalW/2.0;
+ fChamberUAboxd[iDet][0] = fChamberUAboxd[iDet][0] + fgkCroW;
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
+ // The Wacosit frames
+ sprintf(cTagV,"UB%02d",iDet);
+ parCha[0] = fCwidth[iplan]/2.0 - fgkCalT;
+ parCha[1] = -1.0;
+ parCha[2] = -1.0;
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
+ // The glue around the radiator
+ sprintf(cTagV,"UX%02d",iDet);
+ parCha[0] = fCwidth[iplan]/2.0 - fgkCalT - fgkCclsT;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0 - fgkCclfT;
+ parCha[2] = fgkCraH/2.0;
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1311-1],parCha,kNparCha);
+ // The inner part of radiator (air)
+ sprintf(cTagV,"UC%02d",iDet);
+ parCha[0] = fCwidth[iplan]/2.0 - fgkCalT - fgkCclsT - fgkCglT;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0 - fgkCclfT - fgkCglT;
+ parCha[2] = -1.0;
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
+
+ // The upper part of the readout chambers (amplification volume)
+ // The Wacosit frames
+ sprintf(cTagV,"UD%02d",iDet);
+ parCha[0] = fCwidth[iplan]/2.0 + fgkCroW;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0;
+ parCha[2] = fgkCamH/2.0;
+ fChamberUDboxd[iDet][0] = parCha[0];
+ fChamberUDboxd[iDet][1] = parCha[1];
+ fChamberUDboxd[iDet][2] = parCha[2];
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
+ // The inner part of the Wacosit frame (air)
+ sprintf(cTagV,"UE%02d",iDet);
+ parCha[0] = fCwidth[iplan]/2.0 + fgkCroW - fgkCcuT;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0 - fgkCcuT;
+ parCha[2] = -1.;
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
+
+ // The support structure (pad plane, back panel, readout boards)
+ // The aluminum frames
+ sprintf(cTagV,"UF%02d",iDet);
+ parCha[0] = fCwidth[iplan]/2.0 + fgkCroW;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0;
+ parCha[2] = fgkCroH/2.0;
+ fChamberUFboxd[iDet][0] = parCha[0];
+ fChamberUFboxd[iDet][1] = parCha[1];
+ fChamberUFboxd[iDet][2] = parCha[2];
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
+ // The inner part of the aluminum frames
+ sprintf(cTagV,"UG%02d",iDet);
+ parCha[0] = fCwidth[iplan]/2.0 + fgkCroW - fgkCauT;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0 - fgkCauT;
+ parCha[2] = -1.0;
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
+
+ // The material layers inside the chambers
+ // Rohacell layer (radiator)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkRaThick/2.0;
+ sprintf(cTagV,"UH%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
+ // Xe/Isobutane layer (drift volume)
+ parCha[0] = fCwidth[iplan]/2.0 - fgkCalT - fgkCclsT;
+ parCha[1] = fClength[iplan][icham]/2.0 - fgkHspace/2.0 - fgkCclfT;
+ parCha[2] = fgkDrThick/2.0;
+ sprintf(cTagV,"UJ%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
+ // Xe/Isobutane layer (amplification volume)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkAmThick/2.0;
+ sprintf(cTagV,"UK%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
+ // Cu layer (wire plane)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkWrThick/2.0;
+ sprintf(cTagV,"UW%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1303-1],parCha,kNparCha);
+ // Cu layer (pad plane)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkCuThick/2.0;
+ sprintf(cTagV,"UL%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
+ // Epoxy layer (glue)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkGlThick/2.0;
+ sprintf(cTagV,"UY%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1311-1],parCha,kNparCha);
+ // G10 layer (support structure / honeycomb)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkSuThick/2.0;
+ sprintf(cTagV,"UM%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1310-1],parCha,kNparCha);
+ // G10 layer (PCB readout board)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkRpThick/2;
+ sprintf(cTagV,"UN%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
+ // Cu layer (traces in readout board)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkRcThick/2.0;
+ sprintf(cTagV,"UO%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1306-1],parCha,kNparCha);
+ // Cu layer (other material on in readout board)
+ parCha[0] = -1.0;
+ parCha[1] = -1.0;
+ parCha[2] = fgkRoThick/2.0;
+ sprintf(cTagV,"UV%02d",iDet);
+ gMC->Gsvolu(cTagV,"BOX ",idtmed[1304-1],parCha,kNparCha);
+
+ // Position the layers in the chambers
+ xpos = 0.0;
+ ypos = 0.0;
+ // Lower part
+ // Rohacell layer (radiator)
+ zpos = fgkRaZpos;
+ sprintf(cTagV,"UH%02d",iDet);
+ sprintf(cTagM,"UC%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // Xe/Isobutane layer (drift volume)
+ zpos = fgkDrZpos;
+ sprintf(cTagV,"UJ%02d",iDet);
+ sprintf(cTagM,"UB%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // Upper part
+ // Xe/Isobutane layer (amplification volume)
+ zpos = fgkAmZpos;
+ sprintf(cTagV,"UK%02d",iDet);
+ sprintf(cTagM,"UE%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // Cu layer (wire plane inside amplification volume)
+ zpos = fgkWrZpos;
+ sprintf(cTagV,"UW%02d",iDet);
+ sprintf(cTagM,"UK%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // Readout part + support plane
+ // Cu layer (pad plane)
+ zpos = fgkCuZpos;
+ sprintf(cTagV,"UL%02d",iDet);
+ sprintf(cTagM,"UG%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // Epoxy layer (glue)
+ zpos = fgkGlZpos;
+ sprintf(cTagV,"UY%02d",iDet);
+ sprintf(cTagM,"UG%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // G10 layer (support structure)
+ zpos = fgkSuZpos;
+ sprintf(cTagV,"UM%02d",iDet);
+ sprintf(cTagM,"UG%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // G10 layer (PCB readout board)
+ zpos = fgkRpZpos;
+ sprintf(cTagV,"UN%02d",iDet);
+ sprintf(cTagM,"UG%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // Cu layer (traces in readout board)
+ zpos = fgkRcZpos;
+ sprintf(cTagV,"UO%02d",iDet);
+ sprintf(cTagM,"UG%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // Cu layer (other materials on readout board)
+ zpos = fgkRoZpos;
+ sprintf(cTagV,"UV%02d",iDet);
+ sprintf(cTagM,"UG%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+
+ // Position the inner volumes of the chambers in the frames
+ xpos = 0.0;
+ ypos = 0.0;
+ // The inner part of the radiator
+ zpos = 0.0;
+ sprintf(cTagV,"UC%02d",iDet);
+ sprintf(cTagM,"UX%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // The glue around the radiator
+ zpos = fgkCraH/2.0 - fgkCdrH/2.0 - fgkCraH/2.0;
+ sprintf(cTagV,"UX%02d",iDet);
+ sprintf(cTagM,"UB%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // The lower Wacosit frame inside the aluminum frame
+ zpos = 0.0;
+ sprintf(cTagV,"UB%02d",iDet);
+ sprintf(cTagM,"UA%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // The inside of the upper Wacosit frame
+ zpos = 0.0;
+ sprintf(cTagV,"UE%02d",iDet);
+ sprintf(cTagM,"UD%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+ // The inside of the upper aluminum frame
+ zpos = 0.0;
+ sprintf(cTagV,"UG%02d",iDet);
+ sprintf(cTagM,"UF%02d",iDet);
+ gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
+
+ // Position the frames of the chambers in the TRD mother volume
+ xpos = 0.0;
+ ypos = fClength[iplan][0] + fClength[iplan][1] + fClength[iplan][2]/2.0;
+ for (Int_t ic = 0; ic < icham; ic++) {
+ ypos -= fClength[iplan][ic];
+ }
+ ypos -= fClength[iplan][icham]/2.0;
+ zpos = fgkVrocsm + fgkSMpltT + fgkCraH/2.0 + fgkCdrH/2.0 - fgkSheight/2.0
+ + iplan * (fgkCH + fgkVspace);
+ // The lower aluminum frame, radiator + drift region
+ sprintf(cTagV,"UA%02d",iDet);
+ fChamberUAorig[iDet][0] = xpos;
+ fChamberUAorig[iDet][1] = ypos;
+ fChamberUAorig[iDet][2] = zpos;
+ // The upper G10 frame, amplification region
+ sprintf(cTagV,"UD%02d",iDet);
+ zpos += fgkCamH/2.0 + fgkCraH/2.0 + fgkCdrH/2.0;
+ fChamberUDorig[iDet][0] = xpos;
+ fChamberUDorig[iDet][1] = ypos;
+ fChamberUDorig[iDet][2] = zpos;
+ // The upper aluminum frame
+ sprintf(cTagV,"UF%02d",iDet);
+ zpos += fgkCroH/2.0 + fgkCamH/2.0;
+ fChamberUForig[iDet][0] = xpos;
+ fChamberUForig[iDet][1] = ypos;
+ fChamberUForig[iDet][2] = zpos;
+
+ }
+ }
+
+ // Create the volumes of the super module frame
+ CreateFrame(idtmed);
+
+ // Create the volumes of the services
+ CreateServices(idtmed);
+
+ for (Int_t icham = 0; icham < kNcham; icham++) {
+ for (Int_t iplan = 0; iplan < kNplan; iplan++) {
+ GroupChamber(iplan,icham,idtmed);
+ }
+ }
+
+ xpos = 0.0;
+ ypos = 0.0;
+ zpos = 0.0;
+ gMC->Gspos("UTI1",1,"UTS1",xpos,ypos,zpos,0,"ONLY");
+
+ xpos = 0.0;
+ ypos = 0.0;
+ zpos = 0.0;
+ gMC->Gspos("UTS1",1,"UTR1",xpos,ypos,zpos,0,"ONLY");
+
+ // Put the TRD volumes into the space frame mother volumes
+ // if enabled via status flag
+ xpos = 0.0;
+ ypos = 0.0;
+ zpos = 0.0;
+ for (Int_t isect = 0; isect < kNsect; isect++) {
+ if (fSMstatus[isect]) {
+ sprintf(cTagV,"BTRD%d",isect);
+ gMC->Gspos("UTR1",1,cTagV,xpos,ypos,zpos,0,"ONLY");
+ }
+ }