X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PWGLF%2FSTRANGENESS%2FHypernuclei%2FAliAODMCNuclExReplicator.cxx;h=59126274022d23ee92d689dfa0f4051bcd228644;hb=e8f0a920a93fe1a5eeede1802653686d806d87a4;hp=d1951cd3fe6ae6758121acc1b6f55bf7e13333b4;hpb=abed0cc7329c0b7819fcd92fc87e2c1af933468f;p=u%2Fmrichter%2FAliRoot.git

diff --git a/PWGLF/STRANGENESS/Hypernuclei/AliAODMCNuclExReplicator.cxx b/PWGLF/STRANGENESS/Hypernuclei/AliAODMCNuclExReplicator.cxx
index d1951cd3fe6..59126274022 100644
--- a/PWGLF/STRANGENESS/Hypernuclei/AliAODMCNuclExReplicator.cxx
+++ b/PWGLF/STRANGENESS/Hypernuclei/AliAODMCNuclExReplicator.cxx
@@ -284,12 +284,15 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
   TClonesArray *arrayMC = 0;
   AliAODMCHeader *mcHeader=0;
   Int_t mumpdg=-100;
- 
+  Int_t mumpdgNeg=-100;
+  
+
   arrayMC = (TClonesArray*) source.GetList()->FindObject(AliAODMCParticle::StdBranchName());
   if (!arrayMC) {
     Printf("Error: MC particles branch not found!\n");
     return;
   }
+  
   // if(arrayMC)
   //   cout<<"Ho caricato array mc"<<endl;
   
@@ -298,9 +301,10 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
     printf("AliAnalysisTaskSEDs::UserExec: MC header branch not found!\n");
     return;
   } 
+
   // if(mcHeader)
   //   cout<<"Ho caricato MC header"<<endl;
-
+  
   //  cout<<"Centrality AOD source: "<<source.GetHeader()->GetCentrality()<<endl;
 
   Int_t nsv(0);
@@ -316,12 +320,12 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
   Double_t dispersion;
 
   // cout<<"Qui"<<endl;
-  //cout<<source.GetMagneticField()<<endl;
+  // cout<<source.GetMagneticField()<<endl;
 
   AliAODVertex *vtx = source.GetPrimaryVertex();
 						
-  //  cout<<"Source "<<source<<endl;
-  //cout<<"vtx: "<<vtx<<endl;
+  // cout<<"Source "<<source<<endl;
+  // cout<<"vtx: "<<vtx<<endl;
 
   // A Set of very loose cut for a weak strange decay
   
@@ -342,7 +346,7 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
   vtx->GetXYZ(pos);
   vtx->GetCovarianceMatrix(cov);
   fV1 = new AliESDVertex(pos,cov,100.,100,vtx->GetName());
-  //  cout<<"fV1 pricipal loop: "<<fV1<<endl;
+  cout<<"fV1 pricipal loop: "<<fV1<<endl;
   
   if(entries<=0) return;
   indices = new UShort_t[entries];
@@ -424,13 +428,19 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
 
     if(mumpdg == -1010010030){
 
-      //	if(PDGCode==-211 || PDGCode==+211){  
-      if(PDGCode==+211){  
+      if(PDGCode==-211 || PDGCode==+211){  
+	//      if(PDGCode==-211){  
+	
+	//	if(aodtrack->Charge()>0)continue;
+	
 	Track0[nTrack0++]=j;
       }
       
-      //	if(PDGCode==1000020030 ||PDGCode==-1000020030 ){
-      if(PDGCode==-1000020030){
+      if(PDGCode==1000020030 ||PDGCode==-1000020030 ){
+      //if(PDGCode==1000020030){
+	
+	//	if(aodtrack->Charge()<0)continue;
+	
 	Track1[nTrack1++]=j;
 	// 	new((*fNuclei)[nnuclei++]) AliAODTrack(*aodtrack);
       }
@@ -470,6 +480,8 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
     AliAODMCParticle *partPos = (AliAODMCParticle*) arrayMC->At(labelpos);
     Int_t mumidPos = partPos->GetMother();
 
+    //    cout<<"\n\nmumidPos: "<<mumidPos <<endl;
+
     //------------------------------
     
     for (Int_t k=0; k <nTrack0 ; k++) {                           //! Pion Tracks Loop
@@ -486,11 +498,16 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
       AliAODMCParticle *partNeg = (AliAODMCParticle*) arrayMC->At(labelneg);
       Int_t mumidNeg = partNeg->GetMother();
       
-      
+      //      cout<<"mumidNeg: "<<mumidNeg <<endl;
+
+      // AliAODMCParticle *motherNeg = (AliAODMCParticle*) arrayMC->At(mumidNeg);
+      // mumpdgNeg = motherNeg->GetPdgCode();
+
       //------------------------------
       //  if(mumidPos == mumidNeg && mumidNeg > 0){
       isOk=kFALSE;
       
+      // if(mumidPos == mumidNeg && mumidNeg > 0 && mumpdgNeg  == 1010010030)
       if(mumidPos == mumidNeg && mumidNeg > 0)
 	isOk = kTRUE;
       
@@ -525,6 +542,7 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
 	
       if(io2Prong->CosPointingAngle()<fCosMin)continue;
 	
+      //      cout<<"CPA: "<<io2Prong->CosPointingAngle()<<endl;
       
       // AliAODTrack *trk0 = (AliAODTrack*)io2Prong->GetDaughter(0);
       // AliAODTrack *trk1 = (AliAODTrack*)io2Prong->GetDaughter(1);
@@ -537,6 +555,8 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source)
       // cout<<"**********************************************"<<endl;
 	
       //      rd =  new((*fSecondaryVerices)[nsv++]) AliAODRecoDecayLF2Prong(*io2Prong);
+
+      //      cout<<"is ok??? "<<isOk<<endl;
       if(isOk){
 	new((*fSecondaryVerices)[nsv++]) AliAODRecoDecayLF2Prong(*io2Prong);
 	new((*fDaughterTracks)[ntracksD++]) AliAODTrack(*negtrackAOD);