Skip to content

RICH hit associations are not valid #2418

@veprbl

Description

@veprbl

Maybe I missed originally, @DraTeots that you said, this is still not resolved for DRICH. Shouldn't have closed, and nobody pointed out.

I'm able to reproduce with

#include <podio/Frame.h>
#include <podio/ROOTReader.h>
#include <edm4eic/MCRecoTrackerHitAssociationCollection.h>

void test(std::string input_file="root://dtn-eic.jlab.org//volatile/eic/EPIC//RECO/25.10.4/epic_craterlake/DIS/BeAGLE1.03.02-1.0/eH2/10x130/q2_1to1000/BeAGLE1.03.02-1.0_DIS_eH2_10x130_q2_1to1000_ab.0000.eicrecon.edm4eic.root") {
    auto reader = podio::ROOTReader();
    reader.openFile(input_file);
    for (size_t i = 0; i < reader.getEntries(podio::Category::Event); i++) {
        auto frame = podio::Frame(reader.readNextEntry(podio::Category::Event));
        std::cerr << "event " << i << std::endl;
	auto& rawAssoc = frame.get<edm4eic::MCRecoTrackerHitAssociationCollection>("DRICHRawHitsAssociations");
	for (const auto& assoc : rawAssoc) 
	{
		std::cout << "rec: " << assoc.getRawHit() << std::endl;
	      std::cout << "sim: " << assoc.getSimHit() << std::endl;
	      std::cout << "sim: " << assoc.getSimHit().getObjectID() << std::endl;
	}
    }
}
sim: [not available]
sim: ffffffff|-1
rec:  id: b5b1501c|165
 cellID : 1970354927316344
 charge : 299
 timeStamp : 218

The issue must be that we don't store DRICH sim hits:

"DRICHRawHits",
"DRICHRawHitsAssociations",

same for pfRICH
https://github.com/eic/EICrecon/blob/ee80688cda5d0e4260adf5d9ec6c75718580e333/src/services/io/podio/JEventProcessorPODIO.cc#L114C8-L115

we could add those if it doesn't take too much space.

Originally posted by @veprbl in #2324

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions