StMuRpsUtil - Roman Pot data analysis utilities (afterburner)

Updated on Wed, 2017-01-25 14:43. Originally created by rafal_s on 2016-03-20 04:31.

computing

StMuRpsUtil (under development, for testing purposes only!!!)

Should you have any questions/comments/remarks regarding this module please contact rafal.sikora@fis.agh.edu.pl.

1. What is StMuRpsUtil
2. Structure
3. Utilities
4. How to use
5. Useful links

What is StMuRpsUtil
StMuRpsUtil is user-friendly utility class which provides a set of post-processing corrections (afterburner) to Roman Pot data stored in StMuRpsCollection class. It has built-in functionalities which expand standard Roman Pot data collection.

Structure
StMuRpsUtil is a ROOT-based class intended to work in the STAR computation environment, as well as local environments e.g. standalone machines. A typical STAR "Maker" format (inheritance from StMaker class) was abandoned in order to gain possibility to run the same code on MuDST files and other storage formats e.g. private picoDST files. The only limitation/requirement is, that Roman Pot data has to be stored in the StMuRpsCollection class.

Usage of StMuRpsUtil invloves creation of single instance of the class at the beginning of analysis, and invocation of StMuRpsUtil::process() and StMuRpsUtil::clear() methods at the beginning and ending of event analysis, respectively. StMuRpsUtil::process() returns pointer to StMuRpsCollection2 class, a mirror class of standard StMuRpsCollection, which contains RP data post-processed using final calibartions. All elements of StMuRpsCollection2 can be accessed in the very same way as of StMuRpsCollection class.

Utlities
StMuRpsUtil provides the following corrections to data:

run-based alignment calibration
(to be implemented) run-based time slewing corrections
(to be implemented) hot strips removal

The following functionalities are available to user:

user can set position of the vertex that is used in reconstruction of proton kinematics
```
StMuRpsUtil::updateVertex(double x, double y, double z)
```
This method should be invoked before StMuRpsUtil::process(). The unit of arguments is meter.

(to be implemented) user can select type of selection criteria (loose, medium, tight): only proton tracks passing cuts at selected level are present in the tracks collection

How to use
► MuDST analysis (working example: /star/u/rafal_s/StMuRpsUtil_tutorial/)

Setup environment to SL16c or newer.
```
starver SL16c
```
Make sure you have the latest definitions of Roman Pot data classes in your StRoot.

Download StMuRpsUtil package from repository.
```
cvs co offline/users/rafal_s/StMuRpsUtil
```
Put downloaded StMuRpsUtil catalogue under StRoot path in your analysis directory.
```
mv offline/users/rafal_s/StMuRpsUtil myAnalysisDir/StRoot/.
```
Edit setup.h file (myAnalysisDir/StRoot/StMuRpsUtil/setup.h) so that only the following line is uncommented.
```
#define RUN_ON_MUDST // set if afterburner is used on MuDST
```
Edit the header file of your analysis maker class.
Add declaration of StMuRpsUtil class before definition of your analysis maker class, and add pointer to StMuRpsUtil object as the element of your analysis maker class.
```
/*...*/
class StMuRpsUtil;
/*...*/

class MyAnalysisMakerClass: public StMaker{
  /*...*/
  StMuRpsUtil*  mAfterburner;
  /*...*/
}
```

Edit the implementation file of your analysis maker class.
Include StMuRpsUtil and StMuRpsCollection2 headers at the beginning.

/*...*/
#include "StMuRpsUtil/StMuRpsUtil.h"
#include "StMuRpsUtil/StMuRpsCollection2.h"
/*...*/

In the analysis maker constructor, create StMuRpsUtil object passing as an argument pointer to StMuDstMaker.

MyAnalysisMakerClass::MyAnalysisMakerClass(StMuDstMaker* maker): StMaker("MyAnalysisMakerClass"){
  /*...*/
  mAfterburner = new StMuRpsUtil(maker);
}

At the beginning of MyAnalysisMaker::Make() method in your analysis maker class, invoke StMuRpsUtil::process() which will provide you post-processed RP data collection. Don't forger to call StMuRpsUtil::clear() at the end of MyAnalysisMaker::Make().

Int_t MyAnalysisMaker::Make( ){
   /*...*/
   StMuRpsCollection2* muRpsColl = mAfterburner->process(); // <-- muRpsColl can be used to get corrected proton tracks etc. 

   /* here analysis of an event */

   mAfterburner->clear(); // <-- critical!!!
   return kStOK;
}

Download RP data calibration files from http://www.star.bnl.gov/~rafal_s/protected/rpsCalibrations2015.tar.gz, unpack it and put exatracted catalogues under myAnalysisDir (you should then have myAnalysisDir/Alignment etc.).
Add the following line:
```
gSystem->Load("StMuRpsUtil.so");
```
to the macro which starts the analysis chain. It ensures that afterburner libraries are accessible.

Edit you job submission XML file so that directories with calibration files extracted from rpsCalibrations2015.tar.gz are included into sandbox.

<SandBox installer="ZIP">
        <Package>
                <!-- Any other files.... -->
                <File>file:./Alignment</File>
                <File>file:./HotStrips</File>
                <!--        etc.         -->
        </Package>
</SandBox>

After the above steps your code should be compilable and should make use of StMuRpsUtil afterburner in MuDST analysis.

If you find any problems using StMuRpsUtil (code does not compile or crashes at execution) you are kindly requested to report it to developers. We kindly ask to not edit the StMuRpsUtil code on your own.

► picoDST analysis (Krakow format) (description to be added)

Useful links
StMuRpsUtil repository in STAR CVS: http://www.star.bnl.gov/cgi-bin/protected/cvsweb.cgi/offline/users/rafal_s/StMuRpsUtil/
Working example of analysis using StMuRpsUtil: /star/u/rafal_s/StMuRpsUtil_tutorial/
Roman Pot data calibration files (run 2015): http://www.star.bnl.gov/~rafal_s/protected/rpsCalibrations2015.tar.gz
StMuRpsCollection documentation (write-up): https://drupal.star.bnl.gov/STAR/system/files/RomanPotsInStEvent_0.pdf
StMuRpsCollection documentation (doxygen): http://www.star.bnl.gov/webdata/dox/html/classStMuRpsCollection.html
Roman Pot alignment description: to be added

Groups:

Printer-friendly version
Login or register to post comments

The STAR experiment

STAR Protected

Detector Sub-systems

User login

Navigation

Group notifications

StMuRpsUtil - Roman Pot data analysis utilities (afterburner)