Giters

c-dilks / dispin

dihadron spin asymmetry analysis

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

dispin

Dihadron spin asymmetry analysis

See also c-dilks/diagrams for several diagrams related to this analysis

Dependencies

coatjava

  • if on ifarm, you have access to coatjava via module load; see documentation on CLAS12 environment setup for details
  • if local, you can obtain coatjava by downloading a release tarball from the clas12-offline-software repository, then setup your environment: 
    export COATJAVA="/path/to/coatjava"
export PATH="${PATH}:${COATJAVA}/bin"
export CLASSPATH="${CLASSPATH}:${COATJAVA}/lib/clas/*"

j2root

  • if you clone this repository with the option --recurse-submodules, it will clone a copy of j2root to the directory deps/j2root
  • build with buildJ2root.sh, otherwise follow directions below
  • cd to this directory and follow the README.md within for compilation
  • this repository provides the ability to write data to root files with groovy
  • notes:
    • added export JAVA_HOME="/usr/lib/jvm/default-java" to ~/.zshrc
    • troubleshooting:
      • scons unable to find JavaH; solution: roll back to java JDK 1.8.0
        • this is because javah is not in 1.11.0 (openjdk-11-jdk), as it was removed in a previous version
        • install openjdk-8-jdk, and optionally uninstall openjdk-11-jdk
        • then ran sudo update-java-alternatives --list to find paths to JDK, followed by sudo update-java-alternatives --set /path/to/java-1.8.0
        • fix symlink usr/lib/jvm/default-java to point to 1.8.0

QADB

  • if you clone this repository with the option --recurse-submodules, it will clone a copy of clas12-qadb to the directory deps/clas12-qadb
  • this repository contains the file quality assurance (QA) database, and uses the groovy accessors to query information in an event loop

brufit

  • build with buildBrufit.sh, otherwise follow directions below
  • build instructions: see deps/brufit readme
  • source ./env.sh will set $BRUFIT, which is needed to compile brufit
    • if you build on ifarm, you may need to run module load [a newer version of cmake] and build with -DCMAKE_C_COMPILER=/path/to/newer/gcc

StringSpinner

  • this is optional
  • build with buildStringSpinner.sh
    • needs Pythia 8.3 (on Arch: pythia8 and python-pythia8)
    • needs gfortran (on Arch: gcc-fortran)
    • see StringSpinner documentation for usage

Ruby

  • Ruby version 3+ is recommended; you may need to build this locally if it's not available on ifarm
  • various Ruby gems are needed; see Gemfile (e.g., run bundle install)
    • the pycall gem is used extensively to access Python tools, such as matplotlib and numpy
      • the corresponding Python packages are also needed
      • calculation of the systematic uncertainties was used as a learning opportunity to try different tools and techniques; unfortunately this means the list of their dependencies has gotten large

RubyROOT

  • RubyROOT is used by some scripts before I switched to using the pycall gem with PyROOT
    • RubyROOT needs to be built separately
    • be sure that its built CMAKE_INSTALL_PREFIX/lib/ruby is in $RUBYLIB

Building

Once all the dependencies are built, run:

make

Usage

  • first call source env.sh
  • create or link diskim and outroot directories
  • now follow the flowchart below; for example, here is the first step:
    • produce outroot files
      • run runDiskim.sh on skim files
        • arguments: skimFile dataStream
          • set datastream to mcrec to analyze MC
        • use slurm or condor wrapper
        • runs skimDihadrons.groovy followed by calcKinematics.cpp
          • output of skimDihadrons.groovy are diskim/*.root, used as input to calcKinematics.cpp
          • output of calcKinematics.cpp are outroot/*.root
          • intermediate files diskim/*.root are automatically deleted

Flowchart

  • These are summary flow charts, and are the most up-to-date version
  • cf. docDiagram.pdf and its Legend for an older (but more complete) version of the flowchart

Main Analysis

flowchart TB
  classDef alg fill:#ff8888,color:black
  classDef data fill:#00aaaa,color:black
  classDef bigdata fill:#00dddd,color:black
  classDef op fill:#00aa00,color:black

  Skims[(skim files)]:::bigdata
  RunDiskim[<strong>slurm.diskim.sh</strong><br/>** runDiskim.sh<br/>**** skimDihadrons.groovy<br/>**** calcKinematics.cpp]:::alg
  Outroots[(outroot files<br/>Data / MC<br/>Inbending / Outbending)]:::bigdata
  Skims --> RunDiskim --> Outroots

  Diagnostics[diagnostics.cpp]:::alg
  DiagnosticsPlots(plots.root):::data
  Outroots --> Diagnostics --> DiagnosticsPlots

  CatTreeMaker[loop_catTreeMaker.rb<br/>catTreeMaker.sh]:::alg
  CatTreesIBOB(catTree files<br/>Data / MC<br/>Inbending / Outbending):::data
  HaddBibending[hadd_makeBibendingCatTrees.rb<br/>ReweightCatTree.C<br/>IndexCatTree.C<br/>StripTweights.C]:::alg
  CatTreesBB(catTree files<br/>Data / MC<br/>Bibending):::data
  CatTreeWeights(catTreeWeights/ TWeights files):::data
  CatTrees[(catTree files<br/>Data / MC<br/>Inbending / Outbending / Bibending)]:::bigdata
  Outroots --> CatTreeMaker
  CatTreeMaker --> CatTreesIBOB --> HaddBibending
  HaddBibending --> CatTreeWeights
  HaddBibending --> CatTreesBB
  CatTreesIBOB --- CatTrees
  CatTreesBB --- CatTrees

  StackCatTreeDists[loop_StackCatTreeDists.rb<br/>StackCatTreeDists.C<br/>latex_StackCatTreeDists.rb]:::alg
  CatTreeStacks(cattreestacks/ images):::data
  CatTrees --> StackCatTreeDists --> CatTreeStacks

  MeanVsMean[loop_TabulateBinMeans.rb<br/>TabulateBinMeans.cpp]:::alg
  MeanVsMeanRoot(meanvmean/ root files):::data
  MeanVsMeanDraw[loop_drawBinMeans.rb<br/>drawBinMeans.rb<br/>latex_BinMeans.rb]:::alg
  MeanVsMeanImages(meanvmean/ images):::data
  MeanVsMeanInputOR{OR}:::op
  Outroots --> MeanVsMeanInputOR
  CatTrees --> MeanVsMeanInputOR
  MeanVsMeanInputOR --> MeanVsMean --> MeanVsMeanRoot --> MeanVsMeanDraw --> MeanVsMeanImages

  FinalFit[finalFit.rb<br/>asymBruFit.C]:::alg
  BruSpin[(bruspin/ files)]:::data
  CatTrees --> FinalFit --> BruSpin
  CatTreeWeights --> FinalFit
Loading

Systematics

flowchart TB
  classDef alg fill:#ff8888,color:black
  classDef data fill:#00aaaa,color:black
  classDef bigdata fill:#00dddd,color:black
  classDef op fill:#00aa00,color:black

  CatTrees[(catTree files<br/>Data / MC<br/>Inbending / Outbending / Bibending)]:::bigdata
  CatTreesOnlyData{{Data files only}}:::op
  CatTrees --> CatTreesOnlyData

  Outroots[(outroot files<br/>Data / MC<br/>Inbending / Outbending)]:::bigdata
  OutrootsOnlyMC{{MC files only}}:::op
  Outroots --> OutrootsOnlyMC

  PolAlg[systematicPolarization.rb]:::alg
  CatTreesOnlyData --> PolAlg

  subgraph systematicBaryonDecayLoop.rb
      BaryonAlg1[systematicBaryonDecay1.cpp]:::alg
      BaryonTrees(sysunc/baryonDecay/*.root):::data
      BaryonAlg2[systematicBaryonDecay2.rb]:::alg
  end
  OutrootsOnlyMC --> BaryonAlg1 --> BaryonTrees --> BaryonAlg2

  subgraph systematicBinMigrationLoop.rb
      MigrationAlg1[systematicBinMigration1.cpp]:::alg
      MigrationRoots(sysunc/binMigration/*.root):::data
      MigrationAlg2[systematicBinMigration2.rb]:::alg
  end
  OutrootsOnlyMC --> MigrationAlg1 --> MigrationRoots --> MigrationAlg2

  RunDiskimRad[runDiskim.sh with datastream<br/>datarad or mcrecrad]:::alg
  RadRoots[(radroot files)]:::bigdata
  RadAlg1[systematicRadiative1.cpp]:::alg
  RadTrees(sysunc/radiative/tree*.root):::data
  RadAlg2[systematicRadiative2.rb]:::alg
  Outroots --> RadAlg1
  RunDiskimRad --> RadRoots --> RadAlg1
  RadAlg1 --> RadTrees --> RadAlg2

  MisPID[systematicsMisPID.nb]:::alg

  RunDiskimPositron[runDiskim.sh with datastream<br/>datapositron]:::alg
  PositronRoots[(positron outroot files<br/>Inbending / Outbending)]:::bigdata
  PositronAlg[systematicNonDIS.cpp]:::alg
  RunDiskimPositron --> PositronRoots
  Outroots --> PositronAlg
  PositronRoots --> PositronAlg

  subgraph Results sysunc/
      PolResult(sysunc/polarization/polarization.dat<br/>sysunc/polarization/polarization.tex):::data
      BaryonResult(sysunc/baryonDecay/*.png<br/>sysunc/baryonDecay/*.json):::data
      MigrationResult(sysunc/binMigration/*.png<br/>sysunc/binMigration/*.npz):::data
      RadResult(sysunc/radiative/*.png<br/>sysunc/radiative/*.json):::data
  end
  PolAlg        --> PolResult
  BaryonAlg2    --> BaryonResult
  MigrationAlg2 --> MigrationResult
  RadAlg2       --> RadResult
Loading

About

dihadron spin asymmetry analysis

Languages

Language:C++ 63.2%Language:C 17.3%Language:Ruby 9.8%Language:Groovy 3.9%Language:Shell 3.0%Language:Python 1.4%Language:Mathematica 1.2%Language:Makefile 0.2%