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