Commit 3c4e1a3d authored by Warren's avatar Warren
Browse files

TACTIC Garfield electron drift functionality

parent 21007d9b
option(WITH_GARFIELD "Build with Garfield for electron drift in non-uniform field" OFF)
#IMPORTANT::IF CHANGING ABOVE OPTION DELETE CMakeCache IN nps FOR IT TO TAKE AFFECT. IF UPLOADING TO GIT SET THIS OPTION TO OFF, OTHERWISE A FULL COMPILATION OFF ALL DETECTORS IN NPTOOLS WILL REQUIRE GARFIELD!!!
if(WITH_GARFIELD)
find_package(Garfield REQUIRED)
include_directories($ENV{GARFIELD_HOME}/Include)
link_directories($ENV{GARFIELD_HOME}/Library)
file(GLOB headers $ENV{GARFIELD_HOME}/Include/*.hh)
endif()
add_library(NPSTACTIC SHARED TACTIC.cc TACTICScorer.cc)
if(WITH_GARFIELD)
target_link_libraries(NPSTACTIC NPSCore ${ROOT_LIBRARIES} ${Geant4_LIBRARIES} ${NPLib_LIBRARIES} -lNPTACTIC /Users/warren/Garfield/install/lib/libGarfield.0.3.0.dylib /Users/warren/Garfield/install/lib/libGarfield.dylib)
else()
target_link_libraries(NPSTACTIC NPSCore ${ROOT_LIBRARIES} ${Geant4_LIBRARIES} ${NPLib_LIBRARIES} -lNPTACTIC)
endif()
#include "Garfield/GeometrySimple.hh"
#include "Garfield/Sensor.hh"
#include "Garfield/ComponentAnalyticField.hh"
#include "Garfield/MediumMagboltz.hh"
#include "Garfield/TrackHeed.hh"
#include "Garfield/SolidTube.hh"
#include "Garfield/AvalancheMC.hh"
#include "Garfield/ViewSignal.hh"
#include "Garfield/ViewMedium.hh"
void GARFDRIFT(double energy, double t_start, G4ThreeVector start, G4ThreeVector delta_pos, double R, int Pad_start, double ID,
double ScoLength, double SegLength) {
const double rWire = 1.2;
const double rTube = 5.;
const double lTube = 25.19;
//For R6 = 2, R7 = 1
//const double vWire = -375.; //For V_total = 750V
//const double vWire = -645.3312; //For V_total = 1200V //SHOULD BE 583.741 !?
//const double vWire = -535.0959; //For V_total = 1100V
//const double vWire = -875.6116; //For V_total = 1800V
const double vWire = -681.0312; //For V_total = 1400V
const double vTube = 0.;
double x_start, y_start, z_start;
ofstream file;
double production_bias = 1.e-01;
double t_end, x_end, y_end, z_end;
double Pad_end;
//double w_value = 35.; //for HeCO2 --guess
double w_value = 26.31; //for argon
int electrons = (int)(energy/w_value*production_bias);//*production_bias; //Need to carry excess energy over later.
Garfield::MediumMagboltz* gas = new Garfield::MediumMagboltz();
Garfield::ViewMedium* mediumView = new Garfield::ViewMedium();
Garfield::GeometrySimple* geo = new Garfield::GeometrySimple();
Garfield::SolidTube* tube = new Garfield::SolidTube(0, 0, 0, 0, rTube, lTube/2);
Garfield::ComponentAnalyticField* cmp = new Garfield::ComponentAnalyticField();
Garfield::Sensor* sensor = new Garfield::Sensor();
Garfield::AvalancheMC* drift = new Garfield::AvalancheMC();
cout << "\n" << electrons << "\n" << endl;
if(R>rWire && electrons>0) {
gas->LoadGasFile("ar_90_ch4_10.gas"); //atomic fraction in this case (P10).
//gas->LoadGasFile("he_90_co2_10_500mbar.gas"); //mass fraction in this case
gas->Initialise(true);
mediumView->SetMedium(gas);
geo->AddSolid(tube, gas);
cmp->SetGeometry(geo);
cmp->AddWire(0, 0, 2*rWire, vWire, "c");
cmp->AddTube(rTube, vTube, 0, "a");
sensor->AddComponent(cmp);
drift->SetSensor(sensor);
for(int e=0;e<electrons;e++) {
double randomize = (double)std::rand() / (double)RAND_MAX;
x_start = start.x() + delta_pos.x()*randomize;
y_start = start.y() + delta_pos.y()*randomize;
z_start = start.z() + delta_pos.z()*randomize;
drift->DriftElectron(x_start, y_start, z_start, t_start);
int np = drift->GetNumberOfDriftLinePoints();
drift->GetDriftLinePoint(np-1, x_end, y_end, z_end, t_end); //np-1 is last DriftLineEndPoint
Pad_end = (int)((z_end + ScoLength / 2.) / SegLength ) + 1; //new Pad number
file.open("signal.dat", std::ios::app);
file << Pad_end << "\t" << t_end << endl;
file.close();
}
}
delete gas; delete mediumView; delete geo; delete tube; delete cmp; delete sensor; delete drift; //otherwise these are held in memory and cause a killed: 9 crash.
}
......@@ -123,7 +123,7 @@ G4LogicalVolume* TACTIC::BuildCylindricalDetector(){
G4Material* Cu = MaterialManager::getInstance()->GetMaterialFromLibrary("Cu");
G4Material* Mylar = MaterialManager::getInstance()->GetMaterialFromLibrary("Mylar");
G4Material* Vacuum = MaterialManager::getInstance()->GetMaterialFromLibrary("Vacuum");
unsigned const int NumberOfGasMix = m_GasMaterial.size();
double density=0;
......@@ -131,6 +131,17 @@ G4LogicalVolume* TACTIC::BuildCylindricalDetector(){
for(unsigned int i=0; i<NumberOfGasMix; i++){
if(m_GasMaterial[i] == "CO2") GasComponent.push_back(MaterialManager::getInstance()->GetGasFromLibrary(m_GasMaterial[i], 1.0/bar, m_Temperature));
if(m_GasMaterial[i] == "P10_gas") {
//G4Material *P10_gas = new G4Material("P10_gas", 0.00156 * g /cm3, 3, kStateGas, m_Temperature, 1.0/bar); //density for SRIM (1 atm);
G4Material *P10_gas = new G4Material("P10_gas", 0.00156*g/cm3, 3);
G4Element *elAr = new G4Element("Argon","Ar",18.,39.948*g/mole);
G4Element *elC = new G4Element("Carbon","C",6.,12.0107*g/mole);
G4Element *elH = new G4Element("Hydrogen","H",1.,1.00784*g/mole);
P10_gas->AddElement(elAr,90);
P10_gas->AddElement(elC,2);
P10_gas->AddElement(elH,8);
GasComponent.push_back(P10_gas);
}
else GasComponent.push_back(MaterialManager::getInstance()->GetMaterialFromLibrary(m_GasMaterial[i]));
}
......@@ -153,6 +164,7 @@ G4LogicalVolume* TACTIC::BuildCylindricalDetector(){
m_CylindricalDetector = new G4LogicalVolume(anode, Cu, "anode_log",0,0,0);
gas_volume_log = new G4LogicalVolume(gas_volume, TACTIC_gas, "gas_volume_log",0,0,0);
window_log = new G4LogicalVolume(window, Mylar, "window_log",0,0,0);
//window_log = new G4LogicalVolume(window, TACTIC_gas, "window_log",0,0,0); //windows removed (effectively)
vacuum_log = new G4LogicalVolume(vacuum, Vacuum, "vacuum_log",0,0,0);
new G4PVPlacement(0,G4ThreeVector(0,0,0),gas_volume_log,"gas_volume_phys",m_CylindricalDetector,false,0);
......@@ -251,8 +263,8 @@ void TACTIC::ConstructDetector(G4LogicalVolume* world){
G4FastSimulationManager* mng = m_ReactionRegion->GetFastSimulationManager();
unsigned int size = m_ReactionModel.size();
for(unsigned int i = 0 ; i < size ; i++){
mng->RemoveFastSimulationModel(m_ReactionModel[i]);
for(unsigned int j = 0 ; j < size ; j++){
mng->RemoveFastSimulationModel(m_ReactionModel[j]);
}
m_ReactionModel.clear();
......@@ -287,16 +299,22 @@ void TACTIC::ReadSensitive(const G4Event* event ){
G4int BeamCollectionID = G4SDManager::GetSDMpointer()->GetCollectionID("TACTICScorer/BeamScorer");
BeamHitMap = (G4THitsMap<G4double*>*)(event->GetHCofThisEvent()->GetHC(BeamCollectionID));
file.open("/Users/warren/nptool/NPSimulation/Detectors/TACTIC/out.dat",std::ios::app);
file.open("signal.dat", std::ios::app);
file << "Event" << endl;
file.close();
file.open("out.dat",std::ios::app);
for (Beam_itr = BeamHitMap->GetMap()->begin(); Beam_itr != BeamHitMap->GetMap()->end(); Beam_itr++) {
G4double* Info = *(Beam_itr->second);
//file << floor(((Info[3]+TACTIC_NS::active_length*0.5)/(TACTIC_NS::active_length/TACTIC_NS::NumberOfStrips))) << "\t"; // To get PAD number
file << event->GetEventID() << "\t";
for(int s = 0; s<12; s++) {
if(s==11) file << Info[s] << endl;
else file << Info[s] << "\t";
for(int s = 0; s<11; s++) {
//if(s==12) file << Info[s] << endl;
//else
file << Info[s] << "\t";
}
file << Info[11] << endl;
}
BeamHitMap->clear();
......@@ -310,10 +328,12 @@ void TACTIC::ReadSensitive(const G4Event* event ){
G4double* Info = *(Eject_itr->second);
//file << floor(((Info[3]+TACTIC_NS::active_length*0.5)/(TACTIC_NS::active_length/TACTIC_NS::NumberOfStrips))) << "\t"; // To get PAD number
file << event->GetEventID() << "\t";
for(int s = 0; s<12; s++) {
if(s==11) file << Info[s] << endl;
else file << Info[s] << "\t";
for(int s = 0; s<11; s++) {
//if(s==12) file << Info[s] << endl;
//else
file << Info[s] << "\t";
}
file << Info[11] << endl;
}
EjectHitMap->clear();
......@@ -334,7 +354,8 @@ void TACTIC::InitializeScorers() {
// Otherwise the scorer is initialised
G4VPrimitiveScorer* EjectScorer = new TACTICScorer::Gas_Scorer("EjectScorer",1,TACTIC_NS::active_length,(int)TACTIC_NS::NumberOfStrips);
G4VPrimitiveScorer* BeamScorer = new TACTICScorer::Gas_Scorer("BeamScorer",1,TACTIC_NS::active_length,(int)TACTIC_NS::NumberOfStrips);
G4SDParticleFilter* EjectFilter = new G4SDParticleFilter("EjectFilter","proton");
//G4SDParticleFilter* EjectFilter = new G4SDParticleFilter("EjectFilter","proton");
G4SDParticleFilter* EjectFilter = new G4SDParticleFilter("EjectFilter","alpha"); //For studying alpha source data
G4SDParticleFilter* BeamFilter = new G4SDParticleFilter("BeamFilter");
BeamFilter->addIon(11,21);
BeamFilter->addIon(10,18);
......
......@@ -2,11 +2,11 @@
TACTIC
POS= 0 0 0 mm
Shape= Cylindrical
% GasMaterial_1 = P10_1atm
GasMaterial_1 = He_gas
GasMaterial_1 = P10_gas
% GasMaterial_1 = He_gas
GasMaterial_2 = CO2
GasFraction_1 = 90
GasFraction_2 = 10
GasFraction_1 = 100
GasFraction_2 = 0
Temperature = 293.15 kelvin
Pressure = 0.5 bar
% Pressure = 1.0 bar
......@@ -91,7 +91,7 @@ public: // Scorer
void InitializeScorers() ;
// Associated Scorer
G4MultiFunctionalDetector* m_Scorer ;
G4MultiFunctionalDetector* m_Scorer ;
////////////////////////////////////////////////////
///////////Event class to store Data////////////////
////////////////////////////////////////////////////
......
//#define USE_Garfield //only use if compiled with Garfield
#include "TACTICScorer.hh"
#include "G4UnitsTable.hh"
#ifdef USE_Garfield
#include "GARFDRIFT.h"
#endif
using namespace TACTICScorer;
Gas_Scorer::Gas_Scorer(G4String name,G4int Level,G4double ScorerLength,G4int NumberOfSegments, G4int depth) //what do level and depth do?
......@@ -8,7 +15,6 @@ Gas_Scorer::Gas_Scorer(G4String name,G4int Level,G4double ScorerLength,G4int Num
m_NumberOfSegments = NumberOfSegments;
m_SegmentLength = m_ScorerLength / m_NumberOfSegments;
m_Level = Level;
m_Position = G4ThreeVector(-1000,-1000,-1000);
m_SegmentNumber = -1;
m_Index = -1;
......@@ -19,6 +25,7 @@ Gas_Scorer::~Gas_Scorer(){}
G4bool Gas_Scorer::ProcessHits(G4Step* aStep, G4TouchableHistory*){
G4double* Infos = new G4double[12];
//G4double w_value = 26.31*eV;
m_Position = aStep->GetPreStepPoint()->GetPosition();
Infos[0] = aStep->GetTrack()->GetTrackID();
......@@ -49,16 +56,55 @@ G4bool Gas_Scorer::ProcessHits(G4Step* aStep, G4TouchableHistory*){
m_DetectorNumber = aStep->GetPreStepPoint()->GetTouchableHandle()->GetCopyNumber(m_Level);
m_Index = m_DetectorNumber * 1e3 + m_SegmentNumber * 1e6;
//Infos[12] = 0.;
//Infos[13] = 0.;
//Infos[14] = 0.;
//cout << "\n" << "PAD " << Infos[5] << " first_step = " << first_step << endl;
//double TOA_min, TOA_max;
//if(PAD==1000) PAD = Infos[5];
//if(PAD!=Infos[5]) TOA_PAD.clear(); // if new pad clear the TOA vec
#ifdef USE_Garfield
G4ThreeVector delta_Position = aStep->GetDeltaPosition();
// bool last_step;
//if(aStep->GetTrack()->GetTrackStatus()!=fAlive) last_step = 1; else last_step = 0;
//vector<double> TOA_vec = GARFDRIFT(Infos[4]/eV, Infos[2], m_Position/cm, delta_Position/cm, Infos[7]/cm, Infos[5], Infos[1], m_ScorerLength/cm, m_SegmentLength/cm, last_step); //Ionization Energy Deposit, Global Time, x, y, z,R, ParticleID //Garfield works in cm, both G4 and Garfield work in ns.
GARFDRIFT(Infos[4]/eV, Infos[2], m_Position/cm, delta_Position/cm, Infos[7]/cm, Infos[5], Infos[1], m_ScorerLength/cm, m_SegmentLength/cm);
#endif
/*
for(int t=0; t<TOA_vec.size(); t++) TOA_PAD.push_back(TOA_vec[t]);
sort(TOA_PAD.begin(), TOA_PAD.end());
reverse(TOA_PAD.begin(), TOA_PAD.end());
if(TOA_PAD.size() > 42) TOA_max = TOA_PAD[42], TOA_min = *min_element(TOA_PAD.begin(), TOA_PAD.end()); //threshold of 42 (420 e- ~10 mV threshold)
else TOA_min = 1.e06, TOA_max = 0.;
cout << "\n" << "PAD: " << Infos[5] << " vec size: " << TOA_PAD.size() << " TOA_min: " << TOA_min << " TOA_max " << TOA_max << endl;
#else
TOA_min =1.e06;
TOA_max = 0.;
#endif
*/
map<G4int, G4double**>::iterator it;
it= EvtMap->GetMap()->find(m_Index);
if(it!=EvtMap->GetMap()->end()){ //accumulate ionisation energy deposit
if(it!=EvtMap->GetMap()->end()){
G4double* dummy = *(it->second);
Infos[4]+=dummy[4];
/*
if(TOA_min < 1.e06) Infos[12] = TOA_max - TOA_min;
if(Infos[12] < dummy[12]) Infos[12] = dummy[12]; //Ensures that max risetime is maintained
PAD = Infos[5];
*/
Infos[4]+=dummy[4]; //accumulate ionisation energy deposit to get total accross pad
delete dummy;
}
EvtMap->set(m_Index, Infos);
return TRUE;
}
......@@ -92,4 +138,3 @@ void Gas_Scorer::PrintAll(){
G4cout << " PrimitiveScorer " << GetName() << G4endl ;
G4cout << " Number of entries " << EvtMap->entries() << G4endl ;
}
......@@ -36,6 +36,10 @@ namespace TACTICScorer {
G4double m_SegmentLength;
// Level at which to find the copy number linked to the detector number
G4int m_Level;
//bool first_step = 1;
//vector<double> TOA_PAD;
//int PAD = 1000;
private: // inherited from G4VPrimitiveScorer
G4int HCID;
......
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%% An Isotropic Source to be used as EventGenerator %%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Energy are given in MeV , Position in mm %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Isotropic
EnergyLow= 3.99 MeV
EnergyHigh= 5.25 MeV
EnergyDistribution = 2.5859*exp(-0.5*((x-4.61924)/0.212526)*((x-4.61924)/0.212526))
HalfOpenAngleMin= 0 deg
HalfOpenAngleMax= 90 deg
x0= 0 mm
y0= 0 mm
z0= 0 mm
Particle= alpha
ExcitationEnergy= 0 MeV
% Supported particle type: proton, neutron, deuton, triton, He3 , alpha
*----.----1----.----2----.----3----.----4----.----5----.----6----.----7----.----8----.----9----.---10----.---11----.---12----.---13--
% Created 03/05/20 at 06.47.47 < none > GAS "none "
Version : 12
GASOK bits: TFTTFTTTTTTFFFTFFFFF
Identifier: Ar 90%, CH4 10%, T=293.15 K, p=0.5 atm
Clusters :
Dimension : F 20 1 1 1 0
E fields
2.63157895E-03 3.78539444E-03 5.44510021E-03 7.83250379E-03 1.12666642E-02
1.62065319E-02 2.33122840E-02 3.35335523E-02 4.82363345E-02 6.93855500E-02
9.98076366E-02 1.43568284E-01 2.06515782E-01 2.97062603E-01 4.27309668E-01
6.14663544E-01 8.84162707E-01 1.27182375E+00 1.82945473E+00 2.63157895E+00
E-B angles
1.57079633E+00
B fields
0.00000000E+00
Mixture:
0.00000000E+00 9.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 1.00000000E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
Excitation 1: "CH4 EXCITATION TRIPLET DISSOCIATION ELOSS= 7.5 " 7.49974354E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
The gas tables follow:
6.73041747E-02 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 5.03945490E+00 0.00000000E+00
4.19725506E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 2.53961057E+05 0.00000000E+00 1.78834769E+05
0.00000000E+00 1.73504231E+05 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1.07729021E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 4.45483495E+00 0.00000000E+00
3.70785335E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 1.98455545E+05 0.00000000E+00 1.35720387E+05
0.00000000E+00 1.39243142E+05 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1.69143411E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 3.90591956E+00 0.00000000E+00
3.21758934E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 1.52562076E+05 0.00000000E+00 1.04434329E+05
0.00000000E+00 1.02623294E+05 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
2.59721209E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 3.52874606E+00 0.00000000E+00
2.77592734E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 1.24520488E+05 0.00000000E+00 7.64911791E+04
0.00000000E+00 7.76242729E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
4.07272466E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 3.05279517E+00 0.00000000E+00
2.39944125E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 9.31955833E+04 0.00000000E+00 5.80973811E+04
0.00000000E+00 5.70489851E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
6.46116484E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.72984306E+00 0.00000000E+00
2.07594583E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 7.45204314E+04 0.00000000E+00 4.26903078E+04
0.00000000E+00 4.35007138E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1.04347315E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.55930458E+00 0.00000000E+00
1.83312095E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 6.55003992E+04 0.00000000E+00 3.35618885E+04
0.00000000E+00 3.36447598E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1.68372363E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.27663142E+00 0.00000000E+00
1.75201487E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 5.18305064E+04 0.00000000E+00 3.10186763E+04
0.00000000E+00 3.03724461E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
2.64033987E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.07885519E+00 0.00000000E+00
1.64401879E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 4.32163892E+04 0.00000000E+00 2.69444837E+04
0.00000000E+00 2.71114720E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
3.81014290E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.72588805E+00 0.00000000E+00
1.62357123E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 2.97868956E+04 0.00000000E+00 2.72752830E+04
0.00000000E+00 2.54443879E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
4.87744032E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.32623216E+00 0.00000000E+00
1.56473624E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 1.75889175E+04 0.00000000E+00 2.41543701E+04
0.00000000E+00 2.48136196E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
5.47489127E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.06122227E+00 0.00000000E+00
1.63262780E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 1.12619270E+04 0.00000000E+00 2.73900097E+04
0.00000000E+00 2.59194611E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
5.48144361E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 8.95874146E-01 0.00000000E+00
1.61021889E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 8.02590485E+03 0.00000000E+00 2.57090430E+04
0.00000000E+00 2.61470543E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
5.03005153E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 8.16509446E-01 0.00000000E+00
1.62205396E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 6.66687675E+03 0.00000000E+00 2.55913609E+04
0.00000000E+00 2.70298202E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
4.33459998E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 6.87021884E-01 0.00000000E+00
1.61659155E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 4.71999068E+03 0.00000000E+00 2.57070149E+04
0.00000000E+00 2.65603501E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
3.60028353E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 6.34991303E-01 0.00000000E+00
1.61913123E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 4.03213955E+03 0.00000000E+00 2.59601702E+04
0.00000000E+00 2.64715485E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
2.99342239E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 5.91569608E-01 0.00000000E+00
1.67571767E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 3.49954601E+03 0.00000000E+00 2.89889805E+04
0.00000000E+00 2.71716139E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
2.59131814E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 5.65781288E-01 0.00000000E+00
1.66416573E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 3.20108465E+03 0.00000000E+00 2.76265955E+04
0.00000000E+00 2.77623563E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
2.38115105E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 5.79188977E-01 0.00000000E+00
1.62948795E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 3.35459871E+03 0.00000000E+00 2.68351809E+04
0.00000000E+00 2.62694389E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
2.30574815E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 5.75344601E-01 0.00000000E+00
1.51690722E+00 0.00000000E+00-3.59401713E+01 0.00000000E+00-3.59401713E+01-3.59401713E+01 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00-3.00000000E+01 0.00000000E+00 3.31021410E+03 0.00000000E+00 2.39849111E+04
0.00000000E+00 2.20352392E+04 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
H Extr: 1 1 1 1 1 1 1 1 1 1 1 1 1
L Extr: 0 0 0 0 0 0 0 0 0 0 0 0 0
Thresholds: 1 1 1
Interp: 2 2 2 2 2 2 2 2 2 2 2 2 2
A = 0.00000000E+00, Z = 0.00000000E+00, EMPROB= 0.00000000E+00, EPAIR = 0.00000000E+00
Ion diffusion: 0.00000000E+00 0.00000000E+00
CMEAN = 0.00000000E+00, RHO = 0.00000000E+00, PGAS = 3.80000000E+02, TGAS = 2.93150000E+02
CLSTYP : NOT SET
FCNCLS :
NCLS : 0
Average : 0.000000000000000000E+00
Heed initialisation done: F
SRIM initialisation done: F
This diff is collapsed.
/run/beamOn 100
/run/beamOn 10
/exit
\ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment