coverity
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 26 Mar 2011 13:04:10 +0000 (13:04 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 26 Mar 2011 13:04:10 +0000 (13:04 +0000)
PWG4/PartCorrDep/AliAnaShowerParameter.cxx

index 98dba8f..f8c01d8 100644 (file)
@@ -394,7 +394,7 @@ void  AliAnaShowerParameter::MakeAnalysisFillHistograms()
   \r
   //Do analysis and fill histograms\r
   if ( GetDebug() > 0) printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() - Starting analysis.\n");\r
-    \r
+  \r
   // Access MC information in stack if requested, check that it exists.        \r
   AliStack * stack = 0x0;\r
   //TParticle * primary = 0x0;   \r
@@ -402,8 +402,8 @@ void  AliAnaShowerParameter::MakeAnalysisFillHistograms()
   //AliAODMCParticle * aodprimary = 0x0; \r
   //TObjArray * pl = 0x0;\r
   Int_t NClusters = 0 ;\r
-    TLorentzVector momCluster ;\r
-\r
+  TLorentzVector momCluster ;\r
+  \r
   \r
   //Check if the stack is available when analysing MC data.\r
   if(IsDataMC()){\r
@@ -426,17 +426,17 @@ void  AliAnaShowerParameter::MakeAnalysisFillHistograms()
     }// is data and MC\r
   }    \r
   //Loop on stored AOD photons\r
-\r
+  \r
   TClonesArray *  clustArray = 0x0 ; \r
   clustArray = GetAODCaloClusters() ;\r
   NClusters = clustArray->GetEntriesFast() ;    \r
   fhNClusters->Fill(NClusters) ;\r
-\r
+  \r
   if(GetDebug() > 0) printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() - aod branch entries %d\n", NClusters);\r
   \r
   for(Int_t iCluster = 0 ; iCluster < NClusters ; iCluster++){    \r
     Int_t input = 0 ;\r
-\r
+    \r
     AliAODCaloCluster * caloCluster =  (AliAODCaloCluster*) (clustArray->At(iCluster)) ;\r
     caloCluster->GetMomentum(momCluster,GetVertex(0)) ;\r
     AliAODPWG4Particle * aodCluster = new AliAODPWG4Particle(momCluster) ;\r
@@ -446,7 +446,7 @@ void  AliAnaShowerParameter::MakeAnalysisFillHistograms()
     aodCluster->SetCaloLabel(caloCluster->GetID(),-1) ;\r
     aodCluster->SetDetector(fCalorimeter) ;\r
     aodCluster->SetTag(GetMCAnalysisUtils()->CheckOrigin(caloCluster->GetLabels(),caloCluster->GetNLabels(),GetReader(), aodCluster->GetInputFileIndex())) ;\r
-\r
+    \r
     if (GetDebug() > 3) printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() - Cluster particle label = %d\n",labelCluster) ;\r
     \r
     //Fill Cluster histograms \r
@@ -463,49 +463,49 @@ void  AliAnaShowerParameter::MakeAnalysisFillHistograms()
       fhEtaPhiPtCluster->Fill(ptcluster,phicluster,etacluster) ;\r
       fhLambdaPtCluster->Fill(ptcluster,lambdaMainCluster) ;\r
       \r
-    //Play with the MC data if available\r
+      //Play with the MC data if available\r
       if(IsDataMC()){\r
-       \r
-       if(GetReader()->ReadStack() && !stack) return;\r
-       if(GetReader()->ReadAODMCParticles() && !mcparticles0) return;\r
-       \r
-       //Get the tag from AliMCAnalysisUtils for PID\r
-       Int_t tag = aodCluster->GetTag();\r
-       if (GetDebug() > 3) printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() - Cluster particle tag= %d\n",tag) ;\r
-       if ( aodCluster->GetLabel() < 0){\r
-         if(GetDebug() > 0) \r
-           printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() - Label is -1; problem in the MC ESD? ");\r
-         continue;\r
-       }\r
-       \r
-       //Check if the tag matches one of the different particle types and fill the corresponding histograms\r
-       if( GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPhoton)&&!(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPi0))) {\r
-         fhLambdaPtPhoton->Fill(ptcluster,lambdaMainCluster) ;\r
-       }//kMCPhoton\r
-       \r
-       if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPi0)) {     \r
-         fhLambdaPtPi0->Fill(ptcluster,lambdaMainCluster) ;        \r
-       }\r
-       if( GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPion)) {\r
-         fhLambdaPtPion->Fill(ptcluster,lambdaMainCluster) ;\r
-       }\r
-       \r
-       // Access MC information in stack if requested, check that it exists.\r
-       Int_t label = aodCluster->GetLabel();\r
-       if(label < 0) {\r
-         printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() *** bad label ***:  label %d \n", label);\r
-         continue;\r
-       }\r
+        \r
+        //if(GetReader()->ReadStack() && !stack) return;\r
+        //if(GetReader()->ReadAODMCParticles() && !mcparticles0) return;\r
+        \r
+        //Get the tag from AliMCAnalysisUtils for PID\r
+        Int_t tag = aodCluster->GetTag();\r
+        if (GetDebug() > 3) printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() - Cluster particle tag= %d\n",tag) ;\r
+        if ( aodCluster->GetLabel() < 0){\r
+          if(GetDebug() > 0) \r
+            printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() - Label is -1; problem in the MC ESD? ");\r
+          continue;\r
+        }\r
+        \r
+        //Check if the tag matches one of the different particle types and fill the corresponding histograms\r
+        if( GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPhoton)&&!(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPi0))) {\r
+          fhLambdaPtPhoton->Fill(ptcluster,lambdaMainCluster) ;\r
+        }//kMCPhoton\r
+        \r
+        if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPi0)) {     \r
+          fhLambdaPtPi0->Fill(ptcluster,lambdaMainCluster) ;        \r
+        }\r
+        if( GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPion)) {\r
+          fhLambdaPtPion->Fill(ptcluster,lambdaMainCluster) ;\r
+        }\r
+        \r
+        // Access MC information in stack if requested, check that it exists.\r
+        Int_t label = aodCluster->GetLabel();\r
+        if(label < 0) {\r
+          printf("AliAnaShowerParameter::MakeAnalysisFillHistograms() *** bad label ***:  label %d \n", label);\r
+          continue;\r
+        }\r
       }\r
     }\r
   }\r
   \r
   Float_t fPtGener = 0 ;\r
-  if(IsDataMC()){\r
+  if(IsDataMC() && stack){\r
     TParticle* pGener = stack->Particle(0) ;\r
     fPtGener = pGener->Pt() ;\r
   }\r
-\r
+  \r
   \r
 }\r
 \r