Included check for errors in the stdout

[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliAnalysisTaskSEBkgLikeSignJPSI.cxx

diff --git a/PWG3/vertexingHF/AliAnalysisTaskSEBkgLikeSignJPSI.cxx b/PWG3/vertexingHF/AliAnalysisTaskSEBkgLikeSignJPSI.cxx
index 19402db5e7e4b261f6a4f08e1bcb84d0a48a6db1..ca8a816fc3c1e16c9634d4cf37b85ef5dd959403 100644 (file)
--- a/PWG3/vertexingHF/AliAnalysisTaskSEBkgLikeSignJPSI.cxx
+++ b/PWG3/vertexingHF/AliAnalysisTaskSEBkgLikeSignJPSI.cxx
@@ -241,7 +241,7 @@ void AliAnalysisTaskSEBkgLikeSignJPSI::UserExec(Option_t */*option*/)
   // loop over Like sign candidates
   Int_t nPosPairs=0,nNegPairs=0;
   Int_t nLikeSign = arrayLikeSign->GetEntriesFast();
-  printf("+++\n+++Number of like sign pairs ---> %d \n+++\n", nLikeSign);
+  if(fDebug>1) printf("+++\n+++Number of like sign pairs ---> %d \n+++\n", nLikeSign);
 
   for(Int_t iLikeSign = 0; iLikeSign < nLikeSign; iLikeSign++) {
     AliAODRecoDecayHF2Prong *d = (AliAODRecoDecayHF2Prong*)arrayLikeSign->UncheckedAt(iLikeSign);
@@ -251,42 +251,42 @@ void AliAnalysisTaskSEBkgLikeSignJPSI::UserExec(Option_t */*option*/)
         unsetvtx=kTRUE;
     }
     Int_t okBtoJPSIls=0;
-    Float_t ddddd=d->CosPointingAngle();
     if(d->SelectBtoJPSI(fVHF->GetBtoJPSICuts(),okBtoJPSIls)) {
+       AliAODTrack *trk0 = (AliAODTrack*)d->GetDaughter(0);
        fHistMassLS->Fill(d->InvMassJPSIee());
        fHistCPtaLS->Fill(d->CosPointingAngle());
        fHistd0d0LS->Fill(1e8*d->Prodd0d0());
-       fHistCtsLS->Fill(d->CosThetaStarJPSI());
        fHistDCALS->Fill(100*d->GetDCA());
-       PostData(1,fOutput);
-       AliAODTrack *trk0 = (AliAODTrack*)d->GetDaughter(0);
+       //PostData(1,fOutput);
        if(!trk0) {
           trk0=aod->GetTrack(trkIDtoEntry[d->GetProngID(0)]);
           printf("references to standard AOD not available \n");
        }
        if((trk0->Charge())==1) {
           nPosPairs++;
-          fHistCtsLSpos->Fill(d->CosThetaStarJPSI());
-          PostData(1,fOutput);
+          fHistCtsLS->Fill(d->CosThetaStar(0,443,11,11));
+          fHistCtsLSpos->Fill(d->CosThetaStar(0,443,11,11));
+          //PostData(1,fOutput);
         } else {
           nNegPairs++;
+          fHistCtsLS->Fill(d->CosThetaStarJPSI());
           fHistCtsLSneg->Fill(d->CosThetaStarJPSI());
-          PostData(1,fOutput);
+          //PostData(1,fOutput);
         }
        PostData(1,fOutput);
     }
     if(unsetvtx) d->UnsetOwnPrimaryVtx();
   }
 
-  printf("------------ N. of positive pairs in Event ----- %d \n", nPosPairs);
-  printf("------------ N. of negative pairs in Event ----- %d \n", nNegPairs);
+  if(fDebug>1) printf("------------ N. of positive pairs in Event ----- %d \n", nPosPairs);
+  if(fDebug>1) printf("------------ N. of negative pairs in Event ----- %d \n", nNegPairs);
 
   fTotPosPairs += nPosPairs;
   fTotNegPairs += nNegPairs;
 
   // loop over JPSI candidates
   Int_t nBtoJpsiToEle = arrayJPSItoEle->GetEntriesFast();
-  printf("Number of like JPSI -> ee candidates ---> %d \n", nBtoJpsiToEle);
+  if(fDebug>1) printf("Number of like JPSI -> ee candidates ---> %d \n", nBtoJpsiToEle);
 
   for (Int_t iBtoJpsiToEle = 0; iBtoJpsiToEle < nBtoJpsiToEle; iBtoJpsiToEle++) {
     AliAODRecoDecayHF2Prong *d = (AliAODRecoDecayHF2Prong*)arrayJPSItoEle->UncheckedAt(iBtoJpsiToEle);
@@ -338,13 +338,15 @@ void AliAnalysisTaskSEBkgLikeSignJPSI::Terminate(Option_t */*option*/)
   fHistDCAJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistDCAJPSI"));
   fHistDCALS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistDCALS"));
 
-  fHistMassLS->Scale((1/fLsNormalization)*fHistMassLS->GetEntries());
-  fHistCtsLS->Scale((1/fLsNormalization)*fHistCtsLS->GetEntries());
-  fHistCtsLSpos->Scale((1/fLsNormalization)*fHistCtsLSpos->GetEntries());
-  fHistCtsLSneg->Scale((1/fLsNormalization)*fHistCtsLSneg->GetEntries());
-  fHistCPtaLS->Scale((1/fLsNormalization)*fHistCPtaLS->GetEntries());
-  fHistd0d0LS->Scale((1/fLsNormalization)*fHistd0d0LS->GetEntries());
-  fHistDCALS->Scale((1/fLsNormalization)*fHistDCALS->GetEntries());
+  if(fLsNormalization>0.) {
+    fHistMassLS->Scale((1/fLsNormalization)*fHistMassLS->GetEntries());
+    fHistCtsLS->Scale((1/fLsNormalization)*fHistCtsLS->GetEntries());
+    fHistCtsLSpos->Scale((1/fLsNormalization)*fHistCtsLSpos->GetEntries());
+    fHistCtsLSneg->Scale((1/fLsNormalization)*fHistCtsLSneg->GetEntries());
+    fHistCPtaLS->Scale((1/fLsNormalization)*fHistCPtaLS->GetEntries());
+    fHistd0d0LS->Scale((1/fLsNormalization)*fHistd0d0LS->GetEntries());
+    fHistDCALS->Scale((1/fLsNormalization)*fHistDCALS->GetEntries());
+  }
 
   return;
 }