(r6) HiggsBb < ATLAS

ATLAS Web>HiggsBb (revision 6) (raw view) ~~Edit~~ Raw edit ~~Attach~~ Print version
%TOC%

---+ %MAKETEXT{"Introduction"}%

Documenting the setup of Higgsbb analysis code used by [[https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/VHbbBaseLine2012][cern.]]

An overview of the software framework can be found [[https://svnweb.cern.ch/cern/wsvn/atlasusr/mbellomo/ElectroweakBosons/trunk/doc/userguide.pdf][here.]]

---+%MAKETEXT{"Running Mode"}%

The code can be run in two different modes. Locally or using PROOF(Parallel ROOT Facility). PROOF allows analysis of a large number of ROOT files in parallel using multiple machines or processor cores.

---+%MAKETEXT{"Running Locally"}%

---++ %MAKETEXT{"Code Checkout"}%

The [[https://svnweb.cern.ch/cern/wsvn/atlasusr/mbellomo/ElectroweakBosons/trunk/AnalysisWZorHbb][code]] can be [[https://twiki.cern.ch/twiki/bin/viewauth/Atlas/SoftwareDevelopmentWorkBookSVN#SVN_Repositories][checked out]] from the correct repository by doing:

<verbatim>
mkdir myVH_research
cd myVH_research
svn co svn+ssh://svn.cern.ch/reps/atlasusr/mbellomo/ElectroweakBosons/trunk/AnalysisWZorHbb/
</verbatim>

Note: you may want to set up [[http://www.farm.particle.cz/twiki/bin/view/ATLAS/Kerberos][kerberos]] first.(I did not require the GSSAPITrustDNS yes line).

---++ %MAKETEXT{"Compiling And Using The Code"}%

(%RED%Note:%ENDCOLOR% for more complete instructions see the README file in the trunk subdirectory) 

To setup the environment:

<verbatim>
cd myVH_research/ElectroweakBosons/trunk
source scripts/setup_lxplus.sh
</verbatim>

This will setup the required environment variables, the Gnu C++ compiler and ROOTCORE.

To download and compile and extra software required do

(%RED%Note:%ENDCOLOR% If the extra code fails to compile due to the MET package you need to remove the code and run the 2012 versions of all the scripts) 


<verbatim>
./scripts/get_allcode.sh
cd SFrame
source setup.sh
cd ..
</verbatim>

Next compile the software by doing

<verbatim>
make
</verbatim>

(You can also clean up by doing =make clean= and =make distclean=)

For the purposes of analysing H->bb decays the code is stored in the =AnalysisWZorHbb/= directory.We now need to compile this code before we can do the analysis.

<verbatim>
cd AnalysisWZorHbb
make
</verbatim>


Now we can run our analysis by doing =sframe_main config/WZorHbb_config_mc_nu_2011.xml= for example. The .xml file is where the input files and settings are all defined(e.g. what corrections to apply to the data) and can be modified or written to suit the needs of the analysis.

(%RED%Note:%ENDCOLOR% every time you open a shell to run the code you need to do steps 1 and 3 again. i.e. setup the environment and make the code. It is also a good idea to run =make clean= before you recompile anything.)


---++ %MAKETEXT{"Producing Plots"}%

The output from running sframe is stored in .root files as common ntuples. To produce plots from these output files we need to use =HistMaker= by doing

<verbatim>
cd macros
make
</verbatim>

to make the required code. To run the code we use the command 

<verbatim>
./RunHistMaker  <path to config file> <runMJ> <path to files output by sframe> <output directory>      -- explain better?
</verbatim>
 
This will output a .root file containing histograms.

If RunHistMaker crashes saying it did not find the cross-section(xsec) for a given number(correlating to a data set) you may need to edit the function =InitCrossSections2011()=  in the file =macros/src/process.cpp=. You will need to add the following code

<verbatim>

allXSecs.push_back(xsection(dataset, cross-section, k factor, 1, "sample name"));  -- find out what 1 is.

e.g.

allXSecs.push_back(xsection(109351, 4.0638E-05, 1, 1,"ZHnunu"));
</verbatim>

Atlas twiki pages can be used to find the cross-section for [[https://twiki.cern.ch/twiki/bin/view/LHCPhysics/CERNYellowReportPageAt7TeV][7TeV data,]] [[https://twiki.cern.ch/twiki/bin/view/LHCPhysics/CERNYellowReportPageAt8TeV][8 TeV data]] and for [[https://ami.in2p3.fr/AMI/servlet/net.hep.atlas.Database.Bookkeeping.AMI.Servlet.Command][general datasets.]] 


For a 2012 event =InitCrossSection2012()= pulls the information from a file called =configs/2012_cross_sections.txt= which can be edited to add the information for your particular dataset
Topic revision: r6 - 2012-10-19 - PaulMullen
ATLAS