| 1 |
#include "PamVMCDetCat.h" |
| 2 |
|
| 3 |
ClassImp(PamVMCDetCat) |
| 4 |
void PamVMCDetCat::DefineGeometry(){ |
| 5 |
|
| 6 |
TGeoMaterial *m=0; |
| 7 |
TGeoMixture *x=0; |
| 8 |
TString s = "ALUMINIUM"; |
| 9 |
m = GetMat(s); |
| 10 |
if(!m){ |
| 11 |
m = new TGeoMaterial(s,26.98,13,2.7); |
| 12 |
m->SetUniqueID( 9); |
| 13 |
} |
| 14 |
|
| 15 |
|
| 16 |
s = "SCINT"; |
| 17 |
x = (TGeoMixture*)GetMat(s); |
| 18 |
if(!x){ |
| 19 |
x = new TGeoMixture(s,2, 1.03200 ); |
| 20 |
x->SetUniqueID( 17); |
| 21 |
x->DefineElement(0,12.01,6,0.922427); |
| 22 |
x->DefineElement(1,1.01,1,0.7757296E-01); |
| 23 |
} |
| 24 |
|
| 25 |
s = "N2GAS"; |
| 26 |
m = GetMat(s); |
| 27 |
if (!m){ |
| 28 |
m = new TGeoMaterial(s,14.01,7,0.1250000E-02); |
| 29 |
m->SetUniqueID( 19); |
| 30 |
} |
| 31 |
|
| 32 |
s = "PLAST"; |
| 33 |
x = (TGeoMixture*)GetMat(s); |
| 34 |
if(!x){ |
| 35 |
x = new TGeoMixture(s,2, 1.03200 ); |
| 36 |
x->SetUniqueID( 28); |
| 37 |
x->DefineElement(0,12.01,6,0.9929774); |
| 38 |
x->DefineElement(1,1.01,1,0.7022570E-02); |
| 39 |
} |
| 40 |
|
| 41 |
s = "TIALVA"; |
| 42 |
x = (TGeoMixture*)GetMat(s); |
| 43 |
if(!x){ |
| 44 |
x = new TGeoMixture(s,3, 4.49234 ); |
| 45 |
x->SetUniqueID( 32); |
| 46 |
x->DefineElement(0,47.867,22,0.9); |
| 47 |
x->DefineElement(1,26.9815,13,0.6000000E-01); |
| 48 |
x->DefineElement(2,50.9415,23,0.4000000E-01); |
| 49 |
} |
| 50 |
|
| 51 |
TGeoMedium *n =0; |
| 52 |
|
| 53 |
s="ALUMINIUM"; |
| 54 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03); |
| 55 |
|
| 56 |
s = "N2 GAS"; |
| 57 |
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); |
| 58 |
|
| 59 |
s = "SCINTILLATOR"; |
| 60 |
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); |
| 61 |
|
| 62 |
s = "PLASTIC"; |
| 63 |
n = (GetMed(s))?GetMed(s): new TGeoMedium(s,18,28,-1,1,40,1,100,0.5000000E-01,0.1000000E-01,0.5000000E-03); |
| 64 |
|
| 65 |
s = "TITANIUM Mixture"; |
| 66 |
n = (GetMed(s))?GetMed(s):new TGeoMedium(s,24,32,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.3200000E-01,0.5000000E-03); |
| 67 |
|
| 68 |
|
| 69 |
TGeoVolume *v = 0; |
| 70 |
|
| 71 |
TPLAfirst = kFALSE; |
| 72 |
s = "TPLA"; |
| 73 |
v = GetVol(s); |
| 74 |
if(!v) { |
| 75 |
v = gGeoManager->MakeBox("TPLA",GetMed("N2 GAS"),28.,27.165,1.6); |
| 76 |
TPLAfirst = kTRUE; |
| 77 |
} |
| 78 |
|
| 79 |
s = "TPTU"; |
| 80 |
v = (GetVol(s))?GetVol(s):gGeoManager->MakeBox(s,GetMed("TITANIUM Mixture"),2.55,1.95,0.6); |
| 81 |
s = "TPTM"; |
| 82 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("TITANIUM Mixture"),3.45,3.3,0.75); |
| 83 |
s = "TPTL"; |
| 84 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("TITANIUM Mixture"),21.9,21.1,0.25); |
| 85 |
s = "TPCV"; |
| 86 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),8.5,7.,0.25); |
| 87 |
s = "CATA"; |
| 88 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),28.,0.6,27.165); |
| 89 |
s = "ATBB"; |
| 90 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),28.,0.6,27.165); |
| 91 |
s = "PLA2"; |
| 92 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("PLASTIC"),27.65,0.4,26.815); |
| 93 |
s = "TOP1"; |
| 94 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SCINTILLATOR"),26.85,0.4,26.015); |
| 95 |
s = "PTBM"; |
| 96 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("PLASTIC"),10.11,0.4,8.44); |
| 97 |
s = "ATBM"; |
| 98 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),9.41,0.4,7.74); |
| 99 |
s = "LTBM"; |
| 100 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),9.11,0.6,7.44); |
| 101 |
s = "PTBC"; |
| 102 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("PLASTIC"),4.975,0.4,5.5575); |
| 103 |
s = "ATBC"; |
| 104 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),4.975,0.4,5.5575); |
| 105 |
s = "LTBC"; |
| 106 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),4.975,0.6,5.5575); |
| 107 |
s = "PTRT"; |
| 108 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("PLASTIC"),0.,16.9,0.4,2.405); |
| 109 |
s = "ATRT"; |
| 110 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("ALUMINIUM"),0.,17.7,0.4,2.405); |
| 111 |
s = "LTRT"; |
| 112 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("N2 GAS"),0.,18.05,0.6,2.405); |
| 113 |
s = "PTRB"; |
| 114 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("PLASTIC"),16.9,0.,0.4,2.405); |
| 115 |
s = "ATRB"; |
| 116 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("ALUMINIUM"),17.7,0.,0.4,2.405); |
| 117 |
s = "LTRB"; |
| 118 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("N2 GAS"),18.05,0.,0.6,2.405); |
| 119 |
s = "PTRL"; |
| 120 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("PLASTIC"),0.,14.9,0.4,2.375); |
| 121 |
s = "ATRL"; |
| 122 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("ALUMINIUM"),0.,15.7,0.4,2.375); |
| 123 |
s = "LTRL"; |
| 124 |
v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("N2 GAS"),0.,16.05,0.6,2.375); |
| 125 |
|
| 126 |
GetVol("TPTL")->AddNode(GetVol("TPCV"),1,gGeoIdentity); |
| 127 |
GetVol("CATA")->AddNode(GetVol("ATBB"),1,gGeoIdentity); //only |
| 128 |
GetVol("ATBB")->AddNodeOverlap(GetVol("PLA2"),1,gGeoIdentity); //only |
| 129 |
GetVol("PLA2")->AddNodeOverlap(GetVol("TOP1"),1,gGeoIdentity); //only |
| 130 |
GetVol("TOP1")->AddNode(GetVol("PTBM") ,1,gGeoIdentity);//only |
| 131 |
GetVol("TOP1")->AddNode(GetVol("PTBC"),1,new TGeoTranslation(-21.875,0.,20.4575)); |
| 132 |
GetVol("TOP1")->AddNode(GetVol("PTBC"),2,new TGeoTranslation(21.875,0.,20.4575)); |
| 133 |
GetVol("TOP1")->AddNode(GetVol("PTBC"),3,new TGeoTranslation(-21.875,0.,-20.4575)); |
| 134 |
GetVol("TOP1")->AddNode(GetVol("PTBC"),4,new TGeoTranslation(21.875,0.,-20.4575)); |
| 135 |
GetVol("TOP1")->AddNode(GetVol("PTRT"),1,new TGeoTranslation(0.,0.,23.61)); |
| 136 |
GetVol("TOP1")->AddNode(GetVol("PTRB"),1,new TGeoTranslation(0.,0.,-23.61)); |
| 137 |
GetVol("TOP1")->AddNode(GetVol("PTRL"),1,new TGeoCombiTrans(-24.475,0.,0.,GetRot("rot7"))); |
| 138 |
GetVol("TOP1")->AddNode(GetVol("PTRL"),2,new TGeoCombiTrans(24.475,0.,0.,GetRot("rot8"))); |
| 139 |
GetVol("PLA2")->AddNode(GetVol("ATBM"),1,gGeoIdentity); |
| 140 |
GetVol("PLA2")->AddNode(GetVol("ATBC"),1,new TGeoTranslation(-22.675,0.,21.2575)); |
| 141 |
GetVol("PLA2")->AddNode(GetVol("ATBC"),2,new TGeoTranslation(22.675,0.,21.2575)); |
| 142 |
GetVol("PLA2")->AddNode(GetVol("ATBC"),3,new TGeoTranslation(-22.675,0.,-21.2575)); |
| 143 |
GetVol("PLA2")->AddNode(GetVol("ATBC"),4,new TGeoTranslation(22.675,0.,-21.2575)); |
| 144 |
GetVol("PLA2")->AddNode(GetVol("ATRT"),1,new TGeoTranslation(0.,0.,24.41)); |
| 145 |
GetVol("PLA2")->AddNode(GetVol("ATRB"),1,new TGeoTranslation(0.,0.,-24.41)); |
| 146 |
GetVol("PLA2")->AddNode(GetVol("ATRL"),1,new TGeoCombiTrans(-25.275,0.,0.,GetRot("rot7"))); |
| 147 |
GetVol("PLA2")->AddNode(GetVol("ATRL"),2,new TGeoCombiTrans(25.275,0.,0.,GetRot("rot8"))); |
| 148 |
GetVol("ATBB")->AddNode(GetVol("LTBM"),1,gGeoIdentity); |
| 149 |
GetVol("ATBB")->AddNode(GetVol("LTBC"),1,new TGeoTranslation(-23.025,0.,21.6075)); |
| 150 |
GetVol("ATBB")->AddNode(GetVol("LTBC"),2,new TGeoTranslation(23.025,0.,21.6075)); |
| 151 |
GetVol("ATBB")->AddNode(GetVol("LTBC"),3,new TGeoTranslation(-23.025,0.,-21.6075)); |
| 152 |
GetVol("ATBB")->AddNode(GetVol("LTBC"),4,new TGeoTranslation(23.025,0.,-21.6075)); |
| 153 |
GetVol("ATBB")->AddNode(GetVol("LTRT"),1,new TGeoTranslation(0.,0.,24.76)); |
| 154 |
GetVol("ATBB")->AddNode(GetVol("LTRB"),1,new TGeoTranslation(0.,0.,-24.76)); |
| 155 |
GetVol("ATBB")->AddNode(GetVol("LTRL"),1,new TGeoCombiTrans(-25.625,0.,0.,GetRot("rot7"))); |
| 156 |
GetVol("ATBB")->AddNode(GetVol("LTRL"),2,new TGeoCombiTrans(25.625,0.,0.,GetRot("rot8"))); |
| 157 |
GetVol("TPLA")->AddNodeOverlap(GetVol("TPTU"),1,new TGeoTranslation(-20.65,-20.15,1.)); |
| 158 |
GetVol("TPLA")->AddNodeOverlap(GetVol("TPTU"),2,new TGeoTranslation(-20.65,20.15,1.)); |
| 159 |
GetVol("TPLA")->AddNodeOverlap(GetVol("TPTU"),3,new TGeoTranslation(20.65,-20.15,1.)); |
| 160 |
GetVol("TPLA")->AddNodeOverlap(GetVol("TPTU"),4,new TGeoTranslation(20.65,20.15,1.)); |
| 161 |
|
| 162 |
|
| 163 |
if(TPLAfirst) SetMotherProp(GetVol("TPLA"),1,new TGeoTranslation(0.,0.,73.489)); |
| 164 |
|
| 165 |
GetVol("TPLA")->AddNode(GetVol("TPTL"),1,new TGeoTranslation(0.,0.,-1.35)); |
| 166 |
GetVol("TPLA")->AddNode(GetVol("TPTM"),1,new TGeoTranslation(-19.75,-18.8,-0.35)); |
| 167 |
GetVol("TPLA")->AddNode(GetVol("TPTM"),2,new TGeoTranslation(-19.75,18.8,-0.35)); |
| 168 |
GetVol("TPLA")->AddNode(GetVol("TPTM"),3,new TGeoTranslation(19.75,-18.8,-0.35)); |
| 169 |
GetVol("TPLA")->AddNode(GetVol("TPTM"),4,new TGeoTranslation(19.75,18.8,-0.35)); |
| 170 |
GetVol("TPLA")->AddNode(GetVol("CATA"),1,new TGeoCombiTrans(0.,0.,1.,GetRot("rot2"))); |
| 171 |
}; |
| 172 |
|
| 173 |
void PamVMCDetCat::DefineCuts(){ |
| 174 |
|
| 175 |
TString s ="ALUMINIUM"; // default GPAMELA parameters |
| 176 |
if (GetMed(s) && !GetCC(s)) |
| 177 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 178 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 179 |
s ="N2 GAS"; |
| 180 |
if (GetMed(s) && !GetCC(s)) |
| 181 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 182 |
0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.)); |
| 183 |
s ="SCINTILLATOR";// default GPAMELA parameters |
| 184 |
if (GetMed(s) && !GetCC(s)) |
| 185 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 186 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 187 |
s = "PLASTIC";// default GPAMELA parameters |
| 188 |
if (GetMed(s) && !GetCC(s)) |
| 189 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 190 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 191 |
s = "TITANIUM Mixture"; // default GPAMELA parameters |
| 192 |
if (GetMed(s) && !GetCC(s)) |
| 193 |
SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, |
| 194 |
0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.)); |
| 195 |
|
| 196 |
} |
Parent Directory
| 
Revision Log