message: Removing dummy charge; variables not needed for monopoles moved to else part
diff: @@ -2411,9 +2411,6 @@ void AliTPCv2::StepManager()
  if(TMath::Abs(charge)<=0. && pdg!=kPdgMonopole) return; // take only charged particles (except magnetic monopoles)
  - if(pdg==kPdgMonopole && TMath::Abs(charge)<=0.) // need to take care of monopoles with 0 charge (safety for later used formulae)
  - charge = 0.000001;
  -
  // check the sensitive volume
  id = mc->CurrentVolID(copy); // vol ID and copy number (starts from 1!)
  @@ -2516,6 +2513,16 @@ void AliTPCv2::StepManager()
  // Calculate number of ionization electrons nel
  Int_t nel=0;
  + // parameters needed for magnetic monopoles as well
  + const Double_t Emin = fTPCParam->GetFpot();
  + const Double_t wIon = fTPCParam->GetWmean();
  +
  + // special treatmeant of magnetic monopoles
  + if(pdg==kPdgMonopole){
  + nel = (Int_t)((mc->Edep()-Emin)/wIon) + 1;
  + }
  + else{
  +
  // Stepsize in cm
  Double_t stepSize = mc->TrackStep();
  // printf("Stepsize %f\n", stepSize);
  @@ -2538,17 +2545,10 @@ void AliTPCv2::StepManager()
  // Variables needed to generate random powerlaw distributed energy loss
  const Double_t alpha = -1 * fTPCParam->GetExp(); // NA49/G3 value
  const Double_t alpha_p1 = alpha+1;
  - const Double_t Emin = fTPCParam->GetFpot();
  const Double_t Emax = fTPCParam->GetEend();
  const Double_t Kmin = TMath::Power(Emin, alpha_p1);
  const Double_t Kmax = TMath::Power(Emax, alpha_p1);
  - const Double_t wIon = fTPCParam->GetWmean();
  - // special treatmeant of magnetic monopoles
  - if(pdg==kPdgMonopole){
  - nel = (Int_t)((mc->Edep()-Emin)/wIon) + 1;
  - }
  - else{
  for(Int_t n = 0; n < nColl; n++) {
  // Use GEANT3 / NA49 expression:
  // P(Edep) ~ k * Edep^alpha
  