[u/mrichter/AliRoot.git] / TPC / AliTPCCalibRawBase.cxx

 /**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id: AliTPCCalibRawBase.cxx */

/////////////////////////////////////////////////////////////////////////////////////////
//                                                                                     //
//          Base class for the calibration algorithms using raw data as input          //
//                                                                                     //
//   Origin: Jens Wiechula   J.Wiechula@gsi.de                                         //
//                                                                                     //
/////////////////////////////////////////////////////////////////////////////////////////

//Root includes
#include <TDirectory.h>
#include <TFile.h>

//Aliroot includes
#include "AliRawReaderDate.h"
#include "AliRawReader.h"
#include "AliRawEventHeaderBase.h"
#include "AliAltroMapping.h"
#include "AliAltroRawStream.h"
#include "AliTPCROC.h"
#include "AliTPCRawStreamFast.h"
#include "AliTPCRawStreamV3.h"
#include "AliTPCRawStream.h"
#include "AliLog.h"
#include "TTreeStream.h"
#include "event.h"

#include "AliTPCCalibRawBase.h"

ClassImp(AliTPCCalibRawBase)

AliTPCCalibRawBase::AliTPCCalibRawBase() :
  TNamed(),
  fFirstTimeBin(0),
  fLastTimeBin(1000),
  fNevents(0),
  fDebugLevel(0),
  fStreamLevel(0),
  fRunNumber(0),
  fFirstTimeStamp(0),
  fLastTimeStamp(0),
  fTimeStamp(0),
  fEventType(0),
  fAltroL1Phase(0),
  fAltroL1PhaseTB(0),
  fCurrRCUId(-1),
  fPrevRCUId(-1),
  fCurrDDLNum(-1),
  fPrevDDLNum(-1),
  fUseL1Phase(kTRUE),
  fDebugStreamer(0x0),
  fAltroRawStream(0x0),
  fMapping(0x0),
  fROC(AliTPCROC::Instance())
{
    //
    // default ctor
    //

}
//_____________________________________________________________________
AliTPCCalibRawBase::AliTPCCalibRawBase(const AliTPCCalibRawBase &calib) :
  TNamed(calib),
  fFirstTimeBin(calib.fFirstTimeBin),
  fLastTimeBin(calib.fLastTimeBin),
  fNevents(calib.fNevents),
  fDebugLevel(calib.fDebugLevel),
  fStreamLevel(calib.fStreamLevel),
  fRunNumber(calib.fRunNumber),
  fFirstTimeStamp(calib.fFirstTimeStamp),
  fLastTimeStamp(calib.fLastTimeStamp),
  fTimeStamp(0),
  fEventType(0),
  fAltroL1Phase(0),
  fAltroL1PhaseTB(0),
  fCurrRCUId(-1),
  fPrevRCUId(-1),
  fCurrDDLNum(-1),
  fPrevDDLNum(-1),
  fUseL1Phase(kTRUE),
  fDebugStreamer(0x0),
  fAltroRawStream(0x0),
  fMapping(0x0),
  fROC(AliTPCROC::Instance())
{
    //
    // copy ctor
    //
  
}
//_____________________________________________________________________
AliTPCCalibRawBase::~AliTPCCalibRawBase()
{
  //
  // dtor
  //
  if (fDebugStreamer) delete fDebugStreamer;
}
//_____________________________________________________________________
  AliTPCCalibRawBase& AliTPCCalibRawBase::operator = (const  AliTPCCalibRawBase &source)
  {
    //
    // assignment operator
    //
    if (&source == this) return *this;
    new (this) AliTPCCalibRawBase(source);
    
    return *this;
  }
//_____________________________________________________________________
Bool_t AliTPCCalibRawBase::ProcessEventFast(AliTPCRawStreamFast * const rawStreamFast)
{
  //
  // Event Processing loop - AliTPCRawStreamFast
  //
  ResetEvent();
  Bool_t withInput = kFALSE;
  while ( rawStreamFast->NextDDL() ){
    while ( rawStreamFast->NextChannel() ){
      Int_t isector  = rawStreamFast->GetSector();                       //  current sector
      Int_t iRow     = rawStreamFast->GetRow();                          //  current row
      Int_t iPad     = rawStreamFast->GetPad();                          //  current pad
      
      while ( rawStreamFast->NextBunch() ){
        Int_t startTbin = (Int_t)rawStreamFast->GetStartTimeBin();
        Int_t endTbin = (Int_t)rawStreamFast->GetEndTimeBin();
        for (Int_t iTimeBin = startTbin; iTimeBin < endTbin; iTimeBin++){
          Float_t signal=(Float_t)rawStreamFast->GetSignals()[iTimeBin-startTbin];
            Update(isector,iRow,iPad,iTimeBin+1,signal);
          withInput = kTRUE;
        }
      }
    }
  }
  if (withInput){
    EndEvent();
  }
  return withInput;
}
//_____________________________________________________________________
Bool_t AliTPCCalibRawBase::ProcessEventFast(AliRawReader * const rawReader)
{
  //
  //  Event processing loop - AliRawReader
  //
  AliRawEventHeaderBase* eventHeader = (AliRawEventHeaderBase*)rawReader->GetEventHeader();
  if (eventHeader){
    fTimeStamp   = eventHeader->Get("Timestamp");
    fRunNumber = eventHeader->Get("RunNb");
    fEventType = eventHeader->Get("Type");
  }
  if (!fFirstTimeStamp) fFirstTimeStamp=fTimeStamp;
  fLastTimeStamp=fTimeStamp;
  
  AliTPCRawStreamFast *rawStreamFast = new AliTPCRawStreamFast(rawReader, (AliAltroMapping**)fMapping);
  Bool_t res=ProcessEventFast(rawStreamFast);
  delete rawStreamFast;
  return res;
}
//_____________________________________________________________________
Bool_t AliTPCCalibRawBase::ProcessEvent(AliTPCRawStreamV3 * const rawStreamV3)
{
  //
  // Event Processing loop - AliTPCRawStreamV3
  //
  ResetEvent();
  Bool_t withInput = kFALSE;
  fAltroL1Phase=0;
  fAltroL1PhaseTB=0;
//   fAltroRawStream = static_cast<AliAltroRawStream*>(rawStreamV3);
  while ( rawStreamV3->NextDDL() ){
    if (AliLog::GetGlobalDebugLevel()>2) rawStreamV3->PrintRCUTrailer();
    fPrevDDLNum=-1;
    fCurrRCUId=rawStreamV3->GetRCUId();
    fCurrDDLNum=rawStreamV3->GetDDLNumber();
    if (fUseL1Phase){
//         fAltroL1Phase  = fAltroRawStream->GetL1Phase();
      fAltroL1Phase  = rawStreamV3->GetL1Phase();
      fAltroL1PhaseTB = (fAltroL1Phase*1e09/100.);
      AliDebug(1, Form("L1Phase: %.2e (%03d)\n",fAltroL1PhaseTB,fCurrDDLNum));
    }
    UpdateDDL();
    while ( rawStreamV3->NextChannel() ){
      Int_t isector  = rawStreamV3->GetSector();                       //  current sector
      Int_t iRow     = rawStreamV3->GetRow();                          //  current row
      Int_t iPad     = rawStreamV3->GetPad();                          //  current pad
      while ( rawStreamV3->NextBunch() ){
        UInt_t  startTbin    = rawStreamV3->GetStartTimeBin();
//         Int_t  endTbin      = (Int_t)rawStreamV3->GetEndTimeBin();
        Int_t  bunchlength  = rawStreamV3->GetBunchLength();
        const UShort_t *sig = rawStreamV3->GetSignals();
        ProcessBunch(isector,iRow,iPad,bunchlength,startTbin,sig);
        for (Int_t iTimeBin = 0; iTimeBin<bunchlength; iTimeBin++){
          Float_t signal=(Float_t)sig[iTimeBin];
//            printf("%02d - %03d - %03d - %04d: %.1f\n",isector,iRow,iPad,startTbin,signal);
          Update(isector,iRow,iPad,startTbin--,signal);
          fPrevRCUId=fCurrRCUId;
          fPrevDDLNum=fCurrDDLNum;
          withInput = kTRUE;
        }
      }
    }
  }
  if (withInput){
    EndEvent();
  }
  return withInput;
}
//_____________________________________________________________________
Bool_t AliTPCCalibRawBase::ProcessEvent(AliRawReader * const rawReader)
{
  //
  //  Event processing loop - AliRawReader
  //
  AliRawEventHeaderBase* eventHeader = (AliRawEventHeaderBase*)rawReader->GetEventHeader();
  if (eventHeader){
    fTimeStamp   = eventHeader->Get("Timestamp");
    fRunNumber = eventHeader->Get("RunNb");
    fEventType = eventHeader->Get("Type");
  }
  if (!fFirstTimeStamp) fFirstTimeStamp=fTimeStamp;

  AliTPCRawStreamV3 *rawStreamV3 = new AliTPCRawStreamV3(rawReader, (AliAltroMapping**)fMapping);
  Bool_t res=ProcessEvent(rawStreamV3);
  
  fLastTimeStamp=fTimeStamp;
  
  delete rawStreamV3;
  return res;
}
//_____________________________________________________________________
Bool_t AliTPCCalibRawBase::ProcessEvent(AliTPCRawStream * const rawStream)
{
  //
  // Event Processing loop - AliTPCRawStream
  //

  ResetEvent();

  Bool_t withInput = kFALSE;
  fAltroL1Phase=0;
  fAltroL1PhaseTB=0;
  fAltroRawStream = static_cast<AliAltroRawStream*>(rawStream);
  while (rawStream->Next()) {
    if (fUseL1Phase){
      fAltroL1Phase  = fAltroRawStream->GetL1Phase();
      fAltroL1PhaseTB = (fAltroL1Phase*1e09/100.);
    }
    fPrevRCUId=fCurrRCUId;
    fCurrRCUId=rawStream->GetRCUId();
    fPrevDDLNum=fCurrDDLNum;
    fCurrDDLNum=rawStream->GetDDLNumber();
    Int_t isector  = rawStream->GetSector();                       //  current sector
    Int_t iRow     = rawStream->GetRow();                          //  current row
    Int_t iPad     = rawStream->GetPad();                          //  current pad
    Int_t iTimeBin = rawStream->GetTime();                         //  current time bin
    Float_t signal = rawStream->GetSignal();                       //  current ADC signal
    
    Update(isector,iRow,iPad,iTimeBin,signal);
    withInput = kTRUE;
  }
  fAltroRawStream=0x0;
  if (withInput){
    EndEvent();
  }
  return withInput;
}
//_____________________________________________________________________
Bool_t AliTPCCalibRawBase::ProcessEventOld(AliRawReader * const rawReader)
{
  //
  //  Event processing loop - AliRawReader
  //
  AliRawEventHeaderBase* eventHeader = (AliRawEventHeaderBase*)rawReader->GetEventHeader();
  if (eventHeader){
    fTimeStamp   = eventHeader->Get("Timestamp");
    fRunNumber = eventHeader->Get("RunNb");
    fEventType = eventHeader->Get("Type");
  }

  AliTPCRawStream rawStream(rawReader, (AliAltroMapping**)fMapping);
  rawReader->Select("TPC");
  return ProcessEvent(&rawStream);
}
//_____________________________________________________________________
Bool_t AliTPCCalibRawBase::ProcessEvent(eventHeaderStruct * const event)
{
  //
  //  Event processing loop - date event
  //

  fRunNumber=event->eventRunNb;
  fTimeStamp=event->eventTimestamp;
  if (!fFirstTimeStamp) fFirstTimeStamp=fTimeStamp;
  fLastTimeStamp=fTimeStamp;
  fEventType=event->eventType;
  AliRawReader *rawReader = new AliRawReaderDate((void*)event);
  AliTPCRawStreamV3 *rawStreamV3 = new AliTPCRawStreamV3(rawReader, (AliAltroMapping**)fMapping);
  Bool_t result=ProcessEvent(rawStreamV3);
  delete rawStreamV3;
  delete rawReader;
  return result;

}
//_____________________________________________________________________
void AliTPCCalibRawBase::DumpToFile(const Char_t *filename, const Char_t *dir, Bool_t append)
{
    //
    //  Write class to file
    //
  
  TString sDir(dir);
  TString option;
  
  if ( append )
    option = "update";
  else
    option = "recreate";
  
  TDirectory *backup = gDirectory;
  TFile f(filename,option.Data());
  f.cd();
  if ( !sDir.IsNull() ){
    f.mkdir(sDir.Data());
    f.cd(sDir);
  }
  this->Write();
  f.Close();
  
  if ( backup ) backup->cd();
}
//_____________________________________________________________________
TTreeSRedirector *AliTPCCalibRawBase::GetDebugStreamer(){
  //
  // Get Debug streamer
  // In case debug streamer not yet initialized and StreamLevel>0 create new one
  //
  if (fStreamLevel==0) return 0;
  if (fDebugStreamer) return fDebugStreamer;
  TString dsName;
  dsName=GetName();
  dsName+="Debug.root";
  dsName.ReplaceAll(" ","");
  fDebugStreamer = new TTreeSRedirector(dsName.Data());
  return fDebugStreamer;
}
//_____________________________________________________________________
void AliTPCCalibRawBase::MergeBase(const AliTPCCalibRawBase *calib)
{
  //
  // merge this with base
  //
  if (calib->fFirstTimeStamp<fFirstTimeStamp) fFirstTimeStamp=calib->fFirstTimeStamp;
  if (calib->fLastTimeStamp>fLastTimeStamp)   fLastTimeStamp =calib->fLastTimeStamp;
}
Commit	Line	Data
880c3382	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/* $Id: AliTPCCalibRawBase.cxx */
	17
	18	/////////////////////////////////////////////////////////////////////////////////////////
	19	// //
	20	// Base class for the calibration algorithms using raw data as input //
	21	// //
	22	// Origin: Jens Wiechula J.Wiechula@gsi.de //
	23	// //
	24	/////////////////////////////////////////////////////////////////////////////////////////
	25
	26	//Root includes
	27	#include <TDirectory.h>
	28	#include <TFile.h>
	29
	30	//Aliroot includes
	31	#include "AliRawReaderDate.h"
	32	#include "AliRawReader.h"
	33	#include "AliRawEventHeaderBase.h"
	34	#include "AliAltroMapping.h"
	35	#include "AliAltroRawStream.h"
	36	#include "AliTPCROC.h"
	37	#include "AliTPCRawStreamFast.h"
c3066940	38	#include "AliTPCRawStreamV3.h"
880c3382	39	#include "AliTPCRawStream.h"
c3066940	40	#include "AliLog.h"
880c3382	41	#include "TTreeStream.h"
	42	#include "event.h"
	43
	44	#include "AliTPCCalibRawBase.h"
	45
	46	ClassImp(AliTPCCalibRawBase)
	47
	48	AliTPCCalibRawBase::AliTPCCalibRawBase() :
	49	TNamed(),
	50	fFirstTimeBin(0),
	51	fLastTimeBin(1000),
	52	fNevents(0),
	53	fDebugLevel(0),
	54	fStreamLevel(0),
880c3382	55	fRunNumber(0),
78f17711	56	fFirstTimeStamp(0),
78f17711	57	fLastTimeStamp(0),
c3066940	58	fTimeStamp(0),
880c3382	59	fEventType(0),
	60	fAltroL1Phase(0),
	61	fAltroL1PhaseTB(0),
c3066940	62	fCurrRCUId(-1),
	63	fPrevRCUId(-1),
	64	fCurrDDLNum(-1),
	65	fPrevDDLNum(-1),
880c3382	66	fUseL1Phase(kTRUE),
	67	fDebugStreamer(0x0),
	68	fAltroRawStream(0x0),
	69	fMapping(0x0),
	70	fROC(AliTPCROC::Instance())
	71	{
	72	//
	73	// default ctor
	74	//
	75
	76	}
	77	//_____________________________________________________________________
	78	AliTPCCalibRawBase::AliTPCCalibRawBase(const AliTPCCalibRawBase &calib) :
	79	TNamed(calib),
	80	fFirstTimeBin(calib.fFirstTimeBin),
	81	fLastTimeBin(calib.fLastTimeBin),
	82	fNevents(calib.fNevents),
	83	fDebugLevel(calib.fDebugLevel),
	84	fStreamLevel(calib.fStreamLevel),
78f17711	85	fRunNumber(calib.fRunNumber),
	86	fFirstTimeStamp(calib.fFirstTimeStamp),
	87	fLastTimeStamp(calib.fLastTimeStamp),
c3066940	88	fTimeStamp(0),
880c3382	89	fEventType(0),
	90	fAltroL1Phase(0),
	91	fAltroL1PhaseTB(0),
c3066940	92	fCurrRCUId(-1),
	93	fPrevRCUId(-1),
	94	fCurrDDLNum(-1),
	95	fPrevDDLNum(-1),
880c3382	96	fUseL1Phase(kTRUE),
	97	fDebugStreamer(0x0),
	98	fAltroRawStream(0x0),
	99	fMapping(0x0),
	100	fROC(AliTPCROC::Instance())
	101	{
	102	//
	103	// copy ctor
	104	//
	105
	106	}
	107	//_____________________________________________________________________
	108	AliTPCCalibRawBase::~AliTPCCalibRawBase()
	109	{
	110	//
	111	// dtor
	112	//
	113	if (fDebugStreamer) delete fDebugStreamer;
	114	}
	115	//_____________________________________________________________________
	116	AliTPCCalibRawBase& AliTPCCalibRawBase::operator = (const AliTPCCalibRawBase &source)
	117	{
	118	//
	119	// assignment operator
	120	//
	121	if (&source == this) return *this;
	122	new (this) AliTPCCalibRawBase(source);
	123
	124	return *this;
	125	}
	126	//_____________________________________________________________________
7442bceb	127	Bool_t AliTPCCalibRawBase::ProcessEventFast(AliTPCRawStreamFast * const rawStreamFast)
880c3382	128	{
	129	//
	130	// Event Processing loop - AliTPCRawStreamFast
	131	//
	132	ResetEvent();
	133	Bool_t withInput = kFALSE;
	134	while ( rawStreamFast->NextDDL() ){
	135	while ( rawStreamFast->NextChannel() ){
	136	Int_t isector = rawStreamFast->GetSector(); // current sector
	137	Int_t iRow = rawStreamFast->GetRow(); // current row
	138	Int_t iPad = rawStreamFast->GetPad(); // current pad
	139
	140	while ( rawStreamFast->NextBunch() ){
	141	Int_t startTbin = (Int_t)rawStreamFast->GetStartTimeBin();
	142	Int_t endTbin = (Int_t)rawStreamFast->GetEndTimeBin();
	143	for (Int_t iTimeBin = startTbin; iTimeBin < endTbin; iTimeBin++){
	144	Float_t signal=(Float_t)rawStreamFast->GetSignals()[iTimeBin-startTbin];
c3066940	145	Update(isector,iRow,iPad,iTimeBin+1,signal);
880c3382	146	withInput = kTRUE;
	147	}
	148	}
	149	}
	150	}
	151	if (withInput){
	152	EndEvent();
	153	}
	154	return withInput;
	155	}
	156	//_____________________________________________________________________
7442bceb	157	Bool_t AliTPCCalibRawBase::ProcessEventFast(AliRawReader * const rawReader)
880c3382	158	{
	159	//
	160	// Event processing loop - AliRawReader
	161	//
	162	AliRawEventHeaderBase* eventHeader = (AliRawEventHeaderBase*)rawReader->GetEventHeader();
	163	if (eventHeader){
	164	fTimeStamp = eventHeader->Get("Timestamp");
	165	fRunNumber = eventHeader->Get("RunNb");
	166	fEventType = eventHeader->Get("Type");
	167	}
78f17711	168	if (!fFirstTimeStamp) fFirstTimeStamp=fTimeStamp;
	169	fLastTimeStamp=fTimeStamp;
	170
880c3382	171	AliTPCRawStreamFast rawStreamFast = new AliTPCRawStreamFast(rawReader, (AliAltroMapping*)fMapping);
	172	Bool_t res=ProcessEventFast(rawStreamFast);
	173	delete rawStreamFast;
	174	return res;
	175	}
	176	//_____________________________________________________________________
7442bceb	177	Bool_t AliTPCCalibRawBase::ProcessEvent(AliTPCRawStreamV3 * const rawStreamV3)
c3066940	178	{
	179	//
	180	// Event Processing loop - AliTPCRawStreamV3
	181	//
	182	ResetEvent();
	183	Bool_t withInput = kFALSE;
	184	fAltroL1Phase=0;
	185	fAltroL1PhaseTB=0;
	186	// fAltroRawStream = static_cast<AliAltroRawStream*>(rawStreamV3);
	187	while ( rawStreamV3->NextDDL() ){
	188	if (AliLog::GetGlobalDebugLevel()>2) rawStreamV3->PrintRCUTrailer();
5312f439	189	fPrevDDLNum=-1;
	190	fCurrRCUId=rawStreamV3->GetRCUId();
	191	fCurrDDLNum=rawStreamV3->GetDDLNumber();
c3066940	192	if (fUseL1Phase){
	193	// fAltroL1Phase = fAltroRawStream->GetL1Phase();
	194	fAltroL1Phase = rawStreamV3->GetL1Phase();
	195	fAltroL1PhaseTB = (fAltroL1Phase*1e09/100.);
5312f439	196	AliDebug(1, Form("L1Phase: %.2e (%03d)\n",fAltroL1PhaseTB,fCurrDDLNum));
c3066940	197	}
5312f439	198	UpdateDDL();
c3066940	199	while ( rawStreamV3->NextChannel() ){
	200	Int_t isector = rawStreamV3->GetSector(); // current sector
	201	Int_t iRow = rawStreamV3->GetRow(); // current row
	202	Int_t iPad = rawStreamV3->GetPad(); // current pad
	203	while ( rawStreamV3->NextBunch() ){
96bf9029	204	UInt_t startTbin = rawStreamV3->GetStartTimeBin();
c3066940	205	// Int_t endTbin = (Int_t)rawStreamV3->GetEndTimeBin();
96bf9029	206	Int_t bunchlength = rawStreamV3->GetBunchLength();
c3066940	207	const UShort_t *sig = rawStreamV3->GetSignals();
96bf9029	208	ProcessBunch(isector,iRow,iPad,bunchlength,startTbin,sig);
c3066940	209	for (Int_t iTimeBin = 0; iTimeBin<bunchlength; iTimeBin++){
	210	Float_t signal=(Float_t)sig[iTimeBin];
	211	// printf("%02d - %03d - %03d - %04d: %.1f\n",isector,iRow,iPad,startTbin,signal);
	212	Update(isector,iRow,iPad,startTbin--,signal);
	213	fPrevRCUId=fCurrRCUId;
	214	fPrevDDLNum=fCurrDDLNum;
	215	withInput = kTRUE;
	216	}
	217	}
	218	}
	219	}
	220	if (withInput){
	221	EndEvent();
	222	}
	223	return withInput;
	224	}
	225	//_____________________________________________________________________
7442bceb	226	Bool_t AliTPCCalibRawBase::ProcessEvent(AliRawReader * const rawReader)
c3066940	227	{
	228	//
	229	// Event processing loop - AliRawReader
	230	//
	231	AliRawEventHeaderBase* eventHeader = (AliRawEventHeaderBase*)rawReader->GetEventHeader();
	232	if (eventHeader){
	233	fTimeStamp = eventHeader->Get("Timestamp");
	234	fRunNumber = eventHeader->Get("RunNb");
	235	fEventType = eventHeader->Get("Type");
	236	}
78f17711	237	if (!fFirstTimeStamp) fFirstTimeStamp=fTimeStamp;
78f17711	238
c3066940	239	AliTPCRawStreamV3 rawStreamV3 = new AliTPCRawStreamV3(rawReader, (AliAltroMapping*)fMapping);
c3066940	240	Bool_t res=ProcessEvent(rawStreamV3);
78f17711	241
	242	fLastTimeStamp=fTimeStamp;
	243
c3066940	244	delete rawStreamV3;
	245	return res;
	246	}
	247	//_____________________________________________________________________
7442bceb	248	Bool_t AliTPCCalibRawBase::ProcessEvent(AliTPCRawStream * const rawStream)
880c3382	249	{
	250	//
	251	// Event Processing loop - AliTPCRawStream
	252	//
	253
	254	ResetEvent();
	255
	256	Bool_t withInput = kFALSE;
	257	fAltroL1Phase=0;
	258	fAltroL1PhaseTB=0;
	259	fAltroRawStream = static_cast<AliAltroRawStream*>(rawStream);
	260	while (rawStream->Next()) {
	261	if (fUseL1Phase){
	262	fAltroL1Phase = fAltroRawStream->GetL1Phase();
	263	fAltroL1PhaseTB = (fAltroL1Phase*1e09/100.);
	264	}
c3066940	265	fPrevRCUId=fCurrRCUId;
	266	fCurrRCUId=rawStream->GetRCUId();
	267	fPrevDDLNum=fCurrDDLNum;
	268	fCurrDDLNum=rawStream->GetDDLNumber();
880c3382	269	Int_t isector = rawStream->GetSector(); // current sector
	270	Int_t iRow = rawStream->GetRow(); // current row
	271	Int_t iPad = rawStream->GetPad(); // current pad
	272	Int_t iTimeBin = rawStream->GetTime(); // current time bin
	273	Float_t signal = rawStream->GetSignal(); // current ADC signal
	274
	275	Update(isector,iRow,iPad,iTimeBin,signal);
	276	withInput = kTRUE;
	277	}
	278	fAltroRawStream=0x0;
	279	if (withInput){
	280	EndEvent();
	281	}
	282	return withInput;
	283	}
	284	//_____________________________________________________________________
7442bceb	285	Bool_t AliTPCCalibRawBase::ProcessEventOld(AliRawReader * const rawReader)
880c3382	286	{
	287	//
	288	// Event processing loop - AliRawReader
	289	//
	290	AliRawEventHeaderBase* eventHeader = (AliRawEventHeaderBase*)rawReader->GetEventHeader();
	291	if (eventHeader){
	292	fTimeStamp = eventHeader->Get("Timestamp");
	293	fRunNumber = eventHeader->Get("RunNb");
	294	fEventType = eventHeader->Get("Type");
	295	}
	296
	297	AliTPCRawStream rawStream(rawReader, (AliAltroMapping**)fMapping);
	298	rawReader->Select("TPC");
	299	return ProcessEvent(&rawStream);
	300	}
	301	//_____________________________________________________________________
7442bceb	302	Bool_t AliTPCCalibRawBase::ProcessEvent(eventHeaderStruct * const event)
880c3382	303	{
	304	//
	305	// Event processing loop - date event
	306	//
5312f439	307
	308	fRunNumber=event->eventRunNb;
	309	fTimeStamp=event->eventTimestamp;
78f17711	310	if (!fFirstTimeStamp) fFirstTimeStamp=fTimeStamp;
78f17711	311	fLastTimeStamp=fTimeStamp;
5312f439	312	fEventType=event->eventType;
	313	AliRawReader rawReader = new AliRawReaderDate((void)event);
	314	AliTPCRawStreamV3 rawStreamV3 = new AliTPCRawStreamV3(rawReader, (AliAltroMapping*)fMapping);
	315	Bool_t result=ProcessEvent(rawStreamV3);
	316	delete rawStreamV3;
	317	delete rawReader;
	318	return result;
880c3382	319
	320	}
	321	//_____________________________________________________________________
	322	void AliTPCCalibRawBase::DumpToFile(const Char_t filename, const Char_t dir, Bool_t append)
	323	{
	324	//
	325	// Write class to file
	326	//
	327
	328	TString sDir(dir);
	329	TString option;
	330
	331	if ( append )
	332	option = "update";
	333	else
	334	option = "recreate";
	335
	336	TDirectory *backup = gDirectory;
	337	TFile f(filename,option.Data());
	338	f.cd();
	339	if ( !sDir.IsNull() ){
	340	f.mkdir(sDir.Data());
	341	f.cd(sDir);
	342	}
	343	this->Write();
	344	f.Close();
	345
	346	if ( backup ) backup->cd();
	347	}
	348	//_____________________________________________________________________
	349	TTreeSRedirector *AliTPCCalibRawBase::GetDebugStreamer(){
	350	//
	351	// Get Debug streamer
	352	// In case debug streamer not yet initialized and StreamLevel>0 create new one
	353	//
	354	if (fStreamLevel==0) return 0;
	355	if (fDebugStreamer) return fDebugStreamer;
	356	TString dsName;
	357	dsName=GetName();
	358	dsName+="Debug.root";
	359	dsName.ReplaceAll(" ","");
	360	fDebugStreamer = new TTreeSRedirector(dsName.Data());
	361	return fDebugStreamer;
	362	}
78f17711	363	//_____________________________________________________________________
	364	void AliTPCCalibRawBase::MergeBase(const AliTPCCalibRawBase *calib)
	365	{
	366	//
	367	// merge this with base
	368	//
	369	if (calib->fFirstTimeStamp<fFirstTimeStamp) fFirstTimeStamp=calib->fFirstTimeStamp;
	370	if (calib->fLastTimeStamp>fLastTimeStamp) fLastTimeStamp =calib->fLastTimeStamp;
	371	}
	372