EMCAL/AliEMCALSpaceFrame.cxx

   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+2*fCrossBottomRadThick-1, fBeginRadius+2*fCrossBottomRadThick+
 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 + i*fCrossBeamArcLength)*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