for (Int_t iTr=0; iTr<nESDTracks; iTr++)
{
+ Printf("tracj %d", iTr);
TESDTrackInfo &info = ESDTrackInfo[iTr];
info.fOK = 0;
info.fPrimUsedFlag = 0;
//* Check chi^2 for a case
- if( TMath::Sqrt( TMath::Abs(V0.GetChi2()/V0.GetNDF()) >3. )) continue;
+ if( TMath::Sqrt( TMath::Abs(V0.GetChi2()/V0.GetNDF())) >3. ) continue;
//* Get V0 decay length with estimated error
-Int_t AliKFParticleTest(Int_t n1=0,Int_t n2=1000,char *dire="/d/alice10/sma/my_v4-05-Release/pp/"){
+//Int_t AliKFParticleTest(Int_t n1=0,Int_t n2=1000,char *dire="/d/alice10/sma/my_v4-05-Release/pp/"){
+Int_t AliKFParticleTest(Int_t n1=993,Int_t n2=999,char *dire="~/SOFT/AliRoot_git/AliRoot_20140223/test/genkine/TestKF/Events/"){
//* Main macro
// LOOP OVER SERIES OF DIRECTORIES
cout <<" Opening "<<filename<<"\nand ESD "<<esdfile<<endl;
+ TFile* ftemp = TFile::Open(filename);
+ Printf("ftemp = %p", ftemp);
+ if (!ftemp) continue;
+
if (gAlice) {
delete AliRunLoader::Instance();
delete gAlice;
continue;
}
rl->LoadKinematics();
- AliTracker::SetFieldMap(gAlice->Field(),1);
- AliKFParticle::SetField( AliTracker::GetBz() );
//---------------------------------------//
// //
if (!ef || !ef->IsOpen()) {cerr<<"Error open AliESDs.root !\n"; continue ;}
//create event object
- AliESD *event = new AliESDEvent;
+ AliESDEvent *event = new AliESDEvent;
//Set pointer to the esd tree in the file
TTree* tree = (TTree*) ef->Get("esdTree");
if (!tree) {cerr<<"no ESD tree found\n"; continue;};
//Set pointer to the esd object in the tree
- event->ReadFromTree(event);
+ event->ReadFromTree(tree);
//Number of events
Int_t nevents=tree->GetEntriesFast();
for (Int_t iev=0; iev<nevents; iev++){
tree->GetEvent(iev);
rl->GetEvent(iev);
+ AliKFParticle::SetField(event->GetMagneticField());
RunV0(event);
}
EndV0();