- Removed 3D plot, used to much memory
- Added documentation
- New plot: phi vs eta
#include "AliCDBEntry.h"
#include "AliCDBManager.h"
#include "AliITSRecPoint.h"
#include "AliCDBEntry.h"
#include "AliCDBManager.h"
#include "AliITSRecPoint.h"
#include <TFile.h>
#include <TString.h>
#include "TObjString.h"
#include <TFile.h>
#include <TString.h>
#include "TObjString.h"
AliHLTITSClusterHistoComponent::AliHLTITSClusterHistoComponent()
:
fXY(NULL),
AliHLTITSClusterHistoComponent::AliHLTITSClusterHistoComponent()
:
fXY(NULL),
fCharge(NULL),
fPlotCharge(kFALSE),
fPlotXY(kTRUE),
fCharge(NULL),
fPlotCharge(kFALSE),
fPlotXY(kTRUE),
{
// see header file for class documentation
// or
{
// see header file for class documentation
// or
{
fPlotCharge=kFALSE;
fPlotXY=kTRUE;
{
fPlotCharge=kFALSE;
fPlotXY=kTRUE;
- if(fPlotCharge){fCharge = new TH1F("fCharge","Total Charge of clusters",4000,0,4000);}
+ if(fPlotCharge){fCharge = new TH1F("fCharge","Total Charge of clusters",2000,0,2000);}
if(fPlotXY){fXY = new TH2F("fXY","Global XY of ITS clusters",1600,-80,80,1600,-80,80);}
if(fPlotXY){fXY = new TH2F("fXY","Global XY of ITS clusters",1600,-80,80,1600,-80,80);}
- if(fPlotXYZ){fXYZ = new TH3F("fXYZ","Global XYZ of ITS clusters",1600,-80,80,1600,-80,80,2000,-100,100);}
+ if(fPlotPhieta){fPhieta = new TH2F("fPhieta","Global Phieta of ITS clusters",30,-1.5,1.5,60,0,2*TMath::Pi());}
int iResult=0;
TString configuration="";
int iResult=0;
TString configuration="";
// see header file for class documentation
if(fCharge!=NULL) delete fCharge;
if(fXY!=NULL) delete fXY;
// see header file for class documentation
if(fCharge!=NULL) delete fCharge;
if(fXY!=NULL) delete fXY;
- if(fXYZ!=NULL) delete fXYZ;
+ if(fPhieta!=NULL) delete fPhieta;
if(fPlotXY){
fXY->Fill(xyz[0],xyz[1]);
}
if(fPlotXY){
fXY->Fill(xyz[0],xyz[1]);
}
- if(fPlotXYZ){
- fXYZ->Fill(xyz[0],xyz[1],xyz[2]);
+ if(fPlotPhieta){
+ 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->Fill(theta,phi);
+ fPhieta->Fill(eta,phi);
}
if(fPlotCharge){
fCharge->Fill(recpoint.GetQ());
}
if(fPlotCharge){
fCharge->Fill(recpoint.GetQ());
AliHLTUInt32_t fSpecification = 0x0;
PushBack( (TObject*) fXY,kAliHLTDataTypeHistogram,fSpecification);
}
AliHLTUInt32_t fSpecification = 0x0;
PushBack( (TObject*) fXY,kAliHLTDataTypeHistogram,fSpecification);
}
AliHLTUInt32_t fSpecification = 0x0;
AliHLTUInt32_t fSpecification = 0x0;
- PushBack( (TObject*) fXYZ,kAliHLTDataTypeHistogram,fSpecification);
+ PushBack( (TObject*) fPhieta,kAliHLTDataTypeHistogram,fSpecification);
}
HLTInfo("ITSClusterHisto found %d Total Spacepoints", TotalSpacePoint);
}
HLTInfo("ITSClusterHisto found %d Total Spacepoints", TotalSpacePoint);
if (argument.CompareTo("-plot-all")==0) {
HLTInfo("Ploting all historgams");
fPlotXY = kTRUE;
if (argument.CompareTo("-plot-all")==0) {
HLTInfo("Ploting all historgams");
fPlotXY = kTRUE;
fPlotCharge = kTRUE;
continue;
}
fPlotCharge = kTRUE;
continue;
}
- else if (argument.CompareTo("-plot-xyz")==0) {
- HLTInfo("Ploting Global XYZ");
- //fPlotXYZ = kTRUE;
+ else if (argument.CompareTo("-plot-phieta")==0) {
+ HLTInfo("Ploting Global Phieta");
+ fPlotPhieta = kTRUE;
continue;
}
else if (argument.CompareTo("-plot-charge")==0) {
continue;
}
else if (argument.CompareTo("-plot-charge")==0) {
- if(!fCharge && fPlotCharge){fCharge = new TH1F("fCharge","Total Charge of clusters",4000,0,4000);}
+ if(!fCharge && fPlotCharge){fCharge = new TH1F("fCharge","Total Charge of clusters",2000,0,2000);}
if(!fXY && fPlotXY){fXY = new TH2F("fXY","Global XY of ITS clusters",1600,-80,80,1600,-80,80);}
if(!fXY && fPlotXY){fXY = new TH2F("fXY","Global XY of ITS clusters",1600,-80,80,1600,-80,80);}
- if(!fXYZ && fPlotXYZ){fXYZ = new TH3F("fXYZ","Global XYZ of ITS clusters",1600,-80,80,1600,-80,80,2000,-100,100);}
+ if(!fPhieta && fPlotPhieta){fPhieta = new TH2F("fPhieta","Global Phieta of ITS clusters",30,-1.5,1.5,60,0,2*TMath::Pi());}
//* ALICE Experiment at CERN, All rights reserved. *
//* See cxx source for full Copyright notice *
//* ALICE Experiment at CERN, All rights reserved. *
//* See cxx source for full Copyright notice *
-/** @file AliHLTITSQHistoComponent.h
+/** @file AliHLTITSClusterHistoComponent.h
@author Gaute Ovrebekk
@brief Component for ploting clusters
*/
@author Gaute Ovrebekk
@brief Component for ploting clusters
*/
#include "AliHLTProcessor.h"
#include "TH1F.h"
#include "TH2F.h"
#include "AliHLTProcessor.h"
#include "TH1F.h"
#include "TH2F.h"
#include "AliHLTITSSpacePointData.h"
#include "TClonesArray.h"
#include "AliITSRecPoint.h"
#include "AliHLTITSSpacePointData.h"
#include "TClonesArray.h"
#include "AliITSRecPoint.h"
-class AliHLTTPCConfMapper;
-
/**
* @class AliHLTITSQHistoComponent
* Component for ploting charge in clusters
*
/**
* @class AliHLTITSQHistoComponent
* Component for ploting charge in clusters
*
- * Component ID: \b ITSQHisto <br>
- * Library: \b libAliHLTITS.
+ * Component ID: \b ITSClusterHisto <br>
+ * Library: \b libAliHLTITS.so
*
* Mandatory arguments: <br>
*
*
* Optional arguments: <br>
*
* Mandatory arguments: <br>
*
*
* Optional arguments: <br>
+ * \li -plot-all <br>
+ * will plot all histograms in the component.
+ * \li -plot-xy <br>
+ * will plot a historgam of the x and y projection of the clusters.
+ * \li -plot-charge <i> multiplicity </i> <br>
+ * will plot the charge of the clusters.
+ * \li -plot-phieta <br>
+ * will plot the phi vs eta distrubution of the clusters.
* @ingroup alihlt_tpc_components
*/
class AliHLTITSClusterHistoComponent : public AliHLTProcessor
* @ingroup alihlt_tpc_components
*/
class AliHLTITSClusterHistoComponent : public AliHLTProcessor
int Configure(const char* arguments);
TH2F * fXY; //! transient
int Configure(const char* arguments);
TH2F * fXY; //! transient
- TH3F * fXYZ; //! transient
+ TH2F * fPhieta; //! transient
TH1F * fCharge; //! transient
Bool_t fPlotCharge; //! transient
Bool_t fPlotXY; //! transient
TH1F * fCharge; //! transient
Bool_t fPlotCharge; //! transient
Bool_t fPlotXY; //! transient
- Bool_t fPlotXYZ; //! transient
+ Bool_t fPlotPhieta; //! transient
- ClassDef(AliHLTITSClusterHistoComponent, 0);
+ ClassDef(AliHLTITSClusterHistoComponent, 1);