- genevb's home page
- Posts
- 2024
- 2023
- 2022
- September (1)
- 2021
- 2020
- 2019
- 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
- 2013
- 2012
- 2011
- January (3)
- 2010
- February (4)
- 2009
- 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