/[PAMELA software]/PamVMC/src/PamVMCDetPamela.cxx
ViewVC logotype

Annotation of /PamVMC/src/PamVMCDetPamela.cxx

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.5 - (hide annotations) (download)
Fri Jun 12 18:39:29 2009 UTC (15 years, 5 months ago) by pam-rm2
Branch: MAIN
CVS Tags: v1r0, HEAD
Changes since 1.1: +0 -0 lines 
- Introduced user-defined names of output files and random seeds number.
Users can do it use options of PamVMCApplication constructor:
PamVMCApplication(const char* name,  const char *title, const char*
filename="pamtest", Int_t seed=0).
The Random object that I use is TRandom3 object which has astronomical
large period (in case of default initialization 0). All random generators
in the code use this object by calling of gRandom singleton which keeps
it.

- Corrected TOF digitization routine. No problems with TDC hits due to
hadronic interactions anymore.

- Some small changes was done to compile code under Root 5.23. +
geant4_vmc v. 2.6 without any warnings

- Some classes of PamG4RunConfiguartion was changed for geant4_vmc v.
2.6.Some obsolete classes was deleted as soon as developers implemented
regions.

- Navigation was changed from "geomRootToGeant4" to "geomRoot", because on
VMC web page written that as soon as Geant4 has no option ONLY/MANY
translation of overlapped geometry to Geant4 through VGM could be wrong.
I'd like to stay with Root navigation:
http://root.cern.ch/root/vmc/Geant4VMC.html. This should be default
option.

- New Tracker digitization routine written by Sergio was implemented

- PamVMC again became compatible with geant4_vmc v.2.5 and ROOT 5.20.
 The problem was that ROOT developers introduced in TVirtualMC class a new
method SetMagField and new base class:TVirtualMagField from which
user-defined classes shoukd be derived
1 nikolas 1.1 #include "PamVMCDetPamela.h"
2    
3     ClassImp(PamVMCDetPamela)
4    
5     void PamVMCDetPamela::DefineMother(){
6    
7     TGeoMaterial *m =0;
8     TString s = "N2GAS";
9     m = GetMat(s);
10     if (!m){
11     m = new TGeoMaterial(s,14.01,7,0.1250000E-02);
12     m->SetUniqueID( 19);
13     }
14     TGeoMedium *n = 0;
15     s = "N2 GAS";
16     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);
17    
18     TGeoVolume *v = 0;
19    
20     s = "PAME";
21     v = (GetVol(s))?GetVol(s): gGeoManager->MakeTube(s,GetMed("N2 GAS"),0.,53.,130.);
22     gGeoManager->SetTopVolume(v);
23    
24     };
25    
26     void PamVMCDetPamela::DefineShell(){
27    
28     TGeoMaterial *m = 0;
29     TGeoMixture *x = 0;
30     TString s = "ALUMINIUM";
31     m = GetMat(s);
32     if(!m){
33     m = new TGeoMaterial(s,26.98,13,2.7);
34     m->SetUniqueID( 9);
35     }
36    
37     s = "MYLAR";
38     x = (TGeoMixture*) GetMat(s);
39     if(!x){
40     x = new TGeoMixture(s,3, 1.39000 );
41     x->SetUniqueID( 29);
42     x->DefineElement(0,12.011,6,0.6250159);
43     x->DefineElement(1,1.00794,1,0.4196011E-01);
44     x->DefineElement(2,15.9994,8,0.3330241);
45     }
46    
47     s = "NYLON6";
48     x = (TGeoMixture*)GetMat(s);
49     if(!x){
50     x = new TGeoMixture(s,4, 1.14000 );
51     x->SetUniqueID( 35);
52     x->DefineElement(0,12.011,6,0.6368534);
53     x->DefineElement(1,1.00794,1,0.9797978E-01);
54     x->DefineElement(2,15.9994,8,0.1413881);
55     x->DefineElement(3,14.00674,7,0.1237788);
56     }
57    
58     s = "AMG6M";
59     x = (TGeoMixture*)GetMat(s);
60     if(!x){
61     x = new TGeoMixture(s,2, 2.61000 );
62     x->SetUniqueID( 36);
63     x->DefineElement(0,26.9815,13,0.94);
64     x->DefineElement(1,24.305,12,0.6000000E-01);
65     }
66    
67     TGeoMedium *n = 0;
68    
69    
70     s = "ALUMINIUM";
71     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03);
72    
73     s = "MYLAR";
74     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,19,29,-1,1,40,1,100,0.5000000E-01,0.2300000E-04,0.5000000E-03);
75    
76     s = "Nylon6 Mixture";
77     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,25,35,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.1495361E-04,0.5000000E-03);
78    
79     s = "AMg6M Mixture";
80     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,26,36,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.1699982E-02,0.5000000E-03);
81    
82    
83     TGeoVolume *v = 0;
84    
85     s = "SHEL";
86     v = (GetVol(s))?GetVol(s): gGeoManager->MakeCone(s,GetMed("ALUMINIUM"),38.8717,52.5,52.651,43.6423,43.7932);
87     s = "TSPH";
88     v = (GetVol(s))?GetVol(s): gGeoManager->MakeSphere(s,GetMed("AMg6M Mixture"),240.,240.17,0.,10.5062,0.,360.);
89     s = "COV1";
90     v = (GetVol(s))?GetVol(s): gGeoManager->MakeSphere(s,GetMed("Nylon6 Mixture"),240.17,240.1715,0.,10.5062,0.,360.);
91     s = "COV2";
92     v = (GetVol(s))?GetVol(s): gGeoManager->MakeSphere(s,GetMed("MYLAR"),240.1715,240.2165,0.,10.5062,0.,360.);
93     s = "COV3";
94     v = (GetVol(s))?GetVol(s): gGeoManager->MakeSphere(s,GetMed("ALUMINIUM"),240.2165,240.2315,0.,10.5062,0.,360.);
95     s = "MSHE";
96     v = (GetVol(s))?GetVol(s): gGeoManager->MakeTube(s,GetMed("ALUMINIUM"),52.5,52.75,19.95);
97     s = "BSPH";
98     v = (GetVol(s))?GetVol(s): gGeoManager->MakeSphere(s,GetMed("ALUMINIUM"),52.5,52.75,90.,180.,0.,360.);
99    
100     gGeoManager->GetTopVolume()->AddNodeOverlap(GetVol("SHEL"),1,new TGeoTranslation(0.,0.,66.5559)); //66.5607 Idid it because I need to check... This volume overlaps with TSPH and COV1
101     SetMotherProp(GetVol("TSPH"),1,new TGeoTranslation(0.,0.,-130.716));
102     SetMotherProp(GetVol("COV1"),1,new TGeoTranslation(0.,0.,-130.716)); //-130.7112
103     SetMotherProp(GetVol("COV2"),1,new TGeoTranslation(0.,0.,-130.716));
104     SetMotherProp(GetVol("COV3"),1,new TGeoTranslation(0.,0.,-130.716));
105     SetMotherProp(GetVol("MSHE"),1,new TGeoTranslation(0.,0.,7.7342)); //7.739
106     SetMotherProp(GetVol("BSPH"),1,new TGeoTranslation(0.,0.,-12.2158));//-12.211
107    
108     };
109    
110     void PamVMCDetPamela::DefineCuts(){
111    
112     TString s ="ALUMINIUM"; // default GPAMELA parameters
113     if (GetMed(s) && !GetCC(s))
114     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
115     0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
116     s ="N2 GAS";
117     if (GetMed(s) && !GetCC(s))
118     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
119     0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.));
120     s ="MYLAR";
121     if (GetMed(s) && !GetCC(s))
122     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
123     0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
124     s ="Nylon6 Mixture";
125     if (GetMed(s) && !GetCC(s))
126     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
127     0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
128     s ="AMg6M Mixture";
129     if (GetMed(s) && !GetCC(s))
130     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
131     0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
132    
133     }
134    
135    
  ViewVC Help
Powered by ViewVC 1.1.23