2006 Ecal H4 simulation notes

Recipe to install the software

using "ECALTBH4, version: 0_1_0 + Configuration/EcalTB V00-00-06"

  scramv1 project -n ECALTBH4_0_1_0 CMSSW CMSSW_0_8_1
  cd ECALTBH4_0_1_0
  PackageManagement.pl --rel ECALTBH4_0_1_0
  cd src
  project CMSSW
  cvs update -r V00-00-06 Configuration/EcalTB
  scramv1 b

Get other files

  wget 'http://ecaltbh4sim.web.cern.ch/ecaltbh4sim/prod_run.csh'

Running the simulation

usage:
   ./prod_run.csh config_file beam_energy number_of_events
The config file is searched for in $CMSSW_BASE/src/SimG4CMS/EcalTestBeam/test/data/

e.g.

./prod_run.csh crystal1104.cff 120 5

Useful root macro snippets for SimHits

Events->SetAlias("simhit","PCaloHits_SimG4Object_EcalHitsEB_ECALTBH4.obj");
Events->SetAlias("simhit_iphi","simhit.id() & 0x1ff");
Events->SetAlias("simhit_ieta","(simhit.id()  >> 9) & 0x7f");
Events->SetAlias("genpos","PEcalTBInfo_SimEcalTBG4Object__ECALTBH4.obj")

For Crystal 1104:

Events->Draw("simhit.energy()","simhit_iphi == 17 && simhit_ieta == 56")

Drawing the simulated beam profile:

Events->Draw("genpos.evYbeam():genpos.evXbeam()")

Energy and positions (for mips):

Events->Draw("simhit.energy():genpos.evXbeam()","simhit_iphi == 17 && simhit_ieta == 56 && simhit.energy() < 1");
Events->Draw("simhit.energy():genpos.evYbeam():genpos.evXbeam()","simhit_iphi == 17 && simhit_ieta == 56 && simhit.energy() < 1");
Events->Draw("simhit.energy()","simhit_iphi == 17 && simhit_ieta == 56 && simhit.energy() < 1 && abs(genpos.evXbeam()) < 5 && abs(genpos.evYbeam()) < 5");

-- AndreHolzner - 12 Sep 2006

Edit | Attach | Watch | Print version | History: r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r3 - 2006-09-14 - AndreHolzner
 
    • 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-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