fQRange(0.06),
fQNBins(60),
fGenerator(0x0),
+ fStackBuffer(0),
fBufferSize(5),
fNBinsToScale(Int_t(fQNBins*0.1)),
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(0),
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()
}
fGenerator->Generate();
- Int_t j = 0, ntr = 0;
+ Int_t ntr = 0;
if ( genstack->GetNtrack() < fNpart/2)
{
Warning("Generate","************************************************************");
qlong = gRandom->Uniform(qlongc-binwdh, qlongc+binwdh);
TParticle* first = 0;
- while (j < genstack->GetNtrack())
+ Int_t jj = 0;
+
+ while (jj < genstack->GetNtrack())
{
- TParticle* tmpp = genstack->Particle(j++);
+ TParticle* tmpp = genstack->Particle(jj++);
if (tmpp->GetPdgCode() == fPID)
{
if (CheckParticle(tmpp,0x0,stack) == kFALSE)
for (Int_t k = 1; k<=fQNBins; k++)
{
- for (Int_t j = 1; j<=fQNBins; j++)
+ for (j = 1; j<=fQNBins; j++)
{
for (Int_t i = 1; i<=fQNBins; i++)
{
}
Double_t beta = pzsum/esum;
- Double_t gamma = 1.0/TMath::Sqrt(1.0 - beta*beta);
+ Double_t gamma = 1.0/TMath::Sqrt((1.-beta)*(1.+beta));
lon = gamma * ( pzdiff - beta*ediff );
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");
+}
+
+