def ArgoMatch():
palgo_array = cmtool.CPAlgoArray()
palgo1 = cmtool.CPAlgoNHits()
# palgo1.SetMinHits(45)
palgo_array.AddAlgo(palgo1)
algo_array = cmtool.CFAlgoArray()
# wireAlg = argomatch.CFAlgoShowerWireMatch()
timeAlg = cmtool.CFAlgoTimeOverlap()
timeAlg.RequireThreePlanes(False)
# algo_array.AddAlgo(wireAlg)
algo_array.AddAlgo(timeAlg)
return palgo_array, algo_array
def getMatchManager():
nplanes = larutil.Geometry.GetME().Nplanes()
manager = cmtool.CMatchManager(nplanes)
priority_array = cmtool.CPAlgoArray()
palgo1 = cmtool.CPAlgoNHits()
palgo1.SetMinHits(20)
priority_array.AddAlgo(palgo1)
match_array = cmtool.CFAlgoArray()
timeAlg = cmtool.CFAlgoTimeOverlap()
timeAlg.RequireThreePlanes(False)
match_array.AddAlgo(timeAlg)
manager.AddPriorityAlgo(priority_array)
manager.AddMatchAlgo(match_array)
return manager
ana_unit = fmwk.ShowerReco3D() # Attach shower reco alg sralg = showerreco.ShowerRecoTest() sralg.Verbose(False) # sralg.Verbose(True) sralg.SetUseArea(True) #sralg.SetUseModBox(True) ana_unit.SetShowerAlgo(sralg) # # Attach Matching algorithm # palgo_array = cmtool.CPAlgoArray() palgo1 = cmtool.CPAlgoNHits() palgo1.SetMinHits(25) palgo2 = cmtool.CPAlgoIgnoreTracks() palgo_array.AddAlgo(palgo1) palgo_array.AddAlgo(palgo2) ana_unit.GetManager().AddPriorityAlgo(palgo_array) #Andrzej: The algorithms below are ranked by their effectiveness-- TimeOverlap is best, #then 3DAngle, then StartPoint . Right now, only TimeOverlap is called. algo_array = cmtool.CFAlgoArray() #algo_array.SetMode(cmtool.CFAlgoArray.kPositiveAddition) algo_array.AddAlgo(cmtool.CFAlgoTimeOverlap()) #algo_array.AddAlgo(cmtool.CFAlgoTimeProf()) #algo_array.AddAlgo(cmtool.CFAlgo3DAngle())
# Specify analysis output root file name
my_proc.set_ana_output_file("from_test_ana_you_can_remove_me.root")
# Specify data output root file name
my_proc.set_output_file("out.root")
# Create analysis unit
ana_unit = fmwk.ShowerReco3D()
# Attach shower reco alg
ana_unit.SetShowerAlgo(showerreco.ShowerRecoAlg())
# Specify cluster type
ana_unit.SetClusterProducer("mergedfuzzy")
priority_algo = cmtool.CPAlgoNHits()
priority_algo.SetMinHits(20)
ana_unit.GetManager().AddPriorityAlgo(priority_algo)
#
# Attach Matching algorithm
#
#Andrzej: The algorithms below are ranked by their effectiveness-- TimeOverlap is best,
#then 3DAngle, then StartPoint . Right now, only TimeOverlap is called.
algo_array = cmtool.CFAlgoArray()
#algo_array.SetMode(cmtool.CFAlgoArray.kPositiveAddition)
algo_array.AddAlgo(cmtool.CFAlgoTimeOverlap())
#algo_array.AddAlgo(cmtool.CFAlgo3DAngle())
#algo_array.AddAlgo(cmtool.CFAlgoStartPointMatch())