[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()
{
  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	{
139	AliDebug(1,"Create CalFrame Geometry");
140	//////////////////////////////////////Setup/////////////////////////////////////////
141	TGeoVolume* top = gGeoManager->GetVolume("ALIC");
142	TGeoMedium *steel = gGeoManager->GetMedium("EMCAL_S steel$");
143	TGeoMedium *air = gGeoManager->GetMedium("EMCAL_Air$");
144
145
146	//////////////////////////////////// Volumes ///////////////////////////////////////
147	TGeoVolume *calFrameMO =
148	gGeoManager->MakeTubs("CalFrame", air, fBeginRadius-2.1,fEndRadius,
149	fTotalHalfWidth*3,fBeginPhi-3,fEndPhi+3); // Mother Volume
150
151	calFrameMO->SetVisibility(kFALSE);
152
153	// Half Frame Mother Volume
154	TGeoVolume *calHalfFrameMO =
155	gGeoManager->MakeTubs("HalfFrame", air, fBeginRadius-2,fEndRadius,
156	fTotalHalfWidth,fBeginPhi-2.9,fEndPhi+2.9);
157
158	calHalfFrameMO->SetVisibility(kFALSE);
159
160	TGeoVolume *endBeams =
161	gGeoManager->MakeBox("End Beams", steel, fEndBeamRadThick, fCrossTopHeight, fTotalHalfWidth); // End Beams
162
163	TGeoVolume *skin =
164	gGeoManager->MakeTubs("skin", steel, fRibHeight+0.15, fEndRadius,
165	fTotalHalfWidth, fBeginPhi, fEndPhi);// back frame
166
167	TGeoVolume *flangeVolume =
168	gGeoManager->MakeTubs("supportBottom", steel, fBeginRadius, fFlangeHeight,
169	fFlangeWidth, fBeginPhi, fEndPhi); // FlangeVolume Beams
170
171	TGeoVolume *ribVolume =
172	gGeoManager->MakeTubs("RibVolume", steel, fFlangeHeight, fRibHeight, fRibWidth, fBeginPhi, fEndPhi);
173
174	TGeoVolume *subSetCross =
175	gGeoManager->MakeTubs("subSetCross", air, fBeginRadius-1, fBeginRadius+2*fCrossBottomRadThick+
176	2*fCrossTopRadThick+0.15, fCrossBottomWidth, fBeginPhi, fEndPhi); // Cross Beam Containers
177	subSetCross->SetVisibility(kFALSE);
178	/* // Obsolete for now
179	TGeoVolume *subSetCrossTop =
180	gGeoManager->MakeTubs("SubSetCrossTop", air, fBeginRadius+2fCrossBottomRadThick-1, fBeginRadius+2fCrossBottomRadThick+
181	2*fCrossTopRadThick+1, fCrossTopWidth, fBeginPhi, fEndPhi); // Cross
182	subSetCrossTop->SetVisibility(kFALSE);
183	*/
184	TGeoVolume *crossBottomBeams =
185	gGeoManager->MakeBox("crossBottom", steel, fCrossBottomRadThick, fCrossBottomHeight, fCrossBottomWidth); // Cross Beams
186
187	TGeoVolume *crossTopBeams =
188	gGeoManager->MakeBox("crossTop", steel, fCrossTopRadThick, fCrossTopHeight, fCrossTopWidth); // Cross Beams
189
190	TGeoTranslation *trTEST = new TGeoTranslation();
191	TGeoRotation *rotTEST = new TGeoRotation();
192
193	Double_t conv = TMath::Pi()/180.;
194	Double_t radAngle = 0;
195	Double_t endBeamParam=.4;
196	//cout<<"\nfCrossBottomStartRadius: "<<fCrossBottomStartRadius<<"\n";
197
198	for(Int_t i = 0; i < fNumCross; i++){
199
200	Double_t loopPhi = fBeginPhi + 1.8;
201
202	// Cross Bottom Beams
203
204	radAngle = (loopPhi + ifCrossBeamArcLength)conv;
205
206	rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); // SetTranslation(Double_t dx, Double_t dy, Double_t dz);
207	trTEST->SetTranslation(cos(radAngle)fCrossBottomStartRadius, sin(radAngle)fCrossBottomStartRadius,0);
208
209	TGeoCombiTrans combo = new TGeoCombiTrans(trTEST, *rotTEST); // TGeoTranslation &tr, const TGeoRotation &rot);
210	combo->RegisterYourself();
211	crossBottomBeams->SetVisibility(1);
212	subSetCross->AddNode(crossBottomBeams, i+1, combo);
213	if (i != 0 && i!=fNumCross-1){
214	// Cross Bottom Beams
215	rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); // SetTranslation(Double_t dx, Double_t dy, Double_t dz);
216	trTEST->SetTranslation(cos(radAngle)fCrossTopStart, sin(radAngle)fCrossTopStart,0);
217	crossTopBeams->SetVisibility(1);
218	subSetCross->AddNode(crossTopBeams, i+1, new TGeoCombiTrans(trTEST, rotTEST));
219	}
220
221
222	else if(i ==0){
223	rotTEST->SetAngles(fBeginPhi + i*fCrossBeamArcLength, 0, 0); // SetTranslation(Double_t dx, Double_t dy, Double_t dz);
224	trTEST->SetTranslation(cos((77-endBeamParam)conv)(fEndBeamBeginRadius), sin((77-endBeamParam)conv)(fEndBeamBeginRadius),0);
225	endBeams->SetVisibility(1);
226	calHalfFrameMO->AddNode(endBeams, 1, new TGeoCombiTrans(trTEST, rotTEST));
227	}
228	else{
229	rotTEST->SetAngles(193.03, 0, 0); // SetTranslation(Double_t dx, Double_t dy, Double_t dz);
230	trTEST->SetTranslation(cos((193.03+endBeamParam)conv)(fEndBeamBeginRadius)/more duct tape/, sin((193.03+endBeamParam)conv)(fEndBeamBeginRadius),0);
231	endBeams->SetVisibility(1);
232	calHalfFrameMO->AddNode(endBeams, 2, new TGeoCombiTrans(trTEST, rotTEST));
233	}
234	}
235
236	//Beam Containers
237
238	// Translations
239
240	TGeoTranslation *origin1 = new TGeoTranslation(0,0,0); // Equivalent to gGeoIdentity
241	TGeoTranslation origin2 = new TGeoTranslation(0,0,2(fCrossBottomWidth+fFlangeWidth));
242	TGeoTranslation origin3 = new TGeoTranslation(0,0,-2(fCrossBottomWidth+fFlangeWidth));
243
244	// FlangeVolume translations
245	TGeoTranslation str1 = new TGeoTranslation(0,0,-3(fCrossBottomWidth+fFlangeWidth));
246	TGeoTranslation *str2 = new TGeoTranslation(0,0,-(fCrossBottomWidth+fFlangeWidth));
247	TGeoTranslation *str3 = new TGeoTranslation(0,0,(fCrossBottomWidth+fFlangeWidth));
248	TGeoTranslation str4 = new TGeoTranslation(0,0,3(fCrossBottomWidth+fFlangeWidth));
249
250	// Half Frame Translations
251	TGeoTranslation *halfTrans1 = new TGeoTranslation(0,0,fHalfFrameTrans);
252	TGeoTranslation *halfTrans2 = new TGeoTranslation(0,0,-fHalfFrameTrans);
253
254	// Beams Volume
255	calHalfFrameMO->AddNode(flangeVolume, 1, str1);
256	calHalfFrameMO->AddNode(flangeVolume, 2, str2);
257	calHalfFrameMO->AddNode(flangeVolume, 3, str3);
258	calHalfFrameMO->AddNode(flangeVolume, 4, str4);
259
260	calHalfFrameMO->AddNode(ribVolume, 1, str1);
261	calHalfFrameMO->AddNode(ribVolume, 2, str2);
262	calHalfFrameMO->AddNode(ribVolume, 3, str3);
263	calHalfFrameMO->AddNode(ribVolume, 4, str4);
264
265	// Cross Beams
266	calHalfFrameMO->AddNode(subSetCross, 1, origin1);
267	calHalfFrameMO->AddNode(subSetCross, 2, origin2);
268	calHalfFrameMO->AddNode(subSetCross, 3, origin3);
269	/* // Obsolete for now
270	calHalfFrameMO->AddNode(subSetCrossTop, 1, origin1);
271	calHalfFrameMO->AddNode(subSetCrossTop, 2, origin2);
272	calHalfFrameMO->AddNode(subSetCrossTop, 3, origin3);
273	*/
274
275	calHalfFrameMO->AddNode(skin, 1, gGeoIdentity);
276
277	calFrameMO->AddNode(calHalfFrameMO, 1, halfTrans1);
278	calFrameMO->AddNode(calHalfFrameMO, 2, halfTrans2);
279
280	top->AddNode(calFrameMO,1,gGeoIdentity);
281	// cout<<"**********************************\nfEndRadius:\t"<<fEndRadius;
282	}
283