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 190fd5cf4086190612c2bbca4f9be2122ee14c16..ca8a816fc3c1e16c9634d4cf37b85ef5dd959403 100644 (file)
--- a/PWG3/vertexingHF/AliAnalysisTaskSEBkgLikeSignJPSI.cxx
+++ b/PWG3/vertexingHF/AliAnalysisTaskSEBkgLikeSignJPSI.cxx
@@ -24,11 +24,15 @@
 #include <TSystem.h>
 #include <TROOT.h>
 #include <TClonesArray.h>
+#include <TNtuple.h>
+#include <TList.h>
+#include <TH1F.h>
 
 #include "AliAODEvent.h"
 #include "AliAODVertex.h"
 #include "AliAODTrack.h"
 #include "AliAODRecoDecayHF2Prong.h"
+#include "AliAnalysisVertexingHF.h"
 #include "AliAnalysisTaskSE.h"
 #include "AliAnalysisTaskSEBkgLikeSignJPSI.h"
 
@@ -116,8 +120,6 @@ void AliAnalysisTaskSEBkgLikeSignJPSI::Init()
   fVHF->PrintStatus();
 
   return;
-
-  return;
 }
 
 //________________________________________________________________________
@@ -239,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);
@@ -250,40 +252,41 @@ void AliAnalysisTaskSEBkgLikeSignJPSI::UserExec(Option_t */*option*/)
     }
     Int_t okBtoJPSIls=0;
     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);
-      }
+       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);
@@ -301,7 +304,7 @@ void AliAnalysisTaskSEBkgLikeSignJPSI::UserExec(Option_t */*option*/)
       fHistDCAJPSI->Fill(100*d->GetDCA());
       PostData(1,fOutput);
     }
-   if(unsetvtx) d->UnsetOwnPrimaryVtx();
+    if(unsetvtx) d->UnsetOwnPrimaryVtx();
   }
 
   return;
@@ -335,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;
 }