Added the angular distributions for W1->cs and W2->ud

PostCutCode/CMakeLists.txt

@@ -1,48 +1,49 @@
 cmake_minimum_required(VERSION 3.8)
 
-
 project(PostCutCode)
 message("announcing PostCutCode")
 set(CMAKE_BUILD_TYPE Debug)
 
 add_compile_options(-std=c++17 -fPIC -Wall -Wextra -Wno-unused-parameter -O2)
 
-# You need to tell CMake where to find the ROOT installation. This can be done in a number of ways:
-#   - ROOT built with classic configure/make use the provided $ROOTSYS/etc/cmake/FindROOT.cmake
-#   - ROOT built with CMake. Add in CMAKE_PREFIX_PATH the installation prefix for ROOT
 list(APPEND CMAKE_PREFIX_PATH $ENV{ROOTSYS})
 
-#---Locate the ROOT package and defines a number of variables (e.g. ROOT_INCLUDE_DIRS)
-find_package(ROOT REQUIRED COMPONENTS RIO EG Minuit RooFit RooFitCore RooStats HistFactory Core RIO Net Hist Graf Graf3d Gpad Tree Rint Postscript Matrix Physics MathCore Thread MultiProc TMVA)
-# find_package(nlohmann_json REQUIRED)
-
-#---Define useful ROOT functions and macros (e.g. ROOT_GENERATE_DICTIONARY)
+find_package(ROOT REQUIRED COMPONENTS RIO EG Minuit RooFit RooFitCore RooStats HistFactory Core Net Hist Graf Graf3d Gpad Tree Rint Postscript Matrix Physics MathCore Thread MultiProc TMVA)
 include(${ROOT_USE_FILE})
 
-# compile options
-
 include_directories(${CMAKE_CURRENT_SOURCE_DIR})
-# link_directories(/cvmfs/sw.hsf.org/spackages7/intel-tbb/2020.3/x86_64-centos7-gcc11.2.0-opt/ey3ft/lib /cvmfs/sw.hsf.org/spackages7/zlib/1.2.13/x86_64-centos7-gcc11.2.0-opt/2wmsk/lib)
 
-#---Create a shared library with generated dictionary
+# Collect source files
 file(GLOB sources "Root/*.cxx")
 
-add_library(PostCutCodeLib SHARED ${sources})
+# Generate the dictionary
+set(DICT_HEADERS PostCutCode/RVecDict.h)
+set(DICT_LINKDEF PostCutCode/LinkDef.h)
+
+ROOT_GENERATE_DICTIONARY(RVecDict
+    ${DICT_HEADERS}
+    LINKDEF ${DICT_LINKDEF}
+)
+
+# Add the shared library and include the dictionary source variable
+add_library(PostCutCodeLib SHARED ${sources} ${RVecDict_dictionary})
+target_include_directories(PostCutCodeLib PUBLIC ${CMAKE_CURRENT_SOURCE_DIR})
 target_link_libraries(PostCutCodeLib ${ROOT_LIBRARIES})
 
-#---Link the exe files
+# Build executables
 file(GLOB_RECURSE exe_files "utils/*.cxx")
 foreach(_exe_file ${exe_files})
     get_filename_component(_theExec ${_exe_file} NAME_WE)
     add_executable(${_theExec} ${_exe_file})
     target_link_libraries(${_theExec} PostCutCodeLib)
 endforeach()
 
-set(ENV{WSCANNER_DATADIR} ${CMAKE_CURRENT_SOURCE_DIR}/data/)
-
-
-# ---------- Generating the setup.sh file
+# Setup script
 set(SETUP ${CMAKE_CURRENT_BINARY_DIR}/../setup.sh)
 file(WRITE ${SETUP} "#!/bin/bash\n")
-file(APPEND ${SETUP} "# this is an auto-generated setup script\n" )
+file(APPEND ${SETUP} "# this is an auto-generated setup script\n")
 file(APPEND ${SETUP} "export PATH=\${PATH}:${CMAKE_CURRENT_BINARY_DIR}/\n")
+
+# Optional: export data env
+set(ENV{WSCANNER_DATADIR} ${CMAKE_CURRENT_SOURCE_DIR}/data/)
+

PostCutCode/PostCutCode/AnalysisWWCR.h

@@ -21,7 +21,46 @@ class AnalysisWWCR: public AnalysisBase
       void run() override;
 
    protected:
+      std::vector<double> ee_corr;
+      std::vector<double> theta_corr;
+      std::vector<double> phi_corr;
+      std::vector<double> chi_corr;
+      std::vector<double> delR_corr;
+      std::vector<double> p_jet_const_1;
+      std::vector<double> p_jet_const_2;
 
+      std::vector<double> p_reco_Jet1;
+      std::vector<double> p_reco_Jet2;
+      std::vector<double> p_reco_Jet3;
+      std::vector<double> p_reco_Jet4;
+      std::vector<double> theta_Jet1_Jet2;
+      std::vector<double> theta_Jet3_Jet4;
+      std::vector<double> theta_Jet1_Jet3;
+      std::vector<double> theta_Jet2_Jet4;
+      std::vector<double> theta_Jet1_Jet4;
+      std::vector<double> theta_Jet2_Jet3;
+      std::vector<double> phi_Jet1_Jet2;
+      std::vector<double> phi_Jet3_Jet4;
+      std::vector<double> phi_Jet1_Jet3;
+      std::vector<double> phi_Jet2_Jet4;
+      std::vector<double> phi_Jet1_Jet4;
+      std::vector<double> phi_Jet2_Jet3;
+      std::vector<double> chi_Jet1_Jet2;
+      std::vector<double> chi_Jet3_Jet4;
+      std::vector<double> chi_Jet1_Jet3;
+      std::vector<double> chi_Jet2_Jet4;
+      std::vector<double> chi_Jet1_Jet4;
+      std::vector<double> chi_Jet2_Jet3;
+
+      std::vector<double> mass_W1;
+      std::vector<double> mass_W2;
+      std::vector<double> p_W1;
+      std::vector<double> p_W2;
+      std::vector<double> e_W1;
+      std::vector<double> e_W2;
+      // std::vector<double> delta_R_W_boson;
+
+      double mc_weight;
 
 };
 

PostCutCode/PostCutCode/LinkDef.h

@@ -0,0 +1,5 @@
+#ifdef __CLING__
+#pragma link C++ class ROOT::VecOps::RVec<float>+;
+#pragma link C++ class ROOT::VecOps::RVec<ROOT::VecOps::RVec<float>>+;
+#endif
+

PostCutCode/PostCutCode/RVecDict.h

@@ -0,0 +1,2 @@
+#pragma once
+#include <ROOT/RVec.hxx>