shower::TCShowerElectronLikelihood Class Reference

Inheritance diagram for shower::TCShowerElectronLikelihood:

Public Member Functions
	TCShowerElectronLikelihood (fhicl::ParameterSet const &pset)
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Member Functions
void	beginJob () override
void	analyze (const art::Event &evt) override
void	getShowerProfile (detinfo::DetectorClocksData const &clockdata, detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit >> showerhits, TVector3 shwvtx, TVector3 shwdir)
void	findEnergyBin ()
double	getLongLikelihood ()
double	getTranLikelihood ()
void	resetProfiles ()

Private Attributes
std::string	fTemplateFile
std::string	fROOTfile
TH3F *	longTemplate
TH3F *	tranTemplate
TH3F *	tranTemplate_1
TH3F *	tranTemplate_2
TH3F *	tranTemplate_3
TH3F *	tranTemplate_4
TH3F *	tranTemplate_5
TProfile2D *	longTemplateProf2D
TProfile2D *	tranTemplateProf2D
TProfile2D *	tranTemplateProf2D_1
TProfile2D *	tranTemplateProf2D_2
TProfile2D *	tranTemplateProf2D_3
TProfile2D *	tranTemplateProf2D_4
TProfile2D *	tranTemplateProf2D_5
TH1F *	energyDist
TH1F *	longLikelihoodHist
TH1F *	tranLikelihoodHist
TH1F *	longProfHist
TH1F *	tranProfHist_1
TH1F *	tranProfHist_2
TH1F *	tranProfHist_3
TH1F *	tranProfHist_4
TH1F *	tranProfHist_5
TH1F *	longProfile
TH1F *	tranProfile
TH1F *	tranProfile_1
TH1F *	tranProfile_2
TH1F *	tranProfile_3
TH1F *	tranProfile_4
TH1F *	tranProfile_5
int	energyGuess
double	energyChi2
int	maxt
const int	LBINS = 20
const int	LMIN = 0
const int	LMAX = 5
const int	TBINS = 20
const int	TMIN = -5
const int	TMAX = 5
const double	X0 = 14
std::string	fHitModuleLabel
std::string	fShowerModuleLabel
std::string	fTemplateModuleLabel
std::string	fGenieGenModuleLabel
std::string	fDigitModuleLabel
calo::CalorimetryAlg	fCalorimetryAlg

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 37 of file TCShowerElectronLikelihood_module.cc.

Constructor & Destructor Documentation

shower::TCShowerElectronLikelihood::TCShowerElectronLikelihood ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 122 of file TCShowerElectronLikelihood_module.cc.

  : EDAnalyzer(pset)
  , fHitModuleLabel(pset.get<std::string>("HitModuleLabel", "trajcluster"))
  , fShowerModuleLabel(pset.get<std::string>("ShowerModuleLabel", "tcshower"))
  , fGenieGenModuleLabel(pset.get<std::string>("GenieGenModuleLabel", "generator"))
  , fCalorimetryAlg(pset.get<fhicl::ParameterSet>("CalorimetryAlg"))
{
  fTemplateFile = pset.get<std::string>("TemplateFile");
  cet::search_path sp("FW_SEARCH_PATH");
  if (!sp.find_file(fTemplateFile, fROOTfile))
    throw cet::exception("TCShowerElectronLikelihood")
      << "cannot find the root template file: \n"
      << fTemplateFile << "\n bail ungracefully.\n";

  std::unique_ptr<TFile> file{TFile::Open(fROOTfile.c_str())};

  longTemplate = file->Get<TH3F>("tcshowertemplate/fLongitudinal");
  tranTemplate = file->Get<TH3F>("tcshowertemplate/fTransverse");
  tranTemplate_1 = file->Get<TH3F>("tcshowertemplate/fTransverse_1");
  tranTemplate_2 = file->Get<TH3F>("tcshowertemplate/fTransverse_2");
  tranTemplate_3 = file->Get<TH3F>("tcshowertemplate/fTransverse_3");
  tranTemplate_4 = file->Get<TH3F>("tcshowertemplate/fTransverse_4");
  tranTemplate_5 = file->Get<TH3F>("tcshowertemplate/fTransverse_5");

  longTemplateProf2D = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoLong2D");
  tranTemplateProf2D = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D");
  tranTemplateProf2D_1 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_1");
  tranTemplateProf2D_2 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_2");
  tranTemplateProf2D_3 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_3");
  tranTemplateProf2D_4 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_4");
  tranTemplateProf2D_5 = file->Get<TProfile2D>("tcshowertemplate/fShowerProfileRecoTrans2D_5");

  longProfile = new TH1F("longProfile", "longitudinal shower profile;t;Q", LBINS, LMIN, LMAX);
  tranProfile = new TH1F("tranProfile", "transverse shower profile;dist (cm);Q", TBINS, TMIN, TMAX);

  tranProfile_1 = new TH1F(
    "tranProfile_1", "transverse shower profile [0 <= t < 1];dist (cm);Q", TBINS, TMIN, TMAX);
  tranProfile_2 = new TH1F(
    "tranProfile_2", "transverse shower profile [1 <= t < 2];dist (cm);Q", TBINS, TMIN, TMAX);
  tranProfile_3 = new TH1F(
    "tranProfile_3", "transverse shower profile [2 <= t < 3];dist (cm);Q", TBINS, TMIN, TMAX);
  tranProfile_4 = new TH1F(
    "tranProfile_4", "transverse shower profile [3 <= t < 4];dist (cm);Q", TBINS, TMIN, TMAX);
  tranProfile_5 = new TH1F(
    "tranProfile_5", "transverse shower profile [4 <= t < 5];dist (cm);Q", TBINS, TMIN, TMAX);
} // TCShowerElectronLikelihood

Member Function Documentation

void shower::TCShowerElectronLikelihood::analyze ( const art::Event &  evt )

overrideprivate

Definition at line 200 of file TCShowerElectronLikelihood_module.cc.

{
  resetProfiles();

  art::Handle<std::vector<recob::Hit>> hitListHandle;
  std::vector<art::Ptr<recob::Hit>> hitlist;
  if (evt.getByLabel(fHitModuleLabel, hitListHandle)) art::fill_ptr_vector(hitlist, hitListHandle);

  art::Handle<std::vector<recob::Shower>> showerListHandle;
  std::vector<art::Ptr<recob::Shower>> showerlist;
  if (evt.getByLabel(fShowerModuleLabel, showerListHandle))
    art::fill_ptr_vector(showerlist, showerListHandle);

  art::Handle<std::vector<simb::MCTruth>> mctruthListHandle;
  std::vector<art::Ptr<simb::MCTruth>> mclist;
  if (evt.getByLabel(fGenieGenModuleLabel, mctruthListHandle))
    art::fill_ptr_vector(mclist, mctruthListHandle);

  art::FindManyP<recob::Hit> shwfm(showerListHandle, evt, fShowerModuleLabel);

  if (showerlist.size()) {
    auto const clock_data =
      art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
    auto const det_prop =
      art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clock_data);
    std::vector<art::Ptr<recob::Hit>> showerhits = shwfm.at(0);
    getShowerProfile(
      clock_data, det_prop, showerhits, showerlist[0]->ShowerStart(), showerlist[0]->Direction());

    maxt = std::ceil((90 - showerlist[0]->ShowerStart().Z()) / X0);

    findEnergyBin();

    longLikelihoodHist->Fill(getLongLikelihood());
    tranLikelihoodHist->Fill(getTranLikelihood());

    // check true shower energy
    if (mclist.size()) {
      art::Ptr<simb::MCTruth> mctruth = mclist[0];
      if (mctruth->NeutrinoSet()) {
        if (std::abs(mctruth->GetNeutrino().Nu().PdgCode()) == 12 &&
            mctruth->GetNeutrino().CCNC() == 0) {
          double elep = mctruth->GetNeutrino().Lepton().E();
          energyDist->Fill(elep - energyGuess);
        } // cc nue
      }   // neutrinoset
    }     // mclist
  }       // showerlist

  longProfHist = longProfile;
  tranProfHist_1 = tranProfile_1;
  tranProfHist_2 = tranProfile_2;
  tranProfHist_3 = tranProfile_3;
  tranProfHist_4 = tranProfile_4;
  tranProfHist_5 = tranProfile_5;
} // analyze

void shower::TCShowerElectronLikelihood::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 172 of file TCShowerElectronLikelihood_module.cc.

{

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

  energyDist = tfs->make<TH1F>("energyDist", "true energy - guess energy", 41, -20.5, 20.5);
  longLikelihoodHist = tfs->make<TH1F>("longLikelihoodHist", "longitudinal likelihood", 20, 0, 2);
  tranLikelihoodHist = tfs->make<TH1F>("tranLikelihoodHist", "transverse likelihood", 20, 0, 3);

  // just for printing purposes
  longProfHist =
    tfs->make<TH1F>("longProfHist", "longitudinal e- profile (reco);t;Q", LBINS, LMIN, LMAX);
  tranProfHist_1 = tfs->make<TH1F>(
    "tranProfHist", "transverse e- profile (reco) [0 <= t < 1];dist (cm);Q", TBINS, TMIN, TMAX);
  tranProfHist_2 = tfs->make<TH1F>(
    "tranProfHist", "transverse e- profile (reco) [1 <= t < 2];dist (cm);Q", TBINS, TMIN, TMAX);
  tranProfHist_3 = tfs->make<TH1F>(
    "tranProfHist", "transverse e- profile (reco) [2 <= t < 3];dist (cm);Q", TBINS, TMIN, TMAX);
  tranProfHist_4 = tfs->make<TH1F>(
    "tranProfHist", "transverse e- profile (reco) [3 <= t < 4];dist (cm);Q", TBINS, TMIN, TMAX);
  tranProfHist_5 = tfs->make<TH1F>(
    "tranProfHist", "transverse e- profile (reco) [4 <= t < 5];dist (cm);Q", TBINS, TMIN, TMAX);

} // beginJob

void shower::TCShowerElectronLikelihood::findEnergyBin ( )

private

Definition at line 353 of file TCShowerElectronLikelihood_module.cc.

{

  if (longProfile->GetNbinsX() != longTemplate->GetNbinsX())
    throw cet::exception("TCShowerElectronLikelihood")
      << "Bin mismatch in longitudinal profile template \n";

  if (tranProfile->GetNbinsX() != tranTemplate->GetNbinsX())
    throw cet::exception("TCShowerElectronLikelihood")
      << "Bin mismatch in transverse profile template \n";

  double chi2min = 999999;
  double bestbin = -1;

  int ebins = longTemplate->GetNbinsY();
  int lbins = longTemplate->GetNbinsX();
  int tbins = tranTemplate->GetNbinsX();

  for (int i = 0; i < ebins; ++i) {
    double thischi2 = 0;

    TProfile* ltemp = longTemplateProf2D->ProfileX("_x", i + 1, i + 1);
    TProfile* ttemp_1 = tranTemplateProf2D_1->ProfileX("_x_1", i + 1, i + 1);
    TProfile* ttemp_2 = tranTemplateProf2D_2->ProfileX("_x_2", i + 1, i + 1);
    TProfile* ttemp_3 = tranTemplateProf2D_3->ProfileX("_x_3", i + 1, i + 1);
    TProfile* ttemp_4 = tranTemplateProf2D_4->ProfileX("_x_4", i + 1, i + 1);
    TProfile* ttemp_5 = tranTemplateProf2D_5->ProfileX("_x_5", i + 1, i + 1);

    int nlbins = 0;
    int ntbins = 0;

    for (int j = 0; j < lbins; ++j) {
      double obs = longProfile->GetBinContent(j + 1);
      double exp = ltemp->GetBinContent(j + 1);
      if (obs != 0) {
        thischi2 += cet::square(obs - exp) / exp;
        ++nlbins;
      }
    } // loop through longitudinal bins

    for (int j = 0; j < tbins; ++j) {
      double obs = tranProfile_1->GetBinContent(j + 1);
      double exp = ttemp_1->GetBinContent(j + 1);
      if (obs != 0) {
        thischi2 += cet::square(obs - exp) / exp;
        ++ntbins;
      }

      obs = tranProfile_2->GetBinContent(j + 1);
      exp = ttemp_2->GetBinContent(j + 1);
      if (obs != 0) {
        thischi2 += cet::square(obs - exp) / exp;
        ++ntbins;
      }

      obs = tranProfile_3->GetBinContent(j + 1);
      exp = ttemp_3->GetBinContent(j + 1);
      if (obs != 0) {
        thischi2 += cet::square(obs - exp) / exp;
        ++ntbins;
      }

      obs = tranProfile_4->GetBinContent(j + 1);
      exp = ttemp_4->GetBinContent(j + 1);
      if (obs != 0) {
        thischi2 += cet::square(obs - exp) / exp;
        ++ntbins;
      }

      obs = tranProfile_5->GetBinContent(j + 1);
      exp = ttemp_5->GetBinContent(j + 1);
      if (obs != 0) {
        thischi2 += cet::square(obs - exp) / exp;
        ++ntbins;
      }
    } // loop through longitudinal bins

    thischi2 /= (nlbins + ntbins);

    if (thischi2 < chi2min) {
      chi2min = thischi2;
      bestbin = i;
    }

  } // loop through energy bins

  energyChi2 = chi2min;
  energyGuess = bestbin + 1;
} // findEnergyBin

double shower::TCShowerElectronLikelihood::getLongLikelihood ( )

private

Definition at line 446 of file TCShowerElectronLikelihood_module.cc.

{
  if (energyGuess < 0) return -9999.;
  int energyBin = energyGuess;

  double longLikelihood = 0.;
  int nbins = 0;

  for (int i = 0; i < LBINS; ++i) {
    double qval = longProfile->GetBinContent(i + 1);
    int qbin = longTemplate->GetZaxis()->FindBin(qval);
    int binentries = longTemplate->GetBinContent(i + 1, energyBin, qbin);
    int totentries = longTemplate->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
    if (qval > 0) {
      ++nbins;
      double prob = (double)binentries / totentries * 100;
      if (binentries > 0) longLikelihood += log(prob);
    }
  } // loop through

  longLikelihood /= nbins;

  std::cout << longLikelihood << std::endl;
  return longLikelihood;
} // getLongLikelihood

void shower::TCShowerElectronLikelihood::getShowerProfile ( detinfo::DetectorClocksData const &  clockdata, detinfo::DetectorPropertiesData const &  detProp, std::vector< art::Ptr< recob::Hit >>  showerhits, TVector3  shwvtx, TVector3  shwdir  )

private

Definition at line 279 of file TCShowerElectronLikelihood_module.cc.

{
  art::ServiceHandle<geo::Geometry const> geom;

  auto collectionPlane = geo::PlaneID(0, 0, 1);

  double shwVtxTime = detProp.ConvertXToTicks(shwvtx[0], collectionPlane);
  double shwVtxWire = geom->WireCoordinate(shwvtx[1], shwvtx[2], collectionPlane);

  double shwTwoTime = detProp.ConvertXToTicks(shwvtx[0] + shwdir[0], collectionPlane);
  double shwTwoWire =
    geom->WireCoordinate(shwvtx[1] + shwdir[1], shwvtx[2] + shwdir[2], collectionPlane);

  for (size_t i = 0; i < showerhits.size(); ++i) {
    if (showerhits[i]->WireID().Plane != collectionPlane.Plane) continue;

    double wirePitch = geom->WirePitch(showerhits[i]->WireID());
    double tickToDist = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
    tickToDist *= 1.e-3 * sampling_rate(clockdata); // 1e-3 is conversion of 1/us to 1/ns

    double xvtx = shwVtxTime * tickToDist;
    double yvtx = shwVtxWire * wirePitch;

    double xtwo = shwTwoTime * tickToDist;
    double ytwo = shwTwoWire * wirePitch;

    double xtwoorth = (ytwo - yvtx) + xvtx;
    double ytwoorth = -(xtwo - xvtx) + yvtx;

    double xhit = showerhits[i]->PeakTime() * tickToDist;
    double yhit = showerhits[i]->WireID().Wire * wirePitch;

    double ldist = std::abs((ytwoorth - yvtx) * xhit - (xtwoorth - xvtx) * yhit + xtwoorth * yvtx -
                            ytwoorth * xvtx) /
                   std::hypot(ytwoorth - yvtx, xtwoorth - xvtx);
    double tdist = ((ytwo - yvtx) * xhit - (xtwo - xvtx) * yhit + xtwo * yvtx - ytwo * xvtx) /
                   std::hypot(ytwo - yvtx, xtwo - xvtx);

    double to3D = 1. / std::hypot(xvtx - xtwo,
                                  yvtx - ytwo); // distance between two points in 3D space is one
    ldist *= to3D;
    tdist *= to3D;
    double t = ldist / X0;

    double Q = showerhits[i]->Integral() *
               fCalorimetryAlg.LifetimeCorrection(clockdata, detProp, showerhits[i]->PeakTime());

    longProfile->Fill(t, Q);
    tranProfile->Fill(tdist, Q);

    if (t < 1)
      tranProfile_1->Fill(tdist, Q);
    else if (t < 2)
      tranProfile_2->Fill(tdist, Q);
    else if (t < 3)
      tranProfile_3->Fill(tdist, Q);
    else if (t < 4)
      tranProfile_4->Fill(tdist, Q);
    else if (t < 5)
      tranProfile_5->Fill(tdist, Q);

  } // loop through showerhits

  return;

} // getShowerProfile

double shower::TCShowerElectronLikelihood::getTranLikelihood ( )

private

Definition at line 475 of file TCShowerElectronLikelihood_module.cc.

{
  if (energyGuess < 0) return -9999.;
  int energyBin = energyGuess;

  double tranLikelihood_1 = 0;
  double tranLikelihood_2 = 0;
  double tranLikelihood_3 = 0;
  double tranLikelihood_4 = 0;
  double tranLikelihood_5 = 0;

  double qval;
  int qbin, binentries, totentries;

  int nbins = 0;

  for (int i = 0; i < TBINS; ++i) {
    qval = tranProfile_1->GetBinContent(i + 1);
    qbin = tranTemplate_1->GetZaxis()->FindBin(qval);
    binentries = tranTemplate_1->GetBinContent(i + 1, energyBin, qbin);
    totentries = tranTemplate_1->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
    if (qval > 0) {
      ++nbins;
      double prob = (double)binentries / totentries * 100;
      if (binentries > 0) tranLikelihood_1 += log(prob);
    }

    qval = tranProfile_2->GetBinContent(i + 1);
    qbin = tranTemplate_2->GetZaxis()->FindBin(qval);
    binentries = tranTemplate_2->GetBinContent(i + 1, energyBin, qbin);
    totentries = tranTemplate_2->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
    if (qval > 0) {
      ++nbins;
      double prob = (double)binentries / totentries * 100;
      if (binentries > 0) tranLikelihood_2 += log(prob);
    }

    qval = tranProfile_3->GetBinContent(i + 1);
    qbin = tranTemplate_3->GetZaxis()->FindBin(qval);
    binentries = tranTemplate_3->GetBinContent(i + 1, energyBin, qbin);
    totentries = tranTemplate_3->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
    if (qval > 0) {
      ++nbins;
      double prob = (double)binentries / totentries * 100;
      if (binentries > 0) tranLikelihood_3 += log(prob);
    }

    qval = tranProfile_4->GetBinContent(i + 1);
    qbin = tranTemplate_4->GetZaxis()->FindBin(qval);
    binentries = tranTemplate_4->GetBinContent(i + 1, energyBin, qbin);
    totentries = tranTemplate_4->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
    if (qval > 0) {
      ++nbins;
      double prob = (double)binentries / totentries * 100;
      if (binentries > 0) tranLikelihood_4 += log(prob);
    }

    qval = tranProfile_5->GetBinContent(i + 1);
    qbin = tranTemplate_5->GetZaxis()->FindBin(qval);
    binentries = tranTemplate_5->GetBinContent(i + 1, energyBin, qbin);
    totentries = tranTemplate_5->Integral(i + 1, i + 1, energyBin, energyBin, 0, 100);
    if (qval > 0) {
      ++nbins;
      double prob = (double)binentries / totentries * 100;
      if (binentries > 0) tranLikelihood_5 += log(prob);
    }

  } // loop through

  double const tranLikelihood =
    (tranLikelihood_1 + tranLikelihood_2 + tranLikelihood_3 + tranLikelihood_4 + tranLikelihood_5) /
    nbins;
  std::cout << tranLikelihood << std::endl;
  return tranLikelihood;
} // getTranLikelihood

void shower::TCShowerElectronLikelihood::resetProfiles ( )

private

Definition at line 260 of file TCShowerElectronLikelihood_module.cc.

{

  longProfile->Reset();
  tranProfile->Reset();
  tranProfile_1->Reset();
  tranProfile_2->Reset();
  tranProfile_3->Reset();
  tranProfile_4->Reset();
  tranProfile_5->Reset();

  energyGuess = -9999;
  energyChi2 = -9999;
  maxt = -9999;
} // resetProfiles

Member Data Documentation

double shower::TCShowerElectronLikelihood::energyChi2

private

Definition at line 95 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::energyDist

private

Definition at line 75 of file TCShowerElectronLikelihood_module.cc.

int shower::TCShowerElectronLikelihood::energyGuess

private

Definition at line 94 of file TCShowerElectronLikelihood_module.cc.

calo::CalorimetryAlg shower::TCShowerElectronLikelihood::fCalorimetryAlg

private

Definition at line 114 of file TCShowerElectronLikelihood_module.cc.

std::string shower::TCShowerElectronLikelihood::fDigitModuleLabel

private

Definition at line 112 of file TCShowerElectronLikelihood_module.cc.

std::string shower::TCShowerElectronLikelihood::fGenieGenModuleLabel

private

Definition at line 111 of file TCShowerElectronLikelihood_module.cc.

std::string shower::TCShowerElectronLikelihood::fHitModuleLabel

private

Definition at line 108 of file TCShowerElectronLikelihood_module.cc.

std::string shower::TCShowerElectronLikelihood::fROOTfile

private

Definition at line 57 of file TCShowerElectronLikelihood_module.cc.

std::string shower::TCShowerElectronLikelihood::fShowerModuleLabel

private

Definition at line 109 of file TCShowerElectronLikelihood_module.cc.

std::string shower::TCShowerElectronLikelihood::fTemplateFile

private

Definition at line 56 of file TCShowerElectronLikelihood_module.cc.

std::string shower::TCShowerElectronLikelihood::fTemplateModuleLabel

private

Definition at line 110 of file TCShowerElectronLikelihood_module.cc.

const int shower::TCShowerElectronLikelihood::LBINS = 20

private

Definition at line 98 of file TCShowerElectronLikelihood_module.cc.

const int shower::TCShowerElectronLikelihood::LMAX = 5

private

Definition at line 100 of file TCShowerElectronLikelihood_module.cc.

const int shower::TCShowerElectronLikelihood::LMIN = 0

private

Definition at line 99 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::longLikelihoodHist

private

Definition at line 76 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::longProfHist

private

Definition at line 79 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::longProfile

private

Definition at line 86 of file TCShowerElectronLikelihood_module.cc.

TH3F* shower::TCShowerElectronLikelihood::longTemplate

private

Definition at line 59 of file TCShowerElectronLikelihood_module.cc.

TProfile2D* shower::TCShowerElectronLikelihood::longTemplateProf2D

private

Definition at line 66 of file TCShowerElectronLikelihood_module.cc.

int shower::TCShowerElectronLikelihood::maxt

private

Definition at line 96 of file TCShowerElectronLikelihood_module.cc.

const int shower::TCShowerElectronLikelihood::TBINS = 20

private

Definition at line 102 of file TCShowerElectronLikelihood_module.cc.

const int shower::TCShowerElectronLikelihood::TMAX = 5

private

Definition at line 104 of file TCShowerElectronLikelihood_module.cc.

const int shower::TCShowerElectronLikelihood::TMIN = -5

private

Definition at line 103 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranLikelihoodHist

private

Definition at line 77 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_1

private

Definition at line 80 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_2

private

Definition at line 81 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_3

private

Definition at line 82 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_4

private

Definition at line 83 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfHist_5

private

Definition at line 84 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfile

private

Definition at line 87 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfile_1

private

Definition at line 88 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfile_2

private

Definition at line 89 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfile_3

private

Definition at line 90 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfile_4

private

Definition at line 91 of file TCShowerElectronLikelihood_module.cc.

TH1F* shower::TCShowerElectronLikelihood::tranProfile_5

private

Definition at line 92 of file TCShowerElectronLikelihood_module.cc.

TH3F* shower::TCShowerElectronLikelihood::tranTemplate

private

Definition at line 60 of file TCShowerElectronLikelihood_module.cc.

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_1

private

Definition at line 61 of file TCShowerElectronLikelihood_module.cc.

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_2

private

Definition at line 62 of file TCShowerElectronLikelihood_module.cc.

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_3

private

Definition at line 63 of file TCShowerElectronLikelihood_module.cc.

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_4

private

Definition at line 64 of file TCShowerElectronLikelihood_module.cc.

TH3F* shower::TCShowerElectronLikelihood::tranTemplate_5

private

Definition at line 65 of file TCShowerElectronLikelihood_module.cc.

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D

private

Definition at line 67 of file TCShowerElectronLikelihood_module.cc.

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_1

private

Definition at line 68 of file TCShowerElectronLikelihood_module.cc.

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_2

private

Definition at line 69 of file TCShowerElectronLikelihood_module.cc.

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_3

private

Definition at line 70 of file TCShowerElectronLikelihood_module.cc.

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_4

private

Definition at line 71 of file TCShowerElectronLikelihood_module.cc.

TProfile2D* shower::TCShowerElectronLikelihood::tranTemplateProf2D_5

private

Definition at line 72 of file TCShowerElectronLikelihood_module.cc.

const double shower::TCShowerElectronLikelihood::X0 = 14

private

Definition at line 106 of file TCShowerElectronLikelihood_module.cc.

---

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

• larreco/larreco/ShowerFinder/TCShowerElectronLikelihood_module.cc