sim::GenericCRTUtility Class Reference

#include <GenericCRT.h>

Public Member Functions
sim::AuxDetIDE	toAuxDetIDE (sim::AuxDetHit const &InputHit) const
std::vector< unsigned int >	GetAuxDetChannels (const std::vector< sim::AuxDetHit > &InputHitVector) const
std::vector< sim::AuxDetSimChannel >	GetAuxDetSimChannels (const std::vector< sim::AuxDetHit > &InputHitVector) const
sim::AuxDetSimChannel	GetAuxDetSimChannelByNumber (const std::vector< sim::AuxDetHit > &InputHitVector, unsigned int inputchannel) const

Private Attributes
art::ServiceHandle< geo::Geometry const >	fGeo

Detailed Description

Definition at line 26 of file GenericCRT.h.

Member Function Documentation

std::vector< unsigned int > sim::GenericCRTUtility::GetAuxDetChannels

(

const std::vector< sim::AuxDetHit > &

InputHitVector

)

const

Definition at line 41 of file GenericCRT.cxx.

{

  std::vector<unsigned int> AuxDetChanNumber;
  AuxDetChanNumber.reserve(size(InputHitVector));


  for(auto const& hit : InputHitVector)
  {

  std::vector<unsigned int>::iterator Chanitr
      = std::find(AuxDetChanNumber.begin(), AuxDetChanNumber.end(), hit.GetID());

   if(Chanitr == AuxDetChanNumber.end()){ //If trackID is already in the map, update it
         //if channel ID is not in the set yet, add it
      AuxDetChanNumber.push_back(hit.GetID());
        }//

  }

return AuxDetChanNumber;

}

sim::AuxDetSimChannel sim::GenericCRTUtility::GetAuxDetSimChannelByNumber	(	const std::vector< sim::AuxDetHit > &	InputHitVector,
		unsigned int	inputchannel
	)		const

Definition at line 67 of file GenericCRT.cxx.

{
        std::vector<sim::AuxDetIDE> IDEvector;
        //loop over sim::AuxDetHits and assign them to AuxDetSimChannels.

        size_t ad_id_no = 9999;
        size_t ad_sen_id_no = 9999;

        for (auto const& auxDetHit : InputHitVector)
        {

            double xcoordinate = (auxDetHit.GetEntryX() + auxDetHit.GetExitX())/2.0;
            double ycoordinate = (auxDetHit.GetEntryY() + auxDetHit.GetExitY())/2.0;
            double zcoordinate = (auxDetHit.GetEntryZ() + auxDetHit.GetExitZ())/2.0;
            double worldPos[3] = {xcoordinate,ycoordinate,zcoordinate};
            fGeo->FindAuxDetSensitiveAtPosition(worldPos, ad_id_no, ad_sen_id_no, 0.0001);
            if(auxDetHit.GetID() == inputchannel)   // this is the channel we want.
            {
                auto tempIDE = toAuxDetIDE(auxDetHit);

                std::vector<sim::AuxDetIDE>::iterator IDEitr
                = std::find(IDEvector.begin(), IDEvector.end(), tempIDE);

                if(IDEitr != IDEvector.end()){ //If trackID is already in the map, update it
                    //Andrzej's note - following logic from AuxDetReadout in Legacy, but why are the other paremeters getting overwritten like that?
                    IDEitr->energyDeposited += tempIDE.energyDeposited;
                    IDEitr->exitX            = tempIDE.exitX;
                    IDEitr->exitY            = tempIDE.exitY;
                    IDEitr->exitZ            = tempIDE.exitZ;
                    IDEitr->exitT            = tempIDE.exitT;
                    IDEitr->exitMomentumX    = tempIDE.exitMomentumX;
                    IDEitr->exitMomentumY    = tempIDE.exitMomentumY;
                    IDEitr->exitMomentumZ    = tempIDE.exitMomentumZ;
                    }
                    else{  //if trackID is not in the set yet, add it
                        IDEvector.push_back(std::move(tempIDE));
                    }//else

                    break;
            } // end if the AuxDetHit channel checks out.

        } // end main loop on AuxDetHit

        //push back the AuxDetSimChannel Vector.
        //TODO check the last parameter values.
        return sim::AuxDetSimChannel(ad_id_no, std::move(IDEvector), ad_sen_id_no);
}

std::vector< sim::AuxDetSimChannel > sim::GenericCRTUtility::GetAuxDetSimChannels

(

const std::vector< sim::AuxDetHit > &

InputHitVector

)

const

Definition at line 118 of file GenericCRT.cxx.

{
  auto const auxDetChannels = GetAuxDetChannels(InputHitVector);
  std::vector<sim::AuxDetSimChannel> auxDetVector;
  auxDetVector.reserve(size(auxDetChannels));

  for (auto const channelNum : auxDetChannels) {
    auxDetVector.push_back(GetAuxDetSimChannelByNumber(InputHitVector, channelNum));
  }

  return auxDetVector;
}

sim::AuxDetIDE sim::GenericCRTUtility::toAuxDetIDE

(

sim::AuxDetHit const &

InputHit

)

const

Title: GenericCRT Utility Class Author: Andrzej Szelc (andrz.nosp@m.ejs@.nosp@m.fnal..nosp@m.gov)

Description: Class with Algorithms to convert sim::AuxDetHits to sim::AuxDetSimChannels

Definition at line 17 of file GenericCRT.cxx.

  {
    sim::AuxDetIDE outputIDE;

   outputIDE.trackID                = InputHit.GetTrackID();
   outputIDE.energyDeposited    = InputHit.GetEnergyDeposited();
   outputIDE.entryX             = InputHit.GetEntryX();
   outputIDE.entryY             = InputHit.GetEntryY();
   outputIDE.entryZ             = InputHit.GetEntryZ();
   outputIDE.entryT             = InputHit.GetEntryT();
   outputIDE.exitX              = InputHit.GetExitX();
   outputIDE.exitY              = InputHit.GetExitY();
   outputIDE.exitZ              = InputHit.GetExitZ();
   outputIDE.exitT              = InputHit.GetExitT();
   outputIDE.exitMomentumX      = InputHit.GetExitMomentumX();
   outputIDE.exitMomentumY      = InputHit.GetExitMomentumY();
   outputIDE.exitMomentumZ      = InputHit.GetExitMomentumZ();


   return outputIDE;
}

Member Data Documentation

art::ServiceHandle<geo::Geometry const> sim::GenericCRTUtility::fGeo

private

Definition at line 40 of file GenericCRT.h.

---

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

• larsim/larsim/GenericCRT/GenericCRT.h
• larsim/larsim/GenericCRT/GenericCRT.cxx