// *******************
// Gamma conversion
// *******************
+ /*
if(step==kReconstruction)
if( (TMath::Sqrt(lPtArmV0*lPtArmV0 + lAlphaV0*lAlphaV0) < 0.2) && isNaturalPart ){
fAssocParticles->Add( new AliMiniParticle(centrality, zv, iV0, pt, lPhi, lEta, lMCAssocNegDaug, lMCAssocPosDaug, 2) );
}
+ */
// *******************
// K0s selection
fTriggerParticles->RemoveAt(iArray);
fTriggerParticles->AddAt( new AliMiniParticle(centrality, zv, idTrig, ptTrig, phiTrig, etaTrig, 0, 0, 0), iArray);
+
} // Close isTrigFromV0daug
}// End K0s Mass cut
fK0sMassPtRap[curCentBin]->Fill(massK0s,pt,v0->RapK0Short());
fK0sMassPtPhi->Fill(massK0s,pt,lPhi);
+ /*
if( (pt>kPtBinV0[0]) && (pt<kPtBinV0[kN1]) && isNaturalPart )
fAssocParticles->Add( new AliMiniParticle(centrality, zv, iV0, pt, lPhi, lEta, lMCAssocNegDaug, lMCAssocPosDaug, 3) );
-
+ */
+
}
if( fDoQA && lCheckMcK0Short && isNaturalPart && (pt<10.) ){ // Quality Assurance
fLambdaMassPtRap[curCentBin]->Fill(massLambda,pt,v0->RapLambda());
fLambdaMassPtPhi->Fill(massLambda,pt,lPhi);
+ /*
if( (pt>kPtBinV0[0]) && (pt<kPtBinV0[kN1]) && isNaturalPart )
fAssocParticles->Add( new AliMiniParticle(centrality, zv, iV0, pt, lPhi, lEta, lMCAssocNegDaug, lMCAssocPosDaug, 4) );
-
+ */
+
}
// Invariant Mass cut
} // End check ID
+
fTriggerParticles->RemoveAt(iArray);
fTriggerParticles->AddAt( new AliMiniParticle(centrality, zv, idTrig, ptTrig, phiTrig, etaTrig, 0, 0, 0), iArray);
fAntiLambdaMassPtRap[curCentBin]->Fill(massAntiLambda,pt,v0->RapLambda());
fAntiLambdaMassPtPhi->Fill(massAntiLambda,pt,lPhi);
+ /*
if( (pt>kPtBinV0[0]) && (pt<kPtBinV0[kN1]) && isNaturalPart )
fAssocParticles->Add( new AliMiniParticle(centrality, zv, iV0, pt, lPhi, lEta, lMCAssocNegDaug, lMCAssocPosDaug, 5) );
+ */
}
Double_t phi=t->Phi();
- if( (pt>kPtBinCharged[0]) && (pt<kPtBinCharged[kNc]) && (TMath::Abs(eta)<fTrigEtaMax) )
+ if( (pt>=kPtBinCharged[0]) && (pt<kPtBinCharged[kNc]) && (TMath::Abs(eta)<fTrigEtaMax) )
fChargedAssocParticles->Add( new AliMiniParticle(centrality, zv, i, pt, phi, eta, 0, 0, 1) );
if( (pt>fTrigPtMin) && (pt<fTrigPtMax) && (TMath::Abs(eta)<fTrigEtaMax) ) {
// ----------------------------
// 2. Checking if the trigger particle
// might be a daughter from the V0-candidate
+ /*
for (Int_t i=0; i<(fTriggerParticles->GetEntriesFast()); i++){
AliMiniParticle* trig = (AliMiniParticle*) fTriggerParticles->At(i);
Int_t id = trig->ID();
V0Loop(kTriggerCheck,kFALSE,i,id);
}
-
+ */
}
} // End pt selection
// Xi-
+ /*
if(lPdgcodeCurrentPart == kXiMinus || lPdgcodeCurrentPart == 3322){
if( isNaturalPart )
fAssocPartMC->Add( new AliMiniParticle(centrality, zv, iTrkMC, lPtCurrentPart, lPhiCurrentPart, lEtaCurrentPart, 0, 0, 6) );
} //End Xi
-
+ */
+
} // End loop over MC
// -----------------------------------------
} // End pt bin
// Xi Minus (h-Xi correlations)
+ /*
if( (assocMC->WhichCandidate()==6) ){
fXiMinusPtMCAssoc->Fill(assocMC->Pt(),centrality);
} // End pt bin
} // End Xi selection
-
+ */
} // End loop over trigger particles
} // End loop over trigger particles
// With natural Xi's do loop over the MC trigger particle
// to see if Xi are the leading particle in the event:
// in that case, does the correlations with charged particles as associated particles
-
+ /*
Bool_t isXiLP = kFALSE;
for (Int_t iTrkMC = 0; iTrkMC < nTrkMC; iTrkMC++){
fXiTriggerPartMC->Add( new AliMiniParticle(centrality, zv, iTrkMC, p0->Pt(), p0->Phi(), p0->Eta(), 0, 0, 10) );
}
- }
+ }*/
// ----------------------------------------
// Xi-h: Do correlations with primary particles
+ /*
if(isXiLP){
for(Int_t ii=0; ii<(fXiTriggerPartMC->GetEntriesFast()); ii++){
AliMiniParticle* trigXiMC = (AliMiniParticle*) fXiTriggerPartMC->At(ii);
} // end loop over Xi
} // End condition of Xi is LP
-
+ */
} // End MC condition
// *************************************************
// Correlations
//-------------------------------------------------------------
Float_t ptTrig=0., pxTrig=0., pyTrig=0.;
+ Float_t pt=-100.; // ***
+ Float_t lPhi=0., lEta=0.; // ***
+/*
Float_t massK0s=0., mK0s=0., sK0s=0.;
Float_t massL=0., mL=0., sL=0.;
Float_t massAL=0.; //, mAL=0., sAL=0.;
Float_t lAlphaV0=0., lPtArmV0=0, dcaPos=0., dcaNeg=0.;
Float_t dx=-100., dy=-100., dz=-100., lt=-100., res=-100.;
Float_t dlK=-100., dlL=-100.;
+*/
Float_t dPhi=-100., dEta=-100., radio=-100.;
// --------------------------------
for (Int_t i=0; i<(fTriggerParticles->GetEntriesFast()); i++){
AliMiniParticle* trig = (AliMiniParticle*) fTriggerParticles->At(i);
- for(Int_t j=0; j<(fAssocParticles->GetEntriesFast()); j++){
- AliMiniParticle* tChAssoc = (AliMiniParticle*) (fAssocParticles->At(j));
+ for(Int_t j=0; j<(fChargedAssocParticles->GetEntriesFast()); j++){
+ AliMiniParticle* tChAssoc = (AliMiniParticle*) (fChargedAssocParticles->At(j));
lPhi = tChAssoc->Phi();
lEta = tChAssoc->Eta();
}
- }
+ } // End loop over trigger particles
// --------------------------------
// h-V0 correlations
+ /*
for (Int_t i=0; i<(fTriggerParticles->GetEntriesFast()); i++){
AliMiniParticle* trig = (AliMiniParticle*) fTriggerParticles->At(i);
if( trig->WhichCandidate() == 0 ) continue;
} // End AntiLambda selection
} // End loop over associated particles
-
+
} // End loop over trigger particles
-
+ */
+
//-------------------------------------------------------------
// Mixing
//-------------------------------------------------------------
if( trackAssocME->Pt()>=kPtBinCharged[k] && trackAssocME->Pt()<kPtBinCharged[k+1] )
fChargeddPhidEtaME[curCentBin*kNc*kNVtxZ + k*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);
- }
+ } // End loop over charged partilces
+
// --- V0 associated particles
+ /*
for(Int_t j=0; j<fAssocParticles->GetEntriesFast(); j++){
AliMiniParticle* trackAssocME = (AliMiniParticle*) (fAssocParticles->At(j));
else if( trackAssocME->WhichCandidate() == 5 )
fAntiLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);
- }
-
+ } // End loop over V0's
+ */
}
}