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	#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