Monte Carlo Validation

Goal

To write a package to do streamline Monte Carlo Validation. Action items:

The code is structured as an EDAnalyzer, located at the SPRACE CVS repository.

export CVSROOT=/hdacs/shared/cvsrepository

Check out the code - make sure to be inside a directory CMSSW_X_Y_Z/src/SomeAnalyzer beforehand!

cvs co ValidationAnalyzer

cd CMSSW_X_Y_Z/src/SomeAnalyzer cvs ci ValidationAnalyzer

To reconstruct t and tbar, I located the vertex that has a t (tbar) ingoing, and took the outgoing particles of the vertices. However, some of these top vertices also got particles other than the W and b: uubar, gluon. I've (under Thiago's advice) just ignored it and tried to reconstruct the top mass with just W's and b's. Solved: status 2 particles.

To compile with TFileService, take care to put the lines in the BuildFile in the right place (lines above, not below).

Successfully reconstructed invariant mass of the W's from hadronic and leptonic decays.

Successfully reconstructed invariant mass of t and tbar from hadronic and leptonic decays.

For the jet reconstruction, what must be done is the following:
- In the _cfg.py file: use a module that converts HepMC to GenParticles
- In the _cfg.py file: use a module that makes GenJets.
- In the ValidationAnalyzer file: get the GenJets with something like: edm::Handle<GenJetCollection>. Match those jets with the EXTRA partons - that means, if we are dealing with ttbar+2partons, you want to match those jets to the 2 partons, NOT to the ttbar or the b.
- By "match" I really mean "the deltaR in between the parton and the jet is smaller than some number" - at least the size of the jet (usually 0.5), but preferably smaller.