Skip to content
Snippets Groups Projects
Commit 447212f9 authored by Greg Christian's avatar Greg Christian
Browse files

Fixed error message w/ trying to calculate energy loss for negative values

parent e8ccf3e0
No related branches found
No related tags found
No related merge requests found
......@@ -324,10 +324,12 @@ void Analysis::TreatEvent(){
Si_E_TH = TH->Strip_E[countTiaraHyball];
Energy = Si_E_TH; // calibration for hyball is in MeV
// Correct for energy loss using the thickness of the target and the dead layer
ELab = LightSi.EvaluateInitialEnergy( Energy ,0.61*micrometer , ThetaTHSurface); // equivalent to 0.1 um of Aluminum
if(Energy > 0){
ELab = LightSi.EvaluateInitialEnergy( Energy ,0.61*micrometer , ThetaTHSurface); // equivalent to 0.1 um of Aluminum
//by Shuya 170530
//if(ThetaNormalTarget < halfpi) ELab = LightCBacking.EvaluateInitialEnergy( ELab ,0.044*micrometer , ThetaNormalTarget); //10 ug/cm2 carbon
ELab = LightTarget.EvaluateInitialEnergy( ELab ,TargetThickness/2., ThetaNormalTarget);
//if(ThetaNormalTarget < halfpi) ELab = LightCBacking.EvaluateInitialEnergy( ELab ,0.044*micrometer , ThetaNormalTarget); //10 ug/cm2 carbon
ELab = LightTarget.EvaluateInitialEnergy( ELab ,TargetThickness/2., ThetaNormalTarget);
}
/////////////////////////////
// Part 3 : Excitation Energy Calculation
......
......@@ -6,9 +6,9 @@ ConfigTiaraHyball
STRIP_ENERGY_MATCHING_NUMBER_OF_SIGMA 5
%If you want to use Ring signal instead of Ring-Sector average, comment this out
%To get weighted averaged energy, you need to input these below (default; RingE 18 keV, SectorE=22 keV)
TAKE_E_RING_SECTOR_Average
SIGMA_RING_E 0.018
SIGMA_SECTOR_E 0.022
% TAKE_E_RING_SECTOR_Average
% SIGMA_RING_E 0.018
% SIGMA_SECTOR_E 0.022
TAKE_E_RING
%TAKE_E_SECTOR
%TAKE_T_RING
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment