Kalanand's March 2011 Log

   February 2011   
March 2011
SuMoTuWeThFrSa
  12345
6789101112
13141516171819
20212223242526
2728293031  
   April 2011   

March 1st

EDFilters to write Selector selcted objects

Posted by Sal Rappoccio. There is a new class template for people to play with that can write objects that satisfy the "Selector" for a class to the event record.

https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuidePATSelectors#New_Using_Selector_to_write_sele

This alleviates the burden of having to write your own each time, now you can just declare it to the framework as a plugin and be done.

This is queued for 4.2.x.

March 30th

Re-running HLT on RAW events

From Jim Hirschauer.
I started with the python config for generating QCD MC.

I modified it to do 

process.digitisation_step,
process.L1simulation_step,
process.digi2raw_step,
process.HLTSchedule

The file is here:
/uscms/home/jhirsch/nobackup/IsoTracks_334_data/CMSSW_3_3_4/test/sum09_rehlt.py

1) The cmsDriver command is in the config.

2) This is where you load the HLT config: process.load('HLTrigger/Configuration/HLT_1E31_cff')
   You can presumably put your own file in there.   You may need to recompile!!

3) I only save HLT info:

    outputCommands = cms.untracked.vstring('drop *',
                                           'keep *_hltL1GtObjectMap_*_*',
                                           'keep *_TriggerResults_*_*',
                                           'keep *_hltTriggerSummaryAOD_*_*'),                                  
#    outputCommands = process.FEVTDEBUGHLTEventContent.outputCommands,


To use the "two file solution", you need to specify primary and secondary files:

http://cmslxr.fnal.gov/lxr/search?filestring=&string=secondaryFileNames

Make sure that primary contains the version of the HLT output that you want to use.

March 31st

Neutrino pz calculator

Here is the MET calculator: https://github.com/VPlusJetsAnalyzers/VPlusJets/blob/master/interface/METzCalculator.h

You just run it with something like:
gROOT.ProcessLine('.L METzCalculator.cc+')
METzCalculator = METzCalculator()

# estimate Pz of neutrino
    METzCalculator.SetMET(p4MET)
    METzCalculator.SetLepton(p4lepton)
    if Lepton=="electron":
        METzCalculator.SetLeptonType("electron")
    pzNu = METzCalculator.Calculate()
    p4Nu = TLorentzVector()
    p4OtherNu = TLorentzVector()
    pzOtherNu = METzCalculator.getOther()
You can also check the number of complex solutions with:
if METzCalculator.IsComplex():

Go to February's log


Last modified: Thu Mar 31 12:14:45 CST 2011