| ID |
Date |
Author |
Subject |
|
250
|
Fri Mar 31 08:23:32 2017 |
CG PV JA | Run 38 |
Date: 31/03/2017
Start: 16:21
BRIKEN run: 038
AIDA Run: R16_581
BigRIPS run: 182
DLT: 170331_1524_Mg40_037.dlt
BRIKEN Rate: 40 cps
F11 Rate: 12 cps
Date: 31/03/2017
Stop:
BRIKEN run: 38
AIDA Run: R16_611
BigRIPS run:182
Root file:170331_1222_1721_B038_Mg40.root
OBSERVATIONS: |
|
251
|
Fri Mar 31 09:23:45 2017 |
CG PV JA | Run 39 |
Date: 31/03/2017
Start: 17:22
BRIKEN run: 039
AIDA Run: R16_611
BigRIPS run: 183
DLT: 170331_1722_Mg40_039.dlt
BRIKEN Rate: 40 cps
F11 Rate: 12 cps
Date: 31/03/2017
Stop: 18:33
BRIKEN run:39
AIDA Run: R16_628
BigRIPS run:183
Root file:170331_1222_1833_B039_Mg40.root
OBSERVATIONS: Stopped to change a stripper foil |
|
157
|
Thu Nov 24 02:30:56 2016 |
CG | RIBF123R1 Local ELogs |
PDFs of the local ELog for the RIBF123R1 experiment Nov 2016 |
| Attachment 1: RIBF123_Elog_1
|
| Attachment 2: RIBF123_Elog_2
|
| Attachment 3: RIBF123_Elog_3
|
| Attachment 4: RIBF123_Elog_4
|
| Attachment 5: RIBF123_Elog_5
|
| Attachment 6: RIBF123_Elog_6
|
| Attachment 7: RIBF123_Elog_7
|
| Attachment 8: RIBF123_Elog_8
|
| Attachment 9: RIBF123_Elog_9
|
|
174
|
Sat Nov 26 16:21:08 2016 |
Briken collaboration | 41Al run 2-15 |
13:08 Beam is still OFF. Background acquisition is started.
Degrader:
3 mm Pb ONLY!
13:25 Signal cables for clover G7 are connected, now G7 spectra are visible in ONLINE.
13:35 Beam Off!
14:00 Beam ON
15:08 Beam OFF, still problems with beam!
15:11 Beam ON
15:18 Acquisition is stopped. Files
161126_1308_002.dlt to 161126_1308_015.dlt
161126_1308_002.root (online)
Rates:
F11R : 4.2 & 1.0 & 9.3 Hz
BRIKEN : 150 & 12.9 & 160 Hz |
|
324
|
Tue Jun 13 06:50:12 2017 |
Briken Colaboration | 152Eu calibration |
After experiment 152Eu calibration
Working directory: /data/1706FinalCalibrations
Configuration file: 170604_Conf_BrikenFull_AfterCrisis.xlsx
DLT file name: 170613_1451_152Eu_000.dlt
Start: 14:51
Stop: 15:24
Online rootfile:152Eu_cal_z0_1451_1524_170613.root
COMMENTS:
1. Source centered in both BRIKEN and clovers (faced towards clover D4)
2. No AIDA in
3. D4 clover pulsers are at 253, 260, 265, and 271 keV
4. G7 clover pulsers are beyond 5 MeV. No pulser for crystal G7Green |
| Attachment 1: 170604_Conf_BrikenFull_AfterCrisis.xlsx
|
|
133
|
Mon Nov 7 14:17:23 2016 |
BRIKEN collaboration | Picture of shift schedule |
|
| Attachment 1: IMG_20161107_183105.jpg
|
|
|
168
|
Fri Nov 25 23:18:04 2016 |
BRIKEN collaboration | Fallon Parasitic experiments |
===== START OF PARASITIC EXPERIMENTS ====== |
|
176
|
Sat Nov 26 20:48:16 2016 |
BRIKEN collaboration | run |
===== Start of measurements 41Al ====== |
|
178
|
Sun Nov 27 04:36:16 2016 |
BRIKEN collaboration | Run 020-023 |
======= START =======
Date 2016/11/26
Start time : 20:05
BRIKEN run : 020
Daq Conf. File : 161125Conf_BrikenFull.xlsx
DLT File name : 161126_2004_
AIDA run : R10-0
AIDA start : 20:05
BigRIPS run : 40Mg-1011
BigRIPS start : 19:30
---------------------------
Shifters: AT&AT, AE, GK
======= Setup =======
Central fragment: 41Al
--- Degraders ---
Fixed: 3mm Pb attached to beam hole in BRIKEN PE shielding
Variable:
1=(0.5mm W+1mm Al), 2=(0.3mm W+1mm Al), 3=(1mm W +1mm Al)
1. OUT
2. IN
3. IN
======= COMMENTS =======
2016/11/27
02:48 Beam Off for a while.
07:15 Two peaks appeared in the syncMonitor for Ribf-AIDA. It was after the refilling of Liquid Ni. Not clear if
the possible sync problem affected the measurements.
07:20 Beam OFF
------ Rates ------
BRIKEN : 1.90963e+02 Hz
F11R : 9.75134 Hz
Clovers : 2.03958e+02 Hz
VPLT : 1.87119e+01 Hz
VPLB : 7.76665 Hz
AIDAPL : 1.00896e+01 Hz
========= END =========
Date 2016/11/27
BRIKEN run : 023
BRIKEN stop : 06:25
ROOT online : 161126_2004_020_023.root
AIDA run : R10-211
AIDA stop : 06:25
BigRIPS RUN : 1022 |
|
179
|
Sun Nov 27 04:38:15 2016 |
BRIKEN collaboration | Run 030 |
======= START =======
Date 2016/11/27
Start time : 10:41
BRIKEN run : 030
Daq Conf. File : 161125Conf_BrikenFull.xlsx
DLT File name : 161126_1041
AIDA run : R17_0
AIDA start : 10:41
BigRIPS run : 1025
BigRIPS start : 08:57:33
---------------------------
Shifters:
======= Setup =======
Central fragment: 40Mg
--- Degraders ---
Fixed: 3mm Pb attached to beam hole in BRIKEN PE shielding
Variable:
1=(0.5mm W+1mm Al), 2=(0.3mm W+1mm Al), 3=(1mm W +1mm Al)
1. IN
2. IN
3. OUT
======= COMMENTS =======
13:00 The syncro was lost. The TimeStamp was reset. We did not stop recording data. AIDA was writing
file R18_28, BRIKEN was writing file 161126_1041_030.dlt (@637MB). Notice that BRIKEN online histograms were not
reset, and the rates will not be plane.
------ Rates ------
BRIKEN : 200 Hz
F11R : 100 Hz
Clovers : 275 Hz
VPLT : 18 Hz
VPLB : 8 Hz
AIDAPL : 10 Hz
========= END =========
Date 2016/11/26
BRIKEN run : 030
BRIKEN stop : 13:26
ROOT online : 161126_1326_029_030.root
AIDA run : R28_33
AIDA stop : 13:26
BigRIPS RUN : 1028 |
|
186
|
Mon Nov 28 01:22:50 2016 |
BRIKEN collaboration | Elog RIBF123 & BigRIPS |
Full RIB123 Elog entries for Takechi && Fallon parasitic experiments
Update:03Dec16
Instructions:
1) Download RIBF123FullEntries.tar.gz-part-aa and RIBF123FullEntries.tar.gz-part-ab
2) Join and uncompress the file on the terminal:
# cat RIBF123FullEntries.tar.gz-part-*| tar xz
3) cd RIBF123FullEntries
4) Open with a Browser ELOGRIBF123_161203_2030_fullWeb.html
-----------------------------------------------------------
HTML elogs for BigRIPS are found by unziping BigRIPSRunSummary.zip and EXPDBRunSummary.zip
-----------------------------------------------------------
The full elog is also included in pdf file. See:
ElogRIBF123_161203_2030_Full_part1.pdf
ElogRIBF123_161203_2030_Full_part2.pdf |
| Attachment 1: BigRIPSRunSummary.zip
|
| Attachment 2: EXPDBRunSummary.zip
|
| Attachment 3: RIBF123FullEntries.tar.gz-part-aa
|
| Attachment 4: RIBF123FullEntries.tar.gz-part-ab
|
| Attachment 5: ElogRIBF123_161203_2030_Full_part1.pdf
|
|
| Attachment 6: ElogRIBF123_161203_2030_Full_part2.pdf
|
|
|
328
|
Thu Jun 15 03:40:17 2017 |
BRIKEN collaboration | Z-scaning 252Cf Jun17 |
A scaning with 252Cf along the Z-axis of the neutron detector is carried out. The setup is similar to the scaning carried out in October 2016.
Observations:
- Clovers are put inside the detector without voltage.
- AIDA detector is out of the moderator.
- Position in Z is given with respect to the moderator side opposite to AIDA.
- The shielding for in-beam neutrons is mounted. It corresponds to HDPE 2.5cm + 1mm Cd + 20cm HDPE. The distance moderator to shieldind is 34.6mm
- There are two additional neutron sources in the Hall, aprox 10m from the detector. One 241AmBe (30mC) and 244Cm/O13 (500mC). Both sources produce an increase of the neutron background up to approx1n/s (x10 higher).
- Datafiles are located in the dac computer, path:
/data/1706FinalCalibrations/170615_Zscanning
MEASUREMENTS
Configuration file: 170615_Conf_BrikenFull_AfterCrisis.xlsx
0) A background measurement is carried out previous to the Zscan.
FILE: 170615_1100_Bkg.root
1) Z=25mm; Start 11:13 Stop 11:23
FILES: 170615_1113_Zscan_252Cf_25mm_000.dlt, 170615_1113_1123_Zscan_252Cf_25mm.root
2) Z=50mm; Start 11:25 Stop 11:34
FILES: 170615_1125_Zscan_252Cf_50mm_001.dlt, 170615_11125_1134_Zscan_252Cf_50mm.root
3) Z=75mm; Start 11:37 Stop 11:46
FILES: 170615_1137_Zscan_252Cf_75mm_002.dlt, 170615_1137_1146_Zscan_252Cf_75mm.root
4) Z=100mm; Start 13:25 Stop 13:33
FILES: 170615_1325_Zscan_252Cf_100mm_003.dlt 170615_1325_1333_Zscan_252Cf_100mm.root
5) Z=125mm; Start 13:36 Stop 13:44
FILES: 170615_1335_Zscan_252Cf_125mm_004.dlt 170615_1335_1344_Zscan_252Cf_125mm.root
6) Z=150mm; Start 13:51 Stop 13:59
FILES: 170615_1351_Zscan_252Cf_150mm_005.dlt 170615_1351_1359_Zscan_252Cf_150mm.root
7) Z=175mm; Start 14:01 Stop 14:09
FILES: 170615_1401_Zscan_252Cf_175mm_006.dlt 170615_1401_1409_Zscan_252Cf_175mm.root
8) Z=200mm; Start 14:12 Stop 14:21
FILES: 170615_1412_Zscan_252Cf_200mm_007.dlt 170615_1412_1421_Zscan_252Cf_200mm.root
9) Z=225mm; Start 14:24 Stop 14:31
FILES: 170615_1424_Zscan_252Cf_225mm_008.dlt 170615_1424_1431_Zscan_252Cf_225mm.root
10) Z=250mm; Start 14:33 Stop 14:40
FILES: 170615_1433_Zscan_252Cf_250mm_009.dlt 170615_1433_1440_Zscan_252Cf_250mm.root
11) Z=275mm; Start 14:43 Stop 14:50
FILES: 170615_1443_Zscan_252Cf_275mm_010.dlt 170615_1443_1450_Zscan_252Cf_275mm.root
12) Z=300mm; Start 14:54 Stop 15:01
FILES: 170615_1454_Zscan_252Cf_300mm_011.dlt 170615_1454_1501_Zscan_252Cf_300mm.root
13) Z=325mm; Start 15:04 Stop 15:11
FILES: 170615_1504_Zscan_252Cf_325mm_012.dlt 170615_1504_1511_Zscan_252Cf_325mm.root
14) Z=350mm; Start 15:14 Stop 15:22
FILES: 170615_1514_Zscan_252Cf_350mm_013.dlt 170615_1514_1522_Zscan_252Cf_350mm.root
15) Z=375mm; Start 15:24 Stop 15:31
FILES: 170615_1524_Zscan_252Cf_375mm_014.dlt 170615_1524_1531_Zscan_252Cf_375mm.root
16) Z=400mm; Start 15:34 Stop 15:41
FILES: 170615_1533_Zscan_252Cf_400mm_015.dlt 170615_1533_1541_Zscan_252Cf_400mm.root
17) Z=425mm; Start 15:44 Stop 15:57
FILES: 170615_1543_Zscan_252Cf_425mm_016.dlt 170615_1543_1557_Zscan_252Cf_425mm.root
18) Z=450mm; Start 15:59 Stop 16:08
FILES: 170615_1559_Zscan_252Cf_450mm_017.dlt 170615_1559_Zscan_252Cf_450mm_017.dlt |
| Attachment 1: WhatsApp_Image_2017-06-15_at_18.11.52.jpeg
|
|
|
329
|
Thu Jun 15 10:07:56 2017 |
BRIKEN collaboration | Gamma calibrations |
Clover calibrations.
- Datafiles are located in the dac computer, path:
/data/1706FinalCalibrations/GammaCalibration
MEASUREMENTS
Configuration file: 170615_Conf_BrikenFull_AfterCrisis.xlsx
1) 152Eu. Source at the center of the moderator (0,0,0)
Source ID: 09 7011
Files: 170615_1805_152Eu_SourceNum272_000.dlt 170615_1805_1818_152Eu_SourceNum272.root
2) 137Cs. Source at the center of the moderator (0,0,0)
Source ID: 4763
Files: 170615_1837_137Cs_SourceNum4763_001.dlt 170615_1837_1846_137Cs_SourceNum4763.root
3) 152Eu. Source at the center of the moderator (0,0,0)????
Measurement over night.
Source ID: 09 7011
Files: 170615_1805_1850_152Eu_SourceNum272_002.dlt |
| Attachment 1: WhatsApp_Image_2017-06-15_at_18.00.36.jpeg
|
|
| Attachment 2: WhatsApp_Image_2017-06-15_at_18.00.51.jpeg
|
|
| Attachment 3: WhatsApp_Image_2017-06-15_at_18.37.12.jpeg
|
|
|
330
|
Mon Oct 16 05:52:14 2017 |
BRIKEN collaboration | START preparation |
We start to check every thing for the coming run
New directory: /data/1710Preparation/
Copy last configuration file of RIBF148 experiment
renamed:171016_Conf_BrikenFull.xlsx
Edinburgh clock parameters:
-Square signal.
-Amplitude: 1 Vpp
-Offset: 0.5V
-FreqCh1=FreqCh2
Status:
-F11 (high/low gain), AIDA_PL, DT_Pulser, Sync_Pulse, He tubes are defined and connected. Preliminary
configuration "171016_Conf_BrikenFull.xlsx" attached.
-YSO & clovers are not connected yet. |
| Attachment 1: 171016_Conf_BrikenFull.xlsx
|
|
355
|
Thu Nov 16 02:32:23 2017 |
BRIKEN collaboration | YSO thresholds |
YSO Tests
Thresholds & calibrations
For HV between (500,600)V, the rate of visual scalers do not change. The current thresholds,
Ch Thr
xa 130
xb 130
ya 130
yb 150
dy 90
Attachment 1: dynVSdynGated.
>>Black, dyn signal
>>Red, dyn signal gated by 4 anodes (histo xpos) |
| Attachment 1: dynVSdynGated.png
|
|
|
185
|
Sun Nov 27 23:14:49 2016 |
BRIKEN coll | run 036 |
======= START =======
Date 2016/11/27
Start time : 19:28
BRIKEN run : 36
Daq Conf. File : 161125Conf_BrikenFull.xlsx
DLT File name : 161128_1928
AIDA run : R24
AIDA start : 19:28
BigRIPS run : 1034
BigRIPS start : 18:33
---------------------------
Shifters: PV, KG, TD
======= Setup =======
Central fragment: 41Al @F8
--- Degraders ---
Fixed: 3mm Pb between Veto plastics up and down + 1mm Al
Variable:
1=(0.5mm W+1mm Al), 2=(0.3mm W+1mm Al), 3=(1mm W +1mm Al)
1. IN
2. IN
3. OUT
======= COMMENTS =======
TimeStamp has been reset two times during this run.
ERROR! Data acquisition with errors during this ring because not restarted after Liquid Nitrogen refilling.
------ Rates ------
BRIKEN : Hz
F11R : Hz
Clovers : Hz
VPLT : Hz
VPLB : Hz
AIDAPL : Hz
========= END =========
Date 2016/11/27
BRIKEN run : 036
BRIKEN stop : 20:38
ROOT online : 161127_1928_036.root
AIDA run : R25_0
AIDA stop :
BigRIPS RUN : 1035 |
|
180
|
Sun Nov 27 07:10:29 2016 |
BRIKEN colaboration | Run 31 |
======= START =======
Date 2016/11/27
Start time : 13:53
BRIKEN run : 031
Daq Conf. File : 161125Conf_BrikenFull.xlsx
DLT File name : 161127_1353
AIDA run : R20_0
AIDA start : 13:51
BigRIPS run : 1028
BigRIPS start : 12:32:36
---------------------------
Shifters:
======= Setup =======
Central fragment: 40Mg
--- Degraders ---
Fixed: 3mm Pb attached to beam hole in BRIKEN PE shielding
Variable:
1=(0.5mm W+1mm Al), 2=(0.3mm W+1mm Al), 3=(1mm W +1mm Al)
1. IN
2. IN
3. OUT
======= COMMENTS =======
------ Rates ------
BRIKEN : Hz
F11R : Hz
Clovers : Hz
VPLT : Hz
VPLB :
AIDAPL : Hz
========= END =========
Date 2016/11/26
BRIKEN run : 031
BRIKEN stop : 16:04
ROOT online : 161127_1353_031.root
AIDA run : R20_27
AIDA stop : 16:04
BigRIPS RUN : 1030 |
|
150
|
Fri Nov 11 08:07:18 2016 |
BCR SG NTB | BRIKEN Hit Visualization |
These are the files used to create and visualize the hit patterns of the BRIKEN detector.
There are three files needed plus one input data file. The three files are "HybridBRIKENPhysicalHistogram.hpp", "HybridBRIKENPhysicalHistogram.cpp", and "HybryidBRIKEN_Tube_Info_1.txt". These files are attached to this eLog entry.
In order to run this type in a terminal
root -l HybidBRIKENPhysicalHistogram.cpp+
The plus is important, it tells ROOT to compile this with the gcc compiler and then run that compiled version of the code. It runs slow, if at all, in the ROOT C++ interpreter CINT.
The code has been tested with ROOT v5.32 and v5.34, not version 6.
At the bottom of the HybidBRIKENPhysicalHistogram.cpp file is the "main program". In it you can set the input data file, which should be a BRIKEN online root output file, and the 3He energy cuts. It is easiest if all the files are in the same folder. If there are any questions please feel free to email B. Charles (Charlie) Rasco. |
| Attachment 1: HybridBRIKENPhysicalHistogram.cpp
|
//
// HybridBRIKENPhysicalHistogram.cpp
//
//
// Created by Bertis Charles Rasco on 11/4/16.
//
//
#include "HybridBRIKENPhysicalHistogram.hpp"
#include <cmath>
#include <fstream>
#include <iostream>
#include <sstream>
#include <string>
#include <stdio.h>
//#include <stdlib.h>
#include "TH1F.h"
#include "TRandom3.h"
#include "TFile.h"
// ROOT Needs this for some reason.
HybridBRIKENHistogram::HybridBRIKENHistogram(): TH2Poly()
{
}
HybridBRIKENHistogram::HybridBRIKENHistogram( const char *name, const char *title, Double_t xlow, Double_t xup, Double_t ylow, Double_t yup ):
TH2Poly( name, title, xlow, xup, ylow, yup )
{
this->ReadinTubes();
this->FillTubeCountsRandom();
// this->SetupTubes();
// this->Draw();
}
HybridBRIKENHistogram::HybridBRIKENHistogram( const char *name, const char *title, Double_t xlow, Double_t xup, Double_t ylow, Double_t yup, TString inputFileName ):
TH2Poly( name, title, xlow, xup, ylow, yup )
{
m_InputFileName = inputFileName;
this->ReadinTubes();
this->ReadInDataFromFile();
// this->SetupTubes();
// this->Draw();
}
HybridBRIKENHistogram::~HybridBRIKENHistogram()
{
}
void HybridBRIKENHistogram::ReadinTubes( void )
{
std::ifstream file("HybridBRIKEN_Tube_Info_1.txt");
std::string str;
Int_t lineCount = 0;
Int_t tubeCount = 0;
Double_t rad;
Double_t pos[2];
Int_t tubeID;
Int_t ringID;
Char_t tubeType[20];
// std::stringstream ss;
// std::string::size_type sz;
while ( std::getline(file, str) )
{
// Process str
std::stringstream ss( str);
if( lineCount > 0 )
{
tubeType[0] = 'X';
ss >> tubeID >> tubeType >> rad >> pos[0] >> pos[1] >> ringID;
rad *= 0.5;// 2R is what is in file.
// tubeID = std::stoi( str, &sz );
if( tubeID > 0 && (tubeType[0]=='R' || tubeType[0]=='U' || tubeType[0]=='O' ) )
{
if( tubeType[0] != 'O' ) m_Tubes.push_back( new HeTube( rad, pos[0], pos[1], tubeID, ringID, tubeType ) );
else if( tubeType[0] == 'O' )
{
tubeType[1] = tubeType[4];
m_Tubes.push_back( new HeTube( rad, pos[0], pos[1], tubeID, ringID, tubeType ) );
}
tubeCount++;
}
}
// if( lineCount < 11 )
// {
// std::cout << str << std::endl;
// std::cout << tubeID << ", " << tubeType[0] << tubeType[1] << ", " << rad << ", " << pos[0] << ", " << pos[1] << ", " << ringID << std::endl;
// ss >> tubeID >> tubeType >> rad >> pos[0] >> pos[1] >> ringID;
// }
lineCount++;
}
Int_t binID;
// Double_t minX = m_Tubes[0]->pos[0] - m_Tubes[0]->radius;
// Double_t maxX = m_Tubes[0]->pos[0] + m_Tubes[0]->radius;
// Double_t minY = m_Tubes[0]->pos[1] - m_Tubes[0]->radius;
// Double_t maxY = m_Tubes[0]->pos[1] + m_Tubes[0]->radius;
// Fix
// binID = this->AddBin( minX, minY, maxX, maxY );
binID = this->AddCircularBin( 0 );
// std::cout << "X: " << this->GetNbinsX() << std::endl;
// std::cout << "Y: " << this->GetNbinsY() << std::endl;
// std::cout << "Z: " << this->GetNbinsZ() << std::endl;
// std::cout << "BinID: " << binID << std::endl;
// this->SetBinContent( binID, 9);
for( Int_t i = 1; i < (Int_t) m_Tubes.size(); i++ )
{
// minX = m_Tubes[i]->pos[0] - m_Tubes[i]->radius;
// maxX = m_Tubes[i]->pos[0] + m_Tubes[i]->radius;
// minY = m_Tubes[i]->pos[1] - m_Tubes[i]->radius;
// maxY = m_Tubes[i]->pos[1] + m_Tubes[i]->radius;
// Fix AddCircularBin
// binID = this->AddBin( minX, minY, maxX, maxY );
binID = this->AddCircularBin( i );
// std::cout << "BinID: " << binID << std::endl;
// this->SetBinContent( binID, 15 + binID * 11);
}
std::cout << "Number of Lines: " << lineCount << std::endl;// = How many tubes in total?
std::cout << "Number of Tubes: " << binID << std::endl;// = How many tubes in total?
// std::cout << "Number of Bins: " << this->GetNcells() << std::endl;// = How many tubes in total?
}
void HybridBRIKENHistogram::SetupTubes( void )
{
// 1" tubes
const Int_t maxN1 = 8;
for( Int_t i = 0; i < maxN1; i++ )
{
m_Tubes.push_back( new HeTube(0.5, 2.0 * cos( i * 6.28 / maxN1 ), 2.0 * sin(i * 6.28 / maxN1 ), i, 1, "RI" ) );
}
// 2" tube at 14" radius
const Int_t maxN = 3;
for( Int_t i = 0; i < maxN; i++ )
{
m_Tubes.push_back( new HeTube(1.0, 14.0 * cos( i * 6.28 / maxN ), 14.0 * sin(i * 6.28 / maxN ), i, 2, "O2" ) );
}
Int_t binID;
Double_t minX = m_Tubes[0]->pos[0] - m_Tubes[0]->radius;
Double_t maxX = m_Tubes[0]->pos[0] + m_Tubes[0]->radius;
Double_t minY = m_Tubes[0]->pos[1] - m_Tubes[0]->radius;
Double_t maxY = m_Tubes[0]->pos[1] + m_Tubes[0]->radius;
// Fix this
binID = this->AddBin( minX, minY, maxX, maxY );
// std::cout << binID << std::endl;
// this->SetBinContent( binID, 9);
for( Int_t i = 0; i < (Int_t) m_Tubes.size(); i++ )
{
minX = m_Tubes[i]->pos[0] - m_Tubes[i]->radius;
maxX = m_Tubes[i]->pos[0] + m_Tubes[i]->radius;
minY = m_Tubes[i]->pos[1] - m_Tubes[i]->radius;
maxY = m_Tubes[i]->pos[1] + m_Tubes[i]->radius;
// Fix this
binID = this->AddBin( minX, minY, maxX, maxY );
// std::cout << binID << std::endl;
// this->SetBinContent( binID, 15 + binID * 11);
}
}
void HybridBRIKENHistogram::ReadInDataFromFile( void )
{
const Int_t lowHeCut = 600;//20;// Cut is by bin not by energy. But for 1 keV per bin they are the same.
const Int_t highHeCut = 900;//200;// Cut is by bin not by energy. But for 1 keV per bin they are the same.
m_InputFile = new TFile( m_InputFileName );
Int_t binID;
Double_t counts;
char tempCh[20];
for( Int_t i = 0; i < m_Tubes.size(); i++ )
{
binID = i + 1;
TString histName("He");
if( binID < 10 )
{
histName.Append("00");
}
else if( binID < 100 )
{
histName.Append("0");
}
// else if( binID < 1000 )
// {
// }
// histName.Append( TString::Itoa( (Int_t) binID, 10 ) );// Works for ROOT 5.34, but not for 5.32
// itoa(binID,tempCh,10);// Does not work for 5.32 unkown if works for 5.34
sprintf (tempCh, "%d", binID);// works for 5.32
histName.Append( TString( tempCh ) );
// std::cout << histName << std::endl;
TH1F *histogram = (TH1F*)m_InputFile->Get( histName.Data() );
counts = histogram->Integral( lowHeCut, highHeCut );
this->SetBinContent( binID, counts);
}
}
void HybridBRIKENHistogram::FillTubeCountsRandom( void )
{
TRandom3 rand;
// for( Int_t i = 0; i < this->GetNbins() ; i++)
// {
// hist->Fill( rand.Uniform(1.0,140.0) );
// }
for( Int_t i = 1; i <= (Int_t) m_Tubes.size(); i++ )
{
this->SetBinContent(i, rand.Gaus(1000., 50.0) );
}
}
Int_t HybridBRIKENHistogram::AddCircularBin( Int_t tubeID )
{
Int_t binID;
Double_t r = m_Tubes[tubeID]->radius;
Double_t theta;
const Int_t N = 16;
Double_t* x;
Double_t* y;
x = new Double_t[N];
y = new Double_t[N];
for( Int_t i = 0; i < N; i++ )
{
theta = i * 2.0 * 3.1415 / N;
x[i] = m_Tubes[tubeID]->pos[0] + r * cos( theta );
y[i] = m_Tubes[tubeID]->pos[1] + r * sin( theta );
}
binID = this->AddBin( N, x, y );
// Double_t minX = m_Tubes[tubeID]->pos[0] - m_Tubes[tubeID]->radius;
// Double_t maxX = m_Tubes[tubeID]->pos[0] + m_Tubes[tubeID]->radius;
// Double_t minY = m_Tubes[tubeID]->pos[1] - m_Tubes[tubeID]->radius;
// Double_t maxY = m_Tubes[tubeID]->pos[1] + m_Tubes[tubeID]->radius;
//
////
// binID = this->AddBin( minX, minY, maxX, maxY );
return binID;
}
//***************************************************************************************
// Main Program...
void HybridBRIKENPhysicalHistogram( void )
{
Double_t halfHDPEWidth = 500.0;// mm
// TString inputFileName( "161103_1426_152Eu.root" );
// TString inputFileName( "161103_1319_252Cf.root" );
// TString inputFileName( "161107_0833_Overnight.root" );
// TString inputFileName( "161108_1609_ShortBackground.root" );
// TString inputFileName( "161013_1518_137Cs60Co.root" );
TString inputFileName( "161108_2323_717273Ni.root" );
HybridBRIKENHistogram *hist = new HybridBRIKENHistogram("hBR", "HybridBRIKEN", -halfHDPEWidth, halfHDPEWidth, -halfHDPEWidth, halfHDPEWidth, inputFileName );
// HybridBRIKENHistogram *hist = new HybridBRIKENHistogram("hBR", "HybridBRIKEN", -halfHDPEWidth, halfHDPEWidth, -halfHDPEWidth, halfHDPEWidth );
hist->Draw("COLZ");
}
|
| Attachment 2: HybridBRIKENPhysicalHistogram.hpp
|
//
// HybridBRIKENPhysicalHistogram.hpp
//
//
// Created by Bertis Charles Rasco on 11/4/16.
//
//
#ifndef HybridBRIKENPhysicalHistogram_hpp
#define HybridBRIKENPhysicalHistogram_hpp
#include <vector>
#include "TFile.h"
#include "TH2Poly.h"
#include "TObject.h"
#include "TString.h"
//class TFile;
class HeTube
{
public:
HeTube( Double_t r, Double_t x, Double_t y, Int_t tID, Int_t rID, char tubeTypeChars[2] )
{ radius = r; pos[0] = x; pos[1] = y; tubeID = tID; ringID = rID;
tubeType[0] = tubeTypeChars[0]; tubeType[1] = tubeTypeChars[1]; }
//private:
Double_t radius;
Double_t pos[2];
Int_t tubeID;
Int_t ringID;
Char_t tubeType[2]; // RI (RIKEN), UP (UPC), O1 (ORNL1), O2 (ORNL2)
};
class HybridBRIKENHistogram : public TH2Poly
{
public:
HybridBRIKENHistogram();// ROOT needs this to be defined if it is inherited from a ROOT class.
HybridBRIKENHistogram( const char *name, const char *title, Double_t xlow, Double_t xup, Double_t ylow, Double_t yup );
HybridBRIKENHistogram( const char *name, const char *title, Double_t xlow, Double_t xup, Double_t ylow, Double_t yup,
TString inputFileName );
~HybridBRIKENHistogram();
ClassDef(HybridBRIKENHistogram,1);// ROOT also needs this if it is inherited from a ROOT class. Not sure what the 1 means.
private:
void ReadinTubes( void );
void SetupTubes( void );
void ReadInDataFromFile( void );
void FillTubeCountsRandom( void );
Int_t AddCircularBin( Int_t tubeID );
TString m_InputFileName;
TFile *m_InputFile;
std::vector<HeTube*> m_Tubes;
};
#endif /* HybridBRIKENPhysicalHistogram_hpp */
|
| Attachment 3: HybridBRIKEN_Tube_Info_1.txt
|
ID Type 2R(mm) X(mm) Y(mm) RING
1 RIKEN 27.5 -77.45 -39.25 1
2 RIKEN 27.5 -77.45 0 1
3 RIKEN 27.5 -77.45 39.25 1
4 RIKEN 27.5 -106.052 -19.625 2
5 RIKEN 27.5 -106.052 19.625 2
6 RIKEN 27.5 -134.655 0 3
7 RIKEN 27.5 77.45 -39.25 1
8 RIKEN 27.5 77.45 0 1
9 RIKEN 27.5 77.45 39.25 1
10 RIKEN 27.5 106.052 -19.625 2
11 RIKEN 27.5 106.052 19.625 2
12 RIKEN 27.5 134.655 0 3
13 RIKEN 27.5 -177.45 -39.25 1
14 RIKEN 27.5 -177.45 0 1
15 RIKEN 27.5 -177.45 39.25 1
16 RIKEN 27.5 -206.052 -19.625 2
17 RIKEN 27.5 -206.052 19.625 2
18 RIKEN 27.5 -234.655 0 3
19 RIKEN 27.5 177.45 -39.25 1
20 RIKEN 27.5 177.45 0 1
21 RIKEN 27.5 177.45 39.25 1
22 RIKEN 27.5 206.052 -19.625 2
23 RIKEN 27.5 206.052 19.625 2
24 RIKEN 27.5 234.655 0 3
25 UPC 27.5 -118.3 77.45 3
26 UPC 27.5 -69.7 77.45 1
27 UPC 27.5 -34.85 77.45 1
28 UPC 27.5 0 77.45 1
29 UPC 27.5 34.85 77.45 1
30 UPC 27.5 69.7 77.45 1
31 UPC 27.5 118.3 77.45 3
32 UPC 27.5 -97.2 116.7 3
33 UPC 27.5 -60.362 112.3 2
34 UPC 27.5 -17.425 107.631 2
35 UPC 27.5 17.425 107.631 2
36 UPC 27.5 60.362 112.3 2
37 UPC 27.5 97.2 116.7 3
38 UPC 27.5 -73.928 151.378 3
39 UPC 27.5 -34.85 137.812 3
40 UPC 27.5 0 147.15 3
41 UPC 27.5 34.85 137.812 3
42 UPC 27.5 73.928 151.378 3
43 UPC 27.5 -27.0595 178.438 3
44 UPC 27.5 27.0595 178.438 3
45 UPC 27.5 -118.3 -77.45 3
46 UPC 27.5 -69.7 -77.45 1
47 UPC 27.5 -34.85 -77.45 1
48 UPC 27.5 0 -77.45 1
49 UPC 27.5 34.85 -77.45 1
50 UPC 27.5 69.7 -77.45 1
51 UPC 27.5 118.3 -77.45 3
52 UPC 27.5 -97.2 -116.7 3
53 UPC 27.5 -60.362 -112.3 2
54 UPC 27.5 -17.425 -107.631 2
55 UPC 27.5 17.425 -107.631 2
56 UPC 27.5 60.362 -112.3 2
57 UPC 27.5 97.2 -116.7 3
58 UPC 27.5 -73.928 -151.378 3
59 UPC 27.5 -34.85 -137.812 3
60 UPC 27.5 0 -147.15 3
61 UPC 27.5 34.85 -137.812 3
62 UPC 27.5 73.928 -151.378 3
63 UPC 27.5 -27.0595 -178.438 3
64 UPC 27.5 27.0595 -178.438 3
65 ORNL1 27.5 -215.5 77.45 5
66 ORNL1 27.5 -166.9 77.45 4
67 ORNL1 27.5 166.9 77.45 4
68 ORNL1 27.5 215.5 77.45 5
69 ORNL1 27.5 -301.012 214.133 7
70 ORNL1 27.5 -252.393 235.169 7
71 ORNL1 27.5 252.393 235.169 7
72 ORNL1 27.5 301.012 214.133 7
73 ORNL1 27.5 -215.5 -77.45 5
74 ORNL1 27.5 -166.9 -77.45 4
75 ORNL1 27.5 166.9 -77.45 4
76 ORNL1 27.5 215.5 -77.45 5
77 ORNL1 27.5 -301.012 -214.133 7
78 ORNL1 27.5 -252.393 -235.169 7
79 ORNL1 27.5 252.393 -235.169 7
80 ORNL1 27.5 301.012 -214.133 7
81 ORNL2 53 -139.932 142.007 4
82 ORNL2 53 -70.014 198.718 4
83 ORNL2 53 0 217.478 4
84 ORNL2 53 70.014 198.718 4
85 ORNL2 53 139.932 142.007 4
86 ORNL2 53 -214.903 127.3 5
87 ORNL2 53 -133.044 210.494 5
88 ORNL2 53 -48.6976 259.192 5
89 ORNL2 53 48.6976 259.192 5
90 ORNL2 53 133.044 210.494 5
91 ORNL2 53 214.903 127.3 5
92 ORNL2 53 -273.203 102.2 6
93 ORNL2 53 -263.894 172.91 6
94 ORNL2 53 -204.385 193.709 6
95 ORNL2 53 -118.139 282.073 6
96 ORNL2 53 0 313.728 6
97 ORNL2 53 118.139 282.073 6
98 ORNL2 53 204.385 193.709 6
99 ORNL2 53 263.894 172.91 6
100 ORNL2 53 273.203 102.2 6
101 ORNL2 53 -339.453 102.2 7
102 ORNL2 53 -323.857 159.305 7
103 ORNL2 53 -201.103 278.553 7
104 ORNL2 53 -166.264 365.428 7
105 ORNL2 53 -73.6089 352.162 7
106 ORNL2 53 73.6089 352.162 7
107 ORNL2 53 166.264 365.428 7
108 ORNL2 53 201.103 278.553 7
109 ORNL2 53 323.857 159.305 7
110 ORNL2 53 339.453 102.2 7
111 ORNL2 53 -139.932 -142.007 4
112 ORNL2 53 -70.014 -198.718 4
113 ORNL2 53 0 -217.478 4
114 ORNL2 53 70.014 -198.718 4
115 ORNL2 53 139.932 -142.007 4
116 ORNL2 53 -214.903 -127.3 5
117 ORNL2 53 -133.044 -210.494 5
118 ORNL2 53 -48.6976 -259.192 5
119 ORNL2 53 48.6976 -259.192 5
120 ORNL2 53 133.044 -210.494 5
121 ORNL2 53 214.903 -127.3 5
122 ORNL2 53 -273.203 -102.2 6
123 ORNL2 53 -263.894 -172.91 6
124 ORNL2 53 -204.385 -193.709 6
125 ORNL2 53 -118.139 -282.073 6
126 ORNL2 53 0 -313.728 6
127 ORNL2 53 118.139 -282.073 6
128 ORNL2 53 204.385 -193.709 6
129 ORNL2 53 263.894 -172.91 6
130 ORNL2 53 273.203 -102.2 6
131 ORNL2 53 -339.453 -102.2 7
132 ORNL2 53 -323.857 -159.305 7
133 ORNL2 53 -201.103 -278.553 7
134 ORNL2 53 -166.264 -365.428 7
135 ORNL2 53 -73.6089 -352.162 7
136 ORNL2 53 73.6089 -352.162 7
137 ORNL2 53 166.264 -365.428 7
138 ORNL2 53 201.103 -278.553 7
139 ORNL2 53 323.857 -159.305 7
140 ORNL2 53 339.453 -102.2 7
|
|
352
|
Thu Oct 26 09:01:53 2017 |
BCR | Number of Counts vs Expected Efficiency for Proposals |
Here is a doc with an estimate of the statistical uncertainty P1n and P2n as a function of the number of implant-beta correlations for future experiments.
This will hopefully be useful for the new PAC proposals as an estimate for the time needed to measure a certain precision.
If there are any questions please feel free to email me (Charlie Rasco).
|
| Attachment 1: ExpectedBRIKEN_Precision.pdf
|
|
|
333
|
Fri Oct 20 08:27:05 2017 |
BC Rasco | Analysis Slides |
Here are the ORNL/UTK BRIKEN analysis slides presented by BC Rasco on 2017/10/20.
Enjoy the fine reading.
The online seems to be missing some slides, but the pdf has all of the slides. |
| Attachment 1: PreliminaryPnCalcsWithBkgCoin_BRIKENGroup.pdf
|
|