Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:31:59

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 from $ERA_MOD_PATH.Era_$ERA_cff import *
0004 process = cms.Process('CTPPSTest', $ERA)
0005 
0006 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0007 process.load('Validation.CTPPS.ctppsLHCInfoPlotter_cfi')
0008 process.load('Configuration.Generator.randomXiThetaGunProducer_cfi')
0009 process.load("CondCore.CondDB.CondDB_cfi")
0010 
0011 # minimal logger settings
0012 process.MessageLogger = cms.Service("MessageLogger",
0013     statistics = cms.untracked.vstring(),
0014     destinations = cms.untracked.vstring('cout'),
0015     cout = cms.untracked.PSet(
0016         threshold = cms.untracked.string('WARNING')
0017     )
0018 )
0019 
0020 # particle generator
0021 process.generator.xi_max = 0.25
0022 process.generator.theta_x_sigma = 60.e-6
0023 process.generator.theta_y_sigma = 60.e-6
0024 
0025 # default source
0026 process.source = cms.Source("EmptySource",
0027     firstRun = cms.untracked.uint32(1),
0028 )
0029 
0030 process.CondDB.connect = 'frontier://FrontierProd/CMS_CONDITIONS'
0031 process.PoolDBESSource = cms.ESSource("PoolDBESSource",
0032     process.CondDB,
0033     toGet = cms.VPSet(cms.PSet(
0034         record = cms.string('CTPPSPixelAnalysisMaskRcd'),
0035         tag = cms.string("CTPPSPixelAnalysisMask_Run3_v1_hlt"))
0036         ))
0037 
0038 # random seeds
0039 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0040     sourceSeed = cms.PSet(initialSeed = cms.untracked.uint32(98765)),
0041     generator = cms.PSet(initialSeed = cms.untracked.uint32(98766)),
0042     beamDivergenceVtxGenerator = cms.PSet(initialSeed = cms.untracked.uint32(3849)),
0043     ppsDirectProtonSimulation = cms.PSet(initialSeed = cms.untracked.uint32(4981))
0044 )
0045 
0046 # number of events
0047 process.maxEvents = cms.untracked.PSet(
0048     input = cms.untracked.int32(int($N_EVENTS))
0049 )
0050 
0051 # LHCInfo plotter
0052 process.ctppsLHCInfoPlotter.outputFile = "$OUT_LHCINFO"
0053 
0054 # track distribution plotter
0055 process.ctppsTrackDistributionPlotter = cms.EDAnalyzer("CTPPSTrackDistributionPlotter",
0056     tagTracks = cms.InputTag("ctppsLocalTrackLiteProducer"),
0057     outputFile = cms.string("$OUT_TRACKS")
0058 )
0059 
0060 # reconstruction plotter
0061 process.ctppsProtonReconstructionPlotter = cms.EDAnalyzer("CTPPSProtonReconstructionPlotter",
0062     tagTracks = cms.InputTag("ctppsLocalTrackLiteProducer"),
0063     tagRecoProtonsSingleRP = cms.InputTag("ctppsProtons", "singleRP"),
0064     tagRecoProtonsMultiRP = cms.InputTag("ctppsProtons", "multiRP"),
0065     outputFile = cms.string("$OUT_PROTONS")
0066 )
0067 
0068 process.generation = cms.Path(process.generator)
0069 
0070 process.validation = cms.Path(
0071     process.ctppsLHCInfoPlotter
0072     * process.ctppsTrackDistributionPlotter
0073     * process.ctppsProtonReconstructionPlotter
0074 )
0075 
0076 # processing path
0077 process.schedule = cms.Schedule(
0078     process.generation,
0079     process.validation
0080 )
0081 
0082 from SimPPS.Configuration.Utils import setupPPSDirectSim
0083 setupPPSDirectSim(process)
0084 
0085 process.ctppsBeamParametersFromLHCInfoESSource.vtxOffsetX45 = 0.
0086 process.ctppsBeamParametersFromLHCInfoESSource.vtxOffsetY45 = 0.
0087 process.ctppsBeamParametersFromLHCInfoESSource.vtxOffsetZ45 = 0.
0088 process.source.numberEventsInLuminosityBlock = process.ctppsCompositeESSource.generateEveryNEvents
0089 process.ctppsTrackDistributionPlotter.rpId_45_F = process.rpIds.rp_45_F
0090 process.ctppsTrackDistributionPlotter.rpId_45_N = process.rpIds.rp_45_N
0091 process.ctppsTrackDistributionPlotter.rpId_56_N = process.rpIds.rp_56_N
0092 process.ctppsTrackDistributionPlotter.rpId_56_F = process.rpIds.rp_56_F
0093 process.ctppsProtonReconstructionPlotter.rpId_45_F = process.rpIds.rp_45_F
0094 process.ctppsProtonReconstructionPlotter.rpId_45_N = process.rpIds.rp_45_N
0095 process.ctppsProtonReconstructionPlotter.rpId_56_N = process.rpIds.rp_56_N
0096 process.ctppsProtonReconstructionPlotter.rpId_56_F = process.rpIds.rp_56_F