PWGGA/GammaConv/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 "TH2F.h"
  23 #include "AliAnaConvCorrPion.h"
  24 #include "AliAODConversionParticle.h"
  25 #include <iostream>
  26
  27
  28 using namespace std;
  29 ClassImp(AliAnaConvCorrPion)
  30
  31 //________________________________________________________________________________
  32 AliAnaConvCorrPion::AliAnaConvCorrPion() :
  33 AliAnaConvCorrBase("pion_hadron_corr", "Pion dPhi"),
  34   fAxisM()
  35 {
  36   //consctructor
  37   InitMassAxis();
  38 }
  39 //________________________________________________________________________________
  40 AliAnaConvCorrPion::AliAnaConvCorrPion(TString name, TString title = "Pion Corr") :
  41   AliAnaConvCorrBase(name, title),
  42   fAxisM()
  43 {
  44   //consctructor
  45   InitMassAxis();
  46 }
  47
  48
  49 //________________________________________________________________________________
  50 AliAnaConvCorrPion::~AliAnaConvCorrPion() {
  51   //destructor
  52 }
  53
  54 //________________________________________________________________________________
  55 void AliAnaConvCorrPion::InitMassAxis() {
  56   //Init mass axis
  57
  58   Double_t mbins[7] = {0.1, 0.11, 0.12, 0.15, 0.16, 0.18, 0.2};
  59   fAxisM.Set(6, mbins);
  60   fAxisM.SetNameTitle("InvMass", "invariant mass");
  61   GetAxisList().AddAt(&fAxisM, 5);
  62   GetTrackAxisList().AddAt(&fAxisM, 5);
  63   GetTrigAxisList().AddAt(&fAxisM, 4);
  64 }
  65
  66 ///________________________________________________________________________________
  67 void AliAnaConvCorrPion::CreateHistograms() {
  68   //Create histograms
  69   CreateBaseHistograms();
  70
  71   for(Int_t i = 0; i < 3; i++) {
  72         TH2F * hpthist = new TH2F(Form("%s_iso_%d_ptVsmass", fName.Data(), i),
  73                                                           Form("%s_iso_%d_ptVsmass", fName.Data(), i),
  74                                                           400, 0, .400,
  75                                                           GetAxistPt().GetNbins(), GetAxistPt().GetXbins()->GetArray());
  76         GetHistograms()->Add(hpthist);
  77   }
  78 }
  79
  80
  81 ///________________________________________________________________________________
  82 void AliAnaConvCorrPion::FillTriggerCounters(const AliAODConversionParticle * particle, Int_t leading) {
  83   //Fill histograms counting triggers
  84   TH2F *hpt = dynamic_cast<TH2F*>(GetHistograms()->FindObject(Form("%s_iso_%d_ptVsmass", fName.Data(), leading)));
  85   if(hpt) hpt->Fill(particle->M(), particle->Pt());
  86 }