1 /**************************************************************************
2 * This file is property of and copyright by the ALICE HLT Project *
3 * ALICE Experiment at CERN, All rights reserved. *
5 * Primary Author: Svein Lindal <slindal@fys.uio.no> *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
16 /// @file AliAnaConvCorrPion.cxx
17 /// @author Svein Lindal
18 /// @brief Class used to run conversion gamma/pion - hadron/jet analysis
23 #include "AliAnaConvCorrPion.h"
24 //#include "AliAODTrack.h"
25 #include "TClonesArray.h"
26 #include "AliAODConversionParticle.h"
27 //#include "AliAODConversionMother.h"
28 //#include "AliAODConversionPhoton.h"
29 //#include "THnSparse.h"
36 ClassImp(AliAnaConvCorrPion)
38 //________________________________________________________________________________
39 AliAnaConvCorrPion::AliAnaConvCorrPion() :
40 AliAnaConvCorrBase("pion_hadron_corr", "Pion dPhi"),
41 //hTriggerPtvsMass(NULL),
47 //________________________________________________________________________________
48 AliAnaConvCorrPion::AliAnaConvCorrPion(TString name, TString title = "Pion Corr") :
49 AliAnaConvCorrBase(name, title),
50 //hTriggerPtvsMass(NULL),
58 //________________________________________________________________________________
59 AliAnaConvCorrPion::~AliAnaConvCorrPion() {
63 //________________________________________________________________________________
64 void AliAnaConvCorrPion::InitMassAxis() {
65 Double_t mbins[7] = {0.1, 0.11, 0.12, 0.15, 0.16, 0.18, 0.2};
67 fAxisM.SetNameTitle("InvMass", "invariant mass");
68 GetAxisList().AddAt(&fAxisM, 5);
69 GetTrackAxisList().AddAt(&fAxisM, 5);
70 GetTrigAxisList().AddAt(&fAxisM, 4);
73 ///________________________________________________________________________________
74 void AliAnaConvCorrPion::CreateHistograms() {
76 CreateBaseHistograms();
78 hTriggerPtvsMass[0] = new TH2D("hTriggerPtvsMass_all", "Pt vs Mass all pizero", 400, 0, .400, GetAxistPt().GetNbins(), GetAxistPt().GetXbins()->GetArray());
79 hTriggerPtvsMass[1] = new TH2D("hTriggerPtvsMass_leadingcone", "Pt vs Mass leading cone", 400, 0, .400, GetAxistPt().GetNbins(), GetAxistPt().GetXbins()->GetArray());
80 hTriggerPtvsMass[2] = new TH2D("hTriggerPtvsMass_leadingevent", "Pt vs Mass leading event", 400, 0, .400, GetAxistPt().GetNbins(), GetAxistPt().GetXbins()->GetArray());
81 GetHistograms()->Add(hTriggerPtvsMass[0]);
82 GetHistograms()->Add(hTriggerPtvsMass[1]);
83 GetHistograms()->Add(hTriggerPtvsMass[2]);
87 ///________________________________________________________________________________
88 void AliAnaConvCorrPion::FillTriggerCounters(const AliAODConversionParticle * particle, Int_t leading) {
89 //Fill histograms counting triggers
90 //fHNTriggers[leading]->Fill(particle->Pt());
91 hTriggerPtvsMass[leading]->Fill(particle->M(), particle->Pt());
94 //________________________________________________________________________________
95 // void AliAnaConvCorrPion::Process(TClonesArray * pions, TClonesArray * photons, TClonesArray * tracks) {
97 // for(Int_t ip = 0; ip < pions->GetEntriesFast(); ip++) {
99 // AliAODConversionParticle * pion = static_cast<AliAODConversionParticle*>(pions->UncheckedAt(ip));
101 // Int_t tIDs[4] = {-1, -1, -1, -1};
102 // AliAODConversionParticle * photon1 = static_cast<AliAODConversionParticle*>(photons->UncheckedAt(pion->GetLabel(0)));
103 // tIDs[0] = photon1->GetLabel(0);
104 // tIDs[1] = photon1->GetLabel(1);
105 // AliAODConversionParticle * photon2 = static_cast<AliAODConversionParticle*>(photons->UncheckedAt(pion->GetLabel(1)));
106 // tIDs[2] = photon2->GetLabel(0);
107 // tIDs[3] = photon2->GetLabel(1);
109 // CorrelateWithTracks(static_cast<AliAODConversionParticle*>(pion), tracks, tIDs, kFALSE);