fLegPdg1(11),
fLegPdg2(-11),
fAlpha(0.),
- fDebug(0)
+ fDebug(0),
+ fPol(0)
{
//
// Default constructor
//
}
+
+//_____________________________________________________________________________
+AliGenPairFlat::~AliGenPairFlat()
+{
+ //
+ // Destructor
+ //
+ delete fPol;
+
+}
+
//_____________________________________________________________________________
void AliGenPairFlat::Generate()
{
TLorentzVector mother, dau1, dau2;
Float_t random[6];
- TF1* pol = new TF1("Pol","1.+[0]*x*x",-1.,1.);
- pol->SetParameter(0, fAlpha);
+ fPol = new TF1("fPol","1.+[0]*x*x",-1.,1.);
+ fPol->SetParameter(0, fAlpha);
for (j=0;j<3;j++) origin[j]=fOrigin[j];
time = fTime;
}
+ if(fDebug == 2){
printf("\n\n------------------GENERATOR SETTINGS------------------\n\n");
printf("You choosed for the mother the Mass range %f - %f \n",fPairMassMin,fPairMassMax);
printf("You choosed for the mother the transverse Momentum range %f - %f \n",fPairPtMin,fPairPtMax);
printf("polarization factor is alpha == %lf \n",fAlpha);
printf("vertex is at x == %f || y == %f || z == %f \n",origin[0],origin[1],origin[2]);
printf("\n----------------------------------------------------------\n");
-
+ }
for(i=0;i<fPairNpart;i++) {
if (fDebug == 2) printf("p = (%+11.4e,%+11.4e,%+11.4e) GeV || pt = %+11.4e || y = %+11.4e || weight=%+11.4e || E= %+11.4e\n",p[0],p[1],p[2],pt, y,weight,e);
//decay procedure
- if(!Decay(mother,dau1,dau2,pol))continue;
+ if(!Decay(mother,dau1,dau2,fPol))continue;
pt1 = dau1.Pt();
pt2 = dau2.Pt();
}
//_____________________________________________________________________________
-
void AliGenPairFlat::Init()
{
//
//_____________________________________________________________________________
-
-
-Bool_t AliGenPairFlat::Decay(TLorentzVector mother, TLorentzVector &dau1, TLorentzVector &dau2, TF1* polfactor)
+Bool_t AliGenPairFlat::Decay(const TLorentzVector& mother, TLorentzVector &dau1, TLorentzVector &dau2, TF1* polfactor)
{
//
// decay procedure