Python CellPositionsHCalBarrelTool Examples

Example #1

File: preparePileup.py Project: faltovaj/FCCSimJobs

# active material volume name
hcalVolumeName = ["moduleVolume", "wedgeVolume", "layerVolume"]
# ECAL bitfield names & values 
hcalFieldNames=["system"]
hcalFieldValues=[8]

##############################################################################################################                                                                                           #######                                       CELL POSITIONS  TOOLS                              #############                                                                                           ##############################################################################################################                                                                                           # Configure tools for calo cell positions                                                                            
from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelNoSegTool, CellPositionsHCalBarrelTool
ECalBcells = CellPositionsECalBarrelTool("CellPositionsECalBarrel",
                                         readoutName = ecalBarrelReadoutNamePhiEta,
                                         OutputLevel = INFO)
HCalBcellVols = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalBarrelVols",
                                                 readoutName = hcalBarrelReadoutName,
                                                 OutputLevel = INFO)
HCalBsegcells = CellPositionsHCalBarrelTool("CellPositionsHCalSegBarrel",
                                            readoutName = hcalBarrelReadoutNamePhiEta,
                                            radii = [2910.5, 3010.5, 3135.5, 3285.5, 3435.5, 3585.5, 3785.5, 4135.5, 4285.5, 4535.5],
                                            OutputLevel = INFO)

##############################################################################################################                                                                                           #######                                       RESEGMENT HCAL                                   #############                                                                                             ##############################################################################################################                                                                                                                              
from Configurables import CreateVolumeCaloPositions,RedoSegmentation,CreateCaloCells,CreateCellPositions
# Create cells in HCal                                           
# 2. step - rewrite the cellId using the Phi-Eta segmentation    
# 3. step - merge new cells corresponding to eta-phi segmentation
# Hcal barrel cell positions                                                                                                                                                                             
posHcalBarrel = CreateCellPositions("posHcalBarrel",
                                    positionsTool=HCalBcellVols,
                                    hits = prefix+"HCalBarrelCells",
                                    positionedHits = "HCalBarrelPositions")

# Use Phi-Eta segmentation in Hcal barrel                       
resegmentHcalBarrel = RedoSegmentation("ReSegmentationHcal",

Example #2

File: runBarrelCaloSystem_ReconstructionTopoClusters_electrNoise_fullGran_coneSelection.py Project: shehu0/FCCSW

smear = SimG4SmearGenParticles("smear",
                               inParticles = "GenParticles",
                               smearedParticles = "SmearedParticles",
                               smearTool = smeartool,
                               OutputLevel=DEBUG)

#Configure tools for calo reconstruction
from Configurables import ConstNoiseTool 
noiseTool = ConstNoiseTool("ConstNoiseTool")

#Configure tools for calo cell positions
from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelTool, CellPositionsHCalBarrelNoSegTool, CellPositionsCaloDiscsTool, CellPositionsTailCatcherTool 
ECalBcells = CellPositionsECalBarrelTool("CellPositionsECalBarrel", 
                                         readoutName = ecalBarrelReadoutName)
HCalBcells = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel",
                                         readoutName = hcalBarrelReadoutName, #"BarHCal_Readout_phieta",
                                         radii = [291.05, 301.05, 313.55, 328.55, 343.55, 358.55, 378.55, 413.55, 428.55, 453.55],
                                         OutputLevel = INFO)
HCalExtBcells = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalExtBarrel",
                                                readoutName = hcalExtBarrelReadoutName) 
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC",
                                       readoutName = ecalEndcapReadoutName,
                                       OutputLevel = INFO)
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd",
                                          readoutName = ecalFwdReadoutName,
                                          OutputLevel = INFO)
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC",
                                      readoutName = hcalEndcapReadoutName,
                                      OutputLevel = INFO)
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd",
                                          readoutName = hcalFwdReadoutName,
                                          OutputLevel = INFO)

Example #3

File: runFullCaloSystem_ReconstructionSW_noiseFromFile.py Project: simonecid/FCCSW

from Configurables import CellPositionsCaloDiscsTool
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC")
EMECcells.readoutName = ecalBarrelReadoutName

from Configurables import CellPositionsCaloDiscsTool
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd")
ECalFwdcells.readoutName = ecalFwdReadoutName
ECalFwdcells.OutputLevel = INFO

from Configurables import CellPositionsHCalBarrelNoSegTool
HCalBcellVols = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalBarrelVols")
HCalBcellVols.readoutName = "HCalBarrelReadout"
HCalBcellVols.OutputLevel = INFO

from Configurables import CellPositionsHCalBarrelTool
HCalBsegcells = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel")
HCalBsegcells.readoutName = "BarHCal_Readout_phieta"
HCalBsegcells.radii = [291.05, 301.05, 313.55, 328.55, 343.55, 358.55, 378.55, 413.55, 428.55, 453.55]
HCalBsegcells.OutputLevel = INFO

from Configurables import CellPositionsCaloDiscsTool
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC")
HECcells.readoutName = hcalEndcapReadoutName

from Configurables import CellPositionsCaloDiscsTool
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd")
HCalFwdcells.readoutName = hcalFwdReadoutName


# additionally for HCal                                 
from Configurables import RewriteBitfield

Example #4

File: runBarrelCaloSystem_ReconstructionTopoClusters_electrNoise.py Project: tguillemfcc/FCCSW

#Configure tools for calo cell positions
from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelTool, CellPositionsCaloDiscsTool, CellPositionsTailCatcherTool
ECalBcells = CellPositionsECalBarrelTool("CellPositionsECalBarrel",
                                         readoutName=ecalBarrelReadoutName,
                                         OutputLevel=INFO)
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC",
                                       readoutName=ecalEndcapReadoutName,
                                       OutputLevel=INFO)
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd",
                                          readoutName=ecalFwdReadoutName,
                                          OutputLevel=INFO)
HCalBcells = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel",
                                         readoutName=hcalBarrelReadoutName,
                                         radii=[
                                             291.05, 301.05, 313.55, 328.55,
                                             343.55, 358.55, 378.55, 403.55,
                                             428.55, 453.55
                                         ],
                                         OutputLevel=INFO)
HCalExtBcells = CellPositionsHCalBarrelTool(
    "CellPositionsHCalExtBarrel",
    readoutName=hcalExtBarrelReadoutName,
    radii=[
        356.05, 373.55, 398.55, 423.55, 291.05, 301.05, 313.55, 328.55, 348.55,
        373.55, 398.55, 423.55
    ],
    OutputLevel=INFO)
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC",
                                      readoutName=hcalEndcapReadoutName,
                                      OutputLevel=INFO)
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd",

Example #5

File: recPositions.py Project: yenari/FCCSimJobs

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

#Configure tools for calo cell positions
from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelTool, CellPositionsCaloDiscsTool, CellPositionsCaloDiscsTool, CellPositionsTailCatcherTool
ECalBcells = CellPositionsECalBarrelTool(
    "CellPositionsECalBarrel",
    readoutName=ecalBarrelReadoutNamePhiEta,
    OutputLevel=INFO)
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC",
                                       readoutName=ecalEndcapReadoutName,
                                       OutputLevel=INFO)
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd",
                                          readoutName=ecalFwdReadoutName,
                                          OutputLevel=INFO)
HCalBcells = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel",
                                         readoutName=hcalBarrelReadoutName,
                                         OutputLevel=INFO)
HCalExtBcells = CellPositionsHCalBarrelTool(
    "CellPositionsHCalExtBarrel",
    readoutName=hcalExtBarrelReadoutName,
    OutputLevel=INFO)
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC",
                                      readoutName=hcalEndcapReadoutName,
                                      OutputLevel=INFO)
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd",
                                          readoutName=hcalFwdReadoutName,
                                          OutputLevel=INFO)
TailCatchercells = CellPositionsTailCatcherTool(
    "CellPositionsTailCatcher",
    readoutName=tailCatcherReadoutName,
    centralRadius=901.5,

Example #6

                                    readoutName = ecalBarrelReadoutName,
                                    OutputLevel = INFO)
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC",
                                    readoutName = ecalEndcapReadoutName,
                                    OutputLevel = INFO)
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd",
                                        readoutName = ecalFwdReadoutName,
                                        OutputLevel = INFO)
HCalBcells = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalBarrel",
                                              readoutName = hcalBarrelReadoutName,
                                              OutputLevel = INFO)
HCalExtBcells = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalExtBarrel",
                                                 readoutName = hcalExtBarrelReadoutName,
                                                 OutputLevel = INFO)
HCalBsegcells = CellPositionsHCalBarrelTool("CellPositionsHCalSegBarrel",
                                            readoutName = hcalBarrelReadoutNamePhiEta,
                                            radii = [2910.5, 3010.5, 3135.5, 3285.5, 3435.5, 3585.5, 3785.5, 4035.5, 4285.5, 4535.5],
                                            OutputLevel = INFO)
HCalExtBsegcells = CellPositionsHCalBarrelTool("CellPositionsHCalExtSegBarrel",
                                               readoutName = hcalExtBarrelReadoutNamePhiEta,
                                               radii = [ 356.05
                                                         , 373.55
                                                         , 398.55
                                                         , 423.55
                                                         , 291.05
                                                         , 301.05
                                                         , 313.55
                                                         , 328.55
                                                         , 348.55
                                                         , 373.55
                                                         , 398.55
                                                         , 423.55