// --- ROOT system ---
class TFile;
#include <TFolder.h>
-#include <TROOT.h>
#include <TTree.h>
#include <TVirtualMC.h>
// --- Standard library ---
-#include <Riostream.h>
// --- AliRoot header files ---
#include "AliMagF.h"
#include "AliPHOS.h"
#include "AliPHOSGeometry.h"
#include "AliPHOSLoader.h"
-#include "AliPHOSQAChecker.h"
-#include "AliRun.h"
+
ClassImp(AliPHOS)
//____________________________________________________________________________
//____________________________________________________________________________
void AliPHOS::SetTreeAddress()
{
-
+ // Links Hits in the Tree to Hits array
TBranch *branch;
char branchname[20];
sprintf(branchname,"%s",GetName());
//
//*-- Author: Laurent Aphecetche & Yves Schutz (SUBATECH)
-#include <stdlib.h>
// --- ROOT system ---
-
class TString ;
class TTask ;
class TFolder ;
+class TTree ;
// --- AliRoot header files ---
-
#include "AliDetector.h"
class AliPHOSGeometry ;
class AliPHOSQAChecker ;
virtual const TString Version() const {return TString(" ") ; }
virtual void WriteQA() ;
AliPHOS & operator = (const AliPHOS & /*rvalue*/) {
- // assignement operator requested by coding convention
- // but not needed
- Fatal("operator =", "not implemented") ;
- return *this ;
- }
+ Fatal("operator =", "not implemented") ; return *this ; }
virtual AliLoader* MakeLoader(const char* topfoldername);
#include "TH2.h"
#include "TParticle.h"
#include "TClonesArray.h"
-#include "TTree.h"
#include "TMath.h"
#include "TROOT.h"
-#include "TFolder.h"
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliRun.h"
+//#include "AliRun.h"
#include "AliStack.h"
-#include "AliPHOSv1.h"
+#include "AliPHOSGeometry.h"
#include "AliPHOSAnalyze.h"
#include "AliPHOSDigit.h"
#include "AliPHOSSDigitizer.h"
+#include "AliPHOSEmcRecPoint.h"
+#include "AliPHOSCpvRecPoint.h"
#include "AliPHOSTrackSegment.h"
#include "AliPHOSRecParticle.h"
-#include "AliPHOSCpvRecPoint.h"
#include "AliPHOSLoader.h"
//////////////////////////////////////////////////////////////////////////////
// --- ROOT system ---
-#include "TFile.h"
+
// --- Standard library ---
#include <fstream>
#include "TROOT.h"
#include "TF1.h"
#include "TObjString.h"
+
// --- Standard library ---
// --- AliRoot header files ---
cmngr->WriteData(&gains) ;
}
//_____________________________________________________________________________
-void AliPHOSCalibrator::Print(const Option_t *)const
+void AliPHOSCalibrator::Print()const
{
// prints everything
printf("--------------AliPHOSCalibrator-----------------\n") ;
void PlotGain(Int_t channel) ; //Plot histo with gains for a channel
void PlotGains() ; //Plot all gains
- virtual void Print(const Option_t *) const ;
+ virtual void Print() const ;
TH1F * PedestalHisto(Int_t channel)
{return dynamic_cast<TH1F* >(fPedHistos->At(channel)) ;} ;
//////////////////////////////////////////////////////////////////////////////
// --- ROOT system ---
-#include "TGeometry.h"
-#include "TDirectory.h"
-#include "TFile.h"
-#include "TTree.h"
// --- Standard library ---
-#include <stdlib.h>
// --- AliRoot header files ---
-#include "AliRun.h"
#include "AliPHOSClusterizer.h"
-#include "AliHeader.h"
-#include "AliPHOSGetter.h"
-#include "AliPHOSSDigitizer.h"
-#include "AliPHOSDigitizer.h"
ClassImp(AliPHOSClusterizer)
// --- ROOT system ---
-#include "TROOT.h"
-#include "TFile.h"
-#include "TFolder.h"
#include "TMath.h"
#include "TMinuit.h"
#include "TTree.h"
-#include "TSystem.h"
#include "TBenchmark.h"
// --- Standard library ---
-//#include <iostream.h>
// --- AliRoot header files ---
#include "AliPHOSGetter.h"
+#include "AliPHOSGeometry.h"
#include "AliPHOSClusterizerv1.h"
#include "AliPHOSEmcRecPoint.h"
#include "AliPHOSCpvRecPoint.h"
#include "AliPHOSDigit.h"
#include "AliPHOSDigitizer.h"
-#include "AliPHOS.h"
-#include "AliPHOSCalibrationDB.h"
ClassImp(AliPHOSClusterizerv1)
//Now the same for CPV
for(index = 0; index < cpvRecPoints->GetEntries(); index++)
- dynamic_cast<AliPHOSRecPoint *>( cpvRecPoints->At(index) )->EvalAll(fW0CPV,digits) ;
+ dynamic_cast<AliPHOSRecPoint *>( cpvRecPoints->At(index) )->EvalAll(digits) ;
cpvRecPoints->Sort() ;
// --- ROOT system ---
#include "TBenchmark.h"
-#include "TROOT.h"
// --- Standard library ---
#include "AliPHOSEvalRecPoint.h"
#include "AliPHOSRecCpvManager.h"
#include "AliPHOSRecEmcManager.h"
+#include "AliPHOSGeometry.h"
ClassImp(AliPHOSClusterizerv2)
return 0 ;
}
//____________________________________________________________________________
-void AliPHOSConTableDB::Print(Option_t *)const {
+void AliPHOSConTableDB::Print()const {
//prints configuraion
TString message ;
//Transforms AbsId number in PHOS into channel number in prototype
Int_t AbsId2Raw(Int_t AbsId)const ;
- virtual void Print(Option_t *) const ;
+ virtual void Print() const ;
AliPHOSConTableDB& operator=(const AliPHOSConTableDB& cdb) ;
// 30 October 2000
// --- ROOT system ---
-#include "TPad.h"
-#include "TH2.h"
+
#include "TMath.h"
-#include "TCanvas.h"
#include "TClonesArray.h"
// --- Standard library ---
// --- AliRoot header files ---
-
+#include "AliPHOSGeometry.h"
+#include "AliPHOSDigit.h"
#include "AliPHOSCpvRecPoint.h"
#include "AliPHOSLoader.h"
+
ClassImp(AliPHOSCpvRecPoint)
//____________________________________________________________________________
//____________________________________________________________________________
-void AliPHOSCpvRecPoint::Print(Option_t *)
+void AliPHOSCpvRecPoint::Print()
{
// Print the list of digits belonging to the cluster
// 30 October 2000
// --- ROOT system ---
-#include "TObject.h"
-#include "TArrayI.h"
+//#include "TObject.h"
+//#include "TArrayI.h"
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliPHOSDigit.h"
+//#include "AliPHOSDigit.h"
#include "AliPHOSEmcRecPoint.h"
-#include "AliPHOSGeometry.h"
+//#include "AliPHOSGeometry.h"
class AliPHOSCpvRecPoint : public AliPHOSEmcRecPoint {
AliPHOSCpvRecPoint() ;
AliPHOSCpvRecPoint(const char * opt) ;
AliPHOSCpvRecPoint(const AliPHOSCpvRecPoint & rp) : AliPHOSEmcRecPoint(rp) {
- // cpy ctor requested by Coding Convention
- // but not yet needed
- assert(0==1) ;
+ Fatal("cpy ctor", "not implemented") ;
}
virtual ~AliPHOSCpvRecPoint() ;
Bool_t IsEmc(void) const {return kFALSE ; } // tells that this is not a EMC
Bool_t IsCPV(void) const {return kTRUE ; } // true if the recpoint is in CPV
Bool_t IsSortable() const { return kTRUE ; } // tells that this is a sortable object
- void Print(Option_t *) ;
+ void Print() ;
AliPHOSCpvRecPoint & operator = (const AliPHOSCpvRecPoint & /*rvalue*/) {
- // assignement operator requested by coding convention but not needed
- assert(0==1) ;
- return *this ;
+ Fatal("operator =", "not implemented") ; return *this ;
}
protected:
}
//____________________________________________________________________________
-void AliPHOSDigit::Print(Option_t *) const
+void AliPHOSDigit::Print() const
{
+ // Print the digit together with list of primaries
printf("PHOS digit: Amp=%d, Id=%d, Time=%f, NPrim=%d ",fAmp,fId,fTime,fNprimary);
for(Int_t index = 0; index <fNprimary; index ++ )
printf(" %d ",fPrimary[index]);
// --- ROOT system ---
-#include "TObject.h"
+//#include "TObject.h"
// --- Standard library ---
Bool_t operator==(const AliPHOSDigit &rValue) const;
AliPHOSDigit& operator+(AliPHOSDigit const &rValue) ;
AliPHOSDigit& operator*(Float_t factor) ;
-
+ AliPHOSDigit& operator=(const AliPHOSDigit) {
+ Fatal("operator = ", "not implemented") ; return *this ; }
Int_t Compare(const TObject * obj) const ;
Int_t GetNprimary() const { return fNprimary ; }
Int_t GetPrimary(Int_t index) const ;
Float_t GetTime(void) const {return fTime ;}
Bool_t IsSortable() const { return kTRUE ; }
- void Print(Option_t *) const;
+ void Print() const;
void SetAmp(Int_t Amp) { fAmp=Amp ; }
void SetTime(Float_t Time) {fTime = Time ;}
void ShiftPrimary(Int_t shift); // shift to separate different TreeK in merging
//
// --- ROOT system ---
-#include "TFile.h"
#include "TTree.h"
#include "TSystem.h"
-#include "TROOT.h"
-#include "TFolder.h"
-#include "TObjString.h"
-#include "TGeometry.h"
#include "TBenchmark.h"
+#include "TRandom.h"
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliRun.h"
-#include "AliHeader.h"
-#include "AliStream.h"
+
#include "AliRunDigitizer.h"
#include "AliPHOSDigit.h"
-#include "AliPHOS.h"
#include "AliPHOSGetter.h"
#include "AliPHOSDigitizer.h"
#include "AliPHOSSDigitizer.h"
// --- ROOT system ---
-#include "TObjString.h"
+//#include "TObjString.h"
class TArrayI ;
class TClonesArray ;
return *this ;
}
- Float_t * GetStripHalfSize() {return fStripHalfSize ;}
- Float_t GetStripWallWidthOut() {return fStripWallWidthOut ;}
+ Float_t * GetStripHalfSize() {return fStripHalfSize ;}
+ Float_t GetStripWallWidthOut() const {return fStripWallWidthOut ;}
Float_t * GetAirCellHalfSize() {return fAirCellHalfSize ;}
Float_t * GetWrappedHalfSize() {return fWrappedHalfSize ;}
Float_t GetAirGapLed() const {return fAirGapLed ;}
Float_t * GetCrystalHalfSize() {return fCrystalHalfSize ;}
Float_t * GetSupportPlateHalfSize() { return fSupportPlateHalfSize ;}
- Float_t * GetSupportPlateInHalfSize() {return fSupportPlateInHalfSize ;}
+ Float_t * GetSupportPlateInHalfSize() {return fSupportPlateInHalfSize ;}
Float_t GetSupportPlateThickness(void) const { return fSupportPlateThickness ; }
Float_t * GetPreampHalfSize() {return fPreampHalfSize ;}
Float_t * GetAPDHalfSize(void) {return fPinDiodeHalfSize ; }
Float_t * GetOuterThermoParams(void) {return fOuterThermoParams ; }
- Float_t * GetCoolerHalfSize(void) {return fCoolerHalfSize ;}
- Float_t * GetAirGapHalfSize(void) {return fAirGapHalfSize; }
+ Float_t * GetCoolerHalfSize(void) {return fCoolerHalfSize ;}
+ Float_t * GetAirGapHalfSize(void) {return fAirGapHalfSize; }
Float_t * GetInnerThermoHalfSize(void) {return fInnerThermoHalfSize ; }
Float_t * GetAlCoverParams() {return fAlCoverParams ; }
Float_t * GetFiberGlassHalfSize() {return fFiberGlassHalfSize ; }
Float_t * GetTSupport2HalfSize() {return fTSupport2HalfSize ;}
Float_t * GetTCables1HalfSize() {return fTCables1HalfSize ; }
Float_t * GetTCables2HalfSize() {return fTCables2HalfSize ; }
- Float_t GetTSupportDist() {return fTSupportDist ; }
+ Float_t GetTSupportDist() const {return fTSupportDist ; }
Float_t * GetFrameXHalfSize() {return fFrameXHalfSize ;}
Float_t * GetFrameZHalfSize() {return fFrameZHalfSize ;}
Float_t * GetFrameXPosition() {return fFrameXPosition ;}
// --- ROOT system ---
-#include "TPad.h"
#include "TH2.h"
#include "TMath.h"
#include "TCanvas.h"
// --- Standard library ---
// --- AliRoot header files ---
-
+#include "AliPHOSLoader.h"
#include "AliGenerator.h"
#include "AliPHOSGeometry.h"
+#include "AliPHOSDigit.h"
#include "AliPHOSEmcRecPoint.h"
-#include "AliRun.h"
-#include "AliPHOSGetter.h"
-
+
ClassImp(AliPHOSEmcRecPoint)
//____________________________________________________________________________
EvalDispersion(logWeight, digits) ;
EvalCoreEnergy(logWeight, digits);
EvalTime(digits) ;
- AliPHOSRecPoint::EvalAll(logWeight,digits) ;
+ AliPHOSRecPoint::EvalAll(digits) ;
}
//____________________________________________________________________________
void AliPHOSEmcRecPoint::EvalLocalPosition(Float_t logWeight, TClonesArray * digits)
// --- ROOT system ---
-//#include "TObject.h"
-#include "TArrayI.h"
-
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliPHOSDigit.h"
+//#include "AliPHOSDigit.h"
#include "AliPHOSRecPoint.h"
class AliPHOSEmcRecPoint : public AliPHOSRecPoint {
// --- AliRoot header files ---
#include "AliConfig.h"
+#include "AliPHOSDigit.h"
#include "AliPHOSClusterizer.h"
#include "AliPHOSEvalRecPoint.h"
#include "AliRun.h"
#include "AliPHOSRecCpvManager.h"
#include "AliPHOSRecEmcManager.h"
#include "AliPHOSDigitizer.h"
+#include "AliPHOSGeometry.h"
// --- Standard library ---
// fParent=parent;
TObjArray* wPool = (TObjArray*)GetWorkingPool();
if(!wPool) {
- Error("AliPHOSEvalRecPoint", "Couldn't find working pool. Exit.") ;
- exit(1);
+ Fatal("ctor", "Couldn't find working pool") ;
}
fParent = wPool->IndexOf((TObject*)parent);
TFolder* wPoolF = (TFolder*)aliceF->FindObject("WhiteBoard/RecPoints/PHOS/SmP");
AliPHOSClusterizer* clu = (AliPHOSClusterizer*)wPoolF->FindObject("PHOS:clu-v1");
if(!clu) {
- Error("GetClusterizer", "Couldn't find Clusterizer. Exit.") ;
- exit(1);
+ Fatal("GetClusterizer", "Couldn't find Clusterizer") ;
}
return clu;
// initialization
AliPHOSClusterizer* clusterizer = GetClusterizer();
if(!clusterizer) {
- Error("Init", "Cannot get clusterizer. Exit.") ;
- exit(1);
+ Fatal("Init", "Cannot get clusterizer") ;
}
TClonesArray* digits = AliPHOSLoader::GetPHOSLoader(fEventFolderName)->Digits();
newRP->AddDigit(*digit,eDigit);
}
Info("UnfoldTwoMergedPoints", "======= Split: daughter rec point %d =================", iMax) ;
- newRP->Print("");
+ newRP->Print();
}
if(TooClose(rp))
{
Info("MergeClosePoint", "+++++++ Merging point 1: ++++++") ;
- this->Print("");
+ this->Print();
Info("MergeClosePoint", "+++++++ and point 2: ++++++++++") ;
- ((AliPHOSEvalRecPoint*)rp)->Print("");
+ ((AliPHOSEvalRecPoint*)rp)->Print();
//merge two rec. points
TVector3 lpos1;
delete rp;
Info("MergeClosePoint", "++++++ Resulting point: ++++++++") ;
- this->Print("");
+ this->Print();
break;
}
AliPHOSClusterizer* clusterizer = GetClusterizer();
if(!clusterizer) {
- Error("UnfoldLocalMaxima", "Cannot get clusterizer. Exit.") ;
- exit(1);
+ Fatal("UnfoldLocalMaxima", "Cannot get clusterizer") ;
}
if(this->IsEmc()) {
newRP->EvalLocalPosition(logWeight,digits);
Info("UnfoldLocalMaxima", "======= Unfold: daughter rec point %d =================", iMax) ;
- newRP->Print("");
+ newRP->Print();
}
// RemoveFromWorkingPool(this);
}
-void AliPHOSEvalRecPoint::PrintPoint(Option_t* opt)
+void AliPHOSEvalRecPoint::PrintPoint()
{
// print rec.point to stdout
- AliPHOSCpvRecPoint::Print(opt);
+ AliPHOSCpvRecPoint::Print();
TVector3 lpos;
GetLocalPosition(lpos);
TFolder* wPoolF = (TFolder*)aliceF->FindObject("WhiteBoard/RecPoints/PHOS/SmP");
AliPHOSRecManager* recMng = (AliPHOSRecManager*)wPoolF->FindObject("AliPHOSRecManager");
if(!recMng) {
- Error("GetReconstructionManager", "Couldn't find Reconstruction Manager. Exit.") ;
- exit(1);
+ Fatal("GetReconstructionManager", "Couldn't find Reconstruction Manager") ;
}
return recMng;
TFolder* wPoolF = (TFolder*)aliceF->FindObject("WhiteBoard/RecPoints/PHOS/SmP");
TObject* wPool = wPoolF->FindObject("SmartPoints");
if(!wPool) {
- Error("GetWorkingPool", "Couldn't find Working Pool. Exit.") ;
- exit(1);
+ Fatal("GetWorkingPool", "Couldn't find Working Pool") ;
}
return wPool;
void AliPHOSEvalRecPoint::PrintWorkingPool()
{
// print pool of rec.points to stdout
- ((TObjArray*)GetWorkingPool())->Print("");
+ ((TObjArray*)GetWorkingPool())->Print();
}
AliPHOSClusterizer* GetClusterizer();
AliPHOSRecManager* GetReconstructionManager() const;
- void PrintPoint(Option_t* opt);
+ void PrintPoint();
AliPHOSRecPoint* Parent();
Float_t Chi2Dof() const;
//!
// Author: Yuri Kharlov. 17 April 2003
-#include "AliRun.h"
#include "AliGenerator.h"
#include "AliPHOSGetter.h"
#include "AliPHOSFastGlobalReconstruction.h"
// create an instance of the PHOS getter,
// create an array or reconstructed particles.
- gime = AliPHOSGetter::Instance(headerFile);
+ fgime = AliPHOSGetter::Instance(headerFile);
fGenerator = gAlice->Generator();
fParticles = new TClonesArray("TParticle",100);
fNParticles = 0;
TLorentzVector p,v;
Int_t kf,ks,imom1,imom2,idaug1,idaug2;
- gime->Event(event,"X") ;
+ fgime->Event(event,"X") ;
fParticles ->Clear();
fNParticles = 0;
- Int_t nPrimaries = gime->NPrimaries();
- TClonesArray *primaries = gime->Primaries();
+ Int_t nPrimaries = fgime->NPrimaries();
+ TClonesArray *primaries = fgime->Primaries();
for (Int_t iprim=0; iprim<nPrimaries; iprim++) {
primary = (TParticle*)primaries->At(iprim);
void SmearMomentum(TLorentzVector &p);
private:
- AliPHOSGetter *gime; //! Instance of the PHOS getter
+ AliPHOSGetter *fgime; //! Instance of the PHOS getter
AliGenerator *fGenerator; //! MC generator used in simulation
TClonesArray *fParticles; //! Array of reconstructed particles
Int_t fNParticles; //! Number of reconstructed particles
}
//____________________________________________________________________________
-void AliPHOSFastRecParticle::Print(Option_t *)const
+void AliPHOSFastRecParticle::Print()const
{
// Print the type, energy and momentum of the reconstructed particle
TString Name() const ;
virtual void Paint(Option_t * option="");
- virtual void Print(Option_t *) const ;
+ virtual void Print() const ;
void SetType(Int_t type) ;
// --- ROOT system ---
+#include "TRandom.h"
-#include "TString.h"
-#include "TFile.h"
#include "TLorentzVector.h"
#include "TList.h"
-#include "TTree.h"
#include "TParticle.h"
-#include "AliRun.h"
#include "AliPHOSGammaJet.h"
#include "AliPHOSGetter.h"
}
//____________________________________________________________________________
-void AliPHOSGammaJet::Exec(Option_t *)
+void AliPHOSGammaJet::Exec(Option_t * opt)
{
// does the job
+ if (! opt)
+ return ;
+
AliPHOSGetter * gime = AliPHOSGetter::Instance() ;
Int_t nEvent1 = gime->MaxEvent() ;
Int_t iEvent = 0 ;
}//final particle etacut
}//for (iParticle<nParticle)
TLorentzVector gamma,charge,pi0, gammapair;
- Int_t id_g = -1;
- GetGammaJet(particleList,gamma, id_g);
- GetLeadingCharge(particleList,charge, id_g);
+ Int_t idg = -1;
+ GetGammaJet(particleList,gamma, idg);
+ GetLeadingCharge(particleList,charge, idg);
GetLeadingPi0(particleList,pi0);
- Info("Pi0Decay", "Gamma: %f %d", gamma.Energy(), id_g) ;
+ Info("Pi0Decay", "Gamma: %f %d", gamma.Energy(), idg) ;
Info("Pi0Decay", "Charge: %f", charge.Energy()) ;
Info("Pi0Decay", "Pi0: %f", pi0.Energy()) ;
- // GetLeadingGammaPair(particleList, gammapair, id_g,
+ // GetLeadingGammaPair(particleList, gammapair, idg,
// 0.04,0.2, 1.0,0.13,0.14);
Info("Pi0Decay", "Pair: %f", gammapair.Energy()) ;
}//loop: events
//____________________________________________________________________________
void AliPHOSGammaJet::GetGammaJet(TList &particleList, TLorentzVector &gamma, Int_t & id)
{
+ // Get the lists of jet particles and gamma
TParticle *particle = 0x0;
Int_t iPrimary=-1, id_motherg = -1;
//____________________________________________________________________________
void AliPHOSGammaJet::GetLeadingCharge(TList &particleList, TLorentzVector &charge, Int_t & id)
{
+ // Gets the leading particle from the list of charged particles
TParticle *particle = 0x0;
Int_t iPrimary=-1;
- Double_t pt_max = 0, pt_i = 0;
+ Double_t ptmax = 0, pti = 0;
TIter next(&particleList);
while ( (particle = (TParticle*)next()) ) {
iPrimary++;
if((ksCode == 1)&&(id != iPrimary)
&&(particle->GetPDG(0)->Charge()!=0)){
- pt_i = particle->Pt();
- if(pt_i> pt_max){
- pt_max = pt_i;
+ pti = particle->Pt();
+ if(pti> ptmax){
+ ptmax = pti;
particle->Momentum(charge);
- }//pt_max
+ }//ptmax
}// kscode == 1
}// while
}
//____________________________________________________________________________
void AliPHOSGammaJet::GetLeadingPi0(TList &particleList, TLorentzVector &pi0)
{
+ // Gets the leading pi0 from the list of particles
TParticle *particle = 0x0;
Int_t iPrimary=-1;
- Double_t pt_max = 0, pt_i = 0;
+ Double_t ptmax = 0, pti = 0;
TIter next(&particleList);
while ( (particle = (TParticle*)next()) ) {
iPrimary++;
if((ksCode == 2))
{
- pt_i = particle->Pt();
- if(pt_i> pt_max){
- pt_max = pt_i;
+ pti = particle->Pt();
+ if(pti> ptmax){
+ ptmax = pti;
particle->Momentum(pi0);
- }//pt_max
+ }//ptmax
}// kscode == 1
}// while
}
AliPHOSGammaJet() ; // default ctor
AliPHOSGammaJet(const TString inputfilename) ; //ctor
AliPHOSGammaJet(const AliPHOSGammaJet & gj) ; // cpy ctor
- ~AliPHOSGammaJet() ; // dtor
- virtual void Exec(Option_t *);
+ virtual ~AliPHOSGammaJet() ; // dtor
+ virtual void Exec(Option_t * = "");
void GetGammaJet(TList & particleList, TLorentzVector & gamma, Int_t & id) ;
void GetLeadingCharge(TList & particleList, TLorentzVector & charge, Int_t & id) ;
void GetLeadingPi0(TList & particleList, TLorentzVector & pi0) ;
#include "TVector3.h"
#include "TRotation.h"
-#include "TFolder.h"
-#include "TROOT.h"
+#include "TParticle.h"
// --- Standard library ---
-#include <stdlib.h>
-
// --- AliRoot header files ---
#include "AliPHOSGeometry.h"
#include "AliPHOSEMCAGeometry.h"
#include "AliPHOSRecPoint.h"
-#include "AliConst.h"
ClassImp(AliPHOSGeometry) ;
AliPHOSGeometry * AliPHOSGeometry::fgGeom = 0 ;
Bool_t AliPHOSGeometry::fgInit = kFALSE ;
+//____________________________________________________________________________
+AliPHOSGeometry::AliPHOSGeometry() {
+ // default ctor
+ // must be kept public for root persistency purposes, but should never be called by the outside world
+ fPHOSAngle = 0 ;
+ fGeometryEMCA = 0 ;
+ fGeometrySUPP = 0 ;
+ fGeometryCPV = 0 ;
+ fgGeom = 0 ;
+ fRotMatrixArray = 0 ;
+}
+
//____________________________________________________________________________
AliPHOSGeometry::~AliPHOSGeometry(void)
{
}
}
+//____________________________________________________________________________
Bool_t AliPHOSGeometry::Impact(const TParticle * particle) const
{
- Bool_t In=kFALSE;
- Int_t ModuleNumber=0;
+ // Tells if a particle enters PHOS
+ Bool_t in=kFALSE;
+ Int_t moduleNumber=0;
Double_t z,x;
- ImpactOnEmc(particle->Theta(),particle->Phi(),ModuleNumber,z,x);
- if(ModuleNumber) In=kTRUE;
- else In=kFALSE;
- return In;
+ ImpactOnEmc(particle->Theta(),particle->Phi(),moduleNumber,z,x);
+ if(moduleNumber)
+ in=kTRUE;
+ else
+ in=kFALSE;
+ return in;
}
//____________________________________________________________________________
//
//*-- Author: Yves Schutz (SUBATECH)
-#include <assert.h>
-
// --- ROOT system ---
-#include "TString.h"
-#include "TObjArray.h"
-#include "TParticle.h"
-#include "TVector3.h"
-
// --- AliRoot header files ---
#include "AliGeometry.h"
-
#include "AliPHOSEMCAGeometry.h"
#include "AliPHOSCPVGeometry.h"
#include "AliPHOSSupportGeometry.h"
public:
- AliPHOSGeometry() {
- // default ctor
- // must be kept public for root persistency purposes, but should never be called by the outside world
- fPHOSAngle = 0 ;
- fGeometryEMCA = 0;
- fGeometrySUPP = 0;
- fGeometryCPV = 0;
- fgGeom = 0;
- fRotMatrixArray = 0; }
+ AliPHOSGeometry() ;
AliPHOSGeometry(const AliPHOSGeometry & geom) : AliGeometry(geom) {
- // cpy ctor requested by Coding Convention but not yet needed
-
- assert(0==1) ;
+ Fatal("cpy ctor", "not implemented") ;
}
virtual ~AliPHOSGeometry(void) ;
virtual Bool_t Impact(const TParticle * particle) const ;
AliPHOSGeometry & operator = (const AliPHOSGeometry & /*rvalue*/) const {
- // assignement operator requested by coding convention but not needed
- assert(0==1) ;
- return *(GetInstance()) ;
- }
+ Fatal("operator =", "nt implemented") ; return *(GetInstance()) ; }
// General
// --- ROOT system ---
#include "TSystem.h"
-#include "TFile.h"
#include "TROOT.h"
// --- Standard library ---
// --- AliRoot header files ---
-
#include "AliPHOSGetter.h"
+#include "AliPHOS.h"
#include "AliRunLoader.h"
#include "AliStack.h"
#include "AliPHOSLoader.h"
-// #include "AliPHOSRaw2Digits.h"
-//#include "AliPHOSCalibrationDB.h"
#include "AliPHOSBeamTestEvent.h"
#include "AliMC.h"
//____________________________________________________________________________
AliPHOSClusterizer * AliPHOSGetter::Clusterizer()
{
+ // Returns pointer to the Clusterizer task
AliPHOSClusterizer * rv ;
rv = dynamic_cast<AliPHOSClusterizer *>(PhosLoader()->Reconstructioner()) ;
if (!rv) {
//____________________________________________________________________________
AliPHOSDigitizer * AliPHOSGetter::Digitizer()
{
+ // Returns pointer to the Digitizer task
AliPHOSDigitizer * rv ;
rv = dynamic_cast<AliPHOSDigitizer *>(PhosLoader()->Digitizer()) ;
if (!rv) {
}
//____________________________________________________________________________
-AliPHOSTrackSegmentMaker * AliPHOSGetter::TrackSegmentMaker()
+AliPHOSTrackSegmentMaker * AliPHOSGetter::TrackSegmentMaker()
{
+ // Returns pointer to the TrackSegmentMaker task
AliPHOSTrackSegmentMaker * rv ;
rv = dynamic_cast<AliPHOSTrackSegmentMaker *>(PhosLoader()->TrackSegmentMaker()) ;
if (!rv) {
//____________________________________________________________________________
-AliPHOSPID * AliPHOSGetter::PID()
+AliPHOSPID * AliPHOSGetter::PID()
{
+ // Returns pointer to the PID task
AliPHOSPID * rv ;
rv = dynamic_cast<AliPHOSPID *>(PhosLoader()->PIDTask()) ;
if (!rv) {
}
//____________________________________________________________________________
-AliPHOSSDigitizer * AliPHOSGetter::SDigitizer()
+AliPHOSSDigitizer * AliPHOSGetter::SDigitizer()
{
+ // Returns pointer to the SDigitizer task
AliPHOSSDigitizer * rv ;
rv = dynamic_cast<AliPHOSSDigitizer *>(PhosLoader()->SDigitizer()) ;
if (!rv) {
//____________________________________________________________________________
TTree * AliPHOSGetter::TreeD() const
{
+ // Returns pointer to the Digits Tree
TTree * rv = 0 ;
rv = PhosLoader()->TreeD() ;
if ( !rv ) {
//____________________________________________________________________________
TTree * AliPHOSGetter::TreeH() const
{
+ // Returns pointer to the Hits Tree
TTree * rv = 0 ;
rv = PhosLoader()->TreeH() ;
if ( !rv ) {
//____________________________________________________________________________
TTree * AliPHOSGetter::TreeR() const
{
+ // Returns pointer to the RecPoints Tree
TTree * rv = 0 ;
rv = PhosLoader()->TreeR() ;
if ( !rv ) {
//____________________________________________________________________________
TTree * AliPHOSGetter::TreeT() const
{
+ // Returns pointer to the TrackSegments Tree
TTree * rv = 0 ;
rv = PhosLoader()->TreeT() ;
if ( !rv ) {
//____________________________________________________________________________
TTree * AliPHOSGetter::TreeP() const
{
+ // Returns pointer to the RecParticles Tree
TTree * rv = 0 ;
rv = PhosLoader()->TreeP() ;
if ( !rv ) {
//____________________________________________________________________________
TTree * AliPHOSGetter::TreeS() const
-{
+{
+ // Returns pointer to the SDigits Tree
TTree * rv = 0 ;
rv = PhosLoader()->TreeS() ;
if ( !rv ) {
// --- ROOT system ---
#include "TObject.h"
-#include "TClonesArray.h"
-// #include "TFolder.h"
-// #include "TTree.h"
-// #include "TFile.h"
-// class TString ;
- class TParticle ;
-// class TTask ;
+class TParticle ;
+class TTree ;
// --- Standard library ---
#include "AliConfig.h"
// #include "AliRun.h"
-class AliPHOS ;
-#include "AliPHOSHit.h"
-
-class AliPHOSGeometry ;
+#include "AliPHOSLoader.h"
+#include "AliPHOSHit.h"
+#include "AliPHOSDigit.h"
+#include "AliPHOSEmcRecPoint.h"
+#include "AliPHOSCpvRecPoint.h"
+#include "AliPHOSTrackSegment.h"
+#include "AliPHOSRecParticle.h"
#include "AliPHOSDigitizer.h"
-#include "AliPHOSSDigitizer.h"
-// //class AliPHOSCalibrationDB ;
-// class AliPHOSConTableDB ;
+#include "AliPHOSSDigitizer.h"
+class AliPHOS ;
+class AliPHOSGeometry ;
+class AliPHOSClusterizer ;
+class AliPHOSTrackSegmentMaker ;
+class AliPHOSPID ;
class AliPHOSBeamTestEvent ;
-#include "AliPHOSLoader.h"
-
class AliPHOSGetter : public TObject {
public:
//=========== Primaries ============
// TTree * TreeK(TString filename="") ;
- TClonesArray * Primaries(void) ;
+ TClonesArray * Primaries(void) ;
TParticle * Primary(Int_t index) const ;
Int_t NPrimaries()const { return fNPrimaries; }
TParticle * Secondary(const TParticle * p, const Int_t index=1) const ;
TTree * TreeS() const ;
AliPHOSSDigitizer * SDigitizer() ;
- TString GetSDigitsFileName() { return PhosLoader()->GetSDigitsFileName() ; }
- Int_t LoadSDigits(Option_t* opt="") { return PhosLoader()->LoadSDigits(opt) ; }
- Int_t LoadSDigitizer(Option_t* opt=""){ return PhosLoader()->LoadSDigitizer(opt) ; }
- Int_t WriteSDigits(Option_t* opt="") { return PhosLoader()->WriteSDigits(opt) ; }
- Int_t WriteSDigitizer(Option_t* opt=""){
+ TString GetSDigitsFileName() const { return PhosLoader()->GetSDigitsFileName() ; }
+ Int_t LoadSDigits(Option_t* opt="") const { return PhosLoader()->LoadSDigits(opt) ; }
+ Int_t LoadSDigitizer(Option_t* opt="") const { return PhosLoader()->LoadSDigitizer(opt) ; }
+ Int_t WriteSDigits(Option_t* opt="") const { return PhosLoader()->WriteSDigits(opt) ; }
+ Int_t WriteSDigitizer(Option_t* opt="") const {
return PhosLoader()->WriteSDigitizer(opt) ; }
//========== Digits ================
AliPHOSDigit * Digit(const Int_t index) { return static_cast<AliPHOSDigit *>(Digits()->At(index)) ;}
TTree * TreeD() const ;
AliPHOSDigitizer * Digitizer() ;
- TString GetDigitsFileName() { return PhosLoader()->GetDigitsFileName() ; }
- Int_t LoadDigits(Option_t* opt="") { return PhosLoader()->LoadDigits(opt) ; }
- Int_t LoadDigitizer(Option_t* opt=""){
+ TString GetDigitsFileName() const { return PhosLoader()->GetDigitsFileName() ; }
+ Int_t LoadDigits(Option_t* opt="") const { return PhosLoader()->LoadDigits(opt) ; }
+ Int_t LoadDigitizer(Option_t* opt="") const {
return PhosLoader()->LoadDigitizer(opt) ; }
- Int_t WriteDigits(Option_t* opt="") { return PhosLoader()->WriteDigits(opt) ; }
- Int_t WriteDigitizer(Option_t* opt=""){
+ Int_t WriteDigits(Option_t* opt="") const { return PhosLoader()->WriteDigits(opt) ; }
+ Int_t WriteDigitizer(Option_t* opt="") const {
return PhosLoader()->WriteDigitizer(opt) ; }
//========== RecPoints =============
TObjArray * CpvRecPoints() ;
AliPHOSCpvRecPoint * CpvRecPoint(const Int_t index) { return static_cast<AliPHOSCpvRecPoint *>(CpvRecPoints()->At(index)) ;}
TTree * TreeR() const ;
- AliPHOSClusterizer * Clusterizer() ;
- TString GetRecPointsFileName() { return PhosLoader()->GetRecPointsFileName() ; }
- Int_t LoadRecPoints(Option_t* opt="") { return PhosLoader()->LoadRecPoints(opt) ; }
- Int_t LoadClusterizer(Option_t* opt=""){
+ AliPHOSClusterizer * Clusterizer() ;
+ TString GetRecPointsFileName() const { return PhosLoader()->GetRecPointsFileName() ; }
+ Int_t LoadRecPoints(Option_t* opt="") const { return PhosLoader()->LoadRecPoints(opt) ; }
+ Int_t LoadClusterizer(Option_t* opt="") const {
return PhosLoader()->LoadClusterizer(opt) ; }
- Int_t WriteRecPoints(Option_t* opt="") { return PhosLoader()->WriteRecPoints(opt) ; }
- Int_t WriteClusterizer(Option_t* opt=""){
+ Int_t WriteRecPoints(Option_t* opt="") const { return PhosLoader()->WriteRecPoints(opt) ; }
+ Int_t WriteClusterizer(Option_t* opt="") const {
return PhosLoader()->WriteClusterizer(opt) ; }
//========== TrackSegments TClonesArray * TrackSegments(const char * name = 0) {
AliPHOSTrackSegment * TrackSegments(const Int_t index) { return static_cast<AliPHOSTrackSegment *>(TrackSegments()->At(index)) ;}
TTree * TreeT() const ;
AliPHOSTrackSegmentMaker * TrackSegmentMaker() ;
- TString GetTracksFileName() { return PhosLoader()->GetTracksFileName() ; }
- Int_t LoadTracks(Option_t* opt="") { return PhosLoader()->LoadTracks(opt) ; }
- Int_t LoadTrackSegementMaker(Option_t* opt=""){
+ TString GetTracksFileName() const { return PhosLoader()->GetTracksFileName() ; }
+ Int_t LoadTracks(Option_t* opt="") const { return PhosLoader()->LoadTracks(opt) ; }
+ Int_t LoadTrackSegementMaker(Option_t* opt="") const {
return PhosLoader()->LoadTrackSegmentMaker(opt) ; }
- Int_t WriteTracks(Option_t* opt="") { return PhosLoader()->WriteTracks(opt) ; }
- Int_t WriteTrackSegmentMaker(Option_t* opt=""){
+ Int_t WriteTracks(Option_t* opt="") const { return PhosLoader()->WriteTracks(opt) ; }
+ Int_t WriteTrackSegmentMaker(Option_t* opt="") const {
return PhosLoader()->WriteTracker(opt) ; }
//========== RecParticles ===========
AliPHOSRecParticle * RecPaticles(const Int_t index) { return static_cast<AliPHOSRecParticle *>(RecParticles()->At(index)) ;}
TTree * TreeP() const ;
AliPHOSPID * PID() ;
- TString GetRecParticlesFileName() { return PhosLoader()->GetRecParticlesFileName() ; }
- Int_t LoadRecParticles(Option_t* opt="") { return PhosLoader()->LoadRecParticles(opt) ; }
- Int_t LoadPID(Option_t* opt=""){
+ TString GetRecParticlesFileName() const { return PhosLoader()->GetRecParticlesFileName() ; }
+ Int_t LoadRecParticles(Option_t* opt="") const { return PhosLoader()->LoadRecParticles(opt) ; }
+ Int_t LoadPID(Option_t* opt="") const {
return PhosLoader()->LoadPID(opt) ; }
- Int_t WriteRecParticles(Option_t* opt="") { return PhosLoader()->WriteRecParticles(opt) ; }
- Int_t WritePID(Option_t* opt=""){
+ Int_t WriteRecParticles(Option_t* opt="") const { return PhosLoader()->WriteRecParticles(opt) ; }
+ Int_t WritePID(Option_t* opt="") const {
return PhosLoader()->WritePID(opt) ; }
// AliPHOSCalibrationDB * fcdb ; //!
- static AliPHOSLoader * fgPhosLoader ;
+ static AliPHOSLoader * fgPhosLoader ; // the loader for the NewIO
static AliPHOSGetter * fgObjGetter; // pointer to the unique instance of the singleton
enum EDataTypes{kHits,kSDigits,kDigits,kRecPoints,kTracks,kNDataTypes};
//*-- Author: Maxime Volkov (RRC KI) & Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC KI & SUBATECH)
// --- ROOT system ---
+#include "TVector3.h"
// --- Standard library ---
-#include <string.h>
-#include <stdlib.h>
// --- AliRoot header files ---
#include "AliPHOSHit.h"
#include "AliRun.h"
-#include "AliConst.h"
-#include "AliPHOSLoader.h"
#include "AliPHOSGeometry.h"
#include "AliPHOS.h"
//*-- Author: Maxime Volkov (RRC KI) & Yves Schutz (SUBATECH)
// --- ROOT system ---
-#include <TLorentzVector.h>
// --- AliRoot header files ---
#include "AliHit.h"
//____________________________________________________________________________
void AliPHOSImpact::Print()
{
+ // Prints particle info
printf("Impact Pid=%4d, p=(%6.2f,%6.2f,%6.2f,%6.2f) GeV, x=(%6.1f,%6.1f,%6.1f) cm\n",
fPid,fMomentum.Px(),fMomentum.Py(),fMomentum.Pz(),fMomentum.E(),
fX,fY,fZ);
// --- AliRoot header files ---
#include "AliPHOSJet.h"
+#include "AliPHOSGeometry.h"
#include "AliPHOSDigit.h"
#include "AliPHOSGetter.h"
#include "AliPHOSJetFinder.h"
// dtor
}
Int_t Compare(const TObject * obj) const;
- Int_t GetEmc(void) {
- // returns the index of EMC
- return fEmcN;
- }
- Int_t GetPpsd(void) {
- // returns the index of PPSD
- return fPpsdN ;
- }
- Float_t GetR(void) {
- // returns the distance between EMC and PPSD
- return fR ;
- }
+ Int_t GetEmc(void) const { return fEmcN; }
+ Int_t GetPpsd(void) const { return fPpsdN ; }
+ Float_t GetR(void) const { return fR ; }
Bool_t IsSortable() const{
// tells if this is a sortable object
return kTRUE ;
#include "TFile.h"
#include "TTree.h"
#include "TROOT.h"
-#include "TObjString.h"
-#include "TFolder.h"
-#include "TParticle.h"
// --- Standard library ---
-#include <Riostream.h>
// --- AliRoot header files ---
-#include "AliRun.h"
-#include "AliConfig.h"
-#include "AliDataLoader.h"
#include "AliPHOSLoader.h"
#include "AliPHOS.h"
-#include "AliPHOSDigitizer.h"
-#include "AliPHOSSDigitizer.h"
-#include "AliPHOSClusterizer.h"
-#include "AliPHOSClusterizerv1.h"
-#include "AliPHOSTrackSegmentMaker.h"
-#include "AliPHOSTrackSegmentMakerv1.h"
-#include "AliPHOSTrackSegment.h"
-#include "AliPHOSPID.h"
-#include "AliPHOSPIDv1.h"
-#include "AliPHOSGeometry.h"
+#include "AliPHOSHit.h"
#include "AliPHOSCalibrationDB.h"
ClassImp(AliPHOSLoader)
// --- ROOT system ---
#include "TClonesArray.h"
-#include "TFolder.h"
-#include "TTree.h"
class TString ;
class TParticle ;
class TTask ;
// --- Standard library ---
-#include <stdlib.h>
-#include <Riostream.h>
// --- AliRoot header files ---
#include "AliRun.h"
#include "AliLoader.h"
#include "AliRunLoader.h"
-#include "AliPHOS.h"
-#include "AliPHOSHit.h"
-#include "AliPHOSDigit.h"
-#include "AliPHOSEmcRecPoint.h"
-#include "AliPHOSCpvRecPoint.h"
-#include "AliPHOSTrackSegment.h"
-#include "AliPHOSRecParticle.h"
-class AliPHOSGeometry ;
-#include "AliPHOSDigitizer.h"
-#include "AliPHOSSDigitizer.h"
#include "AliPHOSClusterizer.h"
#include "AliPHOSTrackSegmentMaker.h"
#include "AliPHOSPID.h"
+class AliPHOS ;
+class AliPHOSGeometry ;
+class AliPHOSHit ;
+class AliPHOSDigit ;
+class AliPHOSEmcRecPoint ;
+class AliPHOSCpvRecPoint ;
+class AliPHOSTrackSegment ;
+class AliPHOSRecParticle ;
class AliPHOSCalibrationDB ;
+class AliPHOSSDigitizer ;
+class AliPHOSDigitizer ;
+
//
// --- ROOT system ---
-#include "TGeometry.h"
-#include "TDirectory.h"
-#include "TFile.h"
-#include "TTree.h"
-
+
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliRun.h"
+#include "AliConfig.h"
#include "AliPHOSPID.h"
ClassImp(AliPHOSPID)
virtual const char * Version() const { Warning("Version", "not defined" ) ; return 0 ; }
virtual void WriteRecParticles() = 0;
-private:
- virtual void Init() { Warning("Init", "not defined" ) ; }
-
protected:
TString fEventFolderName ; // event folder name
+private:
+ virtual void Init() { Warning("Init", "not defined" ) ; }
+
ClassDef(AliPHOSPID,3) // Particle Identifier algorithm (base class)
} ;
#include "AliGenerator.h"
#include "AliPHOS.h"
#include "AliPHOSPIDv0.h"
+#include "AliPHOSEmcRecPoint.h"
#include "AliPHOSClusterizerv1.h"
#include "AliPHOSTrackSegment.h"
#include "AliPHOSTrackSegmentMakerv1.h"
// --- AliRoot header files ---
#include "AliPHOSQAAlarm.h"
+#include "AliRun.h"
ClassImp(AliPHOSQAAlarm)
// --- Standard library ---
-#include <assert.h>
-
// --- AliRoot header files ---
-#include "AliRun.h"
-
class AliPHOSQAAlarm : public TObject {
public:
//////////////////////////////////////////////////////////////////////////////
// --- ROOT system ---
-#include "TROOT.h"
+#include "TDatime.h"
#include "TFolder.h"
// --- Standard library ---
#include "AliConfig.h"
#include "AliPHOSQAChecker.h"
-#include "AliPHOSQAVirtualCheckable.h"
ClassImp(AliPHOSQAChecker)
// --- Standard library ---
-#include <assert.h>
-
// --- AliRoot header files ---
#include "AliPHOSQAVirtualCheckable.h"
// --- Standard library ---
-#include <unistd.h>
+//#include <unistd.h>
// --- AliRoot header files ---
#include "AliPHOSDigit.h"
#include "AliPHOSConTableDB.h"
#include "AliPHOSBeamTestEvent.h"
#include "AliPHOSRaw2Digits.h"
-#include "AliPHOSv1.h"
-#include "../EVGEN/AliGenBox.h"
#include "AliRun.h"
ClassImp(AliPHOSRaw2Digits)
//____________________________________________________________________________
AliPHOSRaw2Digits::AliPHOSRaw2Digits(AliPHOSRaw2Digits & r2d):TTask(r2d.GetName(), r2d.GetTitle())
{
+ // cpy ctor
fInName=r2d.fInName ;
fMK1 = r2d.fMK1 ;
}
//____________________________________________________________________________
-Bool_t AliPHOSRaw2Digits::StartRootFiles(void ){
+Bool_t AliPHOSRaw2Digits::StartRootFiles(void ) const {
// //Create PHOS geometry, sets magnetic field to zero,
// //create Generator - to store target position,
// //opens out file, creates TreeE
AliPHOSRaw2Digits & operator = ( AliPHOSRaw2Digits & /*r2d*/ ) { return *this ; }
private:
- Bool_t StartRootFiles(void) ;
+ Bool_t StartRootFiles(void) const ;
Bool_t CloseRootFiles(void) ;
Bool_t ProcessRawFile() ;
void Swab4(void *from, void *to, size_t nwords) const ;
// --- AliRoot header files ---
#include "AliPHOSRecCpvManager.h"
-#include "AliPHOSGetter.h"
+#include "AliPHOSGeometry.h"
+#include "AliPHOSLoader.h"
ClassImp(AliPHOSRecCpvManager)
// --- AliRoot header files ---
#include "AliPHOSRecEmcManager.h"
-#include "AliPHOS.h"
-#include "AliRun.h"
ClassImp(AliPHOSRecEmcManager)
return ei;
}
-Float_t AliPHOSRecEmcManager::OneGamChi2(Float_t ai, Float_t ei, Float_t, Float_t& gi) const
+Float_t AliPHOSRecEmcManager::OneGamChi2(Float_t ai, Float_t ei, Float_t fi, Float_t& gi) const
{
// Chi2 used in OneGam (one-gamma fitting).
// gi is d(Chi2)/d(ai).
+ fi = 0 ;
Float_t da = ai - ei;
Float_t d = ei; // we assume that sigma(E) = sqrt(E)
gi = 2.*(ai-ei)/d;
}
-Float_t AliPHOSRecEmcManager::TwoGamChi2(Float_t ai, Float_t ei, Float_t, Float_t& gi) const
+Float_t AliPHOSRecEmcManager::TwoGamChi2(Float_t ai, Float_t ei, Float_t fi, Float_t& gi) const
{
// calculates chi^2
+ fi = 0 ;
Float_t da = ai - ei;
Float_t d = ei; // we assume that sigma(E) = sqrt(E)
gi = 2.*(ai-ei)/d;
// 6 March 2001
#include "AliPHOSRecManager.h"
-#include "AliPHOSGeometry.h"
class AliPHOSRecEmcManager : public AliPHOSRecManager {
// --- AliRoot header files ---
#include "AliPHOSRecParticle.h"
#include "AliPHOSGetter.h"
-#include "TParticle.h"
ClassImp(AliPHOSRecParticle)
//____________________________________________________________________________
const Int_t AliPHOSRecParticle::GetNPrimariesToRecParticles() const
{
-
+ // Get the number of primaries at the origine of the RecParticle
Int_t rv = 0 ;
AliPHOSGetter * gime = AliPHOSGetter::Instance() ;
Int_t emcRPindex = dynamic_cast<AliPHOSTrackSegment*>(gime->TrackSegments()->At(GetPHOSTSIndex()))->GetEmcIndex();
//____________________________________________________________________________
const TParticle * AliPHOSRecParticle::GetPrimary(Int_t index) const
{
+ // Get the list of primary particles at the origine of the RecParticle
if ( index > GetNPrimariesToRecParticles() ) {
if (fDebug)
Warning("GetPrimary", "AliPHOSRecParticle::GetPrimary -> %d is larger that the number of primaries %d",
private:
Int_t fPHOSTrackSegment ; // pointer to the associated track segment in PHOS
- Bool_t fDebug ;
+ Bool_t fDebug ; // to steer debug output
ClassDef(AliPHOSRecParticle,2) // Reconstructed Particle
};
// --- ROOT system ---
#include "TPad.h"
+#include "TGraph.h"
+#include "TPaveText.h"
#include "TClonesArray.h"
// --- Standard library ---
}
}
//____________________________________________________________________________
-void AliPHOSRecPoint::EvalAll(Float_t,TClonesArray * digits)
+void AliPHOSRecPoint::EvalAll(TClonesArray * digits)
{
//evaluates (if necessary) all RecPoint data members
// A recpoint being equivalent to a cluster in encal terminology
//*-- Author: Gines Martinez (SUBATECH)
-#include <assert.h>
-
// --- ROOT system ---
-#include "TMarker.h"
-#include "TGraph.h"
-#include "TPaveText.h"
-
+//#include "TMarker.h"
+//#include "TGraph.h"
+//#include "TPaveText.h"
+ class TClonesArray ;
// --- Standard library ---
// --- AliRoot header files ---
#include "AliRecPoint.h"
-#include "AliPHOSDigit.h"
+class AliPHOSDigit ;
class AliPHOSRecPoint : public AliRecPoint {
AliPHOSRecPoint() ; // ctor
AliPHOSRecPoint(const char * opt) ; // ctor
AliPHOSRecPoint(const AliPHOSRecPoint & rp) : AliRecPoint(rp) {
- // cpy ctor requested by Coding Convention
- // but not yet needed
- assert(0==1) ;
- }
+ Fatal("Cpy ctor", "not implemented") ; }
virtual ~AliPHOSRecPoint(){
// dtor
}
virtual void AddDigit(AliDigitNew &){
- // do not use this definition but the one below
- assert(0==1) ;
+ Fatal("AddDigit", "use AddDigit(AliPHOSDigit & digit, Float_t Energy)") ;
}
virtual void AddDigit(AliPHOSDigit & digit, Float_t Energy) = 0 ;
virtual Int_t Compare(const TObject * obj) const = 0 ;
virtual Int_t DistancetoPrimitive(Int_t px, Int_t py);
virtual void Draw(Option_t * option="") ;
virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py) ;
- virtual void EvalAll(Float_t ,TClonesArray * digits) ;
+ virtual void EvalAll(TClonesArray * digits) ;
virtual void EvalPHOSMod(AliPHOSDigit * digit) ;
virtual void EvalPrimaries(TClonesArray * digits) ;
virtual void GetGlobalPosition(TVector3 & gpos, TMatrix & gmat) const ; // return global position in ALICE
}
AliPHOSRecPoint & operator = (const AliPHOSRecPoint & ) {
- // assignement operator requested by coding convention but not needed
- assert(0==1) ;
- return *this ;
+ Fatal("operator = ", "not implemented") ; return *this ;
}
protected:
// --- ROOT system ---
-#include "TClonesArray.h"
-#include "TROOT.h"
-#include "TTree.h"
-#include "TFile.h"
-
// --- Standard library ---
+#include "Riostream.h"
// --- AliRoot header files ---
-#include "AliRun.h"
#include "AliRunLoader.h"
#include "AliPHOSReconstructioner.h"
#include "AliPHOSClusterizerv1.h"
#include "AliPHOSSDigitizer.h"
#include "AliPHOSTrackSegmentMakerv1.h"
#include "AliPHOSPIDv1.h"
-#include "AliPHOSFastRecParticle.h"
-#include "AliPHOSCpvRecPoint.h"
+
#include "AliPHOSLoader.h"
ClassImp(AliPHOSReconstructioner)
fIsInitialized = kTRUE ;
}
//____________________________________________________________________________
-void AliPHOSReconstructioner::Exec(Option_t)
+void AliPHOSReconstructioner::Exec(Option_t opt)
{
//check, if the names of branches, which should be made conicide with already
//existing
+ if (!opt)
+ return ;
if(!fIsInitialized)
Init() ;
}
//*-- Author: Gines Martinez & Yves Schutz (SUBATECH)
//*-- Dmitri Peressounko (SUBATECH & Kurchatov Institute)
-#include <stdlib.h>
+ //#include <stdlib.h>
// --- ROOT system ---
// --- ROOT system ---
-#include "TFile.h"
-#include "TROOT.h"
#include "TBenchmark.h"
+#include "TRandom.h"
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliRun.h"
+#include "AliPHOSGeometry.h"
#include "AliPHOSDigit.h"
#include "AliPHOSGetter.h"
#include "AliPHOSHit.h"
//____________________________________________________________________________
void AliPHOSSDigitizer::Unload() const
{
+ // Unloads the objects from the folder
AliPHOSGetter * gime = AliPHOSGetter::Instance() ;
AliPHOSLoader * loader = gime->PhosLoader() ;
loader->UnloadHits() ;
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliConfig.h"
class AliPHOSSDigitizer: public TTask {
// Author: Yuri Kharlov (IHEP, Protvino)
// 13 November 2000
-#include <assert.h>
-
// --- ROOT system ---
-#include "TObjArray.h"
+#include "TObject.h"
// --- AliRoot header files ---
AliPHOSSupportGeometry();
AliPHOSSupportGeometry(const AliPHOSSupportGeometry & geom) : TObject(geom) {
- // cpy ctor requested by Coding Convention but not yet needed
- assert(0==1) ;
- }
+ Fatal("cpy ctor", "not implemented") ; }
virtual ~AliPHOSSupportGeometry(void) {}
AliPHOSSupportGeometry & operator = (const AliPHOSSupportGeometry & /*rvalue*/) {
- // assignement operator requested by coding convention but not needed
- assert(0==1) ;
- return *this ;
- }
+ Fatal("operator =", "not implemented") ; return *this ; }
Float_t GetRailOuterSize(Int_t index) const { return fRailOuterSize[index]; }
Float_t GetRailPart1 (Int_t index) const { return fRailPart1[index]; }
Float_t GetRailPart2 (Int_t index) const { return fRailPart2[index]; }
// --- Standard library ---
// --- AliRoot header files ---
-
+#include "AliPHOSEmcRecPoint.h"
#include "AliPHOSTrackSegment.h"
ClassImp(AliPHOSTrackSegment)
}
//____________________________________________________________________________
-void AliPHOSTrackSegment::Print(Option_t *) const
+void AliPHOSTrackSegment::Print() const
{
// Print all information on this track Segment
// --- ROOT system ---
#include "TObject.h"
+class TClonesArray ;
// --- Standard library ---
// --- AliRoot header files ---
-
-#include "AliPHOSEmcRecPoint.h"
-#include "AliPHOSRecPoint.h"
+class AliPHOSRecPoint ;
+class AliPHOSEmcRecPoint ;
+class AliPHOSCpvRecPoint ;
class AliPHOSTrackSegment : public TObject {
Int_t GetEmcIndex() const { return fEmcRecPoint ; }
Int_t GetCpvIndex() const { return fPpsdUpRecPoint; }
- virtual void Print(Option_t *) const;
+ virtual void Print() const;
void SetIndexInList(Int_t val){ fIndexInList = val ; }
void SetCpvRecPoint(AliPHOSRecPoint * PpsdUpRecPoint ); //sets PPSD up Rec Point
// --- ROOT system ---
-#include "TGeometry.h"
-#include "TFile.h"
-#include "TTree.h"
+//#include "TGeometry.h"
+ //#include "TFile.h"
+ //#include "TTree.h"
// --- Standard library ---
// --- AliRoot header files ---
#include "AliPHOSTrackSegmentMaker.h"
-#include "AliHeader.h"
+ //#include "AliHeader.h"
ClassImp( AliPHOSTrackSegmentMaker)
class TFile ;
// --- Standard library ---
-#include <iostream>
+//#include <iostream>
// --- AliRoot header files ---
//
// --- ROOT system ---
-#include "TROOT.h"
-#include "TFile.h"
-#include "TFolder.h"
#include "TTree.h"
-#include "TSystem.h"
#include "TBenchmark.h"
// --- Standard library ---
// --- AliRoot header files ---
-
+#include "AliPHOSGeometry.h"
#include "AliPHOSTrackSegmentMakerv1.h"
-#include "AliPHOSClusterizerv1.h"
#include "AliPHOSTrackSegment.h"
-#include "AliPHOSCpvRecPoint.h"
#include "AliPHOSLink.h"
#include "AliPHOSGetter.h"
-#include "AliPHOS.h"
ClassImp( AliPHOSTrackSegmentMakerv1)
//____________________________________________________________________________
void AliPHOSTrackSegmentMakerv1::InitParameters()
{
+ //Initializes parameters
fRcpv = 10. ;
fEmcFirst = 0 ;
fEmcLast = 0 ;
//____________________________________________________________________________
void AliPHOSTrackSegmentMakerv1::Unload()
{
+ // Unloads the task from the folder
AliPHOSGetter * gime = AliPHOSGetter::Instance() ;
gime->PhosLoader()->UnloadRecPoints() ;
gime->PhosLoader()->UnloadTracks() ;
// --- ROOT system ---
-
-#include <TBRIK.h>
-#include <TNode.h>
#include <TParticle.h>
-#include <TRandom.h>
-#include <TTree.h>
#include <TVirtualMC.h>
// --- Standard library ---
-#include <string.h>
-#include <stdlib.h>
// --- AliRoot header files ---
-
-#include "AliConst.h"
#include "AliPHOSCPVDigit.h"
#include "AliPHOSGeometry.h"
#include "AliPHOSHit.h"
// --- ROOT system ---
#include <TBRIK.h>
-#include <TFile.h>
#include <TGeometry.h>
#include <TNode.h>
#include <TParticle.h>
-#include <TTree.h>
+#include "TClonesArray.h"
#include <TVirtualMC.h>
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliConst.h"
#include "AliPHOSFastRecParticle.h"
#include "AliPHOSGeometry.h"
#include "AliPHOSLoader.h"
//*-- Author: Yves Schutz (SUBATECH)
// --- ROOT system ---
-#include "TClonesArray.h"
+//#include "TClonesArray.h"
#include "TRandom.h"
class TVector3 ;
// --- ROOT system ---
-#include <TTree.h>
+//#include <TTree.h>
#include <TVirtualMC.h>
// --- Standard library ---