#include "AliRawReaderRoot.h"
#include "AliRawReaderDate.h"
#include "AliTPCRawStream.h"
-#include "AliTPCRawStreamFast.h"
#include "AliTPCCalROC.h"
#include "AliTPCCalPad.h"
#include "AliTPCROC.h"
if ( !force || arr->UncheckedAt(sector) )
return (TH2S*)arr->UncheckedAt(sector);
- // if we are forced and histogram doesn't yes exist create it
- Char_t name[255], title[255];
-
- sprintf(name,"hCalib%s%.2d",type,sector);
- sprintf(title,"%s calibration histogram sector %.2d",type,sector);
-
- // new histogram with Q calib information. One value for each pad!
- TH2S* hist = new TH2S(name,title,
+ // if we are forced and histogram doesn't yes exist create it
+ // new histogram with Q calib information. One value for each pad!
+ TH2S* hist = new TH2S(Form("hCalib%s%.2d",type,sector),Form("%s calibration histogram sector %.2d",type,sector),
nbinsY, ymin, ymax,
fROC->GetNChannels(sector),0,fROC->GetNChannels(sector));
hist->SetDirectory(0);
//_____________________________________________________________________
void AliTPCCalibPulser::Merge(AliTPCCalibPulser * const sig)
{
- //
- // Merge reference histograms of sig to the current AliTPCCalibPulser
- //
-
- //merge histograms
+ //
+ // Merge reference histograms of sig to the current AliTPCCalibPulser
+ //
+
+ MergeBase(sig);
+ //merge histograms
for (Int_t iSec=0; iSec<72; ++iSec){
TH2S *hRefQmerge = sig->GetHistoQ(iSec);
TH2S *hRefT0merge = sig->GetHistoT0(iSec);
rocT0->SetValue(iChannel, cogTime0+meanTsec); //offset by mean time of the sector
rocRMS->SetValue(iChannel, cogRMS);
rocOut->SetValue(iChannel, cogOut);
+
+ // in case a channel has no data set the value to 0
+ if (TMath::Abs(cogTime0-fXminT0)<1e-10){
+ rocQ->SetValue(iChannel, 0);
+ rocT0->SetValue(iChannel, 0); //offset by mean time of the sector
+ rocRMS->SetValue(iChannel, 0);
+ }
+
//debug
if ( GetStreamLevel() > 2 ){
TTreeSRedirector *streamer=GetDebugStreamer();