//
// Loop over events
//
-
char tname[30];
+ for (int nev=evNumber1; nev<= evNumber2; nev++) {
+
- for (int nev=0; nev<= evNumber2; nev++) {
+ // for (int nev=0; nev<= evNumber2; nev++) {
sprintf(tname,"TreeT%d",nev);
TTree *tracktree=(TTree*)file->Get(tname);
TBranch *tbranch=tracktree->GetBranch("ITStracks");
+ cout<<" nev = "<<nev<<"\n";
//cout<<" open the file \n";
Int_t nentr=tracktree->GetEntries();
TObjArray tarray(nentr);
- AliITSiotrack *iotrack=0;
+ // AliITSIOTrack *iotrack=0;
printf("nentr %d\n",nentr);
for (Int_t i=0; i<nentr; i++) {
- AliITSiotrack *iotrack=new AliITSiotrack;
- // tarray.AddAt(new AliITSiotrack,i);
+ AliITSIOTrack *iotrack=new AliITSIOTrack;
+ // tarray.AddAt(new AliITSIOTrack,i);
// iotrack=(AliITSiotrack*)tarray.UncheckedAt(i);
tbranch->SetAddress(&iotrack);
tracktree->GetEvent(i);
tarray.AddLast(iotrack);
}
- file->Close();
+ //file->Close();
- AliITSiotrack *iotrack;
- for (Int_t i=0; i<nentr; i++) {
- iotrack=(AliITSiotrack*)tarray.UncheckedAt(i);
+ AliITSIOTrack *iotrack;
+ for (Int_t i=0; i<nentr; i++) {
+ AliITSIOTrack *iotrack=new AliITSIOTrack;
+ iotrack=(AliITSIOTrack*)tarray.UncheckedAt(i);
if(!iotrack) continue;
Int_t label=iotrack->GetLabel();
Double_t phistate=iotrack->GetStatePhi();
Double_t Dr=iotrack->GetStateD();
Double_t C=iotrack->GetStateC();
- Double_t C00,C10,C11,C20,C21,C22,C30,C31,C32,C33,C40,C41,C42,C43,C44;
-
- iotrack->GetCovMatrix(C00,C10,C11,C20,C21,C22,C30,C31,C32,C33,C40,C41,C42,C43,C44);
-
-
- cout<<"\n"<<C00<<" "<<C10<<" "<<C11<<" "<<C20<<" "<<C21<<" "<<C22<<" "
- <<C30<<" "<<C31<<" "<<C32<<" "<<C33<<" "<<C40<<" "<<C41<<" "<<C42<<" "
- <<C43<<" "<<C44<<"\n\n\n";
-
- cout<<" label = "<<label<<"\n";
- cout<< " phi z D tgl C = "<<phistate<<" "<<Zstate<<" "<<Dr<<" "<<tgl<<" "<<C<<"\n"; getchar();
+ // cout<<" track label = "<<label<<"\n";
+ // cout<<" phi z D tanl C = "<<phistate<<" "<<Zstate<<" "<<Dr<<" "<<tgl<<" "<<C<<"\n";
printf(" i %d label x y z %d %f %f %f\n",i,label,x,y,z);
// delete iotrack;
*/
-
+ delete iotrack;
}
- } // event loop
+ } // event loop
+ file->Close();
}