AltTimeHist.cxx

Go to the documentation of this file.

#include "AltTimeHist.h"
#include <Midad/Base/Mint.h>
#include <Midad/Base/MintView.h>
#include <Midad/Base/ZoomPad.h>
#include <Midad/Util/Range.h>

#include <Conventions/Munits.h>

#include <CandDigit/CandDigitHandle.h>
#include <CandDigit/CandDigitListHandle.h>

#include <DataUtil/CDL2STL.h>
using namespace DataUtil;

#include <TH1F.h>

#include <sigc++/sigc++.h>
#include <sigc++/class_slot.h>
using namespace SigC;

#include <vector>
using namespace std;

AltTimeHist::AltTimeHist()
    : fPad(0)
    , fHist1(0)
    , fHist2(0)
{

  fPad = new TPad("Alt Time Hist2", "Alt Time Hist2", 0.42,0.25,0.75,0.44);
  fPad->Draw();
  
  fPad->cd();
  fHist1 = new TH1F("AltTimeHist","Digit Times (ns)",100,0,100);
  fHist2 = new TH1F("AltTimeHistwin","Digit Times (ns)",100,0,100);
  fHist1->Draw();
  fHist2->Draw("same");
  fHist1->SetFillColor(2);
  fHist1->SetLineColor(2);
  fHist1->SetFillStyle(1001);
  fHist2->SetFillColor(3);
  fHist2->SetLineColor(3);
  fHist2->SetFillStyle(1001);
 
}


AltTimeHist::~AltTimeHist()
{
}

void AltTimeHist::Draw(Option_t*)
{
    fPad->Draw();
}
void AltTimeHist::SetMint(Mint* mint)
{
    this->MintView::SetMint(mint);
    mint->GetTimeRange().modified.connect(slot_class(*this,&AltTimeHist::Refresh));
        
}

void AltTimeHist::SetTotTimeRange(RangeDouble& r){
  fTotTimeRange.SetMin(r.Min()); 
  fTotTimeRange.SetMax(r.Max());
}
 
void AltTimeHist::UpdatePad()
{
    fPad->Modified();
    fPad->Update();
}

void AltTimeHist::Refresh()
{
//    cerr << "AltTimeHist::Refresh()\n";

    fHist1->Reset();

    fHist2->Reset();

    const CandDigitListHandle *cdlh = this->GetMint().GetDigits();
    if (! cdlh) return;    

    RangeDouble& r = this->GetMint().GetTimeRange();

    if(fTotTimeRange.Max()-fTotTimeRange.Min()>0){
      fHist1->SetBins(300,fTotTimeRange.Min()/Munits::ns, fTotTimeRange.Max()/Munits::ns);    
      fHist2->SetBins(300,fTotTimeRange.Min()/Munits::ns, fTotTimeRange.Max()/Munits::ns);
      
      vector<const CandDigitHandle*> dc = CDL2STLvector<const CandDigitHandle>(*cdlh);
      
      int siz = dc.size();
      for (int ind=0; ind < siz; ++ind) {
        fHist1->Fill(dc[ind]->GetTime()/Munits::ns,dc[ind]->GetCharge());
      }
      
      for (int ind=0; ind < siz; ++ind) {
        Float_t t = dc[ind]->GetTime();
        if(t>r.Min() && t<r.Max()){
          fHist2->Fill(dc[ind]->GetTime()/Munits::ns,dc[ind]->GetCharge());
        }
      }
    }

    this->UpdatePad();
}

---