[u/mrichter/AliRoot.git] / HLT / ITS / AliHLTITSClusterHistoComponent.cxx

// $Id$
//**************************************************************************
//* This file is property of and copyright by the ALICE HLT Project        * 
//* ALICE Experiment at CERN, All rights reserved.                         *
//*                                                                        *
//* Primary Authors: Gaute Ovrebekk <ovrebekk@ift.uib.no>                  *
//*                  for The ALICE HLT Project.                            *
//*                                                                        *
//* 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.                  *
//**************************************************************************

/** @file   AliHLTITSClusterHistoComponent.cxx
    @author Gaute Ovrebekk
    @brief  Component for plotting charge in clusters
*/

#include "AliHLTITSClusterHistoComponent.h"
#include "AliHLTITSClusterDataFormat.h"
#include "AliCDBEntry.h"
#include "AliCDBManager.h"
#include "AliGeomManager.h"
#include "AliITSRecPoint.h"
#include "TMath.h"
#include <TFile.h>
#include <TString.h>
#include "TObjString.h"
#include "TObjArray.h"

using namespace std;

/** ROOT macro for the implementation of ROOT specific class methods */
ClassImp(AliHLTITSClusterHistoComponent)

AliHLTITSClusterHistoComponent::AliHLTITSClusterHistoComponent()
  :
  fXY(NULL),                     
  fPhieta(NULL),                   
  fCharge(NULL),   
  fPlotCharge(kFALSE),   
  fPlotXYPhiEta(kTRUE),
  fOutputSize(20000)
{
  // see header file for class documentation
  // or
  // refer to README to build package
  // or
  // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
}

AliHLTITSClusterHistoComponent::~AliHLTITSClusterHistoComponent(){
  // see header file for class documentation
}

// Public functions to implement AliHLTComponent's interface.
// These functions are required for the registration process

const char* AliHLTITSClusterHistoComponent::GetComponentID(){
  // see header file for class documentation
  return "ITSClusterHisto";
}

void AliHLTITSClusterHistoComponent::GetInputDataTypes(AliHLTComponentDataTypeList& list){
  // see header file for class documentation
  list.clear();
  list.push_back( kAliHLTDataTypeClusters|kAliHLTDataOriginITSSPD );
  list.push_back( kAliHLTDataTypeClusters|kAliHLTDataOriginITSSDD );
  list.push_back( kAliHLTDataTypeClusters|kAliHLTDataOriginITSSSD );
  list.push_back( kAliHLTDataTypeClusters|kAliHLTDataOriginITS );
}

AliHLTComponentDataType AliHLTITSClusterHistoComponent::GetOutputDataType(){
  // see header file for class documentation
  return kAliHLTDataTypeHistogram;
}

void AliHLTITSClusterHistoComponent::GetOutputDataSize( unsigned long& constBase, double& inputMultiplier ){
  // see header file for class documentation
  // XXX TODO: Find more realistic values.
  constBase = fOutputSize;
  inputMultiplier = 20;
}

AliHLTComponent* AliHLTITSClusterHistoComponent::Spawn(){
  // see header file for class documentation
  return new AliHLTITSClusterHistoComponent;
}

int AliHLTITSClusterHistoComponent::DoInit( int argc, const char** argv ){  
  // see header file for class documentation
     
  if(AliGeomManager::GetGeometry()==NULL){
    AliGeomManager::LoadGeometry();
  }
  
  fPlotCharge   = kFALSE;   
  fPlotXYPhiEta = kTRUE;
     
  if(fPlotCharge) fCharge = new TH1F("fCharge","Total Charge of clusters",2000,0,2000);
  if(fPlotXYPhiEta){
     fXY = new TH2F("fXY","Global XY of ITS clusters",1600,-80,80,1600,-80,80);
     TString name;
     TString title;
     fPhieta = new TH2F*[6];
     for (Int_t iLay=0;iLay<6;iLay++){
     	     name.Form("Phi_vs_Eta_ITS_Layer%d",iLay+1);
     	     title.Form("Phi vs Eta - ITS Layer %d",iLay+1);
     	     fPhieta[iLay]=new TH2F(name,title,60,-1.5,1.5,60,0.,2*TMath::Pi());
     	     fPhieta[iLay]->GetXaxis()->SetTitle("Pseudorapidity");
     	     fPhieta[iLay]->GetYaxis()->SetTitle("#varphi [rad]");
     }
  }
  
  int iResult=0;
  TString configuration="";
  TString argument="";

  for(int i=0; i<argc && iResult>=0; i++){
    argument=argv[i];
    if (!configuration.IsNull()) configuration+=" ";
    configuration+=argument;
  }
  
  if(!configuration.IsNull()){
    iResult=Configure(configuration.Data());
  }  

  return iResult; 
}
  
int AliHLTITSClusterHistoComponent::DoDeinit(){
  // see header file for class documentation
  if(fCharge!=NULL) delete fCharge;
  if(fXY!=NULL) {
     delete fXY;
     for(Int_t i=0;i<6;i++) delete fPhieta[5-i];
  }
  return 0;
}

int AliHLTITSClusterHistoComponent::DoEvent(const AliHLTComponentEventData& /*evtData*/, AliHLTComponentTriggerData& /*trigData*/){
  // see header file for class documentation
   
  static Int_t event = 0;
  event++;
  int TotalSpacePoint = 0;
  
  const AliHLTComponentBlockData* iter = NULL;
  
  AliHLTUInt32_t eventType=gkAliEventTypeUnknown;
  if(!IsDataEvent(&eventType) && eventType!=gkAliEventTypeEndOfRun) {
    // publish the histograms at the end of run
    return 0;
  }
  
  for(iter = GetFirstInputBlock(kAliHLTDataTypeClusters); eventType!=gkAliEventTypeEndOfRun && iter != NULL; iter = GetNextInputBlock()){
  
      if(iter->fDataType!=(kAliHLTAnyDataType|kAliHLTDataOriginITSSPD) && 
         iter->fDataType!=(kAliHLTAnyDataType|kAliHLTDataOriginITSSDD) && 
         iter->fDataType!=(kAliHLTAnyDataType|kAliHLTDataOriginITSSSD) &&
	 iter->fDataType!=(kAliHLTAnyDataType|kAliHLTDataOriginITS))	
         continue;

      const AliHLTITSClusterData* clusterData = (const AliHLTITSClusterData*) iter->fPtr;
      Int_t nSpacepoint = (Int_t) clusterData->fSpacePointCnt;
      TotalSpacePoint += nSpacepoint;
      AliHLTITSSpacePointData *clusters = (AliHLTITSSpacePointData*) clusterData->fSpacePoints;

      for(int i=0;i<nSpacepoint;i++){
         
	  Int_t lab[4]   = {0,0,0,0};
          Float_t hit[6] = {0,0,0,0,0,0};
          Int_t info[3]  = {0,0,0};
          
          lab[0]  = clusters[i].fTracks[0];
          lab[1]  = clusters[i].fTracks[1];
          lab[2]  = clusters[i].fTracks[2];
          lab[3]  = clusters[i].fIndex;
          hit[0]  = clusters[i].fY;
          hit[1]  = clusters[i].fZ;
          hit[2]  = clusters[i].fSigmaY2;
          hit[3]  = clusters[i].fSigmaZ2;
          hit[4]  = clusters[i].fQ;
          hit[5]  = clusters[i].fSigmaYZ;
          info[0] = clusters[i].fNy;
          info[1] = clusters[i].fNz;
          info[2] = clusters[i].fLayer;
 
          Float_t xyz[3];
          AliITSRecPoint recpoint(lab,hit,info);
          recpoint.GetGlobalXYZ(xyz);
    	  Int_t layer = recpoint.GetLayer();
          
	  if(fPlotXYPhiEta){
    	     fXY->Fill(xyz[0],xyz[1]);
    	     Float_t rad=TMath::Sqrt(xyz[0]*xyz[0]+xyz[1]*xyz[1]); 
    	     Float_t theta=TMath::ATan2(rad,xyz[2]);
    	     Float_t eta=-1*TMath::Log(TMath::Tan(theta/2.0));
    	     Float_t phi=TMath::ATan2(xyz[1],xyz[0]);
    	     if(phi<0.0){phi=2 * TMath::Pi() - TMath::Abs(phi);} 
    	     fPhieta[layer]->Fill(eta,phi);
          }
          if(fPlotCharge) fCharge->Fill(recpoint.GetQ());
    }
  }
  
  if(fPlotCharge){
     AliHLTUInt32_t fSpecification = 0x0;
     if (PushBack( (TObject*) fCharge,kAliHLTDataTypeHistogram|kAliHLTDataOriginITS,fSpecification)==-ENOSPC) {
       fOutputSize+=GetLastObjectSize();
     }
  }
  if(fPlotXYPhiEta){
     AliHLTUInt32_t fSpecification = 0x0;
     if (PushBack( (TObject*) fXY,kAliHLTDataTypeHistogram|kAliHLTDataOriginITS,fSpecification)==-ENOSPC) {
       fOutputSize+=GetLastObjectSize();
     }
     for(Int_t ii=0;ii<6;ii++) 
       if (PushBack( (TObject*) fPhieta[ii],kAliHLTDataTypeHistogram|kAliHLTDataOriginITS,fSpecification)==-ENOSPC) {
	 fOutputSize+=GetLastObjectSize();
       }
  }
  
  HLTDebug("ITSClusterHisto found %d Total Spacepoints", TotalSpacePoint);
  
  return 0;
}

int AliHLTITSClusterHistoComponent::Configure(const char* arguments){
  // see header file for class documentation
   
  int iResult=0;
  
  if (!arguments) return iResult;
  
  TString allArgs=arguments;
  TString argument;
  
  TObjArray* pTokens=allArgs.Tokenize(" ");
  
  if (pTokens) {
    for (int i=0; i<pTokens->GetEntries() && iResult>=0; i++) {
      argument=((TObjString*)pTokens->At(i))->GetString();
      if (argument.IsNull()) continue;
      
      if (argument.CompareTo("-plot-all")==0) {
	HLTInfo("Ploting all historgams");
	fPlotXYPhiEta = kTRUE;
	fPlotCharge = kTRUE;
	continue;
      }
      
      else if (argument.CompareTo("-plot-xy")==0) {
	HLTInfo("Ploting Global XY");
	fPlotXYPhiEta = kTRUE;
	continue;
      }

      else if (argument.CompareTo("-plot-charge")==0) {
	HLTInfo("Ploting charge of clusters");
	fPlotCharge = kTRUE;
	continue;
      }

      else {
	HLTError("unknown argument %s", argument.Data());
	iResult=-EINVAL;
	break;
      }
    }
    delete pTokens;
  }
  
  if(!fCharge && fPlotCharge){fCharge = new TH1F("fCharge","Total Charge of clusters",2000,0,2000);}
  if(!fXY && fPlotXYPhiEta){
	fXY = new TH2F("fXY","Global XY of ITS clusters",1600,-80,80,1600,-80,80);
  	TString name;
	TString title;
	fPhieta = new TH2F*[6];
	for (Int_t iLay=0;iLay<6;iLay++) {
	        name.Form("Phi_vs_Eta_ITS_Layer%d",iLay+1);
	        title.Form("Phi vs Eta - ITS Layer %d",iLay+1);
		fPhieta[iLay]=new TH2F(name,title,30,-1.5,1.5,200,0.,2*TMath::Pi());
		fPhieta[iLay]->GetXaxis()->SetTitle("Pseudorapidity");
		fPhieta[iLay]->GetYaxis()->SetTitle("#varphi [rad]");
 	}
  }
  
  return iResult;
}

int AliHLTITSClusterHistoComponent::Reconfigure(const char* cdbEntry, const char* chainId){
  // see header file for class documentation
  int iResult=0;
  
  const char* path="HLT/ConfigITS/HistoComponent";
  const char* defaultNotify="";
  if (cdbEntry) {
    path=cdbEntry;
    defaultNotify=" (default)";
  }
  if (path) {
    HLTInfo("reconfigure from entry %s%s, chain id %s", path, defaultNotify,(chainId!=NULL && chainId[0]!=0)?chainId:"<none>");
    AliCDBEntry *pEntry = AliCDBManager::Instance()->Get(path);
    if (pEntry) {
      TObjString* pString=dynamic_cast<TObjString*>(pEntry->GetObject());
      if (pString) {
	HLTInfo("received configuration object string: \'%s\'", pString->GetString().Data());
	iResult=Configure(pString->GetString().Data());
      } else {
	HLTError("configuration object \"%s\" has wrong type, required TObjString", path);
      }
    } else {
      HLTError("can not fetch object \"%s\" from CDB", path);
    }
  }
  
  return iResult;
}
Commit	Line	Data
c0b68f8e	1	// $Id$
	2	//**************************************************************************
	3	//* This file is property of and copyright by the ALICE HLT Project *
	4	//* ALICE Experiment at CERN, All rights reserved. *
	5	//* *
	6	//* Primary Authors: Gaute Ovrebekk <ovrebekk@ift.uib.no> *
	7	//* for The ALICE HLT Project. *
	8	//* *
	9	//* Permission to use, copy, modify and distribute this software and its *
	10	//* documentation strictly for non-commercial purposes is hereby granted *
	11	//* without fee, provided that the above copyright notice appears in all *
	12	//* copies and that both the copyright notice and this permission notice *
	13	//* appear in the supporting documentation. The authors make no claims *
	14	//* about the suitability of this software for any purpose. It is *
	15	//* provided "as is" without express or implied warranty. *
	16	//**************************************************************************
	17
	18	/** @file AliHLTITSClusterHistoComponent.cxx
	19	@author Gaute Ovrebekk
f44c9362	20	@brief Component for plotting charge in clusters
c0b68f8e	21	*/
c0b68f8e	22
c0b68f8e	23	#include "AliHLTITSClusterHistoComponent.h"
	24	#include "AliHLTITSClusterDataFormat.h"
	25	#include "AliCDBEntry.h"
	26	#include "AliCDBManager.h"
f44c9362	27	#include "AliGeomManager.h"
c0b68f8e	28	#include "AliITSRecPoint.h"
91edcbad	29	#include "TMath.h"
c0b68f8e	30	#include <TFile.h>
	31	#include <TString.h>
	32	#include "TObjString.h"
	33	#include "TObjArray.h"
	34
8e0f6ca3	35	using namespace std;
8e0f6ca3	36
c0b68f8e	37	/** ROOT macro for the implementation of ROOT specific class methods */
	38	ClassImp(AliHLTITSClusterHistoComponent)
	39
	40	AliHLTITSClusterHistoComponent::AliHLTITSClusterHistoComponent()
f44c9362	41	:
f44c9362	42	fXY(NULL),
91edcbad	43	fPhieta(NULL),
c0b68f8e	44	fCharge(NULL),
c0b68f8e	45	fPlotCharge(kFALSE),
59918cb1	46	fPlotXYPhiEta(kTRUE),
618f422f	47	fOutputSize(20000)
c0b68f8e	48	{
	49	// see header file for class documentation
	50	// or
	51	// refer to README to build package
	52	// or
	53	// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
c0b68f8e	54	}
c0b68f8e	55
f44c9362	56	AliHLTITSClusterHistoComponent::~AliHLTITSClusterHistoComponent(){
c0b68f8e	57	// see header file for class documentation
	58	}
	59
	60	// Public functions to implement AliHLTComponent's interface.
	61	// These functions are required for the registration process
	62
f44c9362	63	const char* AliHLTITSClusterHistoComponent::GetComponentID(){
c0b68f8e	64	// see header file for class documentation
c0b68f8e	65	return "ITSClusterHisto";
	66	}
	67
f44c9362	68	void AliHLTITSClusterHistoComponent::GetInputDataTypes(AliHLTComponentDataTypeList& list){
c0b68f8e	69	// see header file for class documentation
c0b68f8e	70	list.clear();
8ce2efb0	71	list.push_back( kAliHLTDataTypeClusters\|kAliHLTDataOriginITSSPD );
	72	list.push_back( kAliHLTDataTypeClusters\|kAliHLTDataOriginITSSDD );
	73	list.push_back( kAliHLTDataTypeClusters\|kAliHLTDataOriginITSSSD );
618f422f	74	list.push_back( kAliHLTDataTypeClusters\|kAliHLTDataOriginITS );
c0b68f8e	75	}
c0b68f8e	76
f44c9362	77	AliHLTComponentDataType AliHLTITSClusterHistoComponent::GetOutputDataType(){
c0b68f8e	78	// see header file for class documentation
c0b68f8e	79	return kAliHLTDataTypeHistogram;
c0b68f8e	80	}
c0b68f8e	81
f44c9362	82	void AliHLTITSClusterHistoComponent::GetOutputDataSize( unsigned long& constBase, double& inputMultiplier ){
c0b68f8e	83	// see header file for class documentation
c0b68f8e	84	// XXX TODO: Find more realistic values.
59918cb1	85	constBase = fOutputSize;
618f422f	86	inputMultiplier = 20;
c0b68f8e	87	}
c0b68f8e	88
f44c9362	89	AliHLTComponent* AliHLTITSClusterHistoComponent::Spawn(){
c0b68f8e	90	// see header file for class documentation
	91	return new AliHLTITSClusterHistoComponent;
	92	}
	93
f44c9362	94	int AliHLTITSClusterHistoComponent::DoInit( int argc, const char** argv ){
	95	// see header file for class documentation
	96
	97	if(AliGeomManager::GetGeometry()==NULL){
	98	AliGeomManager::LoadGeometry();
	99	}
	100
	101	fPlotCharge = kFALSE;
	102	fPlotXYPhiEta = kTRUE;
c0b68f8e	103
f44c9362	104	if(fPlotCharge) fCharge = new TH1F("fCharge","Total Charge of clusters",2000,0,2000);
0b26a583	105	if(fPlotXYPhiEta){
f44c9362	106	fXY = new TH2F("fXY","Global XY of ITS clusters",1600,-80,80,1600,-80,80);
4d4165a5	107	TString name;
4d4165a5	108	TString title;
f44c9362	109	fPhieta = new TH2F*[6];
f44c9362	110	for (Int_t iLay=0;iLay<6;iLay++){
4d4165a5	111	name.Form("Phi_vs_Eta_ITS_Layer%d",iLay+1);
4d4165a5	112	title.Form("Phi vs Eta - ITS Layer %d",iLay+1);
f44c9362	113	fPhieta[iLay]=new TH2F(name,title,60,-1.5,1.5,60,0.,2*TMath::Pi());
	114	fPhieta[iLay]->GetXaxis()->SetTitle("Pseudorapidity");
	115	fPhieta[iLay]->GetYaxis()->SetTitle("#varphi [rad]");
	116	}
0b26a583	117	}
c0b68f8e	118
	119	int iResult=0;
	120	TString configuration="";
	121	TString argument="";
f44c9362	122
f44c9362	123	for(int i=0; i<argc && iResult>=0; i++){
c0b68f8e	124	argument=argv[i];
	125	if (!configuration.IsNull()) configuration+=" ";
	126	configuration+=argument;
	127	}
	128
f44c9362	129	if(!configuration.IsNull()){
c0b68f8e	130	iResult=Configure(configuration.Data());
	131	}
	132
	133	return iResult;
	134	}
	135
f44c9362	136	int AliHLTITSClusterHistoComponent::DoDeinit(){
c0b68f8e	137	// see header file for class documentation
c0b68f8e	138	if(fCharge!=NULL) delete fCharge;
0b26a583	139	if(fXY!=NULL) {
f44c9362	140	delete fXY;
f44c9362	141	for(Int_t i=0;i<6;i++) delete fPhieta[5-i];
0b26a583	142	}
c0b68f8e	143	return 0;
	144	}
	145
f44c9362	146	int AliHLTITSClusterHistoComponent::DoEvent(const AliHLTComponentEventData& /evtData/, AliHLTComponentTriggerData& /trigData/){
	147	// see header file for class documentation
	148
0b26a583	149	static Int_t event = 0;
0b26a583	150	event++;
c0b68f8e	151	int TotalSpacePoint = 0;
	152
	153	const AliHLTComponentBlockData* iter = NULL;
	154
59918cb1	155	AliHLTUInt32_t eventType=gkAliEventTypeUnknown;
	156	if(!IsDataEvent(&eventType) && eventType!=gkAliEventTypeEndOfRun) {
	157	// publish the histograms at the end of run
	158	return 0;
	159	}
f44c9362	160
59918cb1	161	for(iter = GetFirstInputBlock(kAliHLTDataTypeClusters); eventType!=gkAliEventTypeEndOfRun && iter != NULL; iter = GetNextInputBlock()){
c0b68f8e	162
f44c9362	163	if(iter->fDataType!=(kAliHLTAnyDataType\|kAliHLTDataOriginITSSPD) &&
f44c9362	164	iter->fDataType!=(kAliHLTAnyDataType\|kAliHLTDataOriginITSSDD) &&
618f422f	165	iter->fDataType!=(kAliHLTAnyDataType\|kAliHLTDataOriginITSSSD) &&
618f422f	166	iter->fDataType!=(kAliHLTAnyDataType\|kAliHLTDataOriginITS))
f44c9362	167	continue;
8ce2efb0	168
f44c9362	169	const AliHLTITSClusterData* clusterData = (const AliHLTITSClusterData*) iter->fPtr;
	170	Int_t nSpacepoint = (Int_t) clusterData->fSpacePointCnt;
	171	TotalSpacePoint += nSpacepoint;
	172	AliHLTITSSpacePointData clusters = (AliHLTITSSpacePointData) clusterData->fSpacePoints;
c0b68f8e	173
f44c9362	174	for(int i=0;i<nSpacepoint;i++){
	175
	176	Int_t lab[4] = {0,0,0,0};
	177	Float_t hit[6] = {0,0,0,0,0,0};
	178	Int_t info[3] = {0,0,0};
	179
	180	lab[0] = clusters[i].fTracks[0];
	181	lab[1] = clusters[i].fTracks[1];
	182	lab[2] = clusters[i].fTracks[2];
	183	lab[3] = clusters[i].fIndex;
	184	hit[0] = clusters[i].fY;
	185	hit[1] = clusters[i].fZ;
	186	hit[2] = clusters[i].fSigmaY2;
	187	hit[3] = clusters[i].fSigmaZ2;
	188	hit[4] = clusters[i].fQ;
	189	hit[5] = clusters[i].fSigmaYZ;
	190	info[0] = clusters[i].fNy;
	191	info[1] = clusters[i].fNz;
	192	info[2] = clusters[i].fLayer;
	193
	194	Float_t xyz[3];
	195	AliITSRecPoint recpoint(lab,hit,info);
	196	recpoint.GetGlobalXYZ(xyz);
	197	Int_t layer = recpoint.GetLayer();
	198
	199	if(fPlotXYPhiEta){
	200	fXY->Fill(xyz[0],xyz[1]);
	201	Float_t rad=TMath::Sqrt(xyz[0]xyz[0]+xyz[1]xyz[1]);
	202	Float_t theta=TMath::ATan2(rad,xyz[2]);
	203	Float_t eta=-1*TMath::Log(TMath::Tan(theta/2.0));
	204	Float_t phi=TMath::ATan2(xyz[1],xyz[0]);
	205	if(phi<0.0){phi=2 * TMath::Pi() - TMath::Abs(phi);}
	206	fPhieta[layer]->Fill(eta,phi);
	207	}
	208	if(fPlotCharge) fCharge->Fill(recpoint.GetQ());
c0b68f8e	209	}
	210	}
	211
	212	if(fPlotCharge){
f44c9362	213	AliHLTUInt32_t fSpecification = 0x0;
59918cb1	214	if (PushBack( (TObject*) fCharge,kAliHLTDataTypeHistogram\|kAliHLTDataOriginITS,fSpecification)==-ENOSPC) {
	215	fOutputSize+=GetLastObjectSize();
	216	}
c0b68f8e	217	}
0b26a583	218	if(fPlotXYPhiEta){
f44c9362	219	AliHLTUInt32_t fSpecification = 0x0;
59918cb1	220	if (PushBack( (TObject*) fXY,kAliHLTDataTypeHistogram\|kAliHLTDataOriginITS,fSpecification)==-ENOSPC) {
	221	fOutputSize+=GetLastObjectSize();
	222	}
	223	for(Int_t ii=0;ii<6;ii++)
	224	if (PushBack( (TObject*) fPhieta[ii],kAliHLTDataTypeHistogram\|kAliHLTDataOriginITS,fSpecification)==-ENOSPC) {
	225	fOutputSize+=GetLastObjectSize();
	226	}
c0b68f8e	227	}
c0b68f8e	228
f44c9362	229	HLTDebug("ITSClusterHisto found %d Total Spacepoints", TotalSpacePoint);
c0b68f8e	230
	231	return 0;
	232	}
	233
f44c9362	234	int AliHLTITSClusterHistoComponent::Configure(const char* arguments){
	235	// see header file for class documentation
	236
c0b68f8e	237	int iResult=0;
	238
	239	if (!arguments) return iResult;
	240
	241	TString allArgs=arguments;
	242	TString argument;
	243
	244	TObjArray* pTokens=allArgs.Tokenize(" ");
	245
	246	if (pTokens) {
	247	for (int i=0; i<pTokens->GetEntries() && iResult>=0; i++) {
	248	argument=((TObjString*)pTokens->At(i))->GetString();
	249	if (argument.IsNull()) continue;
	250
	251	if (argument.CompareTo("-plot-all")==0) {
	252	HLTInfo("Ploting all historgams");
0b26a583	253	fPlotXYPhiEta = kTRUE;
c0b68f8e	254	fPlotCharge = kTRUE;
	255	continue;
	256	}
	257
	258	else if (argument.CompareTo("-plot-xy")==0) {
	259	HLTInfo("Ploting Global XY");
0b26a583	260	fPlotXYPhiEta = kTRUE;
c0b68f8e	261	continue;
	262	}
	263
c0b68f8e	264	else if (argument.CompareTo("-plot-charge")==0) {
	265	HLTInfo("Ploting charge of clusters");
	266	fPlotCharge = kTRUE;
	267	continue;
	268	}
	269
	270	else {
	271	HLTError("unknown argument %s", argument.Data());
	272	iResult=-EINVAL;
	273	break;
	274	}
	275	}
	276	delete pTokens;
	277	}
	278
91edcbad	279	if(!fCharge && fPlotCharge){fCharge = new TH1F("fCharge","Total Charge of clusters",2000,0,2000);}
0b26a583	280	if(!fXY && fPlotXYPhiEta){
0b26a583	281	fXY = new TH2F("fXY","Global XY of ITS clusters",1600,-80,80,1600,-80,80);
4d4165a5	282	TString name;
4d4165a5	283	TString title;
0b26a583	284	fPhieta = new TH2F*[6];
0b26a583	285	for (Int_t iLay=0;iLay<6;iLay++) {
4d4165a5	286	name.Form("Phi_vs_Eta_ITS_Layer%d",iLay+1);
4d4165a5	287	title.Form("Phi vs Eta - ITS Layer %d",iLay+1);
0b26a583	288	fPhieta[iLay]=new TH2F(name,title,30,-1.5,1.5,200,0.,2*TMath::Pi());
	289	fPhieta[iLay]->GetXaxis()->SetTitle("Pseudorapidity");
	290	fPhieta[iLay]->GetYaxis()->SetTitle("#varphi [rad]");
	291	}
	292	}
c0b68f8e	293
	294	return iResult;
	295	}
	296
f44c9362	297	int AliHLTITSClusterHistoComponent::Reconfigure(const char* cdbEntry, const char* chainId){
c0b68f8e	298	// see header file for class documentation
	299	int iResult=0;
	300
	301	const char* path="HLT/ConfigITS/HistoComponent";
	302	const char* defaultNotify="";
	303	if (cdbEntry) {
	304	path=cdbEntry;
	305	defaultNotify=" (default)";
	306	}
	307	if (path) {
	308	HLTInfo("reconfigure from entry %s%s, chain id %s", path, defaultNotify,(chainId!=NULL && chainId[0]!=0)?chainId:"<none>");
	309	AliCDBEntry *pEntry = AliCDBManager::Instance()->Get(path);
	310	if (pEntry) {
	311	TObjString* pString=dynamic_cast<TObjString*>(pEntry->GetObject());
	312	if (pString) {
	313	HLTInfo("received configuration object string: \'%s\'", pString->GetString().Data());
	314	iResult=Configure(pString->GetString().Data());
	315	} else {
	316	HLTError("configuration object \"%s\" has wrong type, required TObjString", path);
	317	}
	318	} else {
	319	HLTError("can not fetch object \"%s\" from CDB", path);
	320	}
	321	}
	322
	323	return iResult;
	324	}