diff -uNr v01-19-01.ORIG/MarlinReco/v01-17/Analysis/RecoMCTruthLink/src/RecoMCTruthLinker.cc v01-19-01/MarlinReco/v01-17/Analysis/RecoMCTruthLink/src/RecoMCTruthLinker.cc --- v01-19-01.ORIG/MarlinReco/v01-17/Analysis/RecoMCTruthLink/src/RecoMCTruthLinker.cc 2017-01-26 23:00:13.000000000 +0900 +++ v01-19-01/MarlinReco/v01-17/Analysis/RecoMCTruthLink/src/RecoMCTruthLinker.cc 2017-01-28 17:34:32.000000000 +0900 @@ -951,9 +951,9 @@ for ( unsigned kkk=0 ; kkk < mother->getDaughters().size() ; kkk++ ) { MCParticle* sister = dynamic_cast(mother->getDaughters()[kkk]); if ( sister == this_Kid ) continue; - if ( abs(sister->getVertex()[0]-this_Kid->getVertex()[0]) > 0.1 || - abs(sister->getVertex()[1]-this_Kid->getVertex()[1]) > 0.1 || - abs(sister->getVertex()[2]-this_Kid->getVertex()[2]) > 0.1 ) continue; // must check that it is the same vertex: + if ( fabs(sister->getVertex()[0]-this_Kid->getVertex()[0]) > 0.1 || + fabs(sister->getVertex()[1]-this_Kid->getVertex()[1]) > 0.1 || + fabs(sister->getVertex()[2]-this_Kid->getVertex()[2]) > 0.1 ) continue; // must check that it is the same vertex: // several "non-destructive interactions" can // take place (think delta-rays !) if ( sister->isBackscatter()) { diff -uNr v01-19-01.ORIG/MarlinReco/v01-17/Analysis/TrueJet/src/TrueJet.cc v01-19-01/MarlinReco/v01-17/Analysis/TrueJet/src/TrueJet.cc --- v01-19-01.ORIG/MarlinReco/v01-17/Analysis/TrueJet/src/TrueJet.cc 2017-01-26 22:58:31.000000000 +0900 +++ v01-19-01/MarlinReco/v01-17/Analysis/TrueJet/src/TrueJet.cc 2017-01-28 17:37:35.000000000 +0900 @@ -493,9 +493,9 @@ str_tmom[i]=0. ; str_tmomS[i]=0. ; str_mom[i]=0. ; } if ( masssq_t > 0. ) { - if ( abs( sqrt(masssq_t)- p[k[elementon[ijet]][4]][5] )/ p[k[elementon[ijet]][4]][5] > 0.001 ) { + if ( fabs( sqrt(masssq_t)- p[k[elementon[ijet]][4]][5] )/ p[k[elementon[ijet]][4]][5] > 0.001 ) { if ( type[ijet]%100 == 1 && type[ijet-1]%100 == 1 && nfsr[ijet]+nfsr[ijet-1]==0) { - if ( abs( tE[ijet]+tE[ijet-1]- p[k[elementon[ijet]][4]][4] )/ p[k[elementon[ijet]][4]][4] > 0.001 ) { + if ( fabs( tE[ijet]+tE[ijet-1]- p[k[elementon[ijet]][4]][4] )/ p[k[elementon[ijet]][4]][4] > 0.001 ) { streamlog_out(ERROR) << " bad match M (sum/initial) " << " " << sqrt(masssq_t) << " / " << p[k[elementon[ijet]][4]][5] << std::endl; streamlog_out(ERROR) << " E (sum/initial) " << " " << tE[ijet]+tE[ijet-1] << " / " << p[k[elementon[ijet]][4]][4] << std::endl; @@ -537,7 +537,7 @@ for ( int jj= k[kk][4] ; jj <= k[kk][5] ; jj++ ) { e_kid+=p[jj][4]; } - if ( (abs(e_kid - p[kk][4])/ p[kk][4] ) > 0.001 ) { + if ( (fabs(e_kid - p[kk][4])/ p[kk][4] ) > 0.001 ) { streamlog_out(WARNING) << kk << " " << k[kk][4] << " " << k[kk][5] << " " << e_kid << " " << p[kk][4] << std::endl ; } } @@ -1089,7 +1089,7 @@ int parent1=k[line94][3] ; int parent2=k[line94][6] ; int kid1=k[line94][4] ; int kid2=k[line94][5] ; if ( parent1 != 0 && parent2 != 0 ) { - if ( abs((p[kid1][4]+p[kid2][4])-(p[parent1][4]+p[parent2][4]))/(p[parent1][4]+p[parent2][4]) > 0.0001 ) { + if ( fabs((p[kid1][4]+p[kid2][4])-(p[parent1][4]+p[parent2][4]))/(p[parent1][4]+p[parent2][4]) > 0.0001 ) { streamlog_out(WARNING) << " Inconsitent 94 object: Energy of parents = "<< p[kid1][4]+p[kid2][4] << " , of kids = " << p[parent1][4]+p[parent2][4] << std::endl; streamlog_out(WARNING) << " parents: " << parent1 << " " << parent2 << ". Kids: " << kid1 << " " << kid2 << std::endl; streamlog_out(WARNING) << " pdgs: " << k[parent1][2] << " " << k[parent2][2] << ". Kids: " << k[kid1][2] << " " << k[kid2][2] << std::endl; @@ -1919,7 +1919,7 @@ for ( int jj=k[i][4] ; jj <= k[i][5] ; jj++ ) { Ekid+= p[jj][4]; } - if ( abs(Ekid-p[i][4])/p[i][4] > 0.001 ) { + if ( fabs(Ekid-p[i][4])/p[i][4] > 0.001 ) { streamlog_out(WARNING) << " Particle " << i << " has energy " << p[i][4] << " , but the sum of its genstat1 kids is " << Ekid << std::endl; streamlog_out(WARNING) << " indicating that Geant did something fishy and un-documented in MCParticle " @@ -2219,7 +2219,7 @@ k[line][2] == k[k[Candidate_94s[ipi]][5]][2] ) || (k[mothers][2] == k[k[Candidate_94s[ipi]][5]][2] && k[line][2] == k[k[Candidate_94s[ipi]][4]][2] )) && - ( abs(Ekids- Eparents)/Ekids < 0.0001 ) ) { + ( fabs(Ekids- Eparents)/Ekids < 0.0001 ) ) { // Right flavour, right mass : say no more diff -uNr v01-19-01.ORIG/MarlinReco/v01-17/Analysis/TrueJet_Parser/src/TrueJet_Parser.cc v01-19-01/MarlinReco/v01-17/Analysis/TrueJet_Parser/src/TrueJet_Parser.cc --- v01-19-01.ORIG/MarlinReco/v01-17/Analysis/TrueJet_Parser/src/TrueJet_Parser.cc 2017-01-26 23:00:01.000000000 +0900 +++ v01-19-01/MarlinReco/v01-17/Analysis/TrueJet_Parser/src/TrueJet_Parser.cc 2017-01-28 17:38:40.000000000 +0900 @@ -52,7 +52,7 @@ for ( unsigned kk=0 ; kk(mcpvec[kk]); if ( _COUNT_FSR ) { - if ( abs(www[kk]) == 1.0) { + if ( fabs(www[kk]) == 1.0) { E+=mcp->getEnergy(); } } else { @@ -81,7 +81,7 @@ for ( unsigned kk=0 ; kk(mcpvec[kk]); if ( _COUNT_FSR ) { - if ( abs(www[kk]) == 1.0) { + if ( fabs(www[kk]) == 1.0) { const double* mom = mcp->getMomentum(); for (int kk=0 ; kk<3 ; kk++ ) { p3[kk]+=mom[kk];