Ambiente cms

Creazione ambiente cms

  • scram project -n [cartella] CMSSW CMSSW_#versione
mettere al posto di [cartella] il percorso della cartella in cui installare CMSSW; sostituire al posto di #versione il numero della versione desiderata nella forma x_y_z.

p.e.:

scram project -n CMSSW536 CMSSW CMSSW_5_3_6

Inizializzazione cms

Entrare nella cartella src/ dell'ambiente CMSSW e dare il comando

  • cmsenv
p.e.:

cd CMSSW536/src

cmsenv

Creazione MC

Il procedimento composto di 2 parti:

1. Conversione campione lhe in ntpla root

Bisogna creare una ntpla root dal campione lhe eseguendo il seguente script di python.

Formato script python

creare "script lhe-root".py:

# Auto generated configuration file
# using: 
# Revision: 1.353 
# Source: /local/reps/CMSSW.admin/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v 
# with command line options: MCDBtoEDM --filein lhe:5739 --conditions START50_V13::All -s NONE --datatier GEN --no_exec
import FWCore.ParameterSet.Config as cms
process = cms.Process('LHE')

# 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('SimGeneral.MixingModule.mixNoPU_cfi')
process.load('Configuration.StandardSequences.GeometryDB_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.maxEvents = cms.untracked.PSet(

inserire nella seguente riga il numero di eventi da elaborare al posto di #eventi

    input = cms.untracked.int32(#eventi)
)

# Input source
process.source = cms.Source("LHESource",

sostituire a [file input] il percorso del file lhe di cui convertire gli eventi

fileNames = cms.untracked.vstring('file:[file input].lhe')
)
process.options = cms.untracked.PSet(
wantSummary = cms.untracked.bool(True)
)

# Production Info
process.configurationMetadata = cms.untracked.PSet(
    version = cms.untracked.string('$Revision: 1.2 $'),
    annotation = cms.untracked.string('MCDBtoEDM nevts:1'),
    name = cms.untracked.string('PyReleaseValidation')
)

# Output definition
process.RECOSIMoutput = cms.OutputModule("PoolOutputModule",
    splitLevel = cms.untracked.int32(0),
    eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
    outputCommands = process.RECOSIMEventContent.outputCommands,

sostituire al posto di [file output] il nome della ntpla da ottenere

    fileName = cms.untracked.string([file output].root'),
    dataset = cms.untracked.PSet(
        filterName = cms.untracked.string(''),
        dataTier = cms.untracked.string('GEN')
    )
)

# Additional output definition

# Other statements
process.GlobalTag.globaltag = 'START53_V7A::All'

# Path and EndPath definitions
process.RECOSIMoutput_step = cms.EndPath(process.RECOSIMoutput)

# Schedule definition
process.schedule = cms.Schedule(process.RECOSIMoutput_step)

dove "script lhe-root" e il nome dello script di python.

Lanciare lo script per la conversione

eseguire il comando

  • cmsRun "file lhe-root".py
(Ci vogliono circa 1 minuto di preparazione e 1 secondo ogni 1000 eventi per terminare l'operazione).

2. Generazione MC da ntpla root

Si genera infine il MC eseguendo il seguente script di python.

Formato script python

creare "script MC".py:

 # Auto generated configuration file
# using: 
# Revision: 1.381.2.13 
# Source: /local/reps/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v 
# with command line options: Configuration/GenProduction/python/EightTeV/Hadronizer_MgmMatchTuneZ2star_8TeV_madgraph_tauola_cff.py --step GEN,SIM,DIGI,L1,DIGI2RAW,RAW2DIGI,L1Reco,RECO --pileup 2012_Startup_50ns_PoissonOOTPU --conditions START53_V7C::All --beamspot Realistic8TeV2012Collision --eventcontent AODSIM --datatier AODSIM --no_exec --mc --process AODSIMoneStep --filein=file:MCDBtoEDM_NONE_thqhad_1_1_bPh.root
import FWCore.ParameterSet.Config as cms
process = cms.Process('AODSIMoneStep')

# 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('SimGeneral.MixingModule.mix_2012_Startup_50ns_PoissonOOTPU_cfi')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.GeometrySimDB_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeV2012Collision_cfi')
process.load('GeneratorInterface.Core.genFilterSummary_cff')
process.load('Configuration.StandardSequences.SimIdeal_cff')
process.load('Configuration.StandardSequences.Digi_cff')
process.load('Configuration.StandardSequences.SimL1Emulator_cff')
process.load('Configuration.StandardSequences.DigiToRaw_cff')
process.load('Configuration.StandardSequences.RawToDigi_cff')
process.load('Configuration.StandardSequences.L1Reco_cff')
process.load('Configuration.StandardSequences.Reconstruction_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.maxEvents = cms.untracked.PSet(



sostituire al posto di #eventi in numero di eventi da elaborare

    input = cms.untracked.int32(#eventi)
)
# Input source
process.source = cms.Source("PoolSource",
    secondaryFileNames = cms.untracked.vstring(),

mettere il nome della ntpla root di cui generare il MC al posto di [file input]

    fileNames = cms.untracked.vstring('file:[file input].root')
)
process.options = cms.untracked.PSet(
        wantSummary = cms.untracked.bool(True)
)
# Production Info
process.configurationMetadata = cms.untracked.PSet(
    version = cms.untracked.string('$Revision: 1.2 $'),
    annotation = cms.untracked.string('Configuration/GenProduction/python/EightTeV/Hadronizer_MgmMatchTuneZ2star_8TeV_madgraph_tauola_cff.py nevts:1'),
    name = cms.untracked.string('PyReleaseValidation')
)
process.mix.input.fileNames = cms.untracked.vstring(
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/00024E74-D96A-E111-9CE8-0030487D5DC3.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/003F446E-DC67-E111-A0F9-00266CFFA7A8.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/00DDA010-A769-E111-8A12-0030487D5EAF.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/02124719-B769-E111-8971-003048C66BBE.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/0261C34F-C969-E111-B0E6-003048C692E2.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/026B25B2-D967-E111-8561-00266CFFA5E0.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/0289AEA8-6165-E111-8512-0030487D5EB1.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/029CC6E3-9967-E111-95B2-00266CF32684.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/02C1D0AB-F267-E111-B434-003048F0E812.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/02E93ED0-6765-E111-885A-0030487F1F23.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/02F1619A-6E69-E111-A950-0030487E4EC7.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/0413ADCF-DD69-E111-85EC-003048C69414.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/04412411-4468-E111-9826-003048C66BBE.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/044C8BDF-516A-E111-AFFF-003048D436D2.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/044ED2FA-286B-E111-83D9-0030487EAFF9.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/04D7763A-2A6B-E111-A4C4-003048C66184.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/0646CB95-D567-E111-ACED-0025901D4C3C.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/067CF4F5-1568-E111-8B12-00266CF33288.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/06932BAC-BF69-E111-AE28-002481E14F2A.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/06972F6F-116B-E111-BC1A-0030487E0A2D.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0005E496-3661-E111-B31E-003048F0E426.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/003EEBD4-8061-E111-9A23-003048D437F2.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/005825F1-F260-E111-BD97-003048C692DA.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0065594C-B35E-E111-8B8C-003048C693EA.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0091FFD0-6B5E-E111-92FE-003048C693DA.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/00C69AE3-FE60-E111-BC48-0030487D8633.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02079692-AC61-E111-97BB-0025901D4D54.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/021BD915-2D61-E111-8BAD-002481E76052.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02386E4D-DC5E-E111-9413-00266CF1074C.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02446A08-515E-E111-82C7-00266CF330B8.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/024BBC60-CD5D-E111-A41A-0025901D4C32.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/025AB53D-F15E-E111-A7F4-003048C66BBE.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0274007A-A05E-E111-9227-003048D4DFB8.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02AC6853-E65D-E111-A515-0025901D4B04.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02B17485-B861-E111-84FC-003048C6930E.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02CFF160-B85E-E111-B547-003048D3CD92.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0407D67D-F360-E111-A737-0030487F1A73.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0407DB67-2F5E-E111-8115-003048D47976.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0464DA8A-4D61-E111-95F5-003048C64787.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/04670075-3561-E111-AD81-002481E0D50C.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/0015C0BD-0C64-E111-81F8-003048D43942.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/0060F354-3364-E111-A5C7-003048D43960.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/00852170-E061-E111-BA8F-003048D3CDE0.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/00B6AAA8-0D65-E111-966C-0025901D493A.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/00EE834E-3464-E111-9B62-003048D43960.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/00EEF21D-0762-E111-9A59-003048F0E1B2.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02163E93-0A65-E111-B2DE-00266CF32E78.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/021CC150-0A65-E111-90F2-00266CF33340.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/0220793B-9862-E111-A38C-003048CF6338.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02370007-2062-E111-98D4-0030487D5DC3.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/025EC033-EC63-E111-A9D4-0025901D4B22.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/026FA1DC-5F64-E111-8D64-003048D439B4.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/029C9303-0A65-E111-B8DE-0025901D4B22.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02B30A54-C561-E111-AE69-003048D43700.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02B60BCF-DF63-E111-906F-003048C693E6.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02BC9C04-BD64-E111-9B60-003048D436CA.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02CE813D-0A65-E111-9721-0025901D4A58.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02D8803D-EE63-E111-ABE4-0030487F1657.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/0400B93B-E264-E111-B99D-003048C6930E.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/04067303-0A65-E111-BAD0-0025901D4ADE.root',
)

# Output definition
process.AODSIMoutput = cms.OutputModule("PoolOutputModule",
  eventAutoFlushCompressedSize = cms.untracked.int32(15728640),
  outputCommands = process.AODSIMEventContent.outputCommands,



sostituire a [file output] il percorso del file root MC

  fileName = cms.untracked.string('[file MC].root'),
  dataset = cms.untracked.PSet(
    filterName = cms.untracked.string(''),
    dataTier = cms.untracked.string('AODSIM')
  ),
SelectEvents = cms.untracked.PSet(
SelectEvents = cms.vstring('generation_step')
  )
)

# Additional output definition

# Other statements
process.genstepfilter.triggerConditions=cms.vstring("generation_step")
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, 'START53_V7C::All', '')
process.generator = cms.EDFilter("Pythia6HadronizerFilter",
ExternalDecays = cms.PSet(
    Tauola = cms.untracked.PSet(
UseTauolaPolarization = cms.bool(True),
InputCards = cms.PSet(
        mdtau = cms.int32(0),
        pjak2 = cms.int32(0),
        pjak1 = cms.int32(0)
      )
    ),
    parameterSets = cms.vstring('Tauola')
  ),
UseExternalGenerators = cms.untracked.bool(True),
  pythiaPylistVerbosity = cms.untracked.int32(1),
  pythiaHepMCVerbosity = cms.untracked.bool(True),
  comEnergy = cms.double(8000.0),
  jetMatching = cms.untracked.PSet(
        MEMAIN_showerkt = cms.double(0),<span style="font-size: 12pt;">        MEMAIN_maxjets = cms.int32(2),</span>
        MEMAIN_minjets = cms.int32(0),
        MEMAIN_qcut = cms.double(4),
        MEMAIN_excres = cms.string(''),
        MEMAIN_etaclmax = cms.double(5),
        MEMAIN_nqmatch = cms.int32(5),
        outTree_flag = cms.int32(0),
        scheme = cms.string('Madgraph'),
        mode = cms.string('auto')
    ),
    maxEventsToPrint = cms.untracked.int32(0),
    PythiaParameters = cms.PSet(
        pythiaUESettings = cms.vstring('MSTU(21)=1     ! Check on possible errors during program execution', 
            'MSTJ(22)=2     ! Decay those unstable particles', 
            'PARJ(71)=10 .  ! for which ctau  10 mm', 
            'MSTP(33)=0     ! no K factors in hard cross sections', 
            'MSTP(2)=1      ! which order running alphaS', 
            'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)', 
            'MSTP(52)=2     ! work with LHAPDF', 
            'PARP(82)=1.921 ! pt cutoff for multiparton interactions', 
            'PARP(89)=1800. ! sqrts for which PARP82 is set', 
            'PARP(90)=0.227 ! Multiple interactions: rescaling power', 
            'MSTP(95)=6     ! CR (color reconnection parameters)', 
            'PARP(77)=1.016 ! CR', 
            'PARP(78)=0.538 ! CR', 
            'PARP(80)=0.1   ! Prob. colored parton from BBR', 
            'PARP(83)=0.356 ! Multiple interactions: matter distribution parameter', 
            'PARP(84)=0.651 ! Multiple interactions: matter distribution parameter', 
            'PARP(62)=1.025 ! ISR cutoff', 
            'MSTP(91)=1     ! Gaussian primordial kT', 
            'PARP(93)=10.0  ! primordial kT-max', 
            'MSTP(81)=21    ! multiple parton interactions 1 is Pythia default', 
            'MSTP(82)=4     ! Defines the multi-parton model'),

configurare qui come pythia deve effettuare i decadimenti

        processParameters = cms.vstring('MSEL=0         ! User defined processes', 
            'MWID(25) = 2',
            'MDME(210,1)=0           !Higgs decay into dd',
            'MDME(211,1)=0           !Higgs decay into uu',
            'MDME(212,1)=0           !Higgs decay into ss',
            'MDME(213,1)=0           !Higgs decay into cc',
            'MDME(214,1)=0           !Higgs decay into bb',
            'MDME(215,1)=0           !Higgs decay into tt',
            'MDME(216,1)=0           !Higgs decay into',
            'MDME (217,1)= 0            !Higgs decay into Higgs decay',
            'MDME (218,1)= 0            !Higgs decay into e nu e',
            'MDME (219,1)= 0            !Higgs decay into mu nu mu',
            'MDME (220,1)= 1            !Higgs decay into tau nu tau',
            'MDME (221,1)= 0            !Higgs decay into Higgs decay',
            'MDME (222,1)= 0            !Higgs decay into g g',
            'MDME (223,1)= 1            !Higgs decay into gam gam',
            'MDME (224,1)= 0            !Higgs decay into gam Z',
            'MDME (225,1)= 0            !Higgs decay into Z Z',
            'MDME (226,1)= 1    !Higgs decay into W W',
            'PMAS(5,1)=4.8   ! b quark mass', 
            'PMAS(6,1)=173   ! t quark mass', 
            'MSTJ(1)=1       ! Fragmentation/hadronization on or off', 
            'MSTP(61)=1      ! Parton showering on or off',
            'BRAT(210)=0.0     !h decay into d dbar',
            'BRAT(211)=0.0     !h decay into u ubar',
            'BRAT(212)=0.0     !h decay into s sbar',
            'BRAT(213)=0.0     !h decay into c cbar',
            'BRAT(214)=0.0     !h decay into b bbar',
            'BRAT(215)=0.0     !h decay into t tbar',
            'BRAT(216)=0.0     !h decay into b4 b4bar',
            'BRAT(217)=0.0     !h decay into t4 t4bar',
            'BRAT(218)=0.0     !h decay into e+ e-',
            'BRAT(219)=0.0     !h decay into mu+ mu-',
            'BRAT(220)=0.204     !h decay into tau+ tau-',
            'BRAT(221)=0.0     !h decay into tau4+ tau4-',
            'BRAT(222)=0.0     !h decay into g g',
            'BRAT(223)=0.1     !h decay into gamma gamma',
            'BRAT(224)=0.0     !h decay into gamma Z',
            'BRAT(225)=0.0     !h decay into Z Z',
            'BRAT(226)=0.696     !h decay into W+ W-'),
        parameterSets = cms.vstring('pythiaUESettings', 
            'processParameters')
    )
)
# Path and EndPath definitions
process.generation_step = cms.Path(process.pgen)
process.simulation_step = cms.Path(process.psim)
process.digitisation_step = cms.Path(process.pdigi)
process.L1simulation_step = cms.Path(process.SimL1Emulator)
process.digi2raw_step = cms.Path(process.DigiToRaw)
process.raw2digi_step = cms.Path(process.RawToDigi)
process.L1Reco_step = cms.Path(process.L1Reco)
process.reconstruction_step = cms.Path(process.reconstruction)
process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
process.endjob_step = cms.EndPath(process.endOfProcess)
process.AODSIMoutput_step = cms.EndPath(process.AODSIMoutput)
# Schedule definition
process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.simulation_step,process.digitisation_step,process.L1simulation_step,process.digi2raw_step,process.raw2digi_step,process.L1Reco_step,process.reconstruction_step,process.endjob_step,process.AODSIMoutput_step)
# filter all path with the production filter sequence
for path in process.paths:
        getattr(process,path)._seq = process.generator * getattr(process,path)._seq

dove "script MC" il nome dello script python.

Lanciare lo script per la generazione del MC

eseguire il comando

  • cmsRun "script MC".py
(Ci vogliono circa 3 minuti di preparazione e 3 minuti ad evento).

Sistema crab

Con il sistema crab si pu generare il MC dividendo gli eventi da elaborare tra jobs eseguiti da pi macchine.

Creazione del sistema crab

Per poter utilizzare crab necessario disporre delle autorizzazioni necessarie per l'accesso a tutti i file remoti richiesti dal sistema.

Richiesta del certificato CA per GRID

Aprire le impostazioni del proprio CERN Account Management, e richiedere il certificato CA per GRID cliccando su Request a CERN Grid Certification Authority Certificate.

-- MarcoLaRosa - 20 Jan 2014

Edit | Attach | Watch | Print version | History: r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r2 - 2014-01-22 - MarcoLaRosa
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback