48 produces<std::vector<sim::MCHitCollection>>();
49 if (fMakeMCWire) produces<std::vector<sim::MCWireCollection>>();
56 std::unique_ptr<std::vector<sim::MCHitCollection>> hits_v(
57 new std::vector<sim::MCHitCollection>());
58 std::unique_ptr<std::vector<sim::MCWireCollection>> wires_v(
59 new std::vector<sim::MCWireCollection>());
62 wires_v->reserve(nch);
63 for (
size_t ch = 0; ch < nch; ++ch) {
73 <<
"Did not find sim::SimChannel with a label: " <<
fLArG4ModuleName.c_str() << std::endl;
76 for (
size_t simch_index = 0; simch_index < simchArray->size(); ++simch_index) {
80 size_t ch = simch_ptr->
Channel();
82 if (ch >= hits_v->size())
85 <<
"Channel number " << ch <<
" exceeds total # of channels: " << nch << std::endl;
87 auto& mchits = hits_v->at(ch);
88 auto& mcwires = wires_v->at(ch);
90 std::map<sim::MCEnDep, sim::MCWire> edep_wire_map;
91 auto tdc_ide_map = simch_ptr->
TDCIDEMap();
100 if (
fVerbose) std::cout << std::endl <<
"Processing Ch: " << ch << std::endl;
102 for (
auto const& tdc_ide_pair : tdc_ide_map) {
104 auto const& tdc = tdc_ide_pair.first;
105 auto const& ide_v = tdc_ide_pair.second;
107 for (
auto const& ide : ide_v) {
116 auto edep_iter = edep_wire_map.insert(std::make_pair(edep, wire));
119 (edep_iter).first->second.StartTDC() + (edep_iter).first->second.size() - 1;
123 if (edep_iter.second) std::cout << std::endl;
125 std::cout <<
" Track: " << ide.trackID <<
" Vtx: " << ide.x <<
" " << ide.y <<
" " 126 << ide.z <<
" " << ide.energy <<
" ... @ TDC = " << tdc <<
" ... " 127 << ide.numElectrons << std::endl;
130 if (!(edep_iter.second)) {
131 if (last_tdc + 1 != tdc) {
continue; }
133 (edep_iter).first->second.push_back(ide.numElectrons);
140 for (
auto const& edep_wire_pair : edep_wire_map) {
142 auto const&
edep = edep_wire_pair.first;
143 auto const& wire = edep_wire_pair.second;
147 float vtx[3] = {float(
edep.Vertex()[0]),
float(
edep.Vertex()[1]),
float(
edep.Vertex()[2])};
153 for (
size_t wire_index = 0; wire_index < wire.size(); ++wire_index) {
155 auto q = wire.at(wire_index);
161 max_time = wire.StartTDC() + wire_index;
167 mchits.push_back(hit);
171 mcwires.push_back(wire);
177 std::sort((*hits_v).begin(), (*hits_v).end());
178 e.
put(std::move(hits_v));
181 std::sort((*wires_v).begin(), (*wires_v).end());
182 e.
put(std::move(wires_v));
void SetCharge(float qsum, float amp)
Setter function for charge/amplitude.
void SetParticleInfo(const float vtx[], const float energy, const int trackId)
Setter function for partile info.
void SetTime(const float peak, const float width)
Setter function for time.
EDProducer(fhicl::ParameterSet const &pset)
void SetStartTDC(const unsigned int start)
Setter function for time.
void produce(art::Event &e) override
cout<< "Opened file "<< fin<< " ixs= "<< ixs<< endl;if(ixs==0) hhh=(TH1F *) fff-> Get("h1")
void SetTrackId(unsigned int id)
bool isValid() const noexcept
PutHandle< PROD > put(std::unique_ptr< PROD > &&edp, std::string const &instance={})
#define DEFINE_ART_MODULE(klass)
void SetVertex(float x, float y, float z)
std::string fLArG4ModuleName
Detector simulation of raw signals on wires.
raw::ChannelID_t Channel() const
Returns the readout channel this object describes.
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
object containing MC truth information necessary for making RawDigits and doing back tracking ...
TDCIDEs_t const & TDCIDEMap() const
Returns all the deposited energy information as stored.
MCHitFinder(fhicl::ParameterSet const &p)
cet::coded_exception< error, detail::translate > exception