//#pragma link C function GetFloatDialog;
// This class will always be for ITS only
#pragma link C++ class AliITSvtest+;
-
+#pragma link C++ class AliITSTrackerV1+;
+#pragma link C++ class AliITSgeoinfo+;
#pragma link C++ class AliITSclusterV2+;
#pragma link C++ class AliITStrackV2+;
#pragma link C++ class AliITStrackerV2+;
if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
}
-
- AliITS* ITS =(AliITS *)gAlice->GetDetector("ITS");
- if (!ITS) return;
+ AliITS* IITTSS =(AliITS *)gAlice->GetDetector("ITS"); // nuova modificata
+ if (!IITTSS) return; // nuova modificata
+ AliITSTrackerV1* ITStracker = new AliITSTrackerV1(IITTSS); // nuova
//
// Loop over events
TStopwatch timer;
timer.Start();
- ITS->DoTracking(nev,min_t,max_t,file,flagvert);
+ //ITS->DoTracking(nev,min_t,max_t,file,flagvert); // nuova eliminata
+ ITStracker->DoTracking(nev,min_t,max_t,file,flagvert); // nuova
timer.Stop(); timer.Print();
} // event loop
file->Close();
+#include "AliITSgeom.h"
struct GoodTrack {
Int_t lab;
Int_t code;
Int_t numbpoints=0;
AliITSRecPoint *recp=0;
Int_t track=-3;
-
+ // Int_t modmin[]={0,80,240,324,500,1282};
+ // Int_t modmax[]={79,239,323,499,1281,2269};
+
+
+ Int_t modmin[6], modmax[6];
+ Int_t Nladder[6], Ndetector[6];
+ AliITSgeom *g1 = ITS->GetITSgeom();
+
+ for(Int_t i=0; i<6; i++) {
+ Nladder[i]=g1->GetNladders(i+1);
+ Ndetector[i]=g1->GetNdetectors(i+1);
+ //cout<<Nladder[i]<<" "<<Ndetector[i]<<"\n";
+ }
+ for(Int_t i=0; i<6; i++) {
+ modmin[i]=g1->GetModuleIndex(i+1,1,1);
+ modmax[i]=g1->GetModuleIndex(i+1,Nladder[i],Ndetector[i]);
+ //cout<<modmin[i]<<" "<<modmax[i]<<"\n";
+ }
+
+
for (int nev=evNumber1; nev<= evNumber2; nev++) {
Int_t nparticles = gAlice->GetEvent(nev);
// Int_t np=particles->GetEntriesFast();
Int_t np=gAlice->GetNtrack(); //FCA correction
TMatrix goodITS(np,6);
- Int_t modmin[]={0,80,240,324,500,1282};
- Int_t modmax[]={79,239,323,499,1281,2269};
- for(Int_t j=0; j<np; j++){
- for(Int_t j1=0; j1<6; j1++) goodITS(j,j1)=0;
- }
+ //Int_t modmin[]={0,80,240,324,500,1282};
+ //Int_t modmax[]={79,239,323,499,1281,2269};
+ // for(Int_t j=0; j<np; j++){
+ // for(Int_t j1=0; j1<6; j1++) goodITS(j,j1)=0;
+ // }
Int_t nent=gAlice->TreeR()->GetEntries();