ELOG MUGAST

Here is a procedure to convert a run and add data the Kinematic plot.

Every 4Go there will be a file change, convert the file and add it to the chain using the following step:

Step O: Take your favorite ortholinear keyboard

Step 1: Login in the right environment

open a terminal from one of the mugast computer
type lpcgrit to login to lpcgrit and source the right environmt
you should see red lpcgrit invite of command

Step 2: Convert the file using the offline program
This Convert the file from the ganil .dat file and contained only MUGAST and CATS stuff

type npp MUGAST/offline in the terminal to reach the correct folder
./Offline run_com2019/run_0ZZZ-DD-MM-YY.HHhMMmSSs root_com2019/run_ZZZ_0.root
OR
./Offline run_com2019/run_0ZZZ-DD-MM-YY.HHhMMmSSs.R root_com2019/run_ZZZ_R.root

Step 3: Analyse the file using the npanalysis program
This could be used with any file, either from Offline conversion, or the VAMOS Analysis conversion or the Merged Tree from AGATA.

type npp MUGAST in the terminal to reach the correct folder
there is an error below caused by the ortholinear keyboard, I guess, corrected using a normal keyboard:
npanalysis -D mugast.detector -C Calibration.txt -E ./Reaction/16Odp.reaction -T offline/root_com2019/run_0253.root AutoTree -O PhyZZZ.root
The run will be created in the NPOutputs/ folder

Step 4: Add the file to the KineMacro

vim macros/Kine.cxx
Add the new file in the LoadChain method
root macros/Kine.cxx
Enjoy the soothing feeling of a smooth knematic line during a late night shift.

Fixed

com2019

Sat Apr 6 00:29:14 2019

Experimental Kinematic Line from dp + gamma ?????

Triangle is 870 keV

Attachment 1: 20.png

Fixed

com2019

Sat Apr 6 01:06:18 2019

Experimental Kinematic Line from dp

Kinematic Line obtained from run 256+253

Attachment 1: Kine256_0-253.pdf

Fixed

com2019

Sat Apr 6 08:23:09 2019

Kinematic Lines run 256 from 0 to 12 at 8 am

Kinematic lines (E_lab vs Theta_lab) for all angular range (first picture) and zoomed plot for backward angles (second picture).

Attachment 1: Kine256_0-256_12.pdf

Attachment 2: KinematicsLineBackwardAngles.png

Fixed

com2019

Sat Apr 6 08:38:51 2019

We see d,3he more than d,t at forward angle

We see more (d.3he) than dt

there seems to be another reaction going up to 40deg Lab

Attachment 1: Kine256_0-256_3.pdf

Fixed

com2019

Sat Apr 6 09:52:47 2019

*************************
Run : 257
*************************

Start at 2019-04-06 09:52:28.767062

Run Conditions:

Operator: A. Lemasson
Trigger: MUST2 + MG BWd + VAMOS % 100
Target/Source: CD2 1mg
Beam: 16O
Intensity: ~ 40 kHz

Comment:
Second Run

Stop at 2019-04-06 16:54:03.131406
Stop comment: normal stop

Fixed

com2019

Sat Apr 6 10:09:53 2019

Adrien Matta wrote:

We see more (d.3he) than dt

there seems to be another reaction going up to 40deg Lab

The line at 40 deg matched 16O(d,4He)14N g.s.

Fixed

com2019

Sat Apr 6 12:38:17 2019

Tac Si Vamos with gate on GATCONFMASTER > 8

Attachment 1: T_Si_Vamos_GATCONF_inf_8.png

Fixed

com2019

Sat Apr 6 17:22:14 2019

*************************
Run : 258
*************************

Start at 2019-04-06 17:21:23.706857

Run Conditions:

Operator: A. Lemasson
Trigger: MUST2 + MG BWd + VAMOS % 100
Target/Source: CD2 1mg
Beam: 16O
Intensity: ~ 40 kHz

Comment: square back and threshold MG && MM 800keV
Second Run

Stop at 2019-04-06 17:50:11.159270
Stop comment: stop for water leakeage

Fixed

com2019

Sat Apr 6 18:32:26 2019

*************************
Run : 259
*************************

Start at 2019-04-06 18:31:46.303930

Run Conditions:

Operator:
Trigger: MUST2 + MG BWd + VAMOS % 100
Target/Source: CD2 1mg
Beam: no
Intensity:

Comment: Target holder in position

Fixed

com2019

Sat Apr 6 18:37:27 2019

Freddy, Val�rian, Adrien, Antoine

OFFLINE

General

N/A

Excitation Energy and Gammas

This is it!

M.J.
------------------------------------
Measured Statistics (Run 256) about 10h:

GS + 1st excited:

MUGAST(Ex rec) : 2519 (bckg subtracted)
MUGAST(Ex rec + VAMOS) : 1282

Gamma(871keV) : 125 (if subtract bg only 80 cnts)
Gamma(871keV)+MUGAST(Ex) : 68
Gamma(871keV)+MUGAST(Ex)+ VAMOS IC: 27

Deduced gamma efficiency: 68/(2519/2) = 5.4 % ---> to be investigated

------------------------------------

Fit Results Excitation Energy (gated by VAMOS) - sum of two gaussian fit

  EXT PARAMETER                                   STEP         FIRST                                     
  NO.   NAME      VALUE            ERROR          SIZE      DERIVATIVE                                   
   1  p0           1.91342e+02   1.01289e+01   1.78576e-02   3.82264e-06                                 
   2  p1           2.71036e-02   1.14447e-02   2.14103e-05   4.44817e-03                                 
   3  p2           2.13115e-01   1.00329e-02   1.45208e-05   1.03530e-02                                 
   4  p3           1.63729e+02   9.57766e+00   1.64222e-02   1.51782e-05                                 
   5  p4           8.83234e-01   1.59487e-02   3.11205e-05  -8.55421e-03                                 
   6  p5           2.37851e-01   1.88040e-02   2.48042e-05   1.09859e-02

FWHM(gs): 500 keV
FWHM(871keV): 550 keV
Simulation: 500 keV

Attachment 1: run256_allnight_bilan.png

Fixed

com2019

Sun Apr 7 10:26:34 2019

save File KineALL.cxx

1 #include"NPReaction.h"
2
3 TCutG* ETOF=NULL;
4 TCutG* EDE=NULL;
5 TChain* chain=NULL ;
6 TCanvas* c1 = NULL;
7 TCanvas* c2 = NULL;
8
9 ////////////////////////////////////////////////////////////////////////////////
10 void LoadCuts(){
11 /* TFile* File_ETOF = new TFile("cuts/ETOF.root","READ");
12 ETOF = (TCutG*) File_ETOF->FindObjectAny("ETOF");
13
14 TFile* File_EDE = new TFile("cuts/EDE.root","READ");
15 EDE= (TCutG*) File_EDE->FindObjectAny("EDE");
16 */
17 }
18
19 ////////////////////////////////////////////////////////////////////////////////
20 void LoadChainNP(){
21 chain = new TChain("PhysicsTree");
22 chain->Add("../../Outputs/Analysis/Phy253.root");
23 chain->Add("NPOutput/runN_0256.root");
24 chain->Add("NPOutput/runN_0256_1.root");
25 chain->Add("NPOutput/runN_0256_2.root");
26 chain->Add("NPOutput/runN_0256_3.root");
27 chain->Add("NPOutput/runN_0256_4.root");
28 chain->Add("NPOutput/runN_0256_5.root");
29 chain->Add("NPOutput/runN_0256_6.root");
30 chain->Add("NPOutput/runN_0256_7.root");
31 chain->Add("NPOutput/runN_0256_8.root");
32 chain->Add("NPOutput/runN_0256_9.root");
33 chain->Add("NPOutput/runN_0256_10.root");
34 chain->Add("NPOutput/runN_0256_11.root");
35 chain->Add("NPOutput/runN_0256_12.root");
36 chain->Add("NPOutput/runN_0256_13.root");
37 //chain->Add("NPOutput/runN_0256_14.root");
38 //chain->Add("NPOutput/runN_0256_15.root");
39 // chain->Add("NPOutput/runN_0256_16.root");
40 }
41 // void LoadChainNPVAMOS(){
42 // chain = new TChain("PhysicsTree");
43 // chain->Add("NPOutput/runVAMOS_0256_0.root");
44 // }
45 void LoadChainNPVAMOS(){
46 chain = new TChain("PhysicsTree");
47 chain->Add("NPOutput/Phy256_ALL.root");
48 //chain->Add("NPOutput/testFF.root");
49 }
50 //void LoadChainVAMOS(){
51 // chain = new TChain("AD");
52 // chain->Add("../../Outputs/Analysis/Phy253.root");
53 // chain->Add("NPOutput/runN_0256.root");
54 // chain->Add("NPOutput/runN_0256_1.root");
55 // chain->Add("NPOutput/runN_0256_2.root");
56 // chain->Add("NPOutput/runN_0256_3.root");
57 // chain->Add("NPOutput/runN_0256_4.root");
58 // chain->Add("NPOutput/runN_0256_5.root");
59 // chain->Add("NPOutput/runN_0256_6.root");
60 // chain->Add("NPOutput/runN_0256_7.root");
61 // chain->Add("NPOutput/runN_0256_8.root");
62 // chain->Add("NPOutput/runN_0256_9.root");
63 // chain->Add("NPOutput/runN_0256_10.root");
64 // chain->Add("NPOutput/runN_0256_11.root");
65 // chain->Add("NPOutput/runN_0256_12.root");
66 // chain->Add("NPOutput/runN_0256_13.root");
67 // chain->Add("NPOutput/runN_0256_14.root");
68 // chain->Add("NPOutput/runN_0256_15.root");
69 // chain->Add("NPOutput/runN_0256_16.root");
70 // }
71 //////////////////////////////////////////////////////////////////////////////////
72 void LoadEventList(){
73 }
74 ////////////////////////////////////////////////////////////////////////////////
75 void plot_kine(NPL::Reaction r, double Ex,Color_t c,int w, int s){
76 r.SetExcitation4(Ex);
77
78 TGraph* g= r.GetKinematicLine3();
79 g->SetLineColor(c) ;
80 g->SetLineStyle(s) ;
81 g->SetLineWidth(w) ;
82
83 g->Draw("c");
84 }
85 ////////////////////////////////////////////////////////////////////////////////
86 void plot_state(double Ex,double max,Color_t c,int w, int s){
87 TLine* line = new TLine(Ex,0,Ex,max) ;
88 line->SetLineColor(c) ;
89 line->SetLineStyle(s) ;
90 line->SetLineWidth(w) ;
91 line->Draw();
92
93 }
94 ////////////////////////////////////////////////////////////////////////////////
95 void KineALL(){
96 gStyle->SetOptStat(111111111);
97 // Load stuff
98 // LoadChainNP();
99 LoadChainNPVAMOS();
100
101 // for Kinematic calculation
102 NPL::Reaction O17("16O(d,p)17O@96");
103 NPL::Reaction O16dd("16O(d,d)16O@96");
104 NPL::Reaction O16pp("16O(p,p)16O@96");
105 NPL::Reaction O16C12("16O(12C,12C)16O@96");
106 NPL::Reaction O15("16O(d,t)15O@96");
107 NPL::Reaction O16dh("16O(d,3He)15N@96");
108 NPL::Reaction O16d4h("16O(d,4He)15N@96");
109 Color_t color_dp=kOrange+7;
110 Color_t color_dd=kRed;
111 Color_t color_pp=kGreen;
112 Color_t color_dt=kAzure+7;
113 Color_t color_ca=kBlack;
114 Color_t color_dh=kMagenta;
115 Color_t color_d4h=kCyan;
116
117 c1 = new TCanvas("Example1","Example1",0,0,600,600);
118 c2 = new TCanvas("Example2","Example2",0,0,600,600);
119
120 c1->Divide(1,3);
121 c2->Divide(1,4);
122 // Kinematic Plot
123 c1->cd(1);
124 chain->Draw("ELab:ThetaLab>>hKine(720,0,180,1000,0,40)","ELab>0&&ThetaLab>0");
125 TH2F* hKine = (TH2F*) gDirectory->Get("hKine");
126 hKine->Draw("col");
127 hKine->GetXaxis()->SetTitle("#Theta_{Lab} (deg)");
128 hKine->GetYaxis()->SetTitle("E_{Lab} (MeV)");
129
130 // 15N
131 plot_kine(O16dh,0,color_dh,2,1);
132
133 // 17N
134 plot_kine(O16d4h,0,color_d4h,2,1);
135 plot_kine(O16d4h,2.312,color_d4h,2,1);
136
137 // O17
138 plot_kine(O17,0,color_dp,2,1);
139 plot_kine(O17,0.87073,color_dp,2,1);
140 plot_kine(O17,3.05536,color_dp,1,1);
141 // Sn
142 plot_kine(O17,4.143,color_dp,1,2);
143
144 // O15
145 plot_kine(O15,0,color_dt,2,1);
146 //plot_kine(O15,5.183,color_dt,2,1);
147 //plot_kine(O15,5.241,color_dt,1,1);
148 //plot_kine(O15,6.176,color_dt,1,1);
149 // Sp
150 //plot_kine(O15,7.297,color_dt,1,2);
151 // Elastic
152 plot_kine(O16dd,0,color_dd,2,1);
153 plot_kine(O16pp,0,color_pp,1,1);
154 plot_kine(O16C12,0,color_ca,1,2);
155
156 // Ex Plot
157 c1->cd(2);
158 int bin = 1000;
159 double start=-10;
160 double end = 10;
161 double binning = (end-start)/bin;
162 chain->Draw(Form("Ex>>hEx(%d,%f,%f)",1000,start,end),"ThetaLab>90");
163 TH2F* hEx = (TH2F*) gDirectory->Get("hEx");
164
165 chain->Draw(Form("Ex>>hEx2(%d,%f,%f)",1000,start,end),"ThetaLab>90&&IC[1]>2");
166 TH2F* hEx2 = (TH2F*) gDirectory->Get("hEx2");
167
168 hEx->GetXaxis()->SetTitle("E_{17O} (MeV)");
169 hEx->GetYaxis()->SetTitle(Form("counts / %.0f keV",binning*1000));
170
171 hEx->Draw();
172 hEx->Rebin(8);
173 c1->cd(3);
174 hEx2->Draw();
175 hEx2->Rebin(8);
176
177 c2->cd(1);
178 chain->Draw("EDC>>hEDC(2000,0,2)","");
179 TH2F* hEDC = (TH2F*) gDirectory->Get("hEDC");
180 hEDC->Rebin(4);
181
182 c2->cd(2);
183 chain->Draw("EDC>>hEDC1(2000,0,2)","Ex>0.2 && Ex<1.4");
184 TH2F* hEDC1 = (TH2F*) gDirectory->Get("hEDC1");
185 hEDC1->Rebin(4);
186
187 c2->cd(3);
188 chain->Draw("EDC>>hEDC2(2000,0,2)","IC[1]>2");
189 TH2F* hEDC2 = (TH2F*) gDirectory->Get("hEDC2");
190 hEDC2->Rebin(4);
191
192 c2->cd(4);
193 // chain->Draw("EDC>>hEDC3(2000,0,2)","MW_Nr>0");
194 chain->Draw("EDC>>hEDC3(2000,0,2)","IC[1]>2 && Ex>0.2 && Ex<1.4");
195 TH2F* hEDC3 = (TH2F*) gDirectory->Get("hEDC3");
196 hEDC3->Rebin(4);
197
198 double max = hEx->GetMaximum();
199 plot_state(0,max,color_dp,2,1);
200 plot_state(0.87073,max,color_dp,2,1);
201 plot_state(3.05536,max,color_dp,1,1);
202 plot_state(4.143,max,color_dp,1,2);
203 }

Fixed

com2019

Sun Apr 7 16:09:22 2019

Freddy, Val�rian, Adrien, Antoine

OFFLINE

General

N/A

Excitation Energy and Gammas

Freddy, Val�rian, Adrien, Antoine wrote:

  EXT PARAMETER                                   STEP         FIRST                                     
  NO.   NAME      VALUE            ERROR          SIZE      DERIVATIVE                                   
   1  p0           1.91342e+02   1.01289e+01   1.78576e-02   3.82264e-06                                 
   2  p1           2.71036e-02   1.14447e-02   2.14103e-05   4.44817e-03                                 
   3  p2           2.13115e-01   1.00329e-02   1.45208e-05   1.03530e-02                                 
   4  p3           1.63729e+02   9.57766e+00   1.64222e-02   1.51782e-05                                 
   5  p4           8.83234e-01   1.59487e-02   3.11205e-05  -8.55421e-03                                 
   6  p5           2.37851e-01   1.88040e-02   2.48042e-05   1.09859e-02

FWHM(gs): 500 keV
FWHM(871keV): 550 keV
Simulation: 500 keV

Some remarks/ correction :

The top right plot named all gammas in AGATA is in fact already assuming that the multiplicity of charged particle in MUGAST is 1
All the gammas spectra plotted correspond only to nbtrack=1 (fold = 1)

Fixed

com2019

Sun Apr 7 21:35:51 2019

Valerian, Freddy and Antoine

RUN

General

N/A

Overnight run259 crashed during the night

Run 259 which was an overnight run crashed (During the night between the 6 and the 7th of April).

The result of this crash was that the the acquisition started writing a 4gb file each seconds for 80 seconds. We then had to remove the run (~300gb) as the quota allowed by GANIL to the mugast account was reached.
We also transfered the output of the Nptool analysis to the /data/evamosX/com2019/Analysis/nptool_outputs/ and made a symbolic link to this folder in the nptool/Output in order not to ever reach the quota again (Antoine seems to have no quota, or at least multiple Pb of quota...)

Fixed

com2019

Sun Apr 7 21:43:29 2019

*************************
Run : 261
*************************

Start at 2019-04-07 21:43:00.572161

Run Conditions:

Operator: Iulian Marlene
Trigger: MUST2 + MG BWd + VAMOS % 100
Target/Source: CD2 1mg
Beam: no
Intensity:

Comment: Target holder in position

junk

Stop at 2019-04-07 21:44:59.581974
Stop comment: junk

Fixed

com2019

Sun Apr 7 21:45:35 2019

*************************
Run : 262
*************************

Start at 2019-04-07 21:45:06.409673

Run Conditions:

Operator: Iulian Marlene
Trigger: MUST2 + MG BWd + VAMOS % 100
Target/Source: CD2 1mg
Beam: no
Intensity:

Comment: Target holder removed
background check

Stop at 2019-04-08 07:38:35.058779
Stop comment:
end run
Alfa background is continuous in energy particularly in the backward detectors

Attachment 1: bilan_run262_bck_mugast4.png

Attachment 2: bilan_run262_bck_mugast3.png

Attachment 3: bilan_run262_bck_mugast2.png

Attachment 4: bilan_run262_bck_mugast1.png

Goto page Previous 1, 2, 3 ... 67, 68, 69 ... 130, 131, 132 Next

ELOG Home Page