+/**************************************************************************
+ * 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$ */
+
#include "AliGenHBTosl.h"
+#include "AliLog.h"
+
//__________________________________________________________
/////////////////////////////////////////////////////////////
// //
fQRange(0.06),
fQNBins(60),
fGenerator(0x0),
+ fStackBuffer(0),
fBufferSize(5),
fNBinsToScale(Int_t(fQNBins*0.1)),
- fDebug(kFALSE),
+ fDebug(0),
+ fSignalShapeCreated(0),
fMaxIterations(1000),
fMaxChiSquereChange(0.01),
fMaxChiSquerePerNDF(1.5),
fSamplePhiMax(TMath::TwoPi()+0.01),
fSignalRegion(0.0),
fMinFill(1000),
+ fSwapped(0),
fLogFile(0x0)
{
//default constructor
fQRange(0.06),
fQNBins(60),
fGenerator(0x0),
+ fStackBuffer(0),
fBufferSize(5),
fNBinsToScale(Int_t(fQNBins*0.1)),
- fDebug(kFALSE),
+ fDebug(0),
fSignalShapeCreated(kFALSE),
fMaxIterations(1000),
fMaxChiSquereChange(0.01),
fSamplePhiMax(TMath::TwoPi()+0.01),
fSignalRegion(0.0),
fMinFill(1000),
+ fSwapped(0),
fLogFile(0x0)
{
//default constructor
}
+
+AliGenHBTosl::AliGenHBTosl(const AliGenHBTosl & hbt):
+ AliGenerator(-1),
+ fQCoarseBackground(0x0),
+ fQCoarseSignal(0x0),
+ fQSignal(0x0),
+ fQBackground(0x0),
+ fQSecondSignal(0x0),
+ fQSecondBackground(0x0),
+ fQRange(0.06),
+ fQNBins(60),
+ fGenerator(0x0),
+ fStackBuffer(0),
+ fBufferSize(5),
+ fNBinsToScale(Int_t(fQNBins*0.1)),
+ fDebug(0),
+ fSignalShapeCreated(kFALSE),
+ fMaxIterations(1000),
+ fMaxChiSquereChange(0.01),
+ fMaxChiSquerePerNDF(1.5),
+ fQRadius(8.0),
+ fPID(kPiPlus),
+ fSamplePhiMin(-0.01),
+ fSamplePhiMax(TMath::TwoPi()+0.01),
+ fSignalRegion(0.0),
+ fMinFill(1000),
+ fSwapped(0),
+ fLogFile(0x0)
+{
+// Copy constructor
+ hbt.Copy(*this);
+}
/***********************************************************/
AliGenHBTosl::~AliGenHBTosl()
// TParticle(pdg, is, parent, -1, kFirstDaughter, kLastDaughter,
// px, py, pz, e, vx, vy, vz, tof);
- if (GetDebug()) Info("GetOneD","Randomized qinv = %f, obtained = %f",qinv,GetQInv(first,second));
+ AliDebug(1,Form("Randomized qinv = %f, obtained = %f",qinv,GetQInv(first,second)));
}
/***********************************************************/
second->SetMomentum(vector.X(),vector.Y(),vector.Z(),e);
////////////
- if ( GetDebug() > 3 )
+ if ( AliDebugLevel() > 3 )
{
e=TMath::Sqrt(m2+px*px);
TParticle* f = new TParticle(first->GetPdgCode(),0,-1,-1,-1,-1, px , 0.0, 0.0, e,0.0,0.0,0.0,0.0);
}
}
- //if (GetDebug())
+ //if (AliDebugLevel())
TestCoarseSignal();
Info("FillCoarseSignal","DONE");
//Calculates the factor that should be used to scale
//quatience of num and den to 1 at tail
- if (GetDebug()) Info("Scale","Entered Scale()");
+ AliDebug(1,"Entered");
if(!num)
{
- Error("Scale","No numerator");
+ AliError("No numerator");
return 0.0;
}
if(!den)
{
- Error("Scale","No denominator");
+ AliError("No denominator");
return 0.0;
}
if(fNBinsToScale < 1)
{
+ AliError("Number of bins for scaling is smaller than 1");
return 0.0;
- Error("Scale","Number of bins for scaling is smaller thnan 1");
}
Int_t fNBinsToScaleX = fNBinsToScale;
Int_t fNBinsToScaleY = fNBinsToScale;
Int_t nbinsX = num->GetNbinsX();
if (fNBinsToScaleX > nbinsX)
{
- Error("Scale","Number of X bins for scaling is bigger thnan number of bins in histograms");
+ AliError("Number of X bins for scaling is bigger thnan number of bins in histograms");
return 0.0;
}
Int_t nbinsY = num->GetNbinsX();
if (fNBinsToScaleY > nbinsY)
{
- Error("Scale","Number of Y bins for scaling is bigger thnan number of bins in histograms");
+ AliError("Number of Y bins for scaling is bigger thnan number of bins in histograms");
return 0.0;
}
Int_t nbinsZ = num->GetNbinsZ();
if (fNBinsToScaleZ > nbinsZ)
{
- Error("Scale","Number of Z bins for scaling is bigger thnan number of bins in histograms");
+ AliError("Number of Z bins for scaling is bigger thnan number of bins in histograms");
return 0.0;
}
- if (GetDebug()>0) Info("Scale","No errors detected");
+ AliDebug(1,"No errors detected");
Int_t offsetX = nbinsX - fNBinsToScaleX - 1; //bin that we start loop over bins in axis X
Int_t offsetY = nbinsY - fNBinsToScaleY - 1; //bin that we start loop over bins in axis Y
}
}
- if(GetDebug() > 0)
- Info("Scale","numsum=%f densum=%f fNBinsToScaleX=%d fNBinsToScaleY=%d fNBinsToScaleZ=%d",
- numsum,densum,fNBinsToScaleX,fNBinsToScaleY,fNBinsToScaleZ);
+ AliDebug(1,Form("numsum=%f densum=%f fNBinsToScaleX=%d fNBinsToScaleY=%d fNBinsToScaleZ=%d",
+ numsum,densum,fNBinsToScaleX,fNBinsToScaleY,fNBinsToScaleZ));
if (numsum == 0) return 0.0;
Double_t ret = densum/numsum;
- if(GetDebug() > 0) Info("Scale","returning %f",ret);
+ AliDebug(1,Form("returning %f",ret));
return ret;
}
//In other words: To make equations easier
static TVector3 first;
- if (GetDebug())
+ if (AliDebugLevel()>=1)
{
first.SetXYZ(relvector.x(),relvector.y(),relvector.z());
}
vector.RotateY(rotAngleY);
vector.RotateZ(rotAngleZ);
- if (GetDebug()>5)
+ if (AliDebugLevel()>5)
{
TVector3 test(firstPx,0.0,0.0);
test.RotateY(rotAngleY);
test.RotateZ(rotAngleZ);
- ::Info("Rotate","Rotation test: px %f %f",first.x(),test.x());
- ::Info("Rotate","Rotation test: py %f %f",first.y(),test.y());
- ::Info("Rotate","Rotation test: pz %f %f",first.z(),test.z());
+ AliInfo(Form("Rotation test: px %f %f",first.x(),test.x()));
+ AliInfo(Form("Rotation test: py %f %f",first.y(),test.y()));
+ AliInfo(Form("Rotation test: pz %f %f",first.z(),test.z()));
}
}
/***********************************************************/
TVector3 v(x,y,z);
Double_t a1 = -TMath::ATan2(v.Y(),v.X());
- if (GetDebug()>5)
+ if (AliDebugLevel()>5)
{
- ::Info("Rotate","Xpr = %f Ypr = %f",xprime,yprime);
- ::Info("Rotate","Calc sin = %f, and %f",sinphi,TMath::Sin(a1));
- ::Info("Rotate","Calc cos = %f, and %f",cosphi,TMath::Cos(a1));
+ AliInfo(Form("Xpr = %f Ypr = %f",xprime,yprime));
+ AliInfo(Form("Calc sin = %f, and %f",sinphi,TMath::Sin(a1)));
+ AliInfo(Form("Calc cos = %f, and %f",cosphi,TMath::Cos(a1)));
}
Double_t xprimezlength = TMath::Hypot(xprime,z);
Double_t xbis = sintheta*z + costheta*(cosphi*x - sinphi*y);
- ::Info("Rotate","Calculated rot %f, modulus %f",xbis,TMath::Sqrt(x*x+y*y+z*z));
+ AliInfo(Form("Calculated rot %f, modulus %f",xbis,TMath::Sqrt(x*x+y*y+z*z)));
return xbis;
}
/***********************************************************/
{
//Checks if a given particle is falling into signal region with any other particle
//already existing on stack
- return kFALSE;
+ //PH return kFALSE;
if (fSignalRegion <=0) return kFALSE;
fSwapped = kTRUE;
}
+
+AliGenHBTosl& AliGenHBTosl::operator=(const AliGenHBTosl& rhs)
+{
+// Assignment operator
+ rhs.Copy(*this);
+ return *this;
+}
+
+void AliGenHBTosl::Copy(TObject&) const
+{
+ //
+ // Copy
+ //
+ Fatal("Copy","Not implemented!\n");
+}
+
+
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff