Difference between revisions of "NIKHEF PAT workshop Oct2006"
(35 intermediate revisions by one other user not shown) | |||
Line 1: | Line 1: | ||
==Pre-Tutorial== | ==Pre-Tutorial== | ||
− | *This pre-Tutorial is copy of the one | + | *This pre-Tutorial is copy of the one given by Ketevi A. Assamagan, adapted to work in 12.0.2 |
*here is a link to other Tutorials by Ketevi, and other usful PAT links: | *here is a link to other Tutorials by Ketevi, and other usful PAT links: | ||
Line 9: | Line 9: | ||
[http://www.nikhef.nl/pub/experiments/atlaswiki/index.php/NIKHEF_PAT_workshop_Oct2006_pre-Tutorial Pre-Tutorial link] | [http://www.nikhef.nl/pub/experiments/atlaswiki/index.php/NIKHEF_PAT_workshop_Oct2006_pre-Tutorial Pre-Tutorial link] | ||
+ | |||
+ | During the workshop we'll browse over the code that we've used in the pre-Tutorial to try and uderstand it. Later we'll focused on the TAG and AAN selection. | ||
==Event Selection== | ==Event Selection== | ||
Line 19: | Line 21: | ||
cp ~ordonez/public/Workshop_files/UserAnalysis/AnalysisSkeleton.h UserAnalysis/. | cp ~ordonez/public/Workshop_files/UserAnalysis/AnalysisSkeleton.h UserAnalysis/. | ||
cd cmt | cd cmt | ||
− | + | gmake | |
cd .. | cd .. | ||
mkdir run2 | mkdir run2 | ||
Line 34: | Line 36: | ||
*First we generate our own TAG. | *First we generate our own TAG. | ||
− | *Still in the new | + | *Still in the new run2 directory: |
<pre> | <pre> | ||
Line 71: | Line 73: | ||
</pre> | </pre> | ||
*And run again. Note that only will run the events that satisfy the query: | *And run again. Note that only will run the events that satisfy the query: | ||
− | EventSelector.Query = "NLooseMuon>0 && NLooseMuon<3 && abs(LooseMuonEta1)<2.5 && abs(LooseMuonEta2)<2.5 && | + | EventSelector.Query = "NLooseMuon>0 && NLooseMuon<3 && abs(LooseMuonEta1)<2.5 && abs(LooseMuonEta2)<2.5 && LooseMuonPt1>10000 && LooseMuonPt2>10000" |
Line 77: | Line 79: | ||
athena Main_topOptions.py | tee athenaOut.log | athena Main_topOptions.py | tee athenaOut.log | ||
</pre> | </pre> | ||
+ | |||
+ | *As an output you have ZmumuOnAOD_withTAGSelection.aan.root. You can have a look at it in root to make sure that the query has been satisfied | ||
*Now we'll rewrite our own AOD containeing only the events selected. We'll do it for only 50 events to gain some time | *Now we'll rewrite our own AOD containeing only the events selected. We'll do it for only 50 events to gain some time | ||
Line 90: | Line 94: | ||
</pre> | </pre> | ||
− | |||
*Comment out | *Comment out | ||
<pre> | <pre> | ||
Line 107: | Line 110: | ||
</pre> | </pre> | ||
− | You'll see | + | You'll see a new Zmumu_new_AOD_withTAGSelection.root file has been created |
==Event Selection using AAN== | ==Event Selection using AAN== | ||
*Follow as in the previous section but useANN=True instead of useTAG=True. In the end you should have the following two files: | *Follow as in the previous section but useANN=True instead of useTAG=True. In the end you should have the following two files: | ||
− | + | ZmumuOnAOD_withAANSelection.aan.root | |
− | + | Zmumu_new_AOD_withAANSelection.root | |
==Interactive Analysis in ATHENA== | ==Interactive Analysis in ATHENA== | ||
Line 126: | Line 129: | ||
*with this line: | *with this line: | ||
<pre> | <pre> | ||
− | include( "mc11. | + | include ( "mc11.004202.ZmumuJimmy.recon.AOD.v11000401.py" ) |
+ | </pre> | ||
+ | |||
+ | *run | ||
+ | <pre> | ||
+ | athena -i Interactive_topO.py | ||
+ | </pre> | ||
+ | *Initialize application manager | ||
+ | <pre> | ||
+ | athena>theApp.initialize() | ||
+ | </pre> | ||
+ | *Run 1 event | ||
+ | <pre> | ||
+ | athena>theApp.nextEvent() | ||
+ | </pre> | ||
+ | |||
+ | *Retrieve Muons from AOD | ||
+ | <pre> | ||
+ | athena>mcon = PyParticleTools.getMuons("MuidMuonCollection") | ||
+ | </pre> | ||
+ | |||
+ | *Other methods are defined in [http://tmaeno.home.cern.ch/tmaeno/PyParticleTools/index.html PyParticleTools] | ||
+ | |||
+ | *Get the first muon | ||
+ | <pre> | ||
+ | athena>m = mcon[0] | ||
+ | </pre> | ||
+ | |||
+ | *Get a list of methods | ||
+ | <pre> | ||
+ | athena>dir(m) | ||
</pre> | </pre> | ||
− | |||
+ | *See pT | ||
+ | <pre> | ||
+ | athena>m.pt() | ||
+ | </pre> | ||
+ | |||
+ | *Get TrackParticle via ElementLink | ||
+ | <pre> | ||
+ | athena>tp = m.track() | ||
+ | athena>dir(tp) | ||
+ | athena>tp.pt() | ||
+ | </pre> | ||
+ | |||
+ | *Retrieve the MC Truth | ||
+ | <pre> | ||
+ | athena>mcc = PyTruthTools.getMcEvents("TruthEvent") | ||
+ | athena>mc = mcc[0] | ||
+ | athena>mc.alphaQCD() | ||
+ | </pre> | ||
+ | *Dump Data Store | ||
+ | <pre> | ||
+ | athena>dumpSG() | ||
+ | </pre> | ||
− | + | *Run next event | |
− | * | ||
<pre> | <pre> | ||
− | + | athena>theApp.nextEvent() | |
− | |||
− | |||
− | |||
− | |||
</pre> | </pre> | ||
− | + | ||
+ | Rewind event counter | ||
<pre> | <pre> | ||
− | + | athena>theApp.seek(0) | |
− | |||
− | |||
− | |||
− | |||
</pre> | </pre> | ||
− | * | + | |
+ | *See pT distribution of muons | ||
<pre> | <pre> | ||
− | + | athena> plot("MuonContainer#MuidMuonCollection","$x.pt()",nEvent=5) | |
− | |||
− | |||
</pre> | </pre> | ||
− | + | ||
==Analysis with the EventView== | ==Analysis with the EventView== | ||
− | * | + | *Due to the time I won't cover this section during the Workshop. But now you should have the basic knowledge requiered to understand it. You can learn about the [https://twiki.cern.ch/twiki/bin/view/Atlas/EventView Event View] here |
== Analysis on the GRID == | == Analysis on the GRID == | ||
Line 193: | Line 238: | ||
</pre> | </pre> | ||
− | *For job | + | *For job status monitoring visit: [http://gridui02.usatlas.bnl.gov:25880/server/pandamon/query Monitoring] |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | + | ==DQ2 Tools== | |
− | + | *For furhter details visit [https://twiki.cern.ch/twiki/bin/view/Atlas/UsingDQ2 DQ2 wiki] | |
− | + | *This tools make your life easier when working on the GRID | |
− | * | ||
<pre> | <pre> | ||
− | + | dq2_ls *Zmumu*AOD* | |
</pre> | </pre> | ||
− | * | + | *Have a look at the files contained in one of the sets |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
<pre> | <pre> | ||
− | + | dq2_ls -f <datasetname> | |
</pre> | </pre> | ||
− | * | + | *change user.GustavoOrdonezSanz.00008.AAN.root for your own |
<pre> | <pre> | ||
− | + | dq2_get -rv -d . user.GustavoOrdonezSanz.00008.AAN.root | |
</pre> | </pre> | ||
− | + | ==Ganga== | |
− | + | We're not going to cover Ganga in this tutorial. But for reference look at:[https://twiki.cern.ch/twiki/bin/view/Atlas/DistributedAnalysisUsingGanga Ganga link] | |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | == | + | ==Final Exercice== |
− | + | To conclude the workshop we're going to try to make our own Higgs to four muon analysis. Taking as an example the ZmumuOnAOD algorithm, make a H4mu one. The purpose of this, is to dig into the Phisics Analysis Tools used in ZmumuOnAOD, and try to have a realistic approach of how is like to make Analyis in Athena. | |
− | |||
− | |||
− | |||
− | |||
− | + | Good luck! | |
− | |||
− | |||
− |
Latest revision as of 09:08, 25 October 2006
Pre-Tutorial
- This pre-Tutorial is copy of the one given by Ketevi A. Assamagan, adapted to work in 12.0.2
- here is a link to other Tutorials by Ketevi, and other usful PAT links:
- Please follow the following pre-Tutorial before the date of the Workshop. It is a very straight forward, tutorial and will help to set up the athena enviroment for the workshop, so we can focus our time in more interested things
During the workshop we'll browse over the code that we've used in the pre-Tutorial to try and uderstand it. Later we'll focused on the TAG and AAN selection.
Event Selection
First we're going to generate a new AAN that we'll use later for selection.
cd ~scratch0/Workshop/ source setup.sh -tag=12.0.2,opt cd 12.0.2/PhysicsAnalysis/AnalysisCommon/UserAnalysis cp ~ordonez/public/Workshop_files/src/AnalysisSkeleton.cxx src/. cp ~ordonez/public/Workshop_files/UserAnalysis/AnalysisSkeleton.h UserAnalysis/. cd cmt gmake cd .. mkdir run2 cd run2 cp ~ordonez/public/Workshop_files/run2/AnalysisSkeleton_jobOptions.py . cp ~ordonez/public/Workshop_files/data2/PoolFileCatalog.xml . cp ~ordonez/public/Workshop_files/data2/mc11.004202.ZmumuJimmy.recon.AOD.v11000401.py . athena AnalysisSkeleton_jobOptions.py | tee athenaOut_AnalysisSkeleton.log
This will generate the file AnalysisSkeleton.aan.root.
Event Selection using TAG
- First we generate our own TAG.
- Still in the new run2 directory:
cp ~ordonez/public/Workshop_files/run2/Main_topOptions.py .
- Edit Main_topOptions.py and uncomment the line
### Create TAG file doWriteTAG=True
athena Main_topOptions.py | tee athenaOut.log
- This will generate the file called ZmumuOnAOD_TAG.root
- Now we are ready to run our algorithms with a TAG selection on the fly. Edit Main_topOptions.py again comment out:
#doWriteTAG=True
- and uncomment
doRunAODwithSelection=True and useTag=True
- Copy:
cp ~ordonez/public/Workshop_files/run2/ZmumuOnAOD_new_jobOptions.py .
- And run again. Note that only will run the events that satisfy the query:
EventSelector.Query = "NLooseMuon>0 && NLooseMuon<3 && abs(LooseMuonEta1)<2.5 && abs(LooseMuonEta2)<2.5 && LooseMuonPt1>10000 && LooseMuonPt2>10000"
athena Main_topOptions.py | tee athenaOut.log
- As an output you have ZmumuOnAOD_withTAGSelection.aan.root. You can have a look at it in root to make sure that the query has been satisfied
- Now we'll rewrite our own AOD containeing only the events selected. We'll do it for only 50 events to gain some time
- Edit again Main_topOptions.py and change:
TheApp.EvtMax = 500 to TheApp.EvtMax = 50
- Comment out
#doRunAODwithSelection=True
- uncomment:
doRewriteAOD=True
and run once more
athena Main_topOptions.py | tee athenaOut.log
You'll see a new Zmumu_new_AOD_withTAGSelection.root file has been created
Event Selection using AAN
- Follow as in the previous section but useANN=True instead of useTAG=True. In the end you should have the following two files:
ZmumuOnAOD_withAANSelection.aan.root Zmumu_new_AOD_withAANSelection.root
Interactive Analysis in ATHENA
- In this exercise, we will see how to do interactive analysis in ATHENA, how to browse the Raw Data in the AOD. Proceed as follows:
get_files Interactive_topO.py
- Edit Interactive_topO.py and replace this line
EventSelector.InputCollections = [ "AOD.pool.root" ]
- with this line:
include ( "mc11.004202.ZmumuJimmy.recon.AOD.v11000401.py" )
- run
athena -i Interactive_topO.py
- Initialize application manager
athena>theApp.initialize()
- Run 1 event
athena>theApp.nextEvent()
- Retrieve Muons from AOD
athena>mcon = PyParticleTools.getMuons("MuidMuonCollection")
- Other methods are defined in PyParticleTools
- Get the first muon
athena>m = mcon[0]
- Get a list of methods
athena>dir(m)
- See pT
athena>m.pt()
- Get TrackParticle via ElementLink
athena>tp = m.track() athena>dir(tp) athena>tp.pt()
- Retrieve the MC Truth
athena>mcc = PyTruthTools.getMcEvents("TruthEvent") athena>mc = mcc[0] athena>mc.alphaQCD()
- Dump Data Store
athena>dumpSG()
- Run next event
athena>theApp.nextEvent()
Rewind event counter
athena>theApp.seek(0)
- See pT distribution of muons
athena> plot("MuonContainer#MuidMuonCollection","$x.pt()",nEvent=5)
Analysis with the EventView
- Due to the time I won't cover this section during the Workshop. But now you should have the basic knowledge requiered to understand it. You can learn about the Event View here
Analysis on the GRID
- In the following examples we'll run on Zmumu datasets.
- This will be a brief introduction to GRID tools for those of you who already went through the trouble of getting a GRID certificate.
- First we should set up DQ2 tools:
source /afs/usatlas.bnl.gov/Grid/Don-Quijote/dq2_user_client/setup.csh.CERN
PANDA/PATHENA
- For further details visit panda wiki
- Setup PATHENA :
cd ~/scratch0/Workshop/12.0.2/PhysicsAnalysis/AnalysisCommon/UserAnalysis cvs update -r UserAnalysis-00-08-13 -A python cvs update -r UserAnalysis-00-08-13 -A share/pathena cvs update -r UserAnalysis-00-08-13 -A share/pathena_util cvs update -r UserAnalysis-00-08-13 -A share/ConfigExtractor.py cd cmt source setup.sh gmake cd ../run
- run pathena
Note: change user.GustavoOrdonezSanz.00008.AAN.root with an appropiate name!!
pathena ZmumuOnAOD_jobOptions.py --inDS 'csc11.005145.PythiaZmumu.recon.AOD.v11004103' --outDS \ 'user.GustavoOrdonezSanz.00008.AAN.root'
- For job status monitoring visit: Monitoring
DQ2 Tools
- For furhter details visit DQ2 wiki
- This tools make your life easier when working on the GRID
dq2_ls *Zmumu*AOD*
- Have a look at the files contained in one of the sets
dq2_ls -f <datasetname>
- change user.GustavoOrdonezSanz.00008.AAN.root for your own
dq2_get -rv -d . user.GustavoOrdonezSanz.00008.AAN.root
Ganga
We're not going to cover Ganga in this tutorial. But for reference look at:Ganga link
Final Exercice
To conclude the workshop we're going to try to make our own Higgs to four muon analysis. Taking as an example the ZmumuOnAOD algorithm, make a H4mu one. The purpose of this, is to dig into the Phisics Analysis Tools used in ZmumuOnAOD, and try to have a realistic approach of how is like to make Analyis in Athena.
Good luck!