[u/mrichter/AliRoot.git] / EMCAL / AliEMCALSpaceFrame.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

#include <TVirtualMC.h>
#include <TGeoManager.h>
#include <TGeoVolume.h>
#include <TGeoMedium.h>
#include <TGeoMatrix.h>
#include <TGeoVolume.h>
#include <TGeoTube.h>
#include <TGeoCone.h>
#include <TGeoPcon.h>
#include <TGeoCompositeShape.h>

#include "AliConst.h"
#include "AliEMCALSpaceFrame.h"
#include "AliMagF.h"
#include "AliRun.h"
#include "AliLog.h"
 
ClassImp(AliEMCALSpaceFrame)
 
//_____________________________________________________________________________
AliEMCALSpaceFrame::AliEMCALSpaceFrame() 
  : TNamed("EMCALSpaceFrame","Steel Space Frame that supports EMCAL"),
    fNumCross(0),
    fNumSubSets(0),
    fTotalHalfWidth(0.),
    fBeginPhi(0.),
    fEndPhi(0.),
    fTotalPhi(0.),
    fBeginRadius(0.),
    fHalfFrameTrans(0.),
    fFlangeHeight(0.),
    fFlangeWidth(0.),
    fRibHeight(0.),
    fRibWidth(0.),
    fCrossBottomWidth(0.),
    fCrossTopWidth(0.),
    fCrossBottomHeight(0.),
    fCrossBottomRadThick(0.),
    fCrossBeamArcLength(0.),
    fCrossBottomStartRadius(0.),
    fCrossTopHeight(0.),
    fCrossTopRadThick(0.),
    fCrossTopStart(0.),
    fEndRadius(0.),
    fEndBeamRadThick(0),
    fEndBeamBeginRadius(0)
{
  // default constructor for EMCAL Space Frame
  //initialize parameters
  fNumCross = 12; 
  fNumSubSets = 3;
  fTotalHalfWidth = 152.3;   // Half Width of a Half Frame
			     // (CalFrame comes in 2 sections) 
  fBeginPhi = 76.8;
  fEndPhi = 193.03;
  fBeginRadius = 490.;
  
  fHalfFrameTrans = fTotalHalfWidth+57.2/2.;  // Half Frame Connector is 57.2cm wide,
			    // Supermodule is 340cm wide
                            // Sources: HALF-FRAME-CONNECTOR-27E226A.pdf
                            // provided by LBL

  fFlangeWidth = 15.2;
  fRibWidth = 1.5;
  fCrossBottomHeight = 15.2;
  fCrossBottomRadThick = 1.5;
  fCrossTopHeight = 1.5;
  fCrossTopRadThick = 35./2.;

  fTotalPhi = fEndPhi - fBeginPhi;
  fFlangeHeight = fBeginRadius + 3.;
  fRibHeight = fFlangeHeight + 35;
  fCrossBottomWidth = 0.5/(Double_t)fNumSubSets * (2.*fTotalHalfWidth - 8. * fFlangeWidth);
  fCrossTopWidth = fCrossBottomWidth; // fCrossBottomWidth + fFlangeWidth - fRibWidth;
                                      // for future release pending
                                      // overlap correction - new TGeoVolume creation

  fCrossBeamArcLength = (112.62597)/(fNumCross-1)-.001; // To account for shape of TGeoBBox
  fCrossBottomStartRadius = fBeginRadius + fCrossBottomRadThick;
  fCrossTopStart = fBeginRadius + 2.*fCrossBottomRadThick + fCrossTopRadThick+0.015; // 0.015 is a 
										     // bubblegum and duct tape
										     // fix for an overlap problem
										     // will be worked out in future releases 
  fEndRadius = fRibHeight+1.15;
  fEndBeamRadThick = fCrossBottomRadThick+fCrossTopRadThick;
  fEndBeamBeginRadius = fBeginRadius + fEndBeamRadThick;
}

//_____________________________________________________________________________
AliEMCALSpaceFrame::AliEMCALSpaceFrame(const AliEMCALSpaceFrame &frame) 
  : TNamed(frame.GetName(),frame.GetTitle()),
    fNumCross(frame.fNumCross),
    fNumSubSets(frame.fNumSubSets),
    fTotalHalfWidth(frame.fTotalHalfWidth),
    fBeginPhi(frame.fBeginPhi),
    fEndPhi(frame.fEndPhi),
    fTotalPhi(frame.fTotalPhi),
    fBeginRadius(frame.fBeginRadius),
    fHalfFrameTrans(frame.fHalfFrameTrans),
    fFlangeHeight(frame.fFlangeHeight),
    fFlangeWidth(frame.fFlangeWidth),
    fRibHeight(frame.fRibHeight),
    fRibWidth(frame.fRibWidth),
    fCrossBottomWidth(frame.fCrossBottomWidth),
    fCrossTopWidth(frame.fCrossTopWidth),
    fCrossBottomHeight(frame.fCrossBottomHeight),
    fCrossBottomRadThick(frame.fCrossBottomRadThick),
    fCrossBeamArcLength(frame.fCrossBeamArcLength),
    fCrossBottomStartRadius(frame.fCrossBottomStartRadius),
    fCrossTopHeight(frame.fCrossTopHeight),
    fCrossTopRadThick(frame.fCrossTopRadThick),
    fCrossTopStart(frame.fCrossTopStart),
    fEndRadius(frame.fEndRadius),
    fEndBeamRadThick(frame.fEndBeamRadThick),
    fEndBeamBeginRadius(frame.fEndBeamBeginRadius)
{
  // copy constructor for EMCAL Space Frame

}

//_____________________________________________________________________________
void AliEMCALSpaceFrame::CreateGeometry()
{
  // create geometry
  
  AliDebug(1,"Create CalFrame Geometry");
  //////////////////////////////////////Setup/////////////////////////////////////////
  TGeoVolume* top = gGeoManager->GetVolume("ALIC");
  TGeoMedium *steel = gGeoManager->GetMedium("EMCAL_S steel$");
  TGeoMedium *air = gGeoManager->GetMedium("EMCAL_Air$");

	
  //////////////////////////////////// Volumes ///////////////////////////////////////  
  TGeoVolume *calFrameMO = 
    gGeoManager->MakeTubs("CalFrame", air, fBeginRadius-2.1,fEndRadius,
			  fTotalHalfWidth*3,fBeginPhi-3,fEndPhi+3);	// Mother Volume

  calFrameMO->SetVisibility(kFALSE);

  // Half Frame Mother Volume
  TGeoVolume *calHalfFrameMO = 
    gGeoManager->MakeTubs("HalfFrame", air, fBeginRadius-2,fEndRadius,
			  fTotalHalfWidth,fBeginPhi-2.9,fEndPhi+2.9);

  calHalfFrameMO->SetVisibility(kFALSE);

  TGeoVolume *endBeams = 
    gGeoManager->MakeBox("End Beams", steel, fEndBeamRadThick, fCrossTopHeight, fTotalHalfWidth); // End Beams

  TGeoVolume *skin = 
    gGeoManager->MakeTubs("skin", steel, fRibHeight+0.15, fEndRadius, 
  			  fTotalHalfWidth, fBeginPhi, fEndPhi);// back frame

  TGeoVolume *flangeVolume = 
    gGeoManager->MakeTubs("supportBottom", steel, fBeginRadius, fFlangeHeight, 
			  fFlangeWidth, fBeginPhi, fEndPhi);   		// FlangeVolume Beams

  TGeoVolume *ribVolume = 
    gGeoManager->MakeTubs("RibVolume", steel, fFlangeHeight, fRibHeight, fRibWidth, fBeginPhi, fEndPhi);

  TGeoVolume *subSetCross = 
    gGeoManager->MakeTubs("subSetCross", air, fBeginRadius-1,  fBeginRadius+2*fCrossBottomRadThick+
                          2*fCrossTopRadThick+0.15, fCrossBottomWidth, fBeginPhi, fEndPhi); 		// Cross Beam Containers
  subSetCross->SetVisibility(kFALSE);
  /*						// Obsolete for now
  TGeoVolume *subSetCrossTop = 
    gGeoManager->MakeTubs("SubSetCrossTop", air, fBeginRadius+2*fCrossBottomRadThick-1, fBeginRadius+2*fCrossBottomRadThick+
                          2*fCrossTopRadThick+1, fCrossTopWidth, fBeginPhi, fEndPhi);	// Cross 
  subSetCrossTop->SetVisibility(kFALSE);
  */                    
  TGeoVolume *crossBottomBeams = 
    gGeoManager->MakeBox("crossBottom", steel, fCrossBottomRadThick, fCrossBottomHeight, fCrossBottomWidth); // Cross Beams

  TGeoVolume *crossTopBeams = 
    gGeoManager->MakeBox("crossTop", steel, fCrossTopRadThick, fCrossTopHeight, fCrossTopWidth); // Cross Beams
  
  TGeoTranslation *trTEST = new TGeoTranslation();
  TGeoRotation *rotTEST = new TGeoRotation();
  
  Double_t conv = TMath::Pi()/180.;
  Double_t radAngle = 0;
  Double_t endBeamParam=.4;
  //cout<<"\nfCrossBottomStartRadius: "<<fCrossBottomStartRadius<<"\n";
  
  for(Int_t i = 0; i < fNumCross; i++){
    
    Double_t loopPhi = fBeginPhi + 1.8;

    // Cross Bottom Beams
    
    radAngle = (loopPhi + i*fCrossBeamArcLength)*conv; 
    
    rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); //  SetTranslation(Double_t dx, Double_t dy, Double_t dz);
    trTEST->SetTranslation(cos(radAngle)*fCrossBottomStartRadius, sin(radAngle)*fCrossBottomStartRadius,0);

    TGeoCombiTrans *combo = new TGeoCombiTrans(*trTEST, *rotTEST); 	// TGeoTranslation &tr, const TGeoRotation &rot);
    combo->RegisterYourself();
    crossBottomBeams->SetVisibility(1);
    subSetCross->AddNode(crossBottomBeams, i+1, combo);
    if (i != 0 && i!=fNumCross-1){
    // Cross Bottom Beams
    rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); //  SetTranslation(Double_t dx, Double_t dy, Double_t dz);
    trTEST->SetTranslation(cos(radAngle)*fCrossTopStart, sin(radAngle)*fCrossTopStart,0);
    crossTopBeams->SetVisibility(1);
    subSetCross->AddNode(crossTopBeams, i+1,  new TGeoCombiTrans(*trTEST, *rotTEST));
    }
    

    else if(i ==0){
		rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); //  SetTranslation(Double_t dx, Double_t dy, Double_t dz);
    trTEST->SetTranslation(cos((77-endBeamParam)*conv)*(fEndBeamBeginRadius), sin((77-endBeamParam)*conv)*(fEndBeamBeginRadius),0);
    endBeams->SetVisibility(1);
    calHalfFrameMO->AddNode(endBeams, 1,  new TGeoCombiTrans(*trTEST, *rotTEST));
    }
    else{
    rotTEST->SetAngles(193.03, 0, 0); //  SetTranslation(Double_t dx, Double_t dy, Double_t dz);
    trTEST->SetTranslation(cos((193.03+endBeamParam)*conv)*(fEndBeamBeginRadius)/*more duct tape*/, sin((193.03+endBeamParam)*conv)*(fEndBeamBeginRadius),0);
    endBeams->SetVisibility(1);
    calHalfFrameMO->AddNode(endBeams, 2,  new TGeoCombiTrans(*trTEST, *rotTEST));
    }  
  }
  
  //Beam Containers

  // Translations 

  TGeoTranslation *origin1 = new TGeoTranslation(0,0,0); // Equivalent to gGeoIdentity
  TGeoTranslation *origin2 = new TGeoTranslation(0,0,2*(fCrossBottomWidth+fFlangeWidth));
  TGeoTranslation *origin3 = new TGeoTranslation(0,0,-2*(fCrossBottomWidth+fFlangeWidth));

  // FlangeVolume translations  
  TGeoTranslation *str1 = new TGeoTranslation(0,0,-3*(fCrossBottomWidth+fFlangeWidth));
  TGeoTranslation *str2 = new TGeoTranslation(0,0,-(fCrossBottomWidth+fFlangeWidth));
  TGeoTranslation *str3 = new TGeoTranslation(0,0,(fCrossBottomWidth+fFlangeWidth));
  TGeoTranslation *str4 = new TGeoTranslation(0,0,3*(fCrossBottomWidth+fFlangeWidth));

  // Half Frame Translations
  TGeoTranslation *halfTrans1 =  new TGeoTranslation(0,0,fHalfFrameTrans);
  TGeoTranslation *halfTrans2 =  new TGeoTranslation(0,0,-fHalfFrameTrans);

  // Beams Volume 
  calHalfFrameMO->AddNode(flangeVolume, 1, str1);
  calHalfFrameMO->AddNode(flangeVolume, 2, str2);
  calHalfFrameMO->AddNode(flangeVolume, 3, str3);
  calHalfFrameMO->AddNode(flangeVolume, 4, str4);
  
  calHalfFrameMO->AddNode(ribVolume, 1, str1);
  calHalfFrameMO->AddNode(ribVolume, 2, str2);
  calHalfFrameMO->AddNode(ribVolume, 3, str3);
  calHalfFrameMO->AddNode(ribVolume, 4, str4);
  
  // Cross Beams  
  calHalfFrameMO->AddNode(subSetCross, 1, origin1);
  calHalfFrameMO->AddNode(subSetCross, 2, origin2);
  calHalfFrameMO->AddNode(subSetCross, 3, origin3);
/*					// Obsolete for now
  calHalfFrameMO->AddNode(subSetCrossTop, 1, origin1);
  calHalfFrameMO->AddNode(subSetCrossTop, 2, origin2);
  calHalfFrameMO->AddNode(subSetCrossTop, 3, origin3);
*/

  calHalfFrameMO->AddNode(skin, 1, gGeoIdentity);
  
  calFrameMO->AddNode(calHalfFrameMO, 1, halfTrans1);
  calFrameMO->AddNode(calHalfFrameMO, 2, halfTrans2);
 
  top->AddNode(calFrameMO,1,gGeoIdentity);
//  cout<<"**********************************\nfEndRadius:\t"<<fEndRadius;
}
Commit	Line	Data
e3697737	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	#include <TVirtualMC.h>
	17	#include <TGeoManager.h>
	18	#include <TGeoVolume.h>
	19	#include <TGeoMedium.h>
	20	#include <TGeoMatrix.h>
	21	#include <TGeoVolume.h>
	22	#include <TGeoTube.h>
	23	#include <TGeoCone.h>
	24	#include <TGeoPcon.h>
	25	#include <TGeoCompositeShape.h>
	26
	27	#include "AliConst.h"
	28	#include "AliEMCALSpaceFrame.h"
	29	#include "AliMagF.h"
	30	#include "AliRun.h"
c93255fe	31	#include "AliLog.h"
e3697737	32
	33	ClassImp(AliEMCALSpaceFrame)
	34
	35	//_____________________________________________________________________________
	36	AliEMCALSpaceFrame::AliEMCALSpaceFrame()
	37	: TNamed("EMCALSpaceFrame","Steel Space Frame that supports EMCAL"),
	38	fNumCross(0),
	39	fNumSubSets(0),
	40	fTotalHalfWidth(0.),
	41	fBeginPhi(0.),
	42	fEndPhi(0.),
	43	fTotalPhi(0.),
	44	fBeginRadius(0.),
	45	fHalfFrameTrans(0.),
	46	fFlangeHeight(0.),
	47	fFlangeWidth(0.),
	48	fRibHeight(0.),
	49	fRibWidth(0.),
	50	fCrossBottomWidth(0.),
	51	fCrossTopWidth(0.),
	52	fCrossBottomHeight(0.),
	53	fCrossBottomRadThick(0.),
	54	fCrossBeamArcLength(0.),
	55	fCrossBottomStartRadius(0.),
	56	fCrossTopHeight(0.),
	57	fCrossTopRadThick(0.),
	58	fCrossTopStart(0.),
	59	fEndRadius(0.),
	60	fEndBeamRadThick(0),
	61	fEndBeamBeginRadius(0)
	62	{
	63	// default constructor for EMCAL Space Frame
	64	//initialize parameters
	65	fNumCross = 12;
	66	fNumSubSets = 3;
	67	fTotalHalfWidth = 152.3; // Half Width of a Half Frame
	68	// (CalFrame comes in 2 sections)
	69	fBeginPhi = 76.8;
	70	fEndPhi = 193.03;
	71	fBeginRadius = 490.;
	72
	73	fHalfFrameTrans = fTotalHalfWidth+57.2/2.; // Half Frame Connector is 57.2cm wide,
	74	// Supermodule is 340cm wide
	75	// Sources: HALF-FRAME-CONNECTOR-27E226A.pdf
	76	// provided by LBL
	77
	78	fFlangeWidth = 15.2;
	79	fRibWidth = 1.5;
	80	fCrossBottomHeight = 15.2;
	81	fCrossBottomRadThick = 1.5;
	82	fCrossTopHeight = 1.5;
	83	fCrossTopRadThick = 35./2.;
	84
	85	fTotalPhi = fEndPhi - fBeginPhi;
	86	fFlangeHeight = fBeginRadius + 3.;
	87	fRibHeight = fFlangeHeight + 35;
	88	fCrossBottomWidth = 0.5/(Double_t)fNumSubSets * (2.fTotalHalfWidth - 8. fFlangeWidth);
	89	fCrossTopWidth = fCrossBottomWidth; // fCrossBottomWidth + fFlangeWidth - fRibWidth;
	90	// for future release pending
	91	// overlap correction - new TGeoVolume creation
	92
	93	fCrossBeamArcLength = (112.62597)/(fNumCross-1)-.001; // To account for shape of TGeoBBox
	94	fCrossBottomStartRadius = fBeginRadius + fCrossBottomRadThick;
	95	fCrossTopStart = fBeginRadius + 2.*fCrossBottomRadThick + fCrossTopRadThick+0.015; // 0.015 is a
96	// bubblegum and duct tape
97	// fix for an overlap problem
98	// will be worked out in future releases
99	fEndRadius = fRibHeight+1.15;
100	fEndBeamRadThick = fCrossBottomRadThick+fCrossTopRadThick;
101	fEndBeamBeginRadius = fBeginRadius + fEndBeamRadThick;
102	}
103
104	//_____________________________________________________________________________
105	AliEMCALSpaceFrame::AliEMCALSpaceFrame(const AliEMCALSpaceFrame &frame)
106	: TNamed(frame.GetName(),frame.GetTitle()),
107	fNumCross(frame.fNumCross),
108	fNumSubSets(frame.fNumSubSets),
109	fTotalHalfWidth(frame.fTotalHalfWidth),
110	fBeginPhi(frame.fBeginPhi),
111	fEndPhi(frame.fEndPhi),
112	fTotalPhi(frame.fTotalPhi),
113	fBeginRadius(frame.fBeginRadius),
114	fHalfFrameTrans(frame.fHalfFrameTrans),
115	fFlangeHeight(frame.fFlangeHeight),
116	fFlangeWidth(frame.fFlangeWidth),
117	fRibHeight(frame.fRibHeight),
118	fRibWidth(frame.fRibWidth),
119	fCrossBottomWidth(frame.fCrossBottomWidth),
120	fCrossTopWidth(frame.fCrossTopWidth),
121	fCrossBottomHeight(frame.fCrossBottomHeight),
122	fCrossBottomRadThick(frame.fCrossBottomRadThick),
123	fCrossBeamArcLength(frame.fCrossBeamArcLength),
124	fCrossBottomStartRadius(frame.fCrossBottomStartRadius),
125	fCrossTopHeight(frame.fCrossTopHeight),
126	fCrossTopRadThick(frame.fCrossTopRadThick),
127	fCrossTopStart(frame.fCrossTopStart),
128	fEndRadius(frame.fEndRadius),
129	fEndBeamRadThick(frame.fEndBeamRadThick),
130	fEndBeamBeginRadius(frame.fEndBeamBeginRadius)
131	{
132	// copy constructor for EMCAL Space Frame
133
134	}
135
136	//_____________________________________________________________________________
137	void AliEMCALSpaceFrame::CreateGeometry()
138	{
93684a6a	139	// create geometry
93684a6a	140
e3697737	141	AliDebug(1,"Create CalFrame Geometry");
	142	//////////////////////////////////////Setup/////////////////////////////////////////
	143	TGeoVolume* top = gGeoManager->GetVolume("ALIC");
	144	TGeoMedium *steel = gGeoManager->GetMedium("EMCAL_S steel$");
	145	TGeoMedium *air = gGeoManager->GetMedium("EMCAL_Air$");
	146
	147
	148	//////////////////////////////////// Volumes ///////////////////////////////////////
	149	TGeoVolume *calFrameMO =
	150	gGeoManager->MakeTubs("CalFrame", air, fBeginRadius-2.1,fEndRadius,
	151	fTotalHalfWidth*3,fBeginPhi-3,fEndPhi+3); // Mother Volume
	152
	153	calFrameMO->SetVisibility(kFALSE);
	154
	155	// Half Frame Mother Volume
	156	TGeoVolume *calHalfFrameMO =
	157	gGeoManager->MakeTubs("HalfFrame", air, fBeginRadius-2,fEndRadius,
	158	fTotalHalfWidth,fBeginPhi-2.9,fEndPhi+2.9);
	159
	160	calHalfFrameMO->SetVisibility(kFALSE);
	161
	162	TGeoVolume *endBeams =
	163	gGeoManager->MakeBox("End Beams", steel, fEndBeamRadThick, fCrossTopHeight, fTotalHalfWidth); // End Beams
	164
	165	TGeoVolume *skin =
	166	gGeoManager->MakeTubs("skin", steel, fRibHeight+0.15, fEndRadius,
	167	fTotalHalfWidth, fBeginPhi, fEndPhi);// back frame
	168
	169	TGeoVolume *flangeVolume =
	170	gGeoManager->MakeTubs("supportBottom", steel, fBeginRadius, fFlangeHeight,
	171	fFlangeWidth, fBeginPhi, fEndPhi); // FlangeVolume Beams
	172
	173	TGeoVolume *ribVolume =
	174	gGeoManager->MakeTubs("RibVolume", steel, fFlangeHeight, fRibHeight, fRibWidth, fBeginPhi, fEndPhi);
	175
	176	TGeoVolume *subSetCross =
	177	gGeoManager->MakeTubs("subSetCross", air, fBeginRadius-1, fBeginRadius+2*fCrossBottomRadThick+
	178	2*fCrossTopRadThick+0.15, fCrossBottomWidth, fBeginPhi, fEndPhi); // Cross Beam Containers
	179	subSetCross->SetVisibility(kFALSE);
	180	/* // Obsolete for now
	181	TGeoVolume *subSetCrossTop =
	182	gGeoManager->MakeTubs("SubSetCrossTop", air, fBeginRadius+2fCrossBottomRadThick-1, fBeginRadius+2fCrossBottomRadThick+
	183	2*fCrossTopRadThick+1, fCrossTopWidth, fBeginPhi, fEndPhi); // Cross
	184	subSetCrossTop->SetVisibility(kFALSE);
	185	*/
	186	TGeoVolume *crossBottomBeams =
	187	gGeoManager->MakeBox("crossBottom", steel, fCrossBottomRadThick, fCrossBottomHeight, fCrossBottomWidth); // Cross Beams
	188
	189	TGeoVolume *crossTopBeams =
	190	gGeoManager->MakeBox("crossTop", steel, fCrossTopRadThick, fCrossTopHeight, fCrossTopWidth); // Cross Beams
	191
	192	TGeoTranslation *trTEST = new TGeoTranslation();
	193	TGeoRotation *rotTEST = new TGeoRotation();
	194
	195	Double_t conv = TMath::Pi()/180.;
	196	Double_t radAngle = 0;
	197	Double_t endBeamParam=.4;
	198	//cout<<"\nfCrossBottomStartRadius: "<<fCrossBottomStartRadius<<"\n";
	199
	200	for(Int_t i = 0; i < fNumCross; i++){
	201
	202	Double_t loopPhi = fBeginPhi + 1.8;
	203
	204	// Cross Bottom Beams
205
206	radAngle = (loopPhi + ifCrossBeamArcLength)conv;
207
208	rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); // SetTranslation(Double_t dx, Double_t dy, Double_t dz);
209	trTEST->SetTranslation(cos(radAngle)fCrossBottomStartRadius, sin(radAngle)fCrossBottomStartRadius,0);
210
211	TGeoCombiTrans combo = new TGeoCombiTrans(trTEST, *rotTEST); // TGeoTranslation &tr, const TGeoRotation &rot);
212	combo->RegisterYourself();
213	crossBottomBeams->SetVisibility(1);
214	subSetCross->AddNode(crossBottomBeams, i+1, combo);
215	if (i != 0 && i!=fNumCross-1){
216	// Cross Bottom Beams
217	rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); // SetTranslation(Double_t dx, Double_t dy, Double_t dz);
218	trTEST->SetTranslation(cos(radAngle)fCrossTopStart, sin(radAngle)fCrossTopStart,0);
219	crossTopBeams->SetVisibility(1);
220	subSetCross->AddNode(crossTopBeams, i+1, new TGeoCombiTrans(trTEST, rotTEST));
221	}
222
223
224	else if(i ==0){
225	rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); // SetTranslation(Double_t dx, Double_t dy, Double_t dz);
226	trTEST->SetTranslation(cos((77-endBeamParam)conv)(fEndBeamBeginRadius), sin((77-endBeamParam)conv)(fEndBeamBeginRadius),0);
227	endBeams->SetVisibility(1);
228	calHalfFrameMO->AddNode(endBeams, 1, new TGeoCombiTrans(trTEST, rotTEST));
229	}
230	else{
231	rotTEST->SetAngles(193.03, 0, 0); // SetTranslation(Double_t dx, Double_t dy, Double_t dz);
232	trTEST->SetTranslation(cos((193.03+endBeamParam)conv)(fEndBeamBeginRadius)/more duct tape/, sin((193.03+endBeamParam)conv)(fEndBeamBeginRadius),0);
233	endBeams->SetVisibility(1);
234	calHalfFrameMO->AddNode(endBeams, 2, new TGeoCombiTrans(trTEST, rotTEST));
235	}
236	}
237
238	//Beam Containers
239
240	// Translations
241
242	TGeoTranslation *origin1 = new TGeoTranslation(0,0,0); // Equivalent to gGeoIdentity
243	TGeoTranslation origin2 = new TGeoTranslation(0,0,2(fCrossBottomWidth+fFlangeWidth));
244	TGeoTranslation origin3 = new TGeoTranslation(0,0,-2(fCrossBottomWidth+fFlangeWidth));
245
246	// FlangeVolume translations
247	TGeoTranslation str1 = new TGeoTranslation(0,0,-3(fCrossBottomWidth+fFlangeWidth));
248	TGeoTranslation *str2 = new TGeoTranslation(0,0,-(fCrossBottomWidth+fFlangeWidth));
249	TGeoTranslation *str3 = new TGeoTranslation(0,0,(fCrossBottomWidth+fFlangeWidth));
250	TGeoTranslation str4 = new TGeoTranslation(0,0,3(fCrossBottomWidth+fFlangeWidth));
251
252	// Half Frame Translations
253	TGeoTranslation *halfTrans1 = new TGeoTranslation(0,0,fHalfFrameTrans);
254	TGeoTranslation *halfTrans2 = new TGeoTranslation(0,0,-fHalfFrameTrans);
255
256	// Beams Volume
257	calHalfFrameMO->AddNode(flangeVolume, 1, str1);
258	calHalfFrameMO->AddNode(flangeVolume, 2, str2);
259	calHalfFrameMO->AddNode(flangeVolume, 3, str3);
260	calHalfFrameMO->AddNode(flangeVolume, 4, str4);
261
262	calHalfFrameMO->AddNode(ribVolume, 1, str1);
263	calHalfFrameMO->AddNode(ribVolume, 2, str2);
264	calHalfFrameMO->AddNode(ribVolume, 3, str3);
265	calHalfFrameMO->AddNode(ribVolume, 4, str4);
266
267	// Cross Beams
268	calHalfFrameMO->AddNode(subSetCross, 1, origin1);
269	calHalfFrameMO->AddNode(subSetCross, 2, origin2);
270	calHalfFrameMO->AddNode(subSetCross, 3, origin3);
271	/* // Obsolete for now
272	calHalfFrameMO->AddNode(subSetCrossTop, 1, origin1);
273	calHalfFrameMO->AddNode(subSetCrossTop, 2, origin2);
274	calHalfFrameMO->AddNode(subSetCrossTop, 3, origin3);
275	*/
276
277	calHalfFrameMO->AddNode(skin, 1, gGeoIdentity);
278
279	calFrameMO->AddNode(calHalfFrameMO, 1, halfTrans1);
280	calFrameMO->AddNode(calHalfFrameMO, 2, halfTrans2);
281
282	top->AddNode(calFrameMO,1,gGeoIdentity);
283	// cout<<"**********************************\nfEndRadius:\t"<<fEndRadius;
284	}
285