Back to home page

sPhenix code displayed by LXR

 
 

    


Warning, /analysis/EMCal-commissioning/emcalHotTowerFinder/readme.txt is written in an unsupported language. File is not indexed.

0001 These are the instructions for running the code to generate lists of noisy towers 
0002 and generate the corresponding CDBTTrees
0003 
0004 Build code:
0005 
0006 mkdir build
0007 cd build
0008 $PWD/../src/autogen.sh --prefix=${Whatever you've named your install environment variable}
0009 
0010 head to macros directory
0011 
0012 Directory structure for condor
0013      mkdir out out_hadd cdbMaps err err_hadd submission submission_hadd mkdir fileLists log log_hadd output
0014 
0015 cd fileLists
0016 
0017 Create file lists for data of interest
0018 
0019 Create file list of filelists
0020 
0021 ls -1 *.list > fileListofLists.list
0022 
0023 cd ../
0024 
0025 makeSingleFileList.csh
0026 
0027 runHotTowermacroMultistage.csh fileLists/fileListofLists.list sigma
0028 
0029 This will submit a condor job for each segment of each run to be processed individually.
0030 Then it will hadd the output files after pass 0
0031 Then it will submit a job per run, which is actually making the bad towers, don't delete the output files from the previous step!
0032 Then it writes the cdbttrees to the cdbMaps directory. These are immediately ready for upload. 
0033 
0034 Where sigma is the desired cutoff. The code currently employs it symmetrically but it can be asymmetric as well.
0035 
0036 Toweridmacro
0037 - Takes in a run from the infile (supplied as argument)
0038 - Iterates through each event (supplied as argument)
0039 - For each tower that exceeds the adc cut(supplied as argument) in an event, the tower registers a "hit" and its frequency is incremented
0040 - A histogram is filled with the total number of hits for each tower throughout the total number of events
0041 - The mode and Standard Deviation of the histogram are calculated. A cutoff frequency is determined by adding a set number of standard deviations (sigmas: supplied as argument) to the mode
0042 - All towers that exceed (record more hits than) the cutoff frequency are listed by channel number in a TTree 
0043 - The histogram containing the hits for each tower, and the TTree with hot channels are written to an outfile (supplied as argument)
0044 
0045 Fun4All_EMCalBadTowerID.C
0046 nEvents - Number of events to process
0047 inFile - input file
0048 cdbname - name of cdbttree output file
0049 adccut_sg - minimum adc value for blocks containing sg-type fibers
0050 adccut_k - minimum adc value for blocks containing k-type fibers
0051 sigmas_lo - number of standard deviations below  mean before a tower is cold
0052 sigmas_hi - number of standard deviations above mean before tower is noisy
0053 SG_f - fiducial cut on tower hit frequency for sg-type fibers
0054 Kur_f - fiducial cut on tower hit frequency for k-type fibers
0055 region_f - fiducial cut on region hit frequency
0056 pass - iteration
0057