Updated MFT version, some problems from the previous attempt are fixed (Antonio)
[u/mrichter/AliRoot.git] / MFT / AliMuonForwardTrackFinder.C
1 //================================================================================================================================
2
3 void AliMuonForwardTrackFinder(Int_t run=0,
4                                Int_t matching=0,
5                                const Char_t *readDir= ".",
6                                const Char_t *outDir = ".",
7                                Int_t nEventsToAnalyze = -1) {
8   
9   //  TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1, AliMagF::k5kG));
10
11   // AliLog::SetClassDebugLevel("AliMuonForwardTrackFinder", 1);
12
13   AliCDBManager* man = AliCDBManager::Instance();
14   man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
15   man->SetSpecificStorage("GRP/GRP/Data", Form("local://%s",gSystem->pwd()));
16
17   AliMuonForwardTrackFinder *finder = new AliMuonForwardTrackFinder();
18   finder->SetDraw(kFALSE);
19   finder->Init(run, readDir, outDir, nEventsToAnalyze);
20
21   finder -> SetSigmaSpectrometerCut(4.0);
22   finder -> SetSigmaClusterCut(4.0);
23   finder -> SetChi2GlobalCut(2.0);
24   finder -> SetRAbsorberCut(0.0);
25   //  finder -> SetRAbsorberCut(26.4);
26   finder -> SetLowPtCut(0.0);
27   //  finder -> SetLowPtCut(0.5);
28   finder -> SetExtrapOriginTransvError(0.05);   // To be imposed if the gen. vertex is fixed in (0,0,0)
29   finder -> SetGaussianBlurZVert(5.0);          // To be imposed if the gen. vertex is fixed in (0,0,0) 
30   finder -> SetMatchingMode(matching);          // 0 -> real matching   1 -> ideal matching
31   finder -> SetMinResearchRadiusAtLastPlane(0.5);
32
33   while (finder->LoadNextTrack()) continue;
34
35   if (finder->GetNRealTracksAnalyzed()) finder->Terminate();
36
37 }
38
39 //================================================================================================================================
40