DVM analysis, a step-by-step description

First thing to do is to convert the raw files ***.bin into root files

macros/process filename.bin

This script is just a simple way to use the macro

macros/initbatch.C

which use the library

macros/read_C..so

base on the macro macros/read.C

please, be aware that in this initbatch you are loading the "read" library in a line like

gROOT->ProcessLine(".L read.C") ;

you should verify the path.

If you want to made changes, please do it in the read.C macro and then create the library using the command

.L read.C++ in a root4star session

in the macros/read.C macro you will find a flag "Flag_BNL" which you can set to 1 or 0. Probably you are analyzing DVM data taken in BNL so the flag should be equal to 1 (the default value). Just be aware that the flag exits and it set by hand to the value you choose.

At the end of running process you will have the output files ***.bin.raw.root and ***.bin.res.root.

While taking new data, this step should be done from time to time and then you could copy the root files somewhere, in order to avoid virgo to be full. This step is loooong and boring.

You can change the file names with the time_and_date inside in order to create a list of "ordered files" for your ***.bin.res.root

For this use the simple macro

.x macros/Convert_myLabview_oldLabview.C("yourpathin","yourpathout")

(in a root4star session)

You should check the path for your input files and where do you want to move this files once their names changes (the output path).

At the end you will have the output files rawdata00.bin.res.root, rawdata01.bin.res.root, rawdata02.bin.res.root, rawdata03.bin.res.root ... in the path of your choice.

The whole analysis to obtain the Drift Velocity and the Argon content proposed with Magboltz calculations will be done then

.x macros/ReadDV_openCuts.C(your_data_path,first_year,last_year,first_month,last_month,first_day,last_day,first_file_number,last_file_number)

for example : .x macros/ReadDV_openCuts.C("Data/",2002,2002,2,2,1,2,0,3)

Be aware that the macro search the files in the path you choose and the in directories likes DvData_2002_02_03 (for 3 of February 2002)
In this macro you are using a set of classes which are located in StRoot/StDriftVelocityMonitor. Because of this you should compile with the command

cons

and verify the path in the macro gSystem->Load(".i386_linux24/obj/StRoot/StDriftVelocityMonitor/StDriftVelocityMonitor.so") ;

This macro include the Pressure and temperature correction of the drift velocity. Of course, this is not ready for the new data and you need to create a histogram with the pressure values. The data can be obtained from the monitor web page http://www.star.bnl.gov/STAR/html/ftpc_l/monitor/ftpc_monitor.html

open this page from a starp computer. Click on "Pressure (mbar)", then enter the period of time to be retrieved, then "View raw data" and "Confirm selection" Finally, "Click here to view" (hold shift while clicking!!) and save the file. Then you can use a macro like

macros/ConvertPressure2Histo.C(file_with_input_data, Histogram_Output)

Then, the output histogram with the pressure values should be included in the macro "macros/ReadDV_openCuts.C" instead of the actual file

sprintf(PTINFileName,"%s/PT.root",directory) ;

At the end of running ReadDV_openCuts you will have a plot with the Drift velocity calculated before and after pressure corrections, a inverse pressure and a Argon content plot in function of the time.

Stored data

September 2001-January 2002 : root files in BNL on virgo: /mnt/data1/pclusa8/Data/ProcessDV/DvData_2001*,DvData_2002*

pressure information in BNL on virgo: /mnt/data1/pclusa8/Data/ProcessDV/DataPressure_Sep_Jan.root

January 2003-May 2003 : root files in BNL on virgo: /mnt/data2/pclusa8/Store-Rawdata2003.tar.gz

Created by Maria MORA