• VSM cern twiki
• cern twiki BoostedChargeAsymmetry
• VSM personal ific twiki
• Victionary
• VSM - Computing
• VSM - Boosted Charge Asymmetry Group

VictoriaBoosted - IFIC - ... in construction ...

References and links

• RootCore twiki
• TopRootCore twiki
• TopD3PDBoosted twiki
• Top-antitop resonance search twiki
• Top resonances MC samples
• AtlasOff

Unfolding

• Dilepton_CEADESY_casymm_Winter2013 - CEA/DESY ttbar dilepton charge assymetry
• FBUnfolding - download FBUnfolding package
• README - FBUnfolding README
• Top charge asymmetry meeting
  • TopChargeAsymmetryMeeting - 10/09/2013
  • TopChargeAsymmetryMeeting - 24/09/2013
• summary - summary

To download the FBUnfold package we are going to work at lxplus or ui06. The necessary steps are:

~> cd $t3
~> export CERN_USER=vsanchez
~> kinit vsanchez@CERN.CH
~> VERSION=FBUunfolding-00-00-02
~> svn co svn+ssh://$CERN_USER@svn.cern.ch/reps/atlasinst/Institutes/Saclay/TopAnalysis/ChargeAsymm/FBUunfolding/tags/$VERSION $VERSION
###~> svn co svn+ssh://$CERN_USER@svn.cern.ch/reps/atlasinst/Institutes/Saclay/TopAnalysis/ChargeAsymm/FBUunfolding/trunk FBUunfolding
~> cd FBUunfolding
~> source /afs/cern.ch/atlas/software/dist/AtlasSetup/scripts/asetup.sh 17.5.0   # ls /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/
~> make
~> ###link/copy all MC, data, QCD files to  ${td} to avoid the need to pass paths when running
~> DATAtruth=true_117050_top
~> SIGNALtrc=TRC_117050_top
~> ./run_unfold.sh -t outdir -D $DATAtruth -S $SIGNALtrc -N 12 -M 4 -R 0 -x 0

OLD TopRootCore^? & TopD3PDBoosted - 2013 - analysis

To download the TopRootCore package we are going to work at lxplus. The necessary steps are:

ssh -Y vicsanma@ui06.ific.uv.es
~> source empezarROOT53001.sh   ### This script is below
~> source empezarT3.sh
~> cd $T3
~> export CERN_USER=vsanchez   ### beware! This is the CERN username
~> kinit vsanchez@CERN.CH
~> mkdir TopRootCore_2013 && cd TopRootCore_2013
~> TRCR=TopRootCoreRelease-14-00-08
~> svn co svn+ssh://$CERN_USER@svn.cern.ch/reps/atlasoff/PhysicsAnalysis/TopPhys/TopRootCoreRelease/tags/$TRCR/ $TRCR
~> cd $T3
~> cd TopRootCore_2013/$TRCR/share
~> emacs packages.txt &   ### Replace TopD3PDBoosted-12-01-16 by TopD3PDBoosted-12-01-19
~> ./build-all.sh
~> exit

The empezarROOT53001.sh (or setRoot.sh) script is like this:
###---------------------------------------------------------------------
#!/bin/bash/

cd $HOME

source /afs/cern.ch/sw/lcg/contrib/gcc/4.3/x86_64-slc5/setup.sh
cd /afs/cern.ch/sw/lcg/app/releases/ROOT/5.32.00/x86_64-slc5-gcc43-opt/root/
source bin/thisroot.sh
#source /afs/cern.ch/sw/lcg/app/releases/ROOT/5.32.00/x86_64-slc5-gcc43-opt/root/bin/thisroot.sh

export LD_LIBRARY_PATH=$ROOTSYS/lib:$LD_LIBRARY_PATH
export PATH=$ROOTSYS/bin:$PATH

cd $HOME
###---------------------------------------------------------------------

The empezarROOT53407.sh (or setRoot.sh) script is like this:
###---------------------------------------------------------------------
#!/bin/bash/

cd $HOME
echo =============== ROOT version 5.34.07===================
date +"VIKI %x %r"


#CERN VIRTUAL MACHINE FILE SYSTEM
#export ATLAS_LOCAL_ROOT_BASE=/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase
#alias setupATLAS='source ${ATLAS_LOCAL_ROOT_BASE}/user/atlasLocalSetup.sh'
#setupATLAS
#localSetupROOT 5.34.07-i686-slc5-gcc4.3


source /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.07/x86_64-slc5-gcc43-opt/root/bin/thisroot.sh


cd $HOME
###---------------------------------------------------------------------

Every time you want to use TopRootCore, you must set up this one:

ssh -Y vicsanma@ui06.ific.uv.es 
~> source empezarROOT53001.sh
~> source empezarT3.sh
~> cd $t3/TopRootCore_2013/
~> source RootCore/scripts/setup.sh

Every time you download a new package, you must compile it as:

ssh -Y vicsanma@ui06.ific.uv.es 
~> cd TopRootCore_2013/PackageToCompile/
~> cd /cmt
~> make -f Makefile.RootCore

To check TRC works correctly, we are going to run a litle example. I prefer to run code in a separate run directory, which we will need to setup with links to the data files and the binaries:

ssh -Y vicsanma@ui06.ific.uv.es 
~> source empezarROOT53001.sh
~> source empezarT3.sh
~> cd $t3/TopRootCore_2013/
~> mkdir run && cd run
~> ln -s ../RootCore/data .
~> ln -s ../TopD3PDAnalysis/control .
~> ln -s ../TopD3PDAnalysis/bin .
~> cd ../
~> source RootCore/scripts/setup.sh  
~> 
~> 
~> cd $t3/TopRootCore_2013/run/   ###VERY IMPORTANT!!! Be sure you are at run/ directory when you run the cut-flow
~> cp control/file_list.txt .
~> emacs file_list.txt &   ### add the NTUP*.root name and NEW LINE!!!
~> 
~> 
~> ### https://docs.google.com/spreadsheet/ccc?key=0AsQ0llq8Me3adGxpQzF2dnQ4MXlEaVlHcVpRTUpDdWc&gid=1#gid=0
~> cd $t3/SAMPLES/
~> DATASET=mc12_8TeV.117050.PowhegPythia_P2011C_ttbar.merge.NTUP_COMMON.e1728_s1581_s1586_r3658_r3549_p1562/
~> FILE1=NTUP_COMMON.01313638._000400.root.1
~> FILE2=NTUP_COMMON.01313638._000401.root.1
~> FILE3=NTUP_COMMON.01313638._000402.root.1
~> FILE4=NTUP_COMMON.01313638._000403.root.1
~> FILE5=NTUP_COMMON.01313638._000404.root.1 
~> dq2-get -f $FILE1 $DATASET
~> 
~> 
~> cd $t3/TopRootCore_2013/run/
~> mkdir OUT_CutFlow_mc12_$FILE
~> 
~> CutFlow -f file_list.txt -o OutputHistos_CutFlow.root -mcType mc12 -p control/settings.txt > outTable_CutFlow.txt
~> mv OutputHistos_CutFlow.root OUT_CutFlow_mc12_$FILE/
~> mv outTable_CutFlow.txt OUT_CutFlow_mc12_$FILE/
~> 
~> 
~> CutFlowCorr -f file_list.txt -o OutputHistos_CutFlowCorr.root -mcType mc12 -p control/settings.txt > outTable_CutFlowCorr.txt
~> mv OutputHistos_CutFlowCorr.root OUT_CutFlow_mc12_$FILE/
~> mv outTable_CutFlowCorr.txt OUT_CutFlow_mc12_$FILE/

The TopD3PDBoosted package provides a single-lepton cut flow for the Boosted ttbar (%NTUPBOOST% D3PD's) event selection along with the reconstruction of the selected event to build the ttbar system. The package can be used by packing the D3PD information into the mini-EDM (event data model) defined in BoostedElectronSelection, BoostedMuonSelection, BoostedJetSelection and FatJetSelection classes.

Works!!! 13-May-2013 The package TopD3PDBoosted-12-01-16 (see the official twiki page and the releases), is implemented already at TopRootCore-14-00-07. You can use its cut-flow or its D3PD-maker:

ssh -Y vicsanma@ui06.ific.uv.es
~> BoostedCutFlow --help
~> D3PD2MiniSLBoost --help

To run locally an example with TopD3PDBoosted, you can get a file and save it at your AFS' /tmp/. We need to create two new files (mc_ch_file.txt and Data_file_list.txt) with the full path of the files you want to run over, inside the run/ folder. Go at folder where are stored these files and do:

~> cd $t3/SAMPLES
~> DATASET=mc12_8TeV.117050.PowhegPythia_P2011C_ttbar.merge.NTUP_TOP.e1728_s1581_s1586_r3658_r3549_p1400_tid01213734_00
~> FILE1=NTUP_TOP.01213734._000001.root.1
~> FILE2=NTUP_TOP.01213734._000002.root.1
~> FILE3=NTUP_TOP.01213734._000003.root.1
~> FILE4=NTUP_TOP.01213734._000004.root.1
~> FILE5=NTUP_TOP.01213734._000005.root.1 
~> dq2-get -f $FILE1 $DATASET   ### or ### dq2-get -n1 $DATASET
~> 
~> DATASETmc=XXXXXXXXX
~> cd $DATASETmc
~> echo $PWD > mc_ch_file.txt
~> emacs mc_ch_file.txt &   ### add the NTUP*.root name and NEW LINE!!!
~> cp mc_ch_file.txt $HOME/$t3/TopRootCore_2013/run/
~> 
~> DATASETdata=XXXXXXXXX
~> cd $DATASETdata
~> echo $PWD > Data_file_list.txt
~> emacs Data_file_list.txt &   ### add the NTUP*.root name
~> cp Data_file_list.txt $HOME/$t3/TopRootCore_2013/run/

Every time you do some change in the files, you should compile the package: [@lxplus403]~>

~> cd TopRootCore_2013/
~> ./RootCore/scripts/compile.sh   ###to compile all the packages
~> 
~> cd TopRootCore_2013/FOLDERtoCOMPILE/cmt/
~> make -f Makefile.RootCore

Very important:
BoostedCutFlow.cxx is located under util/ directory. Editing this file one will be able to control the correction applied to the output.

BoostedCorrection.h is a wrapper class for StandardCorrection to be able to enable some correction requested for instance in the spreadsheet.

If you wish to run on the StandardCorection i.e. the thing you need to run on for getting the wanted corrections, please be sure that this is what is done in the file. So that, check which corrections are enable by the TopD3PDBoosted tag that you are using.

Running BoostedCutFlow

Be sure you are at run/ folder and you've done the TRC and ROOT set up:

~> ssh -Y vsanchez@lxplus303.cern.ch
~> source startRoot5-32.sh
~> cd TopRootCore-2012
~> source RootCore/scripts/setup.sh
~> cd run/

To get the cutflow numbers listed in the spreadsheet you need to run the next command line from your run directory (you have to create this ditectory and put "file_list.txt" file, where your input file is listed)

BoostedCutFlow -f MC_file_list.txt -mcType fullSim -rndmType debug -p control/settings.txt -boost > OUT_BoostedCutFlow_Corr.txt
###BoostedCutFlow -f MC_file_list.txt -mcType fullSim -rndmType debug -p control/settings.txt -boost -NoCor> OUT_BoostedCutFlow_NoCorr.txt
###BoostedCutFlow -f MC_file_list.txt -mcType fullSim -rndmType debug -p control/settings.txt -boost -useTruthParticles> OUT_BoostedCutFlow_Corr_TruthPart.txt

BoostedCutFlow's options: -rndmType debug flag: for seedTool.

~> BoostedCutFlow --help
 Usage: BoostedCutFlow [OPTION]
  The list of options are:
  -f <value>        input file list file (default: "control/file_list.txt")
  -h                print this message
  -mcType <value>   set the MC type (fullSim, atlfast2) (default: use directory name.)
  -nf <value>       number of input files to consider (default: 1)
  -n <value>        maximum number of events to process (default: all events)
  -nocache          turn off the TTreeCache, reduce read failures with DPM RFIO
  -o <value>        output ROOT file (default "OutputHistos.root")
  -i <value>        index of the output file name (default nothing)
  -p <value>        parameters file (default: "control/settings.txt")
  -rndmType <value> reseeding strategy (once, file, debug) (default: once)
  -skn <value>      number of events to skip (default: 0)
  -sknf <value>     number of files to skip (default: 0)
  -vin              validate input files (default: false)
  -dataStream       data stream name: Egamma,Muons.   
  -useLooseElectrons  use loose lepton definition.   
  -useLooseMuonss     use loose lepton definition.   
  -useEMJESJets       use EM+JES jet collection (default: LCTopo jets).   
  -useTaus          use taus (default: false) 
  -boost            use boosted branches (default: false)
  -NoCor            DO NOT use correction data/MC in file BoostedCutFlow ONLY (for challenges) (default: true)
  -JetTrigger       use Jet Trigger      (default: false)
  -useTruthParticles  use TruthInfo or not (default: false)
  -useClusteris     use (default: false)
  -findEvent <value> use this to run on one event number, this is good for debuging
  -print             use this ONLY with MiniSLBoostRun2012, this will print the histos to a pdf file (default: false)

To enable correction you can edit the BoostedCutFlow.cxx file.

To enable event listing to print to a file you have to create a directory:

~> mkdir ../RootCore/data/eventList

BoostedCutFlow -f file_list.txt  -boost -mcType fullSim -rndmType debug  -p control/settings.txt

Running D3PD2MiniSLBoost

~> D3PD2MiniSLBoost --help
---- Starting D3PD2MiniSL
 Usage: D3PD2MiniSLBoost [OPTION]
  The list of options are:
  -f <value>        input file list file (default: "control/file_list.txt")
  -h                print this message
  -mcType <value>   set the MC type (fullSim, atlfast2) (default: use directory name.)
  -nf <value>       number of input files to consider (default: 1)
  -n <value>        maximum number of events to process (default: all events)
  -nocache          turn off the TTreeCache, reduce read failures with DPM RFIO
  -o <value>        output ROOT file (default "OutputHistos.root")
  -i <value>        index of the output file name (default nothing)
  -p <value>        parameters file (default: "control/settings.txt")
  -rndmType <value> reseeding strategy (once, file, debug) (default: once)
  -skn <value>      number of events to skip (default: 0)
  -sknf <value>     number of files to skip (default: 0)
  -vin              validate input files (default: false)
  -dataStream       data stream name: Egamma,Muons.   
  -useLooseElectrons  use loose lepton definition.   
  -useLooseMuonss     use loose lepton definition.   
  -useEMJESJets       use EM+JES jet collection (default: LCTopo jets).   
  -useTaus          use taus (default: false) 
  -boost            use boosted branches (default: false)
  -NoCor            DO NOT use correction data/MC in file BoostedCutFlow ONLY (for challenges) (default: true)
  -JetTrigger       use Jet Trigger      (default: false)
  -useTruthParticles  use TruthInfo or not (default: false)
  -useClusteris     use (default: false)
  -findEvent <value> use this to run on one event number, this is good for debuging
  -print             use this ONLY with MiniSLBoostRun2012, this will print the histos to a pdf file (default: false)

Running boosted regime for MC and data:

~> D3PD2MiniSLBoost -f MC_file_list.txt -boost -mcType fullSim -rndmType debug > OUT_D3PD2MiniSLBoost_Corr.txt
~> D3PD2MiniSLBoost -f data_file_list.txt -boost -dataStream VAL > OUT_D3PD2MiniSLBoost_Corr_DATA.txt

Running resolved regime for MC and data:

~> D3PD2MiniSLBoost -f MC_file_list.txt -mcType fullSim -rndmType debug  > OUT_D3PD2MiniSLBoost_Corr.txt
~> D3PD2MiniSLBoost -f data_file_list.txt -dataStream VAL  > OUT_D3PD2MiniSLBoost_Corr_DATA.txt  ###(same as above)

TopRootCore^? & TopD3PDBoosted - 2013 - analysis

To download the TopRootCore package we are going to work at lxplus. The necessary steps are:

ssh -Y vicsanma@ui06.ific.uv.es
~> source empezarROOT53407.sh   ### This script is below
~> source empezarT3.sh
~> cd $T3
~> export CERN_USER=vsanchez   ### beware! This is the CERN username
~> kinit vsanchez@CERN.CH
~> mkdir TopRootCore_2013 && cd TopRootCore_2013
~> TRCR=TopRootCoreRelease-14-00-15
~> svn co svn+ssh://$CERN_USER@svn.cern.ch/reps/atlasoff/PhysicsAnalysis/TopPhys/TopRootCoreRelease/tags/$TRCR/ $TRCR
~> cd $T3
~> cd TopRootCore_2013/$TRCR/share
~> emacs packages.txt &   ### Replace TopD3PDBoosted-X-Y-Z by TopD3PDBoosted-14-00-05
~> ./build-all.sh
~> exit

The empezarROOT53407.sh (or setRoot.sh) script is like this:
###---------------------------------------------------------------------
#!/bin/bash/

cd $HOME
echo =============== ROOT version 5.34.07===================
date +"VIKI %x %r"


#CERN VIRTUAL MACHINE FILE SYSTEM
#export ATLAS_LOCAL_ROOT_BASE=/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase
#alias setupATLAS='source ${ATLAS_LOCAL_ROOT_BASE}/user/atlasLocalSetup.sh'
#setupATLAS
#localSetupROOT 5.34.07-i686-slc5-gcc4.3


source /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.07/x86_64-slc5-gcc43-opt/root/bin/thisroot.sh


cd $HOME
###---------------------------------------------------------------------

Every time you want to use TopRootCore, you must set up this one:

ssh -Y vicsanma@ui06.ific.uv.es 
~> source empezarROOT53407.sh
~> source empezarT3.sh
~> cd $t3/TopRootCore_2013/
~> source RootCore/scripts/setup.sh

Every time you download a new package, you must compile it as:

ssh -Y vicsanma@ui06.ific.uv.es 
~> cd TopRootCore_2013/PackageToCompile/
~> cd /cmt
~> make -f Makefile.RootCore

To check TRC works correctly, we are going to run a litle example. I prefer to run code in a separate run directory, which we will need to setup with links to the data files and the binaries:

ssh -Y vicsanma@ui06.ific.uv.es 
~> source empezarROOT53407.sh
~> source empezarT3.sh
~> cd $t3/TopRootCore_2013/
~> mkdir run && cd run
~> ln -s ../RootCore/data .
~> ln -s ../TopD3PDAnalysis/control .
~> ln -s ../TopD3PDAnalysis/bin .
~> cd ../
~> source RootCore/scripts/setup.sh  
~> 
~> 
~> cd $t3/TopRootCore_2013/run/   ###VERY IMPORTANT!!! Be sure you are at run/ directory when you run the cut-flow
~> cp control/file_list.txt .
~> emacs file_list.txt &   ### add the NTUP*.root name and NEW LINE!!!
~> 
~> 
~> ### https://docs.google.com/spreadsheet/ccc?key=0AsQ0llq8Me3adGxpQzF2dnQ4MXlEaVlHcVpRTUpDdWc&gid=1#gid=0
~> cd $t3/SAMPLES/
~> DATASET=mc12_8TeV.117050.PowhegPythia_P2011C_ttbar.merge.NTUP_COMMON.e1728_s1581_s1586_r3658_r3549_p1562/
~> FILE1=NTUP_COMMON.01313638._000400.root.1
~> FILE2=NTUP_COMMON.01313638._000401.root.1
~> FILE3=NTUP_COMMON.01313638._000402.root.1
~> FILE4=NTUP_COMMON.01313638._000403.root.1
~> FILE5=NTUP_COMMON.01313638._000404.root.1 
~> dq2-get -f $FILE1 $DATASET
~> 
~> 
~> cd $t3/TopRootCore_2013/run/
~> mkdir OUT_BoostedCutFlow_mc12
~> mkdir OUT_ResolvedCutFlow_mc12
~> 
~> BoostedCutFlow -f file_list.txt -boost -mcType fullSim -common > outTable_BoostedCutFlow.txt
~> mv outTable_BoostedCutFlow.txt OUT_BoostedCutFlow_mc12/
~> 
~> ResolvedCutFlow -f list.txt -boost -mcType fullSim -common > outTable_ResolvedCutFlow.txt
~> mv outTable_ResolvedCutFlow.txt OUT_ResolvedCutFlow_mc12/
~> 

The TopD3PDBoosted package provides a single-lepton cut flow for the Boosted ttbar (%NTUPBOOST% D3PD's) event selection along with the reconstruction of the selected event to build the ttbar system. The package can be used by packing the D3PD information into the mini-EDM (event data model) defined in BoostedElectronSelection, BoostedMuonSelection, BoostedJetSelection and FatJetSelection classes.

Works!!! 25-Nov-2013 The package TopD3PDBoosted-14-00-05 (see the official twiki page and the releases), is implemented already at TopRootCore-14-00-15. You can use its cut-flow or its D3PD-maker:

ssh -Y vicsanma@ui06.ific.uv.es
~> BoostedCutFlow --help
~> D3PD2MiniSLBoost --help

To run locally an example with TopD3PDBoosted, you can get a file and save it at your AFS' /tmp/. We need to create two new files (mc_ch_file.txt and Data_file_list.txt) with the full path of the files you want to run over, inside the run/ folder. Go at folder where are stored these files and do (see steps before):

~> cd $t3/TopRoorCore_2013
~> 
~> DATASETmc=XXXXXXXXX
~> cd $DATASETmc
~> echo $PWD > mc_ch_file.txt
~> emacs mc_ch_file.txt &   ### add the NTUP*.root name and NEW LINE!!!
~> cp mc_ch_file.txt $HOME/$t3/TopRootCore_2013/run/
~> 
~> DATASETdata=XXXXXXXXX
~> cd $DATASETdata
~> echo $PWD > Data_file_list.txt
~> emacs Data_file_list.txt &   ### add the NTUP*.root name
~> cp Data_file_list.txt $HOME/$t3/TopRootCore_2013/run/

Every time you do some change in the files, you should compile the package: [@lxplus403]~>

~> cd TopRootCore_2013/
~> ./RootCore/scripts/compile.sh   ###to compile all the packages
~> 
~> cd TopRootCore_2013/FOLDERtoCOMPILE/cmt/
~> make -f Makefile.RootCore

Very important:
BoostedCutFlow.cxx is located under util/ directory. Editing this file one will be able to control the correction applied to the output.

BoostedCorrection.h is a wrapper class for StandardCorrection to be able to enable some correction requested for instance in the spreadsheet.

If you wish to run on the StandardCorection i.e. the thing you need to run on for getting the wanted corrections, please be sure that this is what is done in the file. So that, check which corrections are enable by the TopD3PDBoosted tag that you are using.

Running BoostedCutFlow

Be sure you are at run/ folder and you've done the TRC and ROOT set up:

~> ssh -Y vsanchez@lxplus303.cern.ch
~> source startRoot5-32.sh
~> cd TopRootCore-2012
~> source RootCore/scripts/setup.sh
~> cd run/

To get the cutflow numbers listed in the spreadsheet you need to run the next command line from your run directory (you have to create this ditectory and put "file_list.txt" file, where your input file is listed)

BoostedCutFlow -f file_list.txt -boost -mcType fullSim -common > OUT_BoostedCutFlow.txt
ResolvedCutFlow -f file_list.txt -boost -mcType fullSim -common > OUT_ResolvedCutFlow.txt

To enable correction you can edit the BoostedCutFlow.cxx file.

To enable event listing to print to a file you have to create a directory:

~> mkdir ../RootCore/data/eventList

BoostedCutFlow -f file_list.txt  -boost -mcType fullSim -rndmType debug  -p control/settings.txt

Running D3PD2MiniSLBoost

~> D3PD2MiniSLBoost --help
---- Starting D3PD2MiniSL
 Usage: D3PD2MiniSLBoost [OPTION]
  The list of options are:
  -f <value>        input file list file (default: "control/file_list.txt")
  -h                print this message
  -mcType <value>   set the MC type (fullSim, atlfast2) (default: use directory name.)
  -nf <value>       number of input files to consider (default: 1)
  -n <value>        maximum number of events to process (default: all events)
  -nocache          turn off the TTreeCache, reduce read failures with DPM RFIO
  -o <value>        output ROOT file (default: no output file)
  -i <value>        index of the output file name (default nothing)
  -p <value>        parameters file (default: "control/settings.txt")
  -rndmType <value> reseeding strategy (once, file, debug) (default: once)
  -skn <value>      number of events to skip (default: 0)
  -sknf <value>     number of files to skip (default: 0)
  -vin              validate input files (default: false)
  -dataStream       data stream name: Egamma,Muons.   
  -useLooseElectrons  use loose lepton definition.   
  -useLooseMuonss     use loose lepton definition.   
  -useEMJESJets       use EM+JES jet collection (default: LCTopo jets).   
  -useTaus          use taus (default: false) 
  -common           running on NTUP_COMMON instead of NTUP_TOP
  -boost            use boosted branches (default: false)
  -oldOR            use 2012 updated Overlap Removal (default: false)
  -NoCor            DO NOT use correction data/MC in file BoostedCutFlow ONLY (for challenges) (default: true)
  -JetTrigger       use Jet Trigger      (default: false)
  -useTruthParticles  use TruthInfo or not (default: false)
  -useClusters     use (default: false)
  -useLoose        produce both tight and loose ntuple, otherwise only tight
  -sysOn           fills systematic files
  -mcTruth_noCut   save only MC variables for all events [no cuts applied] (default: false)
  -wJetsLoose      use the loose selection for WJets background estimation (default: false)
  -noTrim           use save also no-Trimmed fatJet info. Stops at C11 (default: false)
  -findEvent <value> use this to run on one event number, this is good for debuging
  -print             use this ONLY with MiniSLBoostRun2012, this will print the histos to a pdf file (default: false)

Running boosted regime for MC and data:

~> D3PD2MiniSLBoost -f file_list.txt -boost -mcType fullSim -common > OUT_BoostedCutFlow_MC.txt
#~> D3PD2MiniSLBoost -f MC_file_list.txt -boost -mcType fullSim -rndmType debug  > OUT_D3PD2MiniSLBoost_Corr.txt
#~> D3PD2MiniSLBoost -f data_file_list.txt -boost -dataStream VAL  > OUT_D3PD2MiniSLBoost_Corr_DATA.txt

Running resolved regime for MC and data:

~> D3PD2MiniSLBoost -f MC_file_list.txt -mcType fullSim -rndmType debug  > OUT_D3PD2MiniSLBoost_Corr.txt
~> D3PD2MiniSLBoost -f data_file_list.txt -dataStream VAL  > OUT_D3PD2MiniSLBoost_Corr_DATA.txt  ###(same as above)

My Analysis

• Add m_lep_m and m_lep_y branches (for reconstructed and corrected leptons).
• For signal and SM samples:
  • switch on dumpTruthParticles=true.
  • Add m_mc_parents and m_mc_children branches.
  • Add all TRUTH information.
  • for resonance (KKGluon or Z prime) → mc_status=62.
    • if ( td->mc_status->at(i) == 62){
  • for the remaining partibles → mc_status=3 || mc_status=155 |1mc_status=52.
    • if ( (td->mc_status->at(i) = 3) || (td->mc_status->at(i) = 155) || (td->mc_status->at(i) == 52) ){
  • (3 for Pythia6, 155 for Jimmy, 52 & 62 for Pythia8)
• For background samples:
  • switch off dumpTruthParticles=true.

MiniSLOrig.h:

• Add definition of all new branches used at .cxx

MiniSLBoosted.cxx:

• Add m_ditop_y branch (for ttbar reconstructed system).
• Add m_hadtop_y and m_leptop_y branches (tops rapidities).
• Change the definition of FatJet's branches. Replace "m_ttree->Branch("fjet_ZZZ", &m_fjet_ZZZ, "fjet_ZZZ[fjet_n]/F");" by m_ttree->Branch("fjet_ZZZ", &m_fjet_ZZZ, "fjet_ZZZ/F");.

MiniSLBoosted.h:

• Add definition of all new branches used at .cxx

???.cxx: Enable ALL corrections

???.cxx: Enable ALL Systematics

My Analysis OLD

• Add m_lep_m and m_lep_y branches (for reconstructed and corrected leptons).
• Enable dumpTruthParticles=true.
• Add m_mc_parents and m_mc_children branches.
• Add all TRUTH information.
• Change the definition of Jet's branches. Replace "m_ttree->Branch("jet_ZZZ", &m_jet_ZZZ, "jet_ZZZ[alljet_n]/F");" by m_ttree->Branch("jet_ZZZ", &m_jet_ZZZ, "jet_ZZZ/F");.

MiniSLOrig.h:

• Add definition of all new branches used at .cxx

MiniSLBoosted.cxx:

• Add m_ttbar_m and m_ttbar_y branches (for reconstructed system).
• Change the definition of FatJet's branches. Replace "m_ttree->Branch("fjet_ZZZ", &m_fjet_ZZZ, "fjet_ZZZ[fjet_n]/F");" by m_ttree->Branch("fjet_ZZZ", &m_fjet_ZZZ, "fjet_ZZZ/F");.
• Add m_hadtop_y and m_leptop_y branches (tops rapidities).

MiniSLBoosted.h:

• Add definition of all new branches used at .cxx

BoostedCutFlow.cxx:

• It contains differents corrections, for example:
  • corr = new BoostedCorrections(NOMINAL); // will enable all correction
  • corr->enableLeptonCorrectionOnly();
  • corr->enableLeptonCorrectionPlus();
  • corr->enableLeptonCorrectionAll();

???.cxx: Enable ALL corrections

???.cxx: Enable ALL Systematics

BoostedRealm by Matteo Franchini

• WEBsvn - BoostedRealm package...
• BoostedRealm.pdf - Guide to use the package

The BoostedRealm framework is designed to perform a full analysis, from D3PD to fynal plots, for the top and ttbar differential cross section with ATLAS data in a self-cosistent and guided environment. This is also a good starting point for other top analysis, in fact it is enought flexible to be integrated with other pieces of code; some bonus codes already in and lots of typical analysis issiues also ready.

This framework's core is constituted by 4 main codes described below and a series of scripts to easy exec programs and for different analysis pourposes, also described in following sections.

BoostedRealm is a TopRootCore (TRC) based framework that garantee a quite up to date and reliable environment. If you are new to TopRootCore you can have look at the TRC Twiki. The framework uses an official TRC package called TopD3PDBoosted developed by Avi Gershon for a semi-leptonic boosted and resolved top event processing (D3PD2MiniSLBoosted ). The D3PD2MiniSLBoosted code to transform D3PD to ntuple appling proper corrections and cuts is part of this package. A local and GRID submission script is provided.

The HistoMakerMaster code, to produce histograms starting from the ntuples and to reconstruct physics objects, is also integrated in the TopD3PDBoosted package. This code also performs the particle level selection (thanks to Ana Ovcharova for her contribution ) and a boosted tagging selection using the Top Template Overlap Method.

The final analysis step is made by 2 codes:

• the one to produce data/MC and other control plots.
• the program for the unfolding and cross section evaluation.

Both these codes recognaize the input files and histograms typology by names in order to speed up and automatize the elaboration and future modifies; apply luminosity correction; rebin and set histograms in editorial format; evaluate statistic and systematic errors; are customizable via parameters file.

Installation

~> ssh -Y vicsanma@ui06.ific.uv.es
~> source empezarT3.sh
~> source empezarROOT53407.sh
~> cd $T3
~> export CERN_USER=vsanchez
~> kinit -5 $CERN_USER@CERN.CH
~> TRCrealm=BoostedRealm

~> svn co svn+ssh://$CERN_USER@svn.cern.ch/reps/atlas-matteo/matteo/tags/BoostedRealm-01-00-00 $TRCrealm
~> cd $TRCrealm
~> source boostedRealm.inst
n   --- Are you running on Bologna Tier3 user interface?
/lustre/ific.uv.es/grid/atlas/t3/vos/SAMPLES/samples_BoostedRealm   --- Please type where to store data files:
vsanchez   --- Please type your grid cern username:
v   --- or any type to continue
y   --- I am writing in your bashrc the line...
y   --- I am writing in your bashrc the line...
a   --- Do you prefer afs root setup (recommeded for non bologna tier3 user)...
"grid password"
14-00-15   --- Type TRC version (Default:14-00-14):
y   --- You chose this version: 14-00-15. Ok?
y   --- Do you want to change the pakage list?    atlasoff/PhysicsAnalysis/TopPhys/TopD3PDBoosted/tags/TopD3PDBoosted-13-01-01
~> 
~> TRC_HOME=/lustre/ific.uv.es/grid/atlas/t3/vicsanma/BoostedRealm/TRCs/TopRootCore12
~> MASTER_REALM=/lustre/ific.uv.es/grid/atlas/t3/vicsanma/BoostedRealm/MasterConsolle
###~> export ROOTSYS=/afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.07/x86_64-slc5-gcc43-opt/root
~> 
~>source $TRC_HOME/RootCore/scripts/setup.sh
~> 
~> ### pwd = /lustre/ific.uv.es/grid/atlas/t3/vicsanma/TopRootCore_BoostedRealm/BoostedRealm/MasterConsolle
~> realm_setup
~> source /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.07/x86_64-slc5-gcc43-opt/root/bin/thisroot.sh
cd /lustre/ific.uv.es/grid/atlas/t3/vicsanma/TopRootCore_BoostedRealm/BoostedRealm/moveFast
svn co svn+ssh://$CERN_USER@svn.cern.ch/reps/atlas-matteo/matteo/BoostedRealm/moveFast/FastJet
cp -r FastJet/ FastJet_original
cp -r TtjetsClassifier/ TtjetsClassifier_original
~> ./AdaptToGrid.sh   ### to have locally working version with the "l" option   ### to set the last version of the histo maker programs 
Type s to copy your codes in the moveFast repository synconized with SVN. Do it to save your work.
Type l to run locally with FastJet, TruthSelection package and HistMakerMaster
l
~> ./FastJetCompile.sh   ###to setup FastJet environment... before to launch any TRC program!!!
~> ./BuildTRC.sh   ### to compile TopRootCore


~> 
~> 
~> 
~> realm_setup   ###Every time you enter your system you need to setup the environment

~> ssh -Y vicsanma@ui06.ific.uv.es
~> source empezarT3.sh
~> source empezarROOT53407.sh
~> cd $T3
~> export CERN_USER=vsanchez
~> kinit -5 $CERN_USER@CERN.CH
~> TRCfolder=TopRootCore_BoostedRealm
~> mkdir $TRCfolder && cd $TRCfolder
~> TRCR=TopRootCoreRelease-14-00-15
~> svn co svn+ssh://$CERN_USER@svn.cern.ch/reps/atlasoff/PhysicsAnalysis/TopPhys/TopRootCoreRelease/tags/$TRCR/ $TRCR
~> cd $TRCR/share
~> 
~> emacs packages.txt & 
### go to TopD3PDBoosted line and substitute it with atlasoff/PhysicsAnalysis/TopPhys/TopD3PDBoosted/tags/TopD3PDBoosted-13-01-01
~> ./build-all.sh

Useful macros and files

To slim D3PD and select only MC branches use the package (using the GRID) slimming_selectMCbranch_GRID.

To slim D3PD and select only MC branches use the package (running locally) slimming_selectMCbranch_LOCAL.

Second step in the slimming: to select the truth objects you want (using the GRID) selecting_MCinfo_GRID.

Second step in the slimming: to select the truth objects you want (running locally) selecting_MCinfo_LOCAL.

To take out given systematics inside a file.root.tgz TakeOut_FilesInsideTGZfiles.sh.

.............................

Victoria Sánchez Martínez
--- IFIC ---
Edificio Institutos de Investigación
Paterna, apartado 22085
E-46071 Valencia (Spain)

.............................