- B2->SetVertex(0,xbo[1],ybo[1]);
- B2->SetVertex(1,xbo[2],ybo[2]);
- B2->SetVertex(2,xbi[2],ybi[2]);
- B2->SetVertex(3,xbi[1],ybi[1]);
- //
- B3->SetVertex(0,xbo[5],ybo[5]);
- B3->SetVertex(1,xbo[6],ybo[6]);
- B3->SetVertex(2,xbi[6],ybi[6]);
- B3->SetVertex(3,xbi[5],ybi[5]);
- //--------------------------
- C1->SetVertex(0,xco[0],yco[0]);
- C1->SetVertex(1,xco[1],yco[1]);
- C1->SetVertex(2,xci[1],yci[1]);
- C1->SetVertex(3,xci[0],yci[0]);
- //
- C2->SetVertex(0,xco[1],yco[1]);
- C2->SetVertex(1,xco[2],yco[2]);
- C2->SetVertex(2,xci[2],yci[2]);
- C2->SetVertex(3,xci[1],yci[1]);
- //
- C3->SetVertex(0,xco[5],yco[5]);
- C3->SetVertex(1,xco[6],yco[6]);
- C3->SetVertex(2,xci[6],yci[6]);
- C3->SetVertex(3,xci[5],yci[5]);
- // Defining the hole, filled with air
- Double_t p1,c1,x,y,x7[3],y7[3];
- p1 = (xo[0]-xi[0])/(yo[0]-yi[0]);
- c1 = xo[0]+0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(xo[0]-xi[0])+
- SQ(yo[0]-yi[0]))/(xo[0]-xi[0]);
- y = TSCRoutA-2.*TSCarbonFiberThA;
- x = p1*(y-yo[0])+c1;
- Ah1->SetVertex(0,x,y);
- Bh1->SetVertex(0,x,y);
- Ch1->SetVertex(4,x,y);
- y = TSCRinA+TSCarbonFiberThA;
- x = p1*(y-yo[0])+c1;
- Ah1->SetVertex(3,x,y);
- Bh1->SetVertex(3,x,y);
- x7[0] = x; y7[0] = y; // vortexing done after last point
- //Ch1->SetVertex(7,x,y);
- p1 = (xo[1]-xi[1])/(yo[1]-yi[1]);
- c1 = xo[1]-0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(xo[1]-xi[1])+
- SQ(yo[1]-yi[1]))/(xo[1]-xi[1]);
- y = TSCRoutA-2.*TSCarbonFiberThA;
- x = p1*(y-yo[1])+c1;
- Ah1->SetVertex(1,x,y);
- Bh1->SetVertex(1,x,y);
- Ch1->SetVertex(5,x,y);
- y = TSCRinA+TSCarbonFiberThA;
- x = p1*(y-yo[1])+c1;
- Ah1->SetVertex(2,x,y);
- Bh1->SetVertex(2,x,y);
- Ch1->SetVertex(6,x,y);
- //
- // The easist way to get the points for the hole in volume A2 is to
- // rotate it to the Y axis where the y coordinates are easier to know
- // and then rotate it back.
- Double_t xp,yp,xa,ya,xb,yb;
- th = 0.5*TSCAngle*kRadian;
- xa = TMath::Cos(th)*xo[1]-TMath::Sin(th)*yo[1];
- ya = TMath::Sin(th)*xo[1]+TMath::Cos(th)*yo[1];
- xb = TMath::Cos(th)*xi[1]-TMath::Sin(th)*yi[1];
- yb = TMath::Sin(th)*xi[1]+TMath::Cos(th)*yi[1];
- p1 = (xa-xb)/(ya-yb);
- c1 = xa+0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(xa-xb)+SQ(ya-yb))/(xa-xb);
- y = ya-TSCarbonFiberThA;
- x = p1*(y-ya)+c1;
- xp = TMath::Cos(-th)*x-TMath::Sin(-th)*y;
- yp = TMath::Sin(-th)*x+TMath::Cos(-th)*y;
- Ah2->SetVertex(0,xp,yp);
- Bh2->SetVertex(0,xp,yp);
- Ch2->SetVertex(4,xp,yp);
- y = yb+2.0*TSCarbonFiberThA;
- x = p1*(y-ya)+c1;
- xp = TMath::Cos(-th)*x-TMath::Sin(-th)*y;
- yp = TMath::Sin(-th)*x+TMath::Cos(-th)*y;
- Ah2->SetVertex(3,xp,yp);
- Bh2->SetVertex(3,xp,yp);
- x7[1] = x; y7[1] = y; // vortexing done after last point
- //Ch2->SetVertex(7,xp,yp);
- xa = TMath::Cos(th)*xo[2]-TMath::Sin(th)*yo[2];
- ya = TMath::Sin(th)*xo[2]+TMath::Cos(th)*yo[2];
- xb = TMath::Cos(th)*xi[2]-TMath::Sin(th)*yi[2];
- yb = TMath::Sin(th)*xi[2]+TMath::Cos(th)*yi[2];
- p1 = (xa-xb)/(ya-yb);
- c1 = xa-0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(xa-xb)+SQ(ya-yb))/(xa-xb);
- y = ya-TSCarbonFiberThA;
- x = p1*(y-ya)+c1;
- xp = TMath::Cos(-th)*x-TMath::Sin(-th)*y;
- yp = TMath::Sin(-th)*x+TMath::Cos(-th)*y;
- Ah2->SetVertex(1,xp,yp);
- Bh2->SetVertex(1,xp,yp);
- Ch2->SetVertex(5,xp,yp);
- y = yb+2.0*TSCarbonFiberThA;
- x = p1*(y-ya)+c1;
- xp = TMath::Cos(-th)*x-TMath::Sin(-th)*y;
- yp = TMath::Sin(-th)*x+TMath::Cos(-th)*y;
- Ah2->SetVertex(2,xp,yp);
- Bh2->SetVertex(2,xp,yp);
- Ch2->SetVertex(6,xp,yp);
- //
- p1 = (yo[5]-yi[5])/(xo[5]-xi[5]);
- c1 = yo[5]+0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(yo[5]-yi[5])+
- SQ(xo[5]-xi[5]))/(yo[5]-yi[5]);
- x = xo[5]-TSCarbonFiberThA;
- y = p1*(x-xo[5])+c1;
- Ah3->SetVertex(0,x,y);
- Bh3->SetVertex(0,x,y);
- Ch3->SetVertex(4,x,y);
- x = xi[5]+2.0*TSCarbonFiberThA;
- y = p1*(x-xo[5])+c1;
- Ah3->SetVertex(3,x,y);
- Bh3->SetVertex(3,x,y);
- x7[2] = x; y7[2] = y; // vortexing done after last point
- //Ch3->SetVertex(7,x,y);
- y = 2.0*TSCarbonFiberThA;
- x = xo[5]-TSCarbonFiberThA;
- Ah3->SetVertex(1,x,y);
- Bh3->SetVertex(1,x,y);
- Ch3->SetVertex(5,x,y);
- y = 2.0*TSCarbonFiberThA;
- x = xi[5]+2.0*TSCarbonFiberThA;
- Ah3->SetVertex(2,x,y);
- Bh3->SetVertex(2,x,y);
- Ch3->SetVertex(6,x,y);
- //
- for(i=0;i<4;i++){ // define points at +dz
- A1->SetVertex(i+4,(A1->GetVertices())[2*i],(A1->GetVertices())[1+2*i]);
- A2->SetVertex(i+4,(A2->GetVertices())[2*i],(A2->GetVertices())[1+2*i]);
- A3->SetVertex(i+4,(A3->GetVertices())[2*i],(A3->GetVertices())[1+2*i]);
- //
- B1->SetVertex(i+4,(B1->GetVertices())[2*i],(B1->GetVertices())[1+2*i]);
- B2->SetVertex(i+4,(B2->GetVertices())[2*i],(B2->GetVertices())[1+2*i]);
- B3->SetVertex(i+4,(B3->GetVertices())[2*i],(B3->GetVertices())[1+2*i]);
- // C's are a cone which must match up with B's.
- C1->SetVertex(i+4,(B1->GetVertices())[2*i],(B1->GetVertices())[1+2*i]);
- C2->SetVertex(i+4,(B2->GetVertices())[2*i],(B2->GetVertices())[1+2*i]);
- C3->SetVertex(i+4,(B3->GetVertices())[2*i],(B3->GetVertices())[1+2*i]);
- //
- Ah1->SetVertex(i+4,(Ah1->GetVertices())[2*i],
- (Ah1->GetVertices())[1+2*i]);
- Ah2->SetVertex(i+4,(Ah2->GetVertices())[2*i],
- (Ah2->GetVertices())[1+2*i]);
- Ah3->SetVertex(i+4,(Ah3->GetVertices())[2*i],
- (Ah3->GetVertices())[1+2*i]);
- //
- Bh1->SetVertex(i+4,(Bh1->GetVertices())[2*i],
- (Bh1->GetVertices())[1+2*i]);
- Bh2->SetVertex(i+4,(Bh2->GetVertices())[2*i],
- (Bh2->GetVertices())[1+2*i]);
- Bh3->SetVertex(i+4,(Bh3->GetVertices())[2*i],
- (Bh3->GetVertices())[1+2*i]);
- } // end for
- //
- p1 = (xco[0]-xci[0])/(yco[0]-yci[0]);
- c1 = xco[0]+0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(xco[0]-xci[0])+
- SQ(yco[0]-yci[0]))/(xco[0]-xci[0]);
- y = TSCRoutC-2.*TSCarbonFiberThA;
- x = p1*(y-yco[0])+c1;
- Ch1->SetVertex(0,x,y);
- y = TSCRinC+TSCarbonFiberThA;
- x = p1*(y-yci[0])+c1;
- Ch1->SetVertex(2,x,y);
- p1 = (xco[1]-xci[1])/(yco[1]-yci[1]);
- c1 = xco[1]-0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(xco[1]-xci[1])+
- SQ(yco[1]-yci[1]))/(xco[1]-xci[1]);
- y = TSCRoutC-2.*TSCarbonFiberThA;
- x = p1*(y-yco[1])+c1;
- Ch1->SetVertex(1,x,y);
- y = TSCRinC+TSCarbonFiberThA;
- x = p1*(y-yci[1])+c1;
- Ch1->SetVertex(3,x,y);
- //
- th = 0.5*TSCAngle*kRadian;
- xa = TMath::Cos(th)*xco[1]-TMath::Sin(th)*yco[1];
- ya = TMath::Sin(th)*xco[1]+TMath::Cos(th)*yco[1];
- xb = TMath::Cos(th)*xci[1]-TMath::Sin(th)*yci[1];
- yb = TMath::Sin(th)*xci[1]+TMath::Cos(th)*yci[1];
- p1 = (xa-xb)/(ya-yb);
- c1 = xa+0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(xa-xb)+SQ(ya-yb))/(xa-xb);
- y = ya-TSCarbonFiberThA;
- x = p1*(y-ya)+c1;
- xp = TMath::Cos(-th)*x-TMath::Sin(-th)*y;
- yp = TMath::Sin(-th)*x+TMath::Cos(-th)*y;
- yp = ya-TSCarbonFiberThA;
- xp = p1*(y-ya)+c1;
- Ch2->SetVertex(0,xp,yp);
- y = yb+2.0*TSCarbonFiberThA;
- x = p1*(y-ya)+c1;
- xp = TMath::Cos(-th)*x-TMath::Sin(-th)*y;
- yp = TMath::Sin(-th)*x+TMath::Cos(-th)*y;
- Ch2->SetVertex(2,xp,yp);
- xa = TMath::Cos(th)*xco[2]-TMath::Sin(th)*yco[2];
- ya = TMath::Sin(th)*xco[2]+TMath::Cos(th)*yco[2];
- xb = TMath::Cos(th)*xci[2]-TMath::Sin(th)*yci[2];
- yb = TMath::Sin(th)*xci[2]+TMath::Cos(th)*yci[2];
- p1 = (xa-xb)/(ya-yb);
- c1 = xa-0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(xa-xb)+SQ(ya-yb))/(xa-xb);
- y = ya-TSCarbonFiberThA;
- x = p1*(y-ya)+c1;
- xp = TMath::Cos(-th)*x-TMath::Sin(-th)*y;
- yp = TMath::Sin(-th)*x+TMath::Cos(-th)*y;
- Ch2->SetVertex(1,xp,yp);
- y = yb+2.0*TSCarbonFiberThA;
- x = p1*(y-ya)+c1;
- xp = TMath::Cos(-th)*x-TMath::Sin(-th)*y;
- yp = TMath::Sin(-th)*x+TMath::Cos(-th)*y;
- Ch2->SetVertex(3,xp,yp);
- //
- p1 = (yco[5]-yci[5])/(xco[5]-xci[5]);
- c1 = yco[5]+0.5*TSCarbonFiberThA*TMath::Sqrt(SQ(yco[5]-yci[5])+
- SQ(xco[5]-xci[5]))/(yco[5]-yci[5]);
- x = xco[5]-TSCarbonFiberThA;
- y = p1*(x-xco[5])+c1;
- Ch3->SetVertex(0,x,y);
- x = xci[5]+2.0*TSCarbonFiberThA;
- y = p1*(x-xci[5])+c1;
- Ch3->SetVertex(2,x,y);
- y = 2.0*TSCarbonFiberThA;
- x = xco[5]-TSCarbonFiberThA;
- Ch3->SetVertex(1,x,y);
- y = 2.0*TSCarbonFiberThA;
- x = xci[5]+2.0*TSCarbonFiberThA;
- Ch3->SetVertex(3,x,y);
- Ch1->SetVertex(7,x7[0],y7[0]); // 7th point most be done last ???
- Ch2->SetVertex(7,x7[1],y7[1]); // 7th point most be done last ???
- Ch3->SetVertex(7,x7[2],y7[2]); // 7th point most be done last ???
- printArb8(A1);
- printArb8(Ah1);
- printArb8(A2);
- printArb8(Ah2);
- printArb8(A3);
- printArb8(Ah3);
- printArb8(B1);
- printArb8(Bh1);
- printArb8(B2);
- printArb8(Bh2);
- printArb8(B3);
- printArb8(Bh3);
- printArb8(C1);
- printArb8(Ch1);
- printArb8(C2);
- printArb8(Ch2);
- printArb8(C3);
- printArb8(Ch3);
- //
- // Define Minimal volume to inclose this SPD Thermal Sheald.
- M1 = new TGeoPcon("ITSspdShealdVV",0.0,360.0,9);
- M1->Z(0) = 0.5*TSCLengthA+TSCLengthB;
- M1->Rmin(0) = TSCRinB;
- x = B1->GetVertices()[0]; // [0][0]
- y = B1->GetVertices()[1]; // [0][1]
- M1->Rmax(0) = TMath::Sqrt(x*x+y*y);
- M1->Z(1) = M1->GetZ(0)-TSCLengthB;
- M1->Rmin(1) = M1->GetRmin(0);
- M1->Rmax(1) = M1->GetRmax(0);
- M1->Z(2) = M1->GetZ(1);
- M1->Rmin(2) = TSCRinA;
- x = A1->GetVertices()[0]; // [0]0]
- y = A1->GetVertices()[1]; // [0][1]
- M1->Rmax(2) = TMath::Sqrt(x*x+y*y);
- M1->Z(3) = -(M1->GetZ(0)-TSCLengthB);
- M1->Rmin(3) = M1->GetRmin(2);
- M1->Rmax(3) = M1->GetRmax(2);
- M1->Z(4) = M1->GetZ(3);
- M1->Rmin(4) = M1->GetRmin(1);
- M1->Rmax(4) = M1->GetRmax(1);
- M1->Z(5) = -(M1->GetZ(0));
- M1->Rmin(5) = M1->GetRmin(0);
- M1->Rmax(5) = M1->GetRmax(0);
- M1->Z(6) = M1->GetZ(5) - TSCLengthC;
- M1->Rmin(6) = TSCRinC;
- x = C1->GetVertices()[0]; // [0][0]
- y = C1->GetVertices()[1]; // [0][1]
- M1->Rmax(6) = TMath::Sqrt(x*x+y*y);
- M1->Z(7) = M1->GetZ(6);
- M1->Rmin(7) = D->GetRmin();
- M1->Rmax(7) = D->GetRmax();
- M1->Z(8) = M1->Z(7) - TSCLengthD;
- M1->Rmin(8) = M1->GetRmin(7);
- M1->Rmax(8) = M1->GetRmax(7);
- M2 = new TGeoTubeSeg("ITSspdShealdWingVV",
- M1->GetRmax(8),Dw->GetRmax(),Dw->GetDz(),Dw->GetPhi1(),Dw->GetPhi2());
- printTubeSeg(M2);
- //
- x = 0.5*(M1->GetZ(8) + M1->GetZ(7));
- tranITSspdShealdVVt0 = new TGeoTranslation("ITSspdShealdVVt0",0.0,0.0,x);
- tranITSspdShealdVVt0->RegisterYourself();
- TGeoRotation rotz90("",0.0,0.0,90.0); // never registered.
- rotITSspdShealdVVt1 = new TGeoCombiTrans(*tranITSspdShealdVVt0,rotz90);
- rotITSspdShealdVVt1->SetName("ITSspdShealdVVt1");
- rotITSspdShealdVVt1->RegisterYourself();
- TGeoRotation rotz180("",0.0,0.0,180.0); // never registered
- rotITSspdShealdVVt2 = new TGeoCombiTrans(*tranITSspdShealdVVt0,rotz180);
- rotITSspdShealdVVt2->SetName("ITSspdShealdVVt2");
- rotITSspdShealdVVt2->RegisterYourself();
- TGeoRotation rotz270("",0.0,0.0,270.0); // never registered
- rotITSspdShealdVVt3 = new TGeoCombiTrans(*tranITSspdShealdVVt0,rotz270);
- rotITSspdShealdVVt3->SetName("ITSspdShealdVVt3");
- rotITSspdShealdVVt3->RegisterYourself();
- M = new TGeoCompositeShape("ITS SPD Thermal sheald volume",
- "(((ITSspdShealdVV+"
- "ITSspdShealdWingVV:ITSspdShealdVVt0)+"
- "ITSspdShealdWingVV:ITSspdShealdVVt1)+"
- "ITSspdShealdWingVV:ITSspdShealdVVt2)+"
- "ITSspdShealdWingVV:ITSspdShealdVVt3");
- //
- TGeoManager *mgr = gGeoManager;
- SPDcf = mgr->GetMedium("ITSspdCarbonFiber");
- SPDfs = mgr->GetMedium("ITSspdStaselite4411w");
- SPDfo = mgr->GetMedium("ITSspdRohacell50A");
- SPDss = mgr->GetMedium("ITSspdStainlessSteal");
- SPDair= mgr->GetMedium("ITSspdAir");
- TGeoVolume *A1v,*A2v,*A3v,*Ah1v,*Ah2v,*Ah3v;
- TGeoVolume *B1v,*B2v,*B3v,*Bh1v,*Bh2v,*Bh3v;
- TGeoVolume *C1v,*C2v,*C3v,*Ch1v,*Ch2v,*Ch3v;
- TGeoVolume *Dv,*Dsv,*Dwv,*Dwsv,*Mv;
- Mv = new TGeoVolume("ITSspdThermalSheald",M,SPDair);
- Mv->SetVisibility(kTRUE);
- Mv->SetLineColor(7); // light Blue
- Mv->SetLineWidth(1);
- Mv->SetFillColor(Mv->GetLineColor());
- Mv->SetFillStyle(4090); // 90% transparent
- Moth->AddNode(Mv,1,0); ///////////////////// Virtual Volume ////////
- A1v = new TGeoVolume("ITSspdCentCylA1CF",A1,SPDcf);
- A1v->SetVisibility(kTRUE);
- A1v->SetLineColor(4);
- A1v->SetLineWidth(1);
- A2v = new TGeoVolume("ITSspdCentCylA2CF",A2,SPDcf);
- A2v->SetVisibility(kTRUE);
- A2v->SetLineColor(4);
- A2v->SetLineWidth(1);
- A3v = new TGeoVolume("ITSspdCentCylA3CF",A3,SPDcf);
- A3v->SetVisibility(kTRUE);
- A3v->SetLineColor(4);
- A3v->SetLineWidth(1);
- B1v = new TGeoVolume("ITSspdCentCylB1CF",B1,SPDcf);
- B1v->SetVisibility(kTRUE);
- B1v->SetLineColor(4);
- B1v->SetLineWidth(1);
- B2v = new TGeoVolume("ITSspdCentCylB2CF",B2,SPDcf);
- B2v->SetVisibility(kTRUE);
- B2v->SetLineColor(4);
- B2v->SetLineWidth(1);
- B3v = new TGeoVolume("ITSspdCentCylB3CF",B3,SPDcf);
- B3v->SetVisibility(kTRUE);
- B3v->SetLineColor(4);
- B3v->SetLineWidth(1);
- C1v = new TGeoVolume("ITSspdCentCylC1CF",C1,SPDcf);
- C1v->SetVisibility(kTRUE);
- C1v->SetLineColor(4);
- C1v->SetLineWidth(1);
- C2v = new TGeoVolume("ITSspdCentCylC2CF",C2,SPDcf);
- C2v->SetVisibility(kTRUE);
- C2v->SetLineColor(4);
- C2v->SetLineWidth(1);
- C3v = new TGeoVolume("ITSspdCentCylC3CF",C3,SPDcf);
- C3v->SetVisibility(kTRUE);
- C3v->SetLineColor(4);
- C3v->SetLineWidth(1);
- Ah1v = new TGeoVolume("ITSspdCentCylA1AirA",Ah1,SPDair);
- Ah1v->SetVisibility(kTRUE);
- Ah1v->SetLineColor(5); // Yellow
- Ah1v->SetFillColor(Ah1v->GetLineColor());
- Ah1v->SetFillStyle(4090); // 90% transparent
- Ah2v = new TGeoVolume("ITSspdCentCylA2AirA",Ah2,SPDair);
- Ah2v->SetVisibility(kTRUE);
- Ah2v->SetLineColor(5); // Yellow
- Ah2v->SetFillColor(Ah2v->GetLineColor());
- Ah2v->SetFillStyle(4090); // 90% transparent
- Ah3v = new TGeoVolume("ITSspdCentCylA3AirA",Ah3,SPDair);
- Ah3v->SetVisibility(kTRUE);
- Ah3v->SetLineColor(5); // Yellow
- Ah3v->SetFillColor(Ah3v->GetLineColor());
- Ah3v->SetFillStyle(4090); // 90% transparent
- Bh1v = new TGeoVolume("ITSspdCentCylA1AirB",Bh1,SPDair);
- Bh1v->SetVisibility(kTRUE);
- Bh1v->SetLineColor(5); // Yellow
- Bh1v->SetFillColor(Bh1v->GetLineColor());
- Bh1v->SetFillStyle(4090); // 90% transparent
- Bh2v = new TGeoVolume("ITSspdCentCylA2AirB",Bh2,SPDair);
- Bh2v->SetVisibility(kTRUE);
- Bh2v->SetLineColor(5); // Yellow
- Bh2v->SetFillColor(Bh2v->GetLineColor());
- Bh2v->SetFillStyle(4090); // 90% transparent
- Bh3v = new TGeoVolume("ITSspdCentCylA3AirB",Bh3,SPDair);
- Bh3v->SetVisibility(kTRUE);
- Bh3v->SetLineColor(5); // Yellow
- Bh3v->SetFillColor(Bh3v->GetLineColor());
- Bh3v->SetFillStyle(4090); // 90% transparent
- Ch1v = new TGeoVolume("ITSspdCentCylA1AirC",Ch1,SPDair);
- Ch1v->SetVisibility(kTRUE);
- Ch1v->SetLineColor(5); // Yellow
- Ch1v->SetFillColor(Ch1v->GetLineColor());
- Ch1v->SetFillStyle(4090); // 90% transparent
- Ch2v = new TGeoVolume("ITSspdCentCylA2AirC",Ch2,SPDair);
- Ch2v->SetVisibility(kTRUE);
- Ch2v->SetLineColor(5); // Yellow
- Ch2v->SetFillColor(Ch2v->GetLineColor());
- Ch2v->SetFillStyle(4090); // 90% transparent
- Ch3v = new TGeoVolume("ITSspdCentCylA3AirC",Ch3,SPDair);
- Ch3v->SetVisibility(kTRUE);
- Ch3v->SetLineColor(5); // Yellow
- Ch3v->SetFillColor(Ch3v->GetLineColor());
- Ch3v->SetFillStyle(4090); // 90% transparent
- Dv = new TGeoVolume("ITSspdCentCylA1CD",D,SPDcf);
- Dv->SetVisibility(kTRUE);
- Dv->SetLineColor(4);
- Dv->SetLineWidth(1);
- Dwv = new TGeoVolume("ITSspdCentCylA1CDw",Dw,SPDcf);
- Dwv->SetVisibility(kTRUE);
- Dwv->SetLineColor(4);
- Dwv->SetLineWidth(1);
- Dsv = new TGeoVolume("ITSspdCentCylA1Dfill",Ds,SPDfs);
- Dsv->SetVisibility(kTRUE);
- Dsv->SetLineColor(3); // Green
- Dsv->SetFillColor(Dsv->GetLineColor());
- Dsv->SetFillStyle(4010); // 10% transparent
- Dwsv = new TGeoVolume("ITSspdCentCylA1DwingFill",Dws,SPDfs);
- Dwsv->SetVisibility(kTRUE);
- Dwsv->SetLineColor(3); // Green
- Dwsv->SetFillColor(Dwsv->GetLineColor());
- Dwsv->SetFillStyle(4010); // 10% transparent
- //
- A1v->AddNode(Ah1v,1,0);
- A2v->AddNode(Ah2v,1,0);
- A3v->AddNode(Ah3v,1,0);
- B1v->AddNode(Bh1v,1,0);
- B2v->AddNode(Bh2v,1,0);
- B3v->AddNode(Bh3v,1,0);
- C1v->AddNode(Ch1v,1,0);
- C2v->AddNode(Ch2v,1,0);
- C3v->AddNode(Ch3v,1,0);
- Dv ->AddNode(Dsv ,1,0);
- Dwv->AddNode(Dwsv,1,0);
- //
- Mv->AddNode(A1v,1,0);
- Mv->AddNode(A2v,1,0);
- Mv->AddNode(A3v,1,0);
- tranb = new TGeoTranslation("",0.0,0.0,0.5*(TSCLengthA+TSCLengthB));
- tranbm = new TGeoTranslation("",0.0,0.0,0.5*(-TSCLengthA-TSCLengthB));
- Mv->AddNode(B1v,1,tranb);
- Mv->AddNode(B2v,1,tranb);
- Mv->AddNode(B3v,1,tranb);
- Mv->AddNode(B1v,2,tranbm);
- Mv->AddNode(B2v,2,tranbm);
- Mv->AddNode(B3v,2,tranbm);
- // Muon side (rb26) is at -Z.
- tranc = new TGeoTranslation("",0.0,0.0,
- 0.5*(-TSCLengthA-TSCLengthB-TSCLengthC));
- Mv->AddNode(C1v,1,tranc);
- Mv->AddNode(C2v,1,tranc);
- Mv->AddNode(C3v,1,tranc);
- Mv->AddNode(Dv,1,tranITSspdShealdVVt0);
- Mv->AddNode(Dwv,1,tranITSspdShealdVVt0);
- Mv->AddNode(Dwv,2,rotITSspdShealdVVt1);
- Mv->AddNode(Dwv,3,rotITSspdShealdVVt2);
- Mv->AddNode(Dwv,4,rotITSspdShealdVVt3);
- k=2;
- for(i=1;i<10;i++) {
- th = ((Double_t)i)*TSCAngle*kDegree;
- rot = new TGeoRotation("",0.0,0.0,th);
- Mv->AddNode(A1v,i+1,rot);
- Mv->AddNode(B1v,i+2,new TGeoCombiTrans(*tranb,*rot));
- Mv->AddNode(B1v,i+12,new TGeoCombiTrans(*tranbm,*rot));
- Mv->AddNode(C1v,i+1,new TGeoCombiTrans(*tranc,*rot));
- if(i!=0||i!=2||i!=7){
- Mv->AddNode(A2v,k++,rot);
- Mv->AddNode(B2v,k++,new TGeoCombiTrans(*tranb,*rot));
- Mv->AddNode(B2v,k++,new TGeoCombiTrans(*tranbm,*rot));
- Mv->AddNode(C2v,k++,new TGeoCombiTrans(*tranc,*rot));
- } // end if
- if(i==5) {
- Mv->AddNode(A3v,2,rot);
- Mv->AddNode(B3v,3,new TGeoCombiTrans(*tranb,*rot));
- Mv->AddNode(B3v,4,new TGeoCombiTrans(*tranbm,*rot));
- Mv->AddNode(C3v,2,new TGeoCombiTrans(*tranc,*rot));
- } // end if
- } // end for i
- rot = new TGeoRotation("",180.,0.0,0.0);
- Mv->AddNode(A3v,3,rot);
- Mv->AddNode(B3v,5,new TGeoCombiTrans(*tranb,*rot));
- Mv->AddNode(B3v,6,new TGeoCombiTrans(*tranbm,*rot));
- Mv->AddNode(C3v,3,new TGeoCombiTrans(*tranc,*rot));
- rot = new TGeoRotation("",180.,0.0,180.0);
- Mv->AddNode(A3v,4,rot);
- Mv->AddNode(B3v,7,new TGeoCombiTrans(*tranb,*rot));
- Mv->AddNode(B3v,8,new TGeoCombiTrans(*tranbm,*rot));
- Mv->AddNode(C3v,4,new TGeoCombiTrans(*tranc,*rot));
- if(GetDebug()){
- A1v->PrintNodes();
- Ah1v->PrintNodes();
- A2v->PrintNodes();
- Ah2v->PrintNodes();
- A3v->PrintNodes();
- Ah3v->PrintNodes();
- B1v->PrintNodes();
- Bh1v->PrintNodes();
- B2v->PrintNodes();
- Bh2v->PrintNodes();
- B3v->PrintNodes();
- Bh3v->PrintNodes();
- C1v->PrintNodes();
- Ch1v->PrintNodes();
- C2v->PrintNodes();
- Ch2v->PrintNodes();
- C3v->PrintNodes();
- Ch3v->PrintNodes();
- Dv->PrintNodes();
- Dsv->PrintNodes();
- Dwv->PrintNodes();
- Dwsv->PrintNodes();
- //Mv->PrintNodes();
- } // end if