[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 = 0,
				 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->SetPtDimuRange(10, 0., 5.);
  myAnalysis->SetSingleMuonAnalysis(singleMuonAnalysis);
  myAnalysis->SetMuonPairAnalysis(muonPairAnalysis);
  myAnalysis->SetOption(option);
  myAnalysis->SetMatchTrigger(kTRUE);
  myAnalysis->SetMaxNWrongClustersMC(maxNWrongClusters);
  myAnalysis->SetPtMinSingleMuons(ptMinSingleMuons);

  myAnalysis->UseCutOnOffsetChi2(kFALSE);        // cut on the single muons

  myAnalysis->UseBransonForCut(kFALSE);
  myAnalysis->UseBransonForKinematics(kFALSE);

  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,
7e3dd1af	15	Int_t myRandom = 0,
e806e863	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);
3b6fd486	32	myAnalysis->SetPtDimuRange(10, 0., 5.);
e806e863	33	myAnalysis->SetSingleMuonAnalysis(singleMuonAnalysis);
	34	myAnalysis->SetMuonPairAnalysis(muonPairAnalysis);
	35	myAnalysis->SetOption(option);
	36	myAnalysis->SetMatchTrigger(kTRUE);
	37	myAnalysis->SetMaxNWrongClustersMC(maxNWrongClusters);
	38	myAnalysis->SetPtMinSingleMuons(ptMinSingleMuons);
3b6fd486	39
	40	myAnalysis->UseCutOnOffsetChi2(kFALSE); // cut on the single muons
	41
7e3dd1af	42	myAnalysis->UseBransonForCut(kFALSE);
7e3dd1af	43	myAnalysis->UseBransonForKinematics(kFALSE);
e806e863	44
	45	myAnalysis->Init("MuonGlobalTracks.root");
	46
	47	while (myAnalysis->LoadNextEvent()) continue;
	48
	49	myAnalysis->Terminate(Form("outFiles/outFile.%d.%d.%d.root", myAnalysis->GetFirstEvent(), myAnalysis->GetLastEvent(), myRandom));
	50
	51	}
	52
	53	//================================================================================================================================
	54