Difference: DanielHaydenDataSamples (8 vs. 9)

Revision 918 Feb 2010 - DanielHayden

Line: 1 to 1
 
META TOPICPARENT name="DanielHayden"

My Data Samples

Line: 59 to 59
 

Setup Area

Changed:
<
<
First of all, you need to setup your work area, this is best done on lxplus. Make your cmthome area and download a package to work in such as UserAnalysis. This is explained how to do so on the coding tips and tricks page.
>
>
First of all, you need to setup your work area, this is best done on lxplus. Make your cmthome area and download a package to work in such as UserAnalysis. This is explained how to do so on the coding tips and tricks page. I am writing here assuming the use of 15.6.1, with the setup -tag=AtlasProduction,15.6.1,setup,32.
 

Get joboptions files to use/modify

Line: 67 to 67
  http://alxr.usatlas.bnl.gov/lxr/source/atlas/Generators/EvgenJobOptions/share/
Added:
>
>
Or more importantly for MC9 now, here:

https://svnweb.cern.ch/trac/atlasoff/browser/Generators/MC09JobOptions/trunk/share

 All I did, was find the Particle and Decay I was interested in, take the joboptions file, and modify it for my own needs. Most of the time this just involves changing the Mass or some Coupling, that is well documented in the joboptions file itself. If you want to do something more drastic however, like make a new sample that has not even been roughly coded before, then I have not done this, but I would suggest you start by looking at multiple joboptions files, such as if you wanted to simulate G->ee and G->mumu, where currently there are only seperate files that do this; Take the two seperate files, and copy the lines which correspond to the decaying of each process into a new file, start building your joboptions that way, and after testing with a well known joboptions to see what you would expect, test yours and see if it makes sense.

Making a Test Sample

Line: 75 to 79
  Go to your run directory such as:
Changed:
<
<
> cd testarea/15.3.0
>
>
> cd testarea/15.6.1
 > cd PhysicsAnalysis/AnalysisCommon/UserAnalysis/run

You need to obtain the following four files via, although the latter most is really just an example file:

> get_files PDGTABLE.!MeV

Changed:
<
<
> get_files csc_genAtlfastTwoStep08_trf.py
> get_files csc_evgenTruthJetMet08_trf.py
> get_files jobOptions.pythia.py

>
>
> get_files csc_genAtlfastTwoStep08_trf.py (Only for MC8)
> get_files csc_evgenTruthJetMet08_trf.py (Only for MC8)
> get_files jobOptions.pythia.py
> get_files Evgen_trf.py (Importantly for MC9)

 
Changed:
<
<
Then you are ready to make your first sample, how first you have a choice. If you want truth, jet, and missing et information in your sample you should use:
>
>
Then you are ready to make your first sample, but first you have a choice. If you want truth, jet, and missing et information in your sample you should use:
  > csc_evgenTruthJetMet08_trf.py 1 1 5000 1 MC8.000001.template.py evgen.pool.root sampleB.AOD.pool.root

Line: 93 to 98
  > csc_genAtlfastTwoStep08_trf.py 1 1 5000 1 MC8.000001.template.py evgen.pool.root atlfast.AOD.pool.root genfast.NTUP.pool.root

Added:
>
>
Most importantly for MC9, you should use Evgen_trf.py, you can type this with the -h command after to see what the following tags mean, but to summarise in order they are; CentreOfMassE / RunNumber (Has to match your jobOption) / Starting Event / Number of Events / Random 10 Digit Seed / JobOptions / Output Filename.

> Evgen_trf.py 7000 0700000 1 5000 1234567812 MC9.0700000.Pythia_Gee01_1000.py Output.pool.root

 Note here, that MC8.000001.template.py is your joboptions file, and 5000, is the number of events you are going to make in your sample. The other file names are your output. I would recommend redirecting these to your tmp directory on lxplus as they can be 100mb big or so. Just remember to send them to linappserv1 or somewhere before you next log off so they don't disappear!

Making Validation Plots & Collecting Information

The next step is to take one of your no doubt already existing analysis programs, and strip it down so all it does it create plots of basic things like the, et spectrum, eta spectrum, and invariant mass, as an example. Apply no cuts, as you just want to see what the sample looks like. Save these plots, and include them with all the other information that you will eventually upload.

Changed:
<
<
Important! When the job finishes running, you will be left with a number of log files, namely csc_evgen08.log and csc_genAtlfast.log, you will want to keep these files as they contain all of your cross section information!!! Also remember to take note of the event generator efficiency, which is usually set to about 0.90, and can be found on the last line of your joboptions file usually. You will also want to look through them to see exactly what version of Pythia you are using, if you haven't taken note of it before now.
>
>
Important! When the job finishes running, you will be left with a number of log files, namely csc_evgen08.log and csc_genAtlfast.log, you will want to keep these files as they contain all of your cross section information!!! Also remember to take note of the event generator efficiency, which is usually set to about 0.90, and can be found on the last line of your joboptions file usually. You will also want to look through them to see exactly what version of Pythia you are using, if you haven't taken note of it before now.
 

Finishing up

 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding RHUL Physics Department TWiki? Send feedback