March 2012

Main Tasks for March 2012

1. Prepare Slides for Tau ID Meeitng - DONE on 10th March

1a. including new plot of PFGamma.

2. Prepare Slides for Jet faking lepton - send around. DONE on 11th March

3. Include mmet for electron & eemt for muon fake rate in the Analyzer and ploting script.

4. Understand tag and probe tree filling and structure.

5. Run Shape Analysis

6. Do TnP for electrons.

Table of contents:

9th March 2012

Some of the codes are working fine. But still need to merge codes together such that they work in one go.

Installing ROOT again to run ROOFit. Existing root version can be configured for ROOTFIT using -

./configure --enable -ROOFIT
make
Learning python lists and related stuff - Important links -

What is list and Yield

how to use csv - for reading column instead of row from text file

Note : import csv if you want to use csv. To use Yield later on for further processing typecast it to a list.

Working in jet faking lepton. Plot ready for data/MC comparison for loose and tight isolation working point.

Started working on Tag and Probe foe Electron : Following instructions from Anil.

Tag And Probe

As a first attempt things working fine on MC. Have to edit cfg to calculate efficieny of -

Electron Reconstruction

Electron Tight (CiCTight)

Electron Isolation

Double Electron Trigers

10th March 2012

To do -

1. Prepare Slides for Tau ID Meeitng including new plot of PFGamma.

2. Prepare Slides for Jet faking lepton - send around. DONE

3. Include mmet for electron & eemt for muon fake rate in the Analyzer and ploting script.

11th March 2012

Tag and Probe Follow up from Code -

1. Apply good electron selection (loose) on gsf and store them. User defined - as per understanding

2. On O/P of step one apply Required Cut - for which efficiency is to be estimated (may be wrong). User defined - as per understanding

3. Using O/P of one and two make combination (di-lepton) first one is Tag and second one is Probe. Done using CandViewShallowCloneCombiner from TnP official package.

4. Do some analysis using these dilepton pair and store required information. Done using TagProbeFitTreeProducer from TnP official package.

Above four steps are processed by following lines
process.p = cms.Path(
    process.gsfElectron
    +process.gsfWP90Cut
    +process.tagProbe
    +process.tpTree
)

Made Slides for TES and set for comments.

Script to get Number of events for each sample is almost done.

12th March 2012

Slides presented for TES.

Script to get Number of events for each sample is DONE.

Script to get datacard for each mass point is done if method used is not shape analysis.

Updated numbers for signal and bkg after PU reweighting, cross cleaning and with errors and propagated to the limits.

Tag and Probe cfgs are understood better and trying to understand how tree is filled.

TagProbeFitTreeAnalyzer

13th March 2012

Tried to get certificate from DU and CSIR.

Added pt, eta of leading and subleading e and mu in code, Propagated to plotting script.

checked hits cut for muon and found bug, not accessing pixel hit properly. Now fixed.

Started on P-X - Worked for 2 hrs..

14th March 2012

Tried to get certificate from DU & CSIR

15th March 2012

Travelling to SINP Full day in travelling.

16th March 2012

Half day in travelling.

Trying to get PFGamma plot when gendecay is 1prong0pi0 and reco as 1prong1pi0. DONE.

Adding Pt, eta plot for e1,e2,mu1,mu2.

Trying to add data-MC/MC plots on same canvas with proper axis and other stuff.

17th March 2012

Try to get Strip Pt plot when gendecay is 1prong0pi0 and reco as 1prong1pi0.

Fit ratio histogram with straight line

How to fit

Completed fitting of pt, eta, iso for leading and subleading e,mu but seeing ~20% less data.

Asked Ian to check his plots.

One reason can be loss of storage area at bari which may cause loss of some of data files.

Running pat-tupliZer on data again, Will take 3-4 days.

In the mean time looking at TnP for CiCtight id only and fake rate application realted script and Analyzer code.

Looked at Bhawna's Z peak code, didn't worked fine and couldn't find any bug even. frown

18th March 2012

SUNDAY. Reached ~12.00hrs.

No work after lunch.:)

19th March 2012

For TES - Tried to fix code to get information of strips but can't be fixed in CMSSW_4_2_5. So Set up code in 442 and submitted pat-tuplizer jobs in 442. Before submitting I tried for 100 events it is working fine in 442 atleast. Now waiting for jobs to finish. smile

TES code is not working because to get new fuctions DataFormats/PatCandidates has to be updates and hence have to make new pat tuples. But I was trying on old pat tuples.

For HZZ - Check BR of DYJetsToLL sample, it may be a reason for seeing less fakerate for jet->lepton from MC. TODO

Start writting script to multiply all TnP weights to MC. TODO.

See TnP package how TnP uses cic in official package.

Change PU weights in HZZ analyzer, they seems to be a little different. Pile up Official page

Pile Up Reweight Utilities

For Tau POG meeting - Slides sent to Sascha and these are fine.

20th March 2012

Attended Tuesday Meeting & prepared slides for Sunanda Sir.

Not much work did today. Little bit about TES.

21th March 2012

1. mmtt mass plot when a. No tau iso - No ISo means - iso not passed or it can be anything you just don't check b. Iso on first applied. c. Iso on second applied 2. Do 1 a, 1b. 1c for loose and medium tau isolation for eett, mmtt. 1 & 2 are now added in analyzer

3. Add MC to jet->tau rate plot. NOT DONE YET.

4. mmmt & eemt mass when (for jet->mu)

  1. No Iso on 3rd lepton b. When Loose Iso on third lepton c. When Tight Iso On 3rd letpon NOT DONE YET
5. Repeat 4a,4b,4c for eeet, mmet (for jet faking electron). NOT DONE YET

6. See how to estimate fit uncertainity as done by Abdollah. Abdollah will explain it on Friday. I tried something and that is working don't know how good it is :P

7. Run Analyzer on new Pat-Tuples and send plots to Christian. Make SignalPiZeroCandidate Pt in for the gen. 1Prong0pi0 taus which are reconstructed as 1Prong1pi0 ? It would actually be very useful if you could make the PFGamma Pt and the strip Pt plots also for the combination gendecay = oneProng1pi0 and reco decay mode = oneProng1pi0 --> Then we can see more clearly what would be a good Pt cut (the aim is to make the Pt cut high enough so that we reject low Pt PFGammas/strips in gen. oneProng0pi0, but at the same time low enough that we keep the PFGammas/strips gen. oneProng1pi0) Sent plots to Christian and recieved his reply. Going to reconstruct tau after strip pt > 2.5 GeV and see the decay mode performance plots. Submitted jobs after changing the HPSPFTau_cfi.py. Dificult to submit the jobs :P

22th March 2012

Attended Tracker DQM Shift Tutorial.

Checked crab status but jobs are failing with 8020 again and again. frown

Started writing the script to multiply TNP efficiency scale factors. Half done.

Attended 2l2tau internal meeting.

Attended EGAMMA meeting. Nice tutorial on PF Isolation by Florian.

23th March 2012

Read about TRACKER DQM SHIFT

TNP Meeting at 5.00 GVA IMPORTANT.

24th March 2012

Read about TRACKER DQM SHIFT

Plotted decay mode performance plots after decay mode finding only and we could find that -

25th March 2012

Read about TRACKER DQM SHIFT

And after lunch no work.

26th March 2012

TRACKER DQM SHIFT

Shift and shift tutorial took whole day.

27th March 2012

TRACKER DQM SHIFT

Energy Recovery Algorithm instructions by Christian -

------------------------------------------------------------------------------------------------
in today's Tau id. meeting there was the question for a recipe to run
the tau energy recovery algorithm in CMSSW_4_2_x .

I think that all you need to do is to check-out 2 additional files on
top of the latest RecoTauTag tags for CMSSW_4_2_x.
   V01-02-11      RecoTauTag/Configuration
   V01-02-16      RecoTauTag/RecoTau
   V01-02-01      RecoTauTag/TauTagTools
   cvs up -r 1.13 RecoTauTag/RecoTau/python/
RecoTauCombinatoricProducer_cfi.py
   cvs co -r 1.2 RecoTauTag/RecoTau/plugins/
RecoTauEnergyRecoveryPlugin.cc

If you then rerun the PFTau sequence, the reco::PFTau->p4() will
return the four-vector set by the energy recovery algorithm
(and will be returned by pat::Tau->p4() also)

What you should see after using the tau energy recovery algorithm is
that the tail of low tauJetPt/genVisTauPt gets reduced significantly
and the energy response becomes more Gaussian
------------------------------------------------------------------------------------------------------

28th March 2012

TRACKER DQM SHIFT

Important Links

DQM GUI

Run Registry

White Board

Tutorial Slide

29th March 2012

TRACKER DQM SHIFT

Higgs Production at LHC via Glu Glu fusion

30th March 2012

31th March 2012

Edit | Attach | Watch | Print version | History: r21 | r19 < r18 < r17 < r16 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r17 - 2012-03-27 - RamanKhurana
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

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.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback