[u/mrichter/AliRoot.git] / TPC / TPCbase / 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 "AliTPCRawStreamV3.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::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(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"
c3066940	37	#include "AliTPCRawStreamV3.h"
c3066940	38	#include "AliLog.h"
880c3382	39	#include "TTreeStream.h"
	40	#include "event.h"
	41
	42	#include "AliTPCCalibRawBase.h"
	43
	44	ClassImp(AliTPCCalibRawBase)
	45
	46	AliTPCCalibRawBase::AliTPCCalibRawBase() :
	47	TNamed(),
	48	fFirstTimeBin(0),
	49	fLastTimeBin(1000),
	50	fNevents(0),
	51	fDebugLevel(0),
	52	fStreamLevel(0),
880c3382	53	fRunNumber(0),
78f17711	54	fFirstTimeStamp(0),
78f17711	55	fLastTimeStamp(0),
c3066940	56	fTimeStamp(0),
880c3382	57	fEventType(0),
	58	fAltroL1Phase(0),
	59	fAltroL1PhaseTB(0),
c3066940	60	fCurrRCUId(-1),
	61	fPrevRCUId(-1),
	62	fCurrDDLNum(-1),
	63	fPrevDDLNum(-1),
880c3382	64	fUseL1Phase(kTRUE),
	65	fDebugStreamer(0x0),
	66	fAltroRawStream(0x0),
	67	fMapping(0x0),
	68	fROC(AliTPCROC::Instance())
	69	{
	70	//
	71	// default ctor
	72	//
	73
	74	}
	75	//_____________________________________________________________________
	76	AliTPCCalibRawBase::AliTPCCalibRawBase(const AliTPCCalibRawBase &calib) :
	77	TNamed(calib),
	78	fFirstTimeBin(calib.fFirstTimeBin),
	79	fLastTimeBin(calib.fLastTimeBin),
	80	fNevents(calib.fNevents),
	81	fDebugLevel(calib.fDebugLevel),
	82	fStreamLevel(calib.fStreamLevel),
78f17711	83	fRunNumber(calib.fRunNumber),
	84	fFirstTimeStamp(calib.fFirstTimeStamp),
	85	fLastTimeStamp(calib.fLastTimeStamp),
c3066940	86	fTimeStamp(0),
880c3382	87	fEventType(0),
	88	fAltroL1Phase(0),
	89	fAltroL1PhaseTB(0),
c3066940	90	fCurrRCUId(-1),
	91	fPrevRCUId(-1),
	92	fCurrDDLNum(-1),
	93	fPrevDDLNum(-1),
880c3382	94	fUseL1Phase(kTRUE),
	95	fDebugStreamer(0x0),
	96	fAltroRawStream(0x0),
	97	fMapping(0x0),
	98	fROC(AliTPCROC::Instance())
	99	{
	100	//
	101	// copy ctor
	102	//
	103
	104	}
	105	//_____________________________________________________________________
	106	AliTPCCalibRawBase::~AliTPCCalibRawBase()
	107	{
	108	//
	109	// dtor
	110	//
	111	if (fDebugStreamer) delete fDebugStreamer;
	112	}
	113	//_____________________________________________________________________
	114	AliTPCCalibRawBase& AliTPCCalibRawBase::operator = (const AliTPCCalibRawBase &source)
	115	{
	116	//
	117	// assignment operator
	118	//
	119	if (&source == this) return *this;
	120	new (this) AliTPCCalibRawBase(source);
	121
	122	return *this;
	123	}
	124	//_____________________________________________________________________
7442bceb	125	Bool_t AliTPCCalibRawBase::ProcessEvent(AliTPCRawStreamV3 * const rawStreamV3)
c3066940	126	{
	127	//
	128	// Event Processing loop - AliTPCRawStreamV3
	129	//
	130	ResetEvent();
	131	Bool_t withInput = kFALSE;
	132	fAltroL1Phase=0;
	133	fAltroL1PhaseTB=0;
	134	// fAltroRawStream = static_cast<AliAltroRawStream*>(rawStreamV3);
	135	while ( rawStreamV3->NextDDL() ){
	136	if (AliLog::GetGlobalDebugLevel()>2) rawStreamV3->PrintRCUTrailer();
5312f439	137	fPrevDDLNum=-1;
	138	fCurrRCUId=rawStreamV3->GetRCUId();
	139	fCurrDDLNum=rawStreamV3->GetDDLNumber();
c3066940	140	if (fUseL1Phase){
	141	// fAltroL1Phase = fAltroRawStream->GetL1Phase();
	142	fAltroL1Phase = rawStreamV3->GetL1Phase();
	143	fAltroL1PhaseTB = (fAltroL1Phase*1e09/100.);
5312f439	144	AliDebug(1, Form("L1Phase: %.2e (%03d)\n",fAltroL1PhaseTB,fCurrDDLNum));
c3066940	145	}
5312f439	146	UpdateDDL();
c3066940	147	while ( rawStreamV3->NextChannel() ){
	148	Int_t isector = rawStreamV3->GetSector(); // current sector
	149	Int_t iRow = rawStreamV3->GetRow(); // current row
	150	Int_t iPad = rawStreamV3->GetPad(); // current pad
	151	while ( rawStreamV3->NextBunch() ){
96bf9029	152	UInt_t startTbin = rawStreamV3->GetStartTimeBin();
c3066940	153	// Int_t endTbin = (Int_t)rawStreamV3->GetEndTimeBin();
96bf9029	154	Int_t bunchlength = rawStreamV3->GetBunchLength();
c3066940	155	const UShort_t *sig = rawStreamV3->GetSignals();
96bf9029	156	ProcessBunch(isector,iRow,iPad,bunchlength,startTbin,sig);
c3066940	157	for (Int_t iTimeBin = 0; iTimeBin<bunchlength; iTimeBin++){
	158	Float_t signal=(Float_t)sig[iTimeBin];
	159	// printf("%02d - %03d - %03d - %04d: %.1f\n",isector,iRow,iPad,startTbin,signal);
	160	Update(isector,iRow,iPad,startTbin--,signal);
	161	fPrevRCUId=fCurrRCUId;
	162	fPrevDDLNum=fCurrDDLNum;
	163	withInput = kTRUE;
	164	}
	165	}
	166	}
	167	}
	168	if (withInput){
	169	EndEvent();
	170	}
	171	return withInput;
	172	}
	173	//_____________________________________________________________________
7442bceb	174	Bool_t AliTPCCalibRawBase::ProcessEvent(AliRawReader * const rawReader)
c3066940	175	{
	176	//
	177	// Event processing loop - AliRawReader
	178	//
	179	AliRawEventHeaderBase* eventHeader = (AliRawEventHeaderBase*)rawReader->GetEventHeader();
	180	if (eventHeader){
	181	fTimeStamp = eventHeader->Get("Timestamp");
	182	fRunNumber = eventHeader->Get("RunNb");
	183	fEventType = eventHeader->Get("Type");
	184	}
78f17711	185	if (!fFirstTimeStamp) fFirstTimeStamp=fTimeStamp;
78f17711	186
c3066940	187	AliTPCRawStreamV3 rawStreamV3 = new AliTPCRawStreamV3(rawReader, (AliAltroMapping*)fMapping);
c3066940	188	Bool_t res=ProcessEvent(rawStreamV3);
78f17711	189
	190	fLastTimeStamp=fTimeStamp;
	191
c3066940	192	delete rawStreamV3;
	193	return res;
	194	}
	195	//_____________________________________________________________________
7442bceb	196	Bool_t AliTPCCalibRawBase::ProcessEvent(eventHeaderStruct * const event)
880c3382	197	{
	198	//
	199	// Event processing loop - date event
	200	//
5312f439	201
	202	fRunNumber=event->eventRunNb;
	203	fTimeStamp=event->eventTimestamp;
78f17711	204	if (!fFirstTimeStamp) fFirstTimeStamp=fTimeStamp;
78f17711	205	fLastTimeStamp=fTimeStamp;
5312f439	206	fEventType=event->eventType;
	207	AliRawReader rawReader = new AliRawReaderDate((void)event);
	208	AliTPCRawStreamV3 rawStreamV3 = new AliTPCRawStreamV3(rawReader, (AliAltroMapping*)fMapping);
	209	Bool_t result=ProcessEvent(rawStreamV3);
	210	delete rawStreamV3;
	211	delete rawReader;
	212	return result;
880c3382	213
	214	}
	215	//_____________________________________________________________________
	216	void AliTPCCalibRawBase::DumpToFile(const Char_t filename, const Char_t dir, Bool_t append)
	217	{
	218	//
	219	// Write class to file
	220	//
	221
	222	TString sDir(dir);
	223	TString option;
	224
	225	if ( append )
	226	option = "update";
	227	else
	228	option = "recreate";
	229
	230	TDirectory *backup = gDirectory;
	231	TFile f(filename,option.Data());
	232	f.cd();
	233	if ( !sDir.IsNull() ){
	234	f.mkdir(sDir.Data());
	235	f.cd(sDir);
	236	}
	237	this->Write();
	238	f.Close();
	239
	240	if ( backup ) backup->cd();
	241	}
	242	//_____________________________________________________________________
	243	TTreeSRedirector *AliTPCCalibRawBase::GetDebugStreamer(){
	244	//
	245	// Get Debug streamer
	246	// In case debug streamer not yet initialized and StreamLevel>0 create new one
	247	//
	248	if (fStreamLevel==0) return 0;
	249	if (fDebugStreamer) return fDebugStreamer;
	250	TString dsName;
	251	dsName=GetName();
	252	dsName+="Debug.root";
	253	dsName.ReplaceAll(" ","");
	254	fDebugStreamer = new TTreeSRedirector(dsName.Data());
	255	return fDebugStreamer;
	256	}
78f17711	257	//_____________________________________________________________________
	258	void AliTPCCalibRawBase::MergeBase(const AliTPCCalibRawBase *calib)
	259	{
	260	//
	261	// merge this with base
	262	//
	263	if (calib->fFirstTimeStamp<fFirstTimeStamp) fFirstTimeStamp=calib->fFirstTimeStamp;
	264	if (calib->fLastTimeStamp>fLastTimeStamp) fLastTimeStamp =calib->fLastTimeStamp;
	265	}
	266