/[PAMELA software]/trieste/pamVMC/s4/src/PamVMCDetS4.cxx
ViewVC logotype

Annotation of /trieste/pamVMC/s4/src/PamVMCDetS4.cxx

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (hide annotations) (download)
Wed Mar 4 12:51:17 2009 UTC (15 years, 9 months ago) by pamelats
Branch point for: MAIN, pamVMC
Initial revision

1 pamelats 1.1 #include "PamVMCDetS4.h"
2    
3     ClassImp(PamVMCDetS4)
4    
5     void PamVMCDetS4::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     TGeoMedium *n=0;
27    
28     s="ALUMINIUM";
29     n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03);
30    
31     s = "SCINTILLATOR";
32     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);
33    
34     TGeoVolume *v = 0;
35    
36     s = "S4AL";
37     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),24.2,24.2,0.6);
38     s = "S4";
39     v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SCINTILLATOR"),24.1,24.1,0.5);
40    
41    
42     GetVol("S4AL")->AddNode(GetVol("S4"),1,gGeoIdentity);
43    
44     SetMotherProp(GetVol("S4AL"),1,new TGeoTranslation(0.,0.,1.811));
45    
46     };
47    
48     void PamVMCDetS4::DefineCuts(){
49    
50     TString s ="ALUMINIUM"; // default GPAMELA parameters
51     if (GetMed(s) && !GetCC(s))
52     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
53     0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
54     s ="N2 GAS";
55     if (GetMed(s) && !GetCC(s))
56     SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
57     0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.));
58     }
59    

  ViewVC Help
Powered by ViewVC 1.1.23