Recent caf_util modifications

• p21-br-129 April 23, 2010::
  MuonSelector: move muon corrections and fixing to the new processor MuonPtFix
  MuonPtFix: recalculate the MuonPt vs beamspot, redefine medium... ( all this was is MuonSelector before)
  TauPreprocessing, TauCalcNN, TauSelector:
  New TauNN computation. ApplyTauES & NN recomputation are now carried out directly by the TauPreprocessing.
  
  
• p21-br-125 March 1, 2010:
  PhotonSelector: fix some variables init. Does not affect the results.
  TauSelector: add the possibility to create a branch with the leading tau only.
  
• p21-br-124 February 21, 2010:
  TauSelector: change the tau selector default values and the muon veto logic.
  
• p21-br-123 February 2, 2010:
  Include new Unclustered Energy resolution from mET group in METsigAlgo processor.
  
• p21-br-122 January 22, 2010:
  
• p21-br-121 December 14, 2009:
  Add new muon isolation cuts (muisocuts.txt) Add METSig automatic switch between run2a and ru2b.
• p21-br-118, p21-br-119 December 03, 2009:
  ApplyTauES: - first preliminary version of E/P corrections for tau (F. Couderc) - move E recalculation to before function so the input branch is also corrected - new muon isolation definitions (O. Brandt) - bug fix in p21-br-119 (ApplyTauES, p21-br-118 was not loading in cafe)
• p21-br-117 November 18, 2009:
  Electron/Photon Selectors:
  - new emID definitions
  - enable different emID criteria in CC and EC
  Muon Selector:
  - introduce new muon track definitions (cf caf_util/configs/mutrk.cuts).
  - current track defs are still default (new defs accessible via configs ==> MuonSelector.TrkQualityVersion: 2 )
  
  
• p21-br-116 October 16, 2009:
  Changes in MuonSelector: medium definition is now propagated to the TMBMuon candidate (so the TMBMUon::isMedium() function returns the correct value). Fix the bugged isCosmicT value (also propagated to he TMBMuon candidate), now the flags isCosmicT == isCosmic. Changes in ReComputeMET: apply tau corrections as soon as the tau branch is provided (and remove the config parameter).
• p21-br-114 October 02, 2009 Change in the electron selector: hidden pT cut remove, remove duplicate cuts, change the weird pT calculation to the sandard one.
• p21-br-112 July 17, 2009 Allow access to new EMid variables.
• p21-br-111 July 07, 2009 Bug fix in JetSelector: dR[jet-object] (to remove a jet matching a given object) was incorrectly calculated.
• p21-br-110 June 19, 2009 Bug fix in ApplyTauES processor when dealing wih p17 MC.
• p21-br-109 June 17, 2009 Bug fix in MetSig (smearing of the UE)...
• p21-br-108 June 16, 2009 Add smearing of the inclustered energy in the METsig processor (bug correction)
• p21-br-107 June 16, 2009 Add smearing of the inclustered energy in the METsig processor (do not use, this version is buggy).
• p21-br-106 June 12, 2009 Change the default track selection TauSelector to trackloose_tau
• p21-br-105 April 09, 2009 Autodownload track branch in ApplyTauES to access tracks info.
• p21-br-104 April 06, 2009 Propagate the new tau energy scale architecture to the ReComputeMET processor + ApplyTauES: EM scale is now applied also for runIIa data.
• p21-br-101 March 23, 2009 Add a line to ReComputeMET to save METL1 in output processor.
• p21-br-100 March 20, 2009 Add FakeTrackKiller processor.
• p21-br-99 February 14, 2009 Add protections in ApplyTauES, removes hardcoded nsmt>1 cut in TauSelector (use trk selection instead).
• p21-br-98 February 9, 2009 Update MET significance processor with p20 jet resolution functions.
• p21-br-97 December 23, 2008 ResSelector: add possibility to save selected objects in output branches
• p21-br-96 December 18, 2008Add processor to copy b-tag references. Adapt ResSelector to the new approach for electron track matching use
• p21-br-95 December 12, 2008add ApplyTauES processor
• p21-br-94 December 10, 2008:Set back the actas status in RescomputeMET processor for jets
• p21-br-93 December 3, 2008: Modify trk selector and trk selector so they can use mu trk id definitions.
• p21-br-92 November 29, 2008: correctauogic used for the track match flag in Electron Selector. Needed for the dZ(electron, PV) cut.
• p21-br-91 November 28, 2008: correct ELectronSelector behaiviour for the new emid definitions
• p21-br-90 November 24, 2008: Modified ICRElectronSelector.
• p21-br-89 November 24, 2008: Correct tau correction in ReComputeMET. Moved dZ(lepton,PV) selections to MuonSelector and ElectronSelector.
• p21-br-88 November 18, 2008: Use EOP flag from TMBEMcluster in vertex and resonance selector. Make uniform using of PI and TWOPI in ELectronSelector.
• p21-br-87 November 17, 2008:Added track match info with interface to TMBEMCluster. Add possibility to cut on number of track from PV 0
• p21-br-86 November 16, 2008: ApplyJES: update interfaces for the new jetcorr version p21-br-11. MuonSelector: add possibility to cut on the global fit chisq; update defenition of muon and tracks matched to met muons. ReComputeMET: add jet/em removal to this processor; add 2d muon branch for the final state muons (these muons could have looser definition then the standard "d0correct" ones).
• p18-br-74, p21-br-85 November 08, 2008: modify MuonSelector so it uses the new medium definition (only run2b for now).
• p21-br-82 October 22, 2008: modified ElectronSelector.hpp
• p18-br-73 October 20, 2008: Updated ElectronSelector.
• p21-br-81 October 10, 2008: Updated ElectronSelector to analyze p17 data.
• p18-br-71 October 2, 2008: Correct bug in COMPONENT file. EMLikelihood processor was not compiled.
• p21-br-80 October 1, 2008: update of ICR electron selector
• p18-br-70 September 30, 2008: Remove debugging code from Electron selector.
• p21-br-79 September 23, 2008: correct LumiZeroRemoval so it matches the definition used by central caf counting
• p18-br-69 September 23, 2008: Increase the cut to instlumi < 1E30cm-2s-1 (instead of being ==0) because most zero lumi evts were still there (coherent with new alpgen weights)
• p21-br-78 September 12, 2008: METSelector: add ossibility to apply cut on tau corrected MET from p18
• p18-br-68 September 10, 2008: Modified ElectronSelector to work with new electron ID
• p21-br-77 June 24, 2008: Bug fix in TauSelector
• p21-br-76 June 18, 2008: modified caf_util/PhotonSelector.hpp
• p21-br-75 June 17, 2008: modified PhotonSelector.cpp
• p21-br-74 June 17, 2008: modified PhotonSelector
• p18-br-67, p21-br-73 June 10, 2008: Fix partial read bug in LumiZeroEventRemoval.
• p21-br-72 June 7, 2008: Added ICR ele selector to rootrc
• p21-br-71 June 6, 2008: Modified ICR electron selector, photon selector
• p21-br-69 June 5, 2008: Avoid root-bug in StatBranch
• p21-br-68 June 1, 2008: Added support for 3D histograms
• p21-br-67 May 31, 2008: Remove beam spot file from caf_util. MuonSeletor now use the beamspot file from the beamposition package. This package is updated automatically on the regular basis. The latest version always the best one.
• p21-br-66 May 30, 2008: Update the beam spot file.
• p21-br-65 May 27, 2008: ElectronSelector.cpp: removed EM clusters with |ID|>11
• p21-br-64 May 3, 2008: Added back improved statistics for the THIRD time.
• p21-br-63 Apr 29, 2008: modified electron selection code
• p18-br-66, p21-br-62 Apr 15, 2008: Muon cuts defenition: make various names in these files consistent with how the muon spc files are named
• p21-br-61 Apr 7, 2008: Added back improved statistics AGAIN
• p18-br-65, p21-br-60 Mar 21, 2008: ReComputeMET: fix bug in JES progagation for systematics
• p21-br-59 Mar 19, 2008: modified isolation for electron & photon selectors
• p18-br-64 Mar 19, 2008: modified isolation for photon selectors
• p21-br-58 Mar 18, 2008: added electron ANN-7 and modified isolation to ElectronSelector
• p21-br-57 Mar 12, 2008: update of PhotonSelector: added ANN4&5
• p18-br-63 Mar 12, 2008: Modification of PhotonSelector
• p18-br-62, p21-br-56 Mar 11, 2008: Simple Met branches are now created containing the correct, bug-free MET.
• p21-br-55 Feb 6, 2008:MEtSigAlg: synchronized the p21 branch with latest updates to the p18 branch. The updates are:
  
  • Latest muon resolutions are now used
  • The old muon resolutions had a bug in the p17 post-shutodown resolutions, p17 mu+jets should look bad
  • The old code used data jet resolutions both for data and MC The new one uses data jet reso for data and MC for MC
• p21-br-54 Feb 22, 2008: ElectronSelector: calculate the new EMid variables only in the case when they are really used.
• p18-br-61 Jan 22, 2008: Correct behaviour of METsignAlg: MET is being calculated inside the code rather then taken from caf_util::MET. Secondly unclustered energy is being set to 0.
• p21-br-52 Jan 16, 2008: Update beam spot file
• p21-br-51 Jan 9, 2008: Add missing include in EMLikelihood.cpp
• p18-br-60, p21-br-50 Jan 9, 2007: ApplyJES: added cap on muon correction + bug fix for no output muon branch
• p21-br-49 Dec 29, 2007: Modified ElectronSelector with new electron definitions
• p21-br-47 Dec 18, 2007: Add tau correction to the ReComputMET
• p18-br-59, p21-br-43 Dec 11, 2007: Add classes to save event variables to the branch. ApplyJES: by default muon from ApplyJES will not be saved. Fixed small bug in new part of FindDuplicateEvetns code
• p18-br-58, p21-br-41 Dec 8, 2007: Updated ReComputeMet -- changed jet/muon matching, a muon inside a jet cone can be used for the MET muon correction if wasn't used for the JESMU correction of the jet. Also updated the statistics of all obect used for MET correction. Added stat output for the CleanICRJets processor 2. Added support for automatically recognising the MC type in FindDuplicateEvents and access the right MC duplicate event table from mc_dup_evt.
• p18-br-56, p21-br-39 Nov 10, 2007: throw exception in ApplyJES for zero lumi events
• p18-br-55, p21-br-38 Nov 10, 2007: Added new processor for removal of jets overlapping with muons (used in single top QCD selection), minor updates to the LumiZeroEventRemoval - this new processor was broken in last p21 tag.
  Added more statistics for what obejcts were used to recompute MET.
  MuonSelector, ElectronSelector, JetSelector, METSelector, VertexSelector, : updated the stat handling - now the processor name is specified together with the cut.
  Added some statistics (average number of jet/muons/elecrons used etc) for the ReComputeMET processor.
• p18-br-54, p21-br-36 Nov 7, 2007: Add processor to remove MC events with luminosity 0
• p18-br-53, p21-br-35 Oct 30, 2007: Update Tau eneregy correction
• p18-br-52, p21-br-34 Oct 25, 2007: Syncronized ApplyJES with p18 branch (muons used for JES saved in separate branch) + fixed uninitilized counter. Improvements to the HistoHolder - Sumw2 to all histos + support for X and Y titles.
• p18-br-51, p21-br-33 Oct 24, 2007: Add expert options for JES/smearing propagation in ReComputeMET
• p21-br-31 Sep 27, 2007: Add possibility to use p17 NN for Run2b data in TauSelector
• p18-br-48, p21-br-30 Sep 17, 2007: Update MET significance with new jet and muon resolution
• p18-br-47, p21-br-29 Aug 14, 2007: Updated the default muon IDs to version 2 (again) as agreed on with the MUID group - changed default version number to 2 in the MuonSelector. Minor fixes in ApplyJES. Added processor for removing evetns with jet in the ICR (i.e. remove ICR horns).
• p18-br-46, p21-br-26 July 31, 2007: Small fixes to ApplyJES. Added muon ID defs to the official muon files (from top_cafe). New informative messages in the MuonSelector.
• p18-br-45, p21-br-25 July 25, 2007: Added feature in HistoHolder (possiblity to add comments to the histgrams).
• p18-br-44, p21-br-24 July 4, 2007: Remove doJESMU switch from ApplyJES. This will force JESMU recalculation in all cases. Change GetChargedTrack to getPtrChpSpatial track match in ElectronSelector.
• p21-br-23 June 23, 2007: Update beam spot file
• p18-br-43, p21-br-22 June 7, 2007: ApplyJES modifications. The bugs affect the JESMU correction of jets with muons inside:
  • if using AddDataErrorsToMC, all such MC jets have had the DATA JES applied :(
  • if using AddDataErrorsToMC, the uncertainty of the JES MET correction was negative
  • the hasMU flag was not set correctly when matching jets with muons
  • the JES muon correction was not propagated to the MET (this was fine in jetcorr p18-br-06, but not later)
  New features in ApplyJES:
  • (Since earlier) By default, ApplyJES will automatically identify the input (p17/p20, data/MC) and apply the appropriate JES correction.
  • ApplyJES now creates a summary with statistics of "problematic" events (some events have no inst. lumi stored, and there is currently no JES for events with NPV=0 - for these events these values are "approximated" according to recommendations from the JES group).
  • There are new expert options (provided by Mikko)
• p18-br-42 June 4, 2007:
  • ReComputeMET:
    • propagate tau energy scale correction to MET;
    • correct JES corrections in MET:
      • jet smearing is now correctly propagated to MET, but it looks like MET is now oversmeared in Z->ee+jets and QCD events.
      • the procedure to propagate the JES correction to MET is not really correct. We are still following a RunI prescription which is to propagate only the Response part to MET, but the formula implemented in jetcorr is not valid anymore with the RunIIA MPF method (which is different from the runI MPF method). The effect of showering and offset is however a small effect. We could therefore want to propagate the entire JES correction to MET (not yet implemented).
  • MetSigAlg:
    • Actually, the line const TMBMet* met = event.getMet(); is obsolete . It hasn't been used anywere in the code, so I removed it now in order to avoid futher confusion. The MET that I use is obtained later in the code: if (!event.get("caf_util::METx", metx)) throw runtime_error("METsigAlg [" + name() + "] ERROR. MET Selector must be run before.") ; if (!event.get("caf_util::METy", mety)) throw runtime_error("METsigAlg [" + name() + "] ERROR. MET Selector must be run before.") ; So, like it says one has to run METSelector before running METSigAlg. And METSelector selects the corrected MET by default
    • As for the factor of 2 for muons it seems you are right once again. Changed.
    • The cause of the negative values at the output has been found and fixed.
  • Tau: small modifications to TauSelector to use the new isMC() function in event and added protection against potential crash if TMBTrack branch was not dereferenced prior to using a TRef.
  • MuonSelector: corrected dR max cut behavior
• p18-br-41, p21-br-21 June 4, 2007: ApplyJES: Major update. Several bugs fixed, more security checks and new debug information. Events without valid inst. lumi discared by default. Summary information now displayed. Bugs: hasMU flag wasn't set correctly, several bugs when ApplyDataErrorsToMC was set to true: MC jets with muons had DATA JES applied. JESMU errors weren't calculated. metC errors were negative. Should now all be fixed.
• p18-br-40, p21-br-20 May 24, 2007: Adapted interfaces to new JES
• p18-br-39 May 14, 2007: Added support to store only leading selected object in MuonSelector.
• p18-br-38 May 13, 2007: METsigAlg corrected. Adopt isMC() methode from cafe::Event
• p18-br-37 May 4, 2007: ElectronSelector: correct bug in reading fiducial cut files. Update printout.
• p18-br-36 Apr 12, 2007: MEtsigAlg: MET likelihood cut is only applied if the cut value is greater than 0.0.
• p18-br-35 Apr 11, 2007: MuonSelector: add veto on muons from the specified branch. ReComputeMET: add jet-muon dR cut to muon correction for JESMU.
• p18-br-34 Apr 4, 2007: Removed ApplyMuonSmear - this processor is now in caf_mc_util. Added informative messages to ReComputeMET.
• p18-br-33 Apr 3, 2007: Correct electron correction in ReComputeMET.
• p18-br-32 Apr 3, 2007: Correct coarse hadronic bug in ReComputeMET. Removed ActAs warning messages. Instead the JES/smearing to be used is displayed in the constructor. Tau update
• p18-br-31 Apr 2, 2007: Added caf_util::HistoHolder -- a histogram container class.
• p18-br-30 : TauSelector: Include a control parameter on hether to apply MC or DATA energy scale corrections. Add phi crack selections.
• p18-br-29 : PhotonSelector: bug corrections and update to the certified selection.
• p18-br-28 : Add cuts for EC photons to photon selector.
• p18-br-27 : Add AddSelection processor. Add example for using JetControler.
• p18-br-26 : PhotonSelector: add possibility to cut on number of photons; add Stat printout. TauSelector modified.
• p18-br-25 : Add L2EMSelector
• p18-br-23 : Add MET significance processor.
• p18-br-22 : MuonSelector: correct on selection _noSMThit.
• p18-br-21 : JetSelector: correct pTJES cut. See this mail.
• p18-br-19 : Add TauEnergyCorr. TauSelector modified. MuonSelector modified (add option to choose closest beam spot run).
• p18-br-18 : Photon selector updated
• p18-br-17 : Bug fix for FindDuplicateEvents
• p18-br-16 : Add photon selector. Add Track selector and track smearing processors. The track smearing use the muon smearing parameters for now.
• p18-br-14 : TauSelector: add possibility to cut on NNe
• p18-br-13 : VertexSelector: add option to use other vertices than the first one in the selections. IsoTrackSelection: bug fix.
• p18-br-12 :
  • FindDuplicateEvents processor updated to have extra features when it works with MC eventes. The new version retain its old algorithm, and users who are not using the new features can always use it as before. The new features are:
    1. Capability to generate a table containing information about duplicate events and on which input file they occur.
    2. Capability to read the table from feature (1) and to remove it appropriately.
  • Selector for isolated tracks was added.
• p18-br-11 : Minor changes on the TauJetMatching. Now requires user to supply explicitly the name of the tau branch and throw an exception if not. Sort option is set by default to false. Improvement in the debug message and processor's finish() routine.
• p18-br-10 JetSelector: add dR cut between jet and any TLorentVector. ElectronSelector: add cut on maximal electron pt
• p18-br-09 JetSelector: an additional minimum cut on the JES corrected jet pT. (ignores ActAs) [default: 0.0] this allows you to apply a JSSR-like threshold in data on JESMU corrected jets.
• p18-br-08 add TauJetMathing processor to remove jets matched to taus.
• p18-br-07 Correct reading the EM ID cut file from the release directory
• p18-br-06 Add reading the EM ID cut file from the release directory
• p18-br-05 Put back TFile include in FindDuplicated event processor
• p18-br-04 TauSelector: major changes, mainly to recalculate NN outputs in the selector.
• p18-br-03 New muon smearing.
  - The default for the case that no parameter set is given in a way that all smearing parameters are set to zero and therefore privously smearing are undone the pt will be that of ptcorr.
  - "p17_06" is now one additional parameter set in musmear.txt A set "no_smearing" also added which sets the parameters to zero and therefore undoes the smearing. The description in ApplyMuonSmear.hpp updated.
• p18-br-02 Jet selector updated.
  .UseStandardID -> UseJetIDAndEMRemoval: turn-off default and give error message with instructions on how to properly use or not to use EM removal (the old default, configurable naming and header documentation combined were a trap for people who did not need EM removal and didn't notice the default setting).
  .UseStandardQualityCuts -> UseStandardJetID; default true and only applied if UseJetIDAndEMRemoval false, as before.
  .UseStandardEMRemoval -> UseEMRemoval; default false and only applied if both above false, as before Using old parameters will result in an error message with instructions and code exiting.
• v03-00-04 : BTagJetSelector updated. It now uses the TMBBTag object to get all its information instead of using the quick access methods on the jets. Some crashes were found due to this methods pointing nowhere in some cases.
• v03-00-03 : Remove "ntrack > 2" cut for the verteces in ApplyJES. Imrove printout in MuonSelector.
• v03-00-02 : Fix ReComputeMET and StatBranch (Variable class changed in cafe)
• v03-00-01 : Add Trigger object selectors : L1CalTowerSelector, L3BaseSelector, L3TrackSelector
• v02-04-18 : New muon smearing. This is the original mail from Raimund
  Subject: new smearing for p19_09 Monte Carlo
  Date: Sun, 02 Jul 2006 18:57:59 +0200
  From: Raimund Stroehmer
  Hi,
  I have determine the smearing parameters for p19_09 MC and implemented them into ApplyMuonSmear. The new implementation can be found on ... (see caf_util v02-04-18, tmb_tree p18-br-71)
  The files affected are:
  caf_util/src/ApplyMuonSmear.cpp
  caf_util/ApplyMuonSmear.hpp
  caf_util/config/musmear.txt
  tmb_tree/src/TMBMuon.cpp
  tmb_tree/tmb_tree/TMBMuon.hpp
  An example config file is /work/gandalf-clued0/raimund/cafe_p18_07/my_mc.config
  In TMBMuon.cpp the basic smearing routine has been extended to handle the separation of muons with cft_eta>1.6.
  The pt dependence of the smearing is now implemented as A+B1/pt+B2*pt in the current parameter sets either B1 or B2 is set to zero to have the two parametrization A+B1/pt and A+B2*pt.
  ApplyMuonSmear can now read in the smearing parameters from a config file (musmear.txt)
  For: muon_smear.conf_name: p17_06
  the existing p17_06 parameters (using physics eta > 1.5) are used I have implanted smearing parameters for both parametrization
  p17_09 (A+B1/pt)
  p17_09_pt (A+B2*pt)
  both with systematic variations of the smearing at pt=40 GeV
  p17_09_1s_plus,p17_09_1s_minus (p17_09_pt_1s_plus,p17_09_pt_1s_minus) and of pt pt dependence of the smearing (with fixed smearing at pt=40 GeV):
  p17_09_1s_pt_plus,p17_09_1s_pt_minus
  (p17_09_pt_1s_pt_plus,p17_09_pt_1s_pt_minus)
  I have tested the implementation but if some body else could have an additional look that would be good.
  Raimund