| 1 |
#include "PamVMCDetCas.h" |
| 2 |
|
| 3 |
ClassImp(PamVMCDetCas) |
| 4 |
|
| 5 |
void PamVMCDetCas::DefineGeometry(){ |
| 6 |
|
| 7 |
TGeoMaterial *m=0; |
| 8 |
TGeoMixture *x=0; |
| 9 |
TString s = "ALUMINIUM"; |
| 10 |
m = GetMat(s); |
| 11 |
if(!m){ |
| 12 |
m = new TGeoMaterial(s,26.98,13,2.7); |
| 13 |
m->SetUniqueID( 9); |
| 14 |
} |
| 15 |
|
| 16 |
|
| 17 |
s = "SCINT"; |
| 18 |
x = (TGeoMixture*)GetMat(s); |
| 19 |
if(!x){ |
| 20 |
x = new TGeoMixture(s,2, 1.03200 ); |
| 21 |
x->SetUniqueID( 17); |
| 22 |
x->DefineElement(0,12.01,6,0.922427); |
| 23 |
x->DefineElement(1,1.01,1,0.7757296E-01); |
| 24 |
} |
| 25 |
|
| 26 |
s = "N2GAS"; |
| 27 |
m = GetMat(s); |
| 28 |
if (!m){ |
| 29 |
m = new TGeoMaterial(s,14.01,7,0.1250000E-02); |
| 30 |
m->SetUniqueID( 19); |
| 31 |
} |
| 32 |
|
| 33 |
s = "PLAST"; |
| 34 |
x = (TGeoMixture*)GetMat(s); |
| 35 |
if(!x){ |
| 36 |
x = new TGeoMixture(s,2, 1.03200 ); |
| 37 |
x->SetUniqueID(28); |
| 38 |
x->DefineElement(0,12.01,6,0.9929774); |
| 39 |
x->DefineElement(1,1.01,1,0.7022570E-02); |
| 40 |
} |
| 41 |
|
| 42 |
TGeoMedium *n =0; |
| 43 |
|
| 44 |
s="ALUMINIUM"; |
| 45 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03); |
| 46 |
|
| 47 |
s = "N2 GAS"; |
| 48 |
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); |
| 49 |
|
| 50 |
s = "SCINTILLATOR"; |
| 51 |
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); |
| 52 |
|
| 53 |
s = "PLASTIC"; |
| 54 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,18,28,-1,1,40,1,100,0.5000000E-01,0.1000000E-01,0.5000000E-03); |
| 55 |
|
| 56 |
|
| 57 |
TGeoVolume *v = 0; |
| 58 |
|
| 59 |
|
| 60 |
s = "CASA"; |
| 61 |
v = (GetVol(s))?GetVol(s):gGeoManager->MakeBox("CASA",GetMed("N2 GAS"),21.45,0.6,20.55); |
| 62 |
s = "STRP"; |
| 63 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox("STRP",GetMed("ALUMINIUM"),21.45,0.6,20.55); |
| 64 |
s = "SID1"; |
| 65 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox("SID1",GetMed("SCINTILLATOR"),20.3,0.4,19.4); |
| 66 |
s = "PLA1"; |
| 67 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox("PLA1",GetMed("PLASTIC"),21.1,0.4,20.2); |
| 68 |
s = "STP1"; |
| 69 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STP1",GetMed("PLASTIC"),16.1,7.361412634,0.,0.4,0.,0.,0.,0.4,4.16,4.16,0.); |
| 70 |
s = "STP2"; |
| 71 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STP2",GetMed("PLASTIC"),2.4,20.55604522,0.,0.4,0.,0.,0.,0.4,1.8,1.8,0.); |
| 72 |
s = "STP3"; |
| 73 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STP3",GetMed("PLASTIC"),0.9,58.17255342,0.,0.4,1.8,1.8,0.,0.4,4.7,4.7,0.); |
| 74 |
s = "STA1"; |
| 75 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STA1",GetMed("ALUMINIUM"),16.5,7.355684758,0.,0.4,0.,0.,0.,0.4,4.26,4.26,0.); |
| 76 |
s = "STA2"; |
| 77 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STA2",GetMed("ALUMINIUM"),2.8,21.44773633,0.,0.4,0.,0.,0.,0.4,2.2,2.2,0.); |
| 78 |
s = "STA3"; |
| 79 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STA3",GetMed("ALUMINIUM"),0.9,58.17255342,0.,0.4,2.2,2.2,0.,0.4,5.1,5.1,0.); |
| 80 |
s = "STL1"; |
| 81 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STL1",GetMed("N2 GAS"),16.675,7.355377512,0.,0.6,0.,0.,0.,0.6,4.305,4.305,0.); |
| 82 |
s = "STL2"; |
| 83 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STL2",GetMed("N2 GAS"),2.8,21.44773633,0.,0.6,0.,0.,0.,0.6,2.2,2.2,0.); |
| 84 |
s = "STL3"; |
| 85 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrap("STL3",GetMed("N2 GAS"),1.075,58.18078646,0.,0.6,2.15,2.15,0.,0.6,5.615,5.615,0.); |
| 86 |
|
| 87 |
|
| 88 |
GetVol("CASA")->AddNode(GetVol("STRP"),1,gGeoIdentity); |
| 89 |
GetVol("STRP")->AddNodeOverlap(GetVol("PLA1"),1,gGeoIdentity); |
| 90 |
GetVol("PLA1")->AddNodeOverlap(GetVol("SID1"),1,gGeoIdentity); |
| 91 |
GetVol("SID1")->AddNode(GetVol("STP1"),1,new TGeoCombiTrans(18.22,0.,-3.3,GetRot("rot5"))); |
| 92 |
GetVol("SID1")->AddNode(GetVol("STP2"),1,new TGeoCombiTrans(19.4,0.,15.2,GetRot("rot4"))); |
| 93 |
GetVol("SID1")->AddNode(GetVol("STP3"),1,new TGeoCombiTrans(17.05,0.,18.5,GetRot("rot4"))); |
| 94 |
GetVol("PLA1")->AddNode(GetVol("STA1"),1,new TGeoCombiTrans(18.97,0.,-3.7,GetRot("rot5"))); |
| 95 |
GetVol("PLA1")->AddNode(GetVol("STA2"),1,new TGeoCombiTrans(20.,0.,15.6,GetRot("rot4"))); |
| 96 |
GetVol("PLA1")->AddNode(GetVol("STA3"),1,new TGeoCombiTrans(17.45,0.,19.3,GetRot("rot4"))); |
| 97 |
GetVol("STRP")->AddNode(GetVol("STL1"),1,new TGeoCombiTrans(19.2975,0.,-3.875,GetRot("rot5"))); |
| 98 |
GetVol("STRP")->AddNode(GetVol("STL2"),1,new TGeoCombiTrans(20.35,0.,15.6,GetRot("rot4"))); |
| 99 |
GetVol("STRP")->AddNode(GetVol("STL3"),1,new TGeoCombiTrans(17.5675,0.,19.475,GetRot("rot4"))); |
| 100 |
|
| 101 |
SetMotherProp(GetVol("CASA"),1,new TGeoTranslation(8.872,-17.215,47.174)); |
| 102 |
SetMotherProp(GetVol("CASA"),2,new TGeoCombiTrans(17.9285,6.305,47.174,GetRot("rot1"))); |
| 103 |
SetMotherProp(GetVol("CASA"),3,new TGeoCombiTrans(-5.5095,17.215,47.174,GetRot("rot4"))); |
| 104 |
SetMotherProp(GetVol("CASA"),4,new TGeoCombiTrans(-17.9285,-10.195,47.174,GetRot("rot6"))); |
| 105 |
|
| 106 |
}; |
| 107 |
|
| 108 |
|
| 109 |
void PamVMCDetCas::DefineCuts(){ |
| 110 |
|
| 111 |
TString s ="ALUMINIUM"; // default GPAMELA parameters |
| 112 |
if (GetMed(s) && !GetCC(s)) |
| 113 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 114 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 115 |
s ="N2 GAS"; |
| 116 |
if (GetMed(s) && !GetCC(s)) |
| 117 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 118 |
0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.)); |
| 119 |
s ="SCINTILLATOR";// default GPAMELA parameters |
| 120 |
if (GetMed(s) && !GetCC(s)) |
| 121 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 122 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 123 |
s = "PLASTIC";// default GPAMELA parameters |
| 124 |
if (GetMed(s) && !GetCC(s)) |
| 125 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 126 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 127 |
|
| 128 |
} |
| 129 |
|
Parent Directory
| 
Revision Log