Skip to content

Remove skip channel list and move to calibration database #899

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
asanchezcastillo wants to merge 1 commit into
base: develop
Choose a base branch
from
Open

Remove skip channel list and move to calibration database #899

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 11 additions & 2 deletions sbndcode/TPCPMTBarycenterMatching/TPCPMTBarycenterMatching_module.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -59,6 +59,10 @@
#include "lardataobj/RecoBase/PFParticleMetadata.h"
#include "sbnobj/Common/Reco/TPCPMTBarycenterMatch.h"

// Calibration database includes
#include "sbndcode/Calibration/PDSDatabaseInterface/PMTCalibrationDatabase.h"
#include "sbndcode/Calibration/PDSDatabaseInterface/IPMTCalibrationDatabaseService.h"

//ROOT includes
#include <Eigen/Dense>
#include <vector>
Expand Down Expand Up @@ -231,7 +235,6 @@ class TPCPMTBarycenterMatchProducer : public art::EDProducer {
bool fCollectionOnly; ///< Only use TPC spacepoints from the collection plane
double fDistanceCandidateFlashes; ///< Maximum distance between candidate flashes to be considered for matching (cm)
std::vector<double> fCalAreaConst; /// Calibration area constants for wire plane
std::vector<int> fSkipChannelList;
double fOpDetCoVUVEff; // Efficiencies for PMT detection (Coated PMT VUV)
double fOpDetCoVISEff; // Efficiencies for PMT detection (Coated PMT VIS)
double fOpDetUncoVISEff; // Efficiencies for PMT detection (Uncoated PMT VIS)
Expand Down Expand Up @@ -294,6 +297,9 @@ class TPCPMTBarycenterMatchProducer : public art::EDProducer {
fhicl::ParameterSet _vis_params;
std::shared_ptr<phot::OpticalPath> _optical_path_tool;

sbndDB::PMTCalibrationDatabase const* fPMTCalibrationDatabaseService;


};


Expand Down Expand Up @@ -391,6 +397,9 @@ TPCPMTBarycenterMatchProducer::TPCPMTBarycenterMatchProducer(fhicl::ParameterSet
_vis_params = p.get<fhicl::ParameterSet>("VIVHits");
_optical_path_tool = std::shared_ptr<phot::OpticalPath>(art::make_tool<phot::OpticalPath>(p.get<fhicl::ParameterSet>("OpticalPathTool")));
_semi_model = std::make_unique<phot::SemiAnalyticalModel>(_vuv_params, _vis_params, _optical_path_tool, true, false);

//Load PMT Calibration Database
fPMTCalibrationDatabaseService = lar::providerFrom<sbndDB::IPMTCalibrationDatabaseService const>();
}

void TPCPMTBarycenterMatchProducer::produce(art::Event& e)
Expand Down Expand Up @@ -892,7 +901,7 @@ double TPCPMTBarycenterMatchProducer::GetFlashLight(double flash_pe, std::vector
double tot_visibility=0;

for(size_t ch=0; ch<dir_visibility.size(); ch++){
if (std::find(fSkipChannelList.begin(), fSkipChannelList.end(), ch) != fSkipChannelList.end()) continue;
if(!fPMTCalibrationDatabaseService->getReconstructChannel(ch)) continue; // Skip channels not reconstructed
if(fOpDetType[ch]==0) tot_visibility += fOpDetCoVUVEff*dir_visibility[ch] + fOpDetCoVISEff*ref_visibility[ch];
else if(fOpDetType[ch]==1) tot_visibility += fOpDetUncoVISEff*ref_visibility[ch];
else continue; // skip other types
Expand Down