createRootGeoFileBwEnd_2017.C File Reference

Go to the source code of this file.

Functions
void	createRootGeoFileBwEnd_2017 ()

Variables
static const Double_t	mm = 1.0e-1
static const Double_t	kMPThickness = 30.0*mm
static const Double_t	kMPInsDist = 86.0*mm
static const Double_t	kInsertLength = 45.0*mm
static const Double_t	kCapsuleLength = 42.0*mm
static const Double_t	kCrLength = 200.0*mm
static const Double_t	kAlvLength = kCrLength + kCapsuleLength
static const Double_t	kCrystFrontCoolingDist = 6.0*mm
static const Double_t	kFrontCoolingThickness = 15.0*mm
static const Double_t	kFrontVIPThickness = 10.0*mm
static const Double_t	kCoverThickness = 4.0*mm
static const Double_t	kInnerVolLength
static const Double_t	kBWECzLength = kMPThickness + kInnerVolLength + kCoverThickness
static const Double_t	kInnerVolZpos = 0.5*(kMPThickness-kCoverThickness)
static const Double_t	kModuleZpos
static const Double_t	kFrontVIPZpos = 0.5*(kBWECzLength - kMPThickness - kCoverThickness - kFrontVIPThickness)
static const Double_t	kSideVIPZpos = 0.0
static const Double_t	kFrontCoolingZpos = 0.5*(kBWECzLength - kMPThickness - kCoverThickness - 2.0*kFrontVIPThickness - kFrontCoolingThickness)
static const Double_t	kCrWidth = 24.4*mm
static const Double_t	kCrHeight = 24.4*mm
static const Double_t	kAlvHoleWidth = kCrWidth + 0.5*mm
static const Double_t	kAlvWallThickness = 0.3*mm
static const Double_t	k4x4AlvWidth = 4.*kAlvHoleWidth + 5.*kAlvWallThickness
static const Double_t	kSingleAlvWidth = kAlvHoleWidth + 2.*kAlvWallThickness
static const Double_t	kInsAlvDiff = 0.45*mm
static const Double_t	kInsertWidth = k4x4AlvWidth - 2.*kInsAlvDiff
static const Double_t	kInsHoleWidth1 = 0.25*kInsertWidth - 2.*5.0*mm
static const Double_t	kInsHoleWidth2 = 0.50*kInsertWidth - 2.*5.0*mm
static const Double_t	firstCrystPosX = -1.5*(kAlvHoleWidth + kAlvWallThickness)
static const Double_t	firstCrystPosY = 1.5*(kAlvHoleWidth + kAlvWallThickness)
static const Double_t	crystDPos = kAlvHoleWidth + kAlvWallThickness
const Int_t	kNSubmoduleType = 6
static const Int_t	AlvLayout0 []
static const Int_t	AlvLayout1 []
static const Int_t	AlvLayout2 []
static const Int_t	AlvLayout3 []
static const Int_t	AlvLayout4 []
static const Int_t	AlvLayout5 []
static const Int_t	nVtx0 = 4
static const Double_t	v0X [nVtx0]
static const Double_t	v0Y [nVtx0]
static const Int_t	nVtx1 = 6
static const Double_t	v1X [nVtx1]
static const Double_t	v1Y [nVtx1]
static const Int_t	nVtx2 = 4
static const Double_t	v2X [nVtx2]
static const Double_t	v2Y [nVtx2]
static const Int_t	nVtx3 = 6
static const Double_t	v3X [nVtx3]
static const Double_t	v3Y [nVtx3]
static const Int_t	nVtx4 = 8
static const Double_t	v4X [nVtx4]
static const Double_t	v4Y [nVtx4]
static const Int_t	nVtx5 = 4
static const Double_t	v5X [nVtx5]
static const Double_t	v5Y [nVtx5]
static const Int_t	kNsubmodQuarter = 10
static const Double_t	kOuterSideR = 372.0*mm
static const Double_t	kDistBtwAlvs = 0.55*mm
static const Double_t	kSubmodPosUnit = k4x4AlvWidth + kDistBtwAlvs
static const Double_t	kCornerShift = 16.4*mm
static const Double_t	kSubmodPos0 [3] = {-1.*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth , kModuleZpos}
static const Double_t	kSubmodPos1 [3] = {0.0*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth , kModuleZpos}
static const Double_t	kSubmodPos2 [3] = {+1.*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth , kModuleZpos}
static const Double_t	kSubmodPos3 [3] = {-1.*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs , kModuleZpos}
static const Double_t	kSubmodPos4 [3] = {0.0*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs , kModuleZpos}
static const Double_t	kSubmodPos5 [3] = {+1.*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs , kModuleZpos}
static const Double_t	kSubmodPos6 [3] = {+2.*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs - kCornerShift, kModuleZpos}
static const Double_t	kSubmodPos7 [3] = {+3.*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs - kCornerShift, kModuleZpos}
static const Double_t	kSubmodPos8 [3] = {+2.*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth - kCornerShift, kModuleZpos}
static const Double_t	kSubmodPos9 [3] = {+3.*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth - kCornerShift, kModuleZpos}
static const Double_t	kBWECinnerRadius = 165.0*mm
static const Double_t	kBWECouterRadius = 438.0*mm
static const Double_t	kBWECxyWidth = 401.6*mm
static const Double_t	kSideVIPThickness = 10.0*mm
static const Double_t	kSideVIPRout = kBWECouterRadius-kCoverThickness
static const Double_t	kSideVIPRin = kSideVIPRout - kSideVIPThickness
static const Double_t	kSideVIPxyOuterWidth = kBWECxyWidth - kCoverThickness
static const Double_t	kSideVIPxyInnerWidth = kSideVIPxyOuterWidth - kSideVIPThickness
static const Double_t	kSideVIPLength = kInnerVolLength
static const Double_t	kInnerVIPRin = kBWECinnerRadius+kCoverThickness
static const Double_t	kInnerVIPRout = kInnerVIPRin + kSideVIPThickness
static const Double_t	kSideCoolingThickness = 15.0*mm

Function Documentation

void createRootGeoFileBwEnd_2017

(

)

Definition at line 258 of file createRootGeoFileBwEnd_2017.C.

References AlvLayout0, AlvLayout1, AlvLayout2, AlvLayout3, AlvLayout4, Double_t, geobuild, geoFace, geoLoad, gGeoMan, gGeoManager, i, kAlvLength, kBWECzLength, kCapsuleLength, kCrLength, kCrystFrontCoolingDist, kFrontCoolingThickness, kFrontVIPThickness, kInnerVolLength, kInnerVolZpos, kInsAlvDiff, kInsertLength, kModuleZpos, kNsubmodQuarter, kNSubmoduleType, kSubmodPos0, kSubmodPos1, kSubmodPos2, kSubmodPos3, kSubmodPos4, kSubmodPos5, kSubmodPos6, kSubmodPos7, kSubmodPos8, Media, mm, nmed, nVtx0, nVtx1, nVtx2, nVtx3, nVtx4, outfile, printf(), RegisterYourself(), rotName, SetLineColor(), SetLineWidth(), shape, t1, top, TString, v0X, v0Y, v1X, v1Y, v2X, v2Y, v3X, v3Y, v4X, and v4Y.

                                  {


  //----------------------------- initialisations
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libPndData");
  gSystem->Load("libPassive");
 

  
  //----------------------------- load/create materials/media
  FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader");
  FairGeoInterface *geoFace = geoLoad->getGeoInterface();
  geoFace->setMediaFile("/home/luigi/him/sw/pandaroot6/fairsoft_may16p1/pandaroot/geometry/media_pnd.geo");
  geoFace->readMedia();
  geoFace->print();

  FairGeoMedia *Media =  geoFace->getMedia();
  FairGeoBuilder *geobuild=geoLoad->getGeoBuilder();

  FairGeoMedium *N2      = Media->getMedium("Nitrogen");
  FairGeoMedium *PWO     = Media->getMedium("PWO");
  FairGeoMedium *Prepreg = Media->getMedium("Prepreg");
  FairGeoMedium *Alu     = Media->getMedium("aluminium");
  FairGeoMedium *VIPmat  = Media->getMedium("FumedSilica");
  FairGeoMedium *CSmat   = Media->getMedium("CoolingShellMix");

  Int_t nmed=geobuild->createMedium(N2);
  nmed=geobuild->createMedium(PWO);
  nmed=geobuild->createMedium(Prepreg);
  nmed=geobuild->createMedium(Alu);
  nmed=geobuild->createMedium(VIPmat);
  nmed=geobuild->createMedium(CSmat);

  
  //----------------------------- define mother volume
  TGeoManager* gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom");

  TGeoVolume *top = new TGeoVolumeAssembly("Emc4");
                    
  gGeoMan->SetTopVolume(top);


  //----------------------------- outer/inner BWEC volume
  TGeoTube *outerTube = new TGeoTube("outerTube",kBWECinnerRadius,kBWECouterRadius,0.5*kBWECzLength);
  TGeoBBox *outerBox = new TGeoBBox("outerBox",kBWECxyWidth,kBWECxyWidth,0.5*kBWECzLength+0.1);
  TGeoCompositeShape *outerVol = new TGeoCompositeShape("outerVol","outerTube*outerBox");

  TGeoVolume* BWECouterVol = new TGeoVolume("BWECouterVol",outerVol,gGeoMan->GetMedium("aluminium"));
  BWECouterVol->SetLineColor(kGray);
  BWECouterVol->SetLineWidth(5);
  BWECouterVol->SetLineStyle(kDashed);
  cout<<"transp: "<<BWECouterVol->GetTransparency()<<endl;

  top->AddNode(BWECouterVol,0);


  TGeoTube *innerTube = new TGeoTube("innerTube",
                                     kBWECinnerRadius+kCoverThickness,
                                     kBWECouterRadius-kCoverThickness,
                                     0.5*kInnerVolLength);
                                     
  TGeoBBox *innerBox = new TGeoBBox("innerBox",
                                    kBWECxyWidth-kCoverThickness,
                                    kBWECxyWidth-kCoverThickness,
                                    0.5*kInnerVolLength+0.1);
  
  TGeoCompositeShape *innerVol = new TGeoCompositeShape("innerVol","innerTube*innerBox");
  TGeoTranslation *t1 = new TGeoTranslation(0.0,0.0,kInnerVolZpos);
  
  TGeoVolume* BWECinnerVol = new TGeoVolume("BWECinnerVol",innerVol,gGeoMan->GetMedium("Nitrogen"));
  BWECinnerVol->SetLineColor(kAzure+10);

  BWECouterVol->AddNode(BWECinnerVol,0,t1);




  //----------------------------- submodules with alveoles, crystals, capsules...
  //--- version with composite shapes for the alveoles (holes are subtracted volumes)
  int crystalNo = 0;
  int holeNo = 0;
  int capsuleNo = 0;


  //--- submodule type 0: 4x4 array
  TGeoBBox *alveoleBox = new TGeoBBox("alveoleBox",0.5*k4x4AlvWidth,0.5*k4x4AlvWidth,0.5*kAlvLength);
  const Double_t *vX[kNSubmoduleType] = {v0X,v1X,v2X,v3X,v4X,v5X};
  const Double_t *vY[kNSubmoduleType] = {v0Y,v1Y,v2Y,v3Y,v4Y,v5Y};
  
  //--- submodule type 1
  TGeoXtru *alvXtru1 = new TGeoXtru(2); alvXtru1->SetName("alvXtru1");
  alvXtru1->DefinePolygon(nVtx1,v1X,v1Y);
  alvXtru1->DefineSection(0,-0.5*kAlvLength,0.0,0.0,1.0);
  alvXtru1->DefineSection(1,+0.5*kAlvLength,0.0,0.0,1.0);
  
  //--- submodule type 2
  TGeoXtru *alvXtru2 = new TGeoXtru(2); alvXtru2->SetName("alvXtru2");
  alvXtru2->DefinePolygon(nVtx2,v2X,v2Y);
  alvXtru2->DefineSection(0,-0.5*kAlvLength,0.0,0.0,1.0);
  alvXtru2->DefineSection(1,+0.5*kAlvLength,0.0,0.0,1.0);

  //--- submodule type 3
  TGeoXtru *alvXtru3 = new TGeoXtru(2); alvXtru3->SetName("alvXtru3");
  alvXtru3->DefinePolygon(nVtx3,v3X,v3Y);
  alvXtru3->DefineSection(0,-0.5*kAlvLength,0.0,0.0,1.0);
  alvXtru3->DefineSection(1,+0.5*kAlvLength,0.0,0.0,1.0);

  //--- submodule type 4
  TGeoXtru *alvXtru4 = new TGeoXtru(2); alvXtru4->SetName("alvXtru4");
  alvXtru4->DefinePolygon(nVtx4,v4X,v4Y);
  alvXtru4->DefineSection(0,-0.5*kAlvLength,0.0,0.0,1.0);
  alvXtru4->DefineSection(1,+0.5*kAlvLength,0.0,0.0,1.0);

  //--- "loose" crystals in the corners or submodule type 5
  TGeoXtru *alvXtru5 = new TGeoXtru(2); alvXtru5->SetName("alvXtru5");
  alvXtru5->DefinePolygon(nVtx5,v5X,v5Y);
  alvXtru5->DefineSection(0,-0.5*kAlvLength,0.0,0.0,1.0);
  alvXtru5->DefineSection(1,+0.5*kAlvLength,0.0,0.0,1.0);

  
  //--- shape of alveole hole to be subtracted from the bulk alveole shapes
  TGeoBBox *alvHoleBox = new TGeoBBox("alvHoleBox",0.5*kAlvHoleWidth,0.5*kAlvHoleWidth,0.5*kAlvLength+0.1*mm);

  //--- gathering shapes and layout maps (AlvLayout#) into arrays
  const Int_t *AlvLayout[kNSubmoduleType] = {AlvLayout0,AlvLayout1,AlvLayout2,AlvLayout3,AlvLayout4,AlvLayout5};
  TGeoShape *alveoleBulkShape[]={alveoleBox,alvXtru1,alvXtru2,alvXtru3,alvXtru4,alvXtru5};
  TGeoCompositeShape *alveoleShape[kNSubmoduleType];

  TString tname,holeVolFormula[kNSubmoduleType];
  for(Int_t i=0; i<kNSubmoduleType; ++i) holeVolFormula[i] = "";
  TGeoTranslation* holeInAlvXYtrans;

  //--- definition of subtraction formulas for the alveole types
  for(Int_t i=0; i<16;++i){
    tname = "t"; tname += i;
    holeInAlvXYtrans = new TGeoTranslation(tname,firstCrystPosX + (i%4)*crystDPos, firstCrystPosY - (i/4)*crystDPos, 0.0);
    holeInAlvXYtrans->RegisterYourself();
    for(Int_t k=0; k<kNSubmoduleType; ++k){
      if(AlvLayout[k][i]){
        if(holeVolFormula[k] != "") holeVolFormula[k] += " + ";
        holeVolFormula[k] +=  "alvHoleBox:"+tname;
      }
    }
  }

  //--- subtraction of alveole holes
  for(Int_t i=0; i<kNSubmoduleType; ++i){
    TString shapeName = "alveoleShape"; shapeName += i;
    TString shapeFormula = TString(alveoleBulkShape[i]->GetName()) + " - (" +  holeVolFormula[i] + ")";
    alveoleShape[i] = new TGeoCompositeShape(shapeName,shapeFormula);
    cout<<"Subtraction formula for alveole type "<<i<<": "<<shapeFormula<<endl;
  }

  
  //--- the alveole volumes
  TGeoVolume* BWECalveole[kNSubmoduleType];
  for(Int_t i=0; i<kNSubmoduleType; ++i){
    TString volName = "BWECalveole"; volName += i;
    BWECalveole[i] = new TGeoVolume(volName, alveoleShape[i], gGeoMan->GetMedium("Prepreg"));
    BWECalveole[i]->SetLineColor(kGray+3);
    BWECalveole[i]->SetFillStyle(2001);
  }


  //--- the crystal and Al capsule volumes and attributes
  TGeoBBox *crystalBox = new TGeoBBox("crystalBox",0.5*kCrWidth,0.5*kCrHeight,0.5*kCrLength);
  TGeoBBox *capsuleBox = new TGeoBBox("capsuleBox",0.5*kCrWidth,0.5*kCrHeight,0.5*kCapsuleLength);

  TGeoVolume* PWOcrystal    = new TGeoVolume("PWOCrystal",crystalBox,gGeoMan->GetMedium("PWO")); 
  TGeoVolume* crystCapsule  = new TGeoVolume("crystCapsule",capsuleBox,gGeoMan->GetMedium("aluminium"));

  TGeoTranslation* crystInAlvZtrans   = new TGeoTranslation(0.0,0.0,0.5*kCapsuleLength);
  TGeoTranslation* capsuleInAlvZtrans = new TGeoTranslation(0.0,0.0,-0.5*kCrLength);

  PWOcrystal     -> SetLineColor(kWhite);
  crystCapsule   -> SetLineColor(kGray);


  //--- Al insert vol
  TGeoBBox *insertBox  = new TGeoBBox("insertBox", 0.5*kInsertWidth,0.5*kInsertWidth,0.5*kInsertLength); // 4x4 alveole bulk
  TGeoBBox *insHoleBox1 = new TGeoBBox("insHoleBox1",0.5*kInsHoleWidth1,0.5*kInsHoleWidth1,0.5*kInsertLength+0.1); // 1x1 hole
  TGeoBBox *insHoleBox2 = new TGeoBBox("insHoleBox2",0.5*kInsHoleWidth2,0.5*kInsHoleWidth1,0.5*kInsertLength+0.1); // 2x1 hole
  TGeoBBox *insHoleBox4 = new TGeoBBox("insHoleBox4",0.5*kInsHoleWidth2,0.5*kInsHoleWidth2,0.5*kInsertLength+0.1); // 2x2 hole

  Double_t insVX[kNSubmoduleType][16],insVY[kNSubmoduleType][16];
  Int_t nVtx[kNSubmoduleType] = {nVtx0,nVtx1,nVtx2,nVtx3,nVtx4,nVtx5};
  TGeoXtru* insertXtru[kNSubmoduleType];
  TGeoCompositeShape* insertShape[kNSubmoduleType];

  // Translations for holes and composition formulas
  const Int_t nInsHoleTrans = 8;
  TGeoTranslation* insHoleTrans[nInsHoleTrans];
  insHoleTrans[0] = new TGeoTranslation("insHtr1",-0.25*kInsertWidth,+0.25*kInsertWidth,0.0);   // 1: 2x2 up-left
  insHoleTrans[1] = new TGeoTranslation("insHtr2",+0.25*kInsertWidth,+0.25*kInsertWidth,0.0);   // 2: 2x2 up-right
  insHoleTrans[2] = new TGeoTranslation("insHtr3",-0.25*kInsertWidth,-0.25*kInsertWidth,0.0);   // 3: 2x2 down-left
  insHoleTrans[3] = new TGeoTranslation("insHtr4",+0.25*kInsertWidth,-0.25*kInsertWidth,0.0);   // 4: 2x2 down-right
  insHoleTrans[4] = new TGeoTranslation("insHtr5",-0.25*kInsertWidth,-0.125*kInsertWidth,0.0);  // 5: 2x1 down-left
  insHoleTrans[5] = new TGeoTranslation("insHtr6",+0.25*kInsertWidth,-0.125*kInsertWidth,0.0);  // 6: 2x1 down-right
  insHoleTrans[6] = new TGeoTranslation("insHtr7",+0.125*kInsertWidth,+0.125*kInsertWidth,0.0); // 7: 1x1 up-right (middle)
  insHoleTrans[7] = new TGeoTranslation("insHtr8",-0.375*kInsertWidth,-0.375*kInsertWidth,0.0); // 8: 1x1 down-left (corner)
  for(Int_t i=0;i<nInsHoleTrans;++i) insHoleTrans[i]->RegisterYourself();
  TString insSubFormula[kNSubmoduleType];
  insSubFormula[0] = "insertBox   - insHoleBox4:insHtr1 - insHoleBox4:insHtr2 - insHoleBox4:insHtr3 - insHoleBox4:insHtr4";
  insSubFormula[1] = "insertXtru1 - insHoleBox4:insHtr1 - insHoleBox4:insHtr2 - insHoleBox4:insHtr3 - insHoleBox2:insHtr6";
  insSubFormula[2] = "insertXtru2 - insHoleBox4:insHtr1 - insHoleBox4:insHtr2 - insHoleBox2:insHtr5 - insHoleBox2:insHtr6";
  insSubFormula[3] = "insertXtru3 - insHoleBox4:insHtr1 - insHoleBox4:insHtr2 - insHoleBox2:insHtr5 - insHoleBox4:insHtr4";
  insSubFormula[4] = "insertXtru4 - insHoleBox4:insHtr1 - insHoleBox1:insHtr7 - insHoleBox4:insHtr3 - insHoleBox4:insHtr4";
  insSubFormula[5] = "insertXtru5 - insHoleBox1:insHtr8";


  insertShape[0] = new TGeoCompositeShape("insertShape0",insSubFormula[0]);
  TGeoVolume* insert[kNSubmoduleType];
  insert[0] = new TGeoVolume("insert0",insertShape[0],gGeoMan->GetMedium("aluminium"));
  insert[0]->SetLineColor(kYellow+1);

  for(Int_t subi=1;subi<kNSubmoduleType;++subi){
    for(Int_t i=0; i<nVtx[subi]; ++i){
      insVX[subi][i] = vX[subi][i] + (vX[subi][i]<0 ? 1.0:-1.0)*kInsAlvDiff;
      insVY[subi][i] = vY[subi][i] + (vY[subi][i]<0 ? 1.0:-1.0)*kInsAlvDiff;
    }
    TString shapeName = "insertXtru"; shapeName += subi;
    insertXtru[subi] = new TGeoXtru(2); insertXtru[subi]->SetName(shapeName);
    insertXtru[subi]->DefinePolygon(nVtx[subi],insVX[subi],insVY[subi]);
    insertXtru[subi]->DefineSection(0,-0.5*kInsertLength,0.0,0.0,1.0);
    insertXtru[subi]->DefineSection(1,+0.5*kInsertLength,0.0,0.0,1.0);

    shapeName = "insertShape"; shapeName += subi;
    insertShape[subi] = new TGeoCompositeShape(shapeName,insSubFormula[subi]);
    TString volName = "insert"; volName += subi;
    insert[subi] = new TGeoVolume(volName,insertShape[subi],gGeoMan->GetMedium("aluminium"));
    insert[subi]->SetLineColor(kYellow+1);
  }
  TGeoTranslation* insertZtrans = new TGeoTranslation(0.0,0.0,-0.5*(kAlvLength+kInsertLength));



  //--- volume assemblies describing the BWEC submodules
  TGeoVolumeAssembly* BWECsubmodule[kNSubmoduleType];
  for(Int_t i=0; i<kNSubmoduleType; ++i){
    TString volName = "BWECsubmodule"; volName += i;
    BWECsubmodule[i] = new TGeoVolumeAssembly(volName);

    //--- adding the carbon fibre alveole
    BWECsubmodule[i]->AddNode(BWECalveole[i],i);

    //--- adding the crystals and Al capsules
    for(Int_t icr=0; icr<16;++icr){
      TGeoTranslation *crystTrans = new TGeoTranslation(firstCrystPosX + (icr%4)*crystDPos, firstCrystPosY - (icr/4)*crystDPos, 0.5*kCapsuleLength);
      TGeoTranslation *capslTrans = new TGeoTranslation(firstCrystPosX + (icr%4)*crystDPos, firstCrystPosY - (icr/4)*crystDPos,-0.5*kCrLength);
      if( AlvLayout[i][icr] ){
        BWECsubmodule[i]->AddNode(PWOcrystal  ,icr, crystTrans);
        BWECsubmodule[i]->AddNode(crystCapsule,icr, capslTrans);
      }
    }
    
    //--- adding Al insert
    BWECsubmodule[i]->AddNode(insert[i], 0, insertZtrans);

  }



  //----------------------------- Quarters of the BWEC: side (6 submodules) + corner (4 submodules)
  TGeoVolumeAssembly *BWECquarter = new TGeoVolumeAssembly("BWECquarter");

  // side
  BWECquarter->AddNode(BWECsubmodule[0],0,new TGeoTranslation(kSubmodPos0[0],kSubmodPos0[1],kSubmodPos0[2]));
  BWECquarter->AddNode(BWECsubmodule[0],1,new TGeoTranslation(kSubmodPos1[0],kSubmodPos1[1],kSubmodPos1[2]));
  BWECquarter->AddNode(BWECsubmodule[0],2,new TGeoTranslation(kSubmodPos2[0],kSubmodPos2[1],kSubmodPos2[2]));
  BWECquarter->AddNode(BWECsubmodule[1],3,new TGeoTranslation(kSubmodPos3[0],kSubmodPos3[1],kSubmodPos3[2]));
  BWECquarter->AddNode(BWECsubmodule[2],4,new TGeoTranslation(kSubmodPos4[0],kSubmodPos4[1],kSubmodPos4[2]));
  BWECquarter->AddNode(BWECsubmodule[3],5,new TGeoTranslation(kSubmodPos5[0],kSubmodPos5[1],kSubmodPos5[2]));

  // corner
  BWECquarter->AddNode(BWECsubmodule[0],8,new TGeoTranslation(kSubmodPos6[0],kSubmodPos6[1],kSubmodPos6[2]));
  BWECquarter->AddNode(BWECsubmodule[4],9,new TGeoTranslation(kSubmodPos7[0],kSubmodPos7[1],kSubmodPos7[2]));
  BWECquarter->AddNode(BWECsubmodule[4],6,new TGeoTranslation(kSubmodPos8[0],kSubmodPos8[1],kSubmodPos8[2]));
  BWECquarter->AddNode(BWECsubmodule[5],7,new TGeoTranslation(kSubmodPos9[0],kSubmodPos9[1],kSubmodPos9[2]));

  TGeoRotation *quarterRot[4];
  for(Int_t qi=0; qi<4; ++qi){
    quarterRot[qi] = new TGeoRotation;
    TString rotName("quarterRot"); rotName += qi;
    quarterRot[qi]->SetName(rotName);
    quarterRot[qi]->RotateZ(90.*qi);
    quarterRot[qi]->RegisterYourself();
    BWECinnerVol->AddNode(BWECquarter,qi,quarterRot[qi]);
  }



  //----------------------------- Dead material volumes: VIPs and cooling shells
  //--- front VIP
  TGeoTube *frontVIPTube = new TGeoTube("frontVIPTube",kInnerVIPRout,kSideVIPRin,0.5*kFrontVIPThickness);
                                     
  TGeoBBox *frontVIPrBox = new TGeoBBox("frontVIPBox",kSideVIPxyInnerWidth,kSideVIPxyInnerWidth,
                                        0.5*kFrontVIPThickness+0.1);
  TGeoCompositeShape* frontVIPshape = new TGeoCompositeShape("frontVIPshape","frontVIPTube*frontVIPBox");
  TGeoVolume* frontVIP = new TGeoVolume("frontVIP",frontVIPshape,gGeoMan->GetMedium("FumedSilica"));
  frontVIP->SetLineColor(kGray+1);
  frontVIP->SetLineStyle(kSolid);
  BWECinnerVol->AddNode(frontVIP,0,new TGeoTranslation(0.0,0.0,kFrontVIPZpos));


  //--- front Cooling
  TGeoTube *frontCoolingTube = new TGeoTube("frontCoolingTube",kInnerVIPRout,kSideVIPRin,0.5*kFrontCoolingThickness);
                                     
  TGeoBBox *frontCoolingBox = new TGeoBBox("frontCoolingBox",kSideVIPxyInnerWidth,kSideVIPxyInnerWidth,
                                            0.5*kFrontCoolingThickness+0.1);
  TGeoCompositeShape* frontCoolingshape = new TGeoCompositeShape("frontCoolingshape","frontCoolingTube*frontCoolingBox");
  TGeoVolume* frontCooling = new TGeoVolume("frontCooling",frontCoolingshape,gGeoMan->GetMedium("CoolingShellMix"));
  frontCooling->SetLineColor(kWhite);
  frontCooling->SetLineStyle(kSolid);
  BWECinnerVol->AddNode(frontCooling,0,new TGeoTranslation(0.0,0.0,kFrontCoolingZpos));


  //--- Side VIP
  TGeoTube* sideVIPTubeOut = new TGeoTube("sideVIPTubeOut",1.0,kSideVIPRout,0.5*kSideVIPLength);
  TGeoTube* sideVIPTubeIn  = new TGeoTube("sideVIPTubeIn", 0.0,kSideVIPRin, 0.5*kSideVIPLength+0.1);
  TGeoBBox *sideVIPBoxOut  = new TGeoBBox("sideVIPBoxOut",kSideVIPxyOuterWidth,kSideVIPxyOuterWidth,0.5*kSideVIPLength+0.11);
  TGeoBBox *sideVIPBoxIn   = new TGeoBBox("sideVIPBoxIn", kSideVIPxyInnerWidth,kSideVIPxyInnerWidth,0.5*kSideVIPLength+0.11);
  TGeoCompositeShape* sideVIPshape = 
    new TGeoCompositeShape("sideVIPshape","(sideVIPTubeOut*sideVIPBoxOut)-(sideVIPTubeIn*sideVIPBoxIn)" ); 
  TGeoVolume* sideVIP = new TGeoVolume("sideVIP",sideVIPshape,gGeoMan->GetMedium("FumedSilica"));
  sideVIP->SetLineColor(kGray+1);
  sideVIP->SetLineStyle(kSolid);
  BWECinnerVol->AddNode(sideVIP,0,new TGeoTranslation(0.0,0.0,kSideVIPZpos));

  //--- Inner VIP
  TGeoTube* innerVIPTube = new TGeoTube("sideVIPTubeOut",kInnerVIPRin,kInnerVIPRout,0.5*kSideVIPLength);
  TGeoVolume* innerVIP = new TGeoVolume("innerVIP",innerVIPTube,gGeoMan->GetMedium("FumedSilica"));
  innerVIP->SetLineColor(kGray+1);
  innerVIP->SetLineStyle(kSolid);
  BWECinnerVol->AddNode(innerVIP,0,new TGeoTranslation(0.0,0.0,kSideVIPZpos));



  //--- Side and inner Cooling
  TGeoTube *sideCoolingTube = new TGeoTube("sideCoolingTube",kInnerVIPRout,kSideVIPRin,0.5*(kAlvLength+kInsertLength));
  TGeoBBox *sideCoolingBox = new TGeoBBox("sideCoolingBox",kSideVIPxyInnerWidth,kSideVIPxyInnerWidth,0.5*(kAlvLength+kInsertLength)+0.1);

  //--- Complex shape for the hole. Building it up from the data for the alveole Xtru objects.
  Double_t shVX[kNSubmoduleType][16],shVY[kNSubmoduleType][16];
  TGeoXtru *shellHXtru[kNSubmoduleType];

  for(Int_t subi=0;subi<kNSubmoduleType;++subi){
    for(Int_t i=0; i<nVtx[subi]; ++i){
      shVX[subi][i] = vX[subi][i] + ((vX[subi][i]<0) ? -1.0:1.0)*(kDistBtwAlvs+(subi==0?0.5*mm:0.1*mm));
      shVY[subi][i] = vY[subi][i] + ((vY[subi][i]<0) ? -1.0:1.0)*(kDistBtwAlvs+(subi==0?0.5*mm:0.1*mm));
    }
    if(subi == 5) { // special case for the "loose" crystal
      shVX[subi][1] = vX[subi][1] + kDistBtwAlvs + 0.1*mm;
      shVX[subi][2] = vX[subi][2] + kDistBtwAlvs + 0.1*mm;  shVY[subi][2] = vY[subi][2] + kDistBtwAlvs + 0.1*mm;
      shVY[subi][3] = vY[subi][3] + kDistBtwAlvs + 0.1*mm;
    }
    TString shapeName = "shXtru"; shapeName += subi;
    shellHXtru[subi] = new TGeoXtru(2); shellHXtru[subi]->SetName(shapeName);
    shellHXtru[subi]->DefinePolygon(nVtx[subi],shVX[subi],shVY[subi]);
    shellHXtru[subi]->DefineSection(0,-0.5*(kAlvLength+kInsertLength+2.),0.0,0.0,1.0);
    shellHXtru[subi]->DefineSection(1,+0.5*(kAlvLength+kInsertLength+2.),0.0,0.0,1.0);
  }

  TGeoTranslation *shTransl[kNsubmodQuarter];
  const Double_t *shTrV[kNsubmodQuarter] = {kSubmodPos0,kSubmodPos1,kSubmodPos2,kSubmodPos3,kSubmodPos4,
                                            kSubmodPos5,kSubmodPos6,kSubmodPos7,kSubmodPos8,kSubmodPos9};
  TString shXtruName[kNsubmodQuarter]={"shXtru0","shXtru0","shXtru0","shXtru1","shXtru2",
                                       "shXtru3","shXtru0","shXtru4","shXtru4","shXtru5"};

  TString unionFormula="";
  for(Int_t i=0; i<kNsubmodQuarter; ++i){
    TString trName = "shTr"; trName += i;
    shTransl[i] = new TGeoTranslation(trName,shTrV[i][0],shTrV[i][1],shTrV[i][2]-kInsertLength);
    shTransl[i]->RegisterYourself();
    unionFormula += (unionFormula=="")?"":" + ";
    unionFormula += shXtruName[i]+":"+trName;
  }
  cout<<unionFormula<<endl;
  TGeoCompositeShape* sideCoolingHoleShape = new TGeoCompositeShape("sideCoolHole",unionFormula);

  TString sideCoolingFormula = "sideCoolingTube*sideCoolingBox";
  for(Int_t qi=0; qi<4; ++qi) sideCoolingFormula += TString(" - sideCoolHole:")+quarterRot[qi]->GetName();
  cout<<sideCoolingFormula<<endl;

  TGeoCompositeShape* sideCoolingShape = new TGeoCompositeShape("sideCoolingShape",sideCoolingFormula);
  TGeoVolume* sideCooling = new TGeoVolume("sideCooling",sideCoolingShape,gGeoMan->GetMedium("aluminium"));
  sideCooling->SetLineColor(kWhite);
  sideCooling->SetLineStyle(kSolid);
  BWECinnerVol->AddNodeOverlap(sideCooling,0,new TGeoTranslation(0.0,0.0,kModuleZpos-0.5*kInsertLength));




  // checks:
  cout<<"check: "<<kInnerVolLength - kFrontVIPThickness - kFrontCoolingThickness - kCrystFrontCoolingDist - kCrLength - 173.*mm<<endl;
  cout<<"check: "<<0.5*kInnerVolLength - (kModuleZpos + 0.5*kAlvLength + kFrontVIPThickness + kFrontCoolingThickness)<<endl;;

  // to be used in PndEmc.cxx for calculating the position of the bw emc
  TGeoNode *outer=top->FindNode("BWECouterVol_0");
  if(outer != NULL){
    TGeoCompositeShape* shape=(TGeoCompositeShape*)outer->GetVolume()->GetShape();
    Double_t origin[3]={0,0,0}, zdir[3]={0,0,1};
    Double_t halfLength = shape->DistFromInside(origin,zdir);
    printf("halfLength = %e, fullLength = %e\n",halfLength,2.*halfLength);
  }


  
  gGeoMan->CloseGeometry();


  TString outfile = "/home/luigi/him/sw/pandaroot6/fairsoft_may16p1/pandaroot/macro/emc/dedicated/emc_module4_2017-tmp.root";
  TFile *fout = new TFile(outfile,"RECREATE"); 
  //top->SetDirectory(fout);
  top->Write();
  fout->Close();
  
  //gGeoMan->Export(outfile);

  

  //----------------------------- visibility
  PWOcrystal       -> SetVisibility(kFALSE);
  crystCapsule     -> SetVisibility(kTRUE);
  BWECalveole[0]   -> SetVisibility(kTRUE);
  BWECalveole[1]   -> SetVisibility(kTRUE);
  BWECalveole[2]   -> SetVisibility(kTRUE);
  BWECalveole[3]   -> SetVisibility(kTRUE);
  BWECalveole[4]   -> SetVisibility(kTRUE);
  BWECalveole[5]   -> SetVisibility(kTRUE);
  insert[0]        -> SetVisibility(kTRUE); 
  insert[1]        -> SetVisibility(kTRUE); 
  insert[2]        -> SetVisibility(kTRUE); 
  insert[3]        -> SetVisibility(kTRUE); 
  insert[4]        -> SetVisibility(kTRUE); 
  insert[5]        -> SetVisibility(kTRUE); 
  frontVIP         -> SetVisibility(kFALSE);
  sideVIP          -> SetVisibility(kFALSE);
  innerVIP         -> SetVisibility(kFALSE);
  frontCooling     -> SetVisibility(kFALSE);
  sideCooling      -> SetVisibility(kFALSE);
  BWECinnerVol     -> SetVisibility(kFALSE);
  BWECouterVol     -> SetVisibility(kFALSE);


  crystCapsule     -> SetLineWidth(1);
  BWECalveole[0]   -> SetLineWidth(1);
  BWECalveole[1]   -> SetLineWidth(1);
  BWECalveole[2]   -> SetLineWidth(1);
  BWECalveole[3]   -> SetLineWidth(1);
  BWECalveole[4]   -> SetLineWidth(1);
  BWECalveole[5]   -> SetLineWidth(1);
  insert[0]        -> SetLineWidth(1);
  insert[1]        -> SetLineWidth(1);
  insert[2]        -> SetLineWidth(1);
  insert[3]        -> SetLineWidth(1);
  insert[4]        -> SetLineWidth(1);
  insert[5]        -> SetLineWidth(1);
  frontVIP         -> SetLineWidth(1);
  sideVIP          -> SetLineWidth(1);
  innerVIP         -> SetLineWidth(1);
  frontCooling     -> SetLineWidth(1);
  sideCooling      -> SetLineWidth(1);
  BWECinnerVol     -> SetLineWidth(1);
  BWECouterVol     -> SetLineWidth(1);

  top->Raytrace();
  gStyle->SetCanvasPreferGL(kTRUE);
  gGeoManager->SetTopVisible();
  gGeoManager->SetVisOption(0);
  gGeoManager->SetVisLevel(11);
  //top->Draw("ogl");
  top->Draw("ogle");

   
  new TBrowser;

}

Variable Documentation

const Int_t AlvLayout0[]

static

Initial value:

={1,1,1,1,
                                 1,1,1,1,
                                 1,1,1,1,
                                 1,1,1,1}

Definition at line 84 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t AlvLayout1[]

static

Initial value:

={1,1,1,1,
                                 1,1,1,1,
                                 1,1,1,1,
                                 1,1,0,0}

Definition at line 89 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t AlvLayout2[]

static

Initial value:

={1,1,1,1,
                                 1,1,1,1,
                                 1,1,1,1,
                                 0,0,0,0}

Definition at line 94 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t AlvLayout3[]

static

Initial value:

={1,1,1,1,
                                 1,1,1,1,
                                 1,1,1,1,
                                 0,0,1,1}

Definition at line 99 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t AlvLayout4[]

static

Initial value:

={1,1,0,0,
                                 1,1,1,0,
                                 1,1,1,1,
                                 1,1,1,1}

Definition at line 104 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t AlvLayout5[]

static

Initial value:

={0,0,0,0,
                                 0,0,0,0,
                                 0,0,0,0,
                                 1,0,0,0}

Definition at line 109 of file createRootGeoFileBwEnd_2017.C.

const Double_t crystDPos = kAlvHoleWidth + kAlvWallThickness

static

Definition at line 74 of file createRootGeoFileBwEnd_2017.C.

const Double_t firstCrystPosX = -1.5*(kAlvHoleWidth + kAlvWallThickness)

static

Definition at line 72 of file createRootGeoFileBwEnd_2017.C.

const Double_t firstCrystPosY = 1.5*(kAlvHoleWidth + kAlvWallThickness)

static

Definition at line 73 of file createRootGeoFileBwEnd_2017.C.

const Double_t k4x4AlvWidth = 4.*kAlvHoleWidth + 5.*kAlvWallThickness

static

Definition at line 63 of file createRootGeoFileBwEnd_2017.C.

const Double_t kAlvHoleWidth = kCrWidth + 0.5*mm

static

Definition at line 61 of file createRootGeoFileBwEnd_2017.C.

const Double_t kAlvLength = kCrLength + kCapsuleLength

static

Definition at line 21 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kAlvWallThickness = 0.3*mm

static

Definition at line 62 of file createRootGeoFileBwEnd_2017.C.

const Double_t kBWECinnerRadius = 165.0*mm

static

Definition at line 228 of file createRootGeoFileBwEnd_2017.C.

const Double_t kBWECouterRadius = 438.0*mm

static

Definition at line 229 of file createRootGeoFileBwEnd_2017.C.

const Double_t kBWECxyWidth = 401.6*mm

static

Definition at line 230 of file createRootGeoFileBwEnd_2017.C.

const Double_t kBWECzLength = kMPThickness + kInnerVolLength + kCoverThickness

static

Definition at line 33 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kCapsuleLength = 42.0*mm

static

Definition at line 19 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kCornerShift = 16.4*mm

static

Definition at line 200 of file createRootGeoFileBwEnd_2017.C.

const Double_t kCoverThickness = 4.0*mm

static

Definition at line 25 of file createRootGeoFileBwEnd_2017.C.

const Double_t kCrHeight = 24.4*mm

static

Definition at line 59 of file createRootGeoFileBwEnd_2017.C.

const Double_t kCrLength = 200.0*mm

static

Definition at line 20 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kCrWidth = 24.4*mm

static

Definition at line 58 of file createRootGeoFileBwEnd_2017.C.

const Double_t kCrystFrontCoolingDist = 6.0*mm

static

Definition at line 22 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kDistBtwAlvs = 0.55*mm

static

Definition at line 198 of file createRootGeoFileBwEnd_2017.C.

const Double_t kFrontCoolingThickness = 15.0*mm

static

Definition at line 23 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kFrontCoolingZpos = 0.5*(kBWECzLength - kMPThickness - kCoverThickness - 2.0*kFrontVIPThickness - kFrontCoolingThickness)

static

Definition at line 47 of file createRootGeoFileBwEnd_2017.C.

const Double_t kFrontVIPThickness = 10.0*mm

static

Definition at line 24 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kFrontVIPZpos = 0.5*(kBWECzLength - kMPThickness - kCoverThickness - kFrontVIPThickness)

static

Definition at line 45 of file createRootGeoFileBwEnd_2017.C.

const Double_t kInnerVIPRin = kBWECinnerRadius+kCoverThickness

static

Definition at line 245 of file createRootGeoFileBwEnd_2017.C.

const Double_t kInnerVIPRout = kInnerVIPRin + kSideVIPThickness

static

Definition at line 246 of file createRootGeoFileBwEnd_2017.C.

const Double_t kInnerVolLength

static

Initial value:

=
  kMPInsDist + kInsertLength + kAlvLength + kCrystFrontCoolingDist + kFrontCoolingThickness + kFrontVIPThickness

Definition at line 29 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kInnerVolZpos = 0.5*(kMPThickness-kCoverThickness)

static

Definition at line 37 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kInsAlvDiff = 0.45*mm

static

Definition at line 66 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kInsertLength = 45.0*mm

static

Definition at line 18 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kInsertWidth = k4x4AlvWidth - 2.*kInsAlvDiff

static

Definition at line 67 of file createRootGeoFileBwEnd_2017.C.

const Double_t kInsHoleWidth1 = 0.25*kInsertWidth - 2.*5.0*mm

static

Definition at line 68 of file createRootGeoFileBwEnd_2017.C.

const Double_t kInsHoleWidth2 = 0.50*kInsertWidth - 2.*5.0*mm

static

Definition at line 69 of file createRootGeoFileBwEnd_2017.C.

const Double_t kModuleZpos

static

Initial value:

= 0.5*( (kMPInsDist + kInsertLength) -
                                          (kCrystFrontCoolingDist + kFrontCoolingThickness + kFrontVIPThickness) )

Definition at line 40 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kMPInsDist = 86.0*mm

static

Definition at line 17 of file createRootGeoFileBwEnd_2017.C.

const Double_t kMPThickness = 30.0*mm

static

Definition at line 16 of file createRootGeoFileBwEnd_2017.C.

const Int_t kNsubmodQuarter = 10

static

Definition at line 196 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t kNSubmoduleType = 6

Definition at line 80 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kOuterSideR = 372.0*mm

static

Definition at line 197 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSideCoolingThickness = 15.0*mm

static

Definition at line 248 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSideVIPLength = kInnerVolLength

static

Definition at line 243 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSideVIPRin = kSideVIPRout - kSideVIPThickness

static

Definition at line 240 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSideVIPRout = kBWECouterRadius-kCoverThickness

static

Definition at line 239 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSideVIPThickness = 10.0*mm

static

Definition at line 238 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSideVIPxyInnerWidth = kSideVIPxyOuterWidth - kSideVIPThickness

static

Definition at line 242 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSideVIPxyOuterWidth = kBWECxyWidth - kCoverThickness

static

Definition at line 241 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSideVIPZpos = 0.0

static

Definition at line 46 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSingleAlvWidth = kAlvHoleWidth + 2.*kAlvWallThickness

static

Definition at line 64 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSubmodPos0[3] = {-1.*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth , kModuleZpos}

static

Definition at line 203 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos1[3] = {0.0*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth , kModuleZpos}

static

Definition at line 204 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos2[3] = {+1.*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth , kModuleZpos}

static

Definition at line 205 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos3[3] = {-1.*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs , kModuleZpos}

static

Definition at line 206 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos4[3] = {0.0*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs , kModuleZpos}

static

Definition at line 207 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos5[3] = {+1.*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs , kModuleZpos}

static

Definition at line 208 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos6[3] = {+2.*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs - kCornerShift, kModuleZpos}

static

Definition at line 211 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos7[3] = {+3.*kSubmodPosUnit, kOuterSideR - 1.5*k4x4AlvWidth - kDistBtwAlvs - kCornerShift, kModuleZpos}

static

Definition at line 212 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos8[3] = {+2.*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth - kCornerShift, kModuleZpos}

static

Definition at line 213 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t kSubmodPos9[3] = {+3.*kSubmodPosUnit, kOuterSideR - 0.5*k4x4AlvWidth - kCornerShift, kModuleZpos}

static

Definition at line 214 of file createRootGeoFileBwEnd_2017.C.

const Double_t kSubmodPosUnit = k4x4AlvWidth + kDistBtwAlvs

static

Definition at line 199 of file createRootGeoFileBwEnd_2017.C.

const Double_t mm = 1.0e-1

static

Definition at line 11 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t nVtx0 = 4

static

Definition at line 119 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t nVtx1 = 6

static

Definition at line 129 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t nVtx2 = 4

static

Definition at line 143 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t nVtx3 = 6

static

Definition at line 153 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t nVtx4 = 8

static

Definition at line 166 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Int_t nVtx5 = 4

static

Definition at line 182 of file createRootGeoFileBwEnd_2017.C.

const Double_t v0X[nVtx0]

static

Initial value:

= {
  -0.5*k4x4AlvWidth, +0.5*k4x4AlvWidth, 
  +0.5*k4x4AlvWidth, -0.5*k4x4AlvWidth}

Definition at line 120 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v0Y[nVtx0]

static

Initial value:

= {
  -0.5*k4x4AlvWidth, -0.5*k4x4AlvWidth, 
  +0.5*k4x4AlvWidth, 0.5*k4x4AlvWidth}

Definition at line 123 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v1X[nVtx1]

static

Initial value:

= {
  -0.5*k4x4AlvWidth, 
  +0.5*kAlvWallThickness, 0.5*kAlvWallThickness, 
  +0.5*k4x4AlvWidth, 0.5*k4x4AlvWidth,
  -0.5*k4x4AlvWidth}

Definition at line 130 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v1Y[nVtx1]

static

Initial value:

= {
  -0.5*k4x4AlvWidth, -0.5*k4x4AlvWidth, 
  -0.5*k4x4AlvWidth + kAlvWallThickness + kAlvHoleWidth,
  -0.5*k4x4AlvWidth + kAlvWallThickness + kAlvHoleWidth,
  +0.5*k4x4AlvWidth, 0.5*k4x4AlvWidth}

Definition at line 135 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v2X[nVtx2]

static

Initial value:

= {
  -0.5*k4x4AlvWidth,  0.5*k4x4AlvWidth, 
  +0.5*k4x4AlvWidth, -0.5*k4x4AlvWidth}

Definition at line 144 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v2Y[nVtx2]

static

Initial value:

= {
  -0.5*k4x4AlvWidth + kAlvHoleWidth + kAlvWallThickness, 
  -0.5*k4x4AlvWidth + kAlvHoleWidth + kAlvWallThickness, 
  +0.5*k4x4AlvWidth,  0.5*k4x4AlvWidth}

Definition at line 147 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v3X[nVtx3]

static

Initial value:

= {
  -0.5*k4x4AlvWidth, 
  -0.5*kAlvWallThickness, -0.5*kAlvWallThickness, 
  +0.5*k4x4AlvWidth, 0.5*k4x4AlvWidth,
  -0.5*k4x4AlvWidth}

Definition at line 154 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v3Y[nVtx3]

static

Initial value:

= {
  -0.5*k4x4AlvWidth + kAlvWallThickness + kAlvHoleWidth,
  -0.5*k4x4AlvWidth + kAlvWallThickness + kAlvHoleWidth,
  -0.5*k4x4AlvWidth, -0.5*k4x4AlvWidth, 
  +0.5*k4x4AlvWidth,  0.5*k4x4AlvWidth}

Definition at line 159 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v4X[nVtx4]

static

Initial value:

= {
  -0.5*k4x4AlvWidth, 
  +0.5*k4x4AlvWidth, 0.5*k4x4AlvWidth,
  +0.5*k4x4AlvWidth - kAlvWallThickness - kAlvHoleWidth,
  +0.5*k4x4AlvWidth - kAlvWallThickness - kAlvHoleWidth,
  +0.5*kAlvWallThickness, 0.5*kAlvWallThickness, 
  -0.5*k4x4AlvWidth}

Definition at line 167 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v4Y[nVtx4]

static

Initial value:

= {
  -0.5*k4x4AlvWidth, -0.5*k4x4AlvWidth, 
  +0.5001*kAlvWallThickness, 0.5001*kAlvWallThickness, 
  +0.5*k4x4AlvWidth - kAlvWallThickness - kAlvHoleWidth,
  +0.5*k4x4AlvWidth - kAlvWallThickness - kAlvHoleWidth,
  +0.5*k4x4AlvWidth,  0.5*k4x4AlvWidth}

Definition at line 174 of file createRootGeoFileBwEnd_2017.C.

Referenced by createRootGeoFileBwEnd_2017().

const Double_t v5X[nVtx5]

static

Initial value:

= {
  -0.5*k4x4AlvWidth,  
  -0.5*k4x4AlvWidth + kSingleAlvWidth, 
  -0.5*k4x4AlvWidth + kSingleAlvWidth, 
  -0.5*k4x4AlvWidth}

Definition at line 183 of file createRootGeoFileBwEnd_2017.C.

const Double_t v5Y[nVtx5]

static

Initial value:

= {
  -0.5*k4x4AlvWidth, -0.5*k4x4AlvWidth,  
  -0.5*k4x4AlvWidth + kSingleAlvWidth, 
  -0.5*k4x4AlvWidth + kSingleAlvWidth}

Definition at line 188 of file createRootGeoFileBwEnd_2017.C.