opdet::OpDetAnalyzer Class Reference

Inheritance diagram for opdet::OpDetAnalyzer:

Public Member Functions
	OpDetAnalyzer (fhicl::ParameterSet const &p)
	OpDetAnalyzer (OpDetAnalyzer const &)=delete
	OpDetAnalyzer (OpDetAnalyzer &&)=delete
OpDetAnalyzer &	operator= (OpDetAnalyzer const &)=delete
OpDetAnalyzer &	operator= (OpDetAnalyzer &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
void	endJob () override

Private Attributes
std::vector< std::string >	fInputModule
std::unique_ptr< opdet::PDMapAlg >	fPDSMapPtr
std::vector< std::string >	fPMTMapLabel
std::vector< std::string >	fArapucaMapLabel
std::vector< int >	fNOpDetAll
std::vector< int >	fNOpDetDirect
std::vector< int >	fNOpDetReflected
TTree *	fAllPhotonsTree
TTree *	fTheOpDetTree
TTree *	fTheEventTree
int	fEventID
int	fOpChannel
float	fWavelength
float	fTime
int	fCountOpDetAll
int	fCountOpDetDirect
int	fCountOpDetReflected
float	fOpDetX
float	fOpDetY
float	fOpDetZ
bool	fisPMT
int	fCountEventAll
int	fCountEventDirect
int	fCountEventReflected
int	fVerbosity

Static Private Attributes
static constexpr double	kVUVWavelength = 128.0
	Value used when a typical ultraviolet light wavelength is needed. More...
static constexpr double	kVisibleWavelength = 450.0
	Value used when a typical visible light wavelength is needed. More...

Detailed Description

Definition at line 39 of file OpDetAnalyzer_module.cc.

Constructor & Destructor Documentation

opdet::OpDetAnalyzer::OpDetAnalyzer ( fhicl::ParameterSet const &  p )

explicit

Definition at line 100 of file OpDetAnalyzer_module.cc.

  : EDAnalyzer{pset},
    fPDSMapPtr{art::make_tool<opdet::PDMapAlg>(pset.get<fhicl::ParameterSet>("PDSMapTool"))}
{
  
  fVerbosity = pset.get<int>("Verbosity");
  
  try
  {
     fInputModule = pset.get<std::vector<std::string>>("InputModule");
  }
  catch(...)
  {
     fInputModule.push_back(pset.get<std::string>("InputModule"));
  }
}

opdet::OpDetAnalyzer::OpDetAnalyzer ( OpDetAnalyzer const &  )

delete

opdet::OpDetAnalyzer::OpDetAnalyzer ( OpDetAnalyzer &&  )

delete

Member Function Documentation

void opdet::OpDetAnalyzer::analyze ( art::Event const &  e )

override

Definition at line 150 of file OpDetAnalyzer_module.cc.

{

  // geometry information
  art::ServiceHandle<geo::Geometry> geo;

  // get event number
  fEventID = evt.id().event();
  
  // adapted from standard SimPhotonCounter code

  // Get *ALL* SimPhotonsCollection from Event
  auto photon_handles = evt.getMany<std::vector<sim::SimPhotonsLite>>();
  if (photon_handles.size() == 0)
    throw art::Exception(art::errors::ProductNotFound)<<"sim SimPhotons retrieved and you requested them.";

  // reset counters
  // event
  fCountEventAll=0;
  fCountEventDirect=0;
  fCountEventReflected=0;
  // opdet
  fCountOpDetAll=0;
  fCountOpDetDirect=0;
  fCountOpDetReflected=0;

  // reset vectors
  fNOpDetAll.clear(); fNOpDetAll.resize(geo->NOpChannels(), 0);
  fNOpDetDirect.clear(); fNOpDetDirect.resize(geo->NOpChannels(), 0);
  fNOpDetReflected.clear(); fNOpDetReflected.resize(geo->NOpChannels(), 0);

  // Get SimPhotonsLite from Event
  for(auto const& mod : fInputModule){

    // Loop over direct/reflected photons
    for (auto const& ph_handle: photon_handles) {
      // Do some checking before we proceed
      if (!ph_handle.isValid()) continue;
      if (ph_handle.provenance()->moduleLabel() != mod) continue;   //not the most efficient way of doing this, but preserves the logic of the module. Andrzej

      bool Reflected = (ph_handle.provenance()->productInstanceName() == "Reflected");
      
      if(fVerbosity > 0) std::cout<<"Found OpDet hit collection of size "<< (*ph_handle).size()<<std::endl;

      if((*ph_handle).size()>0)
      {

        for ( auto const& photon : (*ph_handle) )
        {
          //Get data from HitCollection entry
          fOpChannel=photon.OpChannel;
          std::map<int, int> PhotonsMap = photon.DetectedPhotons;

          // do not save if PD is not sensitive to this light
          std::string pd_type=fPDSMapPtr->pdType(fOpChannel);
          if (Reflected && pd_type=="xarapuca_vuv") continue;
          if(!Reflected && (pd_type=="xarapuca_vis" || pd_type=="pmt_uncoated")) continue;

          for(auto it = PhotonsMap.begin(); it!= PhotonsMap.end(); it++)
          {

            // Calculate wavelength in nm
            if (Reflected) fWavelength = kVisibleWavelength;
            else fWavelength= kVUVWavelength;

            // Get arrival time from phot
            fTime= it->first;

            for(int i = 0; i < it->second ; i++)
            {
              // Increment per OpDet counters and fill all photon tree
              fCountOpDetAll++;
              fAllPhotonsTree->Fill();

              // all
              fNOpDetAll[fOpChannel]++;
  
              // direct light
              if (!Reflected){
                fNOpDetDirect[fOpChannel]++;
              }
              // reflected light
              else if (Reflected) {
                fNOpDetReflected[fOpChannel]++;
              }            
            }
          }
        } // eng loop over photons        
      }
    } // end loop over photon handles
  } // end loop over input modules


  // fill information in perOpDet and perEvent trees
  for (unsigned int OpDet = 0; OpDet < geo->NOpChannels(); OpDet++) {
    
    // perOpDet tree
    fOpChannel = OpDet;
    geo::Point_t OpDetCenter = geo->OpDetGeoFromOpDet(OpDet).GetCenter();
    fOpDetX = OpDetCenter.X();
    fOpDetY = OpDetCenter.Y();
    fOpDetZ = OpDetCenter.Z();
    fisPMT =  geo->OpDetGeoFromOpDet(OpDet).isSphere();
    fCountOpDetAll = fNOpDetAll[OpDet];
    fCountOpDetDirect = fNOpDetDirect[OpDet];
    fCountOpDetReflected = fNOpDetReflected[OpDet];
    fTheOpDetTree->Fill();

    // increment counters for perEvent tree
    fCountEventAll+=fCountOpDetAll;
    fCountEventDirect+=fCountOpDetDirect;
    fCountEventReflected+=fCountOpDetReflected;

     if(fVerbosity >2) std::cout<<"OpDetResponseInterface PerOpDet : Event "<<fEventID<<" OpDet " << fOpChannel << " All " << fCountOpDetAll << " Dir " <<fCountOpDetDirect << ", Refl " << fCountOpDetReflected<<std::endl;
  }

  // perEvent tree
  fTheEventTree->Fill();

  if(fVerbosity >1) std::cout<<"OpDetResponseInterface PerEvent : Event "<<fEventID<<" All " << fCountEventAll << " Dir " <<fCountEventDirect << ", Refl " << fCountEventReflected <<std::endl;

}

void opdet::OpDetAnalyzer::beginJob ( )

override

Definition at line 117 of file OpDetAnalyzer_module.cc.

{
  // create output trees
  art::ServiceHandle<art::TFileService> tfs;
  
  // all photons tree
  fAllPhotonsTree = tfs->make<TTree>("AllPhotons", "AllPhotons");
  fAllPhotonsTree->Branch("EventID", &fEventID, "EventID/I");
  fAllPhotonsTree->Branch("Wavelength", &fWavelength, "Wavelength/F");
  fAllPhotonsTree->Branch("OpChannel", &fOpChannel, "OpChannel/I");
  fAllPhotonsTree->Branch("Time", &fTime, "Time/F");
  
  // photons per opdet tree
  fTheOpDetTree = tfs->make<TTree>("PhotonsPerOpDet","PhotonsPerOpDet");
  fTheOpDetTree->Branch("EventID", &fEventID, "EventID/I");
  fTheOpDetTree->Branch("OpChannel", &fOpChannel, "OpChannel/I");
  fTheOpDetTree->Branch("OpDetX", &fOpDetX, "OpDetX/F");
  fTheOpDetTree->Branch("OpDetY", &fOpDetY, "OpDetY/F");
  fTheOpDetTree->Branch("OpDetZ", &fOpDetZ, "OpDetZ/F");
  fTheOpDetTree->Branch("isPMT", &fisPMT);
  fTheOpDetTree->Branch("CountAll", &fCountOpDetAll, "CountAll/I");
  fTheOpDetTree->Branch("CountDirect", &fCountOpDetDirect, "CountDirect/I");
  fTheOpDetTree->Branch("CountReflected", &fCountOpDetReflected, "CountReflected/I");

  // photons per event tree
  fTheEventTree  = tfs->make<TTree>("PhotonsPerEvent","PhotonsPerEvent");
  fTheEventTree->Branch("EventID", &fEventID, "EventID/I");
  fTheEventTree->Branch("CountAll", &fCountEventAll, "CountAll/I");
  fTheEventTree->Branch("CountDirect", &fCountEventDirect, "CountDirect/I");
  fTheEventTree->Branch("CountReflected", &fCountEventReflected, "CountReflected/I");

}

void opdet::OpDetAnalyzer::endJob ( )

override

Definition at line 273 of file OpDetAnalyzer_module.cc.

{
  // Implementation of optional member function here.
}

OpDetAnalyzer& opdet::OpDetAnalyzer::operator= ( OpDetAnalyzer const &  )

delete

OpDetAnalyzer& opdet::OpDetAnalyzer::operator= ( OpDetAnalyzer &&  )

delete

Member Data Documentation

TTree* opdet::OpDetAnalyzer::fAllPhotonsTree

private

Definition at line 70 of file OpDetAnalyzer_module.cc.

std::vector<std::string> opdet::OpDetAnalyzer::fArapucaMapLabel

private

Definition at line 64 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fCountEventAll

private

Definition at line 85 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fCountEventDirect

private

Definition at line 86 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fCountEventReflected

private

Definition at line 87 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fCountOpDetAll

private

Definition at line 79 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fCountOpDetDirect

private

Definition at line 80 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fCountOpDetReflected

private

Definition at line 81 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fEventID

private

Definition at line 74 of file OpDetAnalyzer_module.cc.

std::vector<std::string> opdet::OpDetAnalyzer::fInputModule

private

Definition at line 60 of file OpDetAnalyzer_module.cc.

bool opdet::OpDetAnalyzer::fisPMT

private

Definition at line 83 of file OpDetAnalyzer_module.cc.

std::vector<int> opdet::OpDetAnalyzer::fNOpDetAll

private

Definition at line 66 of file OpDetAnalyzer_module.cc.

std::vector<int> opdet::OpDetAnalyzer::fNOpDetDirect

private

Definition at line 67 of file OpDetAnalyzer_module.cc.

std::vector<int> opdet::OpDetAnalyzer::fNOpDetReflected

private

Definition at line 68 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fOpChannel

private

Definition at line 75 of file OpDetAnalyzer_module.cc.

float opdet::OpDetAnalyzer::fOpDetX

private

Definition at line 82 of file OpDetAnalyzer_module.cc.

float opdet::OpDetAnalyzer::fOpDetY

private

Definition at line 82 of file OpDetAnalyzer_module.cc.

float opdet::OpDetAnalyzer::fOpDetZ

private

Definition at line 82 of file OpDetAnalyzer_module.cc.

std::unique_ptr<opdet::PDMapAlg> opdet::OpDetAnalyzer::fPDSMapPtr

private

Definition at line 63 of file OpDetAnalyzer_module.cc.

std::vector<std::string> opdet::OpDetAnalyzer::fPMTMapLabel

private

Definition at line 64 of file OpDetAnalyzer_module.cc.

TTree* opdet::OpDetAnalyzer::fTheEventTree

private

Definition at line 72 of file OpDetAnalyzer_module.cc.

TTree* opdet::OpDetAnalyzer::fTheOpDetTree

private

Definition at line 71 of file OpDetAnalyzer_module.cc.

float opdet::OpDetAnalyzer::fTime

private

Definition at line 77 of file OpDetAnalyzer_module.cc.

int opdet::OpDetAnalyzer::fVerbosity

private

Definition at line 89 of file OpDetAnalyzer_module.cc.

float opdet::OpDetAnalyzer::fWavelength

private

Definition at line 76 of file OpDetAnalyzer_module.cc.

constexpr double opdet::OpDetAnalyzer::kVisibleWavelength = 450.0

staticprivate

Value used when a typical visible light wavelength is needed.

Definition at line 95 of file OpDetAnalyzer_module.cc.

constexpr double opdet::OpDetAnalyzer::kVUVWavelength = 128.0

staticprivate

Value used when a typical ultraviolet light wavelength is needed.

Definition at line 92 of file OpDetAnalyzer_module.cc.

---

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

• OpDetAnalyzer_module.cc