set stream       $env(STREAM)
set mergefile    /data/harv11/b/lmiller/EventMerge/data/${stream}r0256b7.0001phys
set outputfile   l00spy_10k_$stream.root
set run 153271

# input module for generators
module input GenInputManager

talk GenInputManager
  run_number set $run
exit

# enable generator ( no generators are enabled bu default )
mod enable FAKE_EVENT
talk FAKE_EVENT
  use PT
  use ETA
  use PHI
  RandomNumberMenu
   RandomSeed1 set $env(RANDOM1)
   RandomSeed2 set $env(RANDOM2)
  exit
  generate PT  10. 0. 10. 10. 0. 1.

# Generate mu-
  generate CDFCODE 208
  generate NPARTICLES 1
#  show
exit

talk GeometryManager
# Use the defaults but substitute the detailed for the simple Si detector
# The following set up is to be added to the default GeometryManager set up.
   DetectorMenu
      enableAll set f
      enableCot set t
      enableSvx set t
# central part of the beampipe
      enableBeampipeC set true   
      show
   exit
# Silicon Passive Geometry
   SiliconGeometryMenu
      BuildPassive set true
   exit
exit

talk SimInitManager
  DetectorMenu
     declareSvx set t 
     declareCot set t
#     declarePassive set t
     declareBeampipeC set t
  exit
exit

talk SimulationControlMod
  DebugMenu
    showActiveVolumes set t
  exit
  DetectorMenu
     simulateSvx set t
     simulateCot set t
#     simulateMuon set t
#     simulateTof set t
#     simulateCalor set t
#     simulatePassive set t
  exit
  dumpFactory
  add CdfHalfLadder SvxDigitizer SvxGroup
  add CotSuperLayer CotDigitizer CotGroup
  show add
  ConfigMenus
    CotGroup_CotSuperLayer
      DriftModel set Garfield
    exit
    SvxGroup_CdfHalfLadder
       pick_svx_cdm set PHYSICAL
       svx_cdm_magnetic_field set OFF
       svx_cdm_noise set OFF
#       pick_svx_cdm set PARAMETERIZED
       ZigZagL7 set t
       CreatePropagatedSi set t
       svx_validation set ON
       CreateSIXD set t
     exit
  exit
exit

module talk EventMerger
# This secondary data will be put into the primary
  Filename set $mergefile
# random seed used to add poisson number of events from secondary
  RandomSeed set 3980343
  PoissonFluct set false
  meanNumAdded set 1 
# extra SIXD bank descriptors and
# append secondary dataset...
  WriteSecondSiBanks set true
  WriteMergeSiBanks set true
  SecondarySiBanksDescr set Secondary
  MergeSiBanksDescr set Merged
# To make sure if merging is working, change writeDebug set t
  useNNDataUnpacker set t
  useNNChipThreshold set t
  writeDebug set f
exit

# correct the strips of the merged set
mod talk SiClusteringModule
  StripCorrectorMenu
    RemoveUnderThreshold set false
  exit
  UseNNDataUnpacker set t
  MonitorMode set t
  WriteStripBanks set t
  InputBanksDescr set "Merged"
  WriteStripBanks set t
  StripCorrectorMenu
    PedDBRun set 153214
    PedDBVersion set 4
    PedDBStatus set COMPLETE
    DeadDBRun set 152748
    DeadDBVersion set 1
    DeadDBStatus set COMPLETE
    RemoveUnderThreshold set f
    FlagBad set t
    RemoveBad set f
    SubtractPedestal set t
    DcmsMode set f
    DcmsOffset set 0
#    L00FitOption set 3
    UseL00PedFit set t
    L00FitOrder set 6
  exit
exit

mod talk GeometryManager
  DetectorMenu
   enableMuon set f
  exit
  SiliconGeometryMenu
#    AlignmentSource set  "ofotl_prd_read    100014    2    TEST"
#    AlignmentSource set  "ofotl_prd_read    100014    2    TEST"
  exit
exit

module talk CalibrationManager
   add ofotl_prd OTL cdf_reader/reader@cdfofprd
   DataDB      set ofotl_prd
   Database    set ofotl_prd
exit

mod talk L00Spy
 # descriptors of all the input SIXD banks
 # there is no control over corrected bank descriptors 
 # so you will need to drop any other corrected SIXD banks
 PrimaryBanksDescr set ""
 SecondaryBanksDescr set Secondary
 MergeBanksDescr set Merged
 CorrectedBanksDescr set CORRECTED
 ClusterBanksDescr SiHitSet
 mydebug set f
exit

if [file exists $outputfile] { exec rm $outputfile }
#mod enable HepHbookManager
#mod talk HepHbookManager 
mod enable HepRootManager
mod talk HepRootManager 
  histfile set $outputfile 
exit

#mod enable EventDump
mod disable EventDump
mod talk EventDump
summary set t
#  classList set SIXD_StorableBank ISLD_StorableBank 
#SiHitSet
exit

mod disable CT_TrackingModule
mod disable SiPatternRecModule
path enable AllPath

path list
ev begin -nev 10000
show timer

exit