CaloMonTransforms Namespace Reference

Functions
def	cloneHistFloat (toClone)
def	divideHistCaloMon (inputs, doPercentage=False)
def	divideByOccupancy (inputs, titleToReplace="", replaceTitWith="")
def	divideByAcceptedEvts (inputs)
def	simpleDivideByAcceptedEvts (inputs, titleToReplace="", replaceTitWith="")

Function Documentation

cloneHistFloat()

def CaloMonTransforms.cloneHistFloat

(

toClone

)

If we want to clone the bins from an integer-type histogram, but then wish to fill it with floats (for efficiencies) 

Definition at line 16 of file CaloMonTransforms.py.

def cloneHistFloat(toClone):
    """ If we want to clone the bins from an integer-type histogram, but then wish to fill it with floats (for efficiencies) """
    histype = str(type(toClone))
    newhist = None
    if "TH1I" in histype:
        newhist = R.TH1F()
        toClone.Copy(newhist)
    elif "TH2I" in histype:
        newhist = R.TH2F()
        toClone.Copy(newhist)
    else:
        newhist = toClone.Clone()
    return newhist
 
 

divideByAcceptedEvts()

def CaloMonTransforms.divideByAcceptedEvts

(

inputs

)

This function divides the input histogram by the number of events that passed a given selection

Definition at line 104 of file CaloMonTransforms.py.

def divideByAcceptedEvts(inputs):
    """ This function divides the input histogram by the number of events that passed a given selection"""
    assert len(inputs) == 1  
    assert len(inputs[0][1]) == 2
 
    cl = inputs[0][1][0].Clone()
 
    label = inputs[0][0]['thr']
 
    hnorm = inputs[0][1][1].Clone()
    n0 = hnorm.GetNbinsX()
    fbin = hnorm.GetXaxis().FindBin(label)
 
    if fbin>=n0 or fbin<=0: #threshold not found
        return [cl]
        
    thenorm = hnorm.GetBinContent(fbin)
    
    if thenorm==0:
        return [cl]
    cl.Scale(100.0/thenorm)
    
    return [cl]
 
 

divideByOccupancy()

def CaloMonTransforms.divideByOccupancy	(		inputs,
			titleToReplace = "",
			replaceTitWith = ""
	)

This function returns the ratio of two ROOT histograms 

Definition at line 73 of file CaloMonTransforms.py.

def divideByOccupancy(inputs,titleToReplace="",replaceTitWith=""):
    """ This function returns the ratio of two ROOT histograms """
    assert len(inputs) == 1  
    assert len(inputs[0][1]) == 2
 
    cl = cloneHistFloat(inputs[0][1][0])
    cl.Reset()
 
    if titleToReplace!="":
        tit = cl.GetTitle()
        tit = tit.replace(titleToReplace,replaceTitWith)
        cl.SetTitle(tit)
 
    nCells = cl.GetNcells()
    assert(nCells==inputs[0][1][1].GetNcells())
 
    cppyy.gbl.binbybinDiv(nCells,cl,inputs[0][1][0],inputs[0][1][1])
 #   for i in range(nCells):
  #      t = inputs[0][1][0].GetBinContent(i);
   #     o = inputs[0][1][1].GetBinContent(i);
    #    if o>0:
     #       cl.SetBinContent(i,t/o);
      #      cl.SetBinError(i,TMath.Sqrt(1./o));
       #     pass
        #pass
 
    cl.SetEntries(inputs[0][1][0].GetEntries())
 
    return [cl]
 
 

divideHistCaloMon()

def CaloMonTransforms.divideHistCaloMon	(		inputs,
			doPercentage = False
	)

This function returns the ratio of two ROOT histograms 

Definition at line 31 of file CaloMonTransforms.py.

def divideHistCaloMon(inputs,doPercentage=False):
    """ This function returns the ratio of two ROOT histograms """
    assert len(inputs) == 1  
    assert len(inputs[0][1]) == 3
 
    cl = cloneHistFloat(inputs[0][1][0])
    cl.Divide(inputs[0][1][1])
 
    label = inputs[0][0]['thr']
 
    hnorm = inputs[0][1][2].Clone()
    n0 = hnorm.GetNbinsX()
    fbin = hnorm.GetXaxis().FindBin(label)
 
    if fbin>=n0 or fbin<=0:
        return [cl]
        
    thenorm = hnorm.GetBinContent(fbin)
    cl.SetEntries(thenorm)
    
    if not doPercentage or thenorm==0:
        return [cl]
 
    cl.Scale(100.0/thenorm)
    
    return [cl]
 
 
cppyy.cppdef("""
void binbybinDiv(size_t nCells,TH2* fraction, const TH2* total, const TH2* occupancy) {
for (size_t i=0;i<nCells;++i) {
    const double t = total->GetBinContent(i);
    const double o = occupancy->GetBinContent(i);
    if (o>0) {
      fraction->SetBinContent(i,t/o);
      fraction->SetBinError(i,std::sqrt(1./o));
    }
  }
}
""")
 
 

simpleDivideByAcceptedEvts()

def CaloMonTransforms.simpleDivideByAcceptedEvts	(		inputs,
			titleToReplace = "",
			replaceTitWith = ""
	)

This function divides the input histogram by the number of accepted events

Definition at line 129 of file CaloMonTransforms.py.

def simpleDivideByAcceptedEvts(inputs,titleToReplace="",replaceTitWith=""):
    """ This function divides the input histogram by the number of accepted events"""
    assert len(inputs) == 1 #Expect only one match 
    assert len(inputs[0][1]) == 2 # pair of (regex-match,list-of-hists)
 
    cl = cloneHistFloat(inputs[0][1][0]) # First histogram of the list of histgrams in the pair 
    
    hnorm = inputs[0][1][1]    #Second histogram in the list of histograms is the number of events
    theNorm=hnorm.GetBinContent(1)
    if theNorm>1:
        cl.Scale(100.0/theNorm)
 
    if titleToReplace!="":
        tit = cl.GetTitle()
        tit = tit.replace(titleToReplace,replaceTitWith)
        cl.SetTitle(tit)
 
    return [cl]