readCCLHist Namespace Reference

Functions
def	readInput (filename, nCellsList=20, bcfile=None, cutoff=0.1)

Variables
	parser = argparse.ArgumentParser()
	type
	help
	nargs
	default
	None
	float
	int
	args
	leftover

Function Documentation

readInput()

def readCCLHist.readInput	(		filename,
			nCellsList = 20,
			bcfile = None,
			cutoff = 0.1
	)

Definition at line 10 of file readCCLHist.py.

def readInput(filename, nCellsList=20, bcfile=None, cutoff=0.1):
 
    #Set up LArOnlineID helper class in standalone mode (from xml file)
    from ROOT import IdDictParser
    parser=IdDictParser()
    #Get xml files:
    xmlpath=None
    for dd in os.getenv('XMLPATH').split(os.pathsep):
        d=dd+"/IdDictParser/ATLAS_IDS.xml"
        if os.access(d,os.R_OK):
            xmlpath=dd
            break
    if not xmlpath:
        print ("ERROR, unable to locate identifier dictionaries")
        sys.exit(-1)
               
 
    parser.register_external_entity("LArCalorimeter",xmlpath+"/IdDictParser/IdDictLArCalorimeter_DC3-05-Comm-01.xml")
    idd = parser.parse(xmlpath+"/IdDictParser/ATLAS_IDS.xml")
    from ROOT import LArOnlineID
    larID=LArOnlineID()
    stat=larID.initialize_from_dictionary(idd)
    if stat==1:
        print ("ERROR, failed to init LArOnlineID")
        sys.exit(-1)
 
 
    #Open ROOT HIST file
    f=TFile.Open(filename)
    if not f.IsOpen():
        print("ERROR, failed to open input file",filename)
        sys.exit(-1)
 
    rundir=None
    for d in f.GetListOfKeys():
        mg=re.match("run_([0-9]*)",d.GetName())
        rundir=mg.group(0)
        break
 
    print("Found run dir:",rundir)
    runnbr=rundir[4:]
 
    hist=f.Get(rundir+"/CaloMonitoring/LArClusterCellMon/Summary/cellhashPercent")
    
    freq=[]
    for idx in range(hist.GetNbinsX()):
        freq.append((idx,hist.GetBinContent(idx)))
    
    freq.sort(key=itemgetter(1),reverse=True)
  
    for f in freq[:nCellsList]:
        c=larID.channel_Id(f[0])
        print ("Channel %s constributes to clusters in %.3f %% of events"% (larID.channel_name(c),f[1]))
 
  
    if bcfile:
        bcfile.write("#Bad channel list for run "+runnbr+"\n")
        for (h,f) in freq:
            if f<cutoff: break
            c=larID.channel_Id(h)
            bcfile.write("%i %i %i %i %i 0 highNoiseHG\n"% (larID.barrel_ec(c), larID.pos_neg(c), larID.feedthrough(c),
                                                                  larID.slot(c), larID.channel(c))) 
 
        bcfile.close()
 
 
 

Variable Documentation

args

readCCLHist.args

Definition at line 86 of file readCCLHist.py.

default

readCCLHist.default

Definition at line 82 of file readCCLHist.py.

float

readCCLHist.float

Definition at line 83 of file readCCLHist.py.

help

readCCLHist.help

Definition at line 81 of file readCCLHist.py.

int

readCCLHist.int

Definition at line 84 of file readCCLHist.py.

leftover

readCCLHist.leftover

Definition at line 86 of file readCCLHist.py.

nargs

readCCLHist.nargs

Definition at line 82 of file readCCLHist.py.

None

readCCLHist.None

Definition at line 82 of file readCCLHist.py.

parser

readCCLHist.parser = argparse.ArgumentParser()

Definition at line 80 of file readCCLHist.py.

type

readCCLHist.type

Definition at line 81 of file readCCLHist.py.