1 enum {kNoOption, kPionsKaons, kResonanceOnly};
3 //================================================================================================================================
5 void AliMuonForwardTrackAnalysis(const Char_t *readDir= ".",
6 Int_t option = kNoOption,
9 Double_t massMax = 10.,
10 Bool_t useCutOnOffsetChi2 = kFALSE,
11 Int_t maxNWrongClusters = 999,
12 const Char_t *outDir = ".",
13 Bool_t singleMuonAnalysis = kTRUE,
14 Bool_t muonPairAnalysis = kTRUE,
15 Int_t firstEvent = -1,
18 Double_t ptMinSingleMuons = 0.0) {
20 gROOT -> LoadMacro("./AliMuonForwardTrackAnalysis.cxx+");
21 // AliLog::SetClassDebugLevel("AliMuonForwardTrackPair", 1);
22 // AliLog::SetClassDebugLevel("AliMuonForwardTrackAnalysis", 1);
24 // TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG, AliMagF::kBeamTypeAA, 2750.));
25 // TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG, AliMagF::kBeamTypepp, 7000.));
26 TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
28 AliMuonForwardTrackAnalysis *myAnalysis = new AliMuonForwardTrackAnalysis();
29 myAnalysis->ReadEvents(firstEvent, lastEvent);
30 myAnalysis->SetInputDir(readDir);
31 myAnalysis->SetOutputDir(outDir);
32 myAnalysis->SetMassRange(nMassBin, massMin, massMax);
33 myAnalysis->SetSingleMuonAnalysis(singleMuonAnalysis);
34 myAnalysis->SetMuonPairAnalysis(muonPairAnalysis);
35 myAnalysis->SetOption(option);
36 myAnalysis->SetMaxNWrongClustersMC(maxNWrongClusters);
37 myAnalysis->SetPtMinSingleMuons(ptMinSingleMuons);
38 myAnalysis->UseCutOnOffsetChi2(useCutOnOffsetChi2);
40 myAnalysis->SetPtDimuRange(10, 0., 5.);
41 myAnalysis->SetMatchTrigger(kTRUE);
42 myAnalysis->UseBransonForCut(kFALSE);
43 myAnalysis->UseBransonForKinematics(kFALSE);
45 myAnalysis->Init("MuonGlobalTracks.root");
47 while (myAnalysis->LoadNextEvent()) continue;
49 myAnalysis->Terminate(Form("outFiles/outFile.%d.%d.%d.root", myAnalysis->GetFirstEvent(), myAnalysis->GetLastEvent(), myRandom));
53 //================================================================================================================================