/[PAMELA software]/PamVMC/trk/src/PamVMCDetTrk.cxx
ViewVC logotype

Annotation of /PamVMC/trk/src/PamVMCDetTrk.cxx

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.5 - (hide annotations) (download)
Fri Jun 12 18:39:56 2009 UTC (15 years, 6 months ago) by pam-rm2
Branch: MAIN
CVS Tags: v1r0, HEAD
Changes since 1.1: +12 -2 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 "PamVMCDetTrk.h"
2    
3     ClassImp(PamVMCDetTrk)
4    
5     void PamVMCDetTrk::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     s = "IRON";
17     m = GetMat(s);
18     if(!m){
19     m = new TGeoMaterial(s,55.85,26,7.87);
20     m->SetUniqueID( 10);
21     }
22    
23     s = "SCINT";
24     x = (TGeoMixture*)GetMat(s);
25     if(!x){
26     x = new TGeoMixture(s,2, 1.03200 );
27     x->SetUniqueID( 17);
28     x->DefineElement(0,12.01,6,0.922427);
29     x->DefineElement(1,1.01,1,0.7757296E-01);
30     }
31    
32     s = "N2GAS";
33     m = GetMat(s);
34     if (!m){
35     m = new TGeoMaterial(s,14.01,7,0.1250000E-02);
36     m->SetUniqueID( 19);
37     }
38    
39     s = "SILICON TR";
40     m= GetMat(s);
41     if(!m){
42     m = new TGeoMaterial(s,28.09,14,2.33);
43     m->SetUniqueID( 20);
44     }
45    
46     s = "G10C";
47     x = (TGeoMixture*)GetMat(s);
48     if(!x){
49     x = new TGeoMixture(s,4, 1.70000 );
50     x->SetUniqueID( 27);
51     x->DefineElement(0,28.09,14,0.53);
52     x->DefineElement(1,16,8,0.3);
53     x->DefineElement(2,12.01,6,0.15);
54     x->DefineElement(3,1.01,1,0.2000000E-01);
55     }
56    
57    
58     TGeoMedium * n = 0;
59    
60     s="ALUMINIUM";
61     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03);
62    
63     s = "N2 GAS";
64     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);
65    
66     s = "IRONTR";
67     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,5,10,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.8000000E-01,0.5000000E-03);
68    
69     s = "CARBON FIBER";
70     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,9,17,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.5000000E-02,0.5000000E-03); //CARBON FIBER MADE FROM SCINTILLATOR???
71    
72     s = "SITRACKER";
73     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,12,20,1,1,40,1,0.1000000E-02,0.5000000E-01,0.5000000E-04,0.1569056E-01);
74    
75     s = "G10C";
76     n = (GetMed(s))?GetMed(s):new TGeoMedium(s,17,27,-1,1,40,1,100,0.5000000E-01,0.1200000E-02,0.5000000E-03);
77    
78     TGeoVolume *v = 0;
79    
80     s = "SPEB";
81     v = (GetVol(s))?GetVol(s):gGeoManager->MakeBox(s,GetMed("N2 GAS"),15.,14.4,22.615);
82     s = "MGFR";
83     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),15.,14.4,4.45);
84     s = "MGPL";
85     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("IRONTR"),12.,11.4,4.);
86     s = "MGPI";
87     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),8.07,6.57,4.);
88     s = "TRPB";
89     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),10.49,12.15,0.365);
90     s = "TPAS";
91     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),8.2,9.9,0.365);
92 pam-rm2 1.5 //my new sensitive volume//
93     s = "TPAN";
94     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),8.2,9.9,0.15E-01);
95     //////////////////////////
96 nikolas 1.1 s = "TPAI";
97     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),8.144,1.25,0.5E-01);
98     s = "TRSL";
99     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SITRACKER"),2.6665,3.5,0.15E-01);
100     s = "THBP";
101     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("G10C"),2.6665,2.75,0.15E-01);
102     s = "TSPA";
103     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SITRACKER"),2.596,3.413,0.15E-01);
104     s = "TRCP";
105     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("CARBON FIBER"),0.35E-01,9.9,0.25);
106 pam-rm2 1.5 //my new sensitive volume//
107     s = "TRCN";
108     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("CARBON FIBER"),0.35E-01,9.9,0.15E-01);
109     //////////////////////////
110 nikolas 1.1 s = "TBAL";
111     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),1.49,2.25,0.365);
112     s = "MGPA";
113     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),8.1,6.6,4.);
114     s = "TPGA";
115     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),10.29,10.65,0.4E-01);
116     s = "TPGI";
117     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),8.07,6.57,0.15E-01);
118     s = "TPGU";
119     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),10.49,12.15,0.1E-01);
120     s = "TPGD";
121     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),12.,11.4,0.5E-02);
122    
123     GetVol("SPEB")->AddNode(GetVol("MGFR"),1,new TGeoTranslation(0.,0.,18.165));
124     GetVol("MGFR")->AddNode(GetVol("TPGD"),1,new TGeoTranslation(0.,0.,-4.445));
125     GetVol("MGFR")->AddNode(GetVol("TPGI"),1,new TGeoTranslation(0.,0.,-4.425));
126     GetVol("MGFR")->AddNode(GetVol("MGPL"),1,new TGeoTranslation(0.,0.,-0.41));
127     GetVol("MGPL")->AddNode(GetVol("MGPA"),1,gGeoIdentity);
128     GetVol("MGPA")->AddNode(GetVol("MGPI"),1,gGeoIdentity);
129     GetVol("MGFR")->AddNode(GetVol("TPGI"),2,new TGeoTranslation(0.,0.,3.605));
130     GetVol("MGFR")->AddNode(GetVol("TPGA"),1,new TGeoTranslation(0.,-0.75,3.66));
131     GetVol("MGFR")->AddNode(GetVol("TRPB"),1,new TGeoTranslation(0.,-2.25,4.065));
132 pam-rm2 1.5 GetVol("TRPB")->AddNode(GetVol("TPAS"),4,new TGeoTranslation(0.,-0.75,0.));
133     GetVol("TPAS")->AddNodeOverlap(GetVol("TPAN"),1,new TGeoTranslation(0.,0.,0.15E-01));
134 nikolas 1.1 GetVol("TPAS")->AddNode(GetVol("TRCP"),1,new TGeoTranslation(-8.1075,0.,0.15E-01));
135 pam-rm2 1.5 GetVol("TRCP")->AddNode(GetVol("TRCN"),1,gGeoIdentity);
136 nikolas 1.1 GetVol("TPAS")->AddNode(GetVol("TRCP"),2,new TGeoTranslation(-2.7025,0.,0.15E-01));
137     GetVol("TPAS")->AddNode(GetVol("TRCP"),3,new TGeoTranslation(2.7025,0.,0.15E-01));
138     GetVol("TPAS")->AddNode(GetVol("TRCP"),4,new TGeoTranslation(8.1075,0.,0.15E-01));
139     GetVol("TPAS")->AddNode(GetVol("TRSL"),1,new TGeoTranslation(-5.405,6.252,0.15E-01));
140 pam-rm2 1.5 GetVol("TRSL")->AddNode(GetVol("TSPA"),1,gGeoIdentity);
141 nikolas 1.1 GetVol("TPAS")->AddNode(GetVol("TRSL"),2,new TGeoTranslation(0.,6.252,0.15E-01));
142     GetVol("TPAS")->AddNode(GetVol("TRSL"),3,new TGeoTranslation(5.405,6.252,0.15E-01));
143     GetVol("TPAS")->AddNode(GetVol("TRSL"),4,new TGeoTranslation(-5.405,-0.75,0.15E-01));
144     GetVol("TPAS")->AddNode(GetVol("TRSL"),5,new TGeoTranslation(0.,-0.75,0.15E-01));
145     GetVol("TPAS")->AddNode(GetVol("TRSL"),6,new TGeoTranslation(5.405,-0.75,0.15E-01));
146     GetVol("TPAS")->AddNode(GetVol("THBP"),1,new TGeoTranslation(-5.405,-7.0015,0.15E-01));
147     GetVol("TPAS")->AddNode(GetVol("THBP"),2,new TGeoTranslation(0.,-7.0015,0.15E-01));
148     GetVol("TPAS")->AddNode(GetVol("THBP"),3,new TGeoTranslation(5.405,-7.0015,0.15E-01));
149     GetVol("TPAS")->AddNode(GetVol("TPAI"),4,new TGeoTranslation(0.,-8.65,0.315));
150     GetVol("MGFR")->AddNode(GetVol("TPGU"),2,new TGeoTranslation(0.,-2.25,4.44));
151    
152     GetVol("SPEB")->AddNode(GetVol("MGFR"),2,new TGeoTranslation(0.,0.,9.265));
153     GetVol("SPEB")->AddNode(GetVol("MGFR"),3,new TGeoTranslation(0.,0.,0.365));
154     GetVol("SPEB")->AddNode(GetVol("MGFR"),4,new TGeoTranslation(0.,0.,-8.535));
155     GetVol("SPEB")->AddNode(GetVol("MGFR"),5,new TGeoTranslation(0.,0.,-17.435));
156    
157     GetVol("SPEB")->AddNode(GetVol("TRPB"),2,new TGeoCombiTrans(0.,2.25,-22.25,GetRot("rot10")));
158     GetVol("SPEB")->AddNode(GetVol("TBAL"),1,new TGeoTranslation(5.4,-12.15,-22.25));
159     GetVol("SPEB")->AddNode(GetVol("TBAL"),2,new TGeoTranslation(-5.4,-12.15,-22.25));
160    
161    
162     SetMotherProp(GetVol("SPEB"),1,new TGeoTranslation(0.,0.,49.274));
163     };
164    
165     void PamVMCDetTrk::DefineCuts(){
166    
167     TString s ="ALUMINIUM"; // default GPAMELA parameters
168     if (GetMed(s) && !GetCC(s))
169     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
170     0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
171     s ="N2 GAS";
172     if (GetMed(s) && !GetCC(s))
173     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
174     0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.));
175     s = "IRONTR"; // 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 = "CARBON FIBER"; //default GPAMELA parameters
180     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
181     0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
182     s = "SITRACKER";
183     if (GetMed(s) && !GetCC(s))
184     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.00001, 0.00001, 0.0001, 0.0001, 0.00001, 0.00001,
185     0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 1., 1., 1., 1., 2.));
186    
187     s = "G10C";
188     if (GetMed(s) && !GetCC(s))
189     SetCC(s, new pCutControl(GetMedID(s), 0.00001, 0.00001, 0.00001, 0.00001, 0.001, 0.001, 0.01, 0.00001,
190     0.00001, 0.01, 1., 1., 1., 0., 1., 1., 1., 1., 1., 1., 1., 2.));
191     }
192    
  ViewVC Help
Powered by ViewVC 1.1.23