HMPID/Hgeom.C

   1
   2 Int_t copy; //volume copy number
   3 Double_t dx,dy,dz,r1,r2;//tmp vars for volume dimentions
   4 Double_t cm=1,m=100*cm,mm=0.1*cm,mkm=0.001*cm;//length units
   5 TGeoManager *g=0;
   6
   7 void Hgeom(Bool_t isOnlyChambers=kFALSE)
   8 {
   9
  10   g=new TGeoManager("HMPID","Private HMPID geometry");
  11   Materials();
  12   gGeoManager->MakeBox("ALIC",gGeoManager->GetMedium("Air"),dx=30*m/2,dy=30*m/2,dz=30*m/2); //arbitrary values
  13   gGeoManager->SetTopVolume(gGeoManager->GetVolume("ALIC"));
  14
  15   Hmpid(isOnlyChambers);
  16
  17   gGeoManager->CloseGeometry();
  18
  19   gGeoManager->SetVisOption(0); gGeoManager->SetVisLevel(5);
  20
  21   gGeoManager->GetMasterVolume()->Draw();
  22   Axis();
  23 //  gPad->GetView()->SetView(3,94,-70,0);
  24   new TBrowser;
  25 }
  26 //__________________________________________________________________________________________________
  27 void Materials()
  28 {
  29 //Media for HMPID
  30   Double_t a=0,z=0,den=0,radlen=0,intlen=0;//tmp vars for material parameters
  31   new TGeoMaterial("Air"           ,a=26.98 ,z=13,den=0.4224                             ); new TGeoMedium("Air"          ,1,gGeoManager->GetMaterial("Air"));
  32   new TGeoMaterial("HMPID_CH4"      ,a=26.98 ,z=13,den=0.4224                             ); new TGeoMedium("HMPID_CH4"     ,2,gGeoManager->GetMaterial("HMPID_CH4"));
  33   new TGeoMaterial("HMPID_CsI"      ,a=26.98 ,z=13,den=2.7 ,radlen=24.01*cm,intlen=70.6*cm); new TGeoMedium("HMPID_CsI"     ,3,gGeoManager->GetMaterial("HMPID_CsI"));
  34   new TGeoMaterial("HMPID_Al"       ,a=26.98 ,z=13,den=2.7 ,radlen=24.01*cm,intlen=70.6*cm); new TGeoMedium("HMPID_Al"      ,4,gGeoManager->GetMaterial("HMPID_Al"));
  35   new TGeoMaterial("HMPID_W"        ,a=183.84,z=27,den=19.3,radlen= 9.59*cm,intlen=0.35*cm); new TGeoMedium("HMPID_W"       ,5,gGeoManager->GetMaterial("HMPID_W"));
  36   new TGeoMaterial("HMPID_Cu"       ,a=55.845,z=26,den=7.87,radlen=13.84*cm,intlen=82.8*cm); new TGeoMedium("HMPID_Cu"      ,6,gGeoManager->GetMaterial("HMPID_Cu"));
  37   new TGeoMaterial("HMPID_Rohacell" ,a=12.01 ,z=6 ,den=0.1 ,radlen=18.8,intlen=0);           new TGeoMedium("HMPID_Rohacell",7,gGeoManager->GetMaterial("HMPID_Rohacell"));
  38   new TGeoMaterial("HMPID_SiO2"     ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_SiO2"    ,8,gGeoManager->GetMaterial("HMPID_SiO2"));
  39   new TGeoMaterial("HMPID_C6F14"    ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_C6F14"   ,9,gGeoManager->GetMaterial("HMPID_C6F14"));
  40 //Media for Sr90 source
  41   new TGeoMaterial("HMPID_Perpex"   ,a=55.845,z=26,den=7.87,radlen=13.84*cm,intlen=82.8*cm); new TGeoMedium("HMPID_Perpex"  ,10,gGeoManager->GetMaterial("HMPID_Perpex"));
  42   new TGeoMaterial("HMPID_Steel"    ,a=55.845,z=26,den=7.87,radlen=13.84*cm,intlen=82.8*cm); new TGeoMedium("HMPID_Steel"   ,11,gGeoManager->GetMaterial("HMPID_Steel"));
  43   new TGeoMaterial("HMPID_Mylar"    ,a=55.845,z=26,den=7.87,radlen=13.84*cm,intlen=82.8*cm); new TGeoMedium("HMPID_Mylar"   ,12,gGeoManager->GetMaterial("HMPID_Mylar"));
  44   new TGeoMaterial("HMPID_Sr90"     ,a=87.62 ,z=38,den=7.87,radlen=13.84*cm,intlen=82.8*cm); new TGeoMedium("HMPID_Sr90"    ,13,gGeoManager->GetMaterial("HMPID_Sr90"));
  45 //Media for VHMPID Gas option
  46   new TGeoMaterial("HMPID_CF4"      ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_CF4"     ,14,gGeoManager->GetMaterial("HMPID_CF4"));
  47   new TGeoMaterial("HMPID_C4F10"    ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_C4F10"   ,15,gGeoManager->GetMaterial("HMPID_C4F10"));
  48   new TGeoMaterial("HMPID_Ag"       ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_Ag"      ,16,gGeoManager->GetMaterial("HMPID_Ag"));
  49 //Media for VHMPID aerogel option
  50   new TGeoMaterial("HMPID_Gel24"    ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_Gel24"   ,17,gGeoManager->GetMaterial("HMPID_Gel24"));
  51   new TGeoMaterial("HMPID_Gel26"    ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_Gel26"   ,18,gGeoManager->GetMaterial("HMPID_Gel26"));
  52   new TGeoMaterial("HMPID_Gel28"    ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_Gel28"   ,19,gGeoManager->GetMaterial("HMPID_Gel28"));
  53   new TGeoMaterial("HMPID_Gel30"    ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_Gel30"   ,20,gGeoManager->GetMaterial("HMPID_Gel30"));
  54   new TGeoMaterial("HMPID_Si"       ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_Si"      ,21,gGeoManager->GetMaterial("HMPID_Si"));
  55   new TGeoMaterial("HMPID_Apd"      ,a=0     ,z=0 ,den=0);                                   new TGeoMedium("HMPID_Apd"     ,22,gGeoManager->GetMaterial("HMPID_Apd"));
  56 }//Materials()
  57 //__________________________________________________________________________________________________
  58 void Hmpid(Bool_t isOnlyChambers)
  59 {
  60 //HMPID  chamber
  61 //                                                         z ^
  62 //                                                           |
  63 //                                                           |
  64 //    <-Y0-> X <--Y1--> X <--Y1--> X <--Y1--> X <-Y0->       |            cath wires: r 50mkm; shift Y0=1.05m;, pitch Y1=2.1; center to PC  4.45mm; material Cu
  65 //                                                         x *------->y
  66 //
  67 //    <--Y0--> x <-----------Y1-----------> x <--Y0-->                    anod wires: r 20mkm; shift Y0=2.2mm; pitch Y1=4.0mm; center to PC  2.04mm; material W
  68 //
  69 //
  70 //   |________________________________________________|                   pad size y 8.4mm
  71 //
  72 //
  73   TGeoVolume *pRich=gGeoManager->MakeBox("HMPID",gGeoManager->GetMedium("Air"),dx=(6*mm+1681*mm+6*mm)/2,  //main HMPID volume
  74                                                                               dy=(6*mm+1466*mm+6*mm)/2,
  75                                                                               dz=(80*mm+40*mm)*2/2);     //x,y taken from 2033P1  z from p84 TDR
  76   const Double_t kAngHor=19.5; //  horizontal angle between chambers  19.5 grad
  77   const Double_t kAngVer=20;   //  vertical angle between chambers    20   grad
  78   const Double_t kAngCom=30;   //  common HMPID rotation with respect to x axis  30   grad
  79   const Double_t trans[3]={490,0,0}; //center of the chamber is on window-gap surface
  80   for(Int_t iCh=1;iCh<=7;iCh++){//place 7 chambers
  81     TGeoHMatrix *pMatrix=new TGeoHMatrix;
  82     pMatrix->RotateY(90);           //rotate around y since initial position is in XY plane -> now in YZ plane
  83     pMatrix->SetTranslation(trans); //now plane in YZ is shifted along x
  84     switch(iCh){
  85       case 1:                pMatrix->RotateY(kAngHor);  pMatrix->RotateZ(-kAngVer);  break; //right and down
  86       case 2:                                            pMatrix->RotateZ(-kAngVer);  break; //down
  87       case 3:                pMatrix->RotateY(kAngHor);                               break; //right
  88       case 4:                                                                         break; //no rotation
  89       case 5:                pMatrix->RotateY(-kAngHor);                              break; //left
  90       case 6:                                            pMatrix->RotateZ(kAngVer);   break; //up
  91       case 7:                pMatrix->RotateY(-kAngHor); pMatrix->RotateZ(kAngVer);   break; //left and up
  92     }
  93     pMatrix->RotateZ(kAngCom);     //apply common rotation  in XY plane
  94     gGeoManager->GetVolume("ALIC")->AddNode(pRich,iCh,pMatrix);
  95   }
  96   if(isOnlyChambers) return; //do not construct the detailed geometry
  97 //Pad Panel frame
  98   TGeoVolume *pPpf     =gGeoManager->MakeBox("Rppf"  ,gGeoManager->GetMedium("HMPID_Al")  ,dx=648*mm/2,dy=  411*mm/2 ,dz=40  *mm/2);//PPF 2001P2 inner size of the slab by 1mm more
  99   TGeoVolume *pPpfLarge=gGeoManager->MakeBox("Rppf1" ,gGeoManager->GetMedium("Air")      ,dx=181*mm/2,dy=89.25*mm/2 ,dz=38.3*mm/2);     //large whole
 100   TGeoVolume *pPpfSmall=gGeoManager->MakeBox("Rppf2" ,gGeoManager->GetMedium("Air")      ,dx=114*mm/2,dy=89.25*mm/2 ,dz=38.3*mm/2);//small whole
 101   TGeoVolume *pPc      =gGeoManager->MakeBox("Rpc"   ,gGeoManager->GetMedium("HMPID_CsI") ,dx=644*mm/2,dy=  407*mm/2 ,dz= 1.7*mm/2);//by 0.2 mm more then actual size (PCB 2006P1)
 102
 103   pRich->AddNode(pPpf,copy=1,new TGeoTranslation(-335*mm,-433*mm,8*cm+20*mm));//F1 2040P1 z p.84 TDR
 104   pRich->AddNode(pPpf,copy=2,new TGeoTranslation(+335*mm,-433*mm,8*cm+20*mm));
 105   pRich->AddNode(pPpf,copy=3,new TGeoTranslation(-335*mm,   0*mm,8*cm+20*mm));
 106   pRich->AddNode(pPpf,copy=4,new TGeoTranslation(+335*mm,   0*mm,8*cm+20*mm));
 107   pRich->AddNode(pPpf,copy=5,new TGeoTranslation(-335*mm,+433*mm,8*cm+20*mm));
 108   pRich->AddNode(pPpf,copy=6,new TGeoTranslation(+335*mm,+433*mm,8*cm+20*mm));
 109   pPpf->AddNode( pPc ,copy=1,new TGeoTranslation(   0*mm,   0*mm,-19.15*mm));//PPF 2001P2
 110   pPpf->AddNode(pPpfLarge,copy=1,new TGeoTranslation(-224.5*mm,-151.875*mm,  0.85*mm));
 111   pPpf->AddNode(pPpfLarge,copy=2,new TGeoTranslation(-224.5*mm,- 50.625*mm,  0.85*mm));
 112   pPpf->AddNode(pPpfLarge,copy=3,new TGeoTranslation(-224.5*mm,+ 50.625*mm,  0.85*mm));
 113   pPpf->AddNode(pPpfLarge,copy=4,new TGeoTranslation(-224.5*mm,+151.875*mm,  0.85*mm));
 114   pPpf->AddNode(pPpfSmall,copy=1,new TGeoTranslation(- 65.0*mm,-151.875*mm,  0.85*mm));
 115   pPpf->AddNode(pPpfSmall,copy=2,new TGeoTranslation(- 65.0*mm,- 50.625*mm,  0.85*mm));
 116   pPpf->AddNode(pPpfSmall,copy=3,new TGeoTranslation(- 65.0*mm,+ 50.625*mm,  0.85*mm));
 117   pPpf->AddNode(pPpfSmall,copy=4,new TGeoTranslation(- 65.0*mm,+151.875*mm,  0.85*mm));
 118   pPpf->AddNode(pPpfSmall,copy=5,new TGeoTranslation(+ 65.0*mm,-151.875*mm,  0.85*mm));
 119   pPpf->AddNode(pPpfSmall,copy=6,new TGeoTranslation(+ 65.0*mm,- 50.625*mm,  0.85*mm));
 120   pPpf->AddNode(pPpfSmall,copy=7,new TGeoTranslation(+ 65.0*mm,+ 50.625*mm,  0.85*mm));
 121   pPpf->AddNode(pPpfSmall,copy=8,new TGeoTranslation(+ 65.0*mm,+151.875*mm,  0.85*mm));
 122   pPpf->AddNode(pPpfLarge,copy=5,new TGeoTranslation(+224.5*mm,-151.875*mm,  0.85*mm));
 123   pPpf->AddNode(pPpfLarge,copy=6,new TGeoTranslation(+224.5*mm,- 50.625*mm,  0.85*mm));
 124   pPpf->AddNode(pPpfLarge,copy=7,new TGeoTranslation(+224.5*mm,+ 50.625*mm,  0.85*mm));
 125   pPpf->AddNode(pPpfLarge,copy=8,new TGeoTranslation(+224.5*mm,+151.875*mm,  0.85*mm));
 126 //gap - anod wires
 127   TGeoVolume *pGap =gGeoManager->MakeBox ("Rgap" ,gGeoManager->GetMedium("HMPID_CH4") ,dx=648*mm/2,dy=  411*mm/2 ,dz=4.45*mm/2);//xy as PPF 2001P2 z WP 2099P1
 128   TGeoVolume *pAnod=gGeoManager->MakeTube("Rano" ,gGeoManager->GetMedium("HMPID_W")   ,r1=  0*mm  ,r2=  20*mkm/2 ,dz=648*mm/2); //WP 2099P1 z = gap x PPF 2001P2
 129   TGeoRotation *pAnodRot=new TGeoRotation("RanW",90,90,0);
 130
 131   pRich->AddNode(pGap,copy=1,new TGeoTranslation(-335*mm,-433*mm,8*cm-2.225*mm)); //F1 2040P1 z WP 2099P1
 132   pRich->AddNode(pGap,copy=2,new TGeoTranslation(+335*mm,-433*mm,8*cm-2.225*mm));
 133   pRich->AddNode(pGap,copy=3,new TGeoTranslation(-335*mm,   0*mm,8*cm-2.225*mm));
 134   pRich->AddNode(pGap,copy=4,new TGeoTranslation(+335*mm,   0*mm,8*cm-2.225*mm));
 135   pRich->AddNode(pGap,copy=5,new TGeoTranslation(-335*mm,+433*mm,8*cm-2.225*mm));
 136   pRich->AddNode(pGap,copy=6,new TGeoTranslation(+335*mm,+433*mm,8*cm-2.225*mm));
 137   for(int i=1;i<=96;i++)
 138     pGap->AddNode(pAnod,copy=i,new TGeoCombiTrans( 0*mm, -411/2*mm+i*4*mm, 0.185*mm,pAnodRot)); //WP 2099P1
 139 //frame 3- cathode wires
 140   TGeoVolume *pCath=gGeoManager->MakeTube("RcaW"  ,gGeoManager->GetMedium("Cu")  ,r1=0  ,r2=100*mkm/2,dz=1323*mm/2);//r WP 2099P1 z F3 2041P1
 141   TGeoRotation *pCathRot=new TGeoRotation("CathRot",90,90,0);
 142   for(int i=1;i<=618;i++)
 143     pRich->AddNode(pCath,copy=i,new TGeoCombiTrans( 0*mm, -649.5*mm+i*2.1*mm, 75*mm,pCathRot)); //WP 2099P1
 144 //Frame 4- collection wires
 145   TGeoVolume *pF4  =gGeoManager->MakeBox( "Rfr4"       ,gGeoManager->GetMedium("HMPID_CH4")   ,dx=1407*mm/2 ,dy=1366*mm/2  ,dz=  15*mm/2);//F4 2043P1
 146   TGeoVolume *pF4al=gGeoManager->MakeBox( "Rfr4al"     ,gGeoManager->GetMedium("HMPID_Al")    ,dx=1407*mm/2 ,dy=1366*mm/2  ,dz=  10*mm/2);
 147   TGeoVolume *pF4in=gGeoManager->MakeBox( "Rfr4in"     ,gGeoManager->GetMedium("HMPID_CH4")   ,dx=1323*mm/2 ,dy=1296*mm/2  ,dz=  10*mm/2);
 148   TGeoVolume *pColl=gGeoManager->MakeTube("RcoW"       ,gGeoManager->GetMedium("HMPID_Cu")    ,r1=   0*mm   ,r2=100*mkm/2  ,dz=1323*mm/2);
 149   TGeoRotation *pCollRot=new TGeoRotation("RcoRot",90,90,0);
 150
 151   pRich->AddNode(pF4      ,copy=1,new TGeoTranslation(   0*mm,0*mm,   9*mm)); //F4 to Rich p.84 TDR
 152     pF4  ->AddNode(pF4al    ,copy=1,new TGeoTranslation(   0*mm,0*mm, 2.5*mm)); //F4 al to F4 2043P1
 153     pF4al->AddNode(pF4in    ,copy=1,new TGeoTranslation(   0*mm,0*mm,   0*mm)); //F4 whole F4 al 2043P1
 154     for(int i=1;i<=322;i++)
 155       pF4->AddNode(pColl,copy=i,new TGeoCombiTrans( 0*mm, -1296/2*mm+i*4*mm, -5*mm,pCollRot)); //F4 2043P1
 156 //radiators
 157   TGeoVolume *pRad      =gGeoManager->MakeBox( "Rad"      ,gGeoManager->GetMedium("HMPID_C6F14")    ,dx=1330*mm/2 ,dy= 413*mm/2  ,dz=  24*mm/2); // Rad 2011P1
 158   TGeoVolume *pRadFront =gGeoManager->MakeBox( "RadFront" ,gGeoManager->GetMedium("HMPID_Neoceram") ,dx=1330*mm/2 ,dy= 413*mm/2  ,dz=   4*mm/2);
 159   TGeoVolume *pRadWin   =gGeoManager->MakeBox( "RadWin"   ,gGeoManager->GetMedium("HMPID_SiO2")     ,dx=1330*mm/2 ,dy= 413*mm/2  ,dz=   5*mm/2);
 160   TGeoVolume *pRadLong  =gGeoManager->MakeBox( "RadLong"  ,gGeoManager->GetMedium("HMPID_Neoceram") ,dx=1330*mm/2 ,dy=   5*mm/2  ,dz=  15*mm/2);
 161   TGeoVolume *pRadShort =gGeoManager->MakeBox( "RadShort" ,gGeoManager->GetMedium("HMPID_Neoceram") ,dx=  10*mm/2 ,dy= 403*mm/2  ,dz=  15*mm/2);
 162   TGeoVolume *pRadSpacer=gGeoManager->MakeTube("RadSpacer",gGeoManager->GetMedium("HMPID_SiO2")     ,r1= 0        ,r2=10*mm/2  ,dz=  15*mm/2);
 163
 164   pRich->AddNode(pRad      ,copy=1,new TGeoTranslation(   0*mm,-434*mm,   -12*mm));
 165   pRich->AddNode(pRad      ,copy=2,new TGeoTranslation(   0*mm,   0*mm,   -12*mm));
 166   pRich->AddNode(pRad      ,copy=3,new TGeoTranslation(   0*mm,+434*mm,   -12*mm));
 167
 168   pRad ->AddNode(pRadFront ,copy=1,new TGeoTranslation(   0*mm,   0*mm, -10.0*mm));
 169   pRad ->AddNode(pRadWin   ,copy=1,new TGeoTranslation(   0*mm,   0*mm,   9.5*mm));
 170   pRad ->AddNode(pRadLong  ,copy=1,new TGeoTranslation(   0*mm,-204*mm,  -0.5*mm));
 171   pRad ->AddNode(pRadLong  ,copy=2,new TGeoTranslation(   0*mm,+204*mm,  -0.5*mm));
 172   pRad ->AddNode(pRadShort ,copy=1,new TGeoTranslation(-660*mm,   0*mm,  -0.5*mm));
 173   pRad ->AddNode(pRadShort ,copy=2,new TGeoTranslation(+660*mm,   0*mm,  -0.5*mm));
 174   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));
 175 //sandbox
 176   TGeoVolume *pSandBox  =gGeoManager->MakeBox( "RSandBox"  ,gGeoManager->GetMedium("Air")  ,dx=1419*mm/2 ,dy=1378*mm/2   ,dz=50.5*mm/2);  //2072P1
 177   TGeoVolume *pSandCover=gGeoManager->MakeBox( "RSandCover",gGeoManager->GetMedium("HMPID_Al")       ,dx=1419*mm/2 ,dy=1378*mm/2   ,dz= 0.5*mm/2);
 178   TGeoVolume *pSandComb =gGeoManager->MakeBox( "RSandComb" ,gGeoManager->GetMedium("HMPID_Rohacell") ,dx=1359*mm/2 ,dy=1318*mm/2   ,dz=49.5*mm/2);
 179
 180   pRich->AddNode(pSandBox,copy=1,new TGeoTranslation(  0*mm,0*mm, -73.75*mm)); //p.84 TDR
 181     pSandBox->AddNode(pSandComb  ,copy=1,new TGeoTranslation(  0*mm,0*mm,      0*mm)); //2072P1
 182     pSandBox->AddNode(pSandCover ,copy=1,new TGeoTranslation(  0*mm,0*mm,    +25*mm));
 183     pSandBox->AddNode(pSandCover ,copy=2,new TGeoTranslation(  0*mm,0*mm,    -25*mm));
 184 }//Rich()
 185 //__________________________________________________________________________________________________
 186 void Axis()
 187 {
 188 // Draw axises  on top of geometry
 189   Double_t X[6]={0,0,0,300,0,0};  Double_t Y[6]={0,0,0,0,300,0};  Double_t Z[6]={0,0,0,0,0,300};
 190   TPolyLine3D *pXaxis=new TPolyLine3D(2,X);pXaxis->SetLineColor(kRed);   pXaxis->Draw();
 191   TPolyLine3D *pYaxis=new TPolyLine3D(2,Y);pYaxis->SetLineColor(kGreen); pYaxis->Draw();
 192   TPolyLine3D *pZaxis=new TPolyLine3D(2,Z);pZaxis->SetLineColor(kBlue);  pZaxis->Draw();
 193 }
 194 //__________________________________________________________________________________________________