TV1730RawWaveform.cxx

Go to the documentation of this file.

#include "TV1730RawWaveform.h"
 
#include "TV1730RawData.hxx"
#include "TDirectory.h"          
 
/// Reset the histogram for this canvas
TV1730RawWaveform::TV1730RawWaveform(){
 
  SetSubTabName("V1730 Waveforms");
  SetUpdateOnlyWhenPlotted(true);
  SetNumSamples(64);     
  SetNanosecsPerSample(2); //ADC clock runs at 500Mhz on the v1730 = units of 2 nsecs
 
  CreateHistograms();
}
 
 
void TV1730RawWaveform::CreateHistograms(){
 
    std::cout << "Creating V1730 raw waveforms... " << std::endl;
 
  // check if we already have histograms.
  char tname[100];
  sprintf(tname,"V1730Raw_%i_%i",0,0);
 
  TH1D *tmp = (TH1D*)gDirectory->Get(tname);
  if (tmp) return;
 
  // Otherwise make histogram
  clear();
 
    for(int i = 0; i < 16; i++){ // loop over 16 channels
        
        char name[100];
        char title[100];
        sprintf(name,"V1730Raw_%i",i);
        
        sprintf(title,"V1730 Waveform for channel=%i",i);   
        
        TH1D *tmp = new TH1D(name, title, this->numSamples, 0, this->numSamples);
        tmp->GetXaxis()->SetLimits(0, this->numSamples*this->nanosecsPerSample);
        tmp->SetXTitle("ns");
        tmp->SetYTitle("ADC value");
        
        push_back(tmp);
    }
 
    std::cout << "Finish V1730 raw waveforms... " << std::endl;
 
}
 
 
 
 
void TV1730RawWaveform::UpdateHistograms(TDataContainer& dataContainer){
 
  int eventid = dataContainer.GetMidasData().GetEventId();
  int timestamp = dataContainer.GetMidasData().GetTimeStamp();
    
    TV1730RawData *v1730 = dataContainer.GetEventData<TV1730RawData>("V730");
 
    if(v1730 ){      
 
        std::vector<RawChannelMeasurement> measurements = v1730->GetMeasurements();
 
        for(int i = 0; i < measurements.size(); i++){
            
            int chan = measurements[i].GetChannel();
            
        // Reset the histogram...
            for(int ib = 0; ib < this->numSamples; ib++)
                GetHistogram(chan)->SetBinContent(ib+1,0);
 
            // Hack!
            float offset = 0;
            if(chan == 1)
                offset = 35;
            //std::cout << "Nsamples " <<  measurements[i].GetNSamples() << std::endl;
            for(int ib = 0; ib < measurements[i].GetNSamples(); ib++){
 
 
                GetHistogram(chan)->SetBinContent(ib+1, measurements[i].GetSample(ib)-offset);
                //std::cout << "Setting " << chan << " " << ib << " " << measurements[i].GetSample(ib) << std::endl;
 
            }
 
        }
 
    }
 
}
 
 
void TV1730RawWaveform::Reset(){
 
    // Loop over all the data, reset histos
  for(int iBoard=0; iBoard<32; iBoard++){// Loop over V1730Dpp boards        
        for(int i = 0; i < 8; i++){ // loop over channels           
            int index = iBoard*8 + i;
            
            // Reset the histogram...
            for(int ib = 0; ib < 2500; ib++)
                GetHistogram(index)->SetBinContent(ib,0);
            GetHistogram(index)->Reset();
        }
    }
}
        
/// Take actions at begin run
void TV1730RawWaveform::BeginRun(int transition,int run,int time){
 
  CreateHistograms();
 
}
 
/// Take actions at end run  
void TV1730RawWaveform::EndRun(int transition,int run,int time){
 
}