Analyzer/Comparisons

Setting up Fortran Analyzer (ENGINE) and Replay

• Get copy of Analyzer and Replay directories
  • The ENGINE code has been put in a git repository. If you have not already setup git for the hcana then look at instructions for setting up git and the Git How-to wiki page.
  • Create an "analysis" directory for doing the ENGINE replay.
  • In the "analysis" directory extract the ENGINE code: "git clone git@hallcgit.jlab.org:engine.git" . This creates the subdirectory engine.
  • Use the ENGINE code:
    • go to engine subdirectory
    • git checkout mkj
    • The "mkj" branch has been updated to work on the JLab CentOS 6.2 machines (ifarm). It will work on the CentOS 5.2 machines (ifarml64) but these are being phased out at JLab. The branch also has modifications to make comparisions to HCANA code. The version of the ENGINE is based on the inclusive data analysis.
    • Type "setup cernlib/2005" to set CERN=/site/cernlib/x86_64_rhel6 , CERN_ROOT=/site/cernlib/x86_64_rhel6/2005, CERN_LEVEL=2005 on CentOS 6 machine.
    • Type "make" in engine subdirectory. This will take some time with lots of WARNINGs and should create the executable Linux/bin/engine-replay ( Linux is a subdirectory of the "analysis" directory).

• Get the ENGINE replay files directory
  • Send email to Mark Jones ( jones@jlab.org) to get permission to access the git repository.
  • In the "analysis" directory extract the replay files: "git clone git@hallcgit.jlab.org:jones/replay.git" . This creates the subdirectory replay.

Getting data from the MSS

• Most of examples are for run 52949 which might already be on the cache disk since people are using the file.
• For detailed instructions see new jcache commands)
• An example is : jcache submit default /mss/hallc/daq04/raw/daq04_52949.log.0

Replaying a run with the ENGINE

• Need to have run number 52949 from daq04 experiments on the cache.
• cd to the replay subdirectory in your "analysis" directory
• Type " replay_inter grun=52949 gstop=10000"
  • This will replay 1000 pedestal events and 10000 physics events and create the hbook file ( paw/52949.hbook) and a paw kumac which has the aliases for the histograms ( paw/52949.kumac ).
• To make comparison to HCANA replay, convert the hbook file to root file using h2root.

Replaying a run with the HCANA

• Need to have run number 52949 from daq04 experiments on the cache.
• Assumes that one has checked out the hcana repository and has compile the hcana code.
• Mark Jones created a separate hcana replay directory from the normal hcana/examples subdirectory.
• In the "analysis" directory extract the hcana replay files: "git clone git@hallcgit.jlab.org:jones/hcana_replay.git". This creates the subdirectory replay.
• Replay a run
  • cd to the hcana_replay.
  • If you have compile the hcana code then you should have the environment variable HCANALYZER define, then make symbolic link : " ln -sf $HCANALYZER/hcana hcana"
  • Type "hcana" and get the prompt:"analyzer [0]"
  • For hodo and shower analysis. Type ".x compare_hodo.C(52949,11000)" at the prompt. This replays the 1000 pedestal events and then 10000 physics events. Creates the file compare_hodo_52949.root in the Rootfiles directory.

Compare HCANA to ENGINE for hodoscope and shower

• The root macros assume that the ENGINE root file is in the paw directory that is a subdirectory.
  • Type" ln -sf ../replay/paw paw" to create this subdirectory.
  • Then one can do : "h2root paw/52949.hbook paw/52949.root".
  • Type "root" in the hcana_replay directory
• All plots show the HCANA (red line) plotted over the ENGINE (blue line) in one column. If perfect then just see red line. In other column is the difference between the HCANA and the ENGINE plot. If perfect then zero counts in the histogram.
• For hodoscope ADC hit map:".x Macros/plot_hodo_adc.C(52949)".
  • Scint 1 adc hit map,Scint 2 adc hit map
• For hodoscope TDC hit map:".x Macros/plot_hodo_tdc.C(52949)".
  • Scint 1 tdc hit map,Scint 2 tdc hit map
• For hodoscope focal plane time:".x Macros/plot_fptime.C(52949)"
  • Run 52949 plot of focal plane times
  • One can see that the difference plot shows that the HCANA and ENGINE do not match perfectly at this point. The reason for the difference is explained in the contribution to the April 16th 2013 software meeting.
• For calorimeter ADC hit map:".x Macros/plot_shower_adc.C(52949)".
  • Adc hit map
• For calorimeter individial ADCs :".x Macros/plot_shower_ind_adc.C(52949)".
  • Col A Neg 1-6, Col A Pos 1-6, Col A Neg 7-13, Col A Pos 7-13.
  • Col B Neg 1-6, Col B Pos 1-6, Col B Neg 7-13, Col B Pos 7-13.
  • need to modify the macro to look at columns C and D.

Compare HCANA to ENGINE for drift chambers

• For run 52949, made comparison between HCANA and ENGINE for drift chambers.
  • Plot comparing focal plane positions.
  • Plot comparing focal plane angles.

Hall A Online GUI software

• The code and documentation are available at this public directory.