mergelayersHcalFwd.inhits.Path = "HCalFwdHits"
mergelayersHcalFwd.outhits.Path = "mergedHCalFwdHits"

createHcalFwdCells = CreateCaloCells("CreateHcalFwdCaloCells",
                                     doCellCalibration=True,
                                     calibTool=calibHcalFwd,
                                     addCellNoise=False,
                                     filterCellNoise=False,
                                     OutputLevel=INFO)
createHcalFwdCells.hits.Path = "mergedHCalFwdHits"
createHcalFwdCells.cells.Path = "HCalFwdCells"

out = PodioOutput("out", OutputLevel=INFO)
out.outputCommands = [
    "drop *", "keep ECalBarrelCells", "keep ECalEndcapCells",
    "keep ECalFwdCells", "keep HCalBarrelCells", "keep HCalExtBarrelCells",
    "keep HCalEndcapCells", "keep HCalFwdCells", "keep GenParticles",
    "keep GenVertices"
]
out.filename = "output_fullCalo_SimAndDigi_e50GeV_" + str(
    num_events) + "events.root"

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
#Fill a collection of CaloHitPositions for detailed Cluster analysis
from Configurables import CreateCaloCellPositions 
positionsClusterBarrel =CreateCaloCellPositions("positionsClusterBarrel",
                                                positionsECalBarrelTool = ECalBcells,
                                                positionsHCalBarrelTool = HCalBcells,
                                                positionsHCalExtBarrelTool = HCalExtBcells,
                                                positionsEMECTool = EMECcells,
                                                positionsHECTool = HECcells,
                                                positionsEMFwdTool = ECalFwdcells,
                                                positionsHFwdTool = HCalFwdcells,
                                                hits = "caloClusterBarrelCells",
                                                positionedHits = "caloClusterBarrelCellPositions",
                                                OutputLevel = INFO)

out = PodioOutput("out", filename = "output_BarrelTopo_electrNoise_50GeVe_3ev.root", OutputLevel = DEBUG)
out.outputCommands =["drop *", "keep GenParticles", "keep GenVertices", "keep caloClustersBarrel","keep caloClusterBarrelCells", "keep caloClusterBarrelCellPositions"]

#CPU information
from Configurables import AuditorSvc,ChronoAuditor 
chra = ChronoAuditor() 
audsvc = AuditorSvc() 
audsvc.Auditors =[chra] 
podioinput.AuditExecute = True 
createemptycells.AuditExecute =True 
createEcalBarrelCells.AuditExecute = True 
createHcalBarrelCells.AuditExecute =True 
createTopoClusters.AuditExecute = True 
positionsClusterBarrel.AuditExecute = True

ApplicationMgr(TopAlg =
               [podioinput, 
예제 #3
0
                                          positionedHits = "HCalEndcapCellPositions", 
                                          OutputLevel = INFO)
positionsEcalFwd = CreateCellPositions("positionsEcalFwd", 
                                          positionsTool=ECalFwdcells, 
                                          hits = "ECalFwdCells", 
                                          positionedHits = "ECalFwdCellPositions", 
                                          OutputLevel = INFO)
positionsHcalFwd = CreateCellPositions("positionsHcalFwd", 
                                          positionsTool=HCalFwdcells, 
                                          hits = "HCalFwdCells", 
                                          positionedHits = "HCalFwdCellPositions", 
                                          OutputLevel = INFO)

out = PodioOutput("out", OutputLevel=DEBUG)
out.filename = "digi_cellPositions_50GeVelectrons.root"
out.outputCommands = ["keep *","drop ECalBarrelCells","drop ECalEndcapCells","drop ECalFwdCells","drop HCalBarrelCells", "drop HCalExtBarrelCells", "drop HCalEndcapCells", "drop HCalFwdCells"]

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
positionsEcalEndcap.AuditExecute = True
positionsEcalFwd.AuditExecute = True
positionsHcalBarrel.AuditExecute = True
positionsHcalExtBarrel.AuditExecute = True
positionsHcalEndcap.AuditExecute = True
positionsHcalFwd.AuditExecute = True
out.AuditExecute = True
예제 #4
0
savecaltool.DataOutputs.caloHits.Path = "caloHits"
# next, create the G4 algorithm, giving the list of names of tools ("XX/YY")
particle_converter = SimG4PrimariesFromEdmTool("EdmConverter")
particle_converter.DataInputs.genParticles.Path = "allGenParticles"
geantsim = SimG4Alg("SimG4Alg",
                    outputs = ["SimG4SaveSmearedParticles/saveSmearedParticles", "SimG4SaveCalHits/saveCalHits"],
                    eventProvider=particle_converter)

from Configurables import SimG4FastSimHistograms
hist = SimG4FastSimHistograms("fastHist")
hist.DataInputs.particlesMCparticles.Path = "particleMCparticleAssociation"
THistSvc().Output = ["rec DATAFILE='histFormula.root' TYP='ROOT' OPT='RECREATE'"]
THistSvc().PrintAll=True
THistSvc().AutoSave=True
THistSvc().AutoFlush=True
THistSvc().OutputLevel=INFO

from Configurables import PodioOutput
## PODIO algorithm
out = PodioOutput("out", filename = "out_fast_tracker_formula_calo_gflash.root")
out.outputCommands = ["keep *"]

# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr( TopAlg = [reader, hepmc_converter, geantsim, hist, out],
                EvtSel = 'NONE',
                EvtMax   = 1,
                # order is important, as GeoSvc is needed by SimG4Svc
                ExtSvc = [podioevent, geoservice, geantservice],
                OutputLevel=INFO)
예제 #5
0
createClusters.nEtaDuplicates =  7
createClusters.nPhiDuplicates = 13
createClusters.nEtaFinal =       9
createClusters.nPhiFinal =      17
# Minimal energy to create a cluster in GeV (FCC-ee detectors have to reconstruct low energy particles)
createClusters.energyThreshold = 0.1
createClusters.attachCells = True
createClusters.clusters.Path = "CaloClusters"
createClusters.clusterCells.Path = "CaloClusterCells"
ApplicationMgr().TopAlg += [createClusters]


################ Output
from Configurables import PodioOutput
out = PodioOutput("out")
out.outputCommands = ["keep *", "drop ECalBarrelHits", "drop HCal*", "drop ECalBarrelCellsStep*", "drop ECalBarrelPositionedHits", "drop emptyCaloCells", "drop CaloClusterCells"]
import uuid
out.filename = "output_fullCalo_SimAndDigi.root"
ApplicationMgr().TopAlg += [out]

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
gen.AuditExecute = True
hepmc_converter.AuditExecute = True
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
hepmc_converter.genparticles.Path="genParticles"
hepmc_converter.genvertices.Path="genVertices"

## FCC event-data model output -> define objects to be written out
from Configurables import PodioOutput

out = PodioOutput("out",OutputLevel=messageLevelOut)
out.filename       = "FCCDelphesOutput.root"
if args.outputfile != '':
    out.filename = args.outputfile

#out.outputCommands = ["drop *",
#                      "keep genParticles",
#                      "keep genVertices",
#                      "keep genJets",
#                      "keep genJetsToMC"]
out.outputCommands = ["keep *", "drop genParticles", "drop genVertices"]

############################################################
#
# Run modules
#
############################################################
# Run Pythia + Delphes
ApplicationMgr( TopAlg = [ pythia8gen, hepmc_converter, delphessim, out ],
                EvtSel = 'NONE',
                EvtMax = nEvents,
                ExtSvc = [podioEvent])

ApplicationMgr.EvtMax=10
genfilter.allGenParticles.Path = "allGenParticles"
genfilter.filteredGenParticles.Path = "GenParticles"
from Configurables import SimG4PrimariesFromEdmTool
particle_converter = SimG4PrimariesFromEdmTool("EdmConverter")
particle_converter.genParticles.Path = "GenParticles"
geantsim.eventProvider = particle_converter

# PODIO algorithm
from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput
podioevent = FCCDataSvc("EventDataSvc")
out = PodioOutput("out")
out.outputCommands = [
    "drop *",
    "keep GenParticles",
    "keep GenVertices",
    "keep SimParticles",
    "keep SimVertices",
    "keep TrackerHits",
    "keep TrackerPositionedHits",
    "keep TrackerDigiPostPoint",
]
out.filename = output_name

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
out.AuditExecute = True

list_of_algorithms = [genalg, hepmc_converter, genfilter, geantsim, out]
                   merge = hcalFwdNumberOfLayersToMerge,
                   OutputLevel = INFO)
mergelayersHcalFwd.inhits.Path = "HCalFwdHits"
mergelayersHcalFwd.outhits.Path = "mergedHCalFwdHits"

createHcalFwdCells = CreateCaloCells("CreateHcalFwdCaloCells",
                                 doCellCalibration=True,
                                 calibTool=calibHcalFwd,
                                 addCellNoise=False, filterCellNoise=False,
                                 OutputLevel=INFO)
createHcalFwdCells.hits.Path="mergedHCalFwdHits"
createHcalFwdCells.cells.Path="HCalFwdCells"

out = PodioOutput("out",
                  OutputLevel=INFO)
out.outputCommands = ["drop *", "keep ECalBarrelCells", "keep ECalEndcapCells", "keep ECalFwdCells", "keep HCalBarrelCells", "keep HCalBarrelPositions", "keep HCalExtBarrelCells", "keep HCalEndcapCells", "keep HCalFwdCells", "keep GenParticles","keep GenVertices"]
out.filename = "output_fullCalo_SimAndDigi_pi50GeV_"+str(num_events)+"events.root"

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
createEcalEndcapCells.AuditExecute = True
createEcalFwdCells.AuditExecute = True
createHcalCells.AuditExecute = True
예제 #9
0
hepmc_converter.genparticles.Path="genParticles"
hepmc_converter.genvertices.Path="genVertices"

## FCC event-data model output -> define objects to be written out
from Configurables import PodioOutput

out = PodioOutput("out",OutputLevel=messageLevelOut)
out.filename       = "FCCDelphesOutput.root"
if args.outputfile != '':
    out.filename = args.outputfile

#out.outputCommands = ["drop *",
#                      "keep genParticles",
#                      "keep genVertices",
#                      "keep genJets",
#                      "keep genJetsToMC"]
out.outputCommands = ["keep *", "drop genParticles", "drop genVertices"]

############################################################
#
# Run modules
#
############################################################
# Run Pythia + Delphes
ApplicationMgr( TopAlg = [ pythia8gen, hepmc_converter, delphessim, out ],
                EvtSel = 'NONE',
                EvtMax = nEvents,
                ExtSvc = [podioEvent])

ApplicationMgr.EvtMax=10
예제 #10
0
]
overlay.PileUpTool = pileuptool

from Configurables import CreateExampleEventData
producer = CreateExampleEventData()
producer.magicNumberOffset = 300
producer.trackhits.Path = "trackHits"
producer.positionedtrackhits.Path = "positionedTrackHits"
producer.calohits.Path = "caloHits"
producer.positionedcalohits.Path = "positionedCaloHits"
producer.genparticles.Path = "genParticles"
producer.genvertices.Path = "genVertices"

# PODIO algorithm
from Configurables import PodioOutput
out = PodioOutput("out", OutputLevel=DEBUG)
out.outputCommands = ["drop *", "keep overlaid*"]
out.filename = "dummyEventDataOverlaid.root"

# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(
    TopAlg=[producer, overlay, out],
    EvtSel='NONE',
    EvtMax=300,
    ExtSvc=[
        podioevent,
    ],
    OutputLevel=DEBUG,
)
예제 #11
0
                                          OutputLevel = INFO)
positionsHcalFwd = CreateCellPositions("positionsHcalFwd",
                                          positionsTool=HCalFwdcells,
                                          hits = prefix+"HCalFwdCells",
                                          positionedHits = "HCalFwdCellPositions",
                                          OutputLevel = INFO)
if addMuons:
    positionsTailCatcher = CreateCellPositions("positionsTailCatcher",
                                               positionsTool=TailCatchercells,
                                               hits = "TailCatcherCells",
                                               positionedHits = "TailCatcherCellPositions",
                                               OutputLevel = INFO)

# PODIO algorithm
out = PodioOutput("out", OutputLevel=DEBUG)
out.outputCommands = ["keep *","drop "+prefix+"ECalBarrelCells","drop "+prefix+"ECalEndcapCells","drop "+prefix+"ECalFwdCells","drop "+prefix+"HCalBarrelCells", "drop "+prefix+"HCalExtBarrelCells", "drop "+prefix+"HCalEndcapCells", "drop "+prefix+"HCalFwdCells"]
if addMuons:
    out.outputCommands += ["drop TailCatcherCells"]
out.filename = output_name

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
positionsEcalEndcap.AuditExecute = True
positionsEcalFwd.AuditExecute = True
positionsHcalBarrel.AuditExecute = True
positionsHcalExtBarrel.AuditExecute = True
예제 #12
0
delphessim.DataOutputs.recMETs.Path = "recMETs"
delphessim.DataOutputs.recMuonsToMC.Path = "recMuonsToMC"
delphessim.DataOutputs.recElectronsToMC.Path = "recElectronsToMC"
delphessim.DataOutputs.recChargedToMC.Path = "recChargedToMC"
delphessim.DataOutputs.recNeutralToMC.Path = "recNeutralToMC"
delphessim.DataOutputs.recPhotonsToMC.Path = "recPhotonsToMC"
delphessim.DataOutputs.recJetsToMC.Path = "recJetsToMC"

## FCC event-data model output -> define objects to be written out
from Configurables import PodioOutput

out = PodioOutput("out", OutputLevel=messageLevelOut)
out.filename = "FCCDelphesOutput.root"
out.outputCommands = [
    "drop *", "keep genParticles", "keep genVertices", "keep recMuons",
    "keep recElectrons", "keep recCharged", "keep recPhotons", "keep recJets",
    "keep recMETs", "keep recMuonsToMC", "keep recElectronsToMC",
    "keep recChargedToMC", "keep recPhotonsToMC", "keep recJetsToMC"
]
#out.outputCommands = ["keep *"]

############################################################
#
# Run modules
#
############################################################

# Run Pythia + Delphes
if delphesHepMCInFile == "":
    ApplicationMgr(TopAlg=[pythia8gen, delphessim, out],
                   EvtSel='NONE',
                   EvtMax=nEvents,
예제 #13
0
hepmc_converter2 = HepMCToEDMConverter("Converter2")
hepmc_converter2.hepmc.Path = "hepmc2"
hepmc_converter2.genparticles.Path = "pileupGenParticles"
hepmc_converter2.genvertices.Path = "pileupGenVertices"

from Configurables import GenMerge

mergealg = GenMerge()
mergealg.pileupGenParticles.Path = "pileupGenParticles"
mergealg.pileupGenVertices.Path = "pileupGenVertices"
mergealg.signalGenParticles.Path = "signalGenParticles"
mergealg.signalGenVertices.Path = "signalGenVertices"
mergealg.allGenParticles.Path = "allGenParticles"
mergealg.allGenVertices.Path = "allGenVertices"

from Configurables import PodioOutput
out = PodioOutput("out", filename="output_gen.root")
out.outputCommands = [
    "keep *", "drop pileupGenVertices", "drop pileupGenParticles",
    "drop signalGenVertices", "drop signalGenParticles"
]

ApplicationMgr(
    TopAlg=[gun1, gun2, hepmc_converter1, hepmc_converter2, mergealg, out],
    EvtSel='NONE',
    ExtSvc=[podioevent],
    EvtMax=1,
    OutputLevel=VERBOSE,
)
예제 #14
0
delphessim.DataOutputs.recNeutralToMC.Path     = "recNeutralToMC"
delphessim.DataOutputs.recPhotonsToMC.Path     = "recPhotonsToMC"
delphessim.DataOutputs.recJetsToMC.Path        = "recJetsToMC" 

## FCC event-data model output -> define objects to be written out
from Configurables import PodioOutput

out = PodioOutput("out",OutputLevel=messageLevelOut)
out.filename       = "FCCDelphesOutput.root"
out.outputCommands = ["drop *",
                      "keep genParticles",
                      "keep genVertices",
                      "keep recMuons",
                      "keep recElectrons",
                      "keep recCharged",
                      "keep recPhotons",
                      "keep recJets",
                      "keep recMETs",
                      "keep recMuonsToMC",
                      "keep recElectronsToMC",
                      "keep recChargedToMC",
                      "keep recPhotonsToMC",
                      "keep recJetsToMC"]
#out.outputCommands = ["keep *"]

############################################################
#
# Run modules
#
############################################################

# Run Pythia + Delphes
예제 #15
0
overlay.randomizePileup = False
overlay.mergeTools = ["PileupParticlesMergeTool/MyPileupParticlesMergeTool",
  "PileupTrackHitMergeTool/MyTrackHitsMergeTool"]
overlay.PileUpTool = pileuptool

from Configurables import CreateExampleEventData
producer = CreateExampleEventData()
producer.magicNumberOffset = 300
producer.trackhits.Path = "trackHits"
producer.positionedtrackhits.Path = "positionedTrackHits"
producer.calohits.Path = "caloHits"
producer.positionedcalohits.Path = "positionedCaloHits"
producer.genparticles.Path = "genParticles"
producer.genvertices.Path = "genVertices"

# PODIO algorithm
from Configurables import PodioOutput
out = PodioOutput("out",
                   OutputLevel=DEBUG)
out.outputCommands = ["drop *", "keep overlaid*"]
out.filename = "dummyEventDataOverlaid.root"

# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr( TopAlg=[producer, overlay, out],
                EvtSel='NONE',
                EvtMax=300,
                ExtSvc=[podioevent,],
                OutputLevel=DEBUG,
  )
예제 #16
0
mergelayersHcalFwd.inhits.Path = "HCalFwdHits"
mergelayersHcalFwd.outhits.Path = "mergedHCalFwdHits"

createHcalFwdCells = CreateCaloCells("CreateHcalFwdCaloCells",
                                     doCellCalibration=True,
                                     calibTool=calibHcalFwd,
                                     addCellNoise=False,
                                     filterCellNoise=False,
                                     OutputLevel=INFO)
createHcalFwdCells.hits.Path = "mergedHCalFwdHits"
createHcalFwdCells.cells.Path = "HCalFwdCells"

out = PodioOutput("out", OutputLevel=INFO)
out.outputCommands = [
    "drop *", "keep ECalBarrelCells", "keep ECalEndcapCells",
    "keep ECalFwdCells", "keep newHCalCells", "keep newExtHCalCells",
    "keep HCalEndcapCells", "keep HCalFwdCells"
]
out.filename = "output_fullCalo_SimAndDigi_e50GeV_" + str(
    num_events) + "events.root"

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
예제 #17
0
파일: merge.py 프로젝트: HEP-FCC/FCCSW
hepmc_converter2 = HepMCToEDMConverter("Converter2")
hepmc_converter2.hepmc.Path="hepmc2"
hepmc_converter2.genparticles.Path="pileupGenParticles"
hepmc_converter2.genvertices.Path="pileupGenVertices"


from Configurables import GenMerge

mergealg = GenMerge()
mergealg.pileupGenParticles.Path="pileupGenParticles"
mergealg.pileupGenVertices.Path="pileupGenVertices"
mergealg.signalGenParticles.Path="signalGenParticles"
mergealg.signalGenVertices.Path="signalGenVertices"
mergealg.allGenParticles.Path="allGenParticles"
mergealg.allGenVertices.Path="allGenVertices"


from Configurables import PodioOutput
out = PodioOutput("out", filename = "output_gen.root")
out.outputCommands = ["keep *", "drop pileupGenVertices", "drop pileupGenParticles", "drop signalGenVertices", "drop signalGenParticles"]

ApplicationMgr(
               TopAlg=[gun1, gun2, hepmc_converter1, hepmc_converter2, mergealg, out],
               EvtSel='NONE',
               ExtSvc = [podioevent],
               EvtMax=1,
               OutputLevel=VERBOSE,
)

from Configurables import CreateCaloCells
createcells = CreateCaloCells("CreateCaloCells",
                              geometryTool=barrelHcalGeometry,
                              calibTool=calibHcells,
                              doCellCalibration=True,
                              addCellNoise=True,
                              filterCellNoise=False,
                              noiseTool=noise,
                              OutputLevel=DEBUG)
createcells.hits.Path = "HCalBarrelCellsStep1"
createcells.cells.Path = "HCalCells"

out = PodioOutput("out",
                  filename="output_HCalCells_digitisation_allCells.root",
                  OutputLevel=DEBUG)
out.outputCommands = ["keep *", "drop HCalBarrelCellsStep1"]

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
rewriteHCal.AuditExecute = True
createcells.AuditExecute = True
out.AuditExecute = True

ApplicationMgr(
    TopAlg=[podioinput, rewriteHCal, createcells, out],
    EvtSel='NONE',
    EvtMax=1,
예제 #19
0
pandoralg.ECalToHadGeVCalibrationBarrel = 1.12  #very small effect
pandoralg.ECalToHadGeVCalibrationEndCap = 1.12
pandoralg.HCalToHadGeVCalibration = 1.07
pandoralg.MuonToMipCalibration = 10.0
pandoralg.DigitalMuonHits = 0
pandoralg.MaxHCalHitHadronicEnergy = 1.0
pandoralg.UseOldTrackStateCalculation = 0
pandoralg.AbsorberRadLengthECal = 0.2854
pandoralg.AbsorberIntLengthECal = 0.0101
pandoralg.AbsorberRadLengthHCal = 0.0569
pandoralg.AbsorberIntLengthHCal = 0.006
pandoralg.AbsorberRadLengthOther = 0.0569
pandoralg.AbsorberIntLengthOther = 0.006

##############################################################################

# write PODIO file
from Configurables import PodioOutput
write = PodioOutput("write")
write.filename = "pan_test.root"
write.outputCommands = ["keep *"]

# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(TopAlg=[read, pandoralg],
               EvtSel='NONE',
               EvtMax=10,
               ExtSvc=[dsvc, gearSvc],
               HistogramPersistency="ROOT",
               OutputLevel=INFO)
예제 #20
0
                                         OutputLevel= DEBUG)

from Configurables import CreateCaloCells
createcells = CreateCaloCells("CreateCaloCells",
                              geometryTool = barrelHcalGeometry,
                              calibTool=calibHcells,
                              doCellCalibration = True,
                              addCellNoise = True, filterCellNoise = False,
                              noiseTool = noise,
                              OutputLevel = DEBUG)
createcells.hits.Path="HCalBarrelCellsStep1"
createcells.cells.Path="HCalCells"

out = PodioOutput("out", filename="output_HCalCells_digitisation_allCells.root",
                   OutputLevel = DEBUG)
out.outputCommands = ["keep *", "drop HCalBarrelCellsStep1"]

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
rewriteHCal.AuditExecute = True
createcells.AuditExecute = True
out.AuditExecute = True

ApplicationMgr(
    TopAlg = [podioinput,
              rewriteHCal,
              createcells,
    positionsHCalBarrelTool=HCalBcells,
    positionsHCalExtBarrelTool=HCalExtBcells,
    positionsEMECTool=EMECcells,
    positionsHECTool=HECcells,
    positionsEMFwdTool=ECalFwdcells,
    positionsHFwdTool=HCalFwdcells,
    hits="caloClusterBarrelCells",
    positionedHits="caloClusterBarrelCellPositions",
    OutputLevel=INFO)

out = PodioOutput("out",
                  filename="output_BarrelTopo_50GeVe_3ev.root",
                  OutputLevel=DEBUG)
out.outputCommands = [
    "drop *", "keep GenParticles", "keep GenVertices",
    "keep caloClustersBarrel", "keep caloClusterBarrelCells",
    "keep caloClusterBarrelCellPositions"
]

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
createemptycells.AuditExecute = True
createTopoClusters.AuditExecute = True
positionsClusterBarrel.AuditExecute = True
out.AuditExecute = True

ApplicationMgr(
예제 #22
0
createTailCatcherCells = CreateCaloCells("CreateTailCatcherCells",
                                         doCellCalibration=False,
                                         addCellNoise=False,
                                         filterCellNoise=False,
                                         OutputLevel=INFO,
                                         hits="TailCatcherHits",
                                         cells="TailCatcherCells")

# PODIO algorithm
from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput
podioevent = FCCDataSvc("EventDataSvc")
out = PodioOutput("out")
out.outputCommands = [
    "drop *", "keep GenParticles", "keep GenVertices", "keep TrackerHits",
    "keep TrackerPositionedHits", "keep TrackerDigiPostPoint",
    "keep ECalBarrelCells", "keep ECalEndcapCells", "keep ECalFwdCells",
    "keep HCalBarrelCells", "keep HCalExtBarrelCells", "keep HCalEndcapCells",
    "keep HCalFwdCells", "keep TailCatcherCells"
]
out.filename = output_name

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
예제 #23
0
                                        particleName="e-",
                                        energyMin=energy,
                                        energyMax=energy,
                                        etaMin=-0.36,
                                        etaMax=0.36,
                                        OutputLevel=DEBUG)

geantsim = SimG4Alg("SimG4Alg",
                    outputs=["SimG4SaveCalHits/saveHCalHits"],
                    eventProvider=pgun,
                    OutputLevel=DEBUG)

# PODIO algorithm
from Configurables import PodioOutput
out = PodioOutput("out", OutputLevel=DEBUG)
out.outputCommands = ["keep *"]
out.filename = "output_hcalSim_e" + str(int(
    energy / 1000)) + "GeV_eta036_1events.root"

#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True

# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(
    TopAlg=[geantsim, out],
    EvtSel='NONE',
예제 #24
0
                                        filterCellNoise=False)
createHcalEndcapCells.hits.Path = "pileupHCalEndcapCells"
createHcalEndcapCells.cells.Path = "mergedHCalEndcapCells"
createHcalFwdCells = CreateCaloCells("CreateHcalFwdCells",
                                     doCellCalibration=False,
                                     addCellNoise=False,
                                     filterCellNoise=False)
createHcalFwdCells.hits.Path = "pileupHCalFwdCells"
createHcalFwdCells.cells.Path = "mergedHCalFwdCells"

# PODIO algorithm
from Configurables import PodioOutput
out = PodioOutput("out")
out.outputCommands = [
    "drop *", "keep mergedECalBarrelCells", "keep mergedECalEndcapCells",
    "keep mergedECalFwdCells", "keep mergedHCalBarrelCells",
    "keep mergedHCalExtBarrelCells", "keep mergedHCalEndcapCells",
    "keep mergedHCalFwdCells"
]
out.filename = output_name

list_of_algorithms += [
    overlay, createEcalBarrelCells, createEcalEndcapCells, createEcalFwdCells,
    createHcalBarrelCells, createHcalExtBarrelCells, createHcalEndcapCells,
    createHcalFwdCells, out
]

# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(TopAlg=list_of_algorithms,
               EvtSel='NONE',
               EvtMax=num_events,