7265ce057daee4632bc76f9cd8faf79599db694b
[u/mrichter/AliRoot.git] / PWGGA / GammaConv / ConvCorrelations / AliAnaConvCorrPion.cxx
1 /**************************************************************************
2  * This file is property of and copyright by the ALICE HLT Project        *
3  * ALICE Experiment at CERN, All rights reserved.                         *
4  *                                                                        *
5  * Primary Author: Svein Lindal <slindal@fys.uio.no>                      *
6  *                                                                        *
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  **************************************************************************/
15
16 /// @file   AliAnaConvCorrPion.cxx
17 /// @author Svein Lindal
18 /// @brief  Class used to run conversion gamma/pion - hadron/jet analysis
19
20
21
22 #include "TH2D.h"
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"
30 //#include "TH2F.h"
31
32 #include <iostream>
33
34
35 using namespace std;
36 ClassImp(AliAnaConvCorrPion)
37
38 //________________________________________________________________________________
39 AliAnaConvCorrPion::AliAnaConvCorrPion() :
40 AliAnaConvCorrBase("pion_hadron_corr", "Pion dPhi"),
41   hTriggerPtvsMass(NULL),
42   fAxisM()
43 {
44   //consctructor
45   InitMassAxis();
46 }
47 //________________________________________________________________________________
48 AliAnaConvCorrPion::AliAnaConvCorrPion(TString name, TString title = "Pion Corr") :
49   AliAnaConvCorrBase(name, title),
50   hTriggerPtvsMass(NULL),
51   fAxisM()
52 {
53   //consctructor
54   InitMassAxis();
55 }
56
57
58 //________________________________________________________________________________
59 AliAnaConvCorrPion::~AliAnaConvCorrPion() {
60   //destructor
61 }
62
63 //________________________________________________________________________________
64 void AliAnaConvCorrPion::InitMassAxis() {
65   Double_t mbins[7] = {0.1, 0.11, 0.12, 0.15, 0.16, 0.18, 0.2};
66   fAxisM.Set(6, mbins);
67   fAxisM.SetNameTitle("InvMass", "invariant mass");
68   GetAxisList().AddAt(&fAxisM, 5);
69   GetTrackAxisList().AddAt(&fAxisM, 5);
70   GetTrigAxisList().AddAt(&fAxisM, 4);
71 }
72
73 ///________________________________________________________________________________
74 void AliAnaConvCorrPion::CreateHistograms() {
75   //Create histograms
76   CreateBaseHistograms();
77   hTriggerPtvsMass = new TH2D("hTriggerPtvsMass", "Pt vs Mass", 400, 0, .400, GetAxistPt().GetNbins(), GetAxistPt().GetXbins()->GetArray());
78   GetHistograms()->Add(hTriggerPtvsMass);
79 }
80
81
82 ///________________________________________________________________________________
83 void AliAnaConvCorrPion::FillTriggerCounters(const AliAODConversionParticle * particle, Bool_t isolated) {
84   //Fill histograms counting triggers
85   fHNTriggers[isolated]->Fill(particle->Pt());
86   hTriggerPtvsMass->Fill(particle->M(), particle->Pt());
87 }
88
89 //________________________________________________________________________________
90 // void AliAnaConvCorrPion::Process(TClonesArray * pions, TClonesArray * photons, TClonesArray * tracks) {
91
92 //   for(Int_t ip = 0; ip < pions->GetEntriesFast(); ip++) {
93
94 //      AliAODConversionParticle * pion = static_cast<AliAODConversionParticle*>(pions->UncheckedAt(ip));
95         
96 //      Int_t tIDs[4] = {-1, -1, -1, -1};
97 //      AliAODConversionParticle * photon1 = static_cast<AliAODConversionParticle*>(photons->UncheckedAt(pion->GetLabel(0)));
98 //      tIDs[0] =  photon1->GetLabel(0);
99 //      tIDs[1] =  photon1->GetLabel(1);
100 //      AliAODConversionParticle * photon2 = static_cast<AliAODConversionParticle*>(photons->UncheckedAt(pion->GetLabel(1)));
101 //      tIDs[2] =  photon2->GetLabel(0);
102 //      tIDs[3] =  photon2->GetLabel(1);
103         
104 //      CorrelateWithTracks(static_cast<AliAODConversionParticle*>(pion), tracks, tIDs, kFALSE);
105 //   }
106 // }
107