[u/mrichter/AliRoot.git] / MFT / AliMuonForwardTrackAnalysis.C

enum {kNoOption, kOpenFlavor, kResonanceOnly};

//================================================================================================================================

void AliMuonForwardTrackAnalysis(const Char_t *readDir= ".",
				 Int_t option = kNoOption,
				 Int_t nMassBin = 100, 
				 Double_t massMin = 0.,
				 Double_t massMax = 10.,
				 const Char_t *outDir = ".",
				 Bool_t singleMuonAnalysis = kTRUE,
				 Bool_t muonPairAnalysis = kTRUE,
				 Int_t firstEvent = -1,
				 Int_t lastEvent = -1, 
				 Int_t myRandom,
				 Int_t maxNWrongClusters = 999,
				 Double_t ptMinSingleMuons = 0.0) {
  
  gROOT -> LoadMacro("./AliMuonForwardTrackAnalysis.cxx+");
  //  AliLog::SetClassDebugLevel("AliMuonForwardTrackPair", 1);
  //  AliLog::SetClassDebugLevel("AliMuonForwardTrackAnalysis", 1);

  //  TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG, AliMagF::kBeamTypeAA, 2750.));
  //  TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG, AliMagF::kBeamTypepp, 7000.));
  TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
  
  AliMuonForwardTrackAnalysis *myAnalysis = new AliMuonForwardTrackAnalysis();
  myAnalysis->ReadEvents(firstEvent, lastEvent);
  myAnalysis->SetInputDir(readDir);
  myAnalysis->SetOutputDir(outDir);
  myAnalysis->SetMassRange(nMassBin, massMin, massMax);
  myAnalysis->SetSingleMuonAnalysis(singleMuonAnalysis);
  myAnalysis->SetMuonPairAnalysis(muonPairAnalysis);
  myAnalysis->SetOption(option);
  myAnalysis->SetMatchTrigger(kTRUE);
  myAnalysis->SetMaxNWrongClustersMC(maxNWrongClusters);
  myAnalysis->SetPtMinSingleMuons(ptMinSingleMuons);

  myAnalysis->Init("MuonGlobalTracks.root");

  while (myAnalysis->LoadNextEvent()) continue;

  myAnalysis->Terminate(Form("outFiles/outFile.%d.%d.%d.root", myAnalysis->GetFirstEvent(), myAnalysis->GetLastEvent(), myRandom));

}

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