- genevb's home page
- Posts
- 2024
- 2023
- 2022
- September (1)
- 2021
- 2020
- 2019
- December (1)
- October (4)
- September (2)
- August (6)
- July (1)
- June (2)
- May (4)
- April (2)
- March (3)
- February (3)
- 2018
- 2017
- December (1)
- October (3)
- September (1)
- August (1)
- July (2)
- June (2)
- April (2)
- March (2)
- February (1)
- 2016
- November (2)
- September (1)
- August (2)
- July (1)
- June (2)
- May (2)
- April (1)
- March (5)
- February (2)
- January (1)
- 2015
- December (1)
- October (1)
- September (2)
- June (1)
- May (2)
- April (2)
- March (3)
- February (1)
- January (3)
- 2014
- December (2)
- October (2)
- September (2)
- August (3)
- July (2)
- June (2)
- May (2)
- April (9)
- March (2)
- February (2)
- January (1)
- 2013
- December (5)
- October (3)
- September (3)
- August (1)
- July (1)
- May (4)
- April (4)
- March (7)
- February (1)
- January (2)
- 2012
- December (2)
- November (6)
- October (2)
- September (3)
- August (7)
- July (2)
- June (1)
- May (3)
- April (1)
- March (2)
- February (1)
- 2011
- November (1)
- October (1)
- September (4)
- August (2)
- July (4)
- June (3)
- May (4)
- April (9)
- March (5)
- February (6)
- January (3)
- 2010
- December (3)
- November (6)
- October (3)
- September (1)
- August (5)
- July (1)
- June (4)
- May (1)
- April (2)
- March (2)
- February (4)
- January (2)
- 2009
- November (1)
- October (2)
- September (6)
- August (4)
- July (4)
- June (3)
- May (5)
- April (5)
- March (3)
- February (1)
- 2008
- 2005
- October (1)
- My blog
- Post new blog entry
- All blogs
Cutting events with corrupt RICH scalers in analyses
Following Chris Powell's work, which he You do not have access to view this node at an S&C meeting, to resolve the problem with RICH scalers I have Finding a cut for corrupt scalers in Run 8 dAu, I have implemented some code to remove events with corrupt RICH scalers from the P08ic and P08ie productions of Run 8 dAu. The code can be used two ways:
I hope this is useful. Feedback is welcome. Note that this is relevant only for events with TPC/TPX, but I did not write any code to determine whether these detectors are in the data; that is left to the user if the user cares.
gROOT->Macro("loadMuDst.C");
gSystem->Load("StMuAnalysisMaker");
...
for_loop_over_MuDst_events {
...
if (!StMuScalerCutter::accept(mudst->event())) continue;
...
}
gROOT->Macro("loadMuDst.C");
gSystem->Load("StMuAnalysisMaker");
StChain chain;
StMuDstMaker maker1(...);
StMuScalerCutter maker2();
maker2.SetAttr(".Privilege",1);
...
for_loop_over_MuDst_events {
...
int result = chain.Make();
if (result == kStSkip) continue;
...
}
...and for those who don't use a chain:
gROOT->Macro("loadMuDst.C");
gSystem->Load("StMuAnalysisMaker");
StMuDstMaker maker1(...);
StMuScalerCutter maker2();
...
for_loop_over_MuDst_events {
...
int result1 = maker1.Make();
int result2 = maker2.Make();
if (result2 == kStSkip) continue;
...
}
The code is located in StRoot/StMuAnalysisMaker and should be available in DEV starting February 5, 2009.
Many thanks go to Chris for his work on this topic.
-Gene
- genevb's blog
- Login or register to post comments