D*0 -> D0 gamma Work

BPHY4 Derivation Code

Link for the Analysing with derivation code tutorial:

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 :

Dmeson.png

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

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng Dmeson.png r1 manage 77.8 K 2016-06-15 - 10:44 ZekeriyaUysal First Result for Dmeson
Edit | Attach | Watch | Print version | History: r51 < r50 < r49 < r48 < r47 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r51 - 2016-06-15 - ZekeriyaUysal
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback