PhotonsInData < Main

---+ Studies in data involving photons, for CIT undergrads
<h1>
   * <p>Yousi Ma</p>
   * <p><span style="font-size: 42px; ">Susan Dittmer</span></p>
   * <p><span style="font-size: 42px; ">Ben Faber</span></p>
   * <p><span style="line-height: 29px; font-size: 42px;">Elliot Schneider</span></p>
<span style="font-size: 42px; ">
---+ How to make event displays:

</span></h1>

%NAVY%
<h1><span style="font-size: 42px; ">

f<span style="font-size: 41.76px;">irst filter a RECO file to have only those events you have selected.</span>

e.g. use your analyzer, to pring run:lumi:event triples for the high mass events

Use the Edmtool for this:

https://twiki.cern.ch/twiki/bin/view/CMS/WorkBookPickEvents

if you don't have a cern account:

edmPickEvents.py is part of the =PhysicsTools/Utilities= package with tag =V08-01-01= and works for <span>CMSSW 382</span> or higher. 

How to Run =edmPickEvents.py=

The script =edmPickEvents.py= can be run interactively or via a CRAB job. If you run interactively, before you run it, make sure you have the dataset that you want to pick events from, available locally. If you have few events to pick, run it interactively. In case you do not have the dataset available locally or you have a lot of events to pick, submit a CRAB job.

Run =edmPickEvents.py= Interactively

To run the script interactively do the following, for example

edmPickEvents.py "/EG/Run2010A-PromptReco-v4/RECO" 140160:98:82878561 140160,101,85437704

these instructions are for cern machines, reformat the output so that the files starte with file:/mnt/tier2/store instead of /store

running this command returns a RECO root file, and then all you have to do is:

cmsShow (your root file) 



To setup a CMSSW release look <a target="_top" href="https://twiki.cern.ch/twiki/bin/view/CMS/WorkBookSetComputerNode#RunTime">HERE</a> and then checkout the tag as follows:

cvs co -r V08-01-01 PhysicsTools/Utilities

_Please remember_ to do =scram b= and =rehash= after checking out the package.
</span></h1>

<span style="font-family: 'Courier New', courier, monaco, monospace; color: #7a4707;"><span style="white-space: pre;"><br /></span></span>

---+ This week's data

6807nb are finished and available in /raid1/yma/MPAfiles
---+ Last week's data:

all 4150nb-1 is finished. It had to be run without the bit 40 or 41 that most data up till now has been processed with. please check your selection yields/pb-1 and check these triggers.

/raid1/yma/MPAfiles/MPA_Run2010B_oct22_4140nb.root &lt;- please analyze me

4150nb-1 coming in Run2010B:

now that the files are getting large and precious, they will be placed in hadoop (redundant backup).

/mnt/hadoop/store/user/yma/MPAdata

/mnt/hadoop/store/user/yma/MPAMC
---+ Current open studies/issues

   * Z yields in EE/EE and EB/EE pairs
   * Z mass shift due to Endcap transparency change
   * high mass diphoton pairs (add R9 cut)
   * high mass gammaJet pairs (event displays)
   * effect of preselection (no TT bits, adding a sample with a beam scraping filter)
   * diphoton rate in data, week to week
---+ previous data:

up to oct15th JSON

7875nb-1 in Run2010B.

3068(3049)nb-1 in Run2010A(old partial) file.

---+ previous results
<table style="width: 964px; border: 0px solid #0bf3be;" cellspacing="0" cellpadding="0" border="0"> <col width="215"></col> <col width="187" span="4"></col> <tbody> <tr height="39"> <td style="border: 1px solid #12ec95;">

Runperiod
</td> <td style="border: 1px solid #12ec95;">

Name
</td> <td style="border: 1px solid #12ec95;">

Selection
</td> <td style="border: 1px solid #12ec95;">

Lumi
</td> <td style="border: 1px solid #12ec95;">

Notes
</td> </tr> <tr height="39"> <td style="border: 1px solid #12ec95;">

run2010B &ndash; oct15
</td> <td style="border: 1px solid #12ec95;">

mpa_out_DATA_noPix_Run2010B_7875nb.root
</td> <td style="border: 1px solid #12ec95;">

noTT, noPix, noHLT, pT20, + old tight EgID
</td> <td style="border: 1px solid #12ec95;">

7875nb-1
</td> <td style="border: 1px solid #12ec95;">

Limited use?
</td> </tr> <tr> <td style="border: 1px solid #12ec95;"> </td> <td style="border: 1px solid #12ec95;"> </td> <td style="border: 1px solid #12ec95;"> </td> <td style="border: 1px solid #12ec95;"> </td> <td style="border: 1px solid #12ec95;"> </td> </tr> <tr height="39"> <td style="border: 1px solid #12ec95;">

run2010B &ndash; oct15
</td> <td style="border: 1px solid #12ec95;">

mpa_out_DATA_Run2010B_ExoID_noTT_7875nb.root
</td> <td style="border: 1px solid #12ec95;">

noTT, noPix, noHLT, pt20 + ExoID[1], EgId[2]
</td> <td style="border: 1px solid #12ec95;">

7875nb-1
</td> <td style="border: 1px solid #12ec95;">

EgID for 2nd by accident
</td> </tr> <tr> <td style="border: 1px solid #12ec95;">run2010B - oct15</td> <td style="border: 1px solid #12ec95;">mpa_out_DATA_Run2010B_ExoID_noTT_7875nb.root</td> <td style="border: 1px solid #12ec95;">noTT, noPix, noHLT, pt20 +ExoID</td> <td style="border: 1px solid #12ec95;">7875nb-1</td> <td style="border: 1px solid #12ec95;">FIX of previous line</td> </tr> <tr height="39"> <td style="border: 1px solid #12ec95;">

run2010A
</td> <td style="border: 1px solid #12ec95;">

mpa_out_DATA_Run2010A_ExoID_noTT_3049nb.root
</td> <td style="border: 1px solid #12ec95;">

noTT, noPix, noHLT, pt20 + ExoIDs
</td> <td style="border: 1px solid #12ec95;">

3049nb-1
</td> <td style="border: 1px solid #12ec95;">

Correct ExoID for both
</td> </tr> <tr height="39"> <td style="border: 1px solid #12ec95;">

run2010A
</td> <td style="border: 1px solid #12ec95;">

mpa_out_DATA_Zsel_noTT_Run2010A_3049nb.root
</td> <td style="border: 1px solid #12ec95;">

yesPix, noTT, noHLT, pt20, EgIDs
</td> <td style="border: 1px solid #12ec95;">

3049nb-1
</td> <td style="border: 1px solid #12ec95;">

Tight isolations, Z selection
</td> </tr> <tr height="39"> <td style="border: 1px solid #12ec95;">

run2010B
</td> <td style="border: 1px solid #12ec95;">

mpa_out_DATA_Zsel_noTT_Run2010B_Oct15_7875nb.root
</td> <td style="border: 1px solid #12ec95;">

yesPix, noTT, noHLT, pt20, EgIDs
</td> <td style="border: 1px solid #12ec95;">

7875nb-1
</td> <td style="border: 1px solid #12ec95;">

Tight isolations, Z selection
</td> </tr> <tr> <td style="border: 1px solid #12ec95;">run2010B - Oct22</td> <td style="border: 1px solid #12ec95;">mpa_out_DATA_Zsel_noTT_Run2010B_oct22_4150nb.root</td> <td style="border: 1px solid #12ec95;">yesPix, noTT, noHLT, pt20, ExoID, no '40 OR 41'</td> <td style="border: 1px solid #12ec95;">4150nb-1</td> <td style="border: 1px solid #12ec95;">Zselection</td> </tr> <tr> <td style="border: 1px solid #12ec95;">run2010B - Oct22</td> <td style="border: 1px solid #12ec95;">mpa_out_DATA_Run2010B_ExoID_noTT_Oct22_4150nb.root</td> <td style="border: 1px solid #12ec95;">

noPix, noTT, noHLT, pt20, ExoID

no '40 OR 41'
</td> <td style="border: 1px solid #12ec95;">4150nb-1</td> <td style="border: 1px solid #12ec95;">ExoID</td> </tr> </tbody> </table>

all of Run2010A is available in septReRECO:

/raid1/yma/MPAfiles/MPA_Run2010A_reRECO_3068nb.root

please add your 2.85pb-1 plots here.
---+ Datasets

all paths relative to /raid1/yma/MPAfiles
<table cellspacing="0" cellpadding="0" border="0"> <col width="298"></col><col width="128"></col><col width="133"></col><col width="64"></col> <tbody> <tr height="21"> <td style="border: 1px solid #c33ba4;">

Run Range
</td> <td style="border: 1px solid #c33ba4;">dataset name</td> <td style="border: 1px solid #c33ba4;">

location
</td> <td style="border: 1px solid #c33ba4;">

Lumi
</td> <td style="border: 1px solid #c33ba4;">

comments
</td> <td style="border: 1px solid #c33ba4;"> </td> </tr> <tr height="35"> <td style="border: 1px solid #c33ba4;">

136035-144114

</td> <td style="border: 1px solid #c33ba4;">/EG/Run2010A-Sep17ReRECO</td> <td style="border: 1px solid #c33ba4;">

MPA_Run2010A_reRECO_partial_3049nb.root
</td> <td style="border: 1px solid #c33ba4;">

3049nb-1
</td> <td style="border: 1px solid #c33ba4;">

missing ~20nb-1 from 2/407 jobs
</td> <td style="border: 1px solid #c33ba4;"> </td> </tr> <tr> <td style="border: 1px solid #c33ba4;">146428 - 147454</td> <td style="border: 1px solid #c33ba4;">/Photon/Run2010B-promptreco</td> <td style="border: 1px solid #c33ba4;">MPA_Photon_run2010B_Oct15_7875nb.root</td> <td style="border: 1px solid #c33ba4;">7875nb-1</td> <td style="border: 1px solid #c33ba4;">typo causing missing HLT_double17</td> <td style="border: 1px solid #c33ba4;"> </td> </tr> <tr> <td style="border: 1px solid #c33ba4;">147757 - 148058</td> <td style="border: 1px solid #c33ba4;">/Photon/Run2010B-promptreco</td> <td style="border: 1px solid #c33ba4;">MPA_Run2010B_oct22_4150nb.root</td> <td style="border: 1px solid #c33ba4;">4150nb-1</td> <td style="border: 1px solid #c33ba4;">no bit '40 OR41' requirement vs. before</td> <td style="border: 1px solid #c33ba4;"> </td> </tr> <tr> <td style="border: 1px solid #c33ba4;">148822 - 148864</td> <td style="border: 1px solid #c33ba4;">/PHoton/Run2010B_promptreco</td> <td style="border: 1px solid #c33ba4;">MPA_Run2010B_Oct29_6807nb.root</td> <td style="border: 1px solid #c33ba4;">6807nb-1</td> <td style="border: 1px solid #c33ba4;">also recovered 147754, 147755, 148031</td> <td style="border: 1px solid #c33ba4;"></td> </tr> </tbody> </table>
---+ Monte Carlo

locations relative to /raid1/yma/MPAfiles/HiggsMCSpring2010
<table style="border: 0px solid #400fef;" cellspacing="0" cellpadding="0" border="0"> <col width="298"></col><col width="128"></col><col width="133"></col><col width="64"></col> <tbody> <tr height="21"> <td style="border: 1px solid #f10d32;">

Sample
</td> <td style="border: 1px solid #f10d32;">

Ngen
</td> <td style="border: 1px solid #f10d32;">

Cross section
</td> <td style="border: 1px solid #f10d32;">

filter
</td> <td style="border: 1px solid #f10d32;">MPA_ntuple name</td> <td style="border: 1px solid #f10d32;">Ntuple events</td> </tr> <tr height="35"> <td style="border: 1px solid #f10d32;">

QCD_EMEnriched_Pt20to30
</td> <td style="border: 1px solid #f10d32;">

30 M
</td> <td style="border: 1px solid #f10d32;">

0.2355 mb
</td> <td style="border: 1px solid #f10d32;">

0.007
</td> <td style="border: 1px solid #f10d32;"></td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

QCD_EMEnriched_Pt30to80
</td> <td style="border: 1px solid #f10d32;">

40 M
</td> <td style="border: 1px solid #f10d32;">

0.0593 mb
</td> <td style="border: 1px solid #f10d32;">

0.059
</td> <td style="border: 1px solid #f10d32;"></td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

QCD_EMEnriched_Pt80to170
</td> <td style="border: 1px solid #f10d32;">

5 M
</td> <td style="border: 1px solid #f10d32;">

0.906e-3 mb
</td> <td style="border: 1px solid #f10d32;">

0.148
</td> <td style="border: 1px solid #f10d32;"></td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

<a rel="nofollow" href="https://twiki.hep.caltech.edu/twiki/bin/edit/Main/PhotonJet?topicparent=Main.HiggsToGammaGamma;nowysiwyg=0" title="PhotonJet (this topic does not yet exist; you can create it)">PhotonJet</a> _Pt20to30
</td> <td style="border: 1px solid #f10d32;">

100000
</td> <td style="border: 1px solid #f10d32;">

5.718e+04
</td> <td style="border: 1px solid #f10d32;">

1
</td> <td style="border: 1px solid #f10d32;">MPA_PhotonJet_pt20to30.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

<a rel="nofollow" href="https://twiki.hep.caltech.edu/twiki/bin/edit/Main/PhotonJet?topicparent=Main.HiggsToGammaGamma;nowysiwyg=0" title="PhotonJet (this topic does not yet exist; you can create it)">PhotonJet</a> _Pt30to50
</td> <td style="border: 1px solid #f10d32;">

100000
</td> <td style="border: 1px solid #f10d32;">

1.652e+04
</td> <td style="border: 1px solid #f10d32;">

1
</td> <td style="border: 1px solid #f10d32;">MPA_PhotonJet_pt30to50.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

<a rel="nofollow" href="https://twiki.hep.caltech.edu/twiki/bin/edit/Main/PhotonJet?topicparent=Main.HiggsToGammaGamma;nowysiwyg=0" title="PhotonJet (this topic does not yet exist; you can create it)">PhotonJet</a> _Pt50to80
</td> <td style="border: 1px solid #f10d32;">

100000
</td> <td style="border: 1px solid #f10d32;">

2.723e+03
</td> <td style="border: 1px solid #f10d32;">

1
</td> <td style="border: 1px solid #f10d32;">MPA_PhotonJet_pt50to80.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

<a rel="nofollow" href="https://twiki.hep.caltech.edu/twiki/bin/edit/Main/PhotonJet?topicparent=Main.HiggsToGammaGamma;nowysiwyg=0" title="PhotonJet (this topic does not yet exist; you can create it)">PhotonJet</a> _Pt80to120
</td> <td style="border: 1px solid #f10d32;">

100000
</td> <td style="border: 1px solid #f10d32;">

4.462e+02
</td> <td style="border: 1px solid #f10d32;">

1
</td> <td style="border: 1px solid #f10d32;">MPA_PhotonJet_pt80to120.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

<a rel="nofollow" href="https://twiki.hep.caltech.edu/twiki/bin/edit/Main/PhotonJet?topicparent=Main.HiggsToGammaGamma;nowysiwyg=0" title="PhotonJet (this topic does not yet exist; you can create it)">PhotonJet</a> _Pt120to170
</td> <td style="border: 1px solid #f10d32;">

100000
</td> <td style="border: 1px solid #f10d32;">

8.443e+01
</td> <td style="border: 1px solid #f10d32;">

1
</td> <td style="border: 1px solid #f10d32;">MPA_PhotonJet_pt120to170.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

<a rel="nofollow" href="https://twiki.hep.caltech.edu/twiki/bin/edit/Main/DiPhotonBorn?topicparent=Main.HiggsToGammaGamma;nowysiwyg=0" title="DiPhotonBorn (this topic does not yet exist; you can create it)">DiPhotonBorn</a> _Pt25to250
</td> <td style="border: 1px solid #f10d32;">

500k
</td> <td style="border: 1px solid #f10d32;">

22.37
</td> <td style="border: 1px solid #f10d32;">

1
</td> <td style="border: 1px solid #f10d32;">MPA_box_pt25to250.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

<a rel="nofollow" href="https://twiki.hep.caltech.edu/twiki/bin/edit/Main/DiPhotonBox?topicparent=Main.HiggsToGammaGamma;nowysiwyg=0" title="DiPhotonBox (this topic does not yet exist; you can create it)">DiPhotonBox</a> _Pt25to250
</td> <td style="border: 1px solid #f10d32;">

500k
</td> <td style="border: 1px solid #f10d32;">

12.37
</td> <td style="border: 1px solid #f10d32;">

1
</td> <td style="border: 1px solid #f10d32;">MPA_Born_pt25to250.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

Hgg_gluonfusion_M90,110,120,150
</td> <td style="border: 1px solid #f10d32;">

100k
</td> <td style="border: 1px solid #f10d32;">

</td> <td style="border: 1px solid #f10d32;">

</td> <td style="border: 1px solid #f10d32;">MPA_Hgg_gluonfusion_MXXX.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

Hgg_wbf_M90,110,120,150
</td> <td style="border: 1px solid #f10d32;">

100k
</td> <td style="border: 1px solid #f10d32;">

</td> <td style="border: 1px solid #f10d32;">

</td> <td style="border: 1px solid #f10d32;">MPA_Hgg_wbf_MXXX.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> <tr height="20"> <td style="border: 1px solid #f10d32;">

Hgg_wh_zh_tth_M90,110,120,150
</td> <td style="border: 1px solid #f10d32;">

100k
</td> <td style="border: 1px solid #f10d32;">

</td> <td style="border: 1px solid #f10d32;">

</td> <td style="border: 1px solid #f10d32;">MPA_Hgg_wh_zh_tth_MXXX.root</td> <td style="border: 1px solid #f10d32;"></td> </tr> </tbody> </table>

samples from QCD group/older processing:

/raid1/yma/MPAfiles/MCData
---+ code/CVS

cms.hep.caltech.edu/viewvc - to browse

http://ultralight.caltech.edu/gaewiki/CVS for instructions - ben has used it successfully on friday.

I've made an area under apps repository, /apps/MPAcode for us to share. please post your preselection and selection codes here.
This topic: Main > PhotonsInData
Topic revision: r10 - 2010-11-03 - yma
Copyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback