diff --git a/Examples/options/LCIO_read_pan.py b/Examples/options/LCIO_read_pan.py index 97b5f9e974a0db059bb5a7d4a1da503d07887a16..e29a752b9981cccafc7d379082053fb313b14917 100644 --- a/Examples/options/LCIO_read_pan.py +++ b/Examples/options/LCIO_read_pan.py @@ -10,57 +10,76 @@ read = LCIOInput("read") read.inputs = [ #"/cefs/data/FullSim/CEPC240/CEPC_v4/higgs/smart_final_states/E240.Pffh_invi.e0.p0.whizard195//ffh_inv.e0.p0.00001_1000_sim.slcio" #"/junofs/users/wxfang/CEPC/CEPCOFF/doReco/reco_output/nnh_aa.e0.p0.00010_000000_rec.slcio" -"/cefs/higgs/wxfang/cepc/Pandora/CaloDigi/gamma/Digi_sim_0.slcio" +#"/cefs/higgs/wxfang/cepc/Pandora/CaloDigi/gamma/Digi_sim_0.slcio" + + + ] read.collections = { - #"COILCollection" : "SimTrackerHit", - #"EcalBarrelSiliconCollection" : "SimCalorimeterHit", - "MCParticle" : "MCParticle", - "ECALBarrel" : "CalorimeterHit", - "ECALEndcap" : "CalorimeterHit", - "ECALOther" : "CalorimeterHit", - "HCALBarrel" : "CalorimeterHit", - "HCALEndcap" : "CalorimeterHit", - "HCALOther" : "CalorimeterHit", - "MUON" : "CalorimeterHit", - "LCAL" : "CalorimeterHit", - "LHCAL" : "CalorimeterHit", - "BCAL" : "CalorimeterHit", - #"MarlinTrkTracks" : "Track" - #"TPCCollection" : "SimTrackerHit", - #"VXDCollection" : "SimTrackerHit" + "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", + #"TrackerHit:FTDPixelTrackerHits", + "Track:ClupatraTrackSegments", + "Track:ClupatraTracks", + "Track:ForwardTracks", + "Track:SiTracks", + "Track:SubsetTracks", + "Track:MarlinTrkTracks", + "Vertex:KinkVertices", + "Vertex:ProngVertices", + "Vertex:V0Vertices", + "ReconstructedParticle:KinkRecoParticles", + "ReconstructedParticle:ProngRecoParticles", + "ReconstructedParticle:V0RecoParticles" } ############################################################################## from Configurables import GearSvc gearSvc = GearSvc("GearSvc") -gearSvc.GearXMLFile = "/junofs/users/wxfang/CEPC/CEPCOFF/doSim/fullDet/GearOutput.xml" +gearSvc.GearXMLFile = "../Detector/DetCEPCv4/compact/FullDetGear.xml" ############################################################################## from Configurables import PandoraPFAlg pandoralg = PandoraPFAlg("PandoraPFAlg") -## KEEP same with lcioinput name for the ReadXXX ########### -pandoralg.ReadMCParticle = "MCParticle" -pandoralg.ReadECALBarrel = "ECALBarrel" -pandoralg.ReadECALEndcap = "ECALEndcap" -pandoralg.ReadECALOther = "ECALOther" -pandoralg.ReadHCALBarrel = "HCALBarrel" -pandoralg.ReadHCALEndcap = "HCALEndcap" -pandoralg.ReadHCALOther = "HCALOther" -pandoralg.ReadMUON = "MUON" -pandoralg.ReadLCAL = "LCAL" -pandoralg.ReadLHCAL = "LHCAL" -pandoralg.ReadBCAL = "BCAL" -pandoralg.ReadKinkVertices = "KinkVertices" -pandoralg.ReadProngVertices = "ProngVertices" -pandoralg.ReadSplitVertices = "SplitVertices" -pandoralg.ReadV0Vertices = "V0Vertices" -pandoralg.ReadTracks = "MarlinTrkTracks" +pandoralg.collections = [ + "MCParticle:MCParticle", + "CalorimeterHit:ECALBarrel", + "CalorimeterHit:ECALEndcap", + "CalorimeterHit:ECALOther" , + "CalorimeterHit:HCALBarrel", + "CalorimeterHit:HCALEndcap", + "CalorimeterHit:HCALOther" , + "CalorimeterHit:MUON", + "CalorimeterHit:LCAL", + "CalorimeterHit:LHCAL", + "CalorimeterHit:BCAL", + "Vertex:KinkVertices", + "Vertex:ProngVertices", + "Vertex:SplitVertices", + "Vertex:V0Vertices", + "Track:MarlinTrkTracks", + "MCRecoCaloAssociation:RecoCaloAssociation_ECALBarrel" + ] pandoralg.WriteClusterCollection = "PandoraClusters" pandoralg.WriteReconstructedParticleCollection = "PandoraPFOs" pandoralg.WriteVertexCollection = "PandoraPFANewStartVertices" -pandoralg.AnaOutput = "/cefs/higgs/wxfang/cepc/Pandora/Ana/gamma/Ana_gamma_test.root" +pandoralg.AnaOutput = "Pandora_Ana.root" -pandoralg.PandoraSettingsDefault_xml = "/junofs/users/wxfang/MyGit/MarlinPandora/scripts/PandoraSettingsDefault_wx.xml" +pandoralg.PandoraSettingsDefault_xml = "../Reconstruction/PFA/Pandora/PandoraSettingsDefault.xml" #### Do not chage the collection name, only add or delete ############### pandoralg.TrackCollections = ["MarlinTrkTracks"] pandoralg.ECalCaloHitCollections= ["ECALBarrel", "ECALEndcap", "ECALOther"] diff --git a/Examples/options/tut_detsim_pandora.py b/Examples/options/tut_detsim_pandora.py index b4ceaf07377898130aea5d39e6e304242cda1262..e1ac4a7fa9053fac9cd929609d125fac1ca5f727 100644 --- a/Examples/options/tut_detsim_pandora.py +++ b/Examples/options/tut_detsim_pandora.py @@ -134,24 +134,25 @@ gearSvc.GearXMLFile = "../Detector/DetCEPCv4/compact/FullDetGear.xml" from Configurables import PandoraPFAlg pandoralg = PandoraPFAlg("PandoraPFAlg") -## KEEP same with lcioinput name for the ReadXXX ########### -pandoralg.ReadMCParticle = "MCParticle" -pandoralg.ReadECALBarrel = "ECALBarrel" -pandoralg.ReadECALEndcap = "ECALEndcap" -pandoralg.ReadECALOther = "ECALOther" -pandoralg.ReadHCALBarrel = "HCALBarrel" -pandoralg.ReadHCALEndcap = "HCALEndcap" -pandoralg.ReadHCALOther = "HCALOther" -pandoralg.ReadMUON = "MUON" -pandoralg.ReadLCAL = "LCAL" -pandoralg.ReadLHCAL = "LHCAL" -pandoralg.ReadBCAL = "BCAL" -pandoralg.ReadKinkVertices = "KinkVertices" -pandoralg.ReadProngVertices = "ProngVertices" -pandoralg.ReadSplitVertices = "SplitVertices" -pandoralg.ReadV0Vertices = "V0Vertices" -pandoralg.ReadTracks = "MarlinTrkTracks" -pandoralg.MCRecoCaloAssociation = "RecoCaloAssociation_ECALBarrel" +pandoralg.collections = [ + "MCParticle:MCParticle", + "CalorimeterHit:ECALBarrel", + "CalorimeterHit:ECALEndcap", + "CalorimeterHit:ECALOther" , + "CalorimeterHit:HCALBarrel", + "CalorimeterHit:HCALEndcap", + "CalorimeterHit:HCALOther" , + "CalorimeterHit:MUON", + "CalorimeterHit:LCAL", + "CalorimeterHit:LHCAL", + "CalorimeterHit:BCAL", + "Vertex:KinkVertices", + "Vertex:ProngVertices", + "Vertex:SplitVertices", + "Vertex:V0Vertices", + "Track:MarlinTrkTracks", + "MCRecoCaloAssociation:RecoCaloAssociation_ECALBarrel" + ] pandoralg.WriteClusterCollection = "PandoraClusters" pandoralg.WriteReconstructedParticleCollection = "PandoraPFOs" pandoralg.WriteVertexCollection = "PandoraPFANewStartVertices" diff --git a/Reconstruction/PFA/Pandora/GaudiPandora/include/PandoraPFAlg.h b/Reconstruction/PFA/Pandora/GaudiPandora/include/PandoraPFAlg.h index 0e447790098de21001fb6c07cc3915af81ca8b8a..e357bd34b55532325758d9f45ac9ea0dbfb34868 100644 --- a/Reconstruction/PFA/Pandora/GaudiPandora/include/PandoraPFAlg.h +++ b/Reconstruction/PFA/Pandora/GaudiPandora/include/PandoraPFAlg.h @@ -275,26 +275,12 @@ protected: Gaudi::Property< std::string > m_AnaOutput{ this, "AnaOutput", "/junofs/users/wxfang/MyGit/CEPCSW/Reconstruction/PFA/Pandora/GaudiPandora/Ana.root" }; //###################### + std::map< std::string, std::string > m_collections; + Gaudi::Property<std::vector<std::string>> m_readCols{this, "collections", {}, "Places of collections to read"}; + //the map of collection name to its corresponding DataHandle + std::map<std::string, DataObjectHandleBase*> m_dataHandles; DataHandle<edm4hep::MCParticleCollection> m_mcParCol_r {"MCParticle", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_ECALBarrel_r{"ECALBarrel", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_ECALEndcap_r{"ECALEndcap", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_ECALOther_r {"ECALOther", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_HCALBarrel_r{"HCALBarrel", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_HCALEndcap_r{"HCALEndcap", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_HCALOther_r {"HCALOther", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_MUON_r {"MUON", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_LCAL_r {"LCAL", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_LHCAL_r {"LHCAL", Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::CalorimeterHitCollection> m_BCAL_r {"BCAL", Gaudi::DataHandle::Reader, this}; - - DataHandle<edm4hep::VertexCollection> m_KinkVertices_r {"KinkVertices",Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::VertexCollection> m_ProngVertices_r {"ProngVertices",Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::VertexCollection> m_SplitVertices_r {"SplitVertices",Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::VertexCollection> m_V0Vertices_r {"V0Vertices",Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::TrackCollection> m_MarlinTrkTracks_r {"MarlinTrkTracks",Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::MCRecoCaloAssociationCollection> m_MCRecoCaloAssociation_r {"MCRecoCaloAssociation",Gaudi::DataHandle::Reader, this}; - DataHandle<edm4hep::MCRecoTrackerAssociationCollection> m_MCRecoTrackerAssociation_r {"MCRecoTrackerAssociation",Gaudi::DataHandle::Reader, this}; DataHandle<edm4hep::ClusterCollection> m_ClusterCollection_w {"PandoraClusters",Gaudi::DataHandle::Writer, this}; DataHandle<edm4hep::ReconstructedParticleCollection> m_ReconstructedParticleCollection_w {"PandoraPFOs" ,Gaudi::DataHandle::Writer, this}; diff --git a/Reconstruction/PFA/Pandora/GaudiPandora/include/TrackCreator.h b/Reconstruction/PFA/Pandora/GaudiPandora/include/TrackCreator.h index de380e3f6d4b921a422ad6414d3d772b0f755820..134ff85bc5f08b7643155624e4e82910e7cbdce7 100644 --- a/Reconstruction/PFA/Pandora/GaudiPandora/include/TrackCreator.h +++ b/Reconstruction/PFA/Pandora/GaudiPandora/include/TrackCreator.h @@ -123,6 +123,7 @@ public: */ pandora::StatusCode CreateTrackAssociations(const CollectionMaps& collectionMaps); + const edm4hep::Track* GetTrackAddress(const CollectionMaps& collectionMaps, const edm4hep::ConstTrack& pTrack ); /** * @brief Create tracks, insert user code here * @@ -154,6 +155,7 @@ private: * @param pLCEvent the lcio event */ //pandora::StatusCode ExtractProngsAndSplits(const EVENT::LCEvent *const pLCEvent); + pandora::StatusCode ExtractProngsAndSplits(const CollectionMaps& collectionMaps); /** * @brief Extract v0 information from specified lcio collections @@ -161,6 +163,7 @@ private: * @param pLCEvent the lcio event */ //pandora::StatusCode ExtractV0s(const EVENT::LCEvent *const pLCEvent); + pandora::StatusCode ExtractV0s(const CollectionMaps& collectionMaps); /** * @brief Whether the track vertex conflicts with previously provided relationship information diff --git a/Reconstruction/PFA/Pandora/GaudiPandora/src/CaloHitCreator.cpp b/Reconstruction/PFA/Pandora/GaudiPandora/src/CaloHitCreator.cpp index 9c7fb1da1bcbc516ee597faf3b04d6f8eec33c4a..12120d3f56e3e8b2e8de9a03ddd0ceeb90dab51d 100644 --- a/Reconstruction/PFA/Pandora/GaudiPandora/src/CaloHitCreator.cpp +++ b/Reconstruction/PFA/Pandora/GaudiPandora/src/CaloHitCreator.cpp @@ -192,7 +192,8 @@ pandora::StatusCode CaloHitCreator::CreateECalCaloHits(const CollectionMaps& col caloHitParameters.m_hadronicEnergy = eCalToHadGeVEndCap * pCaloHit->getEnergy(); } - caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * eCalToMip * absorberCorrection; + //caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * eCalToMip * absorberCorrection; + caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * eCalToMip;//FIXME. is absorberCorrection it needed for digi input if (caloHitParameters.m_mipEquivalentEnergy.Get() < eCalMipThreshold) continue; @@ -282,7 +283,8 @@ pandora::StatusCode CaloHitCreator::CreateHCalCaloHits(const CollectionMaps& col this->GetEndCapCaloHitProperties(pCaloHit, endcapLayerLayout, caloHitParameters, absorberCorrection); } - caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_hCalToMip * absorberCorrection; + //caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_hCalToMip * absorberCorrection; + caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_hCalToMip; if (caloHitParameters.m_mipEquivalentEnergy.Get() < m_settings.m_hCalMipThreshold) continue; @@ -456,7 +458,8 @@ pandora::StatusCode CaloHitCreator::CreateLCalCaloHits(const CollectionMaps& col float absorberCorrection(1.); this->GetEndCapCaloHitProperties(pCaloHit, endcapLayerLayout, caloHitParameters, absorberCorrection); - caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_eCalToMip * absorberCorrection; + //caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_eCalToMip * absorberCorrection; + caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_eCalToMip; if (caloHitParameters.m_mipEquivalentEnergy.Get() < m_settings.m_eCalMipThreshold) continue; @@ -528,7 +531,8 @@ pandora::StatusCode CaloHitCreator::CreateLHCalCaloHits(const CollectionMaps& co float absorberCorrection(1.); this->GetEndCapCaloHitProperties(pCaloHit, endcapLayerLayout, caloHitParameters, absorberCorrection); - caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_hCalToMip * absorberCorrection; + //caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_hCalToMip * absorberCorrection; + caloHitParameters.m_mipEquivalentEnergy = pCaloHit->getEnergy() * m_settings.m_hCalToMip; if (caloHitParameters.m_mipEquivalentEnergy.Get() < m_settings.m_hCalMipThreshold) continue; diff --git a/Reconstruction/PFA/Pandora/GaudiPandora/src/PandoraPFAlg.cpp b/Reconstruction/PFA/Pandora/GaudiPandora/src/PandoraPFAlg.cpp index d2e9db69f1a2184b3bbf1f3c6d27d344a2432a2c..14d4ee1a0e82b492706e30bba8cce509823bbbe0 100644 --- a/Reconstruction/PFA/Pandora/GaudiPandora/src/PandoraPFAlg.cpp +++ b/Reconstruction/PFA/Pandora/GaudiPandora/src/PandoraPFAlg.cpp @@ -39,24 +39,6 @@ PandoraPFAlg::PandoraPFAlg(const std::string& name, ISvcLocator* svcLoc) { m_CollectionMaps = new CollectionMaps(); - declareProperty("ReadMCParticle" , m_mcParCol_r, "Handle of the MCParticle input collection" ); - declareProperty("ReadECALBarrel" , m_ECALBarrel_r, "Handle of the ECALBarrel input collection" ); - declareProperty("ReadECALEndcap" , m_ECALEndcap_r, "Handle of the ECALEndcap input collection" ); - declareProperty("ReadECALOther" , m_ECALOther_r, "Handle of the ECALOther input collection" ); - declareProperty("ReadHCALBarrel" , m_HCALBarrel_r, "Handle of the HCALBarrel input collection" ); - declareProperty("ReadHCALEndcap" , m_HCALEndcap_r, "Handle of the HCALEndcap input collection" ); - declareProperty("ReadHCALOther" , m_HCALOther_r, "Handle of the HCALOther input collection" ); - declareProperty("ReadMUON" , m_MUON_r, "Handle of the MUON input collection" ); - declareProperty("ReadLCAL" , m_LCAL_r, "Handle of the LCAL input collection" ); - declareProperty("ReadLHCAL" , m_LHCAL_r, "Handle of the LHCAL input collection" ); - declareProperty("ReadBCAL" , m_BCAL_r, "Handle of the BCAL input collection" ); - declareProperty("ReadKinkVertices" , m_KinkVertices_r, "Handle of the KinkVertices input collection" ); - declareProperty("ReadProngVertices" , m_ProngVertices_r, "Handle of the ProngVertices input collection" ); - declareProperty("ReadSplitVertices" , m_SplitVertices_r, "Handle of the SplitVertices input collection" ); - declareProperty("ReadV0Vertices" , m_V0Vertices_r, "Handle of the V0Vertices input collection" ); - declareProperty("ReadTracks" , m_MarlinTrkTracks_r, "Handle of the Tracks input collection" ); - declareProperty("MCRecoCaloAssociation" , m_MCRecoCaloAssociation_r, "Handle of the MCRecoCaloAssociation input collection" ); - declareProperty("MCRecoTrackerAssociation" , m_MCRecoTrackerAssociation_r, "Handle of the MCRecoTrackerAssociation input collection" ); declareProperty("WriteClusterCollection" , m_ClusterCollection_w, "Handle of the ClusterCollection output collection" ); declareProperty("WriteReconstructedParticleCollection", m_ReconstructedParticleCollection_w, "Handle of the ReconstructedParticleCollection output collection" ); declareProperty("WriteVertexCollection" , m_VertexCollection_w, "Handle of the VertexCollection output collection" ); @@ -123,6 +105,43 @@ StatusCode PandoraPFAlg::initialize() m_tree->Branch("m_mc_charge", &m_mc_charge); m_tree->Branch("m_hasConversion", &m_hasConversion); + for ( const auto& col : m_readCols ) { + auto seperater = col.find(':'); + std::string colType = col.substr(0, seperater); + std::string colName = col.substr(seperater+1); + m_collections[colName] = colType; + + if ( colType == "MCParticle" ) { + m_dataHandles[colName] = + new DataHandle<edm4hep::MCParticleCollection>(colName, Gaudi::DataHandle::Reader, this); + } + else if ( colType == "Track" ) { + m_dataHandles[colName] = + new DataHandle<edm4hep::TrackCollection>(colName, Gaudi::DataHandle::Reader, this); + } + else if ( colType == "CalorimeterHit" ) { + m_dataHandles[colName] = + new DataHandle<edm4hep::CalorimeterHitCollection>(colName, Gaudi::DataHandle::Reader, this); + } + else if ( colType == "Vertex" ) { + m_dataHandles[colName] = + new DataHandle<edm4hep::VertexCollection>(colName, Gaudi::DataHandle::Reader, this); + } + else if ( colType == "MCRecoTrackerAssociation" ) { + m_dataHandles[colName] = + new DataHandle<edm4hep::MCRecoTrackerAssociationCollection>(colName, Gaudi::DataHandle::Reader, this); + } + else if ( colType == "MCRecoCaloAssociation" ) { + m_dataHandles[colName] = + new DataHandle<edm4hep::MCRecoCaloAssociationCollection>(colName, Gaudi::DataHandle::Reader, this); + } + else { + error() << "invalid collection type: " << colType << endmsg; + return StatusCode::FAILURE; + } + } + + // XML file m_settings.m_pandoraSettingsXmlFile = m_PandoraSettingsXmlFile ; // Hadronic energy non-linearity correction @@ -418,160 +437,103 @@ collectionMap_CaloRel.clear(); collectionMap_TrkRel.clear(); } + StatusCode PandoraPFAlg::updateMap() { - const edm4hep::MCParticleCollection* MCParticle = nullptr; - const edm4hep::CalorimeterHitCollection* ECALBarrel = nullptr; - const edm4hep::CalorimeterHitCollection* ECALEndcap = nullptr; - const edm4hep::CalorimeterHitCollection* ECALOther = nullptr; - const edm4hep::CalorimeterHitCollection* HCALBarrel = nullptr; - const edm4hep::CalorimeterHitCollection* HCALEndcap = nullptr; - const edm4hep::CalorimeterHitCollection* HCALOther = nullptr; - const edm4hep::CalorimeterHitCollection* MUON = nullptr; - const edm4hep::CalorimeterHitCollection* LCAL = nullptr; - const edm4hep::CalorimeterHitCollection* LHCAL = nullptr; - const edm4hep::CalorimeterHitCollection* BCAL = nullptr; - const edm4hep::VertexCollection* KinkVertices = nullptr; - const edm4hep::VertexCollection* ProngVertices = nullptr; - const edm4hep::VertexCollection* SplitVertices = nullptr; - const edm4hep::VertexCollection* V0Vertices = nullptr; - const edm4hep::TrackCollection* MarlinTrkTracks = nullptr; - const edm4hep::MCRecoCaloAssociationCollection* mcRecoCaloAssociation = nullptr; - const edm4hep::MCRecoTrackerAssociationCollection* mcRecoTrackerAssociation = nullptr; - StatusCode sc = StatusCode::SUCCESS; - sc = getCol(m_mcParCol_r , MCParticle ); - sc = getCol(m_ECALBarrel_r, ECALBarrel ); - sc = getCol(m_ECALEndcap_r, ECALEndcap ); - sc = getCol(m_ECALOther_r , ECALOther ); - sc = getCol(m_HCALBarrel_r, HCALBarrel ); - sc = getCol(m_HCALEndcap_r, HCALEndcap ); - sc = getCol(m_HCALOther_r , HCALOther ); - sc = getCol(m_MUON_r , MUON ); - sc = getCol(m_LCAL_r , LCAL ); - sc = getCol(m_LHCAL_r , LHCAL ); - sc = getCol(m_BCAL_r , BCAL ); - sc = getCol(m_KinkVertices_r , KinkVertices ); - sc = getCol(m_ProngVertices_r , ProngVertices); - sc = getCol(m_SplitVertices_r , SplitVertices); - sc = getCol(m_V0Vertices_r , V0Vertices ); - sc = getCol(m_MarlinTrkTracks_r , MarlinTrkTracks ); - sc = getCol(m_MCRecoCaloAssociation_r , mcRecoCaloAssociation ); - sc = getCol(m_MCRecoTrackerAssociation_r , mcRecoTrackerAssociation); - - if (NULL != MCParticle ) - { - std::vector<edm4hep::MCParticle> v_mc; - m_CollectionMaps->collectionMap_MC ["MCParticle"] = v_mc; - for(unsigned int i=0 ; i< MCParticle->size(); i++) m_CollectionMaps->collectionMap_MC ["MCParticle"].push_back(MCParticle->at(i)); - } - if (NULL != ECALBarrel ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["ECALBarrel"] = v_cal ; - for(unsigned int i=0 ; i< ECALBarrel->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["ECALBarrel"].push_back(ECALBarrel->at(i)); - } - if (NULL != ECALEndcap ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["ECALEndcap"] = v_cal ; - for(unsigned int i=0 ; i< ECALEndcap->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["ECALEndcap"].push_back(ECALEndcap->at(i)); - } - if (NULL != ECALOther ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["ECALOther"] = v_cal ; - for(unsigned int i=0 ; i< ECALOther->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["ECALOther"].push_back(ECALOther->at(i)); - } - if (NULL != HCALBarrel ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["HCALBarrel"] = v_cal ; - for(unsigned int i=0 ; i< HCALBarrel->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["HCALBarrel"].push_back(HCALBarrel->at(i)); - } - if (NULL != HCALEndcap ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["HCALEndcap"] = v_cal ; - for(unsigned int i=0 ; i< HCALEndcap->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["HCALEndcap"].push_back(HCALEndcap->at(i)); - } - if (NULL != HCALOther ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["HCALOther"] = v_cal ; - for(unsigned int i=0 ; i< HCALOther->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["HCALOther"].push_back(HCALOther->at(i)); - } - if (NULL != MUON ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["MUON"] = v_cal ; - for(unsigned int i=0 ; i< MUON->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["MUON"].push_back(MUON->at(i)); - } - if (NULL != LCAL ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["LCAL"] = v_cal ; - for(unsigned int i=0 ; i< LCAL->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["LCAL"].push_back(LCAL->at(i)); - } - if (NULL != LHCAL ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["LHCAL"] = v_cal ; - for(unsigned int i=0 ; i< LHCAL->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["LHCAL"].push_back(LHCAL->at(i)); - } - if (NULL != BCAL ) - { - std::vector<edm4hep::CalorimeterHit> v_cal; - m_CollectionMaps->collectionMap_CaloHit["BCAL"] = v_cal ; - for(unsigned int i=0 ; i< BCAL->size(); i++) m_CollectionMaps->collectionMap_CaloHit ["BCAL"].push_back(BCAL->at(i)); - } - if (NULL != KinkVertices ) - { - std::vector<edm4hep::Vertex> v_cal; - m_CollectionMaps->collectionMap_Vertex["KinkVertices"] = v_cal ; - for(unsigned int i=0 ; i< KinkVertices->size(); i++) m_CollectionMaps->collectionMap_Vertex ["KinkVertices"].push_back(KinkVertices->at(i)); - } - if (NULL != ProngVertices ) - { - std::vector<edm4hep::Vertex> v_cal; - m_CollectionMaps->collectionMap_Vertex["ProngVertices"] = v_cal ; - for(unsigned int i=0 ; i< ProngVertices->size(); i++) m_CollectionMaps->collectionMap_Vertex ["ProngVertices"].push_back(ProngVertices->at(i)); - } - if (NULL != SplitVertices ) - { - std::vector<edm4hep::Vertex> v_cal; - m_CollectionMaps->collectionMap_Vertex["SplitVertices"] = v_cal ; - for(unsigned int i=0 ; i< SplitVertices->size(); i++) m_CollectionMaps->collectionMap_Vertex ["SplitVertices"].push_back(SplitVertices->at(i)); - } - if (NULL != V0Vertices ) - { - std::vector<edm4hep::Vertex> v_cal; - m_CollectionMaps->collectionMap_Vertex["V0Vertices"] = v_cal ; - for(unsigned int i=0 ; i< V0Vertices->size(); i++) m_CollectionMaps->collectionMap_Vertex ["V0Vertices"].push_back(V0Vertices->at(i)); - } - if (NULL != MarlinTrkTracks ) - { - std::vector<edm4hep::Track> v_cal; - m_CollectionMaps->collectionMap_Track["MarlinTrkTracks"] = v_cal ; - for(unsigned int i=0 ; i< MarlinTrkTracks->size(); i++) m_CollectionMaps->collectionMap_Track ["MarlinTrkTracks"].push_back(MarlinTrkTracks->at(i)); - } - if (NULL != mcRecoCaloAssociation ) - { - std::vector<edm4hep::MCRecoCaloAssociation> v_cal; - m_CollectionMaps->collectionMap_CaloRel["RecoCaloAssociation"] = v_cal ; - for(unsigned int i=0 ; i< mcRecoCaloAssociation->size(); i++) m_CollectionMaps->collectionMap_CaloRel ["RecoCaloAssociation"].push_back(mcRecoCaloAssociation->at(i)); - } - if (NULL != mcRecoTrackerAssociation ) - { - std::vector<edm4hep::MCRecoTrackerAssociation> v_cal; - m_CollectionMaps->collectionMap_TrkRel["RecoTrackerAssociation"] = v_cal ; - for(unsigned int i=0 ; i< mcRecoTrackerAssociation->size(); i++) m_CollectionMaps->collectionMap_TrkRel ["RecoTrackerAssociation"].push_back(mcRecoTrackerAssociation->at(i)); + for(auto &v : m_dataHandles){ + try{ + if(m_collections[v.first]=="MCParticle"){ + auto handle = dynamic_cast<DataHandle<edm4hep::MCParticleCollection>*> (v.second); + auto po = handle->get(); + if(po != NULL){ + std::vector<edm4hep::MCParticle> v_mc; + m_CollectionMaps->collectionMap_MC [v.first] = v_mc; + for(unsigned int i=0 ; i< po->size(); i++) m_CollectionMaps->collectionMap_MC [v.first].push_back(po->at(i)); + std::cout<<"saved col name="<<v.first<<std::endl; + } + else{ + std::cout<<"don't find col name="<<v.first<<std::endl; + } + } + else if(m_collections[v.first]=="CalorimeterHit"){ + auto handle = dynamic_cast<DataHandle<edm4hep::CalorimeterHitCollection>*> (v.second); + auto po = handle->get(); + if(po != NULL){ + std::vector<edm4hep::CalorimeterHit> v_cal; + m_CollectionMaps->collectionMap_CaloHit[v.first] = v_cal ; + for(unsigned int i=0 ; i< po->size(); i++) m_CollectionMaps->collectionMap_CaloHit [v.first].push_back(po->at(i)); + std::cout<<"saved col name="<<v.first<<std::endl; + } + else{ + std::cout<<"don't find col name="<<v.first<<std::endl; + } + } + else if(m_collections[v.first]=="Track"){ + auto handle = dynamic_cast<DataHandle<edm4hep::TrackCollection>*> (v.second); + auto po = handle->get(); + if(po != NULL){ + std::vector<edm4hep::Track> v_cal; + m_CollectionMaps->collectionMap_Track[v.first] = v_cal ; + for(unsigned int i=0 ; i< po->size(); i++) m_CollectionMaps->collectionMap_Track [v.first].push_back(po->at(i)); + std::cout<<"saved col name="<<v.first<<std::endl; + } + else{ + std::cout<<"don't find col name="<<v.first<<std::endl; + } + } + else if(m_collections[v.first]=="Vertex"){ + auto handle = dynamic_cast<DataHandle<edm4hep::VertexCollection>*> (v.second); + auto po = handle->get(); + if(po != NULL){ + std::vector<edm4hep::Vertex> v_cal; + m_CollectionMaps->collectionMap_Vertex[v.first] = v_cal ; + for(unsigned int i=0 ; i< po->size(); i++) m_CollectionMaps->collectionMap_Vertex [v.first].push_back(po->at(i)); + std::cout<<"saved col name="<<v.first<<std::endl; + } + else{ + std::cout<<"don't find col name="<<v.first<<std::endl; + } + } + else if(m_collections[v.first]=="MCRecoCaloAssociation"){ + auto handle = dynamic_cast<DataHandle<edm4hep::MCRecoCaloAssociationCollection>*> (v.second); + auto po = handle->get(); + if(po != NULL){ + std::vector<edm4hep::MCRecoCaloAssociation> v_cal; + m_CollectionMaps->collectionMap_CaloRel[v.first] = v_cal ; + for(unsigned int i=0 ; i< po->size(); i++) m_CollectionMaps->collectionMap_CaloRel [v.first].push_back(po->at(i)); + std::cout<<"saved col name="<<v.first<<std::endl; + } + else{ + std::cout<<"don't find col name="<<v.first<<std::endl; + } + } + else if(m_collections[v.first]=="MCRecoTrackerAssociation"){ + auto handle = dynamic_cast<DataHandle<edm4hep::MCRecoTrackerAssociationCollection>*> (v.second); + auto po = handle->get(); + if(po != NULL){ + std::vector<edm4hep::MCRecoTrackerAssociation> v_cal; + m_CollectionMaps->collectionMap_TrkRel[v.first] = v_cal ; + for(unsigned int i=0 ; i< po->size(); i++) m_CollectionMaps->collectionMap_TrkRel [v.first].push_back(po->at(i)); + std::cout<<"saved col name="<<v.first<<std::endl; + } + else{ + std::cout<<"don't find col name="<<v.first<<std::endl; + } + } + else{ + std::cout<<"wrong type name for col :"<<v.first<<std::endl; + } + }//try + catch(...){ + std::cout<<"don't find "<<v.first<<"in event"<<std::endl; + std::cout<<"don't find col name="<<v.first<<",with type="<<m_collections[v.first]<<" in this event"<<std::endl; } + } return StatusCode::SUCCESS; } - StatusCode PandoraPFAlg::Ana() { int n_current = m_tree->GetEntries()+1; diff --git a/Reconstruction/PFA/Pandora/GaudiPandora/src/TrackCreator.cpp b/Reconstruction/PFA/Pandora/GaudiPandora/src/TrackCreator.cpp index 8dc23c5725546cdb3626cbe3b4c1281648e4da8a..c3b83199123447b624bbafbe38ec0db9d2561080 100644 --- a/Reconstruction/PFA/Pandora/GaudiPandora/src/TrackCreator.cpp +++ b/Reconstruction/PFA/Pandora/GaudiPandora/src/TrackCreator.cpp @@ -138,69 +138,48 @@ TrackCreator::~TrackCreator() //------------------------------------------------------------------------------------------------------------------------------------------ -pandora::StatusCode TrackCreator::CreateTrackAssociations(const CollectionMaps& collectionMaps)//TODO, due to the lack of full Vertex information +pandora::StatusCode TrackCreator::CreateTrackAssociations(const CollectionMaps& collectionMaps) { -// TODO, because the Vertex.getAssociatedParticle() doesn't work for LCIO to plcio transfer -// PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, this->ExtractKinks(collectionMaps)); -// PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, this->ExtractProngsAndSplits(const CollectionMaps& collectionMaps)); -// PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, this->ExtractV0s(const CollectionMaps& collectionMaps)); + PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, this->ExtractKinks(collectionMaps)); + PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, this->ExtractProngsAndSplits(collectionMaps)); + PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, this->ExtractV0s(collectionMaps)); return pandora::STATUS_CODE_SUCCESS; } //------------------------------------------------------------------------------------------------------------------------------------------ -/* pandora::StatusCode TrackCreator::ExtractKinks(const CollectionMaps& collectionMaps) { std::cout<<"start TrackCreator::ExtractKinks:"<<std::endl; for (StringVector::const_iterator iter = m_settings.m_kinkVertexCollections.begin(), iterEnd = m_settings.m_kinkVertexCollections.end(); iter != iterEnd; ++iter) { - if(collectionMaps.CollectionMap_Vertex.find(*iter) == collectionMaps.CollectionMap_Vertex.end()) { std::cout<<"not find "<<(*iter)<<std::endl; continue;} + if(collectionMaps.collectionMap_Vertex.find(*iter) == collectionMaps.collectionMap_Vertex.end()) { std::cout<<"not find "<<(*iter)<<std::endl; continue;} try { - const edm4hep::VertexCollection *pKinkCollection = (collectionMaps.CollectionMap_Vertex.find(*iter))->second; + const std::vector<edm4hep::Vertex>& pKinkCollection = (collectionMaps.collectionMap_Vertex.find(*iter))->second; - for (int i = 0, iMax = pKinkCollection->size(); i < iMax; ++i) + for (int i = 0, iMax = pKinkCollection.size(); i < iMax; ++i) { try { - const edm4hep::Vertex pVertex0 = pKinkCollection->at(i); + const edm4hep::Vertex pVertex0 = pKinkCollection.at(i); const edm4hep::Vertex* pVertex = &(pVertex0); if (NULL == pVertex) throw ("Collection type mismatch"); - std::cout<<"pVertex0 getChi2="<<pVertex0.getChi2()<<std::endl; - std::cout<<"pVertex getChi2="<<pVertex->getChi2()<<std::endl; - std::cout<<"Hi 0 "<<std::endl; - const plcio::ConstReconstructedParticle pReconstructedParticle0 = pVertex0.getAssociatedParticle(); - std::cout<<"Hi 1"<<std::endl; - //EVENT::ReconstructedParticle *pReconstructedParticle = pVertex->getAssociatedParticle(); - const plcio::ConstReconstructedParticle pReconstructedParticle = pVertex->getAssociatedParticle(); - std::cout<<"Hi 2:"<<&pReconstructedParticle<<std::endl; - //const EVENT::TrackVec &trackVec(pReconstructedParticle->getTracks()); - //plcio::ConstTrack trackVec = pReconstructedParticle.getTracks(); - - std::cout<<"pReconstructedParticle0 en="<<pReconstructedParticle0.getEnergy()<<std::endl; - std::cout<<"pReconstructedParticle en="<<pReconstructedParticle.getEnergy()<<std::endl; - //if (this->IsConflictingRelationship(trackVec))continue; + //std::cout<<"pVertex getChi2="<<pVertex->getChi2()<<std::endl; + //std::cout<<"pReconstructedParticle en="<<pReconstructedParticle.getEnergy()<<",type="<<pReconstructedParticle.getType()<<std::endl; + const edm4hep::ConstReconstructedParticle pReconstructedParticle = pVertex->getAssociatedParticle(); if (this->IsConflictingRelationship(pReconstructedParticle))continue; - //const int vertexPdgCode(pReconstructedParticle->getType()); const int vertexPdgCode(pReconstructedParticle.getType()); // Extract the kink vertex information - //for (unsigned int iTrack = 0, nTracks = trackVec.size(); iTrack < nTracks; ++iTrack) - //for (unsigned int iTrack = 0, nTracks = trackVec.tracks_size(); iTrack < nTracks; ++iTrack) for (unsigned int iTrack = 0, nTracks = pReconstructedParticle.tracks_size(); iTrack < nTracks; ++iTrack) { - //EVENT::Track *pTrack = trackVec[iTrack]; - //plcio::ConstTrack pTrack = trackVec.getTracks(iTrack); - plcio::ConstTrack pTrack = pReconstructedParticle.getTracks(iTrack); - //(0 == iTrack) ? m_parentTrackList.insert(pTrack) : m_daughterTrackList.insert(pTrack); + edm4hep::ConstTrack pTrack = pReconstructedParticle.getTracks(iTrack); (0 == iTrack) ? m_parentTrackList.insert(pTrack.id()) : m_daughterTrackList.insert(pTrack.id()); - //std::cout << "KinkTrack " << iTrack << ", nHits " << pTrack.getTrackerHits().size() << std::endl; - std::cout << "KinkTrack " << iTrack << ", nHits " << pTrack.trackerHits_size() << std::endl; int trackPdgCode = pandora::UNKNOWN_PARTICLE_TYPE; @@ -229,7 +208,6 @@ pandora::StatusCode TrackCreator::ExtractKinks(const CollectionMaps& collectionM trackPdgCode = pandora::PI_PLUS; break; default : - //(pTrack->getOmega() > 0) ? trackPdgCode = pandora::PI_PLUS : trackPdgCode = pandora::PI_MINUS; (pTrack.getTrackStates(0).omega > 0) ? trackPdgCode = pandora::PI_PLUS : trackPdgCode = pandora::PI_MINUS; break; } @@ -245,9 +223,7 @@ pandora::StatusCode TrackCreator::ExtractKinks(const CollectionMaps& collectionM { for (unsigned int jTrack = iTrack + 1; jTrack < nTracks; ++jTrack) { - //PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackParentDaughterRelationship(*m_pPandora, pTrack, trackVec[jTrack])); - //PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackParentDaughterRelationship(*m_pPandora, (int*)(pTrack.id()), (int*)(trackVec.getTracks(jTrack).id()) ) ); - PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackParentDaughterRelationship(*m_pPandora, (int*)(pTrack.id()), (int*)(pReconstructedParticle.getTracks(jTrack).id()) ) ); + PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackParentDaughterRelationship(*m_pPandora, GetTrackAddress(collectionMaps, pTrack), GetTrackAddress(collectionMaps, pReconstructedParticle.getTracks(jTrack) ) ) ); } } @@ -256,9 +232,7 @@ pandora::StatusCode TrackCreator::ExtractKinks(const CollectionMaps& collectionM { for (unsigned int jTrack = iTrack + 1; jTrack < nTracks; ++jTrack) { - //PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackSiblingRelationship(*m_pPandora, pTrack, trackVec[jTrack])); - //PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackSiblingRelationship(*m_pPandora, (int*)(pTrack.id()), (int*)(trackVec.getTracks(jTrack).id()) )); - PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackSiblingRelationship(*m_pPandora, (int*)(pTrack.id()), (int*)(pReconstructedParticle.getTracks(jTrack).id()) )); + PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackSiblingRelationship(*m_pPandora, GetTrackAddress(collectionMaps, pTrack), GetTrackAddress(collectionMaps, pReconstructedParticle.getTracks(jTrack) ) ) ); } } } @@ -277,50 +251,44 @@ pandora::StatusCode TrackCreator::ExtractKinks(const CollectionMaps& collectionM return pandora::STATUS_CODE_SUCCESS; } -*/ //------------------------------------------------------------------------------------------------------------------------------------------ -/* -pandora::StatusCode TrackCreator::ExtractProngsAndSplits(const EVENT::LCEvent *const pLCEvent) + +pandora::StatusCode TrackCreator::ExtractProngsAndSplits(const CollectionMaps& collectionMaps) { - for (StringVector::const_iterator iter = m_settings.m_prongSplitVertexCollections.begin(), iterEnd = m_settings.m_prongSplitVertexCollections.end(); - iter != iterEnd; ++iter) + std::cout<<"start TrackCreator::ExtractProngsAndSplits:"<<std::endl; + for (StringVector::const_iterator iter = m_settings.m_prongSplitVertexCollections.begin(), iterEnd = m_settings.m_prongSplitVertexCollections.end(); iter != iterEnd; ++iter) { + if(collectionMaps.collectionMap_Vertex.find(*iter) == collectionMaps.collectionMap_Vertex.end()) { std::cout<<"not find "<<(*iter)<<std::endl; continue;} try { - const EVENT::LCCollection *pProngOrSplitCollection = pLCEvent->getCollection(*iter); + const std::vector<edm4hep::Vertex>& pProngOrSplitCollection = (collectionMaps.collectionMap_Vertex.find(*iter))->second; - for (int i = 0, iMax = pProngOrSplitCollection->getNumberOfElements(); i < iMax; ++i) + for (int i = 0, iMax = pProngOrSplitCollection.size(); i < iMax; ++i) { try { - EVENT::Vertex *pVertex = dynamic_cast<Vertex*>(pProngOrSplitCollection->getElementAt(i)); - - if (NULL == pVertex) - throw EVENT::Exception("Collection type mismatch"); + const edm4hep::Vertex pVertex0 = pProngOrSplitCollection.at(i); + const edm4hep::Vertex* pVertex = &(pVertex0); - EVENT::ReconstructedParticle *pReconstructedParticle = pVertex->getAssociatedParticle(); - const EVENT::TrackVec &trackVec(pReconstructedParticle->getTracks()); + if (NULL == pVertex) throw ("Collection type mismatch"); + const edm4hep::ConstReconstructedParticle pReconstructedParticle = pVertex->getAssociatedParticle(); - if (this->IsConflictingRelationship(trackVec)) - continue; + if (this->IsConflictingRelationship(pReconstructedParticle))continue; // Extract the prong/split vertex information - for (unsigned int iTrack = 0, nTracks = trackVec.size(); iTrack < nTracks; ++iTrack) + for (unsigned int iTrack = 0, nTracks = pReconstructedParticle.tracks_size(); iTrack < nTracks; ++iTrack) { - EVENT::Track *pTrack = trackVec[iTrack]; - (0 == iTrack) ? m_parentTrackList.insert(pTrack) : m_daughterTrackList.insert(pTrack); - streamlog_out(DEBUG) << "Prong or Split Track " << iTrack << ", nHits " << pTrack->getTrackerHits().size() << std::endl; + edm4hep::ConstTrack pTrack = pReconstructedParticle.getTracks(iTrack); + (0 == iTrack) ? m_parentTrackList.insert(pTrack.id()) : m_daughterTrackList.insert(pTrack.id()); - if (0 == m_settings.m_shouldFormTrackRelationships) - continue; + if (0 == m_settings.m_shouldFormTrackRelationships) continue; // Make track parent-daughter relationships if (0 == iTrack) { for (unsigned int jTrack = iTrack + 1; jTrack < nTracks; ++jTrack) { - PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackParentDaughterRelationship(*m_pPandora, - pTrack, trackVec[jTrack])); + PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackParentDaughterRelationship(*m_pPandora, GetTrackAddress(collectionMaps, pTrack), GetTrackAddress(collectionMaps, pReconstructedParticle.getTracks(jTrack) ) ) ); } } @@ -329,21 +297,20 @@ pandora::StatusCode TrackCreator::ExtractProngsAndSplits(const EVENT::LCEvent *c { for (unsigned int jTrack = iTrack + 1; jTrack < nTracks; ++jTrack) { - PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackSiblingRelationship(*m_pPandora, - pTrack, trackVec[jTrack])); + PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackSiblingRelationship(*m_pPandora, GetTrackAddress(collectionMaps, pTrack), GetTrackAddress(collectionMaps, pReconstructedParticle.getTracks(jTrack) ) ) ); } } } } - catch (EVENT::Exception &exception) + catch (...) { - streamlog_out(WARNING) << "Failed to extract prong/split vertex: " << exception.what() << std::endl; + std::cout << "Failed to extract prong/split vertex" <<std::endl; } } } - catch (EVENT::Exception &exception) + catch (...) { - streamlog_out(DEBUG5) << "Failed to extract prong/split vertex collection: " << *iter << ", " << exception.what() << std::endl; + std::cout<< "Failed to extract prong/split vertex collection: " << *iter << std::endl; } } @@ -352,88 +319,84 @@ pandora::StatusCode TrackCreator::ExtractProngsAndSplits(const EVENT::LCEvent *c //------------------------------------------------------------------------------------------------------------------------------------------ -pandora::StatusCode TrackCreator::ExtractV0s(const EVENT::LCEvent *const pLCEvent) +pandora::StatusCode TrackCreator::ExtractV0s(const CollectionMaps& collectionMaps) { - for (StringVector::const_iterator iter = m_settings.m_v0VertexCollections.begin(), iterEnd = m_settings.m_v0VertexCollections.end(); - iter != iterEnd; ++iter) + std::cout<<"start TrackCreator::ExtractV0s:"<<std::endl; + for (StringVector::const_iterator iter = m_settings.m_v0VertexCollections.begin(), iterEnd = m_settings.m_v0VertexCollections.end(); iter != iterEnd; ++iter) { + if(collectionMaps.collectionMap_Vertex.find(*iter) == collectionMaps.collectionMap_Vertex.end()) { std::cout<<"not find "<<(*iter)<<std::endl; continue;} try { - const EVENT::LCCollection *pV0Collection = pLCEvent->getCollection(*iter); + const std::vector<edm4hep::Vertex>& pV0Collection = (collectionMaps.collectionMap_Vertex.find(*iter))->second; - for (int i = 0, iMax = pV0Collection->getNumberOfElements(); i < iMax; ++i) + for (int i = 0, iMax = pV0Collection.size(); i < iMax; ++i) { try { - EVENT::Vertex *pVertex = dynamic_cast<Vertex*>(pV0Collection->getElementAt(i)); + const edm4hep::Vertex pVertex0 = pV0Collection.at(i); + const edm4hep::Vertex* pVertex = &(pVertex0); - if (NULL == pVertex) - throw EVENT::Exception("Collection type mismatch"); + if (NULL == pVertex) throw ("Collection type mismatch"); - EVENT::ReconstructedParticle *pReconstructedParticle = pVertex->getAssociatedParticle(); - const EVENT::TrackVec &trackVec(pReconstructedParticle->getTracks()); + const edm4hep::ConstReconstructedParticle pReconstructedParticle = pVertex->getAssociatedParticle(); - if (this->IsConflictingRelationship(trackVec)) - continue; + if (this->IsConflictingRelationship(pReconstructedParticle))continue; // Extract the v0 vertex information - const int vertexPdgCode(pReconstructedParticle->getType()); + const int vertexPdgCode(pReconstructedParticle.getType()); - for (unsigned int iTrack = 0, nTracks = trackVec.size(); iTrack < nTracks; ++iTrack) + for (unsigned int iTrack = 0, nTracks = pReconstructedParticle.tracks_size(); iTrack < nTracks; ++iTrack) { - EVENT::Track *pTrack = trackVec[iTrack]; - m_v0TrackList.insert(pTrack); - streamlog_out(DEBUG) << "V0Track " << iTrack << ", nHits " << pTrack->getTrackerHits().size() << std::endl; + edm4hep::ConstTrack pTrack = pReconstructedParticle.getTracks(iTrack); + m_v0TrackList.insert(pTrack.id()); int trackPdgCode = pandora::UNKNOWN_PARTICLE_TYPE; switch (vertexPdgCode) { case pandora::PHOTON : - (pTrack->getOmega() > 0) ? trackPdgCode = pandora::E_PLUS : trackPdgCode = pandora::E_MINUS; + (pTrack.getTrackStates(0).omega > 0) ? trackPdgCode = pandora::E_PLUS : trackPdgCode = pandora::E_MINUS; break; case pandora::LAMBDA : - (pTrack->getOmega() > 0) ? trackPdgCode = pandora::PROTON : trackPdgCode = pandora::PI_MINUS; + (pTrack.getTrackStates(0).omega > 0) ? trackPdgCode = pandora::PROTON : trackPdgCode = pandora::PI_MINUS; break; case pandora::LAMBDA_BAR : - (pTrack->getOmega() > 0) ? trackPdgCode = pandora::PI_PLUS : trackPdgCode = pandora::PROTON_BAR; + (pTrack.getTrackStates(0).omega > 0) ? trackPdgCode = pandora::PI_PLUS : trackPdgCode = pandora::PROTON_BAR; break; case pandora::K_SHORT : - (pTrack->getOmega() > 0) ? trackPdgCode = pandora::PI_PLUS : trackPdgCode = pandora::PI_MINUS; + (pTrack.getTrackStates(0).omega > 0) ? trackPdgCode = pandora::PI_PLUS : trackPdgCode = pandora::PI_MINUS; break; default : - (pTrack->getOmega() > 0) ? trackPdgCode = pandora::PI_PLUS : trackPdgCode = pandora::PI_MINUS; + (pTrack.getTrackStates(0).omega > 0) ? trackPdgCode = pandora::PI_PLUS : trackPdgCode = pandora::PI_MINUS; break; } m_trackToPidMap.insert(TrackToPidMap::value_type(pTrack, trackPdgCode)); - if (0 == m_settings.m_shouldFormTrackRelationships) - continue; + if (0 == m_settings.m_shouldFormTrackRelationships) continue; // Make track sibling relationships for (unsigned int jTrack = iTrack + 1; jTrack < nTracks; ++jTrack) { - PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackSiblingRelationship(*m_pPandora, - pTrack, trackVec[jTrack])); + PANDORA_RETURN_RESULT_IF(pandora::STATUS_CODE_SUCCESS, !=, PandoraApi::SetTrackSiblingRelationship(*m_pPandora, GetTrackAddress(collectionMaps, pTrack), GetTrackAddress(collectionMaps, pReconstructedParticle.getTracks(jTrack) ) ) ); } } } - catch (EVENT::Exception &exception) + catch (...) { - streamlog_out(WARNING) << "Failed to extract v0 vertex: " << exception.what() << std::endl; + std::cout<< "Failed to extract v0 vertex" << std::endl; } } } - catch (EVENT::Exception &exception) + catch (...) { - streamlog_out(DEBUG5) << "Failed to extract v0 vertex collection: " << *iter << ", " << exception.what() << std::endl; + std::cout << "Failed to extract v0 vertex collection: " << *iter << std::endl; } } return pandora::STATUS_CODE_SUCCESS; } -*/ + //------------------------------------------------------------------------------------------------------------------------------------------ bool TrackCreator::IsConflictingRelationship(const edm4hep::ConstReconstructedParticle &Particle) const @@ -450,6 +413,20 @@ bool TrackCreator::IsConflictingRelationship(const edm4hep::ConstReconstructedPa return false; } +const edm4hep::Track* TrackCreator::GetTrackAddress(const CollectionMaps& collectionMaps, const edm4hep::ConstTrack& pTrack ) +{ + for (StringVector::const_iterator iter = m_settings.m_trackCollections.begin(), iterEnd = m_settings.m_trackCollections.end(); iter != iterEnd; ++iter) + { + if(collectionMaps.collectionMap_Track.find(*iter) == collectionMaps.collectionMap_Track.end()) { std::cout<<"not find "<<(*iter)<<std::endl; continue;} + const std::vector<edm4hep::Track>& pTrackCollection = (collectionMaps.collectionMap_Track.find(*iter))->second; + for (int i = 0, iMax = pTrackCollection.size(); i < iMax; ++i) + { + const edm4hep::Track& pTrack0 = pTrackCollection.at(i); + if (pTrack.id() == pTrack0.id()) return (const edm4hep::Track*) (&pTrack0); + } + } + return NULL; +} //------------------------------------------------------------------------------------------------------------------------------------------ pandora::StatusCode TrackCreator::CreateTracks(const CollectionMaps& collectionMaps) @@ -493,8 +470,7 @@ pandora::StatusCode TrackCreator::CreateTracks(const CollectionMaps& collectionM const int nTrackHits(static_cast<int>(pTrack->trackerHits_size())); - if ((nTrackHits < minTrackHits) || (nTrackHits > m_settings.m_maxTrackHits)) - continue; + if ((nTrackHits < minTrackHits) || (nTrackHits > m_settings.m_maxTrackHits)) continue; // Proceed to create the pandora track PandoraApi::Track::Parameters trackParameters; @@ -558,7 +534,8 @@ void TrackCreator::GetTrackStates(const edm4hep::Track *const pTrack, PandoraApi this->CopyTrackState(pEndTrack.getTrackStates(3), trackParameters.m_trackStateAtEnd); //FIXME ? LCIO input only has 4 states, so 4 can't be used. - this->CopyTrackState(pEndTrack.getTrackStates(3), trackParameters.m_trackStateAtCalorimeter); + if( pEndTrack.trackStates_size()<5) this->CopyTrackState(pEndTrack.getTrackStates(3), trackParameters.m_trackStateAtCalorimeter); + else this->CopyTrackState(pEndTrack.getTrackStates(4), trackParameters.m_trackStateAtCalorimeter); trackParameters.m_isProjectedToEndCap = ((std::fabs(trackParameters.m_trackStateAtCalorimeter.Get().GetPosition().GetZ()) < m_eCalEndCapInnerZ) ? false : true); diff --git a/Reconstruction/PFA/Pandora/PandoraLikelihoodData9EBin.xml b/Reconstruction/PFA/Pandora/PandoraLikelihoodData9EBin.xml new file mode 100644 index 0000000000000000000000000000000000000000..ab63905438a22a636fc7d58a6853acc4f858533e --- /dev/null +++ b/Reconstruction/PFA/Pandora/PandoraLikelihoodData9EBin.xml @@ -0,0 +1,652 @@ +<NEnergyBins>9</NEnergyBins> +<EnergyBinLowerEdges>0 0.2 0.5 1 1.5 2.5 5 10 20 </EnergyBinLowerEdges> +<NSignalEvents>42213 22561 20497 14219 19767 26770 22039 14547 10186 </NSignalEvents> +<NBackgroundEvents>109034 45109 13323 4171 3447 3033 2077 1206 1006 </NBackgroundEvents> +<PhotonSigPeakRms_0> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0.000213204 0.00485632 0.0255372 0.0994717 0.15481 0.0870111 0.0978845 0.0800701 0.0638429 0.0593656 0.0497714 0.0426172 0.0375951 0.0313884 0.0267216 0.0234999 0.0200649 0.0174117 0.0140478 0.0110393 0.00916779 0.00864663 0.00575652 0.00547225 0.00381399 0.0035771 0.0028901 0.00232156 0.00198991 0.00172933 0.00127923 0.00116078 0.000947575 0.000639613 0.000331651 0.000473788 0.000426409 0.000355341 0.000331651 0.000260583 0.000260583 0.000260583 0.000118447 7.10682e-05 0.000331651 4.73788e-05 4.73788e-05 7.10682e-05 0 </BinContents> +</PhotonSigPeakRms_0> +<PhotonBkgPeakRms_0> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0.000183429 0.00168755 0.0126291 0.0673093 0.147431 0.0799659 0.109259 0.0951355 0.0711154 0.0654383 0.0544968 0.0471596 0.0394831 0.0327054 0.0263221 0.0227177 0.0189482 0.0167562 0.0134729 0.011446 0.00903388 0.00792413 0.00656676 0.0058147 0.00423721 0.00423721 0.00329255 0.00319167 0.00265972 0.00228369 0.0020544 0.00180678 0.00159583 0.00145826 0.00132986 0.00114643 0.00113726 0.000953831 0.000953831 0.000935488 0.000834602 0.000596144 0.000412715 0.000412715 0.000403544 0.000421887 0.000247629 0.000394372 0 </BinContents> +</PhotonBkgPeakRms_0> +<PhotonSigPeakRms_1> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0.000132973 0.00093081 0.00483135 0.0188378 0.0378086 0.0551837 0.0777448 0.1023 0.108949 0.1066 0.0948096 0.0789859 0.0641816 0.0491556 0.0402908 0.0274811 0.0224724 0.0156021 0.0133859 0.0121892 0.00979566 0.00771242 0.00611675 0.00540756 0.00483135 0.0047427 0.00367891 0.00305838 0.00257081 0.00248216 0.00164 0.00234919 0.00195027 0.00203892 0.00132973 0.00115243 0.00101946 0.00093081 0.0012854 0.000531891 0.000398918 0.000398918 0.000664864 0.000398918 0.000576216 0.000487567 0.00031027 0.000265946 0 </BinContents> +</PhotonSigPeakRms_1> +<PhotonBkgPeakRms_1> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 4.43371e-05 0.000665056 0.00181782 0.00589683 0.0195526 0.0389058 0.0514975 0.0672815 0.0762376 0.0795185 0.0771686 0.0732448 0.0666164 0.0605644 0.0513423 0.0425414 0.0407236 0.0338292 0.0288856 0.0250948 0.0215921 0.0178457 0.0153628 0.0141435 0.0114833 0.00935512 0.00791416 0.00773682 0.00649538 0.00476623 0.00418985 0.00474406 0.00399033 0.00336962 0.00272673 0.00270456 0.00266022 0.00195083 0.00192866 0.00172915 0.00157397 0.00175131 0.00164047 0.00168481 0.00104192 0.00113059 0.0011971 0.00115276 0.000709393 0 </BinContents> +</PhotonBkgPeakRms_1> +<PhotonSigPeakRms_2> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 4.87876e-05 0 0.000487876 0.00278089 0.00482998 0.0105869 0.0228814 0.0499097 0.0847929 0.110894 0.131678 0.136605 0.115627 0.0926477 0.0687906 0.0494707 0.0325901 0.021369 0.0157096 0.0120018 0.00785481 0.00658633 0.00390301 0.00336635 0.00219544 0.00195151 0.0015612 0.00160999 0.000926965 0.000926965 0.000390301 0.000878177 0.000634239 0.000390301 0.000487876 0.000439089 0.000243938 0.000292726 0.000243938 4.87876e-05 0.000243938 0.000292726 9.75753e-05 0.000146363 0.000146363 0.000146363 4.87876e-05 0.000243938 0 </BinContents> +</PhotonSigPeakRms_2> +<PhotonBkgPeakRms_2> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0.000225175 0.000600465 0.00150116 0.00360279 0.011559 0.0228177 0.02582 0.0390302 0.0531412 0.0672521 0.0832395 0.0880432 0.0803873 0.0737071 0.0711551 0.0585454 0.0469864 0.0385799 0.0316745 0.029648 0.0231179 0.0178638 0.0179389 0.0138107 0.0114088 0.0101329 0.00900698 0.0082564 0.00683029 0.00637994 0.00547925 0.00450349 0.00382797 0.00397808 0.00390302 0.00330256 0.00240186 0.00307738 0.00225174 0.0023268 0.00172634 0.00217669 0.00150116 0.00120093 0.00112587 0.00142611 0.00127599 0.00120093 0.00105081 0 </BinContents> +</PhotonBkgPeakRms_2> +<PhotonSigPeakRms_3> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0 0.000562627 0.00105493 0.00281314 0.00464168 0.00829875 0.0132217 0.0302412 0.0673043 0.109923 0.138828 0.147338 0.131584 0.100992 0.0720866 0.0485969 0.0335467 0.0244743 0.0164569 0.0103383 0.00689219 0.0056966 0.00506365 0.00344609 0.00295379 0.00225051 0.00168788 0.0014769 0.00126591 0.00091427 0.00105493 0.000773613 0.000632956 0.000562627 0.000492299 0.000492299 0.000492299 0.000281314 0.000140657 7.03284e-05 0.000140657 7.03284e-05 0.000140657 0.000281314 0.000281314 7.03284e-05 7.03284e-05 0 </BinContents> +</PhotonSigPeakRms_3> +<PhotonBkgPeakRms_3> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0.000959003 0.000959003 0.00119875 0.00263726 0.00455526 0.0129465 0.0148645 0.0275713 0.0357228 0.0417166 0.0592184 0.0748022 0.0863102 0.0949413 0.0937425 0.0791177 0.0678494 0.0537042 0.0431551 0.0364421 0.0198993 0.0218173 0.0177415 0.012467 0.0136658 0.00982978 0.00671302 0.00767202 0.00479501 0.00599377 0.00551427 0.00407576 0.00287701 0.00359626 0.00311676 0.00287701 0.00383601 0.00263726 0.00191801 0.00167825 0.00287701 0.00119875 0.00167825 0.00167825 0.0014385 0.00167825 0.00119875 0.00119875 0.00191801 0 </BinContents> +</PhotonBkgPeakRms_3> +<PhotonSigPeakRms_4> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0.000202357 0.000151768 0.00126473 0.00222593 0.00338949 0.00419892 0.00607072 0.0122932 0.0264582 0.0617696 0.105934 0.1459 0.14676 0.137047 0.103455 0.0766429 0.0498305 0.0324278 0.0233217 0.0125968 0.0107249 0.00698133 0.00490717 0.00480599 0.00344008 0.00288359 0.00212475 0.00156827 0.00106238 0.00131532 0.00080943 0.00101179 0.000708251 0.000910609 0.00075884 0.000556483 0.000657662 0.000404715 0.000354126 0.000455304 0.000202357 0.000354126 0.000151768 0.000252947 0.000202357 0.000252947 0.000202357 0 </BinContents> +</PhotonSigPeakRms_4> +<PhotonBkgPeakRms_4> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0.000580215 0.000870322 0.00116043 0.00203075 0.00319118 0.00696258 0.0127647 0.018857 0.0208877 0.0345228 0.0368436 0.05396 0.0791993 0.0771685 0.100667 0.086452 0.0786191 0.0646939 0.0519292 0.0472875 0.0301712 0.0281404 0.0258196 0.0229185 0.0142153 0.00870322 0.00899333 0.0078329 0.00754279 0.00696258 0.00870322 0.00638236 0.00551204 0.00290107 0.00464172 0.00435161 0.00319118 0.00261097 0.00261097 0.00145054 0.00290107 0.00319118 0.00145054 0.00232086 0.00290107 0.00203075 0.000580215 0.00203075 0.000290107 0 </BinContents> +</PhotonBkgPeakRms_4> +<PhotonSigPeakRms_5> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0 7.47105e-05 0.000747105 0.00194247 0.00265222 0.00414643 0.00489354 0.00638775 0.00762047 0.0194247 0.0411655 0.0872619 0.130855 0.147217 0.140082 0.115652 0.0861038 0.0577512 0.0395966 0.0255136 0.0187897 0.0132238 0.0100486 0.0075084 0.00649981 0.00377288 0.00373552 0.00272693 0.00212925 0.00164363 0.00164363 0.0014195 0.0010833 0.00100859 0.000971236 0.000560329 0.000933881 0.000635039 0.000373552 0.000448263 0.000485618 0.000410908 0.000298842 0.000112066 0.000261487 7.47105e-05 0.000112066 0 </BinContents> +</PhotonSigPeakRms_5> +<PhotonBkgPeakRms_5> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0.000329707 0.000329707 0.00131883 0.00098912 0.00395648 0.00758325 0.0102209 0.0181339 0.0214309 0.0276954 0.0412133 0.0626442 0.0791296 0.0919881 0.0900099 0.0824266 0.0787999 0.0646225 0.0586878 0.0412133 0.0306627 0.0253874 0.019123 0.0211012 0.0151665 0.0121991 0.00923178 0.00890208 0.00758325 0.00956149 0.00626442 0.00560501 0.0052753 0.00395648 0.0052753 0.00296736 0.00395648 0.00296736 0.00395648 0.00296736 0.00329707 0.00263765 0.00296736 0.00164853 0.00329707 0.000659413 0.000659413 0 </BinContents> +</PhotonBkgPeakRms_5> +<PhotonSigPeakRms_6> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0 4.53741e-05 0.000453741 0.0012251 0.00235945 0.00349381 0.00725986 0.0077136 0.00830346 0.0107083 0.0142475 0.030809 0.0653387 0.105949 0.14125 0.141113 0.126412 0.0980988 0.0671537 0.045601 0.0323517 0.02305 0.0158356 0.0112528 0.00925632 0.0064885 0.00480966 0.00426517 0.00299469 0.00190571 0.00163347 0.00167884 0.00154272 0.00131585 0.000408367 0.00077136 0.000499115 0.000408367 0.000544489 0.000226871 0.000317619 0.000317619 0.000181496 0.000181496 9.07482e-05 9.07482e-05 4.53741e-05 0 </BinContents> +</PhotonSigPeakRms_6> +<PhotonBkgPeakRms_6> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0 0 0.000481464 0.000962927 0.00385171 0.00674049 0.00818488 0.00770342 0.0129995 0.0187771 0.0187771 0.0438132 0.0370727 0.0597015 0.0746269 0.0702937 0.0794415 0.0871449 0.0645161 0.0654791 0.0452576 0.0409244 0.0365912 0.0284064 0.0327395 0.0202215 0.0231103 0.0216659 0.0144439 0.0129995 0.0110737 0.00914781 0.00433317 0.0052961 0.0052961 0.00481464 0.00337025 0.00288878 0.00337025 0.00144439 0.00192585 0.00337025 0.00337025 0 0.00144439 0.000962927 0.000962927 0 </BinContents> +</PhotonBkgPeakRms_6> +<PhotonSigPeakRms_7> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0.000481199 0.00288719 0.00481199 0.00749295 0.0103801 0.0105864 0.0139548 0.0204853 0.0338901 0.063587 0.10607 0.131986 0.147866 0.122225 0.0924589 0.0685365 0.0504571 0.0290094 0.021104 0.0145047 0.0129924 0.00866158 0.00680553 0.00501822 0.00281845 0.00247474 0.00213102 0.00130611 0.00116863 0.000962398 0.000824912 0.000618684 0.000343713 0.000206228 0.000274971 0.000206228 0.000137485 0.000137485 0 0 6.87427e-05 0 6.87427e-05 0 </BinContents> +</PhotonSigPeakRms_7> +<PhotonBkgPeakRms_7> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0 0 0.000829187 0 0 0.00165837 0.00165837 0.00746269 0.00912106 0.0066335 0.0107794 0.0182421 0.0381426 0.0323383 0.0489221 0.0605307 0.0588723 0.0746269 0.0704809 0.0613599 0.0538972 0.0480929 0.0505804 0.0613599 0.053068 0.0414594 0.0290216 0.0273632 0.0298507 0.0215589 0.0182421 0.00829187 0.0107794 0.00912106 0.0107794 0.0066335 0.00165837 0.00497512 0.00248756 0.00414594 0.00248756 0.000829187 0 0 0.000829187 0 0.000829187 0 </BinContents> +</PhotonBkgPeakRms_7> +<PhotonSigPeakRms_8> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0 9.8174e-05 0.000883566 0.00137444 0.00274887 0.00382878 0.00854114 0.0105046 0.0196348 0.0350481 0.0608679 0.0949342 0.118692 0.125663 0.126743 0.112704 0.0869821 0.060966 0.0416258 0.0243471 0.0175731 0.0121736 0.009032 0.00628313 0.00549774 0.00343609 0.00274887 0.00215983 0.00137444 0.000687218 0.000687218 0.000785392 0.000589044 0.000392696 0.000196348 9.8174e-05 0 0 9.8174e-05 0 0 0 0 0 </BinContents> +</PhotonSigPeakRms_8> +<PhotonBkgPeakRms_8> + <NBinsX>50</NBinsX> + <XLow>0</XLow> + <XHigh>5</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0 0 0 0.000994036 0.000994036 0.00198807 0.00298211 0.00695825 0.0109344 0.0119284 0.0228628 0.0298211 0.0407555 0.05666 0.0685885 0.0854871 0.0755467 0.0805169 0.0785288 0.0675944 0.0606362 0.0536779 0.0447316 0.0377734 0.0357853 0.028827 0.026839 0.0208747 0.0139165 0.00894632 0.00397614 0.00397614 0.00795229 0.00397614 0.00198807 0.000994036 0.000994036 0 0 0.000994036 0 0 0 0 </BinContents> +</PhotonBkgPeakRms_8> +<PhotonSigRmsRatio_0> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.122877 0.0934072 0.085874 0.0783408 0.0719447 0.0632743 0.0581811 0.0520219 0.045223 0.04103 0.0363632 0.0317438 0.0279772 0.0241158 0.021581 0.0195201 0.0181698 0.0155639 0.0136688 0.0120105 0.00952313 0.00933362 0.00840973 0.00774643 0.00696468 0.00559069 0.00570914 0.00495108 0.00483263 0.00405088 0 </BinContents> +</PhotonSigRmsRatio_0> +<PhotonBkgRmsRatio_0> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.128409 0.085056 0.0776363 0.0704001 0.0624851 0.0577893 0.052103 0.0482327 0.0486545 0.0454445 0.0397216 0.0338885 0.0294954 0.0266064 0.0231121 0.0213328 0.0191867 0.0162243 0.015087 0.0132986 0.0117761 0.0114185 0.0105105 0.00974925 0.00872205 0.00805253 0.00725462 0.0067135 0.00589724 0.00574133 0 </BinContents> +</PhotonBkgRmsRatio_0> +<PhotonSigRmsRatio_1> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.122512 0.112273 0.100483 0.0878064 0.0772129 0.0699437 0.0597048 0.0542529 0.0468064 0.0381189 0.0327556 0.0286335 0.0246886 0.0201675 0.019813 0.0155135 0.0132973 0.0117903 0.0107265 0.00899783 0.00802269 0.00589513 0.00567351 0.00465405 0.00429945 0.00421081 0.00367891 0.00279243 0.00314702 0.00212757 0 </BinContents> +</PhotonSigRmsRatio_1> +<PhotonBkgRmsRatio_1> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.11 0.0959232 0.0860139 0.0770135 0.0702077 0.0618724 0.0554657 0.0507438 0.0441154 0.0419428 0.0352923 0.0298832 0.0284644 0.024075 0.0227006 0.0205724 0.0195083 0.0161165 0.0146977 0.0127469 0.0120818 0.0100423 0.0101975 0.00966548 0.00804717 0.00755947 0.00720477 0.00620719 0.00607418 0.0055643 0 </BinContents> +</PhotonBkgRmsRatio_1> +<PhotonSigRmsRatio_2> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.157535 0.138313 0.121725 0.104893 0.0883056 0.0722545 0.0612773 0.0521052 0.0390301 0.0320535 0.0252232 0.0215153 0.0172708 0.013319 0.0115627 0.00756208 0.00751329 0.00648875 0.00502513 0.0035615 0.00224423 0.00224423 0.00180514 0.00200029 0.00126848 0.00107333 0.00112212 0.000926965 0.000390301 0.000390301 0 </BinContents> +</PhotonSigRmsRatio_2> +<PhotonBkgRmsRatio_2> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.124446 0.117316 0.0990768 0.0917961 0.0779104 0.0668768 0.0587705 0.0504391 0.0409818 0.0385799 0.030999 0.0288974 0.0220671 0.0210163 0.0153869 0.0156121 0.0135105 0.0128349 0.010358 0.0101329 0.0082564 0.00698041 0.00668018 0.006455 0.00495384 0.00517901 0.00457855 0.00292727 0.00367785 0.00330256 0 </BinContents> +</PhotonBkgRmsRatio_2> +<PhotonSigRmsRatio_3> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.181096 0.166327 0.131936 0.113369 0.0912863 0.0707504 0.0577396 0.0446586 0.0326324 0.0264435 0.0197623 0.0148393 0.0110416 0.00787678 0.00675153 0.00478233 0.00337576 0.00344609 0.00281314 0.00189887 0.00161755 0.00140657 0.000984598 0.000703284 0.000351642 0.000492299 0.000492299 0.000281314 0.000421971 0.000421971 0 </BinContents> +</PhotonSigRmsRatio_3> +<PhotonBkgRmsRatio_3> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.157996 0.135939 0.110285 0.107169 0.078878 0.0733637 0.0618557 0.0409974 0.0419564 0.0347638 0.0246943 0.0165428 0.0239751 0.0158235 0.011508 0.00982978 0.00863102 0.00551427 0.00839127 0.00455526 0.00503476 0.00287701 0.00407576 0.00263726 0.00335651 0.00191801 0.00311676 0.00191801 0.00119875 0.00119875 0 </BinContents> +</PhotonBkgRmsRatio_3> +<PhotonSigRmsRatio_4> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.203572 0.180199 0.155765 0.116659 0.090302 0.0672839 0.0489199 0.035969 0.0283806 0.0197804 0.0142662 0.00900491 0.00804371 0.00500835 0.0033389 0.00268124 0.00252947 0.00197299 0.00136591 0.00101179 0.000860019 0.000657662 0.000860019 0.000556483 0.000252947 0.000404715 5.05894e-05 0.000101179 5.05894e-05 0.000151768 0 </BinContents> +</PhotonSigRmsRatio_4> +<PhotonBkgRmsRatio_4> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.173194 0.158979 0.120104 0.111981 0.0844212 0.0690455 0.0586017 0.0420656 0.0348129 0.0249492 0.0203075 0.0182768 0.0124746 0.010734 0.010734 0.010734 0.00899333 0.00725268 0.0037714 0.0037714 0.00232086 0.00145054 0.00116043 0.00261097 0.00174064 0.00174064 0.000870322 0.000870322 0.000580215 0.00145054 0 </BinContents> +</PhotonBkgRmsRatio_4> +<PhotonSigRmsRatio_5> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.240605 0.202652 0.159357 0.118379 0.0868136 0.0602167 0.0386627 0.0288383 0.0196489 0.0131117 0.00833022 0.00620097 0.00437056 0.002839 0.00212925 0.00149421 0.00119537 0.00104595 0.000896526 0.000672394 0.000522973 0.000448263 0.000298842 0.000373552 0.000149421 7.47105e-05 0.000224131 0.000112066 0.000261487 7.47105e-05 0 </BinContents> +</PhotonSigRmsRatio_5> +<PhotonBkgRmsRatio_5> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.183976 0.16881 0.149027 0.107814 0.0877019 0.0652819 0.0530828 0.0375865 0.0323112 0.0240686 0.0154962 0.0161556 0.0128586 0.00923178 0.00791296 0.00428619 0.00329707 0.0052753 0.00131883 0.00197824 0.00329707 0.00131883 0.00164853 0.000329707 0.00131883 0.000659413 0.00098912 0.00131883 0.00098912 0.000659413 0 </BinContents> +</PhotonBkgRmsRatio_5> +<PhotonSigRmsRatio_6> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.269931 0.222242 0.163755 0.1182 0.0780435 0.0516357 0.0320795 0.0208267 0.0135215 0.00871183 0.00657925 0.00381143 0.00317619 0.00186034 0.00163347 0.000907482 0.000635238 0.000408367 0.000408367 0.000408367 0.000362993 0.000226871 9.07482e-05 9.07482e-05 0.000136122 9.07482e-05 4.53741e-05 4.53741e-05 4.53741e-05 9.07482e-05 0 </BinContents> +</PhotonSigRmsRatio_6> +<PhotonBkgRmsRatio_6> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.202215 0.189215 0.135773 0.129032 0.0832932 0.0683678 0.0476649 0.0389986 0.0250361 0.0202215 0.0178142 0.00962927 0.00818488 0.0052961 0.00433317 0.0052961 0.00337025 0.000481464 0.000962927 0.00144439 0.000962927 0 0.000962927 0.000962927 0.000481464 0 0 0 0 0 0 </BinContents> +</PhotonBkgRmsRatio_6> +<PhotonSigRmsRatio_7> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.291538 0.231044 0.16842 0.116244 0.0696364 0.0448202 0.0283907 0.0171169 0.0116863 0.00673678 0.00474325 0.00302468 0.00240599 0.000962398 0.000687427 0.000549942 0.000687427 0.000343713 0.000137485 0.000206228 0 6.87427e-05 6.87427e-05 6.87427e-05 0.000137485 6.87427e-05 6.87427e-05 0 6.87427e-05 6.87427e-05 0 </BinContents> +</PhotonSigRmsRatio_7> +<PhotonBkgRmsRatio_7> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.237977 0.217247 0.1534 0.108624 0.0713101 0.0704809 0.0422886 0.0315091 0.0190713 0.0157546 0.00829187 0.00580431 0.00497512 0.00165837 0.00165837 0.00331675 0 0.00248756 0.00165837 0.000829187 0 0 0.000829187 0 0 0.000829187 0 0 0 0 0 </BinContents> +</PhotonBkgRmsRatio_7> +<PhotonSigRmsRatio_8> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.311604 0.245042 0.174357 0.106421 0.0615551 0.0387787 0.0229727 0.0148243 0.00893383 0.00530139 0.00353426 0.002258 0.00147261 0.00098174 0.00049087 0.000196348 0.000294522 0.000392696 9.8174e-05 0.000294522 9.8174e-05 9.8174e-05 0 0 0 0 0 0 0 0 0 </BinContents> +</PhotonSigRmsRatio_8> +<PhotonBkgRmsRatio_8> + <NBinsX>30</NBinsX> + <XLow>1</XLow> + <XHigh>3</XHigh> + <BinContents>0 0.305169 0.223658 0.161034 0.122266 0.0596421 0.0407555 0.0347912 0.0198807 0.0119284 0.0119284 0.00298211 0.00298211 0.00298211 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 </BinContents> +</PhotonBkgRmsRatio_8> +<PhotonSigLongProfileStart_0> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.0520456 0.137233 0.112927 0.123682 0.130979 0.128586 0.106057 0.0800227 0.0637244 0.0455784 0.0191647 0 </BinContents> +</PhotonSigLongProfileStart_0> +<PhotonBkgLongProfileStart_0> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.0566245 0.189996 0.130592 0.128079 0.133674 0.130455 0.0956491 0.0595869 0.0397307 0.0255241 0.0100886 0 </BinContents> +</PhotonBkgLongProfileStart_0> +<PhotonSigLongProfileStart_1> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.148797 0.282878 0.168565 0.116307 0.0837729 0.0670183 0.0447675 0.0332875 0.0282789 0.0185275 0.00780107 0 </BinContents> +</PhotonSigLongProfileStart_1> +<PhotonBkgLongProfileStart_1> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.0692545 0.182469 0.141901 0.13995 0.133122 0.119222 0.0809151 0.0543351 0.0391053 0.0279102 0.0118158 0 </BinContents> +</PhotonBkgLongProfileStart_1> +<PhotonSigLongProfileStart_2> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.234425 0.335122 0.19998 0.102356 0.0572279 0.033322 0.0170269 0.00961116 0.00561058 0.00424452 0.00107333 0 </BinContents> +</PhotonSigLongProfileStart_2> +<PhotonBkgLongProfileStart_2> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.104781 0.189597 0.156571 0.140359 0.123996 0.101554 0.0642498 0.0463859 0.0344517 0.0266456 0.0114088 0 </BinContents> +</PhotonBkgLongProfileStart_2> +<PhotonSigLongProfileStart_3> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.263169 0.352345 0.190027 0.0927632 0.0499332 0.0259512 0.0123778 0.0066812 0.00344609 0.00281314 0.000492299 0 </BinContents> +</PhotonSigLongProfileStart_3> +<PhotonBkgLongProfileStart_3> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.135459 0.198514 0.179094 0.132582 0.120834 0.083673 0.0510669 0.0371614 0.0278111 0.0249341 0.00887077 0 </BinContents> +</PhotonBkgLongProfileStart_3> +<PhotonSigLongProfileStart_4> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.287348 0.345374 0.187181 0.0916173 0.045581 0.0222087 0.0110791 0.00495776 0.00278242 0.00126473 0.000607072 0 </BinContents> +</PhotonSigLongProfileStart_4> +<PhotonBkgLongProfileStart_4> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.156078 0.225413 0.178416 0.125036 0.100087 0.0699159 0.0490281 0.0365535 0.0284305 0.0211778 0.00986365 0 </BinContents> +</PhotonBkgLongProfileStart_4> +<PhotonSigLongProfileStart_5> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.30777 0.33941 0.184535 0.0889055 0.0415764 0.0208069 0.00945088 0.00410908 0.00224131 0.000933881 0.000261487 0 </BinContents> +</PhotonSigLongProfileStart_5> +<PhotonBkgLongProfileStart_5> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.204748 0.23574 0.167491 0.102539 0.0748434 0.060666 0.0487966 0.0346192 0.0323112 0.0250577 0.0131883 0 </BinContents> +</PhotonBkgLongProfileStart_5> +<PhotonSigLongProfileStart_6> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.33114 0.330641 0.172512 0.0880258 0.0432869 0.0202369 0.00825809 0.00372068 0.00131585 0.000544489 0.000317619 0 </BinContents> +</PhotonSigLongProfileStart_6> +<PhotonBkgLongProfileStart_6> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.223881 0.217622 0.142995 0.102552 0.0741454 0.0635532 0.0539239 0.0418873 0.0370727 0.0317766 0.0105922 0 </BinContents> +</PhotonBkgLongProfileStart_6> +<PhotonSigLongProfileStart_7> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.349007 0.319379 0.178525 0.0798103 0.0387021 0.0201416 0.00783667 0.00426205 0.00137485 0.00075617 0.000206228 0 </BinContents> +</PhotonSigLongProfileStart_7> +<PhotonBkgLongProfileStart_7> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.212272 0.197347 0.126866 0.10199 0.0837479 0.0671642 0.0621891 0.0497512 0.0439469 0.039801 0.0149254 0 </BinContents> +</PhotonBkgLongProfileStart_7> +<PhotonSigLongProfileStart_8> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.417632 0.298645 0.151188 0.0719615 0.0337718 0.0147261 0.00736305 0.00323974 0.000785392 0.000392696 0.000294522 0 </BinContents> +</PhotonSigLongProfileStart_8> +<PhotonBkgLongProfileStart_8> + <NBinsX>11</NBinsX> + <XLow>-0.5</XLow> + <XHigh>10.5</XHigh> + <BinContents>0 0.237575 0.150099 0.102386 0.0805169 0.0785288 0.0904573 0.0854871 0.0536779 0.05666 0.0417495 0.0228628 0 </BinContents> +</PhotonBkgLongProfileStart_8> +<PhotonSigLongProfileDiscrepancy_0> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0.000307962 0.000165826 0.00040272 0.000734371 0.0011134 0.000829128 0.00127923 0.00187146 0.00220311 0.00296117 0.00376661 0.00506953 0.00637244 0.00881245 0.011821 0.0150901 0.0186435 0.0240921 0.0282851 0.034018 0.0397271 0.0449625 0.0488475 0.0502926 0.0543671 0.0521403 0.0523772 0.0525667 0.0540592 0.0496766 0.0479473 0.0442518 0.0426646 0.0386137 0.0341127 0.0313174 0.0297302 0.0244474 0.0221496 0.0179092 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_0> +<PhotonBkgLongProfileDiscrepancy_0> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0.000201772 0.000311829 0.000348515 0.000577801 0.000715373 0.000843773 0.00117395 0.00164169 0.0017976 0.00317332 0.00505347 0.00732799 0.0106664 0.0144542 0.0197278 0.0257259 0.0322193 0.0408405 0.0457564 0.0520296 0.0555698 0.0588073 0.0567713 0.0561843 0.0550654 0.051672 0.0497643 0.0474623 0.0426197 0.0395932 0.0371719 0.032577 0.0303116 0.025891 0.0234331 0.019847 0.0166095 0.0146101 0.011877 0.00957499 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_0> +<PhotonSigLongProfileDiscrepancy_1> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 4.43243e-05 8.86486e-05 4.43243e-05 8.86486e-05 0.000132973 0.000753513 0.00168432 0.00332432 0.00558486 0.00944107 0.0123665 0.0179513 0.0242897 0.0314259 0.0373654 0.0428616 0.0487567 0.0504854 0.0563362 0.0557599 0.0583307 0.055627 0.0563362 0.056868 0.0530118 0.0464962 0.0424183 0.0396702 0.0368778 0.0342183 0.0302292 0.0274367 0.0234919 0.0214529 0.0187492 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_1> +<PhotonBkgLongProfileDiscrepancy_1> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0.000110843 0.00015518 0.000731561 0.00126361 0.00219468 0.00396817 0.00605201 0.00793633 0.0109734 0.0134341 0.0167816 0.0216365 0.0263362 0.0294841 0.0341617 0.0377973 0.0412556 0.043539 0.0481722 0.0508768 0.0528276 0.0525616 0.0509433 0.0503226 0.0500787 0.0478175 0.0451129 0.0442705 0.0415438 0.0403911 0.0369328 0.0321 0.0310138 0.027223 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_1> +<PhotonSigLongProfileDiscrepancy_2> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 4.87876e-05 4.87876e-05 4.87876e-05 0.000243938 0.000878177 0.00307362 0.00585452 0.0108796 0.0147339 0.0226375 0.0328341 0.0448846 0.0508855 0.057667 0.0659609 0.069376 0.0702542 0.0671806 0.0664487 0.0584964 0.0569839 0.0494219 0.0436649 0.0370298 0.0342977 0.0274674 0.0243938 0.0201005 0.0169781 0.0133678 0.0112212 0.00931844 0.00780602 0.005513 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_2> +<PhotonBkgLongProfileDiscrepancy_2> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 7.50582e-05 0.000225175 0.000150116 0.000600465 0.000750582 0.00187645 0.00255198 0.00442843 0.00818134 0.0106583 0.0127599 0.0179389 0.0230429 0.0286722 0.0322 0.0370037 0.0390302 0.0454102 0.0457104 0.0506643 0.0499137 0.0558433 0.0454102 0.0497636 0.0475869 0.0499137 0.0484876 0.0478121 0.042633 0.0445845 0.0433086 0.0411319 0.0390302 0.0326503 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_2> +<PhotonSigLongProfileDiscrepancy_3> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0 7.03284e-05 0.000140657 0.00126591 0.0019692 0.00555595 0.0115339 0.0174415 0.0301006 0.0395246 0.0504255 0.0635066 0.076236 0.0784162 0.0793305 0.0786975 0.0744075 0.0677263 0.0590056 0.0495815 0.0414938 0.0380477 0.0316478 0.0239117 0.0189887 0.0152613 0.0120262 0.00893171 0.00696251 0.00555595 0.00436036 0.00337576 0.00260215 0.00189887 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_3> +<PhotonBkgLongProfileDiscrepancy_3> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0 0 0.000239751 0.000479501 0.000479501 0.00263726 0.00671302 0.00743227 0.0112683 0.0210981 0.0227763 0.0345241 0.0388396 0.0445936 0.0465116 0.0484296 0.0501079 0.0515464 0.0407576 0.0469911 0.0469911 0.0383601 0.0462719 0.0409974 0.0393191 0.0383601 0.0393191 0.0326061 0.0357228 0.0414769 0.0299688 0.0309278 0.0342843 0.0299688 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_3> +<PhotonSigLongProfileDiscrepancy_4> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0 0.000303536 0.000455304 0.00212475 0.00440127 0.0130015 0.0211969 0.0348055 0.0498305 0.0649062 0.07897 0.0837254 0.0931856 0.0896949 0.0813477 0.0748217 0.0623767 0.0524612 0.0420904 0.0335913 0.0258512 0.020944 0.0177569 0.0131532 0.0102191 0.00753782 0.00581778 0.00424951 0.00288359 0.00288359 0.00217534 0.00146709 0.00131532 0.000455304 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_4> +<PhotonBkgLongProfileDiscrepancy_4> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0 0.000580215 0 0.00145054 0.00174064 0.00522193 0.0118944 0.0159559 0.0252393 0.0316217 0.0409051 0.0490281 0.0571511 0.05396 0.056861 0.0522193 0.0571511 0.0438062 0.0501886 0.0394546 0.0371337 0.0403249 0.0409051 0.0348129 0.0295909 0.0287206 0.0275602 0.0261097 0.0272701 0.0232086 0.0211778 0.0237888 0.0205976 0.024369 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_4> +<PhotonSigLongProfileDiscrepancy_5> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 3.73552e-05 0.000112066 0.00078446 0.00321255 0.00847964 0.0212178 0.0367949 0.0561076 0.0758685 0.0914456 0.101233 0.0995144 0.0946956 0.083489 0.0738513 0.0589092 0.047217 0.0342921 0.0275682 0.0207695 0.0169966 0.0108704 0.00896526 0.00691072 0.00466941 0.00381024 0.00306313 0.00194247 0.00209189 0.00138214 0.00127008 0.000971236 0.000560329 0.000560329 0.000336197 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_5> +<PhotonBkgLongProfileDiscrepancy_5> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0.000329707 0 0.00197824 0.00428619 0.0108803 0.0201121 0.0339598 0.0369271 0.0520936 0.0639631 0.0679196 0.0698978 0.0712166 0.0573689 0.060666 0.0544016 0.0418727 0.0333004 0.0356083 0.0342895 0.0257171 0.0253874 0.0194527 0.0201121 0.0201121 0.0171447 0.0220903 0.0131883 0.0141774 0.016815 0.0125288 0.0151665 0.0141774 0.0128586 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_5> +<PhotonSigLongProfileDiscrepancy_6> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0.000226871 0.00099823 0.00431054 0.0137484 0.0335315 0.0600753 0.0868007 0.10858 0.119515 0.11103 0.0988702 0.0870729 0.0678343 0.0551295 0.038931 0.0282227 0.0226417 0.0156994 0.0103907 0.00889332 0.00644312 0.00462816 0.00299469 0.00326694 0.00217796 0.00222333 0.00113435 0.00145197 0.000725986 0.000499115 0.000680612 0.000317619 0.000136122 0.000499115 0.000317619 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_6> +<PhotonBkgLongProfileDiscrepancy_6> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0 0.00144439 0.00625903 0.0125181 0.0250361 0.0544054 0.0669234 0.0813674 0.0751083 0.0702937 0.0688493 0.0548869 0.0597015 0.0481464 0.046702 0.0380356 0.033221 0.0240732 0.0327395 0.0207029 0.0221473 0.0192585 0.0221473 0.013481 0.0110737 0.0129995 0.0129995 0.0139624 0.0115551 0.0110737 0.00722195 0.00722195 0.00577756 0.00866635 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_6> +<PhotonSigLongProfileDiscrepancy_7> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0.000137485 0.000962398 0.0065993 0.0227538 0.0514883 0.0906029 0.112257 0.123874 0.11645 0.106139 0.0908778 0.0655805 0.0540318 0.0393896 0.0274971 0.0208978 0.0158108 0.0129924 0.00824912 0.00790541 0.00543067 0.00426205 0.00254348 0.00295594 0.00213102 0.00158108 0.00171857 0.00123737 0.000824912 0.00075617 0.000206228 0.000618684 0.000481199 0.000481199 0.000137485 0.000137485 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_7> +<PhotonBkgLongProfileDiscrepancy_7> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0 0.000829187 0.00414594 0.013267 0.0257048 0.0489221 0.0547264 0.0646766 0.0878939 0.0704809 0.0754561 0.0555556 0.0679934 0.0588723 0.0257048 0.0422886 0.0257048 0.0257048 0.0306799 0.0248756 0.0290216 0.0174129 0.0174129 0.0157546 0.0165837 0.0174129 0.0140962 0.00912106 0.00912106 0.00580431 0.0116086 0.0066335 0.0116086 0.00497512 0.00331675 0.0066335 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_7> +<PhotonSigLongProfileDiscrepancy_8> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0.00049087 0.00628313 0.0170823 0.0405458 0.0650893 0.0862949 0.110151 0.11074 0.114078 0.098763 0.0793246 0.0642058 0.0460436 0.0341645 0.0269978 0.0233654 0.0172786 0.0139407 0.009032 0.0075594 0.00647948 0.00402513 0.00431965 0.00196348 0.00314157 0.00117809 0.000883566 0.000883566 0.000687218 0.000785392 0.00157078 0.000589044 0.000687218 0.000294522 0.000589044 0.000294522 0.000196348 0 0 </BinContents> +</PhotonSigLongProfileDiscrepancy_8> +<PhotonBkgLongProfileDiscrepancy_8> + <NBinsX>42</NBinsX> + <XLow>-0.02</XLow> + <XHigh>0.82</XHigh> + <BinContents>0 0 0 0 0 0 0.00397614 0.0109344 0.0218688 0.0457256 0.0576541 0.0685885 0.0755467 0.083499 0.0626243 0.0616302 0.0636183 0.0487077 0.0347912 0.0397614 0.0328032 0.027833 0.0318091 0.0208747 0.0168986 0.0218688 0.0248509 0.00795229 0.0139165 0.0139165 0.0228628 0.0149105 0.0129225 0.00695825 0.00894632 0.00795229 0.00695825 0.0109344 0.00695825 0.00497018 0.00298211 0.000994036 0 0 </BinContents> +</PhotonBkgLongProfileDiscrepancy_8> +<PhotonSigPeakEnergyFraction_0> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.0229787 0.0358183 0.0323834 0.0279772 0.0252529 0.0217942 0.0179566 0.0172696 0.0155402 0.0154929 0.0144742 0.0133845 0.0126028 0.0116552 0.0111103 0.0106365 0.00992585 0.0104707 0.00978372 0.00862294 0.00916779 0.00791225 0.00805439 0.00734371 0.00698837 0.00656196 0.00686992 0.0064672 0.00594604 0.00596973 0.0053538 0.00521166 0.00497477 0.00466681 0.00360079 0.00338758 0.00281904 0.00258214 0.00189515 0.00175301 0.00151612 0.000923886 0.000686992 0.000497477 4.73788e-05 4.73788e-05 0 0 2.36894e-05 0 0.5436 0 </BinContents> +</PhotonSigPeakEnergyFraction_0> +<PhotonBkgPeakEnergyFraction_0> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.0242585 0.0401343 0.0417943 0.0420694 0.0401618 0.0384376 0.0341912 0.0318616 0.0292661 0.0270008 0.0254966 0.0240384 0.0215621 0.0192692 0.0176 0.016582 0.0153163 0.0139865 0.0134453 0.0116936 0.0110608 0.0103637 0.00957499 0.00946494 0.00808922 0.00772236 0.00766733 0.00679605 0.00663096 0.00602564 0.00559459 0.00526441 0.00477833 0.00393455 0.00391621 0.00298072 0.00294404 0.00227452 0.00212778 0.00143075 0.00118312 0.000779573 0.000394372 0.000210943 7.33716e-05 1.83429e-05 9.17145e-06 9.17145e-06 0 0 0.350515 0 </BinContents> +</PhotonBkgPeakEnergyFraction_0> +<PhotonSigPeakEnergyFraction_1> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.00726918 0.0163557 0.0164443 0.012987 0.0111697 0.0102832 0.0101946 0.00868756 0.00886486 0.00851026 0.00695891 0.00802269 0.00842161 0.00726918 0.0076681 0.00740215 0.00704756 0.00731351 0.00704756 0.00651567 0.00687026 0.00695891 0.00709188 0.00722486 0.00678161 0.00647134 0.00660432 0.00802269 0.0076681 0.00797837 0.00780107 0.0091308 0.00921945 0.00842161 0.0107265 0.0101946 0.0109924 0.0103719 0.0110367 0.0108594 0.0103719 0.00970702 0.00859891 0.00642702 0.00554053 0.00412216 0.00141838 0.00031027 0 0 0.598644 0 </BinContents> +</PhotonSigPeakEnergyFraction_1> +<PhotonBkgPeakEnergyFraction_1> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.00982066 0.0238977 0.0254938 0.0258485 0.0250504 0.0223237 0.0220133 0.0212818 0.0193975 0.0180895 0.0186881 0.0178457 0.0162052 0.0162495 0.0154958 0.0162717 0.015518 0.0145647 0.0143209 0.0140548 0.0138997 0.0146534 0.0141657 0.0127912 0.0129686 0.0123257 0.0115498 0.0117493 0.0118158 0.0106852 0.0115498 0.0112394 0.0105522 0.0100202 0.009998 0.00955463 0.00895608 0.00897825 0.00911126 0.00811368 0.00813585 0.00649538 0.0059855 0.00469973 0.00334745 0.00208384 0.00126361 0.00015518 2.21685e-05 0 0.370702 0 </BinContents> +</PhotonBkgPeakEnergyFraction_1> +<PhotonSigPeakEnergyFraction_2> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.00209787 0.00439089 0.00575694 0.00517149 0.00400059 0.00400059 0.00502513 0.00312241 0.00414695 0.0035615 0.00331756 0.00370786 0.00302483 0.00361028 0.00321998 0.0039518 0.00429331 0.00307362 0.00346392 0.00312241 0.00390301 0.00341513 0.00351271 0.00351271 0.00326877 0.00307362 0.00287847 0.00419574 0.00414695 0.00390301 0.00385422 0.00517149 0.00536664 0.00580573 0.005513 0.00653754 0.0086842 0.0103918 0.0107333 0.0144411 0.0161975 0.0178563 0.0207835 0.0223447 0.0228814 0.0206372 0.0168317 0.00878177 0.00170757 4.87876e-05 0.663561 0 </BinContents> +</PhotonSigPeakEnergyFraction_2> +<PhotonBkgPeakEnergyFraction_2> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.00382797 0.0177888 0.0214666 0.023193 0.025745 0.02582 0.0240937 0.0228177 0.0237934 0.0225175 0.0209412 0.0187645 0.0192899 0.0183893 0.0159123 0.0176387 0.015537 0.0138858 0.0152368 0.0153869 0.0146363 0.0140359 0.0123846 0.0130601 0.0124597 0.0122345 0.0132102 0.0132853 0.0108084 0.0111086 0.0114839 0.0113338 0.0111086 0.0102079 0.010283 0.010358 0.0128349 0.0105081 0.0127599 0.0114088 0.010283 0.011634 0.00998274 0.0111837 0.0121594 0.0114088 0.00840651 0.00442843 0.00120093 0 0.297756 0 </BinContents> +</PhotonBkgPeakEnergyFraction_2> +<PhotonSigPeakEnergyFraction_3> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.000773613 0.00274281 0.00478233 0.00386806 0.00323511 0.00295379 0.00414938 0.00253182 0.00210985 0.00302412 0.00302412 0.00358675 0.00309445 0.00218018 0.00302412 0.00260215 0.00288347 0.00281314 0.00239117 0.00218018 0.00330544 0.00267248 0.00168788 0.00218018 0.00239117 0.00302412 0.00232084 0.00365708 0.00281314 0.00330544 0.00337576 0.00337576 0.00379774 0.00492299 0.00414938 0.004712 0.00696251 0.00647022 0.00696251 0.00858007 0.0128701 0.0160349 0.0175118 0.024193 0.0317884 0.0353049 0.0373444 0.0300302 0.0111822 0.000281314 0.644841 0 </BinContents> +</PhotonSigPeakEnergyFraction_3> +<PhotonBkgPeakEnergyFraction_3> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.0014385 0.00982978 0.0191801 0.0196596 0.0246943 0.0230161 0.0184608 0.0246943 0.0227763 0.0249341 0.0220571 0.0179813 0.0206186 0.0191801 0.0206186 0.0198993 0.0201391 0.0179813 0.0165428 0.0170223 0.0177415 0.016303 0.0141453 0.0139055 0.0122273 0.0131863 0.0129465 0.013426 0.0103093 0.013426 0.0127068 0.0107888 0.0112683 0.0103093 0.010549 0.0131863 0.0112683 0.0119875 0.0122273 0.0119875 0.010549 0.00982978 0.0136658 0.0141453 0.014385 0.0177415 0.015344 0.013426 0.00575402 0 0.250539 0 </BinContents> +</PhotonBkgPeakEnergyFraction_3> +<PhotonSigPeakEnergyFraction_4> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.000354126 0.00278242 0.00389538 0.00374361 0.00344008 0.00308595 0.00323772 0.00298477 0.00258006 0.00278242 0.00308595 0.002833 0.00217534 0.0023777 0.0019224 0.00298477 0.00258006 0.00182122 0.00263065 0.00303536 0.00268124 0.00217534 0.00197299 0.00273183 0.00202357 0.00263065 0.00298477 0.00273183 0.00268124 0.00288359 0.0033389 0.00349067 0.00344008 0.00323772 0.00409774 0.00394597 0.00627308 0.00490717 0.0061719 0.00748723 0.00860019 0.0109779 0.0134568 0.0185663 0.0262559 0.0361208 0.0457328 0.0574695 0.0339961 0.0023777 0.620226 0 </BinContents> +</PhotonSigPeakEnergyFraction_4> +<PhotonBkgPeakEnergyFraction_4> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.000580215 0.00899333 0.0130548 0.0171163 0.0159559 0.0174064 0.0205976 0.0185669 0.0205976 0.0174064 0.0165361 0.0191471 0.0145054 0.0153757 0.0232086 0.0174064 0.0165361 0.0159559 0.0191471 0.0147955 0.0145054 0.0176965 0.0150856 0.0168262 0.0150856 0.013635 0.0156658 0.0156658 0.0145054 0.0124746 0.0150856 0.0116043 0.0118944 0.0191471 0.0165361 0.0124746 0.0133449 0.0121845 0.0127647 0.010734 0.0142153 0.0139252 0.0116043 0.0153757 0.0153757 0.0185669 0.0226284 0.0252393 0.0127647 0.00116043 0.239339 0 </BinContents> +</PhotonBkgPeakEnergyFraction_4> +<PhotonSigPeakEnergyFraction_5> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0 0.00115801 0.00310049 0.00414643 0.00354875 0.00328726 0.0031752 0.0031752 0.00381024 0.00351139 0.0031752 0.00310049 0.00246545 0.00306313 0.00339933 0.00354875 0.00324991 0.00313784 0.00313784 0.00261487 0.00276429 0.00272693 0.00257751 0.00242809 0.00246545 0.0025028 0.00246545 0.00302578 0.002839 0.00324991 0.00332462 0.00298842 0.00339933 0.00354875 0.00351139 0.00377288 0.00496825 0.00429585 0.00545387 0.00612626 0.00668659 0.00870377 0.0112066 0.0146433 0.0182294 0.0270078 0.0429959 0.0679492 0.0812103 0.0189391 0.574188 0 </BinContents> +</PhotonSigPeakEnergyFraction_5> +<PhotonBkgPeakEnergyFraction_5> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0.000329707 0.00230795 0.00725354 0.0125288 0.0108803 0.0118694 0.013518 0.0161556 0.01121 0.0161556 0.0131883 0.0105506 0.0151665 0.0121991 0.0194527 0.0194527 0.016815 0.0118694 0.0161556 0.0145071 0.0181339 0.0158259 0.0187933 0.0187933 0.013518 0.0171447 0.0154962 0.0171447 0.0184636 0.019123 0.0184636 0.016815 0.0145071 0.0164853 0.0148368 0.0174744 0.0154962 0.0201121 0.0148368 0.0164853 0.0148368 0.0197824 0.0161556 0.0141774 0.0178042 0.0184636 0.0243983 0.0339598 0.0425321 0.00791296 0.200462 0 </BinContents> +</PhotonBkgPeakEnergyFraction_5> +<PhotonSigPeakEnergyFraction_6> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0 9.07482e-05 0.00099823 0.00249558 0.00390217 0.00417442 0.00521802 0.00571714 0.00598938 0.00626163 0.00671537 0.00417442 0.00562639 0.00535414 0.00521802 0.00544489 0.00517265 0.00435591 0.00499115 0.00412904 0.00362993 0.00462816 0.00435591 0.00376605 0.00362993 0.00362993 0.00317619 0.00390217 0.00390217 0.0040383 0.00362993 0.00358455 0.0036753 0.00372068 0.00390217 0.00444666 0.0040383 0.00539952 0.00467353 0.0052634 0.00780435 0.00644312 0.00907482 0.0111167 0.0127501 0.0186034 0.0286764 0.0484595 0.0945143 0.0617088 0.533826 0 </BinContents> +</PhotonSigPeakEnergyFraction_6> +<PhotonBkgPeakEnergyFraction_6> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0 0 0.00337025 0.0052961 0.00481464 0.00962927 0.0101107 0.0120366 0.0125181 0.0144439 0.0154068 0.0139624 0.0125181 0.00914781 0.00770342 0.0120366 0.0221473 0.0178142 0.0173327 0.0163698 0.0173327 0.013481 0.0178142 0.0125181 0.0115551 0.0144439 0.0154068 0.0129995 0.0129995 0.0168512 0.0202215 0.0173327 0.0202215 0.0192585 0.0173327 0.01974 0.0226288 0.0178142 0.0216659 0.0211844 0.0149254 0.0250361 0.0250361 0.0235917 0.0178142 0.0235917 0.0240732 0.0385171 0.052961 0.0284064 0.166586 0 </BinContents> +</PhotonBkgPeakEnergyFraction_6> +<PhotonSigPeakEnergyFraction_7> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0 0 0 0.000343713 0.00123737 0.00206228 0.00364336 0.00371211 0.00543067 0.00577439 0.00653056 0.00714924 0.00845535 0.00941775 0.0080429 0.00983021 0.00879906 0.00996769 0.00859284 0.00866158 0.00824912 0.00859284 0.00955523 0.00818038 0.0075617 0.00687427 0.00818038 0.00687427 0.0075617 0.00632433 0.00666804 0.0061181 0.0061181 0.00701175 0.00577439 0.0061181 0.0056369 0.0056369 0.00584313 0.00604936 0.00673678 0.00879906 0.0061181 0.00859284 0.0120987 0.0136111 0.0181481 0.0317591 0.0653743 0.0873032 0.494879 0 </BinContents> +</PhotonSigPeakEnergyFraction_7> +<PhotonBkgPeakEnergyFraction_7> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0 0 0 0 0.000829187 0.00414594 0.00414594 0.00829187 0.0116086 0.00746269 0.00912106 0.0116086 0.00995025 0.013267 0.00497512 0.00912106 0.0107794 0.00995025 0.013267 0.00912106 0.0157546 0.013267 0.0165837 0.0140962 0.0165837 0.0174129 0.0165837 0.00912106 0.013267 0.0165837 0.0190713 0.00995025 0.0149254 0.0257048 0.0174129 0.0165837 0.0232172 0.0232172 0.0257048 0.0257048 0.0331675 0.0298507 0.0323383 0.0323383 0.0447761 0.0339967 0.0339967 0.0348259 0.053068 0.0373134 0.145937 0 </BinContents> +</PhotonBkgPeakEnergyFraction_7> +<PhotonSigPeakEnergyFraction_8> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0.000294522 0.00049087 0.000589044 0.000883566 0.00166896 0.00245435 0.00284704 0.00598861 0.00314157 0.00657766 0.00657766 0.006774 0.00775574 0.00795209 0.00824661 0.00834479 0.0105046 0.0120754 0.00932653 0.0100137 0.013548 0.0144316 0.0123699 0.0127626 0.0145297 0.0126644 0.015806 0.0151188 0.0138425 0.0130571 0.0149224 0.0138425 0.0154133 0.0170823 0.0141371 0.0122717 0.0139407 0.0139407 0.0151188 0.017475 0.0177695 0.0226782 0.0336737 0.0731396 0.453956 0 </BinContents> +</PhotonSigPeakEnergyFraction_8> +<PhotonBkgPeakEnergyFraction_8> + <NBinsX>52</NBinsX> + <XLow>-0.02</XLow> + <XHigh>1.02</XHigh> + <BinContents>0 0 0 0 0 0 0 0 0.000994036 0 0 0.000994036 0.00497018 0.00198807 0 0.00298211 0.00695825 0.00298211 0.00894632 0.00994036 0.0109344 0.00497018 0.0109344 0.00695825 0.00596421 0.0149105 0.0119284 0.0109344 0.0119284 0.0139165 0.00994036 0.0109344 0.0139165 0.0168986 0.0109344 0.0198807 0.0168986 0.0119284 0.026839 0.0208747 0.026839 0.0328032 0.0208747 0.027833 0.027833 0.028827 0.0337972 0.0497018 0.05666 0.0666004 0.0785288 0.0934394 0.153082 0 </BinContents> +</PhotonBkgPeakEnergyFraction_8> +<PhotonSigMinDistanceToTrack_0> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.00172933 0.00206098 0.00151612 0.00187146 0.00139767 0.00130292 0.00132661 0.00106602 0.00130292 0.00127923 0.00130292 0.00113709 0.00132661 0.00137398 0.00132661 0.00149243 0.0015635 0.00144505 0.00158719 0.00168195 0.00144505 0.00198991 0.00168195 0.00258214 0.00198991 0.00232156 0.00227418 0.00260583 0.00244001 0.00213204 0.0028901 0.00281904 0.00307962 0.00305593 0.00329282 0.00293748 0.931372 </BinContents> +</PhotonSigMinDistanceToTrack_0> +<PhotonBkgMinDistanceToTrack_0> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.0364749 0.0241851 0.020434 0.017279 0.0143625 0.0117578 0.011125 0.00979511 0.00868536 0.00851111 0.00896968 0.00886879 0.00886879 0.00831851 0.00877708 0.00855696 0.00810756 0.00819928 0.00778656 0.00751142 0.00762148 0.00792413 0.00791496 0.0076123 0.00805253 0.00827265 0.00818093 0.00819011 0.00824513 0.00779573 0.00800668 0.0076765 0.00853862 0.00822679 0.00796999 0.00804336 0.625154 </BinContents> +</PhotonBkgMinDistanceToTrack_0> +<PhotonSigMinDistanceToTrack_1> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.00616107 0.00390054 0.00328 0.0027481 0.00217189 0.00252648 0.00190594 0.00203892 0.00141838 0.00150703 0.00141838 0.00132973 0.00164 0.00177297 0.00212757 0.00199459 0.00234919 0.0018173 0.00177297 0.00164 0.00208324 0.00212757 0.00208324 0.00190594 0.00203892 0.00168432 0.00212757 0.00243784 0.00208324 0.00203892 0.00226054 0.00265946 0.00270378 0.0027481 0.00248216 0.00239351 0.918621 </BinContents> +</PhotonSigMinDistanceToTrack_1> +<PhotonBkgMinDistanceToTrack_1> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.031834 0.0211044 0.0196413 0.0171806 0.0157618 0.0126361 0.011439 0.0111729 0.0107739 0.0106409 0.0105522 0.0115941 0.0124144 0.0119488 0.0123257 0.0106852 0.011439 0.0117272 0.0105522 0.0106631 0.0108182 0.0111508 0.0112173 0.0114168 0.0106631 0.0111286 0.0107074 0.0100645 0.0102197 0.00968765 0.0090891 0.00886741 0.00964331 0.00937729 0.0103305 0.00982066 0.559711 </BinContents> +</PhotonBkgMinDistanceToTrack_1> +<PhotonSigMinDistanceToTrack_2> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.0035615 0.00175635 0.00136605 0.00131727 0.00107333 0.00082939 0.000780602 0.000487876 0.000634239 0.000683027 0.00126848 0.000780602 0.000731814 0.00151242 0.00131727 0.000731814 0.00180514 0.00151242 0.0011709 0.00131727 0.000878177 0.00190272 0.00141484 0.00234181 0.00239059 0.00219544 0.00209787 0.00190272 0.00180514 0.00165878 0.00146363 0.00180514 0.00185393 0.00224423 0.00200029 0.00239059 0.945016 </BinContents> +</PhotonSigMinDistanceToTrack_2> +<PhotonBkgMinDistanceToTrack_2> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.0443594 0.0261202 0.0249944 0.0190648 0.0199655 0.014186 0.014186 0.0146363 0.0114088 0.0124597 0.0117091 0.0150867 0.0177888 0.0149366 0.0159874 0.0135105 0.0146363 0.0144862 0.0126098 0.0117841 0.0114839 0.0118592 0.0107333 0.011559 0.0127599 0.0119342 0.00960745 0.0096825 0.00990768 0.00893192 0.00983262 0.00863169 0.00960745 0.0091571 0.00840651 0.0073557 0.494633 </BinContents> +</PhotonBkgMinDistanceToTrack_2> +<PhotonSigMinDistanceToTrack_3> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.000984598 0.00119558 0.000773613 0.000351642 0.000703284 0.000562627 0.000632956 0.000773613 0.000843941 0.000492299 0.00091427 0.00091427 0.00133624 0.00140657 0.00119558 0.00105493 0.00189887 0.00182854 0.00161755 0.00232084 0.00182854 0.00182854 0.00253182 0.00203952 0.0019692 0.00218018 0.00239117 0.00274281 0.00302412 0.00253182 0.00154723 0.00182854 0.00302412 0.00210985 0.00232084 0.00302412 0.941276 </BinContents> +</PhotonSigMinDistanceToTrack_3> +<PhotonBkgMinDistanceToTrack_3> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.0366818 0.0249341 0.0246943 0.0230161 0.017262 0.0165428 0.0151043 0.014385 0.012467 0.013426 0.0122273 0.0182211 0.0167825 0.0129465 0.0189403 0.0155838 0.011508 0.0139055 0.0103093 0.0129465 0.00935028 0.0141453 0.00911053 0.0107888 0.010549 0.0100695 0.0100695 0.0112683 0.00647327 0.00671302 0.00791177 0.00599377 0.00791177 0.00863102 0.00959003 0.00647327 0.513066 </BinContents> +</PhotonBkgMinDistanceToTrack_3> +<PhotonSigMinDistanceToTrack_4> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.00131532 0.000556483 0.000657662 0.000505894 0.000708251 0.000303536 0.000708251 0.000657662 0.000505894 0.00106238 0.00121414 0.00116356 0.00131532 0.00151768 0.00172004 0.00182122 0.00242829 0.00166945 0.00202357 0.00202357 0.00207416 0.002833 0.002833 0.00227652 0.00242829 0.00273183 0.00222593 0.00217534 0.00278242 0.00242829 0.00217534 0.002833 0.00258006 0.00318713 0.00303536 0.0033389 0.934183 </BinContents> +</PhotonSigMinDistanceToTrack_4> +<PhotonBkgMinDistanceToTrack_4> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.0272701 0.0194372 0.0203075 0.016246 0.0156658 0.0139252 0.0171163 0.0124746 0.010734 0.0121845 0.0110241 0.0121845 0.016246 0.016246 0.0168262 0.0113142 0.0124746 0.00957354 0.00841311 0.00899333 0.00841311 0.0101538 0.0104439 0.008123 0.00928343 0.00696258 0.00841311 0.00725268 0.00841311 0.00638236 0.0078329 0.00754279 0.00522193 0.00551204 0.00667247 0.00435161 0.590368 </BinContents> +</PhotonBkgMinDistanceToTrack_4> +<PhotonSigMinDistanceToTrack_5> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.000933881 0.000522973 0.000747105 0.000672394 0.000410908 0.0010833 0.000859171 0.0010833 0.000933881 0.00134479 0.0014195 0.00235338 0.00201718 0.002839 0.00268958 0.00246545 0.00287635 0.00276429 0.00257751 0.0031752 0.00276429 0.00332462 0.00351139 0.00339933 0.00306313 0.00384759 0.0042585 0.00369817 0.0039223 0.00324991 0.00373552 0.00474412 0.0035861 0.00362346 0.0039223 0.00366081 0.907919 </BinContents> +</PhotonSigMinDistanceToTrack_5> +<PhotonBkgMinDistanceToTrack_5> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.0171447 0.0145071 0.0201121 0.0148368 0.0148368 0.0171447 0.0121991 0.0121991 0.0105506 0.00923178 0.0118694 0.0158259 0.0121991 0.0131883 0.0105506 0.00857237 0.00824266 0.0108803 0.0115397 0.00560501 0.013518 0.0115397 0.0102209 0.00956149 0.00791296 0.00725354 0.0049456 0.00461589 0.00461589 0.00791296 0.00824266 0.0049456 0.00824266 0.00758325 0.00461589 0.00461589 0.628421 </BinContents> +</PhotonBkgMinDistanceToTrack_5> +<PhotonSigMinDistanceToTrack_6> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.000907482 0.00077136 0.000725986 0.0010436 0.0010436 0.00167884 0.00231408 0.00204183 0.00208721 0.00267707 0.00322156 0.00412904 0.00458279 0.00394755 0.00394755 0.00458279 0.00444666 0.00512727 0.00449204 0.00539952 0.00503653 0.00585326 0.00499115 0.00535414 0.00589863 0.00512727 0.00503653 0.00512727 0.00576251 0.00626163 0.00571714 0.006307 0.00558102 0.00567176 0.00589863 0.00512727 0.85208 </BinContents> +</PhotonSigMinDistanceToTrack_6> +<PhotonBkgMinDistanceToTrack_6> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.0163698 0.0105922 0.0163698 0.0144439 0.0101107 0.0105922 0.0115551 0.0125181 0.00866635 0.00914781 0.0110737 0.0120366 0.0105922 0.0105922 0.0110737 0.00866635 0.00818488 0.00818488 0.00722195 0.0052961 0.00674049 0.0052961 0.00577756 0.0120366 0.00866635 0.00674049 0.00577756 0.00625903 0.0052961 0.00722195 0.00722195 0.00385171 0.00577756 0.00866635 0.00625903 0.00770342 0.677419 </BinContents> +</PhotonBkgMinDistanceToTrack_6> +<PhotonSigMinDistanceToTrack_7> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.00123737 0.00116863 0.00171857 0.00233725 0.00171857 0.00247474 0.00295594 0.00398708 0.00350588 0.00426205 0.00460576 0.00488073 0.00783667 0.00701175 0.00666804 0.00701175 0.00721798 0.00584313 0.00776792 0.00653056 0.00714924 0.00907404 0.00831787 0.00797415 0.00701175 0.00790541 0.0075617 0.00763044 0.00653056 0.0065993 0.0065993 0.00769918 0.00584313 0.00749295 0.00776792 0.0061181 0.791985 </BinContents> +</PhotonSigMinDistanceToTrack_7> +<PhotonBkgMinDistanceToTrack_7> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.0124378 0.00829187 0.0165837 0.00912106 0.00829187 0.00829187 0.00497512 0.00746269 0.00829187 0.00829187 0.00746269 0.00746269 0.00912106 0.0149254 0.00580431 0.00746269 0.0066335 0.0066335 0.00746269 0.0107794 0.00414594 0.00580431 0.00995025 0.00746269 0.00497512 0.00912106 0.000829187 0.00912106 0.013267 0.00497512 0.0066335 0.00746269 0.00331675 0.00331675 0.00829187 0.00331675 0.722222 </BinContents> +</PhotonBkgMinDistanceToTrack_7> +<PhotonSigMinDistanceToTrack_8> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.00196348 0.00284704 0.002258 0.00363244 0.00481052 0.00441783 0.00530139 0.00539957 0.00500687 0.00638131 0.00549774 0.00854114 0.01129 0.00962105 0.0113882 0.0106028 0.00952287 0.0107991 0.0103083 0.0094247 0.0113882 0.009032 0.0109955 0.0106028 0.0104064 0.0113882 0.0102101 0.00932653 0.0106028 0.009032 0.0108973 0.00873748 0.0098174 0.00863931 0.00795209 0.00932653 0.702631 </BinContents> +</PhotonSigMinDistanceToTrack_8> +<PhotonBkgMinDistanceToTrack_8> + <NBinsX>40</NBinsX> + <XLow>0</XLow> + <XHigh>20</XHigh> + <BinContents>0 0 0 0 0 0.00397614 0.00298211 0.0129225 0.00497018 0.00497018 0.00894632 0.00894632 0.00298211 0.00994036 0.00795229 0.00596421 0.00894632 0.0149105 0.00298211 0.0119284 0.00695825 0.0149105 0.00695825 0.00695825 0.00596421 0.0129225 0.00994036 0.0149105 0.0119284 0.00397614 0.00596421 0.0149105 0.00695825 0.00994036 0.00795229 0.00994036 0.00596421 0.0109344 0.00497018 0.00497018 0.00894632 0.699801 </BinContents> +</PhotonBkgMinDistanceToTrack_8> diff --git a/Reconstruction/PFA/Pandora/PandoraSettingsDefault.xml b/Reconstruction/PFA/Pandora/PandoraSettingsDefault.xml index 8bf0e965b235f94400e68629ec565f5f9fd14db2..be26c3643a7017db4f3158da6d6ef338fd30a92e 100644 --- a/Reconstruction/PFA/Pandora/PandoraSettingsDefault.xml +++ b/Reconstruction/PFA/Pandora/PandoraSettingsDefault.xml @@ -69,7 +69,7 @@ <ClusterListName>PhotonClusters</ClusterListName> <ReplaceCurrentClusterList>false</ReplaceCurrentClusterList> <ShouldMakePdfHistograms>false</ShouldMakePdfHistograms> - <HistogramFile>/junofs/users/wxfang/MyGit/MarlinPandora/scripts/PandoraLikelihoodData9EBin.xml</HistogramFile> + <HistogramFile>/junofs/users/wxfang/MyGit/tmp/fork_update_pandora/CEPCSW/Reconstruction/PFA/Pandora/PandoraLikelihoodData9EBin.xml</HistogramFile> </algorithm> <!-- Clustering parent algorithm runs a daughter clustering algorithm -->