Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

HCAL: miscellaneous updates for CMT: all depthes for CapIdErrors #44443

Merged
merged 2 commits into from
Mar 25, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 1 addition & 2 deletions DPGAnalysis/HcalTools/README
Original file line number Diff line number Diff line change
@@ -1,7 +1,6 @@
use
##setenv SCRAM_ARCH slc7_amd64_gcc900

setenv SCRAM_ARCH ec8_amd64_gcc11
setenv SCRAM_ARCH ec8_amd64_gcc12
cmsenv
voms-proxy-init -voms cms -rfc

475 changes: 326 additions & 149 deletions DPGAnalysis/HcalTools/macros/cmt/RemoteMonitoringGLOBAL.cc

Large diffs are not rendered by default.

Original file line number Diff line number Diff line change
Expand Up @@ -118,14 +118,14 @@ int main(int argc, char *argv[]) {
//CUTS: [test][subdetector] CapID(Test=1; ADC amplitude Am(Test= 2); Width for Wm(Test=3); Ratio cut for Rm(Test=4); TS mean for TNm(test=5); TS max for TXm(Test=6);
double MIN_M[7][5] = {{0., 0., 0., 0., 0.},
{0., 0., 0., 0., 0.},
{0, 0., 0., 0., 0.},
{0, 10., 100., 10., 100.},
{0, 1.0, 1.0, 0.2, 0.1},
{0, 0.10, 0.10, 0.18, 0.30},
{0, 0.8, 0.8, 0.8, 0.1},
{0, -0.5, -0.5, -0.5, -0.5}};
double MAX_M[7][5] = {{0., 0., 0., 0., 0.},
{0., 0., 0., 0., 0.},
{0, 900, 900, 9000, 3000},
{0, 300000, 300000, 5000, 100000},
{0, 3.9, 3.9, 4.4, 2.0},
{0, 0.95, 0.98, 0.96, 1.04},
{0, 8.0, 8.0, 8.0, 2.8},
Expand Down
5 changes: 5 additions & 0 deletions DPGAnalysis/HcalTools/plugins/BuildFile.xml
Original file line number Diff line number Diff line change
@@ -1,4 +1,5 @@
<use name="root"/>
<use name="boost"/>
<use name="DataFormats/HcalDetId"/>
<use name="DataFormats/HcalDigi"/>
<use name="DataFormats/HcalRecHit"/>
Expand All @@ -12,8 +13,12 @@
<use name="Geometry/CaloTopology"/>
<use name="Geometry/HcalTowerAlgo"/>
<use name="DataFormats/FEDRawData"/>
<use name="DataFormats/Common"/>
<use name="CalibFormats/HcalObjects"/>
<use name="CalibCalorimetry/HcalAlgos"/>
<use name="FWCore/ServiceRegistry"/>
<use name="FWCore/Framework"/>
<use name="FWCore/ParameterSet"/>
<use name="CommonTools/UtilAlgos"/>
<flags EDM_PLUGIN="1"/>

228 changes: 214 additions & 14 deletions DPGAnalysis/HcalTools/plugins/CMTRawAnalyzer.cc

Large diffs are not rendered by default.

26 changes: 2 additions & 24 deletions DPGAnalysis/HcalTools/scripts/cmt/a.py.beg
Original file line number Diff line number Diff line change
Expand Up @@ -8,50 +8,28 @@
# with command line options: RelVal -s ALCA:HcalCalMinBias --data --scenario=pp -n 100 --conditions auto:run3_data_prompt --datatier USER --customise=HLTrigger/Configuration/CustomConfigs.L1THLT --era Run3 --processName=USER --filein file:RelVal_RECO.root --fileout file:RelVal_ALCA.root --no_exec
import FWCore.ParameterSet.Config as cms

from FWCore.MessageService.MessageLogger_cfi import *

from Configuration.Eras.Era_Run3_cff import Run3

process = cms.Process('USER',Run3)

# import of standard configurations
process.load('Configuration.StandardSequences.Services_cff')
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.MagneticField_cff')
#process.load('Configuration.StandardSequences.AlCaRecoStreams_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')

#process = cms.Process("TEST", eras.Run2_2018)
#process = cms.Process("TEST", eras.Run3)
#process.load("Configuration.StandardSequences.GeometryDB_cff")
#process.load("CondCore.CondDB.CondDB_cfi")
#process.load("EventFilter.L1GlobalTriggerRawToDigi.l1GtUnpack_cfi")
#process.l1GtUnpack.DaqGtInputTag = 'source'
#----------------------------------------------------------------------------
#process.load('Configuration.StandardSequences.Services_cff')
#process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
#process.load('FWCore.MessageService.MessageLogger_cfi')
#process.load('Configuration.EventContent.EventContent_cff')
#process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
#process.load('Configuration.StandardSequences.Reconstruction_Data_cff')
#process.load('Configuration.StandardSequences.EndOfProcess_cff')
#process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
#process.load('Configuration.StandardSequences.EndOfProcess_cff')
#----------------------------------------------------------------------------
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(-1),
output = cms.optional.untracked.allowed(cms.int32,cms.PSet)
)
#process.maxEvents = cms.untracked.PSet(
## input = cms.untracked.int32(1000)
# input = cms.untracked.int32(-1)
#)
#----------------------------------------------------------------------------
#----------------------------------------------------------------------------
# Input source
process.source = cms.Source("PoolSource",
# skipBadFiles=cms.untracked.bool(True),
skipBadFiles=cms.untracked.bool(False),
fileNames = cms.untracked.vstring(
119 changes: 54 additions & 65 deletions DPGAnalysis/HcalTools/scripts/cmt/a.py.end
Original file line number Diff line number Diff line change
Expand Up @@ -2,9 +2,6 @@
secondaryFileNames = cms.untracked.vstring()
)

process.options = cms.untracked.PSet(
SkipEvent = cms.untracked.vstring('ProductNotFound')
)

# Production Info
process.configurationMetadata = cms.untracked.PSet(
Expand Down Expand Up @@ -91,7 +88,7 @@ process.Analyzer = cms.EDAnalyzer("CMTRawAnalyzer",
# -53 for BAD HBHEHF channels from study on shape Ratio
#Verbosity = cms.untracked.int32(-53),
ratioHBMin = cms.double(0.10),
ratioHBMax = cms.double(0.95),
ratioHBMax = cms.double(0.99),
ratioHEMin = cms.double(0.10),
ratioHEMax = cms.double(0.98),
ratioHFMin = cms.double(0.30),
Expand All @@ -100,7 +97,7 @@ process.Analyzer = cms.EDAnalyzer("CMTRawAnalyzer",
ratioHOMax = cms.double(0.96),
# -54 for BAD HBHEHF channels from study on RMS of shapes
#Verbosity = cms.untracked.int32(-54),
rmsHBMin = cms.double(1.0),
rmsHBMin = cms.double(0.6),
rmsHBMax = cms.double(3.9),
rmsHEMin = cms.double(1.0),
rmsHEMax = cms.double(3.9),
Expand All @@ -110,7 +107,7 @@ process.Analyzer = cms.EDAnalyzer("CMTRawAnalyzer",
rmsHOMax = cms.double(4.4),
# -55 for BAD HBHEHF channels from study on TSmean of shapes
#Verbosity = cms.untracked.int32(-55),
TSmeanHBMin = cms.double(1.8),
TSmeanHBMin = cms.double(0.8),
TSmeanHBMax = cms.double(8.0),
TSmeanHEMin = cms.double(1.8),
TSmeanHEMax = cms.double(8.0),
Expand All @@ -130,14 +127,14 @@ process.Analyzer = cms.EDAnalyzer("CMTRawAnalyzer",
TSpeakHOMax = cms.double(6.5),
# -56 for BAD HBHEHOHF channels from study on ADC Amplitude
#Verbosity = cms.untracked.int32(-56),
ADCAmplHBMin = cms.double(10000.),
ADCAmplHBMin = cms.double(100.),
ADCAmplHBMax = cms.double(300000.),
ADCAmplHEMin = cms.double(20000.),
ADCAmplHEMin = cms.double(10.),
ADCAmplHEMax = cms.double(300000.),
ADCAmplHFMin = cms.double(50.),
ADCAmplHFMax = cms.double(9000.),
ADCAmplHOMin = cms.double(50.),
ADCAmplHOMax = cms.double(9000.),
ADCAmplHFMin = cms.double(100.),
ADCAmplHFMax = cms.double(100000.),
ADCAmplHOMin = cms.double(10.),
ADCAmplHOMax = cms.double(5000.),
#
# to see channels w/ PedestalSigma < cut
#Verbosity = cms.untracked.int32(-57),
Expand Down Expand Up @@ -245,23 +242,23 @@ process.Analyzer = cms.EDAnalyzer("CMTRawAnalyzer",
# cuts on Nbadchannels to see LS dependences:
# Verbosity = cms.untracked.int32(-77),
# to select abnormal events,for which Nbcs > this limits
lsdep_cut1_peak_HBdepth1 = cms.int32(20),
lsdep_cut1_peak_HBdepth2 = cms.int32(7),
lsdep_cut1_peak_HEdepth1 = cms.int32(16),
lsdep_cut1_peak_HEdepth2 = cms.int32(13),
lsdep_cut1_peak_HEdepth3 = cms.int32(4),
lsdep_cut1_peak_HFdepth1 = cms.int32(10),
lsdep_cut1_peak_HFdepth2 = cms.int32(5),
lsdep_cut1_peak_HOdepth4 = cms.int32(45),
lsdep_cut1_peak_HBdepth1 = cms.int32(0),
lsdep_cut1_peak_HBdepth2 = cms.int32(0),
lsdep_cut1_peak_HEdepth1 = cms.int32(0),
lsdep_cut1_peak_HEdepth2 = cms.int32(0),
lsdep_cut1_peak_HEdepth3 = cms.int32(0),
lsdep_cut1_peak_HFdepth1 = cms.int32(0),
lsdep_cut1_peak_HFdepth2 = cms.int32(0),
lsdep_cut1_peak_HOdepth4 = cms.int32(0),
# to select events with Nbcs > this limits
lsdep_cut3_max_HBdepth1 = cms.int32(19),
lsdep_cut3_max_HBdepth2 = cms.int32(6),
lsdep_cut3_max_HEdepth1 = cms.int32(15),
lsdep_cut3_max_HEdepth2 = cms.int32(12),
lsdep_cut3_max_HEdepth3 = cms.int32(3),
lsdep_cut3_max_HFdepth1 = cms.int32(9),
lsdep_cut3_max_HFdepth2 = cms.int32(4),
lsdep_cut3_max_HOdepth4 = cms.int32(40),
lsdep_cut3_max_HBdepth1 = cms.int32(0),
lsdep_cut3_max_HBdepth2 = cms.int32(0),
lsdep_cut3_max_HEdepth1 = cms.int32(0),
lsdep_cut3_max_HEdepth2 = cms.int32(0),
lsdep_cut3_max_HEdepth3 = cms.int32(0),
lsdep_cut3_max_HFdepth1 = cms.int32(0),
lsdep_cut3_max_HFdepth2 = cms.int32(0),
lsdep_cut3_max_HOdepth4 = cms.int32(0),
#
#
#old was for runs:
Expand All @@ -270,23 +267,23 @@ process.Analyzer = cms.EDAnalyzer("CMTRawAnalyzer",
# nbadchannels3 = cms.int32(50),
#
#Verbosity = cms.untracked.int32(-79),
# cuts on Estimator1 to see LS dependences: default for normal pp - 2500 everywhere;
lsdep_estimator1_HBdepth1 = cms.double(1500.),
lsdep_estimator1_HBdepth2 = cms.double(1500.),
lsdep_estimator1_HBdepth3 = cms.double(1500.),
lsdep_estimator1_HBdepth4 = cms.double(1500.),
lsdep_estimator1_HEdepth1 = cms.double(1500.),
lsdep_estimator1_HEdepth2 = cms.double(1500.),
lsdep_estimator1_HEdepth3 = cms.double(1500.),
lsdep_estimator1_HEdepth4 = cms.double(1500.),
lsdep_estimator1_HEdepth5 = cms.double(1500.),
lsdep_estimator1_HEdepth6 = cms.double(1500.),
lsdep_estimator1_HEdepth7 = cms.double(1500.),
lsdep_estimator1_HFdepth1 = cms.double(140.),
lsdep_estimator1_HFdepth2 = cms.double(140.),
lsdep_estimator1_HFdepth3 = cms.double(140.),
lsdep_estimator1_HFdepth4 = cms.double(140.),
lsdep_estimator1_HOdepth4 = cms.double(180.),
# cuts on Estimator1 to see LS dependences
lsdep_estimator1_HBdepth1 = cms.double(1700.),
lsdep_estimator1_HBdepth2 = cms.double(2000.),
lsdep_estimator1_HBdepth3 = cms.double(2000.),
lsdep_estimator1_HBdepth4 = cms.double(1800.),
lsdep_estimator1_HEdepth1 = cms.double(10000.),
lsdep_estimator1_HEdepth2 = cms.double(4500.),
lsdep_estimator1_HEdepth3 = cms.double(3000.),
lsdep_estimator1_HEdepth4 = cms.double(3000.),
lsdep_estimator1_HEdepth5 = cms.double(3000.),
lsdep_estimator1_HEdepth6 = cms.double(3000.),
lsdep_estimator1_HEdepth7 = cms.double(6500.),
lsdep_estimator1_HFdepth1 = cms.double(300.),
lsdep_estimator1_HFdepth2 = cms.double(300.),
lsdep_estimator1_HFdepth3 = cms.double(300.),
lsdep_estimator1_HFdepth4 = cms.double(300.),
lsdep_estimator1_HOdepth4 = cms.double(150.),
# cuts on Estimator2 to see LS dependences: default 7.
lsdep_estimator2_HBdepth1 = cms.double(7.),
lsdep_estimator2_HBdepth2 = cms.double(7.),
Expand Down Expand Up @@ -472,31 +469,23 @@ process.load('Configuration.StandardSequences.RawToDigi_Data_cff')
# #InputLabel = cms.InputTag("rawDataCollector"),
#)
###########################################################################
#process.load("Calibration.HcalAlCaRecoProducers.ALCARECOHcalCalMinBias_cff")
###########################################################################
# Schedule definition
# pure CMT:
process.p = cms.Path(process.hcalDigis*process.Analyzer)
# for CMT use with GlobalPSM:
#process.p = cms.Path(process.hcalDigis*process.seqALCARECOHcalCalMinBiasDigiNoHLT*process.seqALCARECOHcalCalMinBias*process.Analyzer)
# see /afs/cern.ch/work/z/zhokin/public/CMSSW_10_4_0_patch1/src/Calibration/HcalAlCaRecoProducers/python/ALCARECOHcalCalMinBias_cff.py
############################################################################
# customisation of the process.
# Automatic addition of the customisation function from HLTrigger.Configuration.CustomConfigs
#from HLTrigger.Configuration.CustomConfigs import L1THLT
#call to customisation function L1THLT imported from HLTrigger.Configuration.CustomConfigs
#process = L1THLT(process)
# End of customisation functions
############################################################################
############################################################################
############################################################################
process.MessageLogger = cms.Service("MessageLogger",
categories = cms.untracked.vstring(''),
destinations = cms.untracked.vstring('cout'),
debugModules = cms.untracked.vstring('*'),
cout = cms.untracked.PSet(
threshold = cms.untracked.string('WARNING'),
WARNING = cms.untracked.PSet(limit = cms.untracked.int32(0))
)
)
MessageLogger.debugModules = cms.untracked.vstring('')
MessageLogger.cerr = cms.untracked.PSet(
FwkReport = cms.untracked.PSet(
limit = cms.untracked.int32(10000000),
reportEvery = cms.untracked.int32(10000)
),
threshold = cms.untracked.string('DEBUG'),
noLineBreaks = cms.untracked.bool(False),
DEBUG = cms.untracked.PSet(limit = cms.untracked.int32(0)),
INFO = cms.untracked.PSet(limit = cms.untracked.int32(0)),
Calibration = cms.untracked.PSet(limit = cms.untracked.int32(-1))
)
############################################################################
25 changes: 1 addition & 24 deletions DPGAnalysis/HcalTools/scripts/cmt/file_lists.csh
Original file line number Diff line number Diff line change
Expand Up @@ -10,34 +10,11 @@ echo " *************** myStart file_list *************** "

foreach i (`cat $1`)
echo ${i}
##
##
#dasgoclient --query="file dataset=/Cosmics/Commissioning2018-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file dataset=/MinimumBias/Commissioning2018-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
##

#dasgoclient --query="file dataset= /HcalNZS/Run2022A-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file dataset= /MinimumBias/Run2022A-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file dataset= /ZeroBias/Run2022A-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /ZeroBias/Run2022B-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2022A-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2022B-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2022C-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2022D-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /ZeroBias/Run2022E-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2022E-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2022F-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2022G-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}
#dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2022H-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}


#
##
dasgoclient --query="file site=T2_CH_CERN dataset= /HcalNZS/Run2023D-v1/RAW run=${i}" --limit=0 | sed "s/\/store/\'\/store/g" | sed "s/root/root\',/g"> TXT_${1}/run_${i}

end

echo "DONE: file_list.csh"

##


Loading