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