hotSpotInTAG Namespace Reference

Functions
def	analyzeTree ()

Variables
	parser = argparse.ArgumentParser()
	type
	int
	dest
	default
	help
	action
	float
	args = parser.parse_args()
	run = args.runNumber
	lowerLumiBlock = args.lowerLB
	upperLumiBlock = args.upperLB
	stream = args.stream
	tag = args.tag
	amiTag = args.amiTag
	etaSpot = args.etaSpot
	phiSpot = args.phiSpot
	thresholdE = args.thresholdE
	deltaSpot = args.deltaSpot
	objectType = args.objectType
	minInLB = args.minInLB
	tagDirectory = args.fileDirectory
	tree = TChain("POOLCollectionTree")
	listOfFiles = pathExtract.returnEosTagPath(run,stream,amiTag,tag)
	entries = tree.GetEntries()
int	nLB = 2500
list	nbHitInHot = [0] * nLB
list	cumEnergInHot = [0] * nLB
list	nbNoiseBurstVeto = [0] * nLB
list	nbLArNoisyRO_Std = [0] * nLB
list	nbLArNoisyRO_SatTight = [0] * nLB
	h0map = TH1D("map","General map of %s with MET > %d MeV"%(objectType,thresholdE),64,-3.14,3.14)
	h0mapClean = TH1D("mapClean","General map of %s with MET > %d MeV - LArFlags != ERROR"%(objectType,thresholdE),64,-3.14,3.14)
	nb = tree.GetEntry( jentry )
list	nLB_offending = []
int	lowerLB = 2500
int	upperLB = 0
string	suffix = "NO LArFlags cut"
	c0 = TCanvas()
	h0Evol = TH1I("h0Evol","Nb of hits in a region of %.2f around (%.2f,%.2f) and Et/Pt > %d MeV - %s"%(deltaSpot,etaSpot,phiSpot,thresholdE,suffix),upperLB-lowerLB+1,lowerLB-0.5,upperLB+0.5)
	h0Evol_E = TH1I("h0Evol_E","Mean E/pt in a region of %.2f around (%.2f,%.2f) and Et/Pt > %d MeV - %s"%(deltaSpot,etaSpot,phiSpot,thresholdE,suffix),upperLB-lowerLB+1,lowerLB-0.5,upperLB+0.5)
	c0_2 = TCanvas()
	c0_3 = TCanvas()
list	canvas = []
int	iCurrent = len(canvas)-1
string	cutC = "1"

Function Documentation

analyzeTree()

def hotSpotInTAG.analyzeTree

(

)

Definition at line 33 of file hotSpotInTAG.py.

def analyzeTree():
  global nbHitInHot
  global cumEnergInHot
  global nbLArNoisyRO_Std
  global nbLArNoisyRO_SatTight
 
  if (tree.LArFlags & 8):nbNoiseBurstVeto[tree.LumiBlockN] = nbNoiseBurstVeto[tree.LumiBlockN] + 1
  # Fill the list of coordinates depending on the object type
  if ("TopoCluster" in objectType):
    coordEnergy = [[tree.TopoClusterEt1,tree.TopoClusterEta1,tree.TopoClusterPhi1]]
  if ("Jet" in objectType):
    coordEnergy = [[tree.JetPt1,tree.JetEta1,tree.JetPhi1],[tree.JetPt2,tree.JetEta2,tree.JetPhi2],[tree.JetPt3,tree.JetEta3,tree.JetPhi3],[tree.JetPt4,tree.JetEta4,tree.JetPhi4],[tree.JetPt5,tree.JetEta5,tree.JetPhi5],[tree.JetPt6,tree.JetEta6,tree.JetPhi6]]
  if ("LoosePhoton" in objectType):
    coordEnergy = [[tree.LoosePhotonPt1,tree.LoosePhotonEta1,tree.LoosePhotonPhi1],[tree.LoosePhotonPt2,tree.LoosePhotonEta2,tree.LoosePhotonPhi2]]
  if ("LooseElectron" in objectType):
    coordEnergy = [[tree.LooseElectronPt1,tree.LooseElectronEta1,tree.LooseElectronPhi1],[tree.LooseElectronPt2,tree.LooseElectronEta2,tree.LooseElectronPhi2],[tree.LooseElectronPt3,tree.LooseElectronEta3,tree.LooseElectronPhi3],[tree.LooseElectronPt4,tree.LooseElectronEta4,tree.LooseElectronPhi4]]
  if ("Tau" in objectType):
    coordEnergy = [[tree.TauJetPt1,tree.TauJetEta1,tree.TauJetPhi1],[tree.TauJetPt2,tree.TauJetEta2,tree.TauJetPhi2]]
  if ("MET" in objectType):
    coordEnergy = [[tree.CellMissingET,0,tree.CellMissingETPhi]]
  # Loop on the list of coordinates to analyze them
  for iCE in coordEnergy:
    if (iCE[0] > thresholdE):
      if ("MET" in objectType): # Fill general map
        h0map.Fill(iCE[2])
        if not (tree.LArFlags & 8):
          h0mapClean.Fill(iCE[2])
      else:
        h0map.Fill(iCE[1],iCE[2])
        if not (tree.LArFlags & 8):
          h0mapClean.Fill(iCE[1],iCE[2])
      if not ((not args.larcleaning) and tree.LArFlags & 8): # Explicit cut LArEventInfo != ERROR to remove noise bursts is larcleaning is false
        if (fabs(iCE[1]-etaSpot)<deltaSpot and fabs(iCE[2]-phiSpot)<deltaSpot): # A candidate has been found in suspicious region - Explict cut LArEventInfo != ERROR to remove noise bursts
#        if (fabs(iCE[1]-etaSpot)<deltaSpot and fabs(iCE[2]-phiSpot)<deltaSpot and not (tree.LArFlags & 8)): # A candidate has been found in suspicious region - Explict cut LArEventInfo != ERROR to remove noise bursts
          nbHitInHot[tree.LumiBlockN]=nbHitInHot[tree.LumiBlockN] + 1
          cumEnergInHot[tree.LumiBlockN]=cumEnergInHot[tree.LumiBlockN]+iCE[0]
          if (tree.LArFlags & 1): nbLArNoisyRO_Std[tree.LumiBlockN] = nbLArNoisyRO_Std[tree.LumiBlockN] + 1
          if (tree.LArFlags & 4): nbLArNoisyRO_SatTight[tree.LumiBlockN] = nbLArNoisyRO_SatTight[tree.LumiBlockN] + 1
 
  return
  
 
# Main===========================================================================================================

Variable Documentation

action

hotSpotInTAG.action

Definition at line 77 of file hotSpotInTAG.py.

amiTag

hotSpotInTAG.amiTag = args.amiTag

Definition at line 102 of file hotSpotInTAG.py.

args

hotSpotInTAG.args = parser.parse_args()

Definition at line 93 of file hotSpotInTAG.py.

c0

hotSpotInTAG.c0 = TCanvas()

Definition at line 192 of file hotSpotInTAG.py.

c0_2

hotSpotInTAG.c0_2 = TCanvas()

Definition at line 211 of file hotSpotInTAG.py.

c0_3

hotSpotInTAG.c0_3 = TCanvas()

Definition at line 217 of file hotSpotInTAG.py.

canvas

list hotSpotInTAG.canvas = []

Definition at line 225 of file hotSpotInTAG.py.

cumEnergInHot

list hotSpotInTAG.cumEnergInHot = [0] * nLB

Definition at line 148 of file hotSpotInTAG.py.

cutC

string hotSpotInTAG.cutC = "1"

Definition at line 231 of file hotSpotInTAG.py.

default

hotSpotInTAG.default

Definition at line 77 of file hotSpotInTAG.py.

deltaSpot

hotSpotInTAG.deltaSpot = args.deltaSpot

Definition at line 106 of file hotSpotInTAG.py.

dest

hotSpotInTAG.dest

Definition at line 77 of file hotSpotInTAG.py.

entries

hotSpotInTAG.entries = tree.GetEntries()

Definition at line 144 of file hotSpotInTAG.py.

etaSpot

int hotSpotInTAG.etaSpot = args.etaSpot

Definition at line 103 of file hotSpotInTAG.py.

float

hotSpotInTAG.float

Definition at line 83 of file hotSpotInTAG.py.

h0Evol

hotSpotInTAG.h0Evol = TH1I("h0Evol","Nb of hits in a region of %.2f around (%.2f,%.2f) and Et/Pt > %d MeV - %s"%(deltaSpot,etaSpot,phiSpot,thresholdE,suffix),upperLB-lowerLB+1,lowerLB-0.5,upperLB+0.5)

Definition at line 193 of file hotSpotInTAG.py.

h0Evol_E

hotSpotInTAG.h0Evol_E = TH1I("h0Evol_E","Mean E/pt in a region of %.2f around (%.2f,%.2f) and Et/Pt > %d MeV - %s"%(deltaSpot,etaSpot,phiSpot,thresholdE,suffix),upperLB-lowerLB+1,lowerLB-0.5,upperLB+0.5)

Definition at line 195 of file hotSpotInTAG.py.

h0map

hotSpotInTAG.h0map = TH1D("map","General map of %s with MET > %d MeV"%(objectType,thresholdE),64,-3.14,3.14)

Definition at line 154 of file hotSpotInTAG.py.

h0mapClean

hotSpotInTAG.h0mapClean = TH1D("mapClean","General map of %s with MET > %d MeV - LArFlags != ERROR"%(objectType,thresholdE),64,-3.14,3.14)

Definition at line 155 of file hotSpotInTAG.py.

help

hotSpotInTAG.help

Definition at line 77 of file hotSpotInTAG.py.

iCurrent

int hotSpotInTAG.iCurrent = len(canvas)-1

Definition at line 228 of file hotSpotInTAG.py.

int

hotSpotInTAG.int

Definition at line 77 of file hotSpotInTAG.py.

listOfFiles

hotSpotInTAG.listOfFiles = pathExtract.returnEosTagPath(run,stream,amiTag,tag)

Definition at line 126 of file hotSpotInTAG.py.

lowerLB

hotSpotInTAG.lowerLB = 2500

Definition at line 175 of file hotSpotInTAG.py.

lowerLumiBlock

hotSpotInTAG.lowerLumiBlock = args.lowerLB

Definition at line 98 of file hotSpotInTAG.py.

minInLB

hotSpotInTAG.minInLB = args.minInLB

Definition at line 108 of file hotSpotInTAG.py.

nb

hotSpotInTAG.nb = tree.GetEntry( jentry )

Definition at line 164 of file hotSpotInTAG.py.

nbHitInHot

list hotSpotInTAG.nbHitInHot = [0] * nLB

Definition at line 147 of file hotSpotInTAG.py.

nbLArNoisyRO_SatTight

list hotSpotInTAG.nbLArNoisyRO_SatTight = [0] * nLB

Definition at line 151 of file hotSpotInTAG.py.

nbLArNoisyRO_Std

list hotSpotInTAG.nbLArNoisyRO_Std = [0] * nLB

Definition at line 150 of file hotSpotInTAG.py.

nbNoiseBurstVeto

list hotSpotInTAG.nbNoiseBurstVeto = [0] * nLB

Definition at line 149 of file hotSpotInTAG.py.

nLB

int hotSpotInTAG.nLB = 2500

Definition at line 146 of file hotSpotInTAG.py.

nLB_offending

list hotSpotInTAG.nLB_offending = []

Definition at line 174 of file hotSpotInTAG.py.

objectType

hotSpotInTAG.objectType = args.objectType

Definition at line 107 of file hotSpotInTAG.py.

parser

hotSpotInTAG.parser = argparse.ArgumentParser()

Definition at line 76 of file hotSpotInTAG.py.

phiSpot

hotSpotInTAG.phiSpot = args.phiSpot

Definition at line 104 of file hotSpotInTAG.py.

run

hotSpotInTAG.run = args.runNumber

Definition at line 97 of file hotSpotInTAG.py.

stream

hotSpotInTAG.stream = args.stream

Definition at line 100 of file hotSpotInTAG.py.

suffix

string hotSpotInTAG.suffix = "NO LArFlags cut"

Definition at line 186 of file hotSpotInTAG.py.

tag

hotSpotInTAG.tag = args.tag

Definition at line 101 of file hotSpotInTAG.py.

tagDirectory

hotSpotInTAG.tagDirectory = args.fileDirectory

Definition at line 109 of file hotSpotInTAG.py.

thresholdE

hotSpotInTAG.thresholdE = args.thresholdE

Definition at line 105 of file hotSpotInTAG.py.

tree

hotSpotInTAG.tree = TChain("POOLCollectionTree")

Definition at line 124 of file hotSpotInTAG.py.

type

hotSpotInTAG.type

Definition at line 77 of file hotSpotInTAG.py.

upperLB

hotSpotInTAG.upperLB = 0

Definition at line 176 of file hotSpotInTAG.py.

upperLumiBlock

hotSpotInTAG.upperLumiBlock = args.upperLB

Definition at line 99 of file hotSpotInTAG.py.