Beispiel #1
0
#!/usr/bin/env python

from Gaudi.Configuration import *

from Configurables import k4DataSvc
dsvc = k4DataSvc("EventDataSvc")

# read LCIO files
from Configurables import LCIOInput
lcioinput = LCIOInput("LCIOInput")

import glob

# inputlists = ["/cefs/higgs/yudan/CEPC240/Reco_tpc_1800/qqh/Reco_qqh__00001.slcio"]
inputlists = ["/cefs/higgs/yudan/CEPC240/Reco_tpc_1800/qqh/Reco_qqh__00003.slcio"]
# inputlists = glob.glob("/cefs/higgs/yudan/CEPC240/Reco_tpc_1800/qqh/Reco_qqh__*.slcio")

lcioinput.inputs = inputlists
lcioinput.collections = [
    "MCParticle:MCParticle",
    "CalorimeterHit:ECALBarrel",
    "CalorimeterHit:ECALEndcap",
    "CalorimeterHit:HCALBarrel",
    "CalorimeterHit:HCALEndcap",
    "CalorimeterHit:HCALOther",
    "ReconstructedParticle:AncientPFOs",
    "ReconstructedParticle:ArborLICHPFOs"
]

from Configurables import TotalInvMass
total_inv_mass = TotalInvMass("TotalInvMass")
Beispiel #2
0
#!/usr/bin/env python

from Gaudi.Configuration import *
from Configurables import k4DataSvc
dsvc = k4DataSvc("EventDataSvc")

# read LCIO files
from Configurables import LCIOInput
read = LCIOInput("read")
read.inputs = [
    #"/cefs/higgs/wxfang/cepc/Pandora/ele_lcio_full_det_rec/reco_sim_0.slcio"
    "/cefs/higgs/wxfang/cepc/Pandora/ele_lcio_full_det_rec/reco_sim_34.slcio"
]
read.collections = [
    "MCParticle:MCParticle",
    #"SimCalorimeterHit:EcalBarrelSiliconCollection",
    "CalorimeterHit:ECALBarrel",
    "CalorimeterHit:ECALEndcap",
    "CalorimeterHit:ECALOther",
    ######### HCAL will effect the reco efficiency close to gap region ######
    "CalorimeterHit:HCALBarrel",
    "CalorimeterHit:HCALEndcap",
    "CalorimeterHit:HCALOther",
    #"TrackerHit:VXDTrackerHits",
    #"TrackerHit:SITTrackerHits",
    "TrackerHit:SITSpacePoints",
    "TrackerHit:TPCTrackerHits",
    #"TrackerHit:SETTrackerHits",
    "TrackerHit:SETSpacePoints",
    #"TrackerHit:FTDStripTrackerHits",
    "TrackerHit:FTDSpacePoints",
Beispiel #3
0
from Configurables import LCIOInput
lcioinput = LCIOInput(
    "LCIOReader",
    collections=[
        "EventHeader", "MCParticle", "TPCCollection"
        #"EventHeader",
        #"MCParticle",
        #"COILCollection",
        #"EcalBarrelSiliconCollection",
        #"EcalBarrelSiliconPreShowerCollection",
        #"EcalEndcapRingCollection",
        #"EcalEndcapRingPreShowerCollection",
        #"EcalEndcapSiliconCollection",
        #"EcalEndcapSiliconPreShowerCollection",
        #"FTD_PIXELCollection",
        #"FTD_STRIPCollection",
        #"HcalBarrelCollection",
        #"HcalEndCapRingsCollection",
        #"HcalEndCapsCollection",
        #"LumiCalCollection",
        #"MuonBarrelCollection",
        #"MuonEndCapCollection",
        #"SETCollection",
        #"SITCollection",
        #"TPCCollection",
        #"TPCSpacePointCollection",
        #"VXDCollection"
    ])
lcioinput.DataSvc = svcname