5252017 Anal
LB_May_2017_Irradiation_Day#5.2F25.2F17
10% Se mixed with Sage leaf ash
Mass of each sample
M_{mix} =
M_{Se} =
M_{Ash}=
Lines for investigation
Found lines from Sage at energies of 147,250,389,1267,1369,1517,2169 and 2755 keV, All seem to have half lines longer than 1 hr.
Compare to Background lines measured at
https://wiki.iac.isu.edu/index.php/PAA_BackGrd_Det_A
Possible sources for the lines
- Mg-25(\gamma,p)Na-24
- :1369, 2755 keV ::T1/2=15 hrs
- Air
- :250, 1267, 1517 keV - Ar-40
147 keV - not sure. Depending on what other lines are there, could be Ag-36, Cl-35, Kr-78
389 keV Sr-88, depending on what else is there
2169 keV - don't know
Time cuts for analysis
0-300 (Sage + Se 10% Se by weight) 300-360 (Co-60) 400-640 ( Pure Se) 680-710 (Co-60) 730-1020 (Mixture) 1030-1080 (Co-60) 1100-1360 ( Pure Se) 1400-1440 (Co-60) 1480-1775 (Mixture) 1800-1840 (Co-60) 1875-2150 (Pure Se) 2190-2220 (Co-60) 2250-2550 (Mixture) 2590-2620 (Co-60) 2650-2930 (Pure Se) 2950-3000 (Co-60) 3050-3300 (Mixture) 3390-3350 (Co-60) 3400-3690 (Pure Se) 3700-3750 (Co-60) 3775-4050 (Mixture) 4060-4100 (Co-60) 4120-4400 (Pure Se) 4420-4470 (Co-60) 4480-4770 (Mixture) 4790-4820 (Co-60) 4840-5130 (Pure Se)
Analysis Code
#define TenPercSeInSage_cxx
#include "TenPercSeInSage.h"
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>
void TenPercSeInSage::Loop()
{
// In a ROOT session, you can do:
// root> .L TenPercSeInSage.C
// root> TenPercSeInSage t
// root> t.GetEntry(12); // Fill t data members with entry number 12
// root> t.Show(); // Show values of entry 12
// root> t.Show(16); // Read and show values of entry 16
// root> t.Loop(); // Loop on all entries
//
// This is the loop skeleton where:
// jentry is the global entry number in the chain
// ientry is the entry number in the current Tree
// Note that the argument to GetEntry must be:
// jentry for TChain::GetEntry
// ientry for TTree::GetEntry and TBranch::GetEntry
//
// To read only selected branches, Insert statements like:
// METHOD1:
// fChain->SetBranchStatus("*",0); // disable all branches
// fChain->SetBranchStatus("branchname",1); // activate branchname
// METHOD2: replace line
// fChain->GetEntry(jentry); //read all branches
//by b_branchname->GetEntry(ientry); //read only this branch
if (fChain == 0) return;
Long64_t nentries = fChain->GetEntriesFast();
Long64_t nbytes = 0, nb = 0;
TH1F *Mix=new TH1F("Mix","Mix",3000,0.5,3000.5);
TH1I *Co60=new TH1I("Co60","Co60",3000,0.5,3000.5);
TH1I *RCo60=new TH1I("RCo60","RCo60",3000,0.5,3000.5);
TH1F *Se=new TH1F("Se","Se",3000,0.5,3000.5);
TH1F *Se103=new TH1F("Se103","Se103",3600,0,60); // Se-81 emmitts 103 keV photon and half life is 57.28 min
TRandom3 R;
Float_t Energy;
Int_t IEnergy;
for (Long64_t jentry=0; jentry<nentries;jentry++) {
Long64_t ientry = LoadTree(jentry);
if (ientry < 0) break;
nb = fChain->GetEntry(jentry); nbytes += nb;
// if (Cut(ientry) < 0) continue;
// Convert Channels to keV
// Energy= (Int_t) (1.02127*evt_Chan-13.2088);
Energy= 1.02127*(evt_Chan)+(R.Rndm()-1)-13.2088; // The ADC may have a broken but so inject a half ADC c\
hannel of noise into the system I never say channel 247 plus others
IEnergy= (Int_t) Energy;
// if(Energy > 234 && Energy < 238)
// cout << "evt_Chan=" << evt_Chan<< " Energy=" << Energy << " IEnergy = " << IEnergy << e\
ndl;
// Mixture Events
if ((evt_Sec>0 && evt_Sec <300) || ( evt_Sec > 730 && evt_Sec<1020 )|| ( evt_Sec > 1480 && evt_Sec<1775\
)|| ( evt_Sec > 2250 && evt_Sec<2550 )|| ( evt_Sec > 3050 && evt_Sec<3300 )|| ( evt_Sec > 3775 && evt_Sec<4\
050 )|| ( evt_Sec > 4480 && evt_Sec<4770))
{
Mix->Fill(Energy,0.0005); // measure Mix sample for 2000 seconds
}
//Co-60 Events lines at 1173 and 1332 keV
if( (evt_Sec > 300 && evt_Sec <360 ) || ( evt_Sec> 680 && evt_Sec <710 )
|| ( evt_Sec> 1030 && evt_Sec <1080 ) || ( evt_Sec> 1400 && evt_Sec <1440 )
|| ( evt_Sec> 1800 && evt_Sec <1840 ) || ( evt_Sec> 2190 && evt_Sec <2220 )
|| ( evt_Sec> 2590 && evt_Sec <2620 ) || ( evt_Sec> 2950 && evt_Sec <3000 )
|| ( evt_Sec> 3390 && evt_Sec <3350 ) || ( evt_Sec> 3700 && evt_Sec <3750 )
|| ( evt_Sec> 4060 && evt_Sec <4100 ) || ( evt_Sec> 4420 && evt_Sec <4470 )
|| ( evt_Sec> 4790 && evt_Sec <4820 ))
{
// Co60->Fill(IEnergy); // measure Mix sample for 2000 seconds
Co60->Fill(evt_Chan); // measure Mix sample for 2000 seconds
RCo60->Fill(Energy); // measure Mix sample for 2000 seconds
}
// Pure Se events measured for 1915 seconds
if( ( evt_Sec> 400 && evt_Sec <640 ) || ( evt_Sec > 1100 && evt_Sec <1360 )
|| ( evt_Sec > 1875 && evt_Sec <2150 ) || ( evt_Sec > 2650 && evt_Sec <2930 )
|| ( evt_Sec > 3400 && evt_Sec <3690 ) || ( evt_Sec > 4120 && evt_Sec <4400 )
|| ( evt_Sec > 4840 && evt_Sec <5130))
{
Se->Fill(Energy,0.0005221932);
if(Energy>101 && Energy < 104) // half life = ln(2)/slope of exponential fit
Se103->Fill(evt_Sec/60); // half life = ln(2)/[0.01172 +/- 0.00013 ] = 59.14 +/- 0.67
// expected half life = 57.28 min
}
}
cout<< "Available Histograms"
<< " Se->Draw();" << endl
<< " Se103->Draw();" << endl
<< " Co60->Draw();" << endl
<< " RCo60->Draw();" << endl
<< " Mix->Draw();" << endl
<< endl;
}