Remove unnecessary print from QA class
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 18 Jan 2011 16:14:10 +0000 (16:14 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 18 Jan 2011 16:14:10 +0000 (16:14 +0000)
AliMCAnalysisUtils: Add to the method to get first common ancestor of 2 particles the possibility to recover the momentum of this ancestor.

PWG4/PartCorrBase/AliMCAnalysisUtils.cxx
PWG4/PartCorrBase/AliMCAnalysisUtils.h
PWG4/PartCorrDep/AliAnaCalorimeterQA.cxx

index d63538d..cfa43e3 100755 (executable)
@@ -107,7 +107,7 @@ Int_t AliMCAnalysisUtils::CheckOrigin(const Int_t * label, const Int_t nlabels,
 
 //_________________________________________________________________________
 Int_t AliMCAnalysisUtils::CheckCommonAncestor(const Int_t index1, const Int_t index2, AliCaloTrackReader* reader, 
-                                              Int_t & ancPDG, Int_t & ancStatus) {
+                                              Int_t & ancPDG, Int_t & ancStatus, TLorentzVector & momentum) {
   //Check the first common ancestor of 2 clusters, given the most likely labels of the primaries generating such clusters.
   Int_t label1[100];
   Int_t label2[100];
@@ -188,6 +188,7 @@ Int_t AliMCAnalysisUtils::CheckCommonAncestor(const Int_t index1, const Int_t in
           if (mom) {
             ancPDG    = mom->GetPdgCode();
             ancStatus = mom->GetStatus();
+            momentum.SetPxPyPzE(mom->Px(),mom->Py(),mom->Pz(),mom->E());
           }
         }
         else {
@@ -195,6 +196,7 @@ Int_t AliMCAnalysisUtils::CheckCommonAncestor(const Int_t index1, const Int_t in
           if (mom) {
             ancPDG    = mom->GetPdgCode();
             ancStatus = mom->GetStatusCode();
+            mom->Momentum(momentum);
           }
         }
         //First ancestor found, end the loops
index 00494e9..958a7ee 100755 (executable)
@@ -45,7 +45,7 @@ class AliMCAnalysisUtils : public TObject {
        
        //Check only the label of the most significant particle
   Int_t   CheckCommonAncestor(const Int_t index1, const Int_t index2, AliCaloTrackReader* reader, 
-                              Int_t & ancPDG, Int_t & ancStatus) ;
+                              Int_t & ancPDG, Int_t & ancStatus, TLorentzVector & momentum) ;
        Int_t   CheckOrigin(const Int_t label, AliCaloTrackReader * reader, const Int_t input) ;
        //Check the label of the most significant particle but do checks on the rest of the contributing labels
        Int_t   CheckOrigin(const Int_t *label, const Int_t nlabels, AliCaloTrackReader * reader, const Int_t input) ;
index 63adeed..1875d71 100755 (executable)
@@ -1526,7 +1526,6 @@ void  AliAnaCalorimeterQA::MakeAnalysisFillHistograms()
           if(imax == ipos) continue;
           absId  = indexList[ipos]; 
           Float_t diff = (tmax-cell->GetCellTime(absId))*1e9;
-          printf("diff time %f %p",diff,fhCellTimeSpreadRespectToCellMax);
           fhCellTimeSpreadRespectToCellMax->Fill(diff);
           if(TMath::Abs(TMath::Abs(diff) > 100)) fhCellIdCellLargeTimeSpread->Fill(absId);
         }// fill cell-cluster histogram loop