D*0 -> D0 gamma Work

BPHY4 Derivation Code

Link for the Analysing with derivation code tutorial:

• https://indico.cern.ch/event/456154/
• https://indico.cern.ch/event/456154/contributions/1122327/attachments/1179820/1707602/ASC_BPhys_DxAOD_Tutorial.pdf

How to run BPHY4 derivation code

First Create your directory

mkdir BPHY4_derivation
cd BPHY4_derivation

Then setup Athena

mkdir 20.1.5.4
cd 20.1.5.4/
asetup 20.1.5.4,here,AtlasDerivation,gcc48

Then setup needed tools

cmt co -r JpsiUpsilonTools-00-02-12 PhysicsAnalysis/JpsiUpsilonTools
cmt co -r TrkVertexAnalysisUtils-00-03-05  Tracking/TrkVertexFitter/TrkVertexAnalysisUtils
cmt co -r TrkVKalVrtFitter-00-07-05        Tracking/TrkVertexFitter/TrkVKalVrtFitter
cmt co PhysicsAnalysis/DerivationFramework/DerivationFrameworkBPhys
cmt co PhysicsAnalysis/DerivationFramework/DerivationFrameworkCore 
cmt co Event/xAOD/xAODBPhys
cmt co -r  xAODMetaDataCnv-00-00-05  Event/xAOD/xAODMetaDataCnv
cmt co -r  xAODMetaData-00-00-06  Event/xAOD/xAODMetaData

Now setup your workarea ..

setupWorkArea.py

Now you need to config and compile all..

cd WorkArea/cmt
cmt bro cmt config
cmt bro gmake

Now you are ready to go , go inside your run directory.

cd ../run

run your derivation

Reco_tf.py --maxEvents=-1 --inputAODFile /disk/lustre/atlas/asc/xAOD/mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_s2601_s2132_r6843_r6264_tid05736188_00/AOD.05736188._003209.pool.root.1 --outputDAODFile output.pool.root --reductionConf BPHY4

How to run on Grid

you can create submit script as below

Content of the submit.sh:

#!/bin/bash


declare -a arr=(
'mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_a782_a787_r6264'
)

for i in "${arr[@]}"
do
   echo "$i"

pathena --allowTaskDuplication  --cmtConfig=x86_64-slc6-gcc48-opt --useNewTRF --trf "Reco_tf.py --inputAODFile %IN --outputDAODFile output.pool.root --reductionConf BPHY4" --inDS=${i}/ --extOutFile DAOD_BPHY4.output.pool.root --outDS user.zuysal.MC3.${i}.BPHY4.r0 --nFiles=-1

done

Ntuple Maker

First Setup the ATLAS software portal and a RootCore release in a working directory;

mkdir WorkDir
cd WorkDir

setupATLAS

rcSetup Base,2.3.22

kinit user@CERN.CH (replace user with you)

Checkout our analysis code from SVN:

svn co svn+ssh://user@svn.cern.ch/reps/atlasinst/Institutes/Birmingham/Run2BPhysics/HadronicDecays (again, replace user with you)

Checkout additional packages:

rc checkout HadronicDecays/share/packages.txt

rc find_packages

Compile the code:

rc compile

Now make a Run directory somewhere and run as RunAnalysis Your_DAOD_datafile_path/ Your_output_Ntuple_filename

mkdir Run
cd Run

RunAnalysis root://eosatlas//eos/atlas/user/z/zuysal/DAOD_MC_Nisan_13TeV/user.zuysal/  charm_mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_a782_a787_r6264_BPHY4.root

Analyzer

Will be updated soon !!!

First Results

First result for 25M event MC15 minbias events :

The left one is the reconstructed mass od D*0 (Kpi gamma) and the right one is the Mass difference between D*0 and D0. The peaks are observed at expected places and they are consistent with the PDG values.

Data Production

Run2 BPHY4 DAOD Data/MC Production Steps

Name	Status	Nevents	Produced	Run	Nfile	Task	Explanation
data15_13TeV.00267358.physics_MinBias.merge.AOD.r7566_p2521	Finished	8665704			533 / 0	http://bigpanda.cern.ch/task/8294852/	Failed ..
data15_13TeV.00267358.physics_MinBias.merge.AOD.r6882_p2358	Finished	8665704	3754546		652 / 347	http://bigpanda.cern.ch/task/8226054/
data15_13TeV.00267358.physics_MinBias.merge.AOD.r6945_p2410	Finished	8665704	795607		418 / 102	http://bigpanda.cern.ch/task/8256197/	DAOD files moved to EOS area.. Waiting for ntup production..
data15_13TeV.00267358.physics_MinBias.merge.AOD.r6849_p2358	Finished	8665704	2953804		831 / 396	http://bigpanda.cern.ch/task/8108517/
mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_a770_a778_r6264	Finished	4981994	4981994		997 / 997	http://bigpanda.cern.ch/task/8266450/	Finished all.. /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set5/ DAOD:361203
mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_a770_a772_r6264	Finished	4988987	4988987		999 / 999	http://bigpanda.cern.ch/task/7956998/	DAOD: 7956998 /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set4/
mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_a771_a772_r6264	Finished	4998981	4953981		1000 / 1000	http://bigpanda.cern.ch/task/8299610/	350 Gb DAOD produced ,Ntuple Production Finished /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set6/ /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set7/
mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_a771_a778_r6264	Finished	4989000	4989000		998/998	http://bigpanda.cern.ch/task/8316192/	/data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set8/ /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set9/
mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_a779_a778_r6264	Ntup prod cont.	4956000	4956000		993/993	http://bigpanda.cern.ch/task/8352944/	Ntuple production contd.. /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set10/
mc15_13TeV.361203.Pythia8_A2_MSTW2008LO_ND_minbias.merge.AOD.e3639_a782_a787_r6264	Grid Prod Finished	19972200	19972200		3995/3995	http://bigpanda.cern.ch/task/8417109/	* Different geometry , must be excluded * 450 of 450 file is moved to eos and Ntuple Production started .. finished sets: /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set11/ /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set12/ /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set13/ /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set14/ /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set15/ /data/data02/userdata/zuysal/13Tev_minbias_Run2_BPHY4_March_set16/

Run2 DAOD Backups

Name	Status	Nevents	produced	task		Explanation	Path on Hep server
user.zuysal.7956998.EXT0._000001.DAOD_BPHY4.output.pool.root	DAOD						/data/data03/userdata/zuysal/13TeV_AOD_Data_Backup/MC

Run1 MC Sample

Name	Status	Nevents	Produced	Run1	Task	Explanation	Path on Hep server
mc11_7TeV.108542.PythiaB_cbmsel1_ChHadr_highpt.merge.AOD.e1039_s1372_s1370_r3044_r2993	HighPt AOD		50 file		http://bigpanda.cern.ch/task/6256567/	AOD	/data/data01/userdata/zuysal/D0_data/high_pt
mc11_7TeV.108532.PythiaB_cbmsel1_ChHadr.merge.AOD.e1039_s1372_s1370_r3044_r2993	LowPt AOD					AOD	/data/data01/userdata/zuysal/D0_data/Low_pt
BPhysLeptonDmeson _new.root	ntuple					ntuple	/data/data01/userdata/zuysal/gladilin
user.zuysal.6212372.EXT0._000001.user.zuysal.Banalysis.root	ntuple					ntuple	/data/data01/userdata/zuysal/gladilin2

Used Samples for Run1

Two Minbias MC sample is used for R1 analysis. Statistics was not enough for us. So we continue with Run2 samples.

• low-pT D mesons (>3 GeV):
  
  mc11_7TeV.108532.PythiaB_cbmsel1_ChHadr.merge.AOD.e1039_s1372_s1370_r3044_r2993/
  
  
• high-pT D mesons (>18 GeV):
  
  mc11_7TeV.108542.PythiaB_cbmsel1_ChHadr_highpt.merge.AOD.e1039_s1372_s1370_r3044_r2993/

Page is Updated by/at : -- ZekeriyaUysal - 2016-06-15 13:30 CET