Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-03-17 11:00:49

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 process = cms.Process("nu")
0004 
0005 # The number of events to be processed.
0006 process.maxEvents = cms.untracked.PSet(
0007     input = cms.untracked.int32(1000000)
0008 )
0009 
0010 
0011 # Include the RandomNumberGeneratorService definition
0012 process.load("IOMC.RandomEngine.IOMC_cff")
0013 
0014 # For histograms
0015 process.load("DQMServices.Core.DQM_cfg")
0016 
0017 # Input
0018 process.source = cms.Source(
0019     "PoolSource",
0020     debugFlag = cms.untracked.bool(True),
0021     debugVebosity = cms.untracked.uint32(10),
0022     fileNames = cms.untracked.vstring(
0023         'rfio:/castor/cern.ch/user/p/pjanot/CMSSW170pre12/SingleParticlePID211-E50/SingleParticlePID211-E50_0.root'
0024     )
0025 )
0026 
0027 process.testNU = cms.EDFilter(
0028     "testNuclearInteractions",
0029     # Full, Fast radii and lengths for plots
0030     BPCylinderRadius = cms.untracked.vdouble(3.2, 3.05),
0031     BPCylinderLength = cms.untracked.vdouble(999.0, 28.3),
0032 
0033     PXB1CylinderRadius = cms.untracked.vdouble(6.0, 6.0),
0034     PXB1CylinderLength = cms.untracked.vdouble(28.2, 28.3),
0035 
0036     PXB2CylinderRadius = cms.untracked.vdouble(8.5, 8.5),
0037     PXB2CylinderLength = cms.untracked.vdouble(28.2, 28.3),
0038 
0039     PXB3CylinderRadius = cms.untracked.vdouble(11.5, 11.5),
0040     PXB3CylinderLength = cms.untracked.vdouble(28.2, 28.3),
0041 
0042     PXBCablesCylinderRadius = cms.untracked.vdouble(16.9, 16.9),
0043     PXBCablesCylinderLength = cms.untracked.vdouble(30.0, 30.0),
0044 
0045     PXD1CylinderRadius = cms.untracked.vdouble(16.0, 17.0),
0046     PXD1CylinderLength = cms.untracked.vdouble(40.0, 40.0),
0047 
0048     PXD2CylinderRadius = cms.untracked.vdouble(16.0, 17.0),
0049     PXD2CylinderLength = cms.untracked.vdouble(52.0, 52.0),
0050 
0051     PXDCablesCylinderRadius = cms.untracked.vdouble(22.0, 20.2),
0052     PXDCablesCylinderLength = cms.untracked.vdouble(999.0, 100.0),
0053 
0054     TIB1CylinderRadius = cms.untracked.vdouble(30.0, 28.0),
0055     TIB1CylinderLength = cms.untracked.vdouble(68.0, 69.0),
0056 
0057     TIB2CylinderRadius = cms.untracked.vdouble(38.0, 37.0),
0058     TIB2CylinderLength = cms.untracked.vdouble(68.0, 69.0),
0059 
0060     TIB3CylinderRadius = cms.untracked.vdouble(46.0, 45.0),
0061     TIB3CylinderLength = cms.untracked.vdouble(68.0, 69.0),
0062 
0063     TIB4CylinderRadius = cms.untracked.vdouble(54.0, 53.0),
0064     TIB4CylinderLength = cms.untracked.vdouble(68.0, 69.0),
0065 
0066     TIBCablesCylinderRadius = cms.untracked.vdouble(53.0, 54.0),
0067     TIBCablesCylinderLength = cms.untracked.vdouble(74.0, 75.0),
0068 
0069     TID1CylinderRadius = cms.untracked.vdouble(52.0, 54.0),
0070     TID1CylinderLength = cms.untracked.vdouble(83.0, 83.0),
0071 
0072     TID2CylinderRadius = cms.untracked.vdouble(52.0, 54.0),
0073     TID2CylinderLength = cms.untracked.vdouble(95.0, 95.0),
0074 
0075     TID3CylinderRadius = cms.untracked.vdouble(52.0, 54.0),
0076     TID3CylinderLength = cms.untracked.vdouble(110.0, 106.0),
0077 
0078     TIDCablesCylinderRadius = cms.untracked.vdouble(58.2, 55.5),
0079     TIDCablesCylinderLength = cms.untracked.vdouble(122.0, 108.5),
0080 
0081     TOB1CylinderRadius = cms.untracked.vdouble(66.6, 65.0),
0082     TOB1CylinderLength = cms.untracked.vdouble(109.0, 109.0),
0083 
0084     TOB2CylinderRadius = cms.untracked.vdouble(75.5, 75.0),
0085     TOB2CylinderLength = cms.untracked.vdouble(109.0, 109.0),
0086 
0087     TOB3CylinderRadius = cms.untracked.vdouble(84.1, 83.0),
0088     TOB3CylinderLength = cms.untracked.vdouble(109.0, 109.0),
0089 
0090     TOB4CylinderLength = cms.untracked.vdouble(109.0, 109.0),
0091     TOB4CylinderRadius = cms.untracked.vdouble(94.0, 92.0),
0092 
0093     TOB5CylinderRadius = cms.untracked.vdouble(105.5, 103.0),
0094     TOB5CylinderLength = cms.untracked.vdouble(109.0, 109.0),
0095 
0096     TOB6CylinderRadius = cms.untracked.vdouble(113.0, 113.0),
0097     TOB6CylinderLength = cms.untracked.vdouble(109.0, 109.0),
0098 
0099     TOBCablesCylinderRadius = cms.untracked.vdouble(113.0, 113.0),
0100     TOBCablesCylinderLength = cms.untracked.vdouble(125.0, 125.0),
0101 
0102     TEC1CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0103     TEC1CylinderLength = cms.untracked.vdouble(138.0, 138.0),
0104 
0105     TEC2CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0106     TEC2CylinderLength = cms.untracked.vdouble(152.0, 152.0),
0107 
0108     TEC3CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0109     TEC3CylinderLength = cms.untracked.vdouble(166.0, 166.0),
0110 
0111     TEC4CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0112     TEC4CylinderLength = cms.untracked.vdouble(180.0, 180.0),
0113 
0114     TEC5CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0115     TEC5CylinderLength = cms.untracked.vdouble(195.0, 195.0),
0116 
0117     TEC6CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0118     TEC6CylinderLength = cms.untracked.vdouble(212.0, 212.0),
0119 
0120     TEC7CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0121     TEC7CylinderLength = cms.untracked.vdouble(232.0, 232.0),
0122 
0123     TEC8CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0124     TEC8CylinderLength = cms.untracked.vdouble(254.0, 254.0),
0125     
0126     TEC9CylinderRadius = cms.untracked.vdouble(110.0, 110.0),
0127     TEC9CylinderLength = cms.untracked.vdouble(272.0, 272.0),
0128 
0129     TrackerCablesCylinderRadius = cms.untracked.vdouble(125.0, 121.0),
0130     TrackerCablesCylinderLength = cms.untracked.vdouble(301.0, 301.0),
0131 
0132     TestParticleFilter = cms.PSet(
0133         # Particles with |eta| > etaMax (momentum direction at primary vertex) 
0134         # are not simulated 
0135         etaMax = cms.double(5.0),
0136         # Charged particles with pT < pTMin (GeV/c) are not simulated
0137         pTMin = cms.double(0.0),
0138         # Particles with energy smaller than EMin (GeV) are not simulated
0139         EMin = cms.double(0.0),
0140         # Protons with energy in excess of this value (GeV) will kept no matter what
0141         EProton = cms.double(99999.0),
0142     )
0143 )
0144 
0145 # Famos SimHits 
0146 process.load("FastSimulation.Configuration.CommonInputsFake_cff")
0147 process.load("FastSimulation.Configuration.FamosSequences_cff")
0148 # Use CMSSW_1_7_0 tuning
0149 process.misalignedTrackerInteractionGeometry.TrackerMaterial.TrackerMaterialVersion = 2
0150 # No vertex smearing
0151 process.famosSimHits.VertexGenerator.BetaStar = 0.0000001
0152 process.famosSimHits.VertexGenerator.SigmaZ = 0.00001
0153 # Magnetic field
0154 process.load("Configuration.StandardSequences.MagneticField_40T_cff")
0155 #process.load("Configuration.StandardSequences.MagneticField_38T_cff")
0156 process.VolumeBasedMagneticFieldESProducer.useParametrizedTrackerField = True
0157 # No SimHits
0158 process.famosSimHits.SimulateCalorimetry = False
0159 process.famosSimHits.SimulateTracking = False
0160 
0161 # Path to run what is needed
0162 process.p = cms.Path(
0163     process.offlineBeamSpot+
0164     process.famosPileUp+
0165     process.famosSimHits+
0166     process.testNU
0167 )
0168 
0169 # Keep the logging output to a nice level #
0170 process.load("FWCore.MessageService.MessageLogger_cfi")
0171 process.MessageLogger.cerr.enable = False
0172 process.MessageLogger.files.test = dict(extension = 'txt')
0173 
0174