#-------------------------------------------------------------------------
# Algorithms
#-------------------------------------------------------------------------
from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()
# Add top algorithms to be run
from SCT2.SCT2Conf import MySCTtscan
job += MySCTtscan("MySCTtscan")   
#--------------------------------------------------------------
# Set output level threshold (DEBUG, INFO, WARNING, ERROR, FATAL)
#--------------------------------------------------------------
job.MySCTtscan.OutputLevel = WARNING
#job.AtRndmGenSvc.OutputLevel = WARNING
#job.AthenaEventLoopMgr.OutputLevel = WARNING
#--------------------------------------------------------------
# Event related parameters
#--------------------------------------------------------------
theApp.EvtMax = 5000
#------------------------------------------------------------------------
#  The properties of the MySCTtscan algorithm
#------------------------------------------------------------------------
job.MySCTtscan.EfieldModel = 2
#job.MySCTtscan.IsSCTDigiModel = FALSE
job.MySCTtscan.IsSCTDigiModel = TRUE
job.MySCTtscan.TimeOffset = 30.0
job.MySCTtscan.TimeOfThreshold_min = -25.0
job.MySCTtscan.TimeOfThreshold_max = 75.0
job.MySCTtscan.IsNewPulseShape = FALSE
#job.MySCTtscan.IsNewPulseShape = TRUE
job.MySCTtscan.Noise = 1600.
job.MySCTtscan.IsLevelSensing = FALSE
job.MySCTtscan.Felectron = 0.00
job.MySCTtscan.Fpion = 1.00
job.MySCTtscan.Fkaon = 0.000
job.MySCTtscan.Fproton = 0.00
job.MySCTtscan.StripNo_min = 0
job.MySCTtscan.StripNo_max = 767
job.MySCTtscan.BarrelLayer = 0  
job.MySCTtscan.LocPhi_min = -90.0
job.MySCTtscan.LocPhi_max = 90.0
job.MySCTtscan.Polarity =  0
job.MySCTtscan.Pt_min = 0.1
job.MySCTtscan.Pt_max = 10.0 
job.MySCTtscan.PtPower = -2.0 
job.MySCTtscan.BarrelLayer = 0  ;
job.MySCTtscan.Zmodule_min = 61.8
job.MySCTtscan.Zmodule_max = 74.2
job.MySCTtscan.DepletionVoltage_VD = 70.
job.MySCTtscan.BiasVoltage_VB = 150.
job.MySCTtscan.MagneticField_B = -2.0
job.MySCTtscan.Temperature_T = 273.
job.MySCTtscan.X_origin_min = 0.0000
job.MySCTtscan.X_origin_max = 0.0080
job.MySCTtscan.StripTilt = 1
job.MySCTtscan.IsdEdX = TRUE
job.MySCTtscan.ChargeTimeStep = 0.5 
job.MySCTtscan.BiasVoltage_VB = 150.
job.MySCTtscan.MagneticField_B = -2.0
job.MySCTtscan.Temperature_T = 273.15
job.MySCTtscan.TransportTimeStep = 0.25
job.MySCTtscan.TransportTimeMax = 25.0
job.MySCTtscan.Threshold = 1.0
job.MySCTtscan.CrossTalk = 1.0
job.MySCTtscan.Ydivision = 500
job.MySCTtscan.eh_pairs = 108.0 
job.MySCTtscan.IsLandau = TRUE
job.MySCTtscan.CoutLevel= 0  # =3(eh), 2(strip), 1(event), 0(no cout)
# save ROOT histograms and Tuple
#------------------------------------------------------------
from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
ServiceMgr.THistSvc.Output = [ "AANT DATAFILE='tscan_A_edge.root' OPT='RECREATE'" ]
#==============================================================
# End of job options file