-- AmnonHarel - 23 Feb 2009


See also AmnonHarelDataMixingRound2 (Mike's "V2" and "V3" files) and AmnonHarelDataMixing.

Testing Mike's DataMixer, circa 20 Feb, with his "V4" files. Goal was to see a plot of HPD discharge noise overlayed near / over MC signal.

Conclusions (on the overlaying process, and reconstruction)


  • work area: /uscms/home/aharel/CMSSW_3_1_0_pre/src/RecoJets/JetAnalyzers
    • Its setup is described below
  • data area: contains both .log and .root files /uscms_data/d1/aharel
  • v4 MC events: /uscms_data/d1/mikeh/QCD_Pt_50_80_cfi_GEN_SIM_DIGI_L1_DIGI2RAW_NoNoise.root
    • examined using gdQCD4new.py --> gdQCD4.newlog * this file contains calorimeter digis, but no reco. calo jets
    • raw2digi, reco & translate to my sc.c format using prepQCD4.py --> prepQCD4.log, prepQCD4events.root (in CMSSW format) & prepQCD4.root (in my format)
      • see comments about raw2digi below
  • Non-collision data events: /uscms_data/d1/mikeh/promptRecoHCalNoise.root (from run 68021)
    • examined using gdData4.py --> gdData4.log
    • this file contains calorimeter digis, but no reco. calo jets
    • raw2digi, reco & translate to my sc.c format using prepData4.py --> prepData4.log, prepData4event.root (in CMSSW format) & prepData4.root (in my format)
      • Finally update to a CRAFT global tag: CRAFT_30X::All. Still pretending the magnetic field is 0 though (don't really care).
      • Looking for anomolous events:
        • #14 (i.e. 645870) looks odd at HB-, 70. Very long low signal (like RBX) but only a few channels in the same HPD.
        • #28 (i.e. 1082063) is notable for having a very short energy deposit (confirmed by ECal, so probably a malformed muon deposit)
        • #364 is the most HPD-discharge like, and there's also no ECal energy. Yet, it just looks like a muon hitting HCal.
        • all the rest look like uncorrelated noise + muons.
      • Not sure why the raw2digi step is needed. Without it complains about:
              Looking for type: edm::SortedCollection<ZDCDataFrame,edm::StrictWeakOrdering<ZDCDataFrame> >
              Looking for module label: hcalDigis
              Looking for productInstanceName:
  • Overlayed events:
    • who's where: /uscms_data/d1/mikeh/Craft_on_QCD_v4_eventlist.txt
    • "processed to Digi stage": /uscms_data/d1/mikeh/CRAFT_on_QCD_NewHcal_Digis_v4.root
      • examined using gdO4Digi.py --> gdO4Digi.log
      • this file contains mixer calorimeter digis and no reconstructed jets (as expected)
    • "Reco version": /uscms_data/d1/mikeh/CRAFT_on_QCD_NewHcal_RECO_v4.root
      • examined using gdO4Reco.py --> gdO4Reco.log
      • this file contains reconstructed jets, and mixer calorimeter digis smile
      • dump in my sc.c format for plots with dumpMix4.py --> dumpMix4.root
        • updated DigiInputTag, it's now "mix"

  • new work area:
    • was set up with:
cmsrel CMSSW_3_1_0_pre1
cd CMSSW_3_1_0_pre1/src
addpkg RecoJets
cd RecoJets/JetAnalyzers/interface/
ln -s ~/test/AHtest.h AHtest.h
cd ../src
ln -s ~/test/AHtest.cc AHtest.cc
cd ..
mkdir python
ln -s ~/test/AHtest_cfi.py python/AHtest_cfi.py
mv BuildFile BuildFile_org
ln -s ~/test/BuildFile_AHtest_CMSSW_3_1_0_pre1 BuildFile
scram b
Edit | Attach | Watch | Print version | History: r4 < r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r2 - 2009-02-23 - AmnonHarel
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

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