+ pRich->AddNode(pF4 ,copy=1,new TGeoTranslation( 0*mm,0*mm, 9*mm)); //F4 to Rich p.84 TDR
+ pF4 ->AddNode(pF4al ,copy=1,new TGeoTranslation( 0*mm,0*mm, 2.5*mm)); //F4 al to F4 2043P1
+ pF4al->AddNode(pF4in ,copy=1,new TGeoTranslation( 0*mm,0*mm, 0*mm)); //F4 whole F4 al 2043P1
+ for(int i=1;i<=322;i++)
+ pF4->AddNode(pColl,copy=i,new TGeoCombiTrans( 0*mm, -1296/2*mm+i*4*mm, -5*mm,pCollRot)); //F4 2043P1
+//Radiators or Sr90 source
+ if(!isSr90){//Radiators
+ TGeoVolume *pRad =g->MakeBox( "Rad" ,g->GetMedium("C6F14") ,dx=1330*mm/2 ,dy= 413*mm/2 ,dz= 24*mm/2); // Rad 2011P1
+ TGeoVolume *pRadFront =g->MakeBox( "RadFront" ,g->GetMedium("Neoceram") ,dx=1330*mm/2 ,dy= 413*mm/2 ,dz= 4*mm/2);
+ TGeoVolume *pRadWin =g->MakeBox( "RadWin" ,g->GetMedium("SiO2") ,dx=1330*mm/2 ,dy= 413*mm/2 ,dz= 5*mm/2);
+ TGeoVolume *pRadLong =g->MakeBox( "RadLong" ,g->GetMedium("Neoceram") ,dx=1330*mm/2 ,dy= 5*mm/2 ,dz= 15*mm/2);
+ TGeoVolume *pRadShort =g->MakeBox( "RadShort" ,g->GetMedium("Neoceram") ,dx= 10*mm/2 ,dy= 403*mm/2 ,dz= 15*mm/2);
+ TGeoVolume *pRadSpacer=g->MakeTube("RadSpacer",g->GetMedium("SiO2") ,r1= 0 ,r2=10*mm/2 ,dz= 15*mm/2);
+
+ pRich->AddNode(pRad,copy=1,new TGeoTranslation( 0*mm,-434*mm, -12*mm));
+ pRich->AddNode(pRad,copy=2,new TGeoTranslation( 0*mm, 0*mm, -12*mm));
+ pRich->AddNode(pRad,copy=3,new TGeoTranslation( 0*mm,+434*mm, -12*mm));
+ pRad->AddNode(g->GetVolume("RadFront") ,copy=1,new TGeoTranslation( 0*mm, 0*mm, -10.0*mm));
+ pRad->AddNode(g->GetVolume("RadWin") ,copy=1,new TGeoTranslation( 0*mm, 0*mm, 9.5*mm));
+ pRad->AddNode(g->GetVolume("RadLong") ,copy=1,new TGeoTranslation( 0*mm,-204*mm, -0.5*mm));
+ pRad->AddNode(g->GetVolume("RadLong") ,copy=2,new TGeoTranslation( 0*mm,+204*mm, -0.5*mm));
+ pRad->AddNode(g->GetVolume("RadShort") ,copy=1,new TGeoTranslation(-660*mm, 0*mm, -0.5*mm));
+ pRad->AddNode(g->GetVolume("RadShort") ,copy=2,new TGeoTranslation(+660*mm, 0*mm, -0.5*mm));
+ for(int i=0;i<3;i++)
+ for(int j=0;j<10;j++)
+ pRad->AddNode(pRadSpacer,copy=10*i+j,new TGeoTranslation(-1330*mm/2+116*mm+j*122*mm,(i-1)*105*mm,-0.5*mm));
+ }else{//Sr90
+ Double_t containerZ=21.8*mm/2,containerR=38*mm/2;
+ Double_t screwLen=15*mm/2 ,screwR=5*mm/2;
+ Double_t srcLen=10*mm/2 ,srcR=2*mm/2;
+ Double_t perpexZ=20*mm/2 ,perpexR=34*mm/2;
+ Double_t perpexWholeLen=10*mm/2,perpexWholeR=4*mm/2;
+ Double_t alBottomLen=containerZ-perpexZ,alWholeR=5*mm/2;
+ g->MakeTube("Sr90" ,pAir ,0 ,containerR ,containerZ);
+ g->MakeTube("Sr90AlWall" ,pAl ,perpexR,containerR ,containerZ);
+ g->MakeTube("Sr90AlBottom" ,pAl ,0 ,perpexR ,alBottomLen);
+ g->MakeTube("Sr90AlWhole" ,pAir ,0 ,alWholeR ,alBottomLen);
+ g->MakeTube("Sr90PerpexPlug" ,g->GetMedium("Perpex"),0 ,perpexR ,perpexZ);
+ g->MakeTube("Sr90PerpexWhole",pAir ,0 ,perpexWholeR,perpexWholeLen);
+ g->MakeTube("Sr90Screw" ,g->GetMedium("Steel") ,0 ,screwR ,screwLen);
+ g->MakeTube("Sr90Source" ,g->GetMedium("Steel") ,0 ,srcR ,srcLen);
+
+ pRich->AddNode(pContainer,1,new TGeoTranslation(30*cm,0,containerZ));
+ pContainer ->AddNode(pAlWall,1);
+ pContainer ->AddNode(pAlBottom,1,new TGeoTranslation(0,0,-containerZ+alBottomLen));
+ pContainer ->AddNode(pPerpexPlug,1,new TGeoTranslation(0,0,containerZ-perpexZ));
+ pAlBottom ->AddNode(pAlWhole,1,new TGeoTranslation(6*mm,0,0));
+ pPerpexPlug->AddNode(pPerpexWhole,1,new TGeoTranslation(6*mm,0,-perpexZ+perpexWholeLen));
+ pPerpexPlug->AddNode(pSource,1 ,new TGeoTranslation(6*mm,0, perpexZ-srcLen));
+ pPerpexPlug->AddNode(pScrew,1,new TGeoTranslation(0,0,perpexZ-screwLen));
+ }//if(!isSr90)
+//Sandbox
+ TGeoVolume *pSandBox=g->MakeBox( "SandBox" ,g->GetMedium("Air") ,dx=1419*mm/2 ,dy=1378*mm/2 ,dz=50.5*mm/2); //2072P1
+ g->MakeBox( "SandCover",g->GetMedium("Al") ,dx=1419*mm/2 ,dy=1378*mm/2 ,dz= 0.5*mm/2);
+ g->MakeBox( "SandComb" ,g->GetMedium("Rohacell") ,dx=1359*mm/2 ,dy=1318*mm/2 ,dz=49.5*mm/2);