• Heavy Ion East meeting

  • Install cygwin on my laptap

  • analysis 14TeV p+p Gen data
    • I produced 14TeV p+p Gen data named pp14TeV_gen.root
    • To see the file I install ROOT and cygwin on my laptap (above method)
    • There is no data, totally empty
    • I guess the reason, data library not loaded when ROOT run.
    • Some of error Warnings are as follows
                                                 :
                                                 :
              Warning in : no dictionary for class HepMC::GenEvent is available
              Warning in : no dictionary for class HepMC::GenVertex is available
              Warning in : no dictionary for class HepMC::GenParticle is available
              Warning in : no dictionary for class HepMC::WeightContainer is available
              Warning in : no dictionary for class HepMC::FourVector is available
              Warning in : no dictionary for class HepMC::Flow is available
                                                 :
                                                 :

  • Contents of pp14TeV_gen.root
http://www.cms-kr.org/twiki/bin/viewfile/Main/20090220?rev=1;filename=2009-02-21_AM_02-31-23.jpg

  • Source code ; GenExample_cfg.py
    import FWCore.ParameterSet.Config as cms
    
    process = cms.Process("Gen")
    # this example configuration offers some minimum
    # annotation, to help users get through; please
    # don't hesitate to read through the comments
    # use MessageLogger to redirect/suppress multiple
    # service messages coming from the system
    #
    # in this config below, we use the replace option to make
    # the logger let out messages of severity ERROR (INFO level
    # will be suppressed), and we want to limit the number to 10
    #
    process.load("FWCore.MessageService.MessageLogger_cfi")
    # in you wish to quiet the logger even more and leave in
    # only info from the Framework, use config example below :
    #
    #service = MessageLogger
    #{
    #   untracked vstring destinations = {"cout"}
    #   untracked vstring categories = { "FwkJob" }
    #   untracked PSet cout =
    #   {
    #      untracked PSet default = { untracked int32 limit = 0 }  # kill all messages in the log
    #      untracked PSet FwkJob  = { untracked int32 limit = -1 } # but FwkJob category - those unlimitted
    #   }
    #}
    process.load("Configuration.StandardSequences.SimulationRandomNumberGeneratorSeeds_cff")
    
    process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
    
    # Vertex smearing
    # There're several types of Vertex smearing module (Gauss, Flat, BeamProfile, Betafunc, BetafuncEarlyCollisions);
    # the cff's are available for each one :
    # Configuration/StandardSequences/data/VtxSmearedGauss.cff
    # Configuration/StandardSequences/data/VtxSmearedFlat.cff
    # Configuration/StandardSequences/data/VtxSmearedBeamProfile.cff
    # Configuration/StandardSequences/data/VtxSmearedBetafuncNominalCollision.cff
    # Configuration/StandardSequences/data/VtxSmearedBetafuncEarlyCollision.cff
    # Currently, we suggest using Betafunc vertex smearing in the desired LHC configuration 
    process.load("Configuration.StandardSequences.VtxSmearedBetafuncEarlyCollision_cff")
    
    # Generator information and event analysis
    process.load("Configuration.StandardSequences.Generator_cff")
    
    # Event output
    process.load("Configuration.EventContent.EventContent_cff")
    
    process.maxEvents = cms.untracked.PSet(
        input = cms.untracked.int32(1)
    )
    process.source = cms.Source("FlatRandomEGunSource",
        PGunParameters = cms.untracked.PSet(
            # you can request more than 1 particle
            # since PartID is a vector, you can place in as many
            # PDG id's as you wish, comma seaparated
            #
            PartID = cms.untracked.vint32(2212),
            MaxEta = cms.untracked.double(2.7),
            MaxPhi = cms.untracked.double(3.14159265359),
            MinEta = cms.untracked.double(-2.7),
            MinE = cms.untracked.double(10.0),
            MinPhi = cms.untracked.double(-3.14159265359), ## in radians
    
            MaxE = cms.untracked.double(7000.0)
        ),
        Verbosity = cms.untracked.int32(0), ## set to 1 (or greater)  for printouts
    
        AddAntiParticle = cms.untracked.bool(True) ## back-to-back particles
    
    )
    
    process.FEVT = cms.OutputModule("PoolOutputModule",
        process.FEVTSIMEventContent,
        fileName = cms.untracked.string('pp14TeV_gen.root')
    )
    
    process.p0 = cms.Path(process.pgen)
    process.outpath = cms.EndPath(process.FEVT)
    process.schedule = cms.Schedule(process.p0,process.outpath) 

  • Result of cmsRun
    [seogoni@gate Python_event_gen]$ cmsenv
    [seogoni@gate Python_event_gen]$ cmsRun GenExample_cfg.py 
    Internal FlatRandomEGun is initialzed
    It is going to generate 1events
    found 539 particles
    Begin processing the 1st record. Run 1, Event 1, LumiSection 1 at 21-Feb-2009 10:43:19 KST
    
    =============================================
    
    MessageLogger Summary
    
    Severity    # Occurrences   Total Occurrences
    --------    -------------   -----------------
    [seogoni@gate Python_event_gen]$  

-- SeoKonKang - 21 Feb 2009

Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2009-02-21 - SeoKonKang
 
KOREA-CMS
This site is powered by the TWiki collaboration platform Powered by Perl This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding KoreaCmsWiki? Send feedback