-
- ////// propagation to the end of TPC //////////////
- Double_t xk=77.415;
- track->PropagateTo(xk, 28.94, 1.204e-3); //Ne
- xk -=0.01;
- track->PropagateTo(xk, 44.77, 1.71); //Tedlar
- xk -=0.04;
- track->PropagateTo(xk, 44.86, 1.45); //kevlar
- xk -=2.0;
- track->PropagateTo(xk, 41.28, 0.029); //Nomex
- xk-=16;
- track->PropagateTo(xk,36.2,1.98e-3); //C02
- xk -=0.01;
- track->PropagateTo(xk, 24.01, 2.7); //Al
- xk -=0.01;
- track->PropagateTo(xk, 44.77, 1.71); //Tedlar
- xk -=0.04;
- track->PropagateTo(xk, 44.86, 1.45); //kevlar
- xk -=0.5;
- track->PropagateTo(xk, 41.28, 0.029); //Nomex
- ////////////////////////////////////////////////////////////////////
-
- AliITSTrackV1 trackITS(*track);
+
+ /// mass definition ////////////////////////
+ Double_t mass=0.13956995;
+ Int_t pcode=211; // a pion by default
+
+ if(realmass) {
+ if(TMath::Abs(pcode)<20443) mass=db->GetParticle(pcode)->Mass();
+ }
+ else {
+ mass = track->GetMass();
+// cout << "Mass = " << mass << endl;
+ }
+
+
+
+ // new propagation to the end of TPC
+ Double_t xk=80.;
+ // track->PropagateTo(xk,0.,0.); //Ne if it's still there //attenzione funziona solo se modifica in TPC
+ // Double_t xk=77.415;
+ track->PropagateTo(xk, 28.94, 1.204e-3);
+ xk-=0.005;
+ track->PropagateTo(xk, 44.77,1.71); //Tedlar
+ xk-=0.02;
+ track->PropagateTo(xk, 44.86, 1.45); //Kevlar
+ xk-=2.0;
+ track->PropagateTo(xk, 41.28, 0.029);//Nomex
+ xk-=0.02;
+ track->PropagateTo(xk, 44.86, 1.45); //Kevlar
+ xk-=0.005;
+ track->PropagateTo(xk, 44.77, 1.71); //Tedlar
+
+ xk=61.;
+ // track->PropagateTo(xk,0.,0.); //C02
+ track->PropagateTo(xk,36.2,1.98e-3); //C02 //attenzione funziona solo se modifica in TPC
+
+ xk -=0.005;
+ track->PropagateTo(xk, 24.01, 2.7); //Al
+ xk -=0.005;
+ track->PropagateTo(xk, 44.77, 1.71); //Tedlar
+ xk -=0.02;
+ track->PropagateTo(xk, 44.86, 1.45); //Kevlar
+ xk -=0.5;
+ track->PropagateTo(xk, 41.28, 0.029); //Nomex
+ xk -=0.02;
+ track->PropagateTo(xk, 44.86, 1.45); //Kevlar
+ xk -=0.005;
+ track->PropagateTo(xk, 44.77, 1.71); //Tedlar
+ xk -=0.005;
+ track->PropagateTo(xk, 24.01, 2.7); //Al
+
+ ////////////////////////////////////////////////////////////////////////////////////////////////////////
+ //AliITSTrackV1 trackITS(*track);
+ AliITSTrackV1 trackITS(*track, fFieldFactor);
+ //cout<<" fFieldFactor = "<<fFieldFactor<<"\n";
+ trackITS.PutMass(mass); //new to add mass to track