diff --git a/exp1016/deadlayers.cxx b/exp1016/deadlayers.cxx new file mode 100644 index 0000000000000000000000000000000000000000..1ba8eaa2c1cee166d099cefd6400473b136ee11c --- /dev/null +++ b/exp1016/deadlayers.cxx @@ -0,0 +1,43 @@ +{ + gSystem->Load("../libTELoss.so"); + + TELoss mSi; +//set parameters for Si and alpha + mSi.SetEL(1, 2.321); // density in g/cm3 + mSi.AddEL(14., 28.086, 1); //Z, mass +// mSi.SetZP(1., 1.); //protons + mSi.SetZP(2., 4.); //alphas, Z, A + mSi.SetEtab(100000, 200.); // ?, MeV calculate ranges + mSi.SetDeltaEtab(300); // calculate delta E + +// TCanvas *c1 = new TCanvas("name", "title which can consist more than one word.", 617, 0, 1058, 972); +// cout << mSi.GetE(50, 1000) << endl; //(input energy E0 in MeV, microns) +// cout << mSi.GetE0(9.04, 1000) << endl; //(output energy E in MeV, microns) + +//calculated energies for run1 - run3 + + Double_t deltaL1 = mSi.GetR(4.789, 4.623); //(doule E0, double E in MeV) calculates layer in mcm for averages + Double_t deltaL2 = mSi.GetR(5.5,5.355); + Double_t deltaL3 = mSi.GetR(6.017,5.880); + Double_t deltaL4 = mSi.GetR(7.693, 7.577); + +//calculated energies for run4 - run5 +/* + Double_t deltaL1 = mSi.GetR(4.782,4.613); //(doule E0, double E in MeV) calculates layer in mcm for averages + Double_t deltaL2 = mSi.GetR(5.503,5.348); + Double_t deltaL3 = mSi.GetR(6.016,5.875); + Double_t deltaL4 = mSi.GetR(7.692,7.571); +*/ + cout << deltaL1 << " mcm" << endl; //MeV, microns - delta layer + cout << deltaL2 << " mcm" << endl; + cout << deltaL3 << " mcm" << endl; + cout << deltaL4 << " mcm" << endl << endl; + + cout << "dead layer 1 is: " << deltaL1/(TMath::Sqrt(2) - 1) << " mcm" << endl; //MeV, microns + cout << "dead layer 2 is: " <