rnc-lbl / auau200gevrun14

LBNL PicoDst HF Analysis Library (PicoHFLib) - Run14 Au+Au 200GeV

PicoDst HF Analysis Library (PicoHFLib)

Contributors:
- Giacomo Contin ([email protected] ),
- Xin Dong ([email protected]),
- Liang He ([email protected]),
- Michael Lomnitz ([email protected]),
- Mustafa Mustafa ([email protected]),
- Kunsu Oh ([email protected]),
- Jochen Thaeder ([email protected]),
- Guannan Xie ([email protected])

Instituations:
- Lawrence Berkeley National Lab.
- Kent State University,
- University of Science and Technology, China
- Purdue University,
- Pusan National University, South Korea

Introduction and how to:

PicoHfLib
StPicoD0Maker
StPicoHFMaker
StPicoMixedEventMaker

Pico production monitors:

Pico production monitor
Pico D0 production monitor
Pico NPE production monitor

auau200gevrun14's Issues

Pt dependent cuts for StHFCuts

Add pt dependent cut sets for secondary and tertiary vertices - using vectors

Fix TPC PID

Need to generalize TPC PID, and please spell calculate typo

Change the code in StPicoPrescales to raise an exception if file is not found.

The current code hangs in there forever!

Error loading pico list

I am following the instructions on slide 9 of Mustafa's presentation on HF pwg April/2/21015 but I am having and error after loading a picoDst.list:

...
StInfo: StPicoD0AnaMaker - Adding :/project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics/P15ic/166/15166010/st_physics_15166010_raw_0000013.picoDst.root
StInfo: StPicoD0AnaMaker - Adding :/project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics/P15ic/166/15166010/st_physics_15166010_raw_1000002.picoDst.root
StInfo: StPicoD0AnaMaker - Adding :/project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics/P15ic/166/15166010/st_physics_15166010_raw_1000012.picoDst.root
Error in TChain::LoadTree: Cannot find tree with name T in file /project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics/P15ic/166/15166010/st_physics_15166010_raw_0000013.picoDst.root
Error in TChain::LoadTree: Cannot find tree with name T in file /project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics/P15ic/166/15166010/st_physics_15166010_raw_1000002.picoDst.root
Error in TChain::LoadTree: Cannot find tree with name T in file /project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics/P15ic/166/15166010/st_physics_15166010_raw_1000012.picoDst.root

*** Break *** segmentation violation
Root > Function runPicoD0AnaMaker() busy flag cleared

isTofX needs to recalculate beta from TOF and pathlength.

In StHfCuts, Beta is calculated for primary tracks only, so that information will not be available for all global tracks. We need to recalculate beta for tracks using TOF and pathlength information.

StHFMaker: public inheritance and false attempt to force inherited methods to be private

@jthaeder , @mlomnitz

https://github.com/rnc-lbl/auau200GeVRun14/blob/master/StRoot/StPicoHFMaker/StPicoHFMaker.h#L121-L127

I am not sure if this was intended. You are doing a public inheritance from StMaker, what does it mean to override inherited methods as private?!
This stuff doesn't work in C++, overriding members of base class does not change the access.

Also, the comment is confusing because this trick won't disallow derived classes from overloading those methods. Beside, "OVERWRITTEN" is the wrong terminology.

I suggest to move them to be public now. Once C++11 is allowed in STAR you can declare these methods as final override.

Create a STAR CVS mirror of this repo.

Test changes in dev for D0AnaMaker

upload submission script

Use of StPicoCutsBase class

Use the StPicoCutsBase class in NPE Maker (see StPicoHFMaker for usage)

• this provides basic cuts for events / bad runs lists / track cuts / TPC and TOF PID

Errors in StPicoD0AnaMaker.cxx -- Stemming from changes in StHfCuts?

I am getting the following errors when trying to run the AnaMaker:

StRoot/StPicoHFMaker/StHFCuts.h:22:43: error: StPicoCutsBase/StPicoCutsBase.h: No such file or directory
In file included from .sl64_gcc447/obj/StRoot/StPicoD0AnaMaker/StPicoD0AnaMaker.cxx:18:
StRoot/StPicoHFMaker/StHFCuts.h:29: error: expected class-name before '{' token
StRoot/StPicoHFMaker/StHFCuts.h: In member function 'virtual void StHFCuts::init()':
StRoot/StPicoHFMaker/StHFCuts.h:38: error: 'initBase' was not declared in this scope
.sl64_gcc447/obj/StRoot/StPicoD0AnaMaker/StPicoD0AnaMaker.cxx: In member function 'virtual Int_t StPicoD0AnaMaker::Make()':
.sl64_gcc447/obj/StRoot/StPicoD0AnaMaker/StPicoD0AnaMaker.cxx:114: error: 'class StHFCuts' has no member named 'isGoodEvent'
.sl64_gcc447/obj/StRoot/StPicoD0AnaMaker/StPicoD0AnaMaker.cxx: In member function 'bool StPicoD0AnaMaker::isGoodPair(const StKaonPion*) const':
.sl64_gcc447/obj/StRoot/StPicoD0AnaMaker/StPicoD0AnaMaker.cxx:173: error: 'class StHFCuts' has no member named 'isGoodTrack'
.sl64_gcc447/obj/StRoot/StPicoD0AnaMaker/StPicoD0AnaMaker.cxx:173: error: 'class StHFCuts' has no member named 'isGoodTrack'
.sl64_gcc447/obj/StRoot/StPicoD0AnaMaker/StPicoD0AnaMaker.cxx:174: error: 'class StHFCuts' has no member named 'isTPCKaon'
.sl64_gcc447/obj/StRoot/StPicoD0AnaMaker/StPicoD0AnaMaker.cxx:174: error: 'class StHFCuts' has no member named 'isTPCPion'

The errors are coming from sections of the maker that are not in the "user analysis" area. I have the updated versions of all of the code. Was there a change in StHFCuts that did not get propagated to the AnaMaker?

Change all runPico*.C macros to compilable versions

I would like to move the command line in the xml to

root only and compiled macro - using a load lib file before

you don't use submitPicoHFMaker.xml right?

How about adding StPicoDstMaker in here?

Do you guys plan to incorporate the StPicoDstMaker/ module into this repository? It would help with automatic builds and pull requests for StPicoDstMaker.

Show D0 signal fits on the production QA page.

HFCuts failure on failing to open the bad runs list

The error message is opposite to what it should be. I also suggest to mention in the message that you are aborting everything. Or, don't abort and issue a warning.

https://github.com/rnc-lbl/auau200GeVRun14/blob/master/StRoot/StPicoHFMaker/StHFCuts.cxx#L134-L137

Update classdef for PicoD0Event

QAInfo:You are using STAR_LEVEL : SL15c, ROOT_LEVEL : 5.34.09_1 and node : pdsf7
root4star [0]
Processing runPicoD0AnaMaker.C("test.list","test_out")...
strangeArrayNames-arrayNames-NARRAYS = 0
loading of shared libraries done
StInfo: Read in picoDst file /project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics2/P15ic/131/15131052/st_physics_15131052_raw_4000026.picoDst.root
StInfo: Total 1 files have been read in.
StInfo: StPicoD0AnaMaker - Adding :/project/projectdirs/starprod/hft/d0tree/Run14/AuAu/200GeV/physics2/P15ic/131/15131052/st_physics_15131052_raw_4000026.picoD0.root
Warning in TStreamerInfo::BuildCheck:
The StreamerInfo of class StPicoD0Event read from file /project/projectdirs/starprod/hft/d0tree/Run14/AuAu/200GeV/physics2/P15ic/131/15131052/st_physics_15131052_raw_4000026.picoD0.root
has the same version (=1) as the active class but a different checksum.
You should update the version to ClassDef(StPicoD0Event,2).
Do not try to write objects with the current class definition,
the files will not be readable.

Warning in TStreamerInfo::CompareContent: The following data member of
the in-memory layout version 1 of class 'StPicoD0Event' is missing from
the on-file layout version 1:
StThreeVector mKfVertex; //
Warning in TStreamerInfo::CompareContent: The following data member of
the in-memory layout version 1 of class 'StPicoD0Event' is missing from
the on-file layout version 1:
int mNKaonPion; //
Warning in TStreamerInfo::CompareContent: The following data member of
the in-memory layout version 1 of class 'StPicoD0Event' is missing from
the on-file layout version 1:
int mNKaons; //
Warning in TStreamerInfo::CompareContent: The following data member of
the in-memory layout version 1 of class 'StPicoD0Event' is missing from
the on-file layout version 1:
int mNPions; //
Warning in TStreamerInfo::CompareContent: The following data member of
the in-memory layout version 1 of class 'StPicoD0Event' is missing from
the on-file layout version 1:
TClonesArray* mKaonPionArray; //
StQAInfo: =================================================================================
StQAInfo: QAInfo:Chain StChain::bfcChain Ast = 21.57 Cpu = 18.79
StQAInfo: QAInfo:Maker StPicoD0AnaMaker::picoD0AnaMaker Ast = 6.83(31.6%) Cpu = 6.60(35.1%)
StQAInfo: QAInfo:Maker StPicoDstMaker::picoDstMaker Ast = 14.75(68.4%) Cpu = 12.19(64.9%)
StInfo: StPicoD0AnaMaker - writing data and closing output file
StInfo: --------------Error Codes-------------------------
StInfo: nStOK nStWarn nStEOF nStErr nStFatal
StInfo: 14540 0 0 0 0
StInfo: --------------------------------------------------

StMemStat::Summary(calls) Min Aver Max RMS

This is the end of ROOT -- Goodbye

How can i get from your repo the daily picoD0 file?

Hi Mustafa,

How can i get from your repo the daily picoD0 file. I did

git pull origin master

but it gives an error that since i made changes to StPicoD0AnaMaker it
aborts and does not merge. Also the picoD0 daily file list i cannot copy

I am sure i am doing something wrong but can you suggest

Anju,

Library README files

Improve the library readme files. Also add institutions logos.

Test issue

Test issue

Clean github history

Purge, prune and clean (gc) github history. It is bloated with duplicated commits and contaminated with history of text files commits"

Nightly build

Make a nightly build for the entire library. Including:
1- cppcheck
2- valgrind (if possible)
3- Production tests

Issue syncing D0 and pico files

Hello Mustafa and Jochen,
I have been testing loading the D0 trees to combine with soft pions. While running tests i often come across a situation where one of the picoDst's doesn't load but it's associated D0 file is. This leads to the error:
StError: StPicoDstarAnaMaker - !!!!!!!!!!!! ATTENTION !!!!!!!!!!!!!
StError: StPicoDstarAnaMaker - SOMETHING TERRIBLE JUST HAPPENED. StPicoEvent and StPicoD0Event are not in sync.
I tried with the following list as input:

/project/projectdirs/starprod/hft/d0tree/Run14/AuAu/200GeV/physics/P15ic/095/15095020/st_physics_adc_15095020_raw_2500005.picoD0.root
/project/projectdirs/starprod/hft/d0tree/Run14/AuAu/200GeV/physics/P15ic/095/15095020/st_physics_adc_15095020_raw_2000035.picoD0.root
/project/projectdirs/starprod/hft/d0tree/Run14/AuAu/200GeV/physics/P15ic/095/15095020/st_physics_adc_15095020_raw_2000056.picoD0.root

And there is no issue, but if I move the first file (2500005) to the end of list it fails to load and I get the crash. Any ideas?

StHfCuts rejecting almost all tracks

You are using primary momentum in the cuts in StHFCuts::isTPCX. You should use global momentum.
Also, I don't think you should apply an Eta cut.

Npe Library bug

Hello Kunsu and all,
I tried compiling the library with the NPE code in and found a minor bug in StElectronPair.cxx lines 81,85 and 89:

StRoot/StPicoNpeEventMaker/StElectronPair.cxx:81: error: 'lowest' is not a member of 'std::numeric_limits'

I don't think lowest is allowed in the version of the compiler and should switch these to 'min' or something along those lines.

Test modified run script for StPicoD0AnaMaker

Hi Mustafa,

I have no setup, to test the run script properly ... could you have a look?

cheers
Jochen

Show total number of events projection on picoProductionMonitor.

Resolve the tagging vs. updating for getting new filelists.

README file for MixedEvent

It will be nice to have a README file inside the StPicoMixedEventMaker package.
Something like:
https://github.com/MustafaMustafa/auau200GeVRun14/tree/master/StRoot/StPicoD0AnaMaker

make StHFCutsBase class

make StHFCutsBase class in extra lib and use it where appropriate

StPicoCutsBase::isGoodEvent default parameter value

Hi Jochen,

The default parameter should be specified at declaration in the header file.
https://github.com/rnc-lbl/auau200GeVRun14/blob/master/StRoot/StPicoCutsBase/StPicoCutsBase.cxx#L137

Using isGoodEvent without specifying aEventCuts result in a compilation error when default parameter is specified in the implementation file.

Macro to run NPE analysis template

Hello @kunsuoh,

Can you also provide a macro to run the NPE analysis code?
Something similar to https://github.com/rnc-lbl/auau200GeVRun14/blob/master/StRoot/macros/runPicoD0AnaMaker.C

Update DoAnaMaker howto

https://github.com/rnc-lbl/auau200GeVRun14/blob/master/StRoot/StPicoD0AnaMaker/README.md

I think you have to add the https git clone for users which don't have a github account and their pub key uploaded.

change usage of StHFCuts to StPicoCutsBase in StPicoD0AnaMaker

Use the StPicoCutsBase class should be enough - removes dependency to HF library