|
Fixed, e744s, Fri Apr 5 22:39:39 2019, OttoRun, RUN, General, Ti-t, Run255
|
*************************
Run : 255
*************************
Start at 2019-04-05 22:39:35.688874
Run Conditions:
- Operator:
- Trigger:
- Target/Source:
- Beam:
- Intensity:
Comment:
Stop at 2019-04-05 22:41:54.607652
Stop comment: |
|
Fixed, com2019, Fri Apr 5 22:58:04 2019, OttoRun, RUN, General, CD2, Run256
|
*************************
Run : 256
*************************
Start at 2019-04-05 22:56:31.606922
Run Conditions:
- Operator: A. Lemasson
- Trigger: MUST2 + MG BWd + VAMOS % 100
- Target/Source: CD2 1mg
- Beam: 16O
- Intensity: ~ 40 kHz
Comment:
First Run With Coupling with AGATA
Stop at 2019-04-06 09:44:30.727058
Stop comment: |
Fixed, com2019, Fri Apr 5 23:19:04 2019, Adrien et Valerian, OFFLINE, General, CD2, Experimental Kinematic Line from dp 
|
We used run 253 to obtain a convincing Kine line and Ex spectrum at backward angle.
Stat is still limited but the distribution is convincing. |
|
Fixed, com2019, Fri Apr 5 23:32:29 2019, Fairouz, CATS, General, CD2, Electronic changes
|
CATS detectors were removed.
For the MUVI stop, we send now HF propre.
The CATS1 signals sent to different TACS are now the HF signal |
|
Pinned, com2019, Fri Apr 5 23:34:46 2019, Adrien et Valerian, OFFLINE, General, N/A, Converting data and adding stat to Kine Plot
|
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, clement, OFFLINE, General, CD2, Experimental Kinematic Line from dp + gamma ?????
|
| Triangle is 870 keV |
|
Fixed, com2019, Sat Apr 6 01:06:18 2019, mugast, OFFLINE, General, CD2, Experimental Kinematic Line from dp
|
| Kinematic Line obtained from run 256+253 |
|
Fixed, com2019, Sat Apr 6 08:23:09 2019, mugast, OFFLINE, General, CD2, 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). |
|
Fixed, com2019, Sat Apr 6 08:38:51 2019, Adrien Matta, OFFLINE, General, CD2, 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 |
|
Fixed, com2019, Sat Apr 6 09:52:47 2019, OttoRun, RUN, General, CD2, Run257
|
*************************
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, OFFLINE, General, CD2, We see d,4He
|
| 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, Marlene et Valerian, RUN, General, CD2, Tac Si Vamos with gate on GATCONFMASTER > 8
|
|
|
Fixed, com2019, Sat Apr 6 17:22:14 2019, OttoRun, RUN, General, CD2, Run258
|
*************************
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, OttoRun, RUN, General, CD2, Run259
|
*************************
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 |
|
Fixed, com2019, Sun Apr 7 10:26:34 2019, mugast, OFFLINE, General, N/A, 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: | 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 |
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, OttoRun, RUN, General, CD2, Run261
|
*************************
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, OttoRun, RUN, General, CD2, Run262
|
*************************
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
|
|