def getShowerRecoAlgModular(): # This function returns the default shower reco module # If you want to extend, customize, or otherwise alter the default # reco (which you should!) it's recommended that you do one of two # things: # 1) Copy this function, and change the lines you need # 2) Use this function as is and replace the modules with the built-in functions # to ShowerRecoAlgModular # Then, use the return value (alg) in place of ShowerRecoAlg in your python script # It should be noted, however, that this "default" set of modules is not yet developed # and if you are developing you ought to be updating it here! alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(False) # 3D Axis Module: axis3D = showerreco.Axis3DModule() axis3D.setMaxIterations(100) axis3D.setNStepsInitial(25) axis3D.setTargetError(0.001) axis3D.setNormalErrorRange(0.01) axis3D.setThetaRangeStart(0.1) axis3D.setThetaRangeMin(0.0005) axis3D.setNStepsStart(4) axis3D.setConvergeRate(0.85) axis3D.setVerbosity(True) axis3D.setSeedVectorErrorCutoff(0.1) angle3D = showerreco.Angle3DFormula() angle3D.setMaxAngleError(0.1) angle3D.setVerbosity(False) energy = showerreco.LinearEnergy() energy.SetUseModBox(True) energy.setVerbosity(False) dqdx = showerreco.dQdxModule() dedx = showerreco.dEdxFromdQdx() dedx.SetUsePitch(False) dedx.setVerbosity(False) #alg.AddShowerRecoModule(axis3D) alg.AddShowerRecoModule(angle3D) alg.AddShowerRecoModule(showerreco.StartPoint3DModule()) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule(dqdx) alg.AddShowerRecoModule(dedx) # alg.AddShowerRecoModule(showerreco.StartPoint2DModule() ) #alg.AddShowerRecoModule(showerreco.OtherStartPoint3D() ) # alg.AddShowerRecoModule(showerreco.ShowerChargeModule() ) alg.AddShowerRecoModule(showerreco.GeoModule()) alg.PrintModuleList() return alg
def getShowerRecoAlgModular(): # This function returns the default shower reco module # If you want to extend, customize, or otherwise alter the default # reco (which you should!) it's recommended that you do one of two # things: # 1) Copy this function, and change the lines you need # 2) Use this function as is and replace the modules with the built-in functions # to ShowerRecoAlgModular # Then, use the return value (alg) in place of ShowerRecoAlg in your python script # It should be noted, however, that this "default" set of modules is not yet developed # and if you are developing you ought to be updating it here! alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(False) # 3D Axis Module: axis3D = showerreco.Axis3DModule() axis3D.setMaxIterations(100) axis3D.setNStepsInitial(25) axis3D.setTargetError(0.001) axis3D.setNormalErrorRange(0.01) axis3D.setThetaRangeStart(0.1) axis3D.setThetaRangeMin(0.0005) axis3D.setNStepsStart(4) axis3D.setConvergeRate(0.85) axis3D.setVerbosity(True) axis3D.setSeedVectorErrorCutoff(0.1) angle3D = showerreco.Angle3DFormula() angle3D.setMaxAngleError(0.1) angle3D.setVerbosity(False) energy = showerreco.LinearEnergy() energy.SetUseModBox(True) energy.setVerbosity(False) dqdx = showerreco.dQdxModule() dedx = showerreco.dEdxFromdQdx() dedx.SetUsePitch(False) dedx.setVerbosity(False) #alg.AddShowerRecoModule(axis3D) alg.AddShowerRecoModule(angle3D) alg.AddShowerRecoModule(showerreco.StartPoint3DModule() ) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule(dqdx) alg.AddShowerRecoModule(dedx) # alg.AddShowerRecoModule(showerreco.StartPoint2DModule() ) #alg.AddShowerRecoModule(showerreco.OtherStartPoint3D() ) # alg.AddShowerRecoModule(showerreco.ShowerChargeModule() ) alg.AddShowerRecoModule(showerreco.GeoModule()) alg.PrintModuleList() return alg
def getShowerRecoAlgModular(): # This function returns the default shower reco module # If you want to extend, customize, or otherwise alter the default # reco (which you should!) it's recommended that you do one of two # things: # 1) Copy this function, and change the lines you need # 2) Use this function as is and replace the modules with the built-in functions # to ShowerRecoAlgModular # Then, use the return value (alg) in place of ShowerRecoAlg in your python script # It should be noted, however, that this "default" set of modules is not yet developed # and if you are developing you ought to be updating it here! alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(True) alg.SetVerbose(False) angle3D = showerreco.Angle3DFormula() angle3D.setMaxAngleError(0.1) angle3D.setValidateDirection(True) angle3D.setVerbosity(False) angle3D = showerreco.Angle3DFromVtxQweighted() energy = showerreco.LinearEnergy() energy.SetElectronLifetime(8000.) # in us MCC7 value #energy.SetElectronLifetime(1e6 ) # in us DATA value energy.SetRecombFactor(0.62) #energy.SetElecGain(198) # e- / ADC DATA value energy.SetElecGain(189) # MCC7 value energy.setVerbosity(False) energy.SetFillTree(True) dqdx = showerreco.dQdxModule() dedx = showerreco.dEdxFromdQdx() dedx.SetUsePitch(False) dedx.setVerbosity(True) geoModule = showerreco.GeoModule() geoModule.setFlipShowerDirection(False) alg.AddShowerRecoModule( showerreco.StartPoint3DModule() ) alg.AddShowerRecoModule(angle3D) alg.AddShowerRecoModule(energy) #alg.AddShowerRecoModule(dqdx) #alg.AddShowerRecoModule(dedx) alg.AddShowerRecoModule( geoModule ) alg.PrintModuleList() return alg
def ShowerReco(): ana_unit = fmwk.ShowerReco3D() alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(False) axis3D = showerreco.Axis3DModule() axis3D.setMaxIterations(100) axis3D.setNStepsInitial(25) axis3D.setTargetError(0.001) axis3D.setNormalErrorRange(0.01) axis3D.setThetaRangeStart(0.1) axis3D.setThetaRangeMin(0.0005) axis3D.setNStepsStart(4) axis3D.setConvergeRate(0.85) axis3D.setVerbosity(False) axis3D.setSeedVectorErrorCutoff(0.1) angle3D = showerreco.Angle3DFormula() angle3D.setMaxAngleError(0.1) angle3D.setVerbosity(False) energy = showerreco.LinearEnergy() energy.SetUseModBox(True) energy.setVerbosity(False) dqdx = showerreco.dQdxModule() dedx = showerreco.dEdxFromdQdx() dedx.SetUsePitch(False) dedx.setVerbosity(False) #alg.AddShowerRecoModule(showerreco.StartPoint2DModule() ) #alg.AddShowerRecoModule(showerreco.SanityCheckInitial() ) alg.AddShowerRecoModule(axis3D) alg.AddShowerRecoModule(showerreco.StartPoint3DModule() ) #alg.AddShowerRecoModule(angle3D) #alg.AddShowerRecoModule(showerreco.StartPoint3DModule() ) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule(dqdx) alg.AddShowerRecoModule(dedx) #alg.AddShowerRecoModule(showerreco.OtherStartPoint3D() ) #alg.AddShowerRecoModule(showerreco.ShowerChargeModule() ) #alg.AddShowerRecoModule(showerreco.GeoModule()) #alg.AddShowerRecoModule(showerreco.SanityCheckFinal() ) alg.PrintModuleList() ana_unit.AddShowerAlgo(alg) return ana_unit
def getShowerRecoAlgModular(): # This function returns the default shower reco module # If you want to extend, customize, or otherwise alter the default # reco (which you should!) it's recommended that you do one of two # things: # 1) Copy this function, and change the lines you need # 2) Use this function as is and replace the modules with the built-in functions # to ShowerRecoAlgModular # Then, use the return value (alg) in place of ShowerRecoAlg in your python script # It should be noted, however, that this "default" set of modules is not yet developed # and if you are developing you ought to be updating it here! alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(False) startPoint = showerreco.StartPointPandora() axis = showerreco.PCA3DAxis() energy = showerreco.LinearEnergy() #energy.SetHitRecoCorrection(1.289, 0) #energy.SetHitRecoCorrection(1.323, 1) #energy.SetHitRecoCorrection(1.336, 2) energy.SetUseModBox(True) energy.setVerbosity(False) dqdx = showerreco.dQdxModule() dedx = showerreco.dEdxFromdQdx() dedx.SetUsePitch(False) dedx.setVerbosity(False) alg.AddShowerRecoModule(showerreco.ToyGeoModule()) alg.AddShowerRecoModule(axis) alg.AddShowerRecoModule(startPoint) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule(dqdx) alg.AddShowerRecoModule(dedx) # alg.AddShowerRecoModule(showerreco.StartPoint2DModule() ) #alg.AddShowerRecoModule(showerreco.OtherStartPoint3D() ) # alg.AddShowerRecoModule(showerreco.ShowerChargeModule() ) # alg.AddShowerRecoModule(showerreco.GeoModule()) alg.PrintModuleList() return alg
def getShowerRecoAlgModular(): # This function returns the default shower reco module # If you want to extend, customize, or otherwise alter the default # reco (which you should!) it's recommended that you do one of two # things: # 1) Copy this function, and change the lines you need # 2) Use this function as is and replace the modules with the built-in functions # to ShowerRecoAlgModular # Then, use the return value (alg) in place of ShowerRecoAlg in your python script # It should be noted, however, that this "default" set of modules is not yet developed # and if you are developing you ought to be updating it here! alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(False) startPoint = showerreco.StartPointPandora() axis = showerreco.PCA3DAxis() energy = showerreco.LinearEnergy() #energy.SetHitRecoCorrection(1.289, 0) #energy.SetHitRecoCorrection(1.323, 1) #energy.SetHitRecoCorrection(1.336, 2) energy.SetUseModBox(True) energy.setVerbosity(False) dqdx = showerreco.dQdxModule() dedx = showerreco.dEdxFromdQdx() dedx.SetUsePitch(False) dedx.setVerbosity(False) alg.AddShowerRecoModule(showerreco.ToyGeoModule()) alg.AddShowerRecoModule(axis) alg.AddShowerRecoModule(startPoint) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule(dqdx) alg.AddShowerRecoModule(dedx) # alg.AddShowerRecoModule(showerreco.StartPoint2DModule() ) #alg.AddShowerRecoModule(showerreco.OtherStartPoint3D() ) # alg.AddShowerRecoModule(showerreco.ShowerChargeModule() ) # alg.AddShowerRecoModule(showerreco.GeoModule()) alg.PrintModuleList() return alg
def getShowerRecoAlgModular(): # This function returns the default shower reco module # If you want to extend, customize, or otherwise alter the default # reco (which you should!) it's recommended that you do one of two # things: # 1) Copy this function, and change the lines you need # 2) Use this function as is and replace the modules with the built-in functions # to ShowerRecoAlgModular # Then, use the return value (alg) in place of ShowerRecoAlg in your python script # It should be noted, however, that this "default" set of modules is not yet developed # and if you are developing you ought to be updating it here! alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(False) angle3D = showerreco.Angle3DFromVtxQweighted() energy = showerreco.LinearEnergy() energy.SetUseModBox(True) energy.SetUseArea(True) energy.setVerbosity(False) energy.SetFillTree(True) energy.SetMinNumHits(20) dqdx = showerreco.dQdxModule() dqdx.setVerbosity(False) dedx = showerreco.dEdxFromdQdx() dedx.setVerbosity(False) alg.AddShowerRecoModule(showerreco.StartPoint3DModule() ) alg.AddShowerRecoModule(angle3D) filter = showerreco.Filter3DDir() filter.SetDotProdCut(0.5) alg.AddShowerRecoModule(filter) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule(dqdx) alg.AddShowerRecoModule(dedx) alg.AddShowerRecoModule(showerreco.GeoModule()) alg.PrintModuleList() return alg
def getShowerRecoAlgModular(): alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(False) alg.SetVerbose(False) angle3D = showerreco.Angle3DFormula() angle3D.setMaxAngleError(0.1) angle3D.setValidateDirection(True) angle3D.setVerbosity(False) angle3D = showerreco.Angle3DFromVtxQweighted() gain = 196.97874008458268 energy = showerreco.LinearEnergy() energy.SetElectronLifetime(1e20) #1e6) #8000.) # in us MCC7 value energy.SetRecombFactor(0.572) energy.SetElecGain(gain) # e- / ADC DATA value, MCC8 #energy.SetElecGain(189) # MCC7 value energy.setVerbosity(False) energy.SetFillTree(True) dqdx = showerreco.dQdxModule() dedx = showerreco.dEdxFromdQdx() dedx.SetUsePitch(False) dedx.setVerbosity(True) geoModule = showerreco.GeoModule() geoModule.setFlipShowerDirection(False) alg.AddShowerRecoModule( showerreco.StartPoint3DModule() ) alg.AddShowerRecoModule(angle3D) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule( geoModule ) alg.PrintModuleList() return alg
def getShowerRecoAlgModular(): # This function returns the default shower reco module # If you want to extend, customize, or otherwise alter the default # reco (which you should!) it's recommended that you do one of two # things: # 1) Copy this function, and change the lines you need # 2) Use this function as is and replace the modules with the built-in functions # to ShowerRecoAlgModular # Then, use the return value (alg) in place of ShowerRecoAlg in your python script # It should be noted, however, that this "default" set of modules is not yet developed # and if you are developing you ought to be updating it here! alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(True) alg.SetVerbose(False) filteralgo = showerreco.FilterPFPart() filteralgo.setMinNHitsAbsolute(5) filteralgo.setMinNHitsLargest(10) angle3D = showerreco.Angle3DFormula() angle3D.setMaxAngleError(0.1) angle3D.setValidateDirection(True) angle3D.setVerbosity(False) angle3D = showerreco.Angle3DFromVtxQweighted() startPt = showerreco.YPlaneStartPoint3D() startPt.setVerbosity(False) energy = showerreco.LinearEnergy() energy.SetElectronLifetime(1e6 ) # in us DATA value energy.SetRecombFactor(0.62) energy.SetElecGain(243.) # MCC8.0 data #energy.SetElecGain(200.) # MCC8.0 value energy.setVerbosity(False) energy.SetFillTree(True) dqdx = showerreco.dQdxModule() dedx = showerreco.dEdxFromdQdx() dedx.SetUsePitch(False) dedx.setVerbosity(True) shrFilter = showerreco.FilterShowers() shrFilter.setAngleCut(15.) shrFilter.setVerbosity(True) alg.AddShowerRecoModule( filteralgo ) alg.AddShowerRecoModule(angle3D) #alg.AddShowerRecoModule( showerreco.StartPoint3DModule() ) #alg.AddShowerRecoModule( showerreco.NearestStartPoint3D() ) alg.AddShowerRecoModule( startPt )#showerreco.YPlaneStartPoint3D() ) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule( showerreco.FillLength() ) alg.AddShowerRecoModule( shrFilter ) alg.PrintModuleList() return alg
def getShowerRecoAlgModular(): # This function returns the default shower reco module # If you want to extend, customize, or otherwise alter the default # reco (which you should!) it's recommended that you do one of two # things: # 1) Copy this function, and change the lines you need # 2) Use this function as is and replace the modules with the built-in functions # to ShowerRecoAlgModular # Then, use the return value (alg) in place of ShowerRecoAlg in your python script # It should be noted, however, that this "default" set of modules is not yet developed # and if you are developing you ought to be updating it here! alg = showerreco.ShowerRecoAlgModular() alg.SetDebug(True) alg.SetVerbose(False) filteralgo = showerreco.FilterPFPart() filteralgo.setMinNHitsAbsolute(5) filteralgo.setMinNHitsLargest(10) angle3D = showerreco.Angle3DFormula() angle3D.setMaxAngleError(0.1) angle3D.setValidateDirection(True) angle3D.setVerbosity(False) angle3D = showerreco.Angle3DFromVtxQweighted() startPt = showerreco.YPlaneStartPoint3D() startPt.setVerbosity(False) energy = showerreco.LinearEnergy() # implement position-dependent calibration energy.CreateResponseMap(15) dQdsAVG = 0. ctr = 0. fin = open('/home/david/Dropbox/Neutrinos/Calorimetry/3DCorrection/MCC82_0903/dqds_0903_data_xyz.txt','r') for line in fin: words = line.split() q = float(words[3]) dQdsAVG += q ctr += 1 dQdsAVG /= ctr print 'dQdsAVG = %.02f'%dQdsAVG fin.seek(0) for line in fin: words = line.split() x = float(words[0]) y = float(words[1]) z = float(words[2]) q = float(words[3]) energy.SetResponse(x,y,z,dQdsAVG/q) energy.SetElectronLifetime(1e6) # in us DATA value energy.SetRecombFactor(0.572) energy.SetElecGain(243.) # MCC8.0 data #energy.SetElecGain(200.) # MCC8.0 value energy.setVerbosity(False) energy.SetFillTree(True) dqdx = showerreco.dQdxModule() dqdx.setTrunkLength(2.) #dqdx.SetFillTree(True) shrFilter = showerreco.FilterShowers() shrFilter.setAngleCut(25.) shrFilter.setVerbosity(False) alg.AddShowerRecoModule( filteralgo ) alg.AddShowerRecoModule(angle3D) #alg.AddShowerRecoModule( showerreco.StartPoint3DModule() ) #alg.AddShowerRecoModule( showerreco.NearestStartPoint3D() ) alg.AddShowerRecoModule( startPt )#showerreco.YPlaneStartPoint3D() ) alg.AddShowerRecoModule( dqdx ) alg.AddShowerRecoModule(energy) alg.AddShowerRecoModule( showerreco.FillLength() ) alg.AddShowerRecoModule( shrFilter ) alg.PrintModuleList() return alg