SimTestPulse Class Reference

Inheritance diagram for SimTestPulse:

Public Member Functions
	SimTestPulse (fhicl::ParameterSet const &p)
	SimTestPulse (SimTestPulse const &)=delete
	SimTestPulse (SimTestPulse &&)=delete
SimTestPulse &	operator= (SimTestPulse const &)=delete
SimTestPulse &	operator= (SimTestPulse &&)=delete
void	produce (art::Event &e) override
void	beginRun (art::Run &run) override
void	beginJob () override
void	endJob () override

Private Attributes
bool	fVerbose
	verbosity More...
double	fTriggerTime
	Trigger timing in electronics clock [us]. More...
std::vector< double >	fSimTime_v
	Charge timing in electronics clock [us]. More...
std::vector< double >	fY_v
	Charge Y position in detector coordinate [cm]. More...
std::vector< double >	fZ_v
	Charge Z position in detector coordinate [cm]. More...
std::vector< double >	fNumElectrons_v
	Charge amount in electron count. More...
std::vector< int >	fTick_v
	Corresponding tick. More...
std::vector< int >	fPlane0Channel_v
	Resulting plane 0 channel number where charge is deposited. More...
std::vector< int >	fPlane1Channel_v
	Resulting plane 1 channel number where charge is deposited. More...
std::vector< int >	fPlane2Channel_v
	Resulting plane 2 channel number where charge is deposited. More...
std::vector< int >	fPlane0Wire_v
	Resulting plane 0 wire number where charge is deposited. More...
std::vector< int >	fPlane1Wire_v
	Resulting plane 1 wire number where charge is deposited. More...
std::vector< int >	fPlane2Wire_v
	Resulting plane 2 wire number where charge is deposited. More...
int	_run
int	_subrun
int	_event
TTree *	fTupleTree
	output analysis tree More...

Detailed Description

Definition at line 39 of file SimTestPulse_module.cc.

Constructor & Destructor Documentation

SimTestPulse::SimTestPulse ( fhicl::ParameterSet const &  p )

explicit

Definition at line 90 of file SimTestPulse_module.cc.

  : EDProducer{p}
// Initialize member data here.
{
    produces< std::vector<sim::SimChannel> >();
    produces< std::vector<sim::SimEnergyDeposit> >();
    produces< std::vector<raw::Trigger> >();
    produces< sumdata::RunData, art::InRun >();
    
    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
    fTriggerTime = clockData.TriggerTime();

//    fTriggerTime    = p.get< double>              ("TriggerTime_us");
    fSimTime_v      = p.get< std::vector<double> >("SimTimeArray_us");
    fY_v            = p.get< std::vector<double> >("YArray_cm");
    fZ_v            = p.get< std::vector<double> >("ZArray_cm");
    fNumElectrons_v = p.get< std::vector<double> >("NumElectronsArray");
    fVerbose        = p.get< bool>                ("Verbose",false);
    
    assert( fSimTime_v.size() == fY_v.size() &&
            fSimTime_v.size() == fZ_v.size() &&
            fSimTime_v.size() == fNumElectrons_v.size() );
}

SimTestPulse::SimTestPulse ( SimTestPulse const &  )

delete

SimTestPulse::SimTestPulse ( SimTestPulse &&  )

delete

Member Function Documentation

void SimTestPulse::beginJob ( )

override

Definition at line 114 of file SimTestPulse_module.cc.

{
    art::ServiceHandle<art::TFileService> tfs;

    fTupleTree = tfs->make<TTree>("simTestPulse", "Tree by SimTestPulse_module");
    fTupleTree->Branch("run",&_run,"run/I");
    fTupleTree->Branch("subrun",&_subrun,"subrun/I");
    fTupleTree->Branch("event",&_event,"event/I");
    fTupleTree->Branch("trigger_time",&fTriggerTime,"trigger_time/D");
    fTupleTree->Branch("charge_time_v","std::vector<double>",&fSimTime_v);
    fTupleTree->Branch("tick_v","std::vector<int>",&fTick_v);
    fTupleTree->Branch("y_v","std::vector<double>",&fY_v);
    fTupleTree->Branch("z_v","std::vector<double>",&fZ_v);
    fTupleTree->Branch("e_v","std::vector<double>",&fNumElectrons_v);
    
    fTupleTree->Branch("ch_plane0","std::vector<int>",&fPlane0Channel_v);
    fTupleTree->Branch("ch_plane1","std::vector<int>",&fPlane1Channel_v);
    fTupleTree->Branch("ch_plane2","std::vector<int>",&fPlane2Channel_v);
    
    fTupleTree->Branch("wire_plane0","std::vector<int>",&fPlane0Wire_v);
    fTupleTree->Branch("wire_plane1","std::vector<int>",&fPlane1Wire_v);
    fTupleTree->Branch("wire_plane2","std::vector<int>",&fPlane2Wire_v);
}

void SimTestPulse::beginRun ( art::Run &  run )

override

Definition at line 143 of file SimTestPulse_module.cc.

{
  // grab the geometry object to see what geometry we are using                                                                   
    art::ServiceHandle<geo::Geometry> geo;
    
    std::unique_ptr<sumdata::RunData> runData(new sumdata::RunData(geo->DetectorName()));
    
    run.put(std::move(runData));
    
    return;
}

void SimTestPulse::endJob ( )

override

Definition at line 138 of file SimTestPulse_module.cc.

{
    alternative::ParamHolder::destroy();
}

SimTestPulse& SimTestPulse::operator= ( SimTestPulse const &  )

delete

SimTestPulse& SimTestPulse::operator= ( SimTestPulse &&  )

delete

void SimTestPulse::produce ( art::Event &  e )

override

track id, keep 0 = invalid

energy, keep -1 = invalid

Definition at line 155 of file SimTestPulse_module.cc.

{
    std::unique_ptr<std::vector<raw::Trigger> > trigger_v(new std::vector<raw::Trigger> );
    trigger_v->push_back(raw::Trigger(0,fTriggerTime,fTriggerTime,1));
    
    std::unique_ptr<std::vector<sim::SimChannel> > simch_v(new std::vector<sim::SimChannel> );

    std::unique_ptr<std::vector<sim::SimEnergyDeposit>> simDep_v(new std::vector<sim::SimEnergyDeposit>);
    
    geo::Point_t chargeDepCoords = {0., 0., 0.};

    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
    art::ServiceHandle<geo::Geometry> geo;

    auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataForJob(clockData);
    
    fPlane0Channel_v.clear();
    fPlane1Channel_v.clear();
    fPlane2Channel_v.clear();
    fPlane0Wire_v.clear();
    fPlane1Wire_v.clear();
    fPlane2Wire_v.clear();
    fTick_v.clear();
    
    auto& pholder = alternative::ParamHolder::get();
    pholder.Clear();

    for(size_t index=0; index < fSimTime_v.size(); ++index) {
        
        int tdc = fSimTime_v[index] / clockData.TPCClock().TickPeriod();
        fTick_v.push_back(clockData.TPCTDC2Tick(tdc));
        
        if(fVerbose)
            std::cout << "[BUFFOON!] Charge injection id " << index << " @ TDC=" << tdc
            << " @ Y=" << fY_v[index] << " @ Z=" << fZ_v[index]
            << " with " << fNumElectrons_v[index] << " electrons " << std::endl;
        
        if(tdc<0) {
            std::cerr << "\033[93m[WARNING BUFFOON!]\033[00m ignoring fSimTime " << fSimTime_v[index]
            << " as it results in negative TDC (invalid)" << std::endl;
            continue;
        }

        const geo::PlaneGeo& planeGeo = geo->Plane(geo::PlaneID(0,0,0)); // Get the coordinates of the first wire plane

        TVector3 planeCoords = planeGeo.GetCenter();
        TVector3 planeNormal = planeGeo.GetNormalDirection();

        // Assume C=0, T=1
        chargeDepCoords = geo::Point_t(planeCoords[0] + 1. * planeNormal[0],fY_v[index],fZ_v[index]);
        
        alternative::TruthHit pulse_record;
        pulse_record.tdc = tdc;
        pulse_record.num_electrons = fNumElectrons_v[index];
        pulse_record.tick = clockData.TPCTDC2Tick(tdc);

        double nElecADC = detProp.ElectronsToADC() * fNumElectrons_v[index];
        double eDeposit = 23.6 * fNumElectrons_v[index] / 1000.;    // This is a guesstimate for now but hopefully not far off, note we assume wirecell will handle recombination

        std::cout << "==> x position of plane 0: " << planeCoords[0] << ", normal: " << planeNormal[0] << ", nElec: " << fNumElectrons_v[index] << ", nElecADC: " << nElecADC << ", edep: " << eDeposit << std::endl;
        std::cout << "    x position of charge deposit: " << chargeDepCoords.X() << std::endl;

        simDep_v->emplace_back(0,
                               fNumElectrons_v[index],
                               0.,
                               eDeposit,
                               geo::Point_t(chargeDepCoords.X(),chargeDepCoords.Y()+0.001,chargeDepCoords.Z()+0.001),
                               geo::Point_t(chargeDepCoords.X(),chargeDepCoords.Y()-0.001,chargeDepCoords.Z()-0.001));

        geo::PlaneID collectionPlaneID(0,0,2);

        for(size_t plane=0; plane<3; ++plane) {
            geo::PlaneID planeID(0,0,plane);
            double       xyz[3] = {chargeDepCoords.X(),chargeDepCoords.Y(),chargeDepCoords.Z()};

            std::cout << "++ Searching for nearest wire id for planeID: " << planeID << ", xyz: "  << xyz[0] << "," << xyz[1] << "," << xyz[2] << std::endl;
            geo::WireID nearestID = geo->NearestWireID(xyz,plane);
            std::cout << "   WireID: " << nearestID << std::endl;
            std::cout << "** Searching for nearest channel with plane: " << plane << ", xyz: " << xyz[0] << "," << xyz[1] << "," << xyz[2] << std::endl;
            auto channel = geo->NearestChannel(xyz,plane);
            std::cout << "   Returned with channel: " << channel << std::endl;
            auto wire = geo->ChannelToWire(channel).front().Wire;
            std::cout << "   which gives wire: " << wire << std::endl;
            if(fVerbose) std::cout << "[BUFFOON!]    plane " << plane << " channel "
                << channel << " ... wire " << wire << std::endl;
            pulse_record.channel_list[plane] = channel;
            switch(plane) {
                case 0: fPlane0Channel_v.push_back(channel); fPlane0Wire_v.push_back(wire); break;
                case 1: fPlane1Channel_v.push_back(channel); fPlane1Wire_v.push_back(wire); break;
                case 2: fPlane2Channel_v.push_back(channel); fPlane2Wire_v.push_back(wire); break;
                default:
                    std::cerr << "[BUFFOON!] unexpected plane! " << plane << std::endl;
                    throw std::exception();
            }
            sim::SimChannel sch(channel);
            unsigned planeTDC = clockData.TPCTick2TDC(pulse_record.tick + detProp.GetXTicksOffset(planeID) - detProp.GetXTicksOffset(collectionPlaneID));
            sch.AddIonizationElectrons(1,   /// track id, keep 0 = invalid
                                       (unsigned int)planeTDC,
                                       fNumElectrons_v[index],
                                       xyz,
                                       100.); /// energy, keep -1 = invalid
            simch_v->emplace_back(std::move(sch));
        }
        pholder.Register(std::move(pulse_record));
    }

    _run = e.id().run();
    _subrun = e.id().subRun();
    _event = e.id().event();
    
    fTupleTree->Fill();
    e.put(std::move(simch_v));
    e.put(std::move(simDep_v));
    e.put(std::move(trigger_v));
    
    return;
}

Member Data Documentation

int SimTestPulse::_event

private

Definition at line 85 of file SimTestPulse_module.cc.

int SimTestPulse::_run

private

Definition at line 85 of file SimTestPulse_module.cc.

int SimTestPulse::_subrun

private

Definition at line 85 of file SimTestPulse_module.cc.

std::vector<double> SimTestPulse::fNumElectrons_v

private

Charge amount in electron count.

Definition at line 71 of file SimTestPulse_module.cc.

std::vector<int> SimTestPulse::fPlane0Channel_v

private

Resulting plane 0 channel number where charge is deposited.

Definition at line 77 of file SimTestPulse_module.cc.

std::vector<int> SimTestPulse::fPlane0Wire_v

private

Resulting plane 0 wire number where charge is deposited.

Definition at line 81 of file SimTestPulse_module.cc.

std::vector<int> SimTestPulse::fPlane1Channel_v

private

Resulting plane 1 channel number where charge is deposited.

Definition at line 78 of file SimTestPulse_module.cc.

std::vector<int> SimTestPulse::fPlane1Wire_v

private

Resulting plane 1 wire number where charge is deposited.

Definition at line 82 of file SimTestPulse_module.cc.

std::vector<int> SimTestPulse::fPlane2Channel_v

private

Resulting plane 2 channel number where charge is deposited.

Definition at line 79 of file SimTestPulse_module.cc.

std::vector<int> SimTestPulse::fPlane2Wire_v

private

Resulting plane 2 wire number where charge is deposited.

Definition at line 83 of file SimTestPulse_module.cc.

std::vector<double> SimTestPulse::fSimTime_v

private

Charge timing in electronics clock [us].

Definition at line 68 of file SimTestPulse_module.cc.

std::vector<int> SimTestPulse::fTick_v

private

Corresponding tick.

Definition at line 76 of file SimTestPulse_module.cc.

double SimTestPulse::fTriggerTime

private

Trigger timing in electronics clock [us].

Definition at line 67 of file SimTestPulse_module.cc.

TTree* SimTestPulse::fTupleTree

private

output analysis tree

Definition at line 87 of file SimTestPulse_module.cc.

bool SimTestPulse::fVerbose

private

verbosity

Definition at line 62 of file SimTestPulse_module.cc.

std::vector<double> SimTestPulse::fY_v

private

Charge Y position in detector coordinate [cm].

Definition at line 69 of file SimTestPulse_module.cc.

std::vector<double> SimTestPulse::fZ_v

private

Charge Z position in detector coordinate [cm].

Definition at line 70 of file SimTestPulse_module.cc.

---

The documentation for this class was generated from the following file:

• SimTestPulse_module.cc