Allowing to reject events with 0 tracks

[u/mrichter/AliRoot.git] / PWGCF / Correlations / DPhi / AliAnalysisTaskMinijet.cxx

diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskMinijet.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskMinijet.cxx
index 2335e73bf839f2e4d842e0d721afca72ca7ff716..44f7562b9beadc6b48474663cbc5e65993da5769 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliAnalysisTaskMinijet.cxx
+++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskMinijet.cxx
@@ -105,7 +105,7 @@ AliAnalysisTaskMinijet::AliAnalysisTaskMinijet(const char *name)
     fChargedPi0(0),
     fVertexCheck(0),
     fPropagateDca(0),
-    fCentralityMethod("V0M")
+    fCentralityMethod("")
 {
     
     //Constructor
@@ -171,8 +171,21 @@ void AliAnalysisTaskMinijet::UserCreateOutputObjects()
     // Called once
     if(fDebug) Printf("In User Create Output Objects.");
     
-    Int_t nbinsCentr = 101;
-    Float_t minbinCentr=-0.5, maxbinCentr=100.5;
+    Int_t nbinsCentr = 0;
+    Float_t minbinCentr=0, maxbinCentr=0;
+
+    if (fCentralityMethod.Length() > 0)
+    {
+        nbinsCentr = 100;
+        minbinCentr=0;
+        maxbinCentr=100;
+    }
+    else
+    {
+        nbinsCentr = 101;
+        minbinCentr=-0.5;
+        maxbinCentr=100.5;
+    }
 
     fStep = new TH1F("fStep", "fStep", 10, -0.5, 9.5);
     fEventStat = new TH1F("fEventStat", "fEventStat", 10, -0.5, 9.5);
@@ -628,7 +641,7 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *)
                 // read tracks
                 if(fESDEvent)     ntracks = ReadEventESD(pt, eta, phi, charge,strangenessWeight, nTracksTracklets, 5);
                 else if(fAODEvent)ntracks = ReadEventAOD(pt, eta, phi, charge,strangenessWeight, nTracksTracklets, 5);
-                else Printf("Fatal Error");
+                else AliInfo("Fatal Error");
                 
                 if (fCentralityMethod.Length() > 0)
                     ntracks = TMath::Nint(centrality);
@@ -645,7 +658,10 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *)
                     // read tracks
                     if(fESDEvent)     ntracks = ReadEventESD(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 7);//stagness version not yet implemented
                     else if(fAODEvent)ntracks = ReadEventAOD(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 7);
-                    else Printf("Fatal Error");
+                    else AliInfo("Fatal Error");
+                    
+                    if (fCentralityMethod.Length() > 0)
+                        ntracks = TMath::Nint(centrality);
                     
                     // analyse
                     if(pt.size()){ //(internally ntracks=fNRecAccept)
@@ -659,7 +675,10 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *)
                     // read tracks
                     if(fESDEvent)       ntracks = ReadEventESDRecMcProp(pt, eta, phi, charge,strangenessWeight, nTracksTracklets, 4);
                     else if(fAODEvent)  ntracks = ReadEventAODRecMcProp(pt, eta, phi, charge,strangenessWeight, nTracksTracklets, 4);
-                    else Printf("Fatal Error");
+                    else AliInfo("Fatal Error");
+                    
+                    if (fCentralityMethod.Length() > 0)
+                        ntracks = TMath::Nint(centrality);
                     
                     //analyse
                     if(pt.size()){//(internally ntracks=fNRecAccept)
@@ -673,7 +692,10 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *)
                         // read tracks
                         if(fESDEvent)       ntracks = ReadEventESDRecMcProp(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 6);//stagness version not yet implemented
                         else if(fAODEvent)  ntracks = ReadEventAODRecMcProp(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 6);
-                        else Printf("Fatal Error");
+                        else AliInfo("Fatal Error");
+                        
+                        if (fCentralityMethod.Length() > 0)
+                            ntracks = TMath::Nint(centrality);
                         
                         //analyse
                         if(pt.size()){//(internally ntracks=fNRecAccept)
@@ -685,7 +707,12 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *)
                     // read tracks
                     if(fESDEvent)       ntracks = ReadEventESDMC(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 3);
                     else if(fAODEvent)  ntracks = ReadEventAODMC(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 3);
-                    else Printf("Fatal Error");
+                    else AliInfo("Fatal Error");
+                    
+                    if (fCentralityMethod.Length() > 0){
+                        fNRecAccept = TMath::Nint(centrality);
+                        fNMcPrimAccept = TMath::Nint(centrality);
+                    }
                     
                     // analyse
                     if(pt.size()){
@@ -709,7 +736,11 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *)
             if(CheckEvent(false)){// all events, with and without reconstucted vertex
                 if(fESDEvent) ntracks       = ReadEventESDMC(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 1);//read tracks
                 else if(fAODEvent) ntracks  = ReadEventAODMC(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 1);//read tracks
-                else Printf("Fatal Error");
+                else AliInfo("Fatal Error");
+                
+                if (fCentralityMethod.Length() > 0)
+                    fNMcPrimAccept = TMath::Nint(centrality);
+                
                 
                 // analyse
                 if(pt.size()){
@@ -730,7 +761,10 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *)
                 //read tracks
                 if(fESDEvent)     ntracks  = ReadEventESDMC(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 0);
                 else if(fAODEvent) ntracks  = ReadEventAODMC(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 0);
-                else Printf("Fatal Error");
+                else AliInfo("Fatal Error");
+                
+                if (fCentralityMethod.Length() > 0)
+                    fNMcPrimAccept = TMath::Nint(centrality);
                 
                 //analyse
                 if(pt.size()){
@@ -745,6 +779,9 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *)
         if(fMode==0)       ntracks  = ReadEventESDMC(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 0);
         else if (fMode==1) ntracks  = ReadEventAODMC(pt, eta, phi, charge, strangenessWeight, nTracksTracklets, 0);
         
+        if (fCentralityMethod.Length() > 0)
+            fNMcPrimAccept = TMath::Nint(centrality);
+        
         // analyse
         if(pt.size()){
             Analyse(pt, eta, phi, charge, strangenessWeight, fNMcPrimAccept, nTracksTracklets[1],nTracksTracklets[2], 0); 
@@ -1306,7 +1343,7 @@ Int_t AliAnalysisTaskMinijet::ReadEventAODRecMcProp( vector<Float_t> &ptArray,
     TClonesArray *mcArray = (TClonesArray*)fAODEvent->
     FindListObject(AliAODMCParticle::StdBranchName());
     if(!mcArray){
-        Printf("No MC particle branch found");
+        AliInfo("No MC particle branch found");
         return kFALSE;
     }
     
@@ -1440,7 +1477,7 @@ Int_t AliAnalysisTaskMinijet::ReadEventAODMC( vector<Float_t> &ptArray,  vector<
     TClonesArray *mcArray = (TClonesArray*)fAODEvent->
     FindListObject(AliAODMCParticle::StdBranchName());
     if(!mcArray){
-        Printf("No MC particle branch found");
+        AliInfo("No MC particle branch found");
         return kFALSE;
     }
     
@@ -1740,16 +1777,16 @@ void AliAnalysisTaskMinijet::Analyse(const vector<Float_t> &pt,
                 fDPhiEventAxis[step]->Fill(dPhi, strangeWeightEventAxis*strangeWeightOthers);
                 
                 //check outliers
-                if(ptEventAxis< 0.4 || ptEventAxis > 100) Printf("particles out of range pt");
-                if(ntracksCharged<0 || ntracksCharged>150) Printf("particles out of range ncharge");
+                if(ptEventAxis< 0.4 || ptEventAxis > 100) AliInfo("particles out of range pt");
+                if(ntracksCharged<-1 || ntracksCharged>1500) AliInfo("particles out of range ncharge");
                 if(TMath::Abs(dEta)>2*fEtaCut) {
-                    Printf("particles out of range dEta");
-                    Printf("eta1=%f, eta2=%f", etaOthers, etaEventAxis);
-                    Printf("step=%d",step);
+                    AliInfo("particles out of range dEta");
+                    AliInfo(Form("eta1=%f, eta2=%f", etaOthers, etaEventAxis));
+                    AliInfo(Form("step=%d",step));
                 }
                 if(dPhi<-0.5*TMath::Pi() || dPhi>1.5*TMath::Pi()){
-                    Printf("particles out of range dPhi");
-                    Printf("phi1=%f, phi2=%f", phiOthers, phiEventAxis);
+                    AliInfo(Form("particles out of range dPhi"));
+                    AliInfo(Form("phi1=%f, phi2=%f", phiOthers, phiEventAxis));
                 }
                 
                 Bool_t isLikeSign = CheckLikeSign(chargeEventAxis, chargeOthers);
@@ -1846,7 +1883,7 @@ Bool_t AliAnalysisTaskMinijet::SelectParticlePlusCharged(const Short_t charge, c
     }
     
     else{
-        Printf("Error: wrong selection of charged/pi0/k0");
+        AliInfo("Error: wrong selection of charged/pi0/k0");
         return 0;
     }
     
@@ -1955,7 +1992,7 @@ Bool_t AliAnalysisTaskMinijet::CheckEvent(const Bool_t recVertex)
             TClonesArray *mcArray = (TClonesArray*)fAODEvent->
             FindListObject(AliAODMCParticle::StdBranchName());
             if(!mcArray){
-                Printf("No MC particle branch found");
+                AliInfo("No MC particle branch found");
                 return false;
             }
             
@@ -2090,11 +2127,11 @@ Bool_t AliAnalysisTaskMinijet::CheckLikeSign(const Short_t chargeEventAxis,
     }
     
     else{
-        Printf("Error: Charge not lower nor higher as zero");
+        AliInfo("Error: Charge not lower nor higher as zero");
         return false;
     }
     
-    Printf("Error: Check values of Charge");
+    AliInfo("Error: Check values of Charge");
     return false;
 }