(Largely DONE) Develop and debug the analysis framework.

Core functionalities are in place and full analysis chain has been tested. In migration to Argonne currently.

(Largely DONE) Trigger performance studies.

Leading two jets 260~270 GeV should do the work. Need some extra checks.

(DONE) Determine the derivation to use.

EXOT2 has all the info needed. If we explore b-trigger in the future we need to make sure trigger matching is available in the derivation.

mkdir ZPrimeto4bSignal

cd ZPrimeto4bSignal

git clone https://:@gitlab.cern.ch:8443/biliu/ANLAnalysis.git

mv ./ANLAnalysis/SignalProduction/ .

rm -rf ANLAnalysis/

cd SignalProduction

source setup_athena.sh

mkdir Zprime4b_500GeV

cd Zprime4b_500GeV

source ../run.sh 320500 123 50 /path/to/your/output/dir 

#argument 1: runumber, runumber - 320000 is the zprime mass
#argument 2: seed, random seed
#argument 3: maxevent
#argument 4: output directory, make sure it exists first

echo "const double n2f=1000000.; grep "cross" Zprime4b_*GeV/log.generate |sed 's/Zprime4b_//g' | sort -n | sed 's?GeV/? ?g'| awk '{print "("$1",n2f*"$6");"}'

# setup athena release

asetup 20.7.8.19,AtlasDerivation,here

#Run reconstruction for TRUTH3

Reco_tf.py --inputEVNTFile=test_Zprime_4b.EVNT.POOL.root 
    --outputDAODFile=test_Zprime_4.AOD.TRUTH3.root --reductionConf=TRUTH3

1   Download tarball from: https://launchpad.net/mg5amcnlo

2   Untar

3   Point the madgraph config file MG5_aMC_v2_6_0/input/mg5_configuration.txt to central versions of FastJet and LHAPDF6:
     lhapdf6 = /share/sl6/lhapdf6/bin/lhapdf-config
     fastjet = /share/sl6/fastjet/bin/fastjet-config

4   Launch MadGraph:
     ./bin/mg5_aMC

5   Install MadAnalysis5 & pythia8: 
     MG5_aMC> install MadAnalysis5
     MG5_aMC> install pythia8

6   Profit.

1   Download tarball from: https://launchpad.net/mg5amcnlo

2   Untar

3   Point the madgraph config file MG5_aMC_v2_6_0/input/mg5_configuration.txt to central versions of FastJet and LHAPDF6:
     lhapdf6 = /afs/cern.ch/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.2.0/x86_64-slc6-gcc48-opt/bin/lhapdf-config
     fastjet = /afs/cern.ch/sw/lcg/external/fastjet/3.2.0/x86_64-slc6-gcc48-opt/bin/fastjet-config

4   Launch MadGraph:
     ./bin/mg5_aMC

5   Install MadAnalysis5 & pythia8: 
     MG5_aMC> install MadAnalysis5
     MG5_aMC> install pythia8

6   Profit.

source /cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase/user/atlasLocalSetup.sh
lsetup rucio

#Replace the p-tag in the example to the latest and greatest one.

rucio ls data17_13TeV:data17_13TeV*EXOT2\.*3372 | grep "CONTAINER" | awk '{print $2}' | sed 's/data17_13TeV://g' > Data13TeV_Main_EXOT2_gridSamples_2017_Dec18th.txt

Note: sometime there are datasets with similar but not exact strings so that please double check the list!

Tip: Do check a given p tag every once a while just in case (Then a diff can tell whether there are new samples.)

The text file generated should go to your MultiBTaggedAlgo/scripts/sampleLists directory as introduced in the framework section.

Luminosity DB tag: OflLumi-13TeV-010

Live fraction trigger: L1_EM24VHI

Physics trigger: HLT_2j225_gsc250_boffperf_split_0eta240_j70_gsc120_boffperf_split_0eta240

LAr EventVeto Tag: LARBadChannelsOflEventVeto-RUN2-UPD4-08

GRL file: data17_13TeV.periodAllYear_DetStatus-v97-pro21-13_Unknown_PHYS_StandardGRL_All_Good_25ns_JetHLT_Normal2017.xml was uploaded successfully to /tmp/lumifiles/392298

--- LumiCalculator      : Livetime     : 23710.3
--- LumiCalculator      : Prescale Weighted Livetime: 23710.3
--- LumiCalculator      : Good LBs     : 406
--- LumiCalculator      : Bad LBs      : 0
--- LumiCalculator      : IntL delivered (ub^-1) : 1.01535e+08
--- LumiCalculator      : IntL after livefraction (ub^-1):  1.00076e+08
--- LumiCalculator      : IntL after LAr fraction (ub^-1):  1.00071e+08
--- LumiCalculator      : IntL recorded after prescale (ub^-1) : 1.00071e+08
--- LumiCalculator      : L1 Prescales: [1107,1195]:1, [1196,1310]:1, [1311,1427]:1, [1428,1562]:1, 
--- LumiCalculator      : HLT Prescales: [1107,1195]:1, [1196,1310]:1, [1311,1427]:1, [1428,1562]:1, 
--- LumiCalculator      : -----------------------------------
--- LumiCalculator      :   LumiCalculator summary
--- LumiCalculator      : -----------------------------------
--- LumiCalculator      :      Total    L1-Acc    L2-Acc    L3-Acc  LiveTime  IntL rec.(ub^-1)  IntL del.(ub^-1)
--- LumiCalculator      :           4157394338         0         04.81438e+06       4.34344e+10       4.54655e+10
--- LumiCalculator      : Total L1/2/3 accept: 41573943380         0         
--- LumiCalculator      : Total L1 livetime trigger accept: 60709860  
--- LumiCalculator      : First Run: 325713    
--- LumiCalculator      : Last Run: 340453    
--- LumiCalculator      : Total L1BeforePresc: 296026634 
--- LumiCalculator      : Total Livetime     : 4.81438e+06
--- LumiCalculator      : Total prescale weighted Livetime: 4.81438e+06
--- LumiCalculator      : Total Good LBs     : 88225
--- LumiCalculator      : Total Bad LBs     : 0
--- LumiCalculator      : Total IntL delivered (ub^-1) : 4.54655e+10
--- LumiCalculator      : Total IntL after livefraction (ub^-1):  4.34513e+10
--- LumiCalculator      : Total IntL after LAr fraction (ub^-1):  4.34344e+10
--- LumiCalculator      : Total IntL recorded (ub^-1) : 4.34344e+10
--- CoolQuery           : Closing database 'CONDBR2
--- CoolQuery           : Closing database 'CONDBR2
--- CoolQuery           : Closing database 'CONDBR2
--- iLumiCalc           : --------------------------------------------
--- TGoodRunsListWriter : GoodRunsList stored as        : ilumicalc_merged_grl_325713.97-340453.1562.xml
--- iLumiCalc           : Real time: 1658.05 s
--- iLumiCalc           : CPU time:  348.2 s
Done
Tue Feb 20 17:37:18 CET 2018

export ATLAS_LOCAL_ROOT_BASE=/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase 
source ${ATLAS_LOCAL_ROOT_BASE}/user/atlasLocalSetup.sh
mkdir FourB
cd FourB
mkdir src build
cd src
asetup AnalysisBase,21.2.10,here
git clone https://github.com/UCATLAS/xAODAnaHelpers.git
git clone https://:@gitlab.cern.ch:8443/biliu/FourBJetAnalysis.git 
git clone https://:@gitlab.cern.ch:8443/biliu/MultiBTaggedAlgo.git 
cd ../build
cmake ../src
make
cd ../
source build/${CMTCONFIG}/setup.sh

export ATLAS_LOCAL_ROOT_BASE=/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase
source ${ATLAS_LOCAL_ROOT_BASE}/user/atlasLocalSetup.sh
kinit your_cern_username@CERN.CH
aklog
mkdir FourB
cd FourB
mkdir src build
cd src
asetup AnalysisBase,21.2.10,here
git clone https://github.com/UCATLAS/xAODAnaHelpers.git
git clone https://:@gitlab.cern.ch:8443/biliu/FourBJetAnalysis.git 
git clone https://:@gitlab.cern.ch:8443/biliu/MultiBTaggedAlgo.git 
cd ../build
cmake ../src
make
cd ../

#Here you may see machine name as centos7 insteed of slac6. This is fine so just ignore it.
source build/x86*/setup.sh

#For test run: you can create a run directory and run analysis code under this directory
#For example mine looks like "build/  runDir/  src/". 
#Under runDir/ you should create an "output/" directory.
#Here is a test command:
run_4b /users/bingxuan.liu/FourB/run/PeriodD/user.biliu.PeriodD.Run331085.physics_Main.Nominal.Part4 /users/bingxuan.liu/FourB/run/HLT_2j225_j70_preselection_noTrigger.txt /users/bingxuan.liu/FourB/run/dataHistograms.txt Inclusive

#(I don’t know if you can read directories from atlas* machine)
#You should ignore ERROR messages related to any TBranches.

cd src
export ATLAS_LOCAL_ROOT_BASE=/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase
source ${ATLAS_LOCAL_ROOT_BASE}/user/atlasLocalSetup.sh
asetup AnalysisBase,21.2.10,here
cd ../
source build/${CMTCONFIG}/setup.sh

source /cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase/user/atlasLocalSetup.sh

lsetup rucio

rucio download mc16_13TeV.361027.Pythia8EvtGen_A14NNPDF23LO_jetjet_JZ7W.deriv.DAOD_EXOT2.e3668_s3126_r9781_r9778_p3340 --nrandom 1

xAH_run.py --config ../src/MultiBTaggedAlgo/data/config_MultiBTaggedAlgo_Nominal_2017.py --files file_you_just_downloaded --isMC --submitDir rel21_test --nevents 100000 direct 

source /cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase/user/atlasLocalSetup.sh
voms-proxy-init -voms atlas
lsetup rucio
lsetup panda
echo data17_13TeV.00333650.physics_Main.deriv.DAOD_EXOT2.f857_m1855_p3372 > test
xAH_run.py --inputRucio --inputList --files test.txt --config $WorkDir_DIR/data/MultiBTaggedAlgo/config_MultiBTaggedAlgo_Nominal_2017.py --force --submitDir gridSubmission/test prun --optGridDestSE=BNL-OSG2_USERDISK --optBatchShellInit 'source $WorkDir_DIR/setup.sh' --optGridOutputSampleName=user.biliu.test.Dec18th.%in:name[1]%.%in:name[2]%.%in:name[3]%.Nominal --optGridNGBPerJob=2

We use TTree (Ntuple) as input format to Histfitter.
Root file format:
       sample_dir/root_file.root:
                                              NomonalTreeName
                                              NomonalTreeName _Sys1Up
                                              NomonalTreeName _Sys1Down
                                              .......

In the Hist Fitter config file define systimatics: jes as an example
jes = Systematic("JES","_NoSys","_JESup","_JESdown","tree","overallSys")