Commit eaab12fe authored by Breno Rilho Lemos's avatar Breno Rilho Lemos 💬
Browse files

Replace AliAnalysisTaskMyTask.cxx

parent 918a691a
......@@ -132,6 +132,7 @@ void AliAnalysisTaskMyTask::UserExec(Option_t *)
Double_t Vz = fESD->GetPrimaryVertex()->GetZ();
Double_t MagneticField = fESD->GetMagneticField();
// Add VERTEZES (x, y, z) and magnetic field to esd-summary.dat file
summary << iTracks << " " << Vx << " " << Vy << " " << Vz << " " << MagneticField << endl;
/*
......@@ -148,16 +149,22 @@ Assumed Units: Mass (GeV/c^2)[CONFIRMED] || Energy (GeV) || Momentum (GeV/c) ||
Double_t Mass = track->M(); // returns the pion mass, if the particle can't be identified properly
Double_t Energy = track->E(); // Returns the energy of the particle given its assumed mass, but assumes the pion mass if the particle can't be identified properly.
Double_t Px = track->Px();
Double_t Py = track->Py();
Double_t Pt = track->Pt();
Double_t PtCal = sqrt(Px*Px+Py*Py); // here, we calculate transversal momentum to make sure of its meaning (then compare it to Pt)
Double_t Pt = track->Pt(); // transversal momentum, in case we need it
Double_t Pz = track->Pz();
Double_t Vel = 1 / sqrt( Mass*Mass / (Pt*Pt + Pz*Pz) + 1 );
Double_t Vx = Px * sqrt( 1 - Vel*Vel ) / Mass;
Double_t Vy = Py * sqrt( 1 - Vel*Vel ) / Mass;
Double_t Vz = Pz * sqrt( 1 - Vel*Vel ) / Mass;
Double_t Charge = track->Charge();
// Add MASS, CHARGE, MOMENTUM (x, y, z) to esd-detail.dat file
// Add MASS, CHARGE and VELOCITIES (x, y, z) to esd-detail.dat file
detail << fixed << Mass << " " << Charge << " ";
detail << fixed << Px << " " << Py << " " << Pz << endl;
detail << fixed << Px << " " << Py << " " << Pz << " " << Vx << " " << Vy << " " << Vz << endl;
//'fixed' fixes the number of decimal places so numbers are vertically aligned
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment