#include "TGeoTrd1.h"
#include "TGeoCompositeShape.h"
#include "TGeoPara.h"
-
+#include "TGeoPhysicalNode.h"
+#include "TGeoHalfSpace.h"
+#include "TTreeStream.h"
+#include "TGeoArb8.h"
using std::ifstream;
using std::ios_base;
ClassImp(AliTPCv0)
TGeoMedium *sm3 = gGeoManager->GetMedium("TPC_Prepreg2");
TGeoMedium *sm4 = gGeoManager->GetMedium("TPC_Nomex");
//
+ //
TGeoVolume *tov1 = new TGeoVolume("TPC_OCV1",to1,sm1);
TGeoVolume *tov2 = new TGeoVolume("TPC_OCV2",to2,sm2);
TGeoVolume *tov3 = new TGeoVolume("TPC_OCV3",to3,sm3);
TGeoMedium *m7=gGeoManager->GetMedium("TPC_Cu");
TGeoMedium *m10 = gGeoManager->GetMedium("TPC_Alumina");
TGeoMedium *m11 = gGeoManager->GetMedium("TPC_Peek");
- // upper and lower rods differ in length!
- Double_t *upar;
- upar=NULL;
- gGeoManager->Volume("TPC_Rod","TUBE",m6->GetId(),upar);
- upar=new Double_t [3];
- upar[0]=1.8;
- upar[1]=2.2;
-
+ TGeoMedium *m13 = gGeoManager->GetMedium("TPC_Brass");
+ //
+ // tpc rod is an assembly of 10 long parts and 2 short parts
+ // connected with alu rings and plagged on both sides.
//
- //HV rods - makrolon + 0.58cm (diameter) Cu
- TGeoTube *hvr = new TGeoTube(0.,2.2,126.64);
- TGeoTube *hvc = new TGeoTube(0.,0.29,126.64);
//
- TGeoVolume *hvrv = new TGeoVolume("TPC_HV_Rod",hvr,m6);
- TGeoVolume *hvcv = new TGeoVolume("TPC_HV_Cable",hvc,m7);
- hvrv->AddNode(hvcv,1);
+// tpc rod long
+//
+ TGeoPcon *rod = new TGeoPcon("rod",0.,360.,6);
+ rod->DefineSection(0,-10.43,1.92,2.08);
+ rod->DefineSection(1,-9.75,1.92,2.08);
+
+ rod->DefineSection(2,-9.75,1.8,2.2);
+ rod->DefineSection(3,9.75,1.8,2.2);
+
+ rod->DefineSection(4,9.75,1.92,2.08);
+ rod->DefineSection(5,10.43,1.92,2.08);
+ //
+ TGeoVolume *mrodl = new TGeoVolume("TPC_mrodl",rod,m6);
+ //
+ // tpc rod short
+ //
+ TGeoPcon *rod1 = new TGeoPcon("rod1",0.,360.,6);
+ rod1->DefineSection(0,-8.93,1.92,2.08);
+ rod1->DefineSection(1,-8.25,1.92,2.08);
+
+ rod1->DefineSection(2,-8.25,1.8,2.2);
+ rod1->DefineSection(3,8.25,1.8,2.2);
+
+ rod1->DefineSection(4,8.25,1.92,2.08);
+ rod1->DefineSection(5,8.93,1.92,2.08);
+ //
+ TGeoVolume *mrods = new TGeoVolume("TPC_mrods",rod1,m6);
+ //
+ // below is for the resistor rod
+ //
+ // hole for the brass connectors
+ //
+
+ new TGeoTube("hhole",0.,0.3,0.3);
+ //
+ //transformations for holes - initialy they
+ // are placed at x=0 and negative y
+ //
+ TGeoRotation *rhole = new TGeoRotation();
+ rhole->RotateX(90.);
+ TGeoCombiTrans *transf[13];
+ Char_t name[30];
+ for(Int_t i=0;i<13;i++){
+ //sprintf(name,"transf%d",i);
+ snprintf(name,30,"transf%d",i);
+ transf[i]= new TGeoCombiTrans(name,0.,-2.,-9.+i*1.5,rhole);
+ transf[i]->RegisterYourself();
+ }
+ // union expression for holes
+ TString operl("hhole:transf0");
+ for (Int_t i=1;i<13;i++){
+ //sprintf(name,"+hhole:transf%d",i);
+ snprintf(name,30,"+hhole:transf%d",i);
+ operl.Append(name);
+ }
+ //
+ TString opers("hhole:transf1");
+ for (Int_t i=2;i<12;i++){
+ //sprintf(name,"+hhole:transf%d",i);
+ snprintf(name,30,"+hhole:transf%d",i);
+ opers.Append(name);
+ }
+ //union of holes
+ new TGeoCompositeShape("hlv",operl.Data());
+ new TGeoCompositeShape("hsv",opers.Data());
+ //
+ TGeoCompositeShape *rodl = new TGeoCompositeShape("rodl","rod-hlv");
+ TGeoCompositeShape *rods = new TGeoCompositeShape("rods","rod1-hsv");
+ //rods - volumes - makrolon rods with holes
+ TGeoVolume *rodlv = new TGeoVolume("TPC_rodl",rodl,m6);
+ TGeoVolume *rodsv = new TGeoVolume("TPC_rods",rods,m6);
+ //brass connectors
+ //connectors
+ TGeoTube *bcon = new TGeoTube(0.,0.3,0.3);//connectors
+ TGeoVolume *bconv = new TGeoVolume("TPC_bcon",bcon,m13);
+ //
+ // hooks holding strips
+ //
+ new TGeoBBox("hk1",0.625,0.015,0.75);
+ new TGeoBBox("hk2",0.625,0.015,0.15);
+ TGeoTranslation *tr21 = new TGeoTranslation("tr21",0.,-0.03,-0.6);
+ TGeoTranslation *tr12 = new TGeoTranslation("tr12",0.,-0.03,0.6);
+ tr21->RegisterYourself();
+ tr12->RegisterYourself();
+
+ TGeoCompositeShape *hook = new TGeoCompositeShape("hook","hk1+hk2:tr21+hk2:tr12");
+ TGeoVolume *hookv = new TGeoVolume("TPC_hook",hook,m13);
+ //
+ // assembly of the short rod with connectors and hooks
+ //
+ //
+ // short rod
+ //
+ TGeoVolumeAssembly *spart = new TGeoVolumeAssembly("TPC_spart");
+ //
+ spart->AddNode( rodsv,1);
+ for(Int_t i=1;i<12;i++){
+ spart->AddNode(bconv,i,transf[i]);
+ }
+ for(Int_t i =0;i<11;i++){
+ spart->AddNode(hookv,i+1,new TGeoTranslation(0.,-2.315,-7.5+i*1.5));
+ }
+ //
+ // long rod
+ //
+ TGeoVolumeAssembly *lpart = new TGeoVolumeAssembly("TPC_lpart");
+ //
+ lpart->AddNode( rodlv,1);
+ for(Int_t i=0;i<13;i++){
+ lpart->AddNode(bconv,i+12,transf[i]);
+ }
+ for(Int_t i =0;i<13;i++){
+ lpart->AddNode(hookv,i+12,new TGeoTranslation(0.,-2.315,-9.+i*1.5));
+ }
+ //
+ // alu ring
+ //
+ new TGeoTube("ring1",2.1075,2.235,0.53);
+ new TGeoTube("ring2",1.7925,1.89,0.43);
+ new TGeoTube("ring3",1.89,2.1075,0.05);
+ TGeoCompositeShape *ring = new TGeoCompositeShape("ring","ring1+ring2+ring3");
+ TGeoVolume *ringv = new TGeoVolume("TPC_ring",ring,m3);
+ //
+ // rod assembly
+ //
+ TGeoVolumeAssembly *tpcrrod = new TGeoVolumeAssembly("TPC_rrod");//rrod
+ TGeoVolumeAssembly *tpcmrod = new TGeoVolumeAssembly("TPC_mrod");//makrolon rod
+ //long pieces
+ for(Int_t i=0;i<11;i++){
+ tpcrrod->AddNode(ringv,i+1,new TGeoTranslation(0.,0.,-105.+i*21));
+ tpcmrod->AddNode(ringv,i+12,new TGeoTranslation(0.,0.,-105.+i*21));
+ }
+ for(Int_t i=0;i<10;i++){
+ tpcrrod->AddNode(lpart,i+1,new TGeoTranslation(0.,0.,-94.5+i*21));//resistor rod
+ tpcmrod->AddNode(mrodl,i+1,new TGeoTranslation(0.,0.,-94.5+i*21));//makrolon rod
+ }
//
- // resistor rods
+ // right plug - identical for all rods
//
- TGeoTube *cri = new TGeoTube(0.,0.45,126.64); //inner
- TGeoTube *cro = new TGeoTube(0.,0.45,126.54); //outer
- TGeoTube *cwi = new TGeoTube(0.,0.15,126.64); // water inner
- TGeoTube *cwo = new TGeoTube(0.,0.15,126.54); // water outer
+ TGeoPcon *tpcrp = new TGeoPcon(0.,360.,6);
//
- TGeoVolume *criv = new TGeoVolume("TPC_CR_I",cri,m10);
- TGeoVolume *crov = new TGeoVolume("TPC_CR_O",cro,m10);
- TGeoVolume *cwiv = new TGeoVolume("TPC_W_I",cwi,m11);
- TGeoVolume *cwov = new TGeoVolume("TPC_W_O",cwo,m11);
+ tpcrp->DefineSection(0,123.05,1.89,2.1075);
+ tpcrp->DefineSection(1,123.59,1.89,2.1075);
//
- // ceramic rod with water
+ tpcrp->DefineSection(2,123.59,1.8,2.2);
+ tpcrp->DefineSection(3,127.,1.8,2.2);
//
- criv->AddNode(cwiv,1);
- crov->AddNode(cwov,1);
+ tpcrp->DefineSection(4,127.,0.,2.2);
+ tpcrp->DefineSection(5,127.5,0.,2.2);
//
- TGeoTube *pri =new TGeoTube(0.2,0.35,126.64); //inner
- TGeoTube *pro = new TGeoTube(0.2,0.35,126.54); //outer
+ TGeoVolume *tpcrpv = new TGeoVolume("TPC_RP",tpcrp,m6);
//
- // peek rod
+ // adding short pieces and right plug
//
- TGeoVolume *priv = new TGeoVolume("TPC_PR_I",pri,m12);
- TGeoVolume *prov = new TGeoVolume("TPC_PR_O",pro,m12);
+ tpcrrod->AddNode(spart,1,new TGeoTranslation(0.,0.,-114.));
+ tpcrrod->AddNode(spart,2,new TGeoTranslation(0.,0.,114.));
+ tpcrrod->AddNode(ringv,23,new TGeoTranslation(0.,0.,-123.));
+ tpcrrod->AddNode(ringv,24,new TGeoTranslation(0.,0.,123.));
+ tpcrrod->AddNode(tpcrpv,1);
//
- // resistor rods assembly
+ tpcmrod->AddNode(mrods,1,new TGeoTranslation(0.,0.,-114.));
+ tpcmrod->AddNode(mrods,2,new TGeoTranslation(0.,0.,114.));
+ tpcmrod->AddNode(ringv,25,new TGeoTranslation(0.,0.,-123.));
+ tpcmrod->AddNode(ringv,26,new TGeoTranslation(0.,0.,123.));
+ tpcmrod->AddNode(tpcrpv,2);
//
- TGeoRotation *rotr = new TGeoRotation("rotr");
- rotr->RotateZ(-21.);
- TGeoTube *rri = new TGeoTube(1.8,2.2,126.64);//inner
- TGeoTube *rro = new TGeoTube(1.8,2.2,126.54);//inner
+ // from the ringv position to the CM is 3.0 cm!
+ //----------------------------------------
//
- TGeoVolume *rriv = new TGeoVolume("TPC_RR_I",rri,m6);
- TGeoVolume *rrov = new TGeoVolume("TPC_RR_O",rro,m6);
//
- TGeoVolumeAssembly *rrin = new TGeoVolumeAssembly("TPC_RROD_I");
- TGeoVolumeAssembly *rrou = new TGeoVolumeAssembly("TPC_RROD_O");
- rrin->AddNode(rriv,1);
- rrin->AddNode(criv,1,new TGeoTranslation(0.5,0.866, 0.));
- rrin->AddNode(criv,2,new TGeoTranslation(0.5,-0.866, 0.));
- rrin->AddNode(priv,1);
+ //HV rods - makrolon + 0.58cm (diameter) Cu ->check the length
+ TGeoTube *hvr = new TGeoTube(0.,1.465,123.);
+ TGeoTube *hvc = new TGeoTube(0.,0.29,123.);
//
- rrou->AddNode(rrov,1);
- rrou->AddNode(crov,1,new TGeoTranslation(0.5,0.866, 0.));
- rrou->AddNode(crov,2,new TGeoTranslation(0.5,-0.866, 0.));
- rrou->AddNode(prov,1);
+ TGeoVolume *hvrv = new TGeoVolume("TPC_HV_Rod",hvr,m6);
+ TGeoVolume *hvcv = new TGeoVolume("TPC_HV_Cable",hvc,m7);
+ hvrv->AddNode(hvcv,1);
+ //
+ //resistor rod
+ //
+ TGeoTube *cr = new TGeoTube(0.,0.45,123.);
+ TGeoTube *cw = new TGeoTube(0.,0.15,123.);
+ TGeoVolume *crv = new TGeoVolume("TPC_CR",cr,m10);
+ TGeoVolume *cwv = new TGeoVolume("TPC_W",cw,m12);
+ //
+ // ceramic rod with water
+ //
+ crv->AddNode(cwv,1);
+ //
+ //peek rod
+ //
+ TGeoTube *pr =new TGeoTube(0.2,0.35,123.);
+ TGeoVolume *prv = new TGeoVolume("TPC_PR",pr,m11);
+ //
+ // copper plates with connectors
+ //
+ new TGeoTube("tub",0.,1.7,0.025);
+ //
+ // half space - points on the plane and a normal vector
+ //
+ Double_t n[3],p[3];
+ Double_t slope = TMath::Tan(22.*TMath::DegToRad());
+ Double_t intp = 1.245;
+ //
+ Double_t b = slope*slope+1.;
+ p[0]=intp*slope/b;
+ p[1]=-intp/b;
+ p[2]=0.;
+ //
+ n[0]=-p[0];
+ n[1]=-p[1];
+ n[2]=0.;
+ Double_t norm;
+ norm=TMath::Sqrt(n[0]*n[0]+n[1]*n[1]);
+ n[0] /= norm;
+ n[1] /=norm;
+ //
+ new TGeoHalfSpace("sp1",p,n);
+ //
+ slope = -slope;
+ //
+ p[0]=intp*slope/b;
+ p[1]=-intp/b;
+ //
+ n[0]=-p[0];
+ n[1]=-p[1];
+ norm=TMath::Sqrt(n[0]*n[0]+n[1]*n[1]);
+ n[0] /= norm;
+ n[1] /=norm;
+ //
+ new TGeoHalfSpace("sp2",p,n);
+ // holes for rods
+ //holes
+ new TGeoTube("h1",0.,0.5,0.025);
+ new TGeoTube("h2",0.,0.35,0.025);
+ //translations:
+ TGeoTranslation *ttr11 = new TGeoTranslation("ttr11",-0.866,0.5,0.);
+ TGeoTranslation *ttr22 = new TGeoTranslation("ttr22",0.866,0.5,0.);
+ ttr11->RegisterYourself();
+ ttr22->RegisterYourself();
+ // elastic connector
+ new TGeoBBox("elcon",0.72,0.005,0.3);
+ TGeoRotation *crr1 = new TGeoRotation();
+ crr1->RotateZ(-22.);
+TGeoCombiTrans *ctr1 = new TGeoCombiTrans("ctr1",-0.36011, -1.09951,-0.325,crr1);
+ctr1->RegisterYourself();
+ TGeoCompositeShape *cs1 = new TGeoCompositeShape("cs1",
+"(((((tub-h1:ttr11)-h1:ttr22)-sp1)-sp2)-h2)+elcon:ctr1");
+ //
+ TGeoVolume *csvv = new TGeoVolume("TPC_RR_CU",cs1,m7);
+ //
+ // resistor rod assembly 2 ceramic rods, peak rod, Cu plates
+ // and resistors
+ //
+ TGeoVolumeAssembly *rrod = new TGeoVolumeAssembly("TPC_RRIN");
+ // rods
+ rrod->AddNode(crv,1,ttr11);
+ rrod->AddNode(crv,2,ttr22);
+ rrod->AddNode(prv,1);
+ //Cu plates
+ for(Int_t i=0;i<165;i++){
+ rrod->AddNode(csvv,i+1,new TGeoTranslation(0.,0.,-122.675+i*1.5));
+ }
+ //resistors
+ TGeoTube *res = new TGeoTube(0.,0.15,0.5);
+ TGeoVolume *resv = new TGeoVolume("TPC_RES",res,m10);
+ TGeoVolumeAssembly *ress = new TGeoVolumeAssembly("TPC_RES_CH");
+ ress->AddNode(resv,1,new TGeoTranslation(0.2,0.,0.));
+ ress->AddNode(resv,2,new TGeoTranslation(-0.2,0.,0.));
+ //
+ TGeoRotation *crr2 = new TGeoRotation();
+ crr2->RotateY(30.);
+ TGeoRotation *crr3 = new TGeoRotation();
+ crr3->RotateY(-30.);
+ //
+ for(Int_t i=0;i<164;i+=2){
+ rrod->AddNode(ress,i+1, new TGeoCombiTrans(0.,1.2,-121.925+i*1.5,crr2));
+ rrod->AddNode(ress,i+2, new TGeoCombiTrans(0.,1.2,-121.925+(i+1)*1.5,crr3));
+ }
+
+ tpcrrod->AddNode(rrod,1,new TGeoCombiTrans(0.,0.,0.5,crr1));
+ //
+ // rod left head with holders - inner
+ //
+ // first element - support for inner holder TPC_IHS
+ Double_t shift1[3] = {0.0,-0.175,0.0};
+
+new TGeoBBox("tpcihs1", 4.7, 0.66, 2.35);
+new TGeoBBox("tpcihs2", 4.7, 0.485, 1.0, shift1);
+new TGeoBBox("tpcihs3", 1.5, 0.485, 2.35, shift1);
+new TGeoTube("tpcihs4", 0.0, 2.38, 0.1);
+//
+Double_t pointstrap[16];
+pointstrap[0]= 0.0;
+pointstrap[1]= 0.0;
+pointstrap[2]= 0.0;
+pointstrap[3]= 1.08;
+pointstrap[4]= 2.3;
+pointstrap[5]= 1.08;
+pointstrap[6]= 3.38;
+pointstrap[7]= 0.0;
+pointstrap[8]= 0.0;
+pointstrap[9]= 0.0;
+pointstrap[10]= 0.0;
+pointstrap[11]= 1.08;
+pointstrap[12]= 2.3;
+pointstrap[13]= 1.08;
+pointstrap[14]= 3.38;
+pointstrap[15]= 0.0;
+//
+TGeoArb8 *tpcihs5 = new TGeoArb8("tpcihs5", 0.6, pointstrap);
+//
+// half space - cutting "legs"
+//
+ p[0]=0.0;
+ p[1]=0.105;
+ p[2]=0.0;
+ //
+n[0] = 0.0;
+n[1] = 1.0;
+n[2] = 0.0;
+
+new TGeoHalfSpace("cutil1", p, n);
+
+//
+// transformations
+//
+TGeoTranslation *trans2 = new TGeoTranslation("trans2", 0.0, 2.84, 2.25);
+trans2->RegisterYourself();
+TGeoTranslation*trans3= new TGeoTranslation("trans3", 0.0, 2.84, -2.25);
+trans3->RegisterYourself();
+//support - composite volume
+//
+TGeoCompositeShape *tpcihs6 = new TGeoCompositeShape("tpcihs6", "tpcihs1-(tpcihs2+tpcihs3)-(tpcihs4:trans2)-(tpcihs4:trans3)-cutil1");
+//
+// volumes - all makrolon
+//
+ TGeoVolume *tpcihss = new TGeoVolume("TPC_IHSS", tpcihs6, m6); //support
+ TGeoVolume *tpcihst = new TGeoVolume("TPC_IHSTR",tpcihs5 , m6); //trapesoid
+ //now assembly
+TGeoRotation *rot111 = new TGeoRotation();
+rot111->RotateY(180.0);
+//
+TGeoVolumeAssembly *tpcihs = new TGeoVolumeAssembly("TPC_IHS"); // assembly of the support
+tpcihs->AddNode(tpcihss, 1);
+tpcihs->AddNode(tpcihst, 1, new TGeoTranslation(-4.7, 0.66, 0.0));
+tpcihs->AddNode(tpcihst, 2, new TGeoCombiTrans(4.7, 0.66, 0.0, rot111));
+//
+// two rod holders (TPC_IRH) assembled with the support
+//
+new TGeoBBox("tpcirh1", 4.7, 1.33, 0.5);
+ shift1[0]=-3.65;
+ shift1[1]=0.53;
+ shift1[2]=0.;
+new TGeoBBox("tpcirh2", 1.05, 0.8, 0.5, shift1);
+ shift1[0]=3.65;
+ shift1[1]=0.53;
+ shift1[2]=0.;
+new TGeoBBox("tpcirh3", 1.05, 0.8, 0.5, shift1);
+ shift1[0]=0.0;
+ shift1[1]=1.08;
+ shift1[2]=0.;
+new TGeoBBox("tpcirh4", 1.9, 0.25, 0.5, shift1);
+new TGeoTube("tpcirh5", 0, 1.9, 5);
+//
+TGeoTranslation *trans4 = new TGeoTranslation("trans4", 0, 0.83, 0.0);
+trans4->RegisterYourself();
+//
+TGeoCompositeShape *tpcirh6 = new TGeoCompositeShape("tpcirh6", "tpcirh1-tpcirh2-tpcirh3-(tpcirh5:trans4)-tpcirh4");
+//
+// now volume
+//
+TGeoVolume *tpcirh = new TGeoVolume("TPC_IRH", tpcirh6, m6);
+//
+// and all together...
+//
+TGeoVolume *tpciclamp = new TGeoVolumeAssembly("TPC_ICLP");
+tpciclamp->AddNode(tpcihs, 1);
+tpciclamp->AddNode(tpcirh, 1, new TGeoTranslation(0, 1.99, 1.1));
+tpciclamp->AddNode(tpcirh, 2, new TGeoTranslation(0, 1.99, -1.1));
+//
+// and now left inner "head"
+//
+TGeoPcon *inplug = new TGeoPcon("inplug", 0.0, 360.0, 14);
+
+inplug->DefineSection(0, 0.3, 0.0, 2.2);
+inplug->DefineSection(1, 0.6, 0.0, 2.2);
+
+inplug->DefineSection(2, 0.6, 0.0, 1.75);
+inplug->DefineSection(3, 0.7, 0.0, 1.75);
+
+inplug->DefineSection(4, 0.7, 1.55, 1.75);
+inplug->DefineSection(5, 1.6, 1.55, 1.75);
+
+inplug->DefineSection(6, 1.6, 1.55, 2.2);
+inplug->DefineSection(7, 1.875, 1.55, 2.2);
+
+inplug->DefineSection(8, 1.875, 1.55, 2.2);
+inplug->DefineSection(9, 2.47, 1.75, 2.2);
+
+inplug->DefineSection(10, 2.47, 1.75, 2.08);
+inplug->DefineSection(11, 2.57, 1.8, 2.08);
+
+inplug->DefineSection(12, 2.57, 1.92, 2.08);
+inplug->DefineSection(13, 2.95, 1.92, 2.08);
+//
+ shift1[0]=0.0;
+ shift1[1]=-2.09;
+ shift1[2]=1.075;
+ //
+new TGeoBBox("pcuti", 1.5, 0.11, 1.075, shift1);
+//
+TGeoCompositeShape *inplleft = new TGeoCompositeShape("inplleft", "inplug-pcuti");
+TGeoVolume *tpcinlplug = new TGeoVolume("TPC_INPLL", inplleft, m6);
+//
+// holder + plugs
+//
+ TGeoVolume *tpcihpl = new TGeoVolumeAssembly("TPC_IHPL"); //holder+2 plugs (reflected)
+ tpcihpl->AddNode(tpcinlplug, 1);
+ tpcihpl->AddNode(tpcinlplug, 2,ref);
+ tpcihpl->AddNode(tpciclamp,1,new TGeoTranslation(0.0, -2.765, 0.0));
+ //
+ // outer holders and clamps
+ //
+
+ // outer membrane holder (between rods)
+pointstrap[0]= 0.0;
+pointstrap[1]= 0.0;
+pointstrap[2]= 0.0;
+pointstrap[3]= 2.8;
+pointstrap[4]= 3.1;
+pointstrap[5]= 2.8-3.1*TMath::Tan(15.*TMath::DegToRad());
+pointstrap[6]= 3.1;
+pointstrap[7]= 0.0;
+pointstrap[8]= 0.0;
+pointstrap[9]= 0.0;
+pointstrap[10]= 0.0;
+pointstrap[11]= 2.8;
+pointstrap[12]= 3.1;
+pointstrap[13]= 2.8-3.1*TMath::Tan(15.*TMath::DegToRad());
+pointstrap[14]= 3.1;
+pointstrap[15]= 0.0;
+//
+TGeoArb8 *tpcomh1 = new TGeoArb8("tpcomh1", 1.05, pointstrap);
+TGeoBBox *tpcomh2 = new TGeoBBox("tpcomh2", 0.8, 1.4, 6);
+//
+TGeoVolume *tpcomh1v = new TGeoVolume("TPC_OMH1", tpcomh1, m7);
+TGeoVolume *tpcomh2v = new TGeoVolume("TPC_OMH2", tpcomh2, m7);
+//
+TGeoVolume *tpcomh3v = new TGeoVolumeAssembly("TPC_OMH3"); // assembly1
+tpcomh3v->AddNode(tpcomh1v, 1, new TGeoTranslation(0.8, -1.4, 4.95));
+ tpcomh3v->AddNode(tpcomh1v, 2, new TGeoTranslation(0.8, -1.4, -4.95));
+tpcomh3v->AddNode(tpcomh2v, 1);
+//
+shift1[0] = 0.9;
+shift1[1] = -1.85;
+shift1[2] = 0.0;
+//
+new TGeoBBox("tpcomh3", 1.65, 1.15, 3.4);
+TGeoBBox *tpcomh4 = new TGeoBBox("tpcomh4", 0.75, 0.7, 3.4, shift1);
+//
+// halfspace 1
+//
+p[0] = 0.0;
+p[1] = -1.05;
+p[2] = -3.4;
+//
+n[0] = 0.0;
+n[1] = -1.0*TMath::Tan(30.*TMath::DegToRad());
+n[2] = 1.0;
+//
+new TGeoHalfSpace("cutomh1", p, n);
+//
+// halfspace 2
+//
+p[0] = 0.0;
+p[1] = -1.05;
+p[2] = 3.4;
+//
+n[0] = 0.0;
+n[1] = -1.0*TMath::Tan(30.*TMath::DegToRad());
+n[2] = -1.0;
+//
+new TGeoHalfSpace("cutomh2", p, n);
+//
+// halfspace 3
+//
+p[0] = -1.65;
+p[1] = 0.0;
+p[2] = -0.9;
+//
+n[0] = 1.0*TMath::Tan(75.*TMath::DegToRad());
+n[1] = 0.0;
+n[2] = 1.0;
+//
+new TGeoHalfSpace("cutomh3", p, n);
+//
+// halfspace 4
+//
+p[0] = -1.65;
+p[1] = 0.0;
+p[2] = 0.9;
+//
+n[0] = 1.0*TMath::Tan(75*TMath::DegToRad());
+n[1] = 0.0;
+n[2] = -1.0;
+//
+new TGeoHalfSpace("cutomh4", p, n);
+//
+// halsfspace 5
+//
+p[0] = 1.65;
+p[1] = -1.05;
+p[2] = 0.0;
+//
+n[0] = -1.0;
+n[1] = -1.0*TMath::Tan(20.*TMath::DegToRad());
+n[2] = 0.0;
+//
+new TGeoHalfSpace("cutomh5", p, n);
+//
+TGeoCompositeShape *tpcomh5 = new TGeoCompositeShape("tpcomh5", "tpcomh3-cutomh1-cutomh2-cutomh3-cutomh4-cutomh5");
+//
+ TGeoVolume *tpcomh5v = new TGeoVolume("TPC_OMH5",tpcomh5,m6);
+ TGeoVolume *tpcomh4v = new TGeoVolume("TPC_OMH6",tpcomh4,m6);
+ //
+ TGeoVolumeAssembly *tpcomh7v = new TGeoVolumeAssembly("TPC_OMH7");
+ tpcomh7v->AddNode(tpcomh5v,1);
+ tpcomh7v->AddNode(tpcomh4v,1);
+ //
+ // full membrane holder - tpcomh3v + tpcomh7v
+ //
+ TGeoVolumeAssembly *tpcomh = new TGeoVolumeAssembly("TPC_OMH");
+ tpcomh->AddNode(tpcomh3v,1,new TGeoTranslation(1.5,0.,0.));
+ tpcomh->AddNode(tpcomh3v,2,new TGeoCombiTrans(-1.5,0.,0.,rot111));
+ tpcomh->AddNode(tpcomh7v,1,new TGeoTranslation(0.65+1.5, 2.55, 0.0));
+ tpcomh->AddNode(tpcomh7v,2,new TGeoCombiTrans(-0.65-1.5, 2.55, 0.0,rot111));
+ //
+ // outer rod holder support
+ //
+new TGeoBBox("tpcohs1", 3.8, 0.675, 2.35);
+//
+shift1[0] = 0.0;
+shift1[1] = 0.175;
+shift1[2] = 0.0;
+//
+new TGeoBBox("tpcohs2", 1.5, 0.5, 2.35, shift1);
+new TGeoBBox("tpcohs3", 3.8, 0.5, 0.85, shift1);
+//
+shift1[0] = 0.0;
+shift1[1] = -1.175;
+shift1[2] = 0.0;
+//
+TGeoBBox *tpcohs4 = new TGeoBBox("tpsohs4", 3.1, 0.5, 0.7, shift1);
+//
+TGeoVolume *tpcohs4v = new TGeoVolume("TPC_OHS4", tpcohs4, m6);
+//
+p[0] = 0.0;
+p[1] = -0.186;
+p[2] = 0.0;
+//
+n[0] = 0.0;
+n[1] = -1.0;
+n[2] = 0.0;
+//
+new TGeoHalfSpace("cutohs1", p, n);
+//
+TGeoCompositeShape *tpcohs5 = new TGeoCompositeShape("tpcohs5", "tpcohs1-tpcohs2-tpcohs3-cutohs1");
+TGeoVolume *tpcohs5v = new TGeoVolume("TPC_OHS5", tpcohs5, m6);
+//
+TGeoVolumeAssembly *tpcohs = new TGeoVolumeAssembly("TPC_OHS");
+tpcohs->AddNode(tpcohs5v, 1);
+tpcohs->AddNode(tpcohs4v, 1);
+//
+// outer rod holder itself
+//
+shift1[0] = 0.0;
+shift1[1] = 1.325;
+shift1[2] = 0.0;
+ new TGeoBBox("tpcorh1", 3.1, 1.825, 0.55); //from this box we cut pieces...
+ //
+shift1[0] = -3.1;
+shift1[1] = -0.5;
+shift1[2] = 0.0;
+//
+new TGeoBBox("tpcorh2", 0.5, 2.75, 1.1, shift1);
+//
+shift1[0] = 3.1;
+shift1[1] = -0.5;
+shift1[2] = 0.0;
+//
+new TGeoBBox("tpcorh3", 0.5, 2.75, 1.1, shift1);
+//
+shift1[0] = 0.0;
+shift1[1] = -0.5;
+shift1[2] = -0.95;
+//
+new TGeoBBox("tpcorh4", 3.9, 2.75, 0.5, shift1);
+//
+shift1[0] = 0.0;
+shift1[1] = -0.5;
+shift1[2] = 0.0;
+//
+new TGeoBBox("tpcorh5", 1.95, 0.5, 1.1, shift1);
+//
+shift1[0] = 0.0;
+shift1[1] = -0.5;
+shift1[2] = 0.55;
+//
+new TGeoBBox("tpcorh6", 2.4, 0.5, 0.6, shift1);
+//
+new TGeoTube("tpcorh7", 0, 1.95, 0.85);
+new TGeoTube("tpcorh8", 0, 2.4, 0.6);
+//
+TGeoTranslation *trans33 = new TGeoTranslation("trans33", 0.0, 0.0, 0.55);
+trans33->RegisterYourself();
+//
+TGeoCompositeShape *tpcorh9 = new TGeoCompositeShape("tpcorh9", "tpcorh1-tpcorh2-tpcorh3-tpcorh4-tpcorh5-tpcorh6-(tpcorh8:trans33)-tpcorh7");
+//
+ TGeoVolume *tpcorh9v = new TGeoVolume("TPC_ORH",tpcorh9,m6); //outer rod holder
+ //
+ // now 2 holders together
+ //
+ TGeoVolumeAssembly *tpcorh = new TGeoVolumeAssembly("TPC_ORH2");
+ //
+ tpcorh->AddNode(tpcorh9v,1,new TGeoTranslation(0.0, 0.0, 1.25));
+ tpcorh->AddNode(tpcorh9v,2,new TGeoCombiTrans(0.0, 0.0, -1.25,rot111));
+//
+// outer rod plug left
+//
+TGeoPcon *outplug = new TGeoPcon("outplug", 0.0, 360.0, 14);
+
+outplug->DefineSection(0, 0.5, 0.0, 2.2);
+outplug->DefineSection(1, 0.7, 0.0, 2.2);
+
+outplug->DefineSection(2, 0.7, 1.55, 2.2);
+outplug->DefineSection(3, 0.8, 1.55, 2.2);
+
+outplug->DefineSection(4, 0.8, 1.55, 1.75);
+outplug->DefineSection(5, 1.2, 1.55, 1.75);
+
+outplug->DefineSection(6, 1.2, 1.55, 2.2);
+outplug->DefineSection(7, 1.875, 1.55, 2.2);
+
+outplug->DefineSection(8, 1.875, 1.55, 2.2);
+outplug->DefineSection(9, 2.47, 1.75, 2.2);
+
+outplug->DefineSection(10, 2.47, 1.75, 2.08);
+outplug->DefineSection(11, 2.57, 1.8, 2.08);
+
+outplug->DefineSection(12, 2.57, 1.92, 2.08);
+outplug->DefineSection(13, 2.95, 1.92, 2.08);
+//
+shift1[0] = 0.0;
+shift1[1] = 2.09;
+shift1[2] = 1.01;
+
+new TGeoBBox("cutout", 2.5, 0.11, 1.01, shift1);
+//
+
+TGeoCompositeShape *outplleft = new TGeoCompositeShape("outplleft", "outplug-cutout");
+TGeoVolume *outplleftv = new TGeoVolume("TPC_OPLL", outplleft, m6);
+//
+// support + holder + plug
+//
+
+
+ TGeoVolumeAssembly *tpcohpl = new TGeoVolumeAssembly("TPC_OHPL");
+ //
+ tpcohpl->AddNode(outplleftv,1); //plug
+ tpcohpl->AddNode(outplleftv,2,ref); //plug reflected
+ tpcohpl->AddNode(tpcorh,1); //rod holder
+ tpcohpl->AddNode(tpcohs,1,new TGeoTranslation(0.0, 3.925, 0)); // support
+ //
+
+ //
+ // main membrane holder
+ //
+pointstrap[0]= 0.0;
+pointstrap[1]= 0.0;
+pointstrap[2]= 0.0;
+pointstrap[3]= 2.8;
+pointstrap[4]= 3.1;
+pointstrap[5]= 1.96;
+pointstrap[6]= 3.1;
+pointstrap[7]= 0.0;
+pointstrap[8]= 0.0;
+pointstrap[9]= 0.0;
+pointstrap[10]= 0.0;
+pointstrap[11]= 2.8;
+pointstrap[12]= 3.1;
+pointstrap[13]= 1.96;
+pointstrap[14]= 3.1;
+pointstrap[15]= 0.0;
+//
+TGeoArb8 *tpcmmh1 = new TGeoArb8("tpcmmh1", 1.75, pointstrap);
+TGeoBBox *tpcmmh2 = new TGeoBBox("tpcmmh2", 0.8, 1.4, 12.5);
+//
+TGeoVolume *tpcmmh1v = new TGeoVolume("TPC_MMH1", tpcmmh1, m6);
+TGeoVolume *tpcmmh2v = new TGeoVolume("TPC_MMH2", tpcmmh2, m6);
+//
+ TGeoVolumeAssembly *tpcmmhs = new TGeoVolumeAssembly("TPC_MMHS");
+ tpcmmhs->AddNode(tpcmmh1v,1,new TGeoTranslation(0.8, -1.4, 10.75));
+ tpcmmhs->AddNode(tpcmmh1v,2,new TGeoTranslation(0.8, -1.4, -10.75));
+ tpcmmhs->AddNode(tpcmmh2v,1);
+ //
+ // main membrahe holder clamp
+ //
+shift1[0] = -0.75;
+shift1[1] = -1.15;
+shift1[2] = 0.0;
+//
+new TGeoBBox("tpcmmhc1", 1.65, 1.85, 8.9);
+new TGeoBBox("tpcmmhc2", 0.9, 0.7, 8.9, shift1);
+//
+// half spaces - cuts
+//
+p[0] = -1.65;
+p[1] = 0.0;
+p[2] = -0.9;
+//
+n[0] = 8.0;
+n[1] = 0.0;
+n[2] = 8.0*TMath::Tan(13.*TMath::DegToRad());
+//
+new TGeoHalfSpace("cutmmh1", p, n);
+//
+p[0] = -1.65;
+p[1] = 0.0;
+p[2] = 0.9;
+//
+n[0] = 8.0;
+n[1] = 0.0;
+n[2] = -8.0*TMath::Tan(13.*TMath::DegToRad());
+//
+new TGeoHalfSpace("cutmmh2", p, n);
+//
+p[0] = 0.0;
+p[1] = 1.85;
+p[2] = -2.8;
+//
+n[0] = 0.0;
+n[1] = -6.1;
+n[2] = 6.1*TMath::Tan(20.*TMath::DegToRad());
+//
+new TGeoHalfSpace("cutmmh3", p, n);
+//
+p[0] = 0.0;
+p[1] = 1.85;
+p[2] = 2.8;
+//
+n[0] = 0.0;
+n[1] = -6.1;
+n[2] = -6.1*TMath::Tan(20*TMath::DegToRad());
+//
+new TGeoHalfSpace("cutmmh4", p, n);
+//
+p[0] = 0.75;
+p[1] = 0.0;
+p[2] = -8.9;
+//
+n[0] = 2.4*TMath::Tan(30*TMath::DegToRad());
+n[1] = 0.0;
+n[2] = 2.4;
+//
+ new TGeoHalfSpace("cutmmh5", p, n);
+//
+p[0] = 0.75;
+p[1] = 0.0;
+p[2] = 8.9;
+//
+n[0] = 2.4*TMath::Tan(30*TMath::DegToRad());
+n[1] = 0.0;
+n[2] = -2.4;
+//
+new TGeoHalfSpace("cutmmh6", p, n);
+
+TGeoCompositeShape *tpcmmhc = new TGeoCompositeShape("TPC_MMHC", "tpcmmhc1-tpcmmhc2-cutmmh1-cutmmh2-cutmmh3-cutmmh4-cutmmh5-cutmmh6");
+
+ TGeoVolume *tpcmmhcv = new TGeoVolume("TPC_MMHC",tpcmmhc,m6);
+ //
+TGeoVolume *tpcmmh = new TGeoVolumeAssembly("TPC_MMH");
+//
+ tpcmmh->AddNode(tpcmmhcv,1,new TGeoTranslation(0.65+1.5, 1.85, 0.0));
+ tpcmmh->AddNode(tpcmmhcv,2,new TGeoCombiTrans(-0.65-1.5, 1.85, 0.0,rot111));
+ tpcmmh->AddNode(tpcmmhs,1,new TGeoTranslation(1.5, 0.0, 0.0));
+ tpcmmh->AddNode(tpcmmhs,2,new TGeoCombiTrans(-1.5, 0.0, 0.0,rot111));
+ //
+
+ //
+
+//--------------------------------------------
+ //
+ // guard ring resistor chain
+ //
+
+ TGeoTube *gres1 = new TGeoTube(0.,0.375,125.);// inside ifc
+ //
+ TGeoVolume *vgres1 = new TGeoVolume("TPC_GRES1",gres1,m10);
+
+ //
+ Double_t xrc,yrc;
+ //
+ xrc=79.3*TMath::Cos(350.*TMath::DegToRad());
+ yrc=79.3*TMath::Sin(350.*TMath::DegToRad());
+ //
+ v9->AddNode(vgres1,1,new TGeoTranslation(xrc,yrc,126.9));
+ v9->AddNode(vgres1,2,new TGeoTranslation(xrc,yrc,-126.9));
+ //
+ xrc=79.3*TMath::Cos(190.*TMath::DegToRad());
+ yrc=79.3*TMath::Sin(190.*TMath::DegToRad());
+ //
+ v9->AddNode(vgres1,3,new TGeoTranslation(xrc,yrc,126.9));
+ v9->AddNode(vgres1,4,new TGeoTranslation(xrc,yrc,-126.9));
+ //------------------------------------------------------------------
+ TGeoRotation refl("refl",90.,0.,90.,90.,180.,0.);
+ TGeoRotation rotrod("rotrod");
+ //
+ TGeoRotation *rotpos[2];
+ //
+ TGeoRotation *rotrod1[2];
+
+
+ //v9 - drift gas
+
+ TGeoRotation rot102("rot102");
+ rot102.RotateY(-90.);
+
for(Int_t i=0;i<18;i++){
Double_t angle,x,y;
Double_t z,r;
angle=TMath::DegToRad()*20.*(Double_t)i;
+ //inner rods
r=81.5;
x=r * TMath::Cos(angle);
y=r * TMath::Sin(angle);
- upar[2]=126.64; //lower
- z= 126.96;
- if(i==3){
- v9->AddNode(rrin,1,new TGeoCombiTrans(x,y,z,rotr)); //A
- v9->AddNode(rrin,2,new TGeoCombiTrans(x,y,-z,rotr)); //C
+ z = 126.;
+ TGeoRotation *rot12 = new TGeoRotation();
+ rot12->RotateZ(-90.0+i*20.);
+ v9->AddNode(tpcihpl,i+1,new TGeoCombiTrans(x, y, 0., rot12));
+ //
+ if(i==11){//resistor rod inner
+ rotrod.RotateZ(-90.+angle);
+ rotrod1[0]= new TGeoRotation();
+ rotpos[0]= new TGeoRotation();
+ //
+ rotrod1[0]->RotateZ(-90.+angle);
+ *rotpos[0] = refl*rotrod; //rotation+reflection
+ v9->AddNode(tpcrrod,1,new TGeoCombiTrans(x,y, z, rotrod1[0])); //A
+ v9->AddNode(tpcrrod,2,new TGeoCombiTrans(x,y,-z, rotpos[0])); //C
}
else {
- gGeoManager->Node("TPC_Rod",i+1,"TPC_Drift",x,y,z,0,kTRUE,upar,3);//shaft
- gGeoManager->Node("TPC_Rod",i+19,"TPC_Drift",x,y,-z,0,kTRUE,upar,3);//muon
+ v9->AddNode(tpcmrod,i+1,new TGeoTranslation(x,y,z));//shaft
+ v9->AddNode(tpcmrod,i+19,new TGeoCombiTrans(x,y,-z,ref));//muon
}
-
-
- //
+ // outer rods
r=254.25;
x=r * TMath::Cos(angle);
y=r * TMath::Sin(angle);
- upar[2]=126.54; //upper
- z=127.06;
+ z=126.;
//
- if(i==15){
- //v9->AddNode(hvrv,1,new TGeoTranslation(x,y,z));//A-side only
- v9->AddNode(hvrv,1,new TGeoTranslation(x,y,127.14));//A-side only
- gGeoManager->Node("TPC_Rod",i+55,"TPC_Drift",x,y,-z,0,kTRUE,upar,3);
- }
- else if(i==11){
- v9->AddNode(rrou,1,new TGeoCombiTrans(x,y,z,rotr)); //A
- v9->AddNode(rrou,2,new TGeoCombiTrans(x,y,-z,rotr)); //C
+ // outer rod holder + outer left plug
+ //
+ TGeoRotation *rot33 = new TGeoRotation();
+ rot33->RotateZ(-90+i*20.);
+ //
+ v9->AddNode(tpcohpl,i+1,new TGeoCombiTrans(x, y, 0., rot33));
+ //
+ Double_t xxx = 256.297*TMath::Cos((i*20.+10.)*TMath::DegToRad());
+ Double_t yyy = 256.297*TMath::Sin((i*20.+10.)*TMath::DegToRad());
+ //
+ TGeoRotation rot101("rot101");
+ rot101.RotateZ(90.+i*20.+10.);
+ TGeoRotation *rot103 = new TGeoRotation("rot103");
+ *rot103 = rot101*rot102;
+ //
+ TGeoCombiTrans *trh100 = new TGeoCombiTrans(xxx,yyy,0.,rot103);
+ //
+ if(i==2) {
+ //main membrane holder
+ v9->AddNode(tpcmmh,1,trh100);
}
else{
+ // "normal" membrane holder
+ v9->AddNode(tpcomh,i+1,trh100);
+ }
+
//
- gGeoManager->Node("TPC_Rod",i+37,"TPC_Drift",x,y,z,0,kTRUE,upar,3);
- gGeoManager->Node("TPC_Rod",i+55,"TPC_Drift",x,y,-z,0,kTRUE,upar,3);
+ if(i==3){//resistor rod outer
+ rotrod.RotateZ(90.+angle);
+ rotrod1[1]= new TGeoRotation();
+ rotpos[1]= new TGeoRotation();
+ rotrod1[1]->RotateZ(90.+angle);
+ *rotpos[1] = refl*rotrod;//rotation+reflection
+ v9->AddNode(tpcrrod,3,new TGeoCombiTrans(x,y, z, rotrod1[1])); //A
+ v9->AddNode(tpcrrod,4,new TGeoCombiTrans(x,y, -z, rotpos[1])); //C
}
+ else {
+ v9->AddNode(tpcmrod,i+37,new TGeoTranslation(x,y,z));//shaft
+ v9->AddNode(tpcmrod,i+55,new TGeoCombiTrans(x,y,-z,ref));//muon
+ }
+ if(i==15){
+ v9->AddNode(hvrv,1,new TGeoTranslation(x,y,z+0.7)); //hv->A-side only
- }
+
+ }
+ } //end of rods positioning
TGeoVolume *alice = gGeoManager->GetVolume("ALIC");
alice->AddNode(v1,1);
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff