//
fLoader = AliRunLoader::Open(fnGalice);
if (gAlice){
- //delete gAlice->GetRunLoader();
+ //delete AliRunLoader::Instance();
delete gAlice;
gAlice = 0x0;
}
fEventNr = nall-firstEvent;
cerr<<"restricted number of events availaible"<<endl;
}
- AliMagF * magf = gAlice->Field();
- AliTracker::SetFieldMap(magf,0);
TGeoManager::Import("geometry.root");
AliESDtrack * trackn = fEvent->GetTrack((v0MI->GetNindex()));
AliESDtrack * trackp = fEvent->GetTrack((v0MI->GetPindex()));
Int_t labels[2]={-1,-1};
- labels[0] = (trackn==0) ? -1 : trackn->GetLabel();
- labels[1] = (trackp==0) ? -1 : trackp->GetLabel();
+ labels[0] = (trackn==0) ? -1 : TMath::Abs(trackn->GetLabel());
+ labels[1] = (trackp==0) ? -1 : TMath::Abs(trackp->GetLabel());
//
for (Int_t i=0;i<2;i++){
- Int_t absLabel = labels[i];
+ Int_t absLabel = TMath::Abs(labels[i]);
if (absLabel < fNParticles) {
if (fMultiRecV0[absLabel]>0){
if (fMultiRecV0[absLabel]<20)
cerr<<"fNParticles, nParticlesTR, fNextTreeGenEntryToRead: "<<fNParticles<<" "
<<nParticlesTR<<" "<<fNextTreeGenEntryToRead<<endl;
TBranch * branch = fTreeCmp->GetBranch("RC");
- TBranch * branchF = fTreeCmp->GetBranch("F");
+ // TBranch * branchF = fTreeCmp->GetBranch("F");
branch->SetAddress(&fRecInfo); // set all pointers
fRecArray = new TObjArray(fNParticles);
while (entry < nParticlesTR) {
fTreeGenV0->GetEntry(entry);
entry++;
+ fRecV0Info->Reset(); //reset all variables
if (eventNr < fGenV0Info->GetMinus().fEventNr) continue;
if (eventNr > fGenV0Info->GetMinus().fEventNr) continue;;
//
AliESDtrack * trackp = fEvent->GetTrack((v0MI2->GetPindex()));
Int_t vlabeln = (trackn==0) ? -1 : trackn->GetLabel();
Int_t vlabelp = (trackp==0) ? -1 : trackp->GetLabel();
+ fRecV0Info->fLab[0]=TMath::Abs(vlabelp);
+ fRecV0Info->fLab[1]=TMath::Abs(vlabeln);
//
if (TMath::Abs(vlabeln)==label &&TMath::Abs(vlabelp)==label2) {
if (v0MI2->GetOnFlyStatus()) {
if (fSignedV0[i]==0){
AliV0 *v0MI = (AliV0*)fEvent->GetV0(i);
if (!v0MI) continue;
+ fRecV0Info->Reset(); //reset all variables
//
new (fRecV0Info->fV0rec) AliV0(*v0MI);
fRecV0Info->fV0Status =-10;
AliESDtrack * trackp = fEvent->GetTrack((v0MI->GetPindex()));
Int_t vlabeln = (trackn==0) ? -1 : trackn->GetLabel();
Int_t vlabelp = (trackp==0) ? -1 : trackp->GetLabel();
-
+ fRecV0Info->fLab[0]=TMath::Abs(vlabelp);
+ fRecV0Info->fLab[1]=TMath::Abs(vlabeln);
+ if (TMath::Abs(fRecV0Info->fLab[0] - fRecV0Info->fLab[1])<2) continue;
AliESDRecInfo* fRecInfo1 = (AliESDRecInfo*)fRecArray->At(TMath::Abs(vlabeln));
AliESDRecInfo* fRecInfo2 = (AliESDRecInfo*)fRecArray->At(TMath::Abs(vlabelp));
if (fRecInfo1 && fRecInfo2){