Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Updated CS2AreaFractionTopExtended for SR1 LowER #178

Open
wants to merge 2 commits into
base: master
Choose a base branch
from
Open

Updated CS2AreaFractionTopExtended for SR1 LowER #178

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
2cf641d
Updated CS2AreaFractionTopExtended for SR1 LowER
DCichon Nov 25, 2019
493c683
Typo fix
DCichon Nov 25, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
78 changes: 78 additions & 0 deletions lax/lichens/postsr1.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -193,6 +193,84 @@ def pre(self, df):
return df


class CS2AreaFractionTopExtendedLowEr(StringLichen):
""""An extension of CS2AreaFractionTop to the entire S2 range
with a designed acceptance of 98% (low ER analysis version).
It is defined in the (cxys2, cs2_aft) space, with:
cxys2 = (cs2_top + cs2_bottom) / s2_lifetime_correction
cs2_aft = cs2_top / (cs2_top + cs2_bottom)
Events where cxys2 > 1713000 PE or cxys2 < 60 PE are not cut.
This cut should be used with a pax version of at least 6.10.0
due to a major update of the S2 desaturation correction.
This version has been designed for SR1 analysis
SR1 S2Width and post-SR1 S2PatternLikelihood have been included in
the preselection for defining the cut, which needs to be taken into
account for data outside their respective regions of validity.
Required minitrees: Corrections
Defined with pax version: 6.10.1
Wiki notes: xenon:xenon1t:double_beta:cs2_aft_extension
xenon:xenon1t:double_beta:cs2_aft_extension_rejection_estimate
xenon:xenon1t:analysis:sciencerun2:cs2_aft
xenon:xenon1t:analysis:sciencerun2:cs2_aft_presel_update
xenon:xenon1t:analysis:sciencerun1:cs2_aft_presel_update_for_lower
Contact: Dominick Cichon <[email protected]>"""

version = 1

top_bound_string = ('(6.233004E-01 + 1.961573E-07 * cxys2 +'
' -5.556886E-13 * cxys2**2 + 7.106575E-19 * cxys2**3 +'
' -4.078553E-25 * cxys2**4 + 8.629125E-32 * cxys2**5 +'
' 2.238235E+00 / sqrt(cxys2) + -7.652169E+00 / cxys2)')
bot_bound_string = ('(6.240076E-01 + -1.454419E-08 * cxys2 +'
' 4.270098E-14 * cxys2**2 + -8.150741E-20 * cxys2**3 +'
' 5.550419E-26 * cxys2**4 + -1.231424E-32 * cxys2**5 +'
' -1.767875E+00 / sqrt(cxys2) + 1.246680E+00 / cxys2)')

string = ('((' + top_bound_string + ' > cs2_aft) & (' + bot_bound_string +
' < cs2_aft)) | (cxys2 > 1713000) | (cxys2 < 60)')

def pre(self, df):
df.loc[:, 'cxys2'] = ((df['cs2_top'] + df['cs2_bottom']) /
df['s2_lifetime_correction'])
df.loc[:, 'cs2_aft'] = df['cs2_top'] / (df['cs2_top'] +
df['cs2_bottom'])
return df


class CS2AreaFractionTopExtendedLowEr99Percent(StringLichen):
""""A version of CS2AreaFractionTopExtendedLowER with a target acceptance of
99%. See the docstring of CS2AreaFractionTopExtendedLowER for more details.
Required minitrees: Corrections
Defined with pax version: 6.10.1
Wiki notes: xenon:xenon1t:double_beta:cs2_aft_extension
xenon:xenon1t:double_beta:cs2_aft_extension_rejection_estimate
xenon:xenon1t:analysis:sciencerun2:cs2_aft
xenon:xenon1t:analysis:sciencerun2:cs2_aft_presel_update
xenon:xenon1t:analysis:sciencerun1:cs2_aft_presel_update_for_lower
Contact: Dominick Cichon <[email protected]>"""

version = 1

top_bound_string = ('(6.226176E-01 + 2.333159E-07 * cxys2 +'
' -6.895548E-13 * cxys2**2 + 9.006261E-19 * cxys2**3 +'
' -5.242889E-25 * cxys2**4 + 1.121020E-31 * cxys2**5 +'
' 2.490721E+00 / sqrt(cxys2) + -8.389662E+00 / cxys2)')
bot_bound_string = ('(6.241466E-01 + -2.600373E-08 * cxys2 +'
' 6.945772E-14 * cxys2**2 + -1.312317E-19 * cxys2**3 +'
' 9.570400E-26 * cxys2**4 + -2.312927E-32 * cxys2**5 +'
' -2.012173E+00 / sqrt(cxys2) + 1.082919E+00 / cxys2)')

string = ('((' + top_bound_string + ' > cs2_aft) & (' + bot_bound_string +
' < cs2_aft)) | (cxys2 > 1713000) | (cxys2 < 60)')

def pre(self, df):
df.loc[:, 'cxys2'] = ((df['cs2_top'] + df['cs2_bottom']) /
df['s2_lifetime_correction'])
df.loc[:, 'cs2_aft'] = df['cs2_top'] / (df['cs2_top'] +
df['cs2_bottom'])
return df


class MisIdS1SingleScatter(Lichen):
"""Cut to target the shoulder on Kr83m data due to mis-identified krypton events.

Expand Down