Conversation
…d neutrinos to particle dictionary
| @@ -0,0 +1,227 @@ | |||
| #! /usr/bin/env python3 | |||
There was a problem hiding this comment.
we don't want to commit this file. You should be able to git rm it
| @@ -0,0 +1,8 @@ | |||
| #!/bin/bash -x | |||
|
|
|||
| git clone -b cleanUp https://github.com/jstupak/aQGC.git | |||
There was a problem hiding this comment.
once this is merged, folks can use master branch again. So you should remove -b cleanUp
| particle={ | ||
| 11:"e", | ||
| 13:"mu", | ||
| 21:"gamma", #using this for all jets |
There was a problem hiding this comment.
I think I introduced this bug, but 21 is a gluon, 22 is the photon
|
|
||
| #------------------------------------------------------------------------------------------------- | ||
| leptons=electrons+muons | ||
|
|
There was a problem hiding this comment.
This is where we want to add some code. First, lets define a new kind of object, a "W jet." W jets should have mass in the 50-100 GeV range (just looking at your leading jet mass plot).
Then, you should check and make sure you found exactly two beam remnant muons and 2 W jets. If not, you can just skip to the next event with a continue statement. This is the "object selection" we want to apply.
| DEBUG=True | ||
| DEBUG=False | ||
|
|
||
| jetType='KTjet' |
There was a problem hiding this comment.
lets switch to the large cone (12?) that mayuri showed here (with N=2)
| #------------------------------------------------------------------------------------------------- | ||
| leptons=electrons+muons | ||
|
|
||
| Zs={} |
There was a problem hiding this comment.
Now that we decided to focus on the case with 2 BR remnants and no additional leptons, you can comment out this chunk of code (through line 332). Just make the Z candidate from your two BR muons
| truthPhotons =selector(event.Particle,'abs(x.Status)==1 and abs(x.PID)==22') | ||
| truthNeutrinos=selector(event.Particle,'x.Status==1 and (abs(x.PID)==12 or abs(x.PID)==14 or abs(x.PID)==16)') | ||
|
|
||
| beamRemnantMuons =selector(truthMuons,'isBeamRemnant(x)') |
There was a problem hiding this comment.
we should probably rename this truthBeamRemnantMuons
1 participant
No description provided.