| 1 |
#include <TMath.h> |
| 2 |
#include "PamVMCDetCard.h" |
| 3 |
|
| 4 |
ClassImp(PamVMCDetCard) |
| 5 |
|
| 6 |
using TMath::Tan; |
| 7 |
using TMath::ATan; |
| 8 |
using TMath::Sin; |
| 9 |
using TMath::Cos; |
| 10 |
using TMath::Pi; |
| 11 |
|
| 12 |
void PamVMCDetCard::DefineGeometry(){ |
| 13 |
|
| 14 |
TGeoMaterial *m=0; |
| 15 |
TGeoMixture *x=0; |
| 16 |
TString s = "ALUMINIUM"; |
| 17 |
m = GetMat(s); |
| 18 |
if(!m){ |
| 19 |
m = new TGeoMaterial(s,26.98,13,2.7); |
| 20 |
m->SetUniqueID( 9); |
| 21 |
} |
| 22 |
|
| 23 |
s = "IRON"; |
| 24 |
m = GetMat(s); |
| 25 |
if(!m){ |
| 26 |
m = new TGeoMaterial(s,55.85,26,7.87); |
| 27 |
m->SetUniqueID( 10); |
| 28 |
} |
| 29 |
|
| 30 |
s = "SCINT"; |
| 31 |
x = (TGeoMixture*)GetMat(s); |
| 32 |
if(!x){ |
| 33 |
x = new TGeoMixture(s,2, 1.03200 ); |
| 34 |
x->SetUniqueID( 17); |
| 35 |
x->DefineElement(0,12.01,6,0.922427); |
| 36 |
x->DefineElement(1,1.01,1,0.7757296E-01); |
| 37 |
} |
| 38 |
|
| 39 |
s = "N2GAS"; |
| 40 |
m = GetMat(s); |
| 41 |
if (!m){ |
| 42 |
m = new TGeoMaterial(s,14.01,7,0.1250000E-02); |
| 43 |
m->SetUniqueID( 19); |
| 44 |
} |
| 45 |
|
| 46 |
s = "PLAST"; |
| 47 |
x = (TGeoMixture*)GetMat(s); |
| 48 |
if(!x){ |
| 49 |
x = new TGeoMixture(s,2, 1.03200 ); |
| 50 |
x->SetUniqueID( 28); |
| 51 |
x->DefineElement(0,12.01,6,0.9929774); |
| 52 |
x->DefineElement(1,1.01,1,0.7022570E-02); |
| 53 |
} |
| 54 |
|
| 55 |
TGeoMedium *n =0; |
| 56 |
|
| 57 |
s="ALUMINIUM"; |
| 58 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03); |
| 59 |
|
| 60 |
s = "N2 GAS"; |
| 61 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,3,19,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.1000000E-02,0.5000000E-03); |
| 62 |
|
| 63 |
s = "SCINTILLATOR"; |
| 64 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,10,17,1,1,40,1,0.5000000E-01,0.5000000E-01,0.1000000E-01,0.5000000E-03); |
| 65 |
|
| 66 |
s = "PLASTIC"; |
| 67 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,18,28,-1,1,40,1,100,0.5000000E-01,0.1000000E-01,0.5000000E-03); |
| 68 |
|
| 69 |
s = "IRON"; |
| 70 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,27,10,-1,1,40,1,100,0.5000000E-01,0.2491200E-02,0.5000000E-03); |
| 71 |
|
| 72 |
|
| 73 |
TGeoVolume *v = 0; |
| 74 |
|
| 75 |
Double_t H1=6.6; |
| 76 |
Double_t H2=11.55; |
| 77 |
Double_t B1=20.95+3.6+5.8+H1/Tan(76./180.*Pi()); |
| 78 |
Double_t C1=H2/Tan(76./180.*Pi()); |
| 79 |
Double_t A1=37.18273-C1; |
| 80 |
Double_t D1=37.18273-20.95-3.6-5.8; |
| 81 |
Double_t ANG=ATan((B1/2.+D1-A1/2.)/18.15)/Pi()*180.; |
| 82 |
|
| 83 |
Double_t B2=17.45+3.6+5.8+H1/Tan(76./180.*Pi()); |
| 84 |
Double_t C2=H2/Tan(76./180.*Pi()); |
| 85 |
Double_t A2=31.98273-C2; |
| 86 |
Double_t D2=31.98273-17.45-3.6-5.8; |
| 87 |
Double_t ANG2=ATan((B2/2.+D2-A2/2.)/18.15)/Pi()*180.; |
| 88 |
|
| 89 |
|
| 90 |
Double_t BB=B1/2.+1.1*(1./Cos(ATan(D1/18.15))-D1/18.15+1./Tan(76./180.*Pi())+1./Sin(76./180.*Pi()))/2.; |
| 91 |
Double_t AA=A1/2.+1.1*(1./Cos(ATan(D1/18.15))+D1/18.15-1./Tan(76./180.*Pi())+1./Sin(76./180.*Pi()))/2.; |
| 92 |
Double_t DD=10.175*2.*D1/18.15; |
| 93 |
Double_t CAR1_2=-ATan((BB+DD-AA)/(2.*10.175))/Pi()*180.; |
| 94 |
Double_t CAR1_9=BB; |
| 95 |
Double_t C1D1_9=B1/2.; |
| 96 |
Double_t C2D1_9=B2/2.; |
| 97 |
|
| 98 |
s = "C1D1"; |
| 99 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("SCINTILLATOR"),9.075,-ANG,0.,0.4,A1/2.,A1/2.,0.,0.4,C1D1_9,C1D1_9,0.); |
| 100 |
|
| 101 |
|
| 102 |
|
| 103 |
s = "C2D1"; |
| 104 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("SCINTILLATOR"),9.075,-ANG2,0.,0.4,A2/2.,A2/2.,0.,0.4,C2D1_9,C2D1_9,0.); |
| 105 |
|
| 106 |
|
| 107 |
s = "CAR1"; |
| 108 |
v = (GetVol(s))?GetVol(s):gGeoManager->MakeTrap(s,GetMed("ALUMINIUM"),10.175,CAR1_2,0.,0.55,AA,AA,0.,0.55,CAR1_9,CAR1_9,0.); |
| 109 |
|
| 110 |
|
| 111 |
BB=B1/2.+0.9*(1./Cos(ATan(D1/18.15))-D1/18.15+1./Tan(76./180.*Pi())+1./Sin(76./180.*Pi()))/2.; |
| 112 |
AA=A1/2.+0.9*(1./Cos(ATan(D1/18.15))+D1/18.15-1./Tan(76./180.*Pi())+1./Sin(76./180.*Pi()))/2.; |
| 113 |
DD=9.975*2.*D1/18.15; |
| 114 |
Double_t CR1P_2=-ATan((BB+DD-AA)/(2.*9.975))/Pi()*180.; |
| 115 |
Double_t CR1P_9=BB; |
| 116 |
|
| 117 |
s = "CR1P"; |
| 118 |
v = (GetVol(s))?GetVol(s):gGeoManager->MakeTrap(s,GetMed("PLASTIC"),9.975,CR1P_2,0.,0.4,AA,AA,0.,0.4,CR1P_9,CR1P_9,0.); |
| 119 |
|
| 120 |
//===========// |
| 121 |
|
| 122 |
BB=B2/2.+1.1*(1./Cos(ATan(D2/18.15))-D2/18.15+1./Tan(76./180.*Pi())+1./Sin(76./180.*Pi()))/2.; |
| 123 |
AA=A2/2.+1.1*(1./Cos(ATan(D2/18.15))+D2/18.15-1./Tan(76./180.*Pi())+1./Sin(76./180.*Pi()))/2.; |
| 124 |
DD=10.175*2.*D2/18.15; |
| 125 |
Double_t CAR2_2=-ATan((BB+DD-AA)/(2.*10.175))/Pi()*180.; |
| 126 |
Double_t CAR2_9=BB; |
| 127 |
|
| 128 |
s = "CAR2"; |
| 129 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("ALUMINIUM"),10.175,CAR2_2,0.,0.55,AA,AA,0.,0.55,CAR2_9,CAR2_9,0.); |
| 130 |
|
| 131 |
//==========// |
| 132 |
|
| 133 |
BB=B2/2.+0.9*(1./Cos(ATan(D2/18.15))-D2/18.15+1./Tan(76./180.*Pi())+1./Sin(76./180.*Pi()))/2.; |
| 134 |
AA=A2/2.+0.9*(1./Cos(ATan(D2/18.15))+D2/18.15-1./Tan(76./180.*Pi())+1./Sin(76./180.*Pi()))/2.; |
| 135 |
DD=9.975*2.*D2/18.15; |
| 136 |
Double_t CR2P_2=-ATan((BB+DD-AA)/(2.*9.975))/Pi()*180.; |
| 137 |
Double_t CR2P_9=BB; |
| 138 |
|
| 139 |
s = "CR2P"; |
| 140 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("PLASTIC"),9.975,CR2P_2,0.,0.4,AA,AA,0.,0.4,CR2P_9,CR2P_9,0.); |
| 141 |
|
| 142 |
//==========// |
| 143 |
|
| 144 |
AA=4.8+1.8; |
| 145 |
BB=AA/Tan(76./180.*Pi())+AA*3.6/4.8; |
| 146 |
Double_t VSN1_2=ATan((BB/2.-AA/Tan(76./180.*Pi()))/AA)/Pi()*180.; |
| 147 |
Double_t VSN1_9=BB/2.; |
| 148 |
|
| 149 |
s = "VSN1"; |
| 150 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("N2 GAS"),3.3,VSN1_2,0.,0.4,0.,0.,0.,0.4,VSN1_9,VSN1_9,0.); |
| 151 |
|
| 152 |
//==========// |
| 153 |
|
| 154 |
AA=1.8; |
| 155 |
BB=AA*(5.8/1.8-3.6/4.8); |
| 156 |
Double_t VSN2_1=0.9; |
| 157 |
Double_t VSN2_2=ATan((BB/2.+AA*3.6/4.8)/AA)/Pi()*180.; |
| 158 |
Double_t VSN2_9=2.225; |
| 159 |
|
| 160 |
s = "VSN2"; |
| 161 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("PLASTIC"),VSN2_1,VSN2_2,0.,0.4,0.,0.,0.,0.4,VSN2_9,VSN2_9,0.); |
| 162 |
|
| 163 |
//=========// |
| 164 |
|
| 165 |
AA=1.8+0.9-Cos(ATan(3.6/4.8))*0.9/Sin(ATan(5.8/1.8)-ATan(3.6/4.8)); |
| 166 |
BB=AA*(5.8/1.8-3.6/4.8); |
| 167 |
Double_t VPN2_1=AA/2.; |
| 168 |
Double_t VPN2_2=ATan((BB/2.+AA*3.6/4.8)/AA)/Pi()*180.; |
| 169 |
Double_t VPN2_9=BB/2.; |
| 170 |
|
| 171 |
|
| 172 |
s = "VPN2"; |
| 173 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("ALUMINIUM"),VPN2_1,VPN2_2,0.,/*0.55*/ 0.4,0.,0.,0.,/*0.55*/0.4,VPN2_9,VPN2_9,0.); |
| 174 |
|
| 175 |
//=========// |
| 176 |
|
| 177 |
AA=(0.9/Sin(76./180.*Pi()))/(1./Tan(76./180.*Pi())+3.6/4.8)+4.8+1.8+0.9; |
| 178 |
BB=AA/Tan(76./180.*Pi())+AA*3.6/4.8; |
| 179 |
Double_t VPN1_1=AA/2.; |
| 180 |
Double_t VPN1_9=BB/2.; |
| 181 |
|
| 182 |
s = "VPN1"; |
| 183 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("N2 GAS"),VPN1_1,VSN1_2,0.,0.4,0.,0.,0.,0.4,VPN1_9,VPN1_9,0.); |
| 184 |
|
| 185 |
//=========// |
| 186 |
|
| 187 |
AA=(1.1/Sin(76./180.*Pi()))/(1./Tan(76./180.*Pi())+3.6/4.8)+4.8+1.8+1.1; |
| 188 |
Double_t VAN1_1=AA/2.; |
| 189 |
BB=AA/Tan(76./180.*Pi())+AA*3.6/4.8; |
| 190 |
Double_t VAN1_9=BB/2.; |
| 191 |
|
| 192 |
s = "VAN1"; |
| 193 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("N2 GAS"),VAN1_1,VSN1_2,0.,0.55,0.,0.,0.,0.55,VAN1_9,VAN1_9,0.); |
| 194 |
|
| 195 |
//=========// |
| 196 |
|
| 197 |
AA=1.8+1.1-Cos(ATan(3.6/4.8))*1.1/Sin(ATan(5.8/1.8)-ATan(3.6/4.8)); |
| 198 |
BB=AA*(5.8/1.8-3.6/4.8); |
| 199 |
Double_t VAN2_1=AA/2.; |
| 200 |
Double_t VAN2_2=ATan((BB/2.+AA*3.6/4.8)/AA)/Pi()*180.; |
| 201 |
Double_t VAN2_9=BB/2.; |
| 202 |
|
| 203 |
s = "VAN2"; |
| 204 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("N2 GAS"),VAN2_1,VAN2_2,0.,0.55,0.,0.,0.,0.55,VAN2_9,VAN2_9,0.); |
| 205 |
|
| 206 |
//=========// |
| 207 |
|
| 208 |
s = "CARD"; |
| 209 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),27.,22.8,12.35); |
| 210 |
s = "LEGB"; |
| 211 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),2.15,3.575,11.85); |
| 212 |
s = "LEG1"; |
| 213 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("N2 GAS"),0.65,0.,0.,0.15,6.85,6.85,0.,0.15,10.85,10.85,0.); |
| 214 |
s = "LEG2"; |
| 215 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),0.65,1.65,11.35); |
| 216 |
s = "LEG3"; |
| 217 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),0.1,1.3,0.325); |
| 218 |
s = "LEG4"; |
| 219 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),0.65,1.625,11.85); |
| 220 |
s = "LEG5"; |
| 221 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),1.125,2.575,11.85); |
| 222 |
s = "LEG6"; |
| 223 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),0.1,0.55,11.2); |
| 224 |
s = "LEG7"; |
| 225 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),1.1,0.75,3.25); |
| 226 |
s = "LEG8"; |
| 227 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),1.1,0.75,4.35); |
| 228 |
s = "LEG9"; |
| 229 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),1.1,0.75,3.25); |
| 230 |
s = "TPLT"; |
| 231 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),27.,22.8,0.5); |
| 232 |
s = "TOPC"; |
| 233 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),20.4,16.,0.5); |
| 234 |
s = "TH11"; |
| 235 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),2.8,3.2,0.5); |
| 236 |
s = "TH12"; |
| 237 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),2.8,5.3,0.5); |
| 238 |
s = "TH21"; |
| 239 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),1.7,3.,0.5); |
| 240 |
s = "TH22"; |
| 241 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),4.9,3.,0.5); |
| 242 |
s = "TOPT"; |
| 243 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap(s,GetMed("N2 GAS"),0.5,0.,0.,2.25,2.25,0.,-26.565,2.25,2.25,0.,-26.565); |
| 244 |
s = "TOHO"; |
| 245 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTube(s,GetMed("N2 GAS"),0.,0.7,0.5); |
| 246 |
s = "MGSC"; |
| 247 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("IRON"),16.,15.4,0.12456); |
| 248 |
s = "MGSH"; |
| 249 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),9.8,8.,0.12456); |
| 250 |
|
| 251 |
|
| 252 |
Double_t X,Y,Z; |
| 253 |
Double_t X1,Y1,Z1,X2,Y2,Z2; |
| 254 |
Double_t ZH=49.229-(13.05+10.639+2.97+2.*22.615+1.6)-1.6; |
| 255 |
|
| 256 |
Y=-16.014+19.45/2.*Sin(17./180.*Pi()); |
| 257 |
Z=47.661+ZH-12.35-19.45*Cos(17./180.*Pi())/2.; |
| 258 |
X=((35.6/2.)-Tan(ANG/180.*Pi())*19.45/2.)-18.67; |
| 259 |
X1=X; Y1=Y; Z1=Z; |
| 260 |
GetVol("CARD")->AddNode(GetVol("CAR1"),1,new TGeoCombiTrans(X1,Y1,Z1,GetRot("rot11"))); |
| 261 |
Z=10.175-VAN1_1; |
| 262 |
X=-(CAR1_9-10.175*Tan(CAR1_2/180.*Pi())-(VAN1_9-VAN1_1*Tan(VSN1_2/180.*Pi()))); |
| 263 |
GetVol("CAR1")->AddNode(GetVol("VAN1"),1,new TGeoTranslation(X,0.,Z));//many |
| 264 |
Z=10.175-VAN2_1; |
| 265 |
X=-(CAR1_9-10.175*Tan(CAR1_2/180.*Pi())-(2.*VAN1_9+VAN2_9-VAN2_1*Tan(VAN2_2/180.*Pi()))); |
| 266 |
GetVol("CAR1")->AddNode(GetVol("VAN2"),1,new TGeoTranslation(X,0.,Z));//many |
| 267 |
X=(1.1-0.9)*(1./Cos(ATan(D1/18.15))-1./Sin(76./180.*Pi()))/2.; |
| 268 |
GetVol("CAR1")->AddNodeOverlap(GetVol("CR1P"),1,new TGeoTranslation(X,0.,0.));//only |
| 269 |
Z=9.975-VPN1_1; |
| 270 |
X=-(CR1P_9-9.975*Tan(CR1P_2/180.*Pi())-(VPN1_9-VPN1_1*Tan(VSN1_2/180.*Pi()))); |
| 271 |
GetVol("CR1P")->AddNode(GetVol("VPN1"),1,new TGeoTranslation(X,0.,Z)); //many |
| 272 |
Z=9.975-VPN2_1; |
| 273 |
X=-(CR1P_9-9.975*Tan(CR1P_2/180.*Pi())-(2.*VPN1_9+VPN2_9-VPN2_1*Tan(VPN2_2/180.*Pi()))); |
| 274 |
GetVol("CR1P")->AddNode(GetVol("VPN2"),1,new TGeoTranslation(X,0.,Z));//many |
| 275 |
X=0.9*(1./Cos(ATan(D1/18.15))-1./Sin(76./180.*Pi()))/2.; |
| 276 |
GetVol("CR1P")->AddNodeOverlap(GetVol("C1D1"),1,new TGeoTranslation(X,0.,0.));//only |
| 277 |
X=-(C1D1_9+9.075*Tan(ANG/180.*Pi())-(VSN1_9-3.3*Tan(VSN1_2/180.*Pi()))); |
| 278 |
GetVol("C1D1")->AddNode(GetVol("VSN1"),1,new TGeoTranslation(X,0.,5.775)); //only |
| 279 |
X=-(C1D1_9+9.075*Tan(ANG/180.*Pi())-(2.*VSN1_9+2.225-0.9*Tan(VSN2_2/180.*Pi()))); |
| 280 |
GetVol("C1D1")->AddNode(GetVol("VSN2"),1,new TGeoTranslation(X,0.,8.175)); //only |
| 281 |
|
| 282 |
|
| 283 |
GetVol("CARD")->AddNode(GetVol("CAR1"),2,new TGeoCombiTrans(-X1,-Y1,Z1,GetRot("rot12"))); |
| 284 |
|
| 285 |
X=19.123-19.45/2.*Sin(20./180.*Pi()); |
| 286 |
Z=47.358+ZH-12.35-19.45*Cos(20./180.*Pi())/2.; |
| 287 |
Y=(30.4/2.-Tan(ANG2/180.*Pi())*19.45/2.)-16.8; |
| 288 |
X2=X; Y2=Y; Z2=Z; |
| 289 |
|
| 290 |
GetVol("CARD")->AddNode(GetVol("CAR2"),1,new TGeoCombiTrans(X2,Y2,Z2,GetRot("rot14"))); |
| 291 |
|
| 292 |
Z=10.175-VAN1_1; |
| 293 |
X=-(CAR2_9-10.175*Tan(CAR2_2/180.*Pi())-(VAN1_9-VAN1_1*Tan(VSN1_2/180.*Pi()))); |
| 294 |
GetVol("CAR2")->AddNode(GetVol("VAN1"),2,new TGeoTranslation(X,0.,Z)); //many |
| 295 |
Z=10.175-VAN2_1; |
| 296 |
X=-(CAR2_9-10.175*Tan(CAR2_2/180.*Pi())-(2.*VAN1_9+VAN2_9-VAN2_1*Tan(VAN2_2/180.*Pi()))); |
| 297 |
GetVol("CAR2")->AddNode(GetVol("VAN2"),2,new TGeoTranslation(X,0.,Z)); //many |
| 298 |
X=(1.1-0.9)*(1./Cos(ATan(D2/18.15))-1./Sin(76./180.*Pi()))/2.; |
| 299 |
GetVol("CAR2")->AddNodeOverlap(GetVol("CR2P"),1,new TGeoTranslation(X,0.,0.)); // only |
| 300 |
Z=9.975-VPN1_1; |
| 301 |
X=-(CR2P_9+9.975*Tan(ANG2/180.*Pi())-(VPN1_9-VPN1_1*Tan(VSN1_2/180.*Pi()))); |
| 302 |
GetVol("CR2P")->AddNode(GetVol("VPN1"),2,new TGeoTranslation(X,0.,Z)); //many |
| 303 |
Z=9.975-VPN2_1; |
| 304 |
X=-(CR2P_9-9.975*Tan(CR2P_2/180.*Pi())-(2.*VPN1_9+VPN2_9-VPN2_1*Tan(VPN2_2/180.*Pi()))); |
| 305 |
GetVol("CR2P")->AddNode(GetVol("VPN2"),2,new TGeoTranslation(X,0.,Z));//many |
| 306 |
X=0.9*(1./Cos(ATan(D2/18.15))-1./Sin(76./180.*Pi()))/2.; |
| 307 |
GetVol("CR2P")->AddNodeOverlap(GetVol("C2D1"),1,new TGeoTranslation(X,0.,0.));//only |
| 308 |
Z=9.075-3.3; |
| 309 |
X=-(C2D1_9+9.075*Tan(ANG2/180.*Pi())-(VSN1_9-3.3*Tan(VSN1_2/180.*Pi()))); |
| 310 |
GetVol("C2D1")->AddNode(GetVol("VSN1"),2,new TGeoTranslation(X,0.,Z)); //only |
| 311 |
Z=9.075-0.9; |
| 312 |
X=-(C2D1_9+9.075*Tan(ANG2/180.*Pi())-(2.*VSN1_9+VSN2_9-VSN2_1*Tan(VSN2_2/180.*Pi()))); |
| 313 |
GetVol("C2D1")->AddNode(GetVol("VSN2"),2,new TGeoTranslation(X,0.,Z)); //only |
| 314 |
|
| 315 |
GetVol("CARD")->AddNode(GetVol("CAR2"),2,new TGeoCombiTrans(-X2,-Y2,Z2,GetRot("rot13"))); |
| 316 |
|
| 317 |
//==Checked==// |
| 318 |
|
| 319 |
GetVol("CARD")->AddNode(GetVol("LEGB"),1,new TGeoTranslation(-22.9,-18.575,-0.5)); |
| 320 |
GetVol("LEGB")->AddNode(GetVol("LEG1"),1,new TGeoCombiTrans(-1.5,-3.425,0.,GetRot("rot7"))); |
| 321 |
GetVol("LEGB")->AddNode(GetVol("LEG1"),2,new TGeoCombiTrans(-1.5,0.175,0.,GetRot("rot7"))); |
| 322 |
GetVol("LEGB")->AddNode(GetVol("LEG2"),1,new TGeoTranslation(-1.5,-1.625,0.)); |
| 323 |
GetVol("LEGB")->AddNode(GetVol("LEG3"),1,new TGeoTranslation(-0.55,2.275,11.525)); |
| 324 |
GetVol("LEGB")->AddNode(GetVol("LEG3"),2,new TGeoTranslation(-0.55,2.275,-11.525)); |
| 325 |
GetVol("LEGB")->AddNode(GetVol("LEG4"),1,new TGeoTranslation(-1.5,1.95,0.)); |
| 326 |
GetVol("LEGB")->AddNode(GetVol("LEG5"),1,new TGeoTranslation(1.025,1.,0.)); |
| 327 |
GetVol("LEGB")->AddNode(GetVol("LEG6"),1,new TGeoTranslation(-0.75,3.025,0.)); |
| 328 |
GetVol("LEGB")->AddNode(GetVol("LEG7"),1,new TGeoTranslation(0.65,-2.825,8.)); |
| 329 |
GetVol("LEGB")->AddNode(GetVol("LEG8"),1,new TGeoTranslation(0.65,-2.825,0.)); |
| 330 |
GetVol("LEGB")->AddNode(GetVol("LEG9"),1,new TGeoTranslation(0.65,-2.825,-8.)); |
| 331 |
GetVol("CARD")->AddNode(GetVol("LEGB"),2,new TGeoCombiTrans(-22.9,18.575,-0.5,GetRot("rot10"))); |
| 332 |
GetVol("CARD")->AddNode(GetVol("LEGB"),3,new TGeoCombiTrans(22.9,18.575,-0.5,GetRot("rot4"))); |
| 333 |
GetVol("CARD")->AddNode(GetVol("LEGB"),4,new TGeoCombiTrans(22.9,-18.575,-0.5,GetRot("rot9"))); |
| 334 |
GetVol("CARD")->AddNode(GetVol("TPLT"),1,new TGeoTranslation(0.,0.,11.85)); |
| 335 |
GetVol("TPLT")->AddNode(GetVol("TOPC"),1,gGeoIdentity); |
| 336 |
GetVol("TPLT")->AddNode(GetVol("TH11"),1,new TGeoTranslation(-23.6,14.4,0.)); |
| 337 |
GetVol("TPLT")->AddNode(GetVol("TH11"),2,new TGeoTranslation(-23.6,-14.4,0.)); |
| 338 |
GetVol("TPLT")->AddNode(GetVol("TH11"),3,new TGeoTranslation(23.6,14.4,0.)); |
| 339 |
GetVol("TPLT")->AddNode(GetVol("TH11"),4,new TGeoTranslation(23.6,-14.4,0.)); |
| 340 |
GetVol("TPLT")->AddNode(GetVol("TH12"),1,new TGeoTranslation(-23.6,5.5,0.)); |
| 341 |
GetVol("TPLT")->AddNode(GetVol("TH12"),2,new TGeoTranslation(-23.6,-5.5,0.)); |
| 342 |
GetVol("TPLT")->AddNode(GetVol("TH12"),3,new TGeoTranslation(23.6,5.5,0.)); |
| 343 |
GetVol("TPLT")->AddNode(GetVol("TH12"),4,new TGeoTranslation(23.6,-5.5,0.)); |
| 344 |
GetVol("TPLT")->AddNode(GetVol("TH21"),1,new TGeoTranslation(-17.2,19.4,0.)); |
| 345 |
GetVol("TPLT")->AddNode(GetVol("TH21"),2,new TGeoTranslation(-17.2,-19.4,0.)); |
| 346 |
GetVol("TPLT")->AddNode(GetVol("TH21"),3,new TGeoTranslation(17.2,19.4,0.)); |
| 347 |
GetVol("TPLT")->AddNode(GetVol("TH21"),4,new TGeoTranslation(17.2,-19.4,0.)); |
| 348 |
GetVol("TPLT")->AddNode(GetVol("TH22"),1,new TGeoTranslation(-10.2,19.4,0.)); |
| 349 |
GetVol("TPLT")->AddNode(GetVol("TH22"),2,new TGeoTranslation(-10.2,-19.4,0.)); |
| 350 |
GetVol("TPLT")->AddNode(GetVol("TH22"),3,new TGeoTranslation(10.2,19.4,0.)); |
| 351 |
GetVol("TPLT")->AddNode(GetVol("TH22"),4,new TGeoTranslation(10.2,-19.4,0.)); |
| 352 |
GetVol("TPLT")->AddNode(GetVol("TH22"),5,new TGeoTranslation(0.,19.4,0.)); |
| 353 |
GetVol("TPLT")->AddNode(GetVol("TH22"),6,new TGeoTranslation(0.,-19.4,0.)); |
| 354 |
GetVol("TPLT")->AddNode(GetVol("TOPT"),1,new TGeoTranslation(-25.875,-20.55,0.)); |
| 355 |
GetVol("TPLT")->AddNode(GetVol("TOPT"),2,new TGeoCombiTrans(-25.875,20.55,0.,GetRot("rot10"))); |
| 356 |
GetVol("TPLT")->AddNode(GetVol("TOPT"),3,new TGeoCombiTrans(25.875,20.55,0.,GetRot("rot4"))); |
| 357 |
GetVol("TPLT")->AddNode(GetVol("TOPT"),4,new TGeoCombiTrans(25.875,-20.55,0.,GetRot("rot9"))); |
| 358 |
GetVol("TPLT")->AddNode(GetVol("TOHO"),1,new TGeoTranslation(20.1,18.,0.)); |
| 359 |
GetVol("TPLT")->AddNode(GetVol("TOHO"),2,new TGeoTranslation(-20.1,18.,0.)); |
| 360 |
GetVol("TPLT")->AddNode(GetVol("TOHO"),3,new TGeoTranslation(-20.1,-18.,0.)); |
| 361 |
GetVol("TPLT")->AddNode(GetVol("TOHO"),4,new TGeoTranslation(20.1,-18.,0.)); |
| 362 |
GetVol("CARD")->AddNode(GetVol("MGSC"),1,new TGeoTranslation(0.,0.,-12.22544)); |
| 363 |
GetVol("MGSC")->AddNode(GetVol("MGSH"),1,gGeoIdentity); |
| 364 |
|
| 365 |
|
| 366 |
SetMotherProp(GetVol("CARD"),1,new TGeoTranslation(0.,0.,87.439)); |
| 367 |
|
| 368 |
}; |
| 369 |
|
| 370 |
void PamVMCDetCard::DefineCuts(){ |
| 371 |
|
| 372 |
TString s ="ALUMINIUM"; // default GPAMELA parameters |
| 373 |
if (GetMed(s) && !GetCC(s)) |
| 374 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 375 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 376 |
s ="N2 GAS"; |
| 377 |
if (GetMed(s) && !GetCC(s)) |
| 378 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 379 |
0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.)); |
| 380 |
s ="SCINTILLATOR";// default GPAMELA parameters |
| 381 |
if (GetMed(s) && !GetCC(s)) |
| 382 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 383 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 384 |
s = "PLASTIC";// default GPAMELA parameters |
| 385 |
if (GetMed(s) && !GetCC(s)) |
| 386 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 387 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 388 |
s = "IRON"; //default GPAMELA parameters |
| 389 |
if (GetMed(s) && !GetCC(s)) |
| 390 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 391 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 392 |
} |
Parent Directory
| 
Revision Log