HOME ESR EXPERIMENTS LABS etc
E121 E127 E132 E125 E143 laser_cooling_2021 E142 Ê128 E146 E0052 E0018
  nuclear two-photon decay, Page 10 of 12  Not logged in ELOG logo
ID Dateup Author Category Subject
  184   Thu Jul 1 23:07:47 2021 IK, KKRunssource
The change of the source happens now
  185   Thu Jul 1 23:38:55 2021 IK, KKRunsChange of gas bottle
The gas bottle has been changed and the beam is back.
  186   Thu Jul 1 23:46:17 2021 IrisGeneralLess beam after gas bottle change?
Looks like there is a bit less beam now in the ring. I just looked at one spectra from before the gas bottle change and after...
Attachment 1: Comparison-21h-2340h.png
Comparison-21h-2340h.png
  187   Thu Jul 1 23:54:43 2021 IrisGeneralLess beam after gas bottle change?
No, everything okay! Must have had a very good earlier run to compare.

> Looks like there is a bit less beam now in the ring. I just looked at one spectra from before the gas bottle change and after...
Attachment 1: Comparison-21h-2350h.png
Comparison-21h-2350h.png
  188   Fri Jul 2 00:08:22 2021 IK, KanikaDAQNTCAP status
The latest recorded tdms files:
 p:\E143\sc\SC_2021-06-30_23-27-45\0000338
 p:\E143\iq\IQ_2021-06-30_23-27-34\0001651
2. The rest space on the drive "p:\" : 28,4TByte
  189   Fri Jul 2 01:31:35 2021 The remote Canadian beavers (Iris, Guy)AnalysisFirst spectra after bottle change
The first good spectra after the gas bottle change is from 23:38:25 (on the 410 MHz Analyzer) 
Attachment 1: 410MHz-2021.07.01.23.38.25.809.tiq.png
410MHz-2021.07.01.23.38.25.809.tiq.png
Attachment 2: 410MHz-2021.07.01.23.38.08.817.tiq.png
410MHz-2021.07.01.23.38.08.817.tiq.png
  190   Fri Jul 2 01:34:03 2021 The remote Canadian beavers (Iris, Guy)AnalysisSummary of spectra 22:00-23:00
The last hour before the gas bottle change is summarized here. 
The first pic shows the 410 MHz spectra taken from 22-23h on the 410 MHz analyzer, in time slices of 0.1 sec from 0s-2s.
The second pic shows the spectra in 20ms slices from 1s-1.16s.

python3 e143_analyzer.py -s 1.0 -t 0.02 /home/dillmann/.../sparc-daq-02.gsi.de/E143-410MHz/410MHz-2021.07.01.22*

The third spectra was taken in May, for comparison.

I don't think we see anything...
Attachment 1: July1-22-23h-410MHz_Spectra_0.1s_steps.png
July1-22-23h-410MHz_Spectra_0.1s_steps.png
Attachment 2: July1-22-23h-410MHz_Spectra_0.02s_steps_start_1s.png
July1-22-23h-410MHz_Spectra_0.02s_steps_start_1s.png
Attachment 3: Sum_410MHz_May10-1-2am_5_0-5_175s.png
Sum_410MHz_May10-1-2am_5_0-5_175s.png
  191   Fri Jul 2 01:39:51 2021 IK, KKRunsdrop of intensity
We see a drop in beam intensity from 7E9 to 2E8 between 01:30 and 01:43
  192   Fri Jul 2 01:43:56 2021 IK, KKDAQNTCAP status
The latest recorded tdms files:
 p:\E143\sc\SC_2021-06-30_23-27-45\0000360
 p:\E143\iq\IQ_2021-06-30_23-27-34\0001759
2. The rest space on the drive "p:\" : 28TByte
  193   Fri Jul 2 02:03:36 2021 The remote Canadian beavers (Iris, Guy)AnalysisSummary of spectra 0:00-1:00
Same as before, for 0-1h on July 2nd


> The last hour before the gas bottle change is summarized here. 
> The first pic shows the 410 MHz spectra taken from 22-23h on the 410 MHz analyzer, in time slices of 0.1 sec from 0s-2s.
> The second pic shows the spectra in 20ms slices from 1s-1.16s.
> 
> python3 e143_analyzer.py -s 1.0 -t 0.02 /home/dillmann/.../sparc-daq-02.gsi.de/E143-410MHz/410MHz-2021.07.01.22*
> 
> The third spectra was taken in May, for comparison.
> 
> I don't think we see anything...
Attachment 1: July2-0-1h-410MHz_Spectra_0.02s_steps_start_1s.png
July2-0-1h-410MHz_Spectra_0.02s_steps_start_1s.png
  194   Fri Jul 2 03:10:03 2021 IK, KKRunsdrop of intensity
We see that intensity stabalized at 9E8 which is smaller then before. 
We neticed this at 2:50, keep observing.
  195   Fri Jul 2 03:19:49 2021 IK, KKRunsdrop of intensity
And the intensity is back again to 5E9 at 03:20
  196   Fri Jul 2 04:27:38 2021 The remote Canadian beavers (Iris and Guy)AnalysisSummary of spectra 0:00-4:00 (410 MHz)
Same as before for 0-4h on July 2nd.

> Same as before, for 0-1h on July 2nd
> 
> 
> > The last hour before the gas bottle change is summarized here. 
> > The first pic shows the 410 MHz spectra taken from 22-23h on the 410 MHz analyzer, in time slices of 0.1 sec from 0s-2s.
> > The second pic shows the spectra in 20ms slices from 1s-1.16s.
> > 
> > python3 e143_analyzer.py -s 1.0 -t 0.02 /home/dillmann/.../sparc-daq-02.gsi.de/E143-410MHz/410MHz-2021.07.01.22*
> > 
> > The third spectra was taken in May, for comparison.
> > 
> > I don't think we see anything...
Attachment 1: July2-0-4h-410MHz_Spectra_0.02s_steps_start_1s.png
July2-0-4h-410MHz_Spectra_0.02s_steps_start_1s.png
  197   Fri Jul 2 04:58:39 2021 IK, KanikaDAQNTCAP status
The latest recorded tdms files:
 p:\E143\sc\SC_2021-06-30_23-27-45\0000405
 p:\E143\iq\IQ_2021-06-30_23-27-34\0001978
2. The rest space on the drive "p:\" : 27.1TByte
  198   Fri Jul 2 05:08:16 2021 Guy Leckenby (Canadian beaver team)Analysise143_analyser.py Summing Script
I've written two scripts that others might find useful for the on-line analysis. 
  1. The first is a simple bash script that runs a for loop to run the e143_analyzer.py script for a sequence of start values (-s). You'll need to change the increment and 
start/end of the for loop to set the time resolution.
  2. The second is a root script that takes the above generated scripts, reads the histograms, and then fills a 2D histogram. You'll need to adjust the TH2 bin number to 
match both the bin number of the TH1s and the number of analysed files. I've included a print line to help identify the TH1 bin number.

To see the first results, I've attached a sum of July 1st 22-23h data. Sadly, I don't see any indication of the isomer, or at least the production is very low.
Attachment 1: e143_sum_ana.bash
#!/usr/bin/env bash

for i in $(seq 0 0.02 2)
do
	python e143_analyzer.py -t 0.02 -s $i -o 410MHz-07.01.22h_s${i}_t0.02 E143_Data/sparc-daq-02.gsi.de/E143-410MHz/410MHz-2021.07.01.22*.tiq
done
Attachment 2: e143_comparison.C
#define e143_comparison_cxx

#include <iostream>
#include <math.h>
#include <vector>
#include "TROOT.h"
#include "TFile.h"
#include "TChain.h"
#include "TTree.h"
#include "TBranch.h"
#include "TCanvas.h"
#include "TGraph.h"
#include "TGraphPainter.h"
#include "TColor.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TH3F.h"
#include "THStack.h"
#include "THistPainter.h"

using namespace std;

void e143_comparison()
   {
	TH2 *h1vt_410MHz_07m01d22h_10ms = new TH2F("h1vt_410MHz_07m01d22h_10ms","410Mhz 07.01.22h Additiion (10 ms res);frequency (Hz);time (s)",256,4.06704e+08,4.06753e+08,200,0,2);
	int nframes = 200;
	vector<TFile*> tf_frames;
	vector<TH1*> h1_frames(nframes,NULL);
	char chname[300];
	for(int i=0; i<nframes; i++)
	   {
		sprintf(chname,"E143_Results/410MHz-07.01.22h_s%.2f_t0.01.root",i*0.01);
		//cout << chname << endl;
		tf_frames.push_back(new TFile(chname));
		tf_frames[i]->GetObject("th1f",h1_frames[i]);
		for(long j=0; j<256; j++)
		   {
			h1vt_410MHz_07m01d22h_10ms->SetBinContent(j+1, i+1, h1_frames[i]->GetBinContent(j));
		   }
	   }
	cout << "Total Bins: " << h1_frames[100]->GetNbinsX() << ", Xlow: " << h1_frames[100]->GetBinCenter(1) << ", Xhigh: " << h1_frames[100]->GetBinCenter(h1_frames[100]->GetNbinsX()) << endl;
	
	TCanvas *c1 = new TCanvas("c1","c1",50,50,1000,800);
	h1vt_410MHz_07m01d22h_10ms->Draw("COLZ");

   	TFile *fhists = new TFile("e143_added_th2s.root","RECREATE");
	h1vt_410MHz_07m01d22h_10ms->Write();
	cout << "Written h1vt_410MHz_07m01d22h_10ms to file." << endl;
	fhists->Close();
   }


int main ()
   {
	e143_comparison();

	return 0;
   }
Attachment 3: 410MHz-07.01.22h-sum_20ms.png
410MHz-07.01.22h-sum_20ms.png
  199   Fri Jul 2 06:33:48 2021 The remote Canadian beavers (Iris and Guy)AnalysisSummary of spectra 4:00-6:00 (410 MHz)
Same as before for 4-5h and 5-6h.
No isomer, I think :-( (Sorry, Wolfram!)

> Same as before for 0-4h on July 2nd.
> 
> > Same as before, for 0-1h on July 2nd
> > 
> > 
> > > The last hour before the gas bottle change is summarized here. 
> > > The first pic shows the 410 MHz spectra taken from 22-23h on the 410 MHz analyzer, in time slices of 0.1 sec from 0s-2s.
> > > The second pic shows the spectra in 20ms slices from 1s-1.16s.
> > > 
> > > python3 e143_analyzer.py -s 1.0 -t 0.02 /home/dillmann/.../sparc-daq-02.gsi.de/E143-410MHz/410MHz-2021.07.01.22*
> > > 
> > > The third spectra was taken in May, for comparison.
> > > 
> > > I don't think we see anything...
Attachment 1: July2-4-6h-410MHz_Spectra_0.02s_steps_start_1s.png
July2-4-6h-410MHz_Spectra_0.02s_steps_start_1s.png
  200   Fri Jul 2 06:51:45 2021 Pavlos, KanikaDAQNTCAP status
The latest recorded tdms files:
 p:\E143\sc\SC_2021-06-30_23-27-45\0000431
 p:\E143\iq\IQ_2021-06-30_23-27-34\0002106
2. The rest space on the drive "p:\" : 26.6TByte
  201   Fri Jul 2 06:52:11 2021 The remote Canadian eager beavers (Iris and Guy)AnalysisSummary of spectra 0:00-6:00 (410 MHz)
Same as before but for 0-6h, 25ms time slices starting at 1.04s.

> Same as before for 4-5h and 5-6h.
> No isomer, I think :-( (Sorry, Wolfram!)
> 
> > Same as before for 0-4h on July 2nd.
> > 
> > > Same as before, for 0-1h on July 2nd
> > > 
> > > 
> > > > The last hour before the gas bottle change is summarized here. 
> > > > The first pic shows the 410 MHz spectra taken from 22-23h on the 410 MHz analyzer, in time slices of 0.1 sec from 0s-2s.
> > > > The second pic shows the spectra in 20ms slices from 1s-1.16s.
> > > > 
> > > > python3 e143_analyzer.py -s 1.0 -t 0.02 /home/dillmann/.../sparc-daq-02.gsi.de/E143-410MHz/410MHz-2021.07.01.22*
> > > > 
> > > > The third spectra was taken in May, for comparison.
> > > > 
> > > > I don't think we see anything...
Attachment 1: July2-4-0-6h-410MHz_AllSpectra_0.025s_steps_start_1.04s.png
July2-4-0-6h-410MHz_AllSpectra_0.025s_steps_start_1.04s.png
  202   Fri Jul 2 07:31:44 2021 The remote Canadian eager beavers (Iris and Guy)Analysis00h to 05h Summed Spectra
Using the scripts I posted previously, I summed 5.5 hours of data from midnight to 5:30 July 2nd CEST. 
I used 10ms time steps for better time resolution, but that did leave questionable frequency resolution. If anyone would like to see higher resolved frequency plots for the 
full 5 hours, let me know.
I've attached the spectra from 0s to 4s  for the 410MHz and 0s to 2s for the 245MHz, including zooming in on the injection timing. 

Given that we are summing 5 hours of data so the background is very smooth, I can only conclude that either the isomer doesn't exit, we cannot resolve it, or we are 
producing incredibly low amounts such that the search needs to be needle in a haystack. 
Attachment 1: 245MHz-07.02.00-05h-sum_10ms.png
245MHz-07.02.00-05h-sum_10ms.png
Attachment 2: 245MHz-07.02.00-05h-sum_10ms_zoom.png
245MHz-07.02.00-05h-sum_10ms_zoom.png
Attachment 3: 410MHz-07.02.00-05h-sum_10ms.png
410MHz-07.02.00-05h-sum_10ms.png
Attachment 4: 410MHz-07.02.00-05h-sum_10ms_zoom.png
410MHz-07.02.00-05h-sum_10ms_zoom.png
  203   Fri Jul 2 08:54:12 2021 Pavlos, AlexandreDAQNTCAP status
The latest recorded tdms files:
 p:\E143\sc\SC_2021-06-30_23-27-45\0000459
 p:\E143\iq\IQ_2021-06-30_23-27-34\0002242
2. The rest space on the drive "p:\" : 26,1TByte
ELOG V3.1.5-fc6679b