Configuring_Hot_Towers

 
A code prepared to find out hot towers.

 This code running on the root files used by the macros which created L2btowW monitorig PDF at http://online.star.bnl.gov/L2algo/l2btowW/

Use the histo:  

W-BTOW-accepted: acc seed tower ET>5.00 GeV

Calculate mean (Xbar) and RMS : 

IF  
Xi -Xbar   >  Threshold  * RMS    ------>  HOT

***************************************************************************************
For ETOW, use the histo
     ETOW: Seed Tower RDO

Calculate mean (Xbar) and RMS : 

IF  
Xi -Xbar   >  Threshold  * RMS    ------>  HOT

Code plot a recap histogram hot towers Vs run number. And also out put text file create which include list of hot towers of each run. 

To test I use runs  from fill 17318.  

F17318, R14092087,430229,430305
 F17318, R14092090,430229,430305
 F17318, R14092091,430229,430305
 F17318, R14092092,430229,430305
 F17318, R14092093,430229,430305
 F17318, R14092097,430229,430305
 F17318, R14092098,430229,430305
 F17318, R14092099,430229,430305
 F17318, R14092100,430229,430305
 F17318, R14092101,430229,430305
 F17318, R14092104,430229,430305
 F17318, R14092105,430229,430305
 F17318, R14092106,430229,430305
 F17318, R14092107,430229,430305
 F17318, R14092108,430229,430305
 F17318, R14092109,430229,430305
 F17318, R14092110,430229,430305
 F17318, R14093001,430229,430305
 F17318, R14093005,430229,430305
 F17318, R14093006,430229,430305
 F17318, R14093007,430229,430305
 F17318, R14093008,430229,430305
 F17318, R14093009,430229,430305
 F17318, R14093010,430229,430305
 F17318, R14093013,430229,430305
 F17318, R14093014,430229,430305
 F17318, R14093015,430229,430305
 F17318, R14093016,430229,430305
 F17318, R14093017,430229,430305
 F17318, R14093018,430229,430305
 F17318, R14093019,430229,430305
 F17318, R14093020,430229,430305
 F17318, R14093021,430229,430305
***********************************************************

Now I need find out a good threshold value.

Threshold == 6

B-  3737     for R 14092097      L2plots ,  J plots      :                 for R 14092098        L2 plots ,  J plots

**********************************************************************************************************************************************************************
Threshold == 10

**********************************************************************************************************************************************************************

Threshold == 12

Seems to me actually hot in this set of runs are soft ID 3737 of Barrel and Soft ID 615 in Endcap.

my feeling is to use threshold = 12

**********************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************

BEMC Hot Towers of P14ia Production
 

Plots of Run vs hot Towers for each fill  :    
pdf threshold = 12                                 pdf threshold = 6

Spread sheet :      Hot Towers Soft ID of each run  : csv