]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
Bug fix (Sergey)
authordainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 9 Sep 2010 15:20:47 +0000 (15:20 +0000)
committerdainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 9 Sep 2010 15:20:47 +0000 (15:20 +0000)
PWG3/vertexingHF/AliAnalysisTaskSESignificance.cxx

index 9cb606905942bfeb2eca97437ac852dbb55de509..e45daf52ebf94029b298c41934a9facf6d707bc9 100644 (file)
@@ -450,6 +450,8 @@ void AliAnalysisTaskSESignificance::UserExec(Option_t */*option*/)
     prongpdg=4122;
     break;
   }
+  Int_t prongPdgPlus=prongpdg;
+  Int_t prongPdgMinus=TMath::Abs(prongPdgPlus)*(-1);
 
   Int_t nHistpermv=((AliMultiDimVector*)fCutList->FindObject("multiDimVectorPtBin0"))->GetNTotCells();
   Int_t nProng = arrayProng->GetEntriesFast();
@@ -490,12 +492,12 @@ void AliAnalysisTaskSESignificance::UserExec(Option_t */*option*/)
        fHistNEvents->Fill(3);
        if(fReadMC){
          Int_t lab=-1;
-         lab = d->MatchToMC(prongpdg,arrayMC,nprongs,pdgdaughters);
+         lab = d->MatchToMC(prongPdgPlus,arrayMC,nprongs,pdgdaughters);
          if(lab>=0){ //signal
            AliAODMCParticle *dMC = (AliAODMCParticle*)arrayMC->At(lab);
            Int_t pdgMC = dMC->GetPdgCode();
            
-           if(pdgMC==+prongpdg) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMass);
+           if(pdgMC==prongPdgPlus) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMass);
            else fRflHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMass);
          }
          else{ //background
@@ -524,11 +526,11 @@ void AliAnalysisTaskSESignificance::UserExec(Option_t */*option*/)
          fMassHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
          if(fReadMC){
            Int_t lab=-1;
-           lab = d->MatchToMC(-prongpdg,arrayMC,nprongs,pdgdaughters);
+           lab = d->MatchToMC(prongPdgMinus,arrayMC,nprongs,pdgdaughters);
            if(lab>=0){ //signal
              AliAODMCParticle *dMC = (AliAODMCParticle*)arrayMC->At(lab);
              Int_t pdgMC = dMC->GetPdgCode();
-             if(pdgMC==-prongpdg) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
+             if(pdgMC==prongPdgMinus) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
              else fRflHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
            }
            else{ //background
@@ -550,12 +552,12 @@ void AliAnalysisTaskSESignificance::UserExec(Option_t */*option*/)
         fMassHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
         if(fReadMC){
          Int_t lab=-1;
-         lab = d->MatchToMC(-prongpdg,arrayMC,nprongs,pdgdaughters);
+         lab = d->MatchToMC(prongPdgMinus,arrayMC,nprongs,pdgdaughters);
                                                                 //cambia il match
         if(lab>=0){ //signal
         AliAODMCParticle *dMC = (AliAODMCParticle*)arrayMC->At(lab);
         Int_t pdgMC = dMC->GetPdgCode();
-        if(pdgMC==-prongpdg) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
+        if(pdgMC==prongPdgMinus) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
         else fRflHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
        }
                                                                                         else{ //background