November 2012
To Be done this Month :
1. Approval slides. DONE
2. Modifying the logic of reconstructing the
PFChargedHadrons in the PFlow algorithm. (working on login03 : ~/TauPOGPatTuples/CMSSW_5_2_3/src)
- Plots Erec/Egen for generated 3-prong 0pi0 decays for allowTracksToBecomeNeutralHadrons = true and allowTracksToBecomeNeutralHadrons = false (will tell us if the allowTracksToBecomeNeutralHadrons flag works as expected and we improve the tau energy reconstruction for 3-prongs.)
- make 1d plots of (recMEtPx - genMEtPx) and (recMetPy - genMEtPy) (to check that our change to the treatment of high Pt tracks by the particle flow algorithm does not degrade the MET resolution.
Since the
WtoTauNu events have genuine MET, the figure of merit is the difference recMet - genMEt.)
- make 2d correlation plots of (recMEtPx - genMEtPx | allowTracksToBecomeNeutralHadrons = true) vs (recMEtPx - genMEtPx | allowTracksToBecomeNeutralHadrons = false)
and (recMEtPy - genMEtPy | allowTracksToBecomeNeutralHadrons = true) vs (recMEtPy - genMEtPy | allowTracksToBecomeNeutralHadrons = false) (to check that our change to the treatment of high Pt tracks by the particle flow algorithm does not degrade the MET resolution.
Since the
WtoTauNu events have genuine MET, the figure of merit is the difference recMet - genMEt.)
3. Installing New Tau ID with _3Hits and checking tau id efficiency (Using DYJets) and fake rate (Using
WtauNu & Data & QCD)
- Code installed at : /cmshome/khurana/NewHZZAnalyzer/Moriond13/CMSSW_5_3_3_patch3/src/Analyzer/TauEfficiencyAnalyzer
4. Installing Electron Isolation, Muon Isolation, Electron MVA ID, with NEW TAU ID (also Electron scale factor, Muon correction ). And prepare 2011, 2012ABC tuples AGAIN.
5. Remove all merged files from /lustre/home/khurana. KEEP THOSE FILES IN MULTIPLE COPIES WHICH ARE USED FOR HCP ITERATION AND DOCUMENT HOW TO RUN THE CODE ON THOSE MERGED FILES.
6. Put all HCP related code in CVS user code area. DONE
7. Check signal - background efficiency for 2012C data set after making following changes -
- Switch to new TAU ID
- Use Looser object definitions
- Use BTag
- Use sumpt of lepton cut (subleading Z)
- Use Z_pt cut
- Relax Z1 isolation
- Include dr(Z1,Z2)
1. Check How RM are installed for HO.
2. Check how to take data for HO,
SiPM.
3. Check how to take various runs and what is the purpose of these different runs.
4. Check how to analyse data and write CODE for analysis in classes format.
5. Check the detector geometry for New HCAL.
6. Get Access to 904 and fix EDH password.
7. Do Biocell training and get access to USC.
8. Once access to USC go with Pooja and see where muTCA is located and check things there.
9. Check what expertise are needed for HG upgrade from my side ( mainly during April when cables will be changed)
10. What about FPGA work ...?
1st to 5th November 2012 :
Working for approval slides, Tau ID efficiency, final limit etc.
6th November 2012 : Tuesday
Analysis Approved. And working on cosmetic changes of the plots.
Plots after modifications are at :
http://cmsdoc.cern.ch/~khurana/PublicPlotsForHCP.html
7th November 2012 : Wednesday
Not much work done today. Added new Tau ID in the efficiency measuring code and submitted jobs for it.
Will get new Tau ID efficiency by 8th November.
Most of the day went in sleeping.
8th, 9th, 10th November.
Preparing plot with proper style. Problem was not in script but it was because of root version.
It seems 4l is preparing plot in 5.29 version of
ROOT (because plot looks fine from Sascha laptop).
But i don't have access to that version. and that version is not available anymore.
11th November
Working on Tau ID code for met resolution before and after changing the PF algo. For MET Info in the treeanalyzer add one new class METInfo.cc
12th November
Tau ID Meeting
HO Burn In meeting
When reading lines from a text file using python, the end-line character often needs to be truncated before processing the text, as in the following example:
f = open("myFile.txt", "r")
for line in f:
line = line[:-1]
# do something with line
13th November
Meeting with Beni & Bobby L Bobby explained his code in some details, I have to learn how to run this code and how exactly is can be written in the Beni framework.
Visit to B 904 : Boring and sleepy
Visit to B28 : Quite intresting and Vladimir explained many things about
SiPM and RM boxes of HO.
Basic terminology :
SiPM card has 18
SiPM each having 3000 pixels, (for HB 30000pixels i guess),
SiPM card is connected to controller card : mainly to regulate the bias voltage from power suppy,
Control card has one bias voltage controller to check the bias voltage (check with Sanjeev)
if there is scratch on
SiPM it doesn't effect its working.
Light from led is coming via fibres (18) and these are fed to light guides (mixtures so that light hit the whole surface of
SiPM not only corresponding to
area of optical fiblre, This light falls on
SiPMs and light is converted into analig current. This current is sent to QIE cards to digitize this. There are three QIE cards
in each RM box and each card has 6 QIE chips (in total 3*6=18 QIE chips) each of these QIE chips are being used to integrate charge from one
SiPm and digitize it and then send to the
HTR card of HO.
14th, 15th November
Worked on New Tau ID efficiency : Plot sent to Savin and looks okay.
Now have to look at fake rate with new Tau ID from data with lltt cuts.
When moved to New Tau ID code started to crash because of new
PhysicsTools/PatAlgos version and
using older version fixed this but have to send on hyupernews why this is happening (crash is because of electron isolations)
Submitted jobs on data for new Tau ID.
16th, 17th November
Checked status of data jobs and almost finished (>90%)
Submitted jobs on MC 2012 for new Tau ID
Fixed code for Customised PF Work. Wrote Analyzer for the same. Pattuples were not available for all events so submitted jobs to get it. Pat tuplisation done now.
Wrote plotting script for Efficiency, METPx_diff, RECMETPt, GENMETPt ,
MEtPx_diff_vs_MEtPy_diff , RECMETPt_vs_GENMETPt.
Efficiency Plots were empty and after some debugging found that there was a check that event should be a Z->tautau at generated level (this is important if I run code on DYJets) but I
was running on
WTauNu Sample so this cut will never pass. Removed this cut and fixed the cut.
18th, 19th November 2012 :
Almost full data MC analyzed now (2012ABCD) with New Tau ID. Will be done by tomorrow , and then I will start producing the numbers.
For Tau ID Found one bug : cloningsnippet works and change the input tag when initial change is in the sequence otherwise not.
Trying to fix this using the patutils.massSearchReplaceAnyInputTag
OR i guess similar can be done using the following snippet
for modulename in process.p.moduleNames():
module = getattr(process, modulename)
if hasattr(module, 'src'):
if module.src == cms.InputTag("gsfElectrons"):
module.src = "correctedGsfElectrons"
if (hasattr(module, 'electronSource') and
module.electronSource == cms.InputTag("gsfElectrons")):
module.electronSource = "correctedGsfElectrons"