#! /usr/bin/env python import sys, math, csv, json, pprint from array import array import ROOT as r import lumiTools as l ### Get Data/MC files ### bcmData = [] bcmData.extend(l.read_csv_file('data_db/BunchByBunch/Fill2663_NoHeader1.csv')) ### Set time frame ### timeStart = 0 timeStop = int(len(bcmData)*1.1) duration = int(timeStop - timeStart) ### Declare histograms ### outFile = r.TFile('histograms/fedHistograms.root', 'RECREATE') outFile.mkdir('OR', 'OR') #outFile.mkdir('XOR-', 'XOR-') #outFile.mkdir('XOR+', 'XOR+') #outFile.mkdir('AND', 'AND') r.TH1.SetDefaultSumw2(r.kTRUE) r.TProfile.SetDefaultSumw2(r.kTRUE) h1_Intensities = r.TH1D('h1_Intensities', 'I_{1}#times I_{2};I_{1}#times I_{2};N_{bx}', 50, 0.8, 1.8) h1_IntensityCorr = r.TH1D('h1_IntensityCorr', 'I_{1}#times I_{2}/<I_{1}#times I_{2}>;BCID;I_{1}#times I_{2}/<I_{1}#times I_{2}>', 3564, 0.5, 3564.5) p1_IntensVsTime = r.TProfile('p1_IntensVsTime', 'I_{1}#times I_{2};BCID;I_{1}#times I_{2}', 3564, 0.5, 3564) outFile.cd('OR') h1_ProbabilityOr = r.TH1D('h1_ProbabilityOr', 'Probability of BCM1F OR;BCID;p(hit)', 3564, 0.5, 3564.5)
#! /usr/bin/env python import sys, math, csv, json, pprint from array import array import ROOT as r import lumiTools as l #fillList = [2629, 2630, 2634, 2535, 2644, 2645, 2646, 2648, 2649, 2651, 2653, 2657, 2658] fillList = [2520] bcmData = [] #bcmData.extend(l.read_csv_file('data_db/Current_bcm.csv')) bcmData.extend(l.read_csv_file('data_db/fill2520bcm.csv')) #pixelData = [] #pixelData.extend(l.read_json_file('data/pixel_lumi_raw_2012.json')) ### Declare histograms ### timeStart = 0 timeStop = int(len(bcmData)*1.1) duration = int(timeStop - timeStart) start = r.TDatime('2012-05-15 12:28:44').Convert() stop = r.TDatime(bcmData[int(len(bcmData))-1]['DT']).Convert() #stop = r.TDatime(bcmData[int(len(bcmData))-1]['DT_UNTIL']).Convert() outFile = r.TFile('histograms/lumiHistograms.root', 'RECREATE') outFile.mkdir('OR', 'OR') outFile.mkdir('XOR-', 'XOR-') outFile.mkdir('XOR+', 'XOR+')
scanResults = [] #-----------# # Fill 2855 # #-----------# if 2855 in fills: ### Scan configuration and corrections ### scans = ['1', '2'] #, '3'] ### Input data files ### bcmData = [ l.read_csv_file('data/July/VDMScanXY1.csv'), l.read_csv_file('data/July/VDMScanXY2.csv'), l.read_csv_file('data/July/VDMScanXY3.csv') ] bunchData = [ l.read_csv_file('data/July/FBCT_2855_Timber/TIMBER_DATA_Fill_2855_BCTFR_A6R4_B1.csv', bunchHeader), l.read_csv_file('data/July/FBCT_2855_Timber/TIMBER_DATA_Fill_2855_BCTFR_A6R4_B2.csv', bunchHeader) ] pltData = r.TFile('histograms/vdmHistograms_2011-07-19_2855.root', 'OPEN') ### Initialize calibrator for fill 2855 ### calibrator = c.VDMCalibrator(2855, bcmData, bunchData, corrections, {'PLT':1342721696 - 1798, 'intensity':1342685132}, algos, dataFormat = 'brm') scanResults.extend(calibrator.multi_scan(scans, bcids, pltData, False, savePath)[0])