s4/src/PamVMCDetS4.cxx

#include "PamVMCDetS4.h"

ClassImp(PamVMCDetS4)

void  PamVMCDetS4::DefineGeometry(){

 TGeoMaterial *m=0;
 TGeoMixture *x=0;
  TString s = "ALUMINIUM"; 
  m = GetMat(s);
  if(!m){
    m = new TGeoMaterial(s,26.98,13,2.7);
    m->SetUniqueID(   9);
  }


 s = "SCINT";
  x = (TGeoMixture*)GetMat(s);  
  if(!x){
    x = new TGeoMixture(s,2,   1.03200    );
    x->SetUniqueID(  17);
    x->DefineElement(0,12.01,6,0.922427);
    x->DefineElement(1,1.01,1,0.7757296E-01);
  }

  TGeoMedium *n=0;

  s="ALUMINIUM";
  n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03);

  s = "SCINTILLATOR";
  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);

  TGeoVolume *v = 0;

  s = "S4AL";
  v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),24.2,24.2,0.6);
  s = "S4";
  v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SCINTILLATOR"),24.1,24.1,0.5);


  GetVol("S4AL")->AddNode(GetVol("S4"),1,gGeoIdentity);

  SetMotherProp(GetVol("S4AL"),1,new TGeoTranslation(0.,0.,1.811));

};

void  PamVMCDetS4::DefineCuts(){

 TString  s ="ALUMINIUM"; // default GPAMELA parameters
  if (GetMed(s) && !GetCC(s))
  SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
                             0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
  s ="N2 GAS";
  if (GetMed(s) && !GetCC(s))
    SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
                             0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.));
}

1	nikolas	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