nlana::SPLifetime Class Reference

Inheritance diagram for nlana::SPLifetime:

Public Member Functions
	SPLifetime (fhicl::ParameterSet const &p)
	SPLifetime (SPLifetime const &)=delete
	SPLifetime (SPLifetime &&)=delete
SPLifetime &	operator= (SPLifetime const &)=delete
SPLifetime &	operator= (SPLifetime &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
void	endJob () override
void	reconfigure (fhicl::ParameterSet const &p)
void	respondToOpenInputFile (art::FileBlock const &infileblock) override
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 Attributes
std::string	fClusterModuleLabel
double	fChiCut
std::vector< float >	fChgCuts
double	fTicksPerBin
unsigned	fMinBins
unsigned	fMinHits
double	fMinDTickSNR
double	fMinDWireSNR
unsigned	fMinHitsSNR
int	lastRun
int	fDebugCluster
std::string	fInFilename
bool	fIsRealData
double	bigLifeInv [12]
double	bigLifeInvErr [12]
double	bigLifeInvCnt [12]
double	signalToNoise [12]
double	signalToNoiseCnt [12]
unsigned int	signalToNoiseClsCnt [12]
TH1F *	fLife
TH2F *	fLifeVTPC
TH1F *	fLifeInv
TH1F *	fFracSelHits
TH1F *	fChiDOF
TH1F *	fLifeInvCA
TH1F *	fLifeInvAC
TProfile *	fLifeInv_E
TProfile *	fLifeInv_Angle
TH1F *	fDriftTime
TH2F *	fDriftTimeVTPC
TH1F *	fSNR
TH2F *	fSNRVTPC
TH1F *	fAmplitudes
TH1F *	fNoise
TH1F *	fNoiseWide

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 67 of file SPLifetime_module.cc.

Constructor & Destructor Documentation

nlana::SPLifetime::SPLifetime ( fhicl::ParameterSet const &  p )

explicit

Definition at line 134 of file SPLifetime_module.cc.

  :
  EDAnalyzer(pset),
  fInFilename("NoInFilenameFound"),
  fIsRealData(true)

 // More initializers here.
{
  reconfigure(pset);
}

nlana::SPLifetime::SPLifetime ( SPLifetime const &  )

delete

nlana::SPLifetime::SPLifetime ( SPLifetime &&  )

delete

Member Function Documentation

void nlana::SPLifetime::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 449 of file SPLifetime_module.cc.

{
  
//  int event  = evt.id().event(); 
  int run    = evt.run();
//  int subrun = evt.subRun();
  fIsRealData = evt.isRealData();
  
  if(lastRun < 0) lastRun = run;
  
  if(run != lastRun) mf::LogVerbatim("LIFE")<<"The run number has changed from "<<lastRun<<" to "<<run<<" This might not be good...";
  
  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
  double msPerTick = 1E-6 * sampling_rate(clockData);


  art::ValidHandle<std::vector<recob::Cluster>> clsVecHandle = evt.getValidHandle<std::vector<recob::Cluster>>(fClusterModuleLabel);
  art::FindManyP<recob::Hit> clsHitsFind(clsVecHandle, evt, fClusterModuleLabel);
  art::ValidHandle<std::vector<recob::Wire>> wireVecHandle = evt.getValidHandle<std::vector<recob::Wire>>("caldata");
  
//  bool prt = false;
  
  
  for(unsigned int icl = 0; icl < clsVecHandle->size(); ++icl) {
    art::Ptr<recob::Cluster> cls = art::Ptr<recob::Cluster>(clsVecHandle, icl);
    // only consider the collection plane
    if(cls->Plane().Plane != 2) continue;
//    prt = (fDebugCluster >= 0 && icl == (unsigned int)fDebugCluster);
    float sTick = cls->StartTick();
    float eTick = cls->EndTick();
    if(sTick > eTick) std::swap(sTick, eTick);
    float dTick = eTick - sTick;
    // Get the hits
    std::vector<art::Ptr<recob::Hit> > clsHits;
    clsHitsFind.get(icl, clsHits);
    if(clsHits.size() == 0) continue;
    unsigned short tpc = clsHits[0]->WireID().TPC;
    fDriftTime->Fill(dTick*msPerTick);
    fDriftTimeVTPC->Fill(tpcMapping.at(tpc),dTick*msPerTick);
    unsigned short nhist = 1 + (unsigned short)(dTick / fTicksPerBin);
    if(nhist < fMinBins) continue;
    if(clsHits.size() < fMinHits) continue;
/*
    if(prt) {
      auto& sht = clsHits[0];
      std::cout<<"Cls "<<icl<<" "<<sht->WireID().TPC<<":"<<sht->WireID().Plane<<":"<<sht->WireID().Wire<<":"<<(int)sht->PeakTime();
      auto& eht = clsHits[clsHits.size()-1];
      std::cout<<"  "<<eht->WireID().TPC<<":"<<eht->WireID().Plane<<":"<<eht->WireID().Wire<<":"<<(int)eht->PeakTime();
      std::cout<<" sTick "<<(int)sTick<<" "<<(int)eTick<<" nhits "<<clsHits.size()<<" nhist "<<nhist<<"\n";
    }
*/
    // Find the average and maximum charge of these histograms
    std::vector<double> tck(nhist), ave(nhist), cnt(nhist), err(nhist);
    std::vector<float> minChg(nhist, 0);
    std::vector<float> maxChg(nhist, 10000);
    for(unsigned short nit = 0; nit < 2; ++nit) {
      for(unsigned short ihist = 0; ihist < nhist; ++ihist) {
        tck[ihist] = 0;
        ave[ihist] = 0;
        cnt[ihist] = 0;
        err[ihist] = 0;
      } // ihist
      // Sum to get the (truncated) average
      for(auto& pht : clsHits) {
        unsigned short ihist = (pht->PeakTime() - sTick) / fTicksPerBin;
        if(ihist > nhist - 1) continue;
        float chg = pht->Integral();
        if(chg < minChg[ihist] || chg > maxChg[ihist]) continue;
        tck[ihist] += pht->PeakTime() - sTick;
        ave[ihist] += chg;
        err[ihist] += chg * chg;
        ++cnt[ihist];
      } // ii
      for(unsigned short ihist = 0; ihist < nhist; ++ihist) if(cnt[ihist] > 0) {
        tck[ihist] /= cnt[ihist];
        ave[ihist] /= cnt[ihist];
      }
      if(nit == 0) {
        // Calculate a max charge cut on the first iteration
        for(unsigned short ihist = 0; ihist < nhist; ++ihist) {
          maxChg[ihist] = 0;
          if(cnt[ihist] < 5) continue;
          maxChg[ihist] = fChgCuts[1] * ave[ihist];
          minChg[ihist] = fChgCuts[0] * ave[ihist];
//          if(prt) std::cout<<ihist<<" Min "<<(int)minChg[ihist]<<" max "<<(int)maxChg[ihist]<<"\n";
        } // ihist
      } else {
        // Calculate the error on the average on the second iteration
        for(unsigned short ihist = 0; ihist < nhist; ++ihist) {
          if(cnt[ihist] < 3) continue;
          double arg = err[ihist] - cnt[ihist] * ave[ihist] * ave[ihist];
          if(arg < 0) {
            cnt[ihist] = 0;
            continue;
          }
          // calculate rms
          err[ihist] = sqrt(arg / (cnt[ihist] - 1));
          // convert to error on the mean
          err[ihist] /= sqrt(cnt[ihist]);
        } // ihist
      } // second iteration
    } // nit
    
    // require a minimum count in each histogram
    unsigned short nok = 0;
    for(auto& icnt : cnt) if(icnt > 4) ++nok;
    if(nok < 5) {
//      if(prt) std::cout<<" failed nok cut "<<nok<<" Need at least 4 \n";
      continue;
    }
/*
    if(prt) {
      for(unsigned short ihist = 0; ihist < nhist; ++ihist) {
        float xx = (ihist + 0.5) * fTicksPerBin + sTick;
        std::cout<<"ihist "<<ihist<<" tick "<<(int)xx<<" ave "<<(int)ave[ihist]<<" +/- "<<(int)err[ihist]<<" cnt "<<(int)cnt[ihist]<<"\n";
      } // ihist
    }
*/
    // fit to find the 1/lifetime
    double sum = 0.;
    double sumx = 0.;
    double sumy = 0.;
    double sumxy = 0.;
    double sumx2 = 0.;
    double sumy2 = 0.;
    double xx, yy, wght, arg;
    double fitcnt = 0;
    
    for(unsigned short ihist = 0; ihist < nhist; ++ihist) {
      if(cnt[ihist] < 3) continue;
      if(err[ihist] == 0) continue;
      xx = (tck[ihist] - tck[0]) * msPerTick;
      yy = log(ave[ihist]);
      // error on log(x) = dx / x
      arg = ave[ihist] / err[ihist];
      wght = arg * arg;
      sum += wght;
      sumx += wght * xx;
      sumy += wght * yy;
      sumx2 += wght * xx * xx;
      sumxy += wght * xx * yy;
      sumy2 += wght * yy * yy;
      ++fitcnt;
    } // ihist
    // calculate coefficients
    double delta = sum * sumx2 - sumx * sumx;
    if(delta == 0.) continue;
    double A = (sumx2 * sumy - sumx * sumxy) / delta;
    // slope = 1 / lifetime
    double B = (sumxy * sum  - sumx * sumy) / delta;
    if(B > 0) continue;
    // calculate the error
    double ndof = fitcnt - 2;
    double varnce = (sumy2 + A*A*sum + B*B*sumx2 - 2 * (A*sumy + B*sumxy - A*B*sumx)) / ndof;
    if(varnce == 0) continue;
//    double BErr = sqrt(varnce * sum / delta);
//    if(prt) std::cout<<"B "<<B<<" "<<BErr;
    
    double lifeInv = -B;
//    double lifeInvErr = BErr / (B * B) ;
    
    // calculate chisq
    double chi = 0;
    for(unsigned short ihist = 0; ihist < nhist; ++ihist) {
      if(cnt[ihist] < 3) continue;
      if(err[ihist] == 0) continue;
      xx = (tck[ihist] - tck[0]) * msPerTick;
      yy = exp(A - xx * lifeInv);
      arg = (yy - ave[ihist]) / err[ihist];
      chi += arg * arg;
//      if(prt) std::cout<<"chk "<<ihist<<" xx "<<xx<<" yy "<<yy<<" ave "<<ave[ihist]<<" arg "<<arg<<"\n";
    }
    chi /= ndof;
//    if(prt) std::cout<<tpc<<" lifeInv "<<lifeInv<<" +/- "<<lifeInvErr<<" chi "<<chi<<"\n";
    fChiDOF->Fill(chi);
    if(chi > fChiCut) continue;
    ++bigLifeInvCnt[tpc];
    bigLifeInv[tpc] += lifeInv;
    bigLifeInvErr[tpc] += lifeInv * lifeInv;
    fLifeInv->Fill(lifeInv);
    fLife->Fill(1./lifeInv);
    fLifeVTPC->Fill(tpcMapping.at(tpc),1./lifeInv);
    
    double selHits = 0;
    for(auto& hcnt : cnt) selHits += hcnt;
    selHits /= (double)(clsHits.size());
    fFracSelHits->Fill(selHits);

    fLifeInv_Angle->Fill(cls->StartAngle(), lifeInv);
  } // icl
  
  std::vector<std::pair<unsigned int, float>> chanPedRMS;
//  std::cout << "n wire data: "<<wireVecHandle->size() << std::endl;
  for(unsigned int witer = 0; witer < wireVecHandle->size(); ++witer) {
    art::Ptr<recob::Wire> thisWire(wireVecHandle, witer);
    const recob::Wire::RegionsOfInterest_t& signalROI = thisWire->SignalROI();
    for(const auto& range : signalROI.get_ranges()) {
      const std::vector<float>& signal = range.data();
//      std::cout << "channel: "<<thisWire->Channel()<<" range begin: " << range.begin_index() << " range size: "<<range.size()<<"signal length: " << signal.size() << std::endl;
      if(signal.size() != 1) continue;
      // protect against unrealistic values
      float rms = signal[0];
      if(rms < 0.001) rms = 0.001;
      chanPedRMS.push_back(std::make_pair(thisWire->Channel(), rms));
//      std::cout << "channel: "<<thisWire->Channel()<<" rms: " << rms << std::endl;
    }
  }// witer

  // calculate S/N using shallow angle clusters
  for(unsigned int icl = 0; icl < clsVecHandle->size(); ++icl) {
    art::Ptr<recob::Cluster> cls = art::Ptr<recob::Cluster>(clsVecHandle, icl);
    // only consider the collection plane
    if(cls->Plane().Plane != 2) continue;
    float dTick = std::abs(cls->StartTick() - cls->EndTick());
    if(dTick < fMinDTickSNR) continue;
    float dWire = std::abs(cls->StartWire() - cls->EndWire());
    if(dWire < fMinDWireSNR) continue;
    // Get the hits
    std::vector<art::Ptr<recob::Hit> > clsHits;
    clsHitsFind.get(icl, clsHits);
    if(clsHits.size() < fMinHitsSNR) continue;
    unsigned short tpc = clsHits[0]->WireID().TPC;
    ++signalToNoiseClsCnt[tpc];
//    float aveSN = 0;
//    float cnt = 0;
    for(auto& hit : clsHits) {
      float pedRMS = -1;
      for(auto& chrms : chanPedRMS) {
        if(chrms.first != hit->Channel()) continue;
        pedRMS = chrms.second;
        break;
      } // chrms
      if(pedRMS < 0) continue;
      float snr = hit->PeakAmplitude() / pedRMS;
      //std::cout<<"SNR "<<(int)snr << " S: " << hit->PeakAmplitude() << " N: "<< pedRMS <<"\n";
      signalToNoise[tpc] += snr;
      ++signalToNoiseCnt[tpc];
      fAmplitudes->Fill(hit->PeakAmplitude());
      fNoise->Fill(pedRMS);
      fNoiseWide->Fill(pedRMS);
      fSNR->Fill(snr);
      fSNRVTPC->Fill(tpcMapping.at(tpc),snr);
//      aveSN += sn;
//      ++cnt;
    } // hit
/*
    if(cnt == 0) continue;
    aveSN /= cnt;
    std::cout<<" aveSN "<<aveSN<<" cnt "<<(int)cnt<<"\n";
*/
  } // icl

} // analyze

void nlana::SPLifetime::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 146 of file SPLifetime_module.cc.

{
  // Implementation of optional member function here.
  
  art::ServiceHandle<art::TFileService> tfs;
  fLife = tfs->make<TH1F>("Life","Electron Lifetime for all APAs", 50, 0, 15);
  setHistTitles(fLife,"Cluster Electron Lifetime [ms]", "Clusters / Bin");
  fLifeVTPC = tfs->make<TH2F>("LifeVTPC","Electron Lifetime v. APA", 6,0.5,6.5,50, 0, 15.);
  setHistTitles(fLifeVTPC,"","Cluster Electron Lifetime [ms]");
  fLifeInv = tfs->make<TH1F>("LifeInv","LifeInv", 50, 0, 1);
  fFracSelHits = tfs->make<TH1F>("FracSelHits","FracSelHits", 50, 0, 1);
  fChiDOF = tfs->make<TH1F>("ChiDOF","ChiDOF", 80, 0, 80);
  fLifeInvCA = tfs->make<TH1F>("LifeInvCA","LifeInv C to A", 50, 0, 1);
  fLifeInvAC = tfs->make<TH1F>("LifeInvAC","LifeInv A to C", 50, 0, 1);
  
  fLifeInv_E = tfs->make<TProfile>("LifeInv_E","LifeInv_E", 20, 0, 40);
  fLifeInv_Angle = tfs->make<TProfile>("LifeInv_Angle","LifeInv_Angle", 15, -1.5, 1.5);

  // initialize an invalid run number
  lastRun = -1;
  for(unsigned short tpc = 0; tpc < 12; ++tpc) {
    bigLifeInv[tpc] = 0;
    bigLifeInvErr[tpc] = 0;
    bigLifeInvCnt[tpc] = 0;
    signalToNoise[tpc] = 0;
    signalToNoiseCnt[tpc] = 0;
    signalToNoiseClsCnt[tpc] = 0;
  }

  fDriftTime = tfs->make<TH1F>("DriftTime","Drift Time", 500, 0.,10.);
  setHistTitles(fDriftTime,"Cluster Drift Time [ms]", "Clusters / Bin");
  fDriftTimeVTPC = tfs->make<TH2F>("DriftTimeVTPC","Drift Time v. APA", 6,0.5,6.5,500,0.,10.);
  setHistTitles(fDriftTimeVTPC,"","Cluster Drift Time [ms]");

  fSNR = tfs->make<TH1F>("SNR","Signal to Noise Ratio", 400, 0.,200);
  setHistTitles(fSNR,"Signal to Noise Ratio", "Hits / Bin");
  fSNRVTPC = tfs->make<TH2F>("SNRVTPC","Signal to Noise Ratio v. APA", 6,0.5,6.5,400,0.,200.);
  setHistTitles(fSNRVTPC,"TPC Number","Signal to Noise Ratio");

  fAmplitudes = tfs->make<TH1F>("Amplitudes","Hit Amplitude", 4100, 0.,4100.);
  setHistTitles(fAmplitudes,"Hit Amplitude [ADC]", "Hits / Bin");
  fNoise = tfs->make<TH1F>("Noise","Wire Noise < 10 ADC", 2000, 0.,10.);
  setHistTitles(fNoise,"RMS Noise [ADC]", "Hit Wires / Bin");
  fNoiseWide = tfs->make<TH1F>("NoiseWide","Wire Noise", 2000, 0.,1000.);
  setHistTitles(fNoiseWide,"RMS Noise [ADC]", "Hit Wires / Bin");

  for(size_t apa = 0; apa < 6; ++apa){
    fLifeVTPC->GetXaxis()->SetBinLabel(apa+1,apaLabels.at(apa).c_str());
    fDriftTimeVTPC->GetXaxis()->SetBinLabel(apa+1,apaLabels.at(apa).c_str());
    fSNRVTPC->GetXaxis()->SetBinLabel(apa+1,apaLabels.at(apa).c_str());
  }

} // beginJob

void nlana::SPLifetime::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 217 of file SPLifetime_module.cc.

{
  std::ofstream purfile;
  purfile.open("Lifetime_Run" + std::to_string(lastRun) + ".txt");
  //purfile<<"Run, tpc, lifetime, error, count, S/N, num S/N clusters\n";
  purfile<<"Run, tpc, lifetime, error, count, S/N, num S/N clusters, drift time\n";

  for(unsigned short tpc = 0; tpc < 12; ++tpc) {
    if(bigLifeInvCnt[tpc] < 2) continue;
    if(bigLifeInv[tpc] <= 0) continue;
    bigLifeInv[tpc] /= bigLifeInvCnt[tpc];
    bigLifeInvErr[tpc] = bigLifeInvErr[tpc] - bigLifeInvCnt[tpc] * bigLifeInv[tpc] * bigLifeInv[tpc];
    if(bigLifeInvErr[tpc] < 0) continue;
    bigLifeInvErr[tpc] = sqrt(bigLifeInvErr[tpc] / (bigLifeInvCnt[tpc] - 1));
    // convert to error on the mean
    bigLifeInvErr[tpc] /= sqrt(bigLifeInvCnt[tpc]);
    float life = 0;
    if(1/bigLifeInv[tpc] != 0) life = 1 / bigLifeInv[tpc];
    float lifeErr = life * bigLifeInvErr[tpc] / bigLifeInv[tpc];
    float sn = -1;
    if(signalToNoiseCnt[tpc] > 100) sn = signalToNoise[tpc] / signalToNoiseCnt[tpc];
    purfile<<lastRun<<", "<<tpc<<", "<<std::fixed<<std::setprecision(2)<<life<<", "<<lifeErr<<", "<<(int)bigLifeInvCnt[tpc];
    purfile<<", "<<std::setprecision(1)<<sn<<", "<<signalToNoiseClsCnt[tpc];

    // Now do drift time
    size_t apa = tpcMapping.at(tpc);
    if (apa == 0) continue;
    size_t nBinsY = fDriftTimeVTPC->GetNbinsY();
    float driftTime = -1.;
    for (int iBin=nBinsY; iBin >=0; iBin--)
    {
      if (fDriftTimeVTPC->GetBinContent(apa,iBin) > 1)
      {
        driftTime = fDriftTimeVTPC->GetYaxis()->GetBinCenter(iBin);
        break;
      }
    }
    purfile<<", "<<std::setprecision(2)<<driftTime;

    purfile<<"\n";
  }
  purfile.close();

  // Now for images
  TCanvas * canvas = new TCanvas("canvas_SPLifetime");
  canvas->SetLogz();
  canvas->SetBottomMargin(0.13);
  canvas->SetRightMargin(0.12);
  std::string imageFileName;

  // Get filename in format p3s wants
  // "run": "run003907_0001_dl05",
  unsigned runNum = lastRun;
  unsigned fileNum = 1;
  unsigned dlNum = 5;
  if (fIsRealData)
  {
    static const std::string patternString = "run(\\d+)_(\\d+)_dl(\\d+)";
    static const std::regex pattern(patternString);
    std::smatch patternMatches;
    const bool foundMatch = std::regex_search(fInFilename,patternMatches,pattern);
    if(foundMatch)
    {
      std::cout << "filename: " << fInFilename << " match: " << patternMatches[0] << " run: " << patternMatches[1] << " file: " << patternMatches[2] << " dl: " << patternMatches[3] << std::endl;
      runNum = std::stoi(patternMatches[1]);
      fileNum = std::stoi(patternMatches[2]);
      dlNum = std::stoi(patternMatches[3]);
    }
    else
    {
      //throw cet::exception("InputFilenameError")
      //  <<"Couldn't parse input filename: '"<<fInFilename<<"' with regex '"<<patternString<<"'";
      std::cout <<"Error: couldn't parse input filename: '"<<fInFilename<<"' with regex '"<<patternString<<"'\n";
    }
  }

  std::stringstream infileNameStrStream;
  infileNameStrStream << "run";
  infileNameStrStream << std::setfill('0') << std::setw(6) << runNum;
  infileNameStrStream << std::setw(0) << '_';
  infileNameStrStream << std::setw(4) << fileNum;
  const std::string identifierNoDL = infileNameStrStream.str();
  infileNameStrStream << std::setw(0) << "_dl";
  infileNameStrStream << std::setw(2) << dlNum;
  const std::string identifierDL = infileNameStrStream.str();

  std::cout << "identifierDL: " << identifierDL << std::endl;
  std::cout << "identifierNoDL: " << identifierNoDL << std::endl;
  const std::string identifierTitle = " " + identifierDL;

  //std::string identifierBase = infilenameStripped;
  //identifierBase.erase(identifierBase.rfind("_dl"));
  //std::cout << "identifierBase: '"<<identifierBase<<"'";
  //std::string dlStr = infilenameStripped;
  //dlStr.erase(0,dlStr.rfind("_dl")+3);
  //std::cout << "dlStr: '"<<dlStr<<"'";

  //// summary json file
  // fn should be like run003907_0001_purity_summary.json
  std::string summary_filename = identifierNoDL+"_purity_summary.json";
  std::string filelist_filename = identifierNoDL+"_purity_FileList.json";

  std::ofstream summaryfile;
  summaryfile.open(summary_filename);
  summaryfile << "[\n  {\n    \"run\": \"" << identifierDL << "\",\n"
              << "    \"Type\": \"purity\"\n  }\n]\n";
  summaryfile.close();

  // file list json file
  std::ofstream filelistfile;
  filelistfile.open(filelist_filename);
  filelistfile << "[\n  {\n    \"Category\": \"Purity Monitor\",\n    \"Files\": {\n      \"Cluster Drift Time\": \"";

  imageFileName = "driftVTPC_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  appendToHistTitle(fDriftTimeVTPC,identifierTitle);
  fDriftTimeVTPC->SetStats(false);
  fDriftTimeVTPC->GetXaxis()->SetLabelSize(0.050);
  fDriftTimeVTPC->Draw("colz");
  canvas->SaveAs(imageFileName.c_str());
  filelistfile << imageFileName<<",";

  imageFileName = "driftVTPC_zoom_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  std::string originalTitle = fDriftTimeVTPC->GetTitle();
  fDriftTimeVTPC->SetTitle((originalTitle+" From 0 to 4 ms").c_str());
  fDriftTimeVTPC->GetYaxis()->SetRangeUser(0,4);
  fDriftTimeVTPC->Draw("colz");
  canvas->SaveAs(imageFileName.c_str());
  fDriftTimeVTPC->SetTitle(originalTitle.c_str());
  filelistfile << imageFileName;

  filelistfile << "\",\n";

  // Now e lifetime
  filelistfile << "      \"Purity\": \"";
  imageFileName = "purity_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  appendToHistTitle(fLife,identifierTitle);
  fLife->Draw();
  canvas->SaveAs(imageFileName.c_str());
  filelistfile << imageFileName<<",";

  canvas->SetLogz(false);
  imageFileName = "purityVTPC_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  appendToHistTitle(fLifeVTPC,identifierTitle);
  fLifeVTPC->SetStats(false);
  fLifeVTPC->GetXaxis()->SetLabelSize(0.050);
  fLifeVTPC->Draw("colz");
  canvas->SaveAs(imageFileName.c_str());
  filelistfile << imageFileName;
  canvas->SetLogz();

  filelistfile << "\",\n";

  // Now SNR
  filelistfile << "      \"SNR\": \"";
  imageFileName = "snr_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  canvas->SetLogy(true);
  appendToHistTitle(fSNR,identifierTitle);
  fLife->Draw();
  fSNR->Draw();
  canvas->SaveAs(imageFileName.c_str());
  canvas->SetLogy(false);
  filelistfile << imageFileName<<",";

  imageFileName = "snrVTPC_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  appendToHistTitle(fSNRVTPC,identifierTitle);
  fSNRVTPC->SetStats(false);
  fSNRVTPC->GetXaxis()->SetLabelSize(0.050);
  fSNRVTPC->Draw("colz");
  canvas->SaveAs(imageFileName.c_str());
  filelistfile << imageFileName<<",";

  canvas->SetLogy(true);
  imageFileName = "noise_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  appendToHistTitle(fNoise,identifierTitle);
  fNoise->Draw();
  canvas->SaveAs(imageFileName.c_str());
  filelistfile << imageFileName<<",";

  imageFileName = "noiseWide_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  appendToHistTitle(fNoiseWide,identifierTitle);
  fNoiseWide->Draw();
  canvas->SaveAs(imageFileName.c_str());
  filelistfile << imageFileName<<",";

  imageFileName = "amplitude_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  appendToHistTitle(fAmplitudes,identifierTitle);
  fAmplitudes->Draw();
  canvas->SaveAs(imageFileName.c_str());
  filelistfile << imageFileName<<",";

  imageFileName = "amplitude_zoom_";
  imageFileName += identifierNoDL;
  imageFileName += ".png";
  originalTitle = fAmplitudes->GetTitle();
  fAmplitudes->SetTitle((originalTitle+" From 0 to 500 ADC").c_str());
  fAmplitudes->GetXaxis()->SetRangeUser(0,500);
  fAmplitudes->Draw();
  canvas->SaveAs(imageFileName.c_str());
  filelistfile << imageFileName;
  canvas->SetLogy(false);

  filelistfile << "\"\n    }\n  }\n]\n";
  filelistfile.close();
  delete canvas;

} // endJob

SPLifetime& nlana::SPLifetime::operator= ( SPLifetime const &  )

delete

SPLifetime& nlana::SPLifetime::operator= ( SPLifetime &&  )

delete

void nlana::SPLifetime::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 201 of file SPLifetime_module.cc.

{
  // Implementation of optional member function here.
  fClusterModuleLabel  = pset.get<std::string>("ClusterModuleLabel");
  fChgCuts             = pset.get<std::vector<float>>("ChgCuts", {0, 5});
  fChiCut              = pset.get<double>("ChiCut", 3);
  fTicksPerBin         = pset.get<double>("TicksPerBin");
  fMinBins             = pset.get<unsigned>("MinBins");
  fMinHits             = pset.get<unsigned>("MinHits");
  fMinDTickSNR         = pset.get<double>("MinDTickSNR");
  fMinDWireSNR         = pset.get<double>("MinDWireSNR");
  fMinHitsSNR          = pset.get<unsigned>("MinHitsSNR");
  fDebugCluster        = pset.get<int>("DebugCluster");
} // reconfigure

void nlana::SPLifetime::respondToOpenInputFile ( art::FileBlock const &  infileblock )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 443 of file SPLifetime_module.cc.

{
    fInFilename = infileblock.fileName();
} // respondToOpenInputFile

Member Data Documentation

double nlana::SPLifetime::bigLifeInv[12]

private

Definition at line 103 of file SPLifetime_module.cc.

double nlana::SPLifetime::bigLifeInvCnt[12]

private

Definition at line 105 of file SPLifetime_module.cc.

double nlana::SPLifetime::bigLifeInvErr[12]

private

Definition at line 104 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fAmplitudes

private

Definition at line 127 of file SPLifetime_module.cc.

std::vector<float> nlana::SPLifetime::fChgCuts

private

Definition at line 92 of file SPLifetime_module.cc.

double nlana::SPLifetime::fChiCut

private

Definition at line 91 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fChiDOF

private

Definition at line 114 of file SPLifetime_module.cc.

std::string nlana::SPLifetime::fClusterModuleLabel

private

Definition at line 90 of file SPLifetime_module.cc.

int nlana::SPLifetime::fDebugCluster

private

Definition at line 100 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fDriftTime

private

Definition at line 121 of file SPLifetime_module.cc.

TH2F* nlana::SPLifetime::fDriftTimeVTPC

private

Definition at line 122 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fFracSelHits

private

Definition at line 113 of file SPLifetime_module.cc.

std::string nlana::SPLifetime::fInFilename

private

Definition at line 101 of file SPLifetime_module.cc.

bool nlana::SPLifetime::fIsRealData

private

Definition at line 102 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fLife

private

Definition at line 110 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fLifeInv

private

Definition at line 112 of file SPLifetime_module.cc.

TProfile* nlana::SPLifetime::fLifeInv_Angle

private

Definition at line 119 of file SPLifetime_module.cc.

TProfile* nlana::SPLifetime::fLifeInv_E

private

Definition at line 118 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fLifeInvAC

private

Definition at line 116 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fLifeInvCA

private

Definition at line 115 of file SPLifetime_module.cc.

TH2F* nlana::SPLifetime::fLifeVTPC

private

Definition at line 111 of file SPLifetime_module.cc.

unsigned nlana::SPLifetime::fMinBins

private

Definition at line 94 of file SPLifetime_module.cc.

double nlana::SPLifetime::fMinDTickSNR

private

Definition at line 96 of file SPLifetime_module.cc.

double nlana::SPLifetime::fMinDWireSNR

private

Definition at line 97 of file SPLifetime_module.cc.

unsigned nlana::SPLifetime::fMinHits

private

Definition at line 95 of file SPLifetime_module.cc.

unsigned nlana::SPLifetime::fMinHitsSNR

private

Definition at line 98 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fNoise

private

Definition at line 128 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fNoiseWide

private

Definition at line 129 of file SPLifetime_module.cc.

TH1F* nlana::SPLifetime::fSNR

private

Definition at line 124 of file SPLifetime_module.cc.

TH2F* nlana::SPLifetime::fSNRVTPC

private

Definition at line 125 of file SPLifetime_module.cc.

double nlana::SPLifetime::fTicksPerBin

private

Definition at line 93 of file SPLifetime_module.cc.

int nlana::SPLifetime::lastRun

private

Definition at line 99 of file SPLifetime_module.cc.

double nlana::SPLifetime::signalToNoise[12]

private

Definition at line 106 of file SPLifetime_module.cc.

unsigned int nlana::SPLifetime::signalToNoiseClsCnt[12]

private

Definition at line 108 of file SPLifetime_module.cc.

double nlana::SPLifetime::signalToNoiseCnt[12]

private

Definition at line 107 of file SPLifetime_module.cc.

---

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

• duneprototypes/duneprototypes/Protodune/singlephase/NearlineMonitor/SPLifetime_module.cc