status of pms

24 Jul 01 Simulation Meeting 1

Status of PMS

Sarah Eno


PMS

A fast simulation package for Run II. Similar to the Run I packages CMS (WZ fast simulation), QSIM, the di-boson group fast MC, and the QCD group’s fast simulations.

Structure: similar to D0Reco, in that there is a “package” for each particle type. Each ID group is supposed to be responsible for the package for its particle (in practice, only electron and muon groups pay any attention).

executables available linked to isajet, herwig, pythia, vecbos to allow generation and simulation in one shot. hope to add cms generator (with help from Josh from Michigan)


Structure: PackagesUtility Packages

pmcsMergeKinePkg: does pileup

pmcs_setup: reads the mckine chunks and puts into a nice format for the rest of the packages

pmcs_output: produces ntuples, root trees, etc

pmcs_util: defines classes, utilities of general use to all packages

Particle Packages:

pmcs_chprt: tracking pmcs_vtx: vertex finding pmcs_em: electrons and photons pmcs_jet: jets pmcs_met: missing Et pmcs_muon: muons pmcs_tau: taus

may need in future a b-tag package and a trigger package


pmcsMergeKinePkgcurrent developers: none

status: stable

Current Functionality:

reads in the mckinechunks from pileup events and adds to event record

Plans:

none, stable


pmcs_setupcurrent developers: Sarah Eno and Marco Verzocchi

status: sort of stable


reads mckinechunks and puts into a “pmcsInputChunk”. pmcsInputChunk has useful functions to get lists of various sorts of particles. create “pmcsOutputChunk”, which is used for communication between packages, and for people who don’t want to access the information via d0phyobj chunks (request of the NP group). does generator-level jet finding.

Plans:

move routine that calculates particle position when hits surface of the calorimeter (used in isolation routines and fiducial routines) from pmcs_chprt::buildEvent to pmcs_setup. this routine may also be buggy. debug.

move routine from pmcs_vtx that can pick the primary vertex position from a histogram to pmcs_setup. make it so it can do this for pileup primary vertices as well


pmcs_outputcurrent developers: Stephanie Baffioni (root tree, part-time), Elemer Nagy (root-tree, part-time), Sarah Eno, Marco Verzocchi

status: under active development


root-tuple, paw-tuple output from pmcsOutputChunk works. output only histograms exists (book and fill a special pmcs class that the user can edit). start of root tree output exists. pmcs_em, pmcs_met, pmcs_chprt can write d0physobj chunks. rest of packages can not write chunks. em_analysis does not work correctly off the electron chunks that are created. other *_analysis ntuple maker code not tried yet.

Plans:

first priority is to develop the d0physobj chunks, and use *_analysis packages to produce ntules. This should also allow us to use the official code to produce thumbnails and root trees. However, we need A LOT of help from the ID conveners to achieve this (as we are having trouble finding documentation on these chunks), and active developers for each of the particle packages.


pmcs_utilcurrent developers: Sarah Eno and Marco Verzocchi

status: stable?


provides utilities, like routines to extrapolate particles in a magnetic field, to do cone isolation calculations at the generator-level.

Plans:

isolation routine needs some debugging


pmcs_chprtcurrent developers: Zhong-Min Wang (reasonable fraction of his time), Frank Filthaut (reasonable fraction of his time)

status: under active development


acceptance and smearing of charged particles according to tracker resolutions. current routine is not well-tested. However, new routines that includes proper correlations of resolutions are under development by Zhong-Min. Will include separate parameterizations for electrons, muons, and pions. Frank has also developed code that does acceptances for toy silicon and fiber trackers, for RunIIb studies. can create a track chunk.

Plans:

Zhong Min has found resolutions in 3 of the 5 tracking variables (pt, phi, eta). He needs to extend to include the other 2, and then code his results. We also need to check that the track object that is created being filled in a reasonable way.


pmcs_vtxcurrent developers: Sarah Eno

status: under some development, but needs a real developer


picks primary vertex position as position of hard scatter. can smear z position according to a gaussian

Plans:

move routine that can choose location of primary vertex according to a root histogram from here to pmcs_setup.

find an active developer

make d0physobj output


pmcs_emcurrent developers: Florencia Canelli (wants to move on to other things, part-time)

status: under active development, but could do with more manpower


most advanced packages. smears energies and angles according to calorimeter resolutions. separate resolutions for electrons and photons, CC, and EC, with PS and without. tuned using single particle MC. however, does not seem to properly give resolutions in Z->ee events (could be because of inconsistent use of track information in Z events between pmcs and d0reco)

Plans:

properly simulate dependence of resolution on energy

develop tools to measure resolution from Z->ee data (not MC)

properly fill all variables in the electron chunk


pmcs_em

Shawn Kwang


pmcs_muoncurrent developers: Bing Hu (very part-time), Dave Hedin (very part-time)

status: under active development, but could use some younger manpower


simulates energy and angular resolution of muon chambers. code exists, but perhaps nobody has checked that its output is reasonable.

Plans:

write a d0physobj chunk (hedin)

add pi/K decay (hedin) (maybe into pmcs_setup, though)


pmcs_metcurrent developers: Sarah Eno

status: under some development


two basic algorithms exists, one similar to the old QSIM algorithm, one similar to the old CMS algorithm. neither is well-tested or well-tuned.

Plans:

Marco will tune the cms-like algorithm using run I wmass methods

check filling of met chunk to see if its reasonable


pmcs_jetcurrent developers: none

status: stagnant


takes generator-level jets and smears with a simple smearing. completely untuned and untested.

Plans:

Sudeshna Banerjee, with advice from Laurent DuFlot, may try to take this one.


pmcs_taucurrent developers: none

status: stagnant


none

Plans:

none


Future

Major effort now to put the PMCS output into d0phyobj chunks, and to make sure we understand what the variables in the chunks really are (so we can simulate them properly)

See the documentation::design button on http://www-d0.fnal.gov/~sceno/pmcs_doc/pmcs.htmlfor a document by Marco and myself that attempts to describe

1) the structure of PMCS2) what we think the variables in each chunk represent, and how they should be simulated3) how the resolutions can be obtained for each particle from MC4) How the resolutions can be obtained from real collider data

status of pms

Documents