-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathBabyIAXO_2024.rml
110 lines (88 loc) · 4.55 KB
/
BabyIAXO_2024.rml
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<BabyIAXO>
<globals>
<constant name="Pi" value="3.1415927" />
<constant name="Bores" value="2" />
<constant name="SpotRadius" value="0.8" />
<constant name="BField" value="1.69" />
<constant name="Years" value="1.5" />
<constant name="MagnetRadius" value="35" />
<constant name="MagnetLength" value="11" />
<constant name="GasLength" value="5" />
<constant name="WindowEff" value="0.81" />
<constant name="OpticsEff" value="0.47*0.9" />
<constant name="DetectorEff" value="0.66" />
<constant name="BckRate" value="1e-7" />
</globals>
<TRestSensitivity name="VacuumPhase" strategy="nodes" >
<TRestExperiment name="Vacuum" exposureTime="Years*365*12hr" useAverage="true">
<!-- Background -->
<TRestComponentFormula name="Flat7" nature="background" >
<parameter name="formulaUnits" value="keV^-1" />
<cVariable name="energy" range="(0,10)keV" bins="20" />
<formula name="bck" expression="Bores*BckRate*Pi*SpotRadius*SpotRadius" />
</TRestComponentFormula>
<!-- Signal -->
<TRestAxionHelioscopeSignal name="Vacuum" nature="signal" parameter="ma"
conversionType="IAXO" bores="Bores"
magnetRadius="MagnetRadius cm" magnetLength="MagnetLength m" magnetStrength="BField T"
opticsEfficiency="OpticsEff" windowEfficiency="WindowEff" >
<cVariable name="energy" range="(0,10)keV" bins="20" />
<parameter name="firstParameterValue" value="0.001" />
<parameter name="lastParameterValue" value="10" />
<parameter name="stepParameterValue" value="1.02" />
<parameter name="exponential" value="true" />
<TRestAxionSolarQCDFlux name="LennertHoofPrimakoff" couplingType="g_ag" couplingStrength="1.e-10"
fluxDataFile="Primakoff_LennertHoof_202203.dat" seed="137" />
</TRestAxionHelioscopeSignal>
</TRestExperiment>
</TRestSensitivity>
<TRestSensitivity name="CombinedPhase" strategy="nodes" >
<TRestExperiment name="Vacuum" exposureTime="1.5*300*12hr" useAverage="true">
<!-- Background -->
<TRestComponentFormula name="Flat7" nature="background" >
<parameter name="formulaUnits" value="keV^-1" />
<cVariable name="energy" range="(0,10)keV" bins="20" />
<formula name="bck" expression="Bores*BckRate*Pi*SpotRadius*SpotRadius" />
</TRestComponentFormula>
<!-- Signal -->
<TRestAxionHelioscopeSignal name="Vacuum" nature="signal" parameter="ma"
conversionType="IAXO" bores="Bores"
magnetRadius="MagnetRadius cm" magnetLength="MagnetLength m" magnetStrength="BField T"
opticsEfficiency="OpticsEff" windowEfficiency="WindowEff" detectorEfficiency="DetectorEff"
gasLength="GasLength m" >
<cVariable name="energy" range="(0,10)keV" bins="20" />
<TRestAxionSolarQCDFlux name="LennertHoofPrimakoff" couplingType="g_ag" couplingStrength="1.e-10"
fluxDataFile="Primakoff_LennertHoof_202203.dat" seed="137" />
<parameter name="firstParameterValue" value="0.001" />
<parameter name="lastParameterValue" value="10" />
<parameter name="stepParameterValue" value="1.02" />
<parameter name="exponential" value="true" />
</TRestAxionHelioscopeSignal>
</TRestExperiment>
<TRestExperimentList name="GasPhase" exposureTime="0"
componentPattern="output/SignalsBabyIAXO_2024.root" experimentsFile="output/KSVZ.settings" useAverage="true">
<TRestComponentFormula name="Flat7" nature="background" >
<parameter name="formulaUnits" value="keV^-1" />
<cVariable name="energy" range="(0,10)keV" bins="20" />
<formula name="bck" expression="Bores*BckRate*Pi*SpotRadius*SpotRadius" />
</TRestComponentFormula>
</TRestExperimentList>
</TRestSensitivity>
<TRestAxionHelioscopeSignal name="GasSignal" nature="signal" parameter="ma"
conversionType="IAXO" bores="Bores"
magnetRadius="MagnetRadius cm" magnetLength="MagnetLength m" magnetStrength="BField T"
opticsEfficiency="OpticsEff" windowEfficiency="WindowEff" detectorEfficiency="DetectorEff"
gasLength="GasLength m" >
<cVariable name="energy" range="(0,10)keV" bins="20" />
<parameter name="firstParameterValue" value="0.001" />
<parameter name="lastParameterValue" value="10" />
<parameter name="stepParameterValue" value="1.02" />
<parameter name="exponential" value="true" />
<TRestAxionSolarQCDFlux name="LennertHoofPrimakoff" couplingType="g_ag" couplingStrength="1.e-10"
fluxDataFile="Primakoff_LennertHoof_202203.dat" seed="137" />
<TRestAxionBufferGas name="helium" >
<gas name="He" density="0.0025e-6g/cm^3"/>
</TRestAxionBufferGas>
</TRestAxionHelioscopeSignal>
</BabyIAXO>