delete fHRapMuon;
delete fHRapResonance;
delete fHPtResonance;
- delete fHInvMassAll_vs_Pt;
+ delete fHInvMassAllvsPt;
}
/*********************************************************/
Int_t AliMuonAnalysis::Init()
{
- //Initilizes anaysis
+ //Initilizes analysis
Info("Init","Histo initialized for MUON Analysis");
fHistoFile = new TFile("MUONmassPlot.root", "RECREATE");
fHRapMuon = new TH1F("hRapMuon"," Muon Rapidity",50,-4.5,-2);
fHRapResonance = new TH1F("hRapResonance"," Resonance Rapidity",50,-4.5,-2);
fHPtResonance = new TH1F("hPtResonance", "Resonance Pt (GeV/c)", 100, 0., 20.);
- fHInvMassAll_vs_Pt = new TH2F("hInvMassAll_vs_Pt","hInvMassAll_vs_Pt",480,0.,12.,80,0.,20.);
+ fHInvMassAllvsPt = new TH2F("hInvMassAll_vs_Pt","hInvMassAll_vs_Pt",480,0.,12.,80,0.,20.);
return 0;
}
Int_t AliMuonAnalysis::ProcessEvent(AliAOD* aodrec, AliAOD* aodsim)
{
-
+ //
+ // process the event
+ //
if (aodrec) {
GetInvMass(aodrec);
// Info("ProcessEvent","Inv Mass Rec");
void AliMuonAnalysis::GetInvMass(AliAOD* aod)
{
-
+ // get the invariant mass distribution
+ // from the oad events
+
TLorentzVector lorV1, lorV2, lorVtot;
Float_t massMin = 9.17;
Float_t massMax = 9.77;
Float_t invMass = lorVtot.M();
fHInvMassAll->Fill(invMass);
- fHInvMassAll_vs_Pt->Fill(invMass,lorVtot.Pt());
+ fHInvMassAllvsPt->Fill(invMass,lorVtot.Pt());
if (invMass > massMin && invMass < massMax) {
fHRapResonance->Fill(lorVtot.Rapidity());
private:
- TFile *fHistoFile;
- TH1F *fHPtMuon;
- TH1F *fHPtMuonPlus;
- TH1F *fHPtMuonMinus;
- TH1F *fHPMuon;
- TH1F *fHInvMassAll;
- TH1F *fHRapMuon;
- TH1F *fHRapResonance;
- TH1F *fHPtResonance;
- TH2F *fHInvMassAll_vs_Pt;
+ TFile *fHistoFile; // histogramm file pointer
+ TH1F *fHPtMuon; // Muon Pt distribution
+ TH1F *fHPtMuonPlus; // Muon Plus Pt distribution
+ TH1F *fHPtMuonMinus; // Muon Minus Pt distribution
+ TH1F *fHPMuon; // Muon momentum distribution
+ TH1F *fHInvMassAll; // Invariant mass distribution
+ TH1F *fHRapMuon; // Muon rapidity distribution
+ TH1F *fHRapResonance; // Muon rapidity distribution around resonance
+ TH1F *fHPtResonance; // Muon Pt distribution around resonance
+ TH2F *fHInvMassAllvsPt; // Invariant mass vs Pt distribution
AliAODParticleCut* fPartCut;//Particle Cut
ClassDef(AliMuonAnalysis,1)