HcalCalibAlgos/test/batchjob_calib_valid.csh

0001 #!/bin/csh
0002 #BSUB -q 1nh
0003 ##BSUB -J jfp
0004
0005 echo 'Start'
0006 #export STAGE_SVCCLASS cmscaf
0007
0008 #cmsrel CMSSW_3_1_0_pre10
0009 #cd CMSSW_3_1_0_pre10/src
0010 cmsrel CMSSW_3_1_4
0011 cd CMSSW_3_1_4/src
0012 cmsenv
0013 cvs co Calibration/HcalCalibAlgos
0014 cp /afs/cern.ch/user/a/andrey/scratch1/CMSSW_3_1_4/src/Calibration/HcalCalibAlgos/plugins/ValidIsoTrkCalib.cc Calibration/HcalCalibAlgos/plugins
0015 scram b
0016 cd Calibration/HcalCalibAlgos/test
0017
0018 #set respcorrdir=/afs/cern.ch/user/a/andrey/scratch1/CMSSW_3_1_0_pre10/src/Calibration/HcalCalibAlgos/data
0019 set respcorrdir=/afs/cern.ch/user/a/andrey/scratch1/CMSSW_3_1_4/src/Calibration/HcalCalibAlgos/data
0020
0021 # if you want to validate your own calibration, copy it to data/ from your local place:
0022 cp $respcorrdir/calibConst_IsoTrk_testCone_26.3cm.txt ../data/response_corrections.txt
0023
0024 cat > validator.py <<@EOF
0025
0026 import FWCore.ParameterSet.Config as cms
0027
0028 process = cms.Process("Validator")
0029 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0030 #process.GlobalTag.globaltag = "IDEAL_31X::All"
0031 process.GlobalTag.globaltag = "MC_31X_V5::All"
0032 process.prefer("GlobalTag")
0033
0034 process.load("Configuration.StandardSequences.VtxSmearedBetafuncEarlyCollision_cff")
0035 process.load("Configuration.StandardSequences.Generator_cff")
0036 process.load("Configuration.StandardSequences.GeometryDB_cff")
0037 process.load("Configuration.StandardSequences.MagneticField_cff")
0038 process.load("Calibration.HcalAlCaRecoProducers.ALCARECOHcalCalIsoTrk_cff")
0039 process.load("Calibration.HcalAlCaRecoProducers.ALCARECOHcalCalIsoTrkNoHLT_cff")
0040 process.load("Calibration.HcalAlCaRecoProducers.alcaisotrk_cfi")
0041 #process.IsoProd.SkipNeutralIsoCheck = cms.untracked.bool(True)
0042 process.IsoProd.MinTrackP = cms.double(5.0)
0043
0044 process.isoHLT.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT")
0045 process.load("Configuration.StandardSequences.Services_cff")
0046 process.load("Configuration.StandardSequences.Reconstruction_cff")
0047 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0048 process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(5000)
0049
0050 process.load("Calibration.HcalCalibAlgos.calib_validator_cfi")
0051 process.ValidationIsoTrk.outputFileName = cms.string("ValidFile_XX.root")
0052 process.ValidationIsoTrk.AxB = cms.string("Cone")
0053 #process.ValidationIsoTrk.AxB = cms.string("3x3")
0054 process.ValidationIsoTrk.calibrationConeSize = cms.double(26.2)
0055 process.es_ascii2 = cms.ESSource("HcalTextCalibrations",
0056  appendToDataLabel = cms.string('recalibrate'),
0057   input = cms.VPSet(
0058   cms.PSet(object = cms.string('RespCorrs'),
0059     file = cms.FileInPath('Calibration/HcalCalibAlgos/data/response_corrections.txt')
0060           ),
0061    )
0062 )
0063
0064 #process.ValidationIsoTrk.outputFileName = cms.string("ValidFile_10_${1}.root")
0065 process.ValidationIsoTrk.outputFileName = cms.string("ValidFile_50_${1}.root")
0066 #process.ValidationIsoTrk.outputFileName = cms.string("ValidFile_100_${1}.root")
0067 #process.ValidationIsoTrk.outputFileName = cms.string("ValidFile_300_${1}.root")
0068
0069 process.load("HLTrigger.Timer.timer_cfi")
0070
0071 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0072 process.source = cms.Source("PoolSource",
0073 fileNames = cms.untracked.vstring(
0074
0075 'rfio:/castor/cern.ch/user/a/abdullin/pi50_fullproduction_312/pi50_${1}.root'
0076
0077 #'rfio:/castor/cern.ch/user/a/abdullin/pi50_fullproduction_310pre10/pi50_${1}.root'
0078 #'rfio:/castor/cern.ch/user/a/abdullin/pi50_fullproduction_310pre10/pi50_HEZS8_${1}.root'
0079 #'rfio:/castor/cern.ch/user/a/abdullin/pi300_fullproduction_310pre10/pi300_${1}.root'
0080
0081      )
0082 )
0083 process.options = cms.untracked.PSet(
0084    wantSummary = cms.untracked.bool(True)
0085 )
0086 process.TimerService = cms.Service("TimerService", useCPUtime = cms.untracked.bool(True))
0087 process.pts = cms.EDFilter("PathTimerInserter")
0088 process.PathTimerService = cms.Service("PathTimerService")
0089 process.p = cms.Path(process.seqALCARECOHcalCalIsoTrkNoHLT*process.ValidationIsoTrk)
0090 #process.p = cms.Path(process.ValidationIsoTrk)
0091
0092 #-----------
0093 @EOF
0094
0095 cmsRun validator.py
0096
0097 set outdir=/afs/cern.ch/user/a/andrey/scratch1/CMSSW_3_1_4/src/Calibration/HcalCalibAlgos/test
0098 #set outdir=/afs/cern.ch/user/a/andrey/scratch1/CMSSW_3_1_0_pre10/src/Calibration/HcalCalibAlgos/test
0099 #set outdir=/castor/cern.ch/user/a/andrey/pi50_310_pre10
0100 #set outdir=/castor/cern.ch/user/a/andrey/pi300_310_pre10
0101
0102
0103 cp ValidFile_*.root $outdir/
0104 #rfcp ValidFile_*.root $outdir/
0105
0106 cd ../../../../../
0107
0108 rm -rf CMSSW*