plots.cxx File Reference

#include <HistMan/HistMan.h>
#include <BeamData/ana/bmnt/ProfMon.h>
#include <BeamData/ana/bmnt/Knot.h>
#include <BeamData/ana/bmnt/Pedestals.h>
#include <BeamData/ana/bmnt/AcnetDevice.h>
#include <TCanvas.h>
#include <TSystem.h>
#include <TTree.h>
#include <TFile.h>
#include <TH1D.h>
#include <TH2D.h>
#include <iostream>
#include <string>
#include <cassert>
#include <cmath>

Go to the source code of this file.

Functions
void	book_once ()
void	fill (Knot &k, ProfMon &, Pedestals &, int)
void	fill_profile (ProfMon &pm, int count, TH1D &xh, TH1D &yh, TH2D &xa, TH2D &ya)
void	book_and_fill (Knot &k, ProfMon &pm, Pedestals &peds, int count)
void	write_plots ()
void	make_plots ()
void	make_prints ()
void	plots ()
Variables
const double	MIN_ETOR101 = 0.15
const char *	toroids []
const char *	profile_monitor []
int	nentries = 0

---

Function Documentation

void book_and_fill (  Knot &  k, ProfMon &  pm, Pedestals &  peds, int  count )

Definition at line 91 of file plots.cxx. References HistMan::Book(), count, fill_profile(), Form(), Knot::GetDevice(), Pedestals::GetPeds(), max, min, profile_monitor, Swic::SetAcnetDevice(), Swic::SetPeds(), and ProfMon::WirePosition(). Referenced by make_plots(). { HistMan hm(Form("bv/spill%03d",count)); double min = ProfMon::WirePosition(ProfMon::MIN_CHANNEL); double max = ProfMon::WirePosition(ProfMon::MAX_CHANNEL); min -= ProfMon::WIRE_SPACING; max -= ProfMon::WIRE_SPACING; int nbins = ProfMon::MAX_CHANNEL - ProfMon::MIN_CHANNEL + 1; int npms = 0; for (; profile_monitor[npms]; ++npms); TH2D* xall = hm.Book<TH2D>(Form("xallpm_%03d",count), Form("All X Profiles for spill %d",count), npms,0,npms,nbins,min,max); TH2D* yall = hm.Book<TH2D>(Form("yallpm_%03d",count), Form("All Y Profiles for spill %d",count), npms,0,npms,nbins,min,max); for (int ind=0; profile_monitor[ind]; ++ind) { const AcnetDevice* ad = k.GetDevice(profile_monitor[ind]); if (!pm.SetPeds(peds.GetPeds(profile_monitor[ind]))) { cerr << "Failed to set peds for " << profile_monitor[ind] << " skipping\n"; continue; } int ret = pm.SetAcnetDevice(*ad); if (! ret) { cerr << "Failed to set acnet device for spill " << count << endl; continue; } string pms = profile_monitor[ind]; pms[1]='_'; TH1D* xhist = hm.Book<TH1D>(Form("%sX_%03d",pms.c_str(),count), Form("ProfMon %s X view spill %d",pms.c_str(),count), nbins,min,max); TH1D* yhist = hm.Book<TH1D>(Form("%sY_%03d",pms.c_str(),count), Form("ProfMon %s Y view spill %d",pms.c_str(),count), nbins,min,max); fill_profile(pm,ind,*xhist,*yhist,*xall,*yall); ++npms; } }

void book_once (   )

Definition at line 53 of file plots.cxx. References HistMan::Book(). Referenced by make_plots(). { HistMan hm("bv"); // ... hm.Book<TH1D>("max_dae_dt","Maximum difference between 2 DAEs in data block",1200,-60,60); }

void fill (  Knot &  k, ProfMon &  , Pedestals &  , int  )

Definition at line 61 of file plots.cxx. References HistMan::Fill1d(), Knot::GetDevice(), profile_monitor, and AcnetDevice::timestamp. Referenced by main(), and make_plots(). { double max_dt = 0; for (int ind1=0; profile_monitor[ind1]; ++ind1) { const AcnetDevice* ad1 = k.GetDevice(profile_monitor[ind1]); for (int ind2=ind1; profile_monitor[ind2]; ++ind2) { const AcnetDevice* ad2 = k.GetDevice(profile_monitor[ind2]); double diff = fabs(ad1->timestamp - ad2->timestamp); if (diff > max_dt) max_dt = diff; } } HistMan hm("bv"); hm.Fill1d("max_dae_dt",max_dt); }

void fill_profile (  ProfMon &  pm, int  count, TH1D &  xh, TH1D &  yh, TH2D &  xa, TH2D &  ya )

Definition at line 78 of file plots.cxx. References count, Swic::GetVoltage(), ProfMon::MAX_CHANNEL, ProfMon::MIN_CHANNEL, ProfMon::WirePosition(), ProfMon::Xindex(), and ProfMon::Yindex(). Referenced by book_and_fill(). { for (int ch=pm.MIN_CHANNEL; ch<=pm.MAX_CHANNEL; ++ch) { xh.Fill(pm.WirePosition(ch),pm.GetVoltage(pm.Xindex(ch))); xa.Fill(count,pm.WirePosition(ch),pm.GetVoltage(pm.Xindex(ch))); yh.Fill(pm.WirePosition(ch),pm.GetVoltage(pm.Yindex(ch))); ya.Fill(count,pm.WirePosition(ch),pm.GetVoltage(pm.Yindex(ch))); } }

void make_plots (   )

Definition at line 152 of file plots.cxx. References book_and_fill(), book_once(), count, AcnetDevice::data, fill(), Knot::GetDevice(), Knot::GetEntry(), Knot::GetSize(), Pedestals::LoadPedFile(), MIN_ETOR101, AcnetDevice::ndata, nentries, and Swic::SetScale(). Referenced by plots(). { //TFile f("first-data.tree.root"); TFile f("second.bmnt.root"); TTree* bd = dynamic_cast<TTree*>(f.FindObjectAny("bd")); assert(bd); Knot k(*bd); book_once(); Pedestals peds; peds.LoadPedFile("peds.txt"); nentries = k.GetSize(); int count = 0; for (int entry=0; entry<nentries; ++entry) { cerr << entry << " "; k.GetEntry(entry); const AcnetDevice* ad = k.GetDevice("E:TOR101"); if (ad->ndata == 0) { cerr << "No E:TOR101!\n"; continue; } double cut = fabs(ad->data[0]); if (cut < MIN_ETOR101) continue; ++count; ProfMon pm; pm.SetScale(cut/MIN_ETOR101); book_and_fill(k,pm,peds,count); fill(k,pm,peds,count); } //write(); }

void make_prints (   )

Definition at line 188 of file plots.cxx. References Form(), and HistMan::Get(). Referenced by plots(). { TCanvas c("c","BeamData Plots",850,1100); c.Print("bdplots.ps["); HistMan hm("bv"); hm.Get<TH1D>("max_dae_dt")->Draw(); c.Print("bdplots.ps"); #if 0 for (int ind=1; ind <= nentries; ++ind) { HistMan hm(Form("bv/spill%03d",ind)); TH2D* hx = hm.Get<TH2D>(Form("xallpm_%03d",ind)); TH2D* hy = hm.Get<TH2D>(Form("yallpm_%03d",ind)); if (!hx || !hy) { cerr << "Failed to get hists for spill " << ind << endl; continue; } c.Clear(); c.Divide(1,2); c.cd(1); hx->Draw("colz"); c.cd(2); hy->Draw("colz"); c.Print("bdplots.ps"); } #endif c.Print("bdplots.ps]"); }

void plots (   )

Definition at line 223 of file plots.cxx. References make_plots(), make_prints(), and write_plots(). Referenced by ChopModule::Ana(), NuCutImps::Bravo::Bravo(), NuCutImps::BravoPrime::BravoPrime(), NuDSTAna::BRevAna(), NuCutImps::CC0325Std::CC0325Std(), NuCutImps::CC0720Test::CC0720Test(), NuCutImps::CCA::CCA(), NuCutImps::CCA_NC::CCA_NC(), NuCutImps::CCAPresel::CCAPresel(), NuCutImps::ChairSound::ChairSound(), NuAnalysis::ChargeSeparationOneSnarl(), NuAnalysis::ChargeSignCut(), NuCutImps::Charlie::Charlie(), NuDSTAna::CoilHoleAna(), NuDSTAna::CSSAna(), NuDSTAna::CSSAnaRashid(), NuCutImps::Delta::Delta(), NuAnalysis::Efficiencies(), NuAnalysis::EnergySpect(), NuDSTAna::FDTestAna(), NuPlots::FillCoilHole(), NuCutter::GetNewCutInstance(), NuDSTAna::IsGoodPreSelectionCuts(), NuDSTAna::IsGoodStdCuts(), NuDSTAna::JeffsTestAna(), NuAnalysis::LIRejectionTest(), NuAnalysis::LoopOverTruthInfo(), NuDSTAna::MakeFCTree(), NuDSTAna::MakeFinalPlots(), NuAnalysis::MakeFullDST(), NuDSTAna::MakeMicroDST(), NuDSTAna::MakeMicroDstFakeData(), NuDSTAna::MakeMicroDstForCSSSystematics(), NuDSTAna::MakeMicroDstHe(), NuDSTAna::MakeMicroDstJJEPresel(), NuDSTAna::MakeMicroDstWithStdCCRecoAndCuts(), NuDSTAna::MakePostPreSelectionPlots(), NuDSTAna::MakeSelMicroDST(), NuDSTAna::MMRereco(), NuDSTAna::MMTransition(), NuAnalysis::N_1(), NuDSTAna::NDTestAna(), NuDSTAna::NewFieldAna(), NuDSTAna::NMBAna(), NuAnalysis::NMBSummaryTreeAna(), NuCut::NuCut(), NuCutImps::NuCutsSelection::NuCutsSelection(), NuCutter::NuCutter(), NuAnalysis::NuMuBarAppearance(), NuDSTAna::QPStudy(), NuDSTAna::SelectorTable(), NuCutter::SetCut(), NuDSTAna::StdCCAna(), NuDSTAna::StdNMBAna(), and NuDSTAna::VsTime(). { make_plots(); make_prints(); write_plots(); }

void write_plots (   )

Definition at line 142 of file plots.cxx. References HistMan::WriteOut(). Referenced by plots(). { HistMan hm("bv"); TFile file("plots.root","RECREATE"); hm.WriteOut(file); file.Close(); }

---

Variable Documentation

const double MIN_ETOR101 = 0.15

Definition at line 29 of file plots.cxx. Referenced by make_plots().

int nentries = 0 [static]

Definition at line 150 of file plots.cxx. Referenced by make_plots().

const char* profile_monitor[]

Initial value: { "E:M101DS", "E:M105DS", "E:M107DS", "E:M108DS", "E:M112DS", "E:M114DS", "E:M115DS", "E:M117DS", "E:M121DS", "E:MTGTDS", 0 } Definition at line 39 of file plots.cxx.

const char* toroids[]

Initial value: { "E:TOR101", "E:TORTGT", 0 } Definition at line 33 of file plots.cxx. Referenced by BDCheckDB::BDCheckDB(), BDCheckDB::CheckSpill(), BeamMonDbuModule::FillPipeline(), BDSwicPeds::IsPedSpill(), and BDDevices::Toroids().

---