-- YetkinYilmaz - 01-Apr-2011 Fix fo the calorimeter conditions in Data mixing.

SIM on Data

Problem

  • Ecal : energy loss by factor of 10 %
  • Hcal : no problem reported.

reasons and approach for fix.

  • Something wrong is going on the the digitization in the data mixer if the customized data condition is used.
  • The basic idea is to use MC global tag for the digitization of MC signal.

solution : Fix on top of Heavy ion event mixing recipe

Pre-processing for heavy ion background : no change.

Overlaying process.

  • run up to DIGI and drop the processes afterward as below.
# 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.datamixing_step = cms.Path(process.pdatamix)                                                                   
#process.L1simulation_step = cms.Path(process.SimL1Emulator)                                                            
#process.digi2raw_step = cms.Path(process.DigiToRaw)                                                                    
process.endjob_step = cms.Path(process.endOfProcess)                                                                   
process.RAWDEBUGHLToutput_step = cms.EndPath(process.RAWDEBUGHLToutput)
...
...

##And, change the schedule as below.
process.schedule = cms.Schedule(process.generation_step,process.simulation_step,process.digitisation_step)
process.schedule.extend([process.endjob_step,process.RAWDEBUGHLToutput_step])



  • Comment out the customise(process) => # customise(process)
    • You don't need to remove the definition of this ( the block of def customise(process): )
  • remove simHcalTTPDigis which will be recovered later. EDIT (4/6/11): is this necessary if DigiDM_cff.py is being used?
 
   process.pdigi.remove(process.simHcalTTPDigis)
  • Save all event content in output file. ( More effort is required to reduce the data size in this step )
   outputCommands = cms.untracked.vstring( "keep *"),

Reconstruction process.

  • Add MIX-DATAMIX-L1 before the Mixing
  • Add crossing_step to create the crossing frame.
    ===== newly added steps =======
    process.crossing_step = cms.Path(process.mix)
    process.datamixing_step = cms.Path(process.pdatamix)
    process.L1simulation_step = cms.Path(process.SimL1Emulator)
    process.digi2raw_step = cms.Path(process.DigiToRaw)
    ===== existing steps out of box ====
    process.raw2digi_step = cms.Path(process.RawToDigi)
    process.L1Reco_step = cms.Path(process.L1Reco)
    process.reconstruction_step = cms.Path(process.reconstructionHeavyIons)
    process.endjob_step = cms.Path(process.endOfProcess)
    process.FEVTDEBUGoutput_step = cms.EndPath(process.FEVTDEBUGoutput)
    
    
  • modify customiseCloneMatchRaw the customisation

def customiseCloneMatchRaw(process):
#    customiseCloneGenerator(process)           <============= comment this out. 
    customiseMatchRecVertex(process)
    customiseSiStripRawDigi(process)

    for outmod in process.outputModules_():
        print "Extending %s content to include 'hiSelectedVertex' from background event" % outmod
        getattr(process,outmod).outputCommands.extend(['keep *_hiSelectedVertex_*_*'])

    return process
  • modify customiseMatchRecVertex
    def customiseMatchRecVertex(process):
    process.load("SimGeneral.MixingModule.MatchRecVtx_cfi")
    process.matchRecVtx.heavyIonLabel = "hiSelectedVertex"
#    process.pgen.replace(process.VertexSmearing,process.matchRecVtx)                                                                     <<==== comment out
    process.mixData.input.sequential = True
    process.mixData.input.fileNames = process.source.fileNames
    return process

  • and then apply this at the very end of the configuration
process = customiseCloneMatchRaw(process)

Data on Data

Problem

  • pedestal noise was not removed in the data + data mixing step

solution : use the fixed code

   
cvs co UserCode/EYKim/dataMixer/DataMixingHcalDigiWorker.cc
mv UserCode/EYKim/dataMixer/DataMixingHcalDigiWorker.cc SimGeneral/DataMixingModule/plugins/
rm -r UserCode/EYKim

 
Edit | Attach | Watch | Print version | History: r6 < r5 < r4 < r3 < r2 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r6 - 2011-04-06 - EdWenger
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    CMSPublic All webs login

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