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.

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

Fix bug in FV z-optimization analysis #145

Merged
merged 4 commits into from
Feb 26, 2018
Merged

Conversation

pdeperio
Copy link
Contributor

@pdeperio pdeperio commented Feb 25, 2018

Re-optimize z-cut based on bug fixed FV analysis (details below) and tighten z cut to -8 cm based on appearance of leakage gas event in SR1 BG data.

After review by @lgrandi, I think I may have found a bug where the ER KDE matrix was transposed when it shouldn't be. Now if we look at "er_kde" in the following figure (black line old cut), it seems to make more sense with the population extending across the top instead of down the side (which is not yet the TPC edge) and the bulge transposed now to the bottom-right corresponding to the other corner with large expected materials contribution:
https://xe1t-wiki.lngs.infn.it/lib/exe/fetch.php?cache=&media=36637720-844071d4-19b0-11e8-99b1-d4bdf0df9e88.png

Re-running the notebook and fit results in the following line:

SR0 | SR1

The new cut is significantly shallower and follows the background shape much better.

@feigaodm @mcfatelin If this is approved, please add to stack prior to next template regeneration.

and tighten z cut based on appearance of leakage gas event
feigaodm
feigaodm previously approved these changes Feb 25, 2018
@pdeperio
Copy link
Contributor Author

pdeperio commented Feb 25, 2018

Fixed another bug in the z_min finder algorithm, which then results in a cut on the bottom-right corner when using WIMP signal region (as expected from ER materials contribution dominating in the signal region instead of surface background), as shown in the following figure:
https://xe1t-wiki.lngs.infn.it/lib/exe/fetch.php?cache=&media=36638549-0c1fc978-19c6-11e8-9c7f-d22a6478e95d.png

Fitting the bottom corner now yields:

10 GeV WIMP | 50 GeV WIMP

Pretty much the same so use 10 GeV which is a bit tighter.

@pdeperio pdeperio dismissed feigaodm’s stale review February 25, 2018 06:01

Fixed another bug right after you approved :(

@pdeperio pdeperio requested a review from feigaodm February 25, 2018 06:01
mcfatelin
mcfatelin previously approved these changes Feb 25, 2018
@mcfatelin mcfatelin dismissed their stale review February 25, 2018 15:28

Have one comment

string = "(-95 < z_3d_nn) & (z_3d_nn < -7) & \
(z_3d_nn < 27.2929 - 0.0276385*r_3d_nn*r_3d_nn)"
version = 1
string = "(-95 < z_3d_nn) & (z_3d_nn < -8) & \
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Where's the R cut?

Copy link
Contributor

@coderdj coderdj Feb 25, 2018

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It's logical 'and' between all 4 lines.

Edit: I think you're right. Max r is missing... also in the last version right?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Added in 9e9eeeb: r_3d_nn < 42.8387

@@ -260,12 +260,12 @@ class FiducialZOptimized(StringLichen):
This first version is a bit rough and may be improved by understanding the underlying BG KDE
and errors better.

https://xe1t-wiki.lngs.infn.it/doku.php?id=xenon:xenon1t:analysis:sciencerun1:summary_fiducial_volume
https://xe1t-wiki.lngs.infn.it/doku.php?id=xenon:xenon1t:analysis:sciencerun1:summary_fiducial_volume_v2
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@jingqiangye
Copy link
Contributor

Hey the fiducial masses for different r bins under current FV version:

R binning(cm) 34.5904 35.0729 35.5489 36.0185 36.4822 36.9400 37.3922 37.8390 38.2806 38.7172 39.1489 39.5758 39.9983 40.4163 40.8300 41.2396 41.6451 42.0468 42.4446 42.8387
Fiducial Mass(kg) 918.3 942.8 967.3 991.5 1015.6 1039.6 1063.3 1086.9 1110.4 1133.6 1156.9 1179.7 1202.3 1224.1 1245.5 1266.1 1286.1 1305.7 1324.6 1343.4

@sreichard sreichard merged commit 20cc025 into master Feb 26, 2018
@sreichard sreichard deleted the FiducialZOptimized_bugfix branch February 26, 2018 09:44
@sreichard sreichard restored the FiducialZOptimized_bugfix branch February 26, 2018 09:45
@feigaodm feigaodm deleted the FiducialZOptimized_bugfix branch February 26, 2018 17:59
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

6 participants