PWGLF/totEt/AliAnalysisHadEtReconstructed.cxx

   1 //_________________________________________________________________________
   2 //  Utility Class for transverse energy studies, charged hadrons
   3 //  Base class for ESD analysis
   4 //  - reconstruction output
   5 // implementation file
   6 //
   7 //Created by Christine Nattrass, Rebecca Scott, Irakli Martashvili
   8 //University of Tennessee at Knoxville
   9 //_________________________________________________________________________
  10
  11 #include <TROOT.h>
  12 #include <TSystem.h>
  13 #include <TInterpreter.h>
  14 #include "AliAnalysisHadEtReconstructed.h"
  15 #include "AliAnalysisEtCuts.h"
  16 #include "AliESDtrack.h"
  17 #include "AliESDCaloCluster.h"
  18 #include "AliVEvent.h"
  19 #include "AliESDEvent.h"
  20 #include "AliESDtrackCuts.h"
  21 #include "AliESDpid.h"
  22 #include "AliVParticle.h"
  23 #include <iostream>
  24 #include "AliAnalysisHadEtCorrections.h"
  25 #include "TFile.h"
  26 #include "TString.h"
  27 #include "AliAnalysisEtCommon.h"
  28 #include "AliAnalysisHadEt.h"
  29 #include "AliCentrality.h"
  30 #include "AliLog.h"
  31 #include "AliPWG0Helper.h"
  32
  33 using namespace std;
  34
  35 ClassImp(AliAnalysisHadEtReconstructed);
  36
  37
  38 AliAnalysisHadEtReconstructed::AliAnalysisHadEtReconstructed() :
  39   AliAnalysisHadEt()
  40   ,fCorrections(0)
  41   ,fConfigFile("ConfigHadEtAnalysis.C")
  42   ,fCorrTotEtFullAcceptanceTPC(0)
  43   ,fCorrTotEtFullAcceptanceITS(0)
  44   ,fCorrHadEtFullAcceptanceTPC(0)
  45   ,fCorrHadEtFullAcceptanceITS(0)
  46   ,fCorrPiKPEtFullAcceptanceTPC(0)
  47   ,fCorrPiKPEtFullAcceptanceITS(0)
  48   ,fCorrTotEtEMCALAcceptanceTPC(0)
  49   ,fCorrTotEtEMCALAcceptanceITS(0)
  50   ,fCorrHadEtEMCALAcceptanceTPC(0)
  51   ,fCorrHadEtEMCALAcceptanceITS(0)
  52   ,fCorrTotEtPHOSAcceptanceTPC(0)
  53   ,fCorrTotEtPHOSAcceptanceITS(0)
  54   ,fCorrHadEtPHOSAcceptanceTPC(0)
  55   ,fCorrHadEtPHOSAcceptanceITS(0)
  56   ,fCorrectedHadEtFullAcceptanceTPCNoPID(0)
  57   ,fCorrectedHadEtFullAcceptanceITSNoPID(0)
  58   ,fCorrectedHadEtEMCALAcceptanceTPCNoPID(0)
  59   ,fCorrectedHadEtEMCALAcceptanceITSNoPID(0)
  60   ,fCorrectedHadEtPHOSAcceptanceTPCNoPID(0)
  61   ,fCorrectedHadEtPHOSAcceptanceITSNoPID(0)
  62   ,fCorrectedHadEtFullAcceptanceTPC(0)
  63   ,fCorrectedHadEtFullAcceptanceITS(0)
  64   ,fCorrectedHadEtFullAcceptanceTPCAssumingPion(0)
  65   ,fCorrectedHadEtFullAcceptanceITSAssumingPion(0)
  66   ,fCorrectedHadEtFullAcceptanceTPCAssumingProton(0)
  67   ,fCorrectedHadEtFullAcceptanceITSAssumingProton(0)
  68   ,fCorrectedHadEtFullAcceptanceTPCAssumingKaon(0)
  69   ,fCorrectedHadEtFullAcceptanceITSAssumingKaon(0)
  70   ,fCorrectedHadEtEMCALAcceptanceTPC(0)
  71   ,fCorrectedHadEtEMCALAcceptanceITS(0)
  72   ,fCorrectedHadEtPHOSAcceptanceTPC(0)
  73   ,fCorrectedHadEtPHOSAcceptanceITS(0)
  74   ,fRawEtFullAcceptanceTPC(0)
  75   ,fRawEtFullAcceptanceITS(0)
  76   ,fRawEtEMCALAcceptanceTPC(0)
  77   ,fRawEtEMCALAcceptanceITS(0)
  78   ,fRawEtPHOSAcceptanceTPC(0)
  79   ,fRawEtPHOSAcceptanceITS(0)
  80   ,fRawEtFullAcceptanceTPCNoPID(0)
  81   ,fRawEtFullAcceptanceITSNoPID(0)
  82   ,fRawEtEMCALAcceptanceTPCNoPID(0)
  83   ,fRawEtEMCALAcceptanceITSNoPID(0)
  84   ,fRawEtPHOSAcceptanceTPCNoPID(0)
  85   ,fRawEtPHOSAcceptanceITSNoPID(0)
  86   ,kIsOfflineV0AND(0)
  87   ,kDoTriggerChecks(0)
  88   ,kDoTriggerChecksOnly(0)
  89 {
  90 }
  91
  92 AliAnalysisHadEtReconstructed::~AliAnalysisHadEtReconstructed()
  93 {
  94   delete fCorrections;
  95 }
  96
  97 Int_t AliAnalysisHadEtReconstructed::AnalyseEvent(AliVEvent* ev, Int_t eventtype)
  98 { // analyse ESD event
  99   if(kDoTriggerChecksOnly){return 1;}//In this case we are just after trigger efficiencies and don't care about the ET reconstructed.
 100   if(kDoTriggerChecks && !kIsOfflineV0AND){return 1;}//In this case we are just after trigger efficiencies and don't care about the ET reconstructed.
 101   ResetEventValues();
 102   if(!ev){
 103     AliFatal("ERROR: Event does not exist");
 104     return 0;
 105   }
 106
 107   AliESDEvent *realEvent = dynamic_cast<AliESDEvent*>(ev);
 108   if(!realEvent){
 109     AliFatal("ERROR: ESD Event does not exist");
 110     return 0;
 111   }
 112   fCentBin= -1;
 113   fGoodEvent = kTRUE;//for p+p collisions if we made it this far we have a good event
 114   if(fDataSet==20100){//If this is Pb+Pb
 115     AliCentrality *centrality = realEvent->GetCentrality();
 116     if(fNCentBins<21) fCentBin= centrality->GetCentralityClass10(fCentralityMethod);
 117     else{ fCentBin= centrality->GetCentralityClass5(fCentralityMethod);}
 118     if(fCentBin ==-1) fGoodEvent = kFALSE;//but for Pb+Pb events we don't want to count events where we did not find a centrality
 119   }
 120   //for PID
 121   AliESDpid *pID = new AliESDpid();
 122   pID->MakePID(realEvent);
 123   TString *strTPC = new TString("TPC");
 124   TString *strITS = new TString("ITS");
 125   TString *strTPCITS = new TString("TPCITS");
 126   for(Int_t cutset=0;cutset<2;cutset++){
 127     bool isTPC = false;
 128     TString *cutName = NULL;
 129     TObjArray* list = NULL;
 130     switch(cutset){
 131     case 0:
 132       cutName = strTPCITS;
 133       list = fEsdtrackCutsITSTPC->GetAcceptedTracks(realEvent);
 134       isTPC = true;
 135       break;
 136     case 1:
 137       cutName = strITS;
 138       list = fEsdtrackCutsITS->GetAcceptedTracks(realEvent);
 139       break;
 140     case 2:
 141       cutName = strTPC;
 142       list = fEsdtrackCutsTPC->GetAcceptedTracks(realEvent);
 143       break;
 144     default:
 145       cerr<<"Error:  cannot fill histograms!"<<endl;
 146       return -1;
 147     }
 148     Int_t nGoodTracks = list->GetEntries();
 149     for (Int_t iTrack = 0; iTrack < nGoodTracks; iTrack++)
 150       {
 151
 152
 153         AliESDtrack *track = dynamic_cast<AliESDtrack*> (list->At(iTrack));
 154         if (!track)
 155           {
 156             Printf("ERROR: Could not get track %d", iTrack);
 157             continue;
 158           }
 159         else{
 160           if(TMath::Abs(track->Eta())>fCorrections->GetEtaCut()) continue;
 161           Float_t nSigmaPion,nSigmaProton,nSigmaKaon,nSigmaElectron;
 162           pID->MakeTPCPID(track);
 163           pID->MakeITSPID(track);
 164           if(cutset!=1){
 165             nSigmaPion = TMath::Abs(pID->NumberOfSigmasTPC(track,AliPID::kPion));
 166             nSigmaProton = TMath::Abs(pID->NumberOfSigmasTPC(track,AliPID::kProton));
 167             nSigmaKaon = TMath::Abs(pID->NumberOfSigmasTPC(track,AliPID::kKaon));
 168             nSigmaElectron = TMath::Abs(pID->NumberOfSigmasTPC(track,AliPID::kElectron));
 169           }
 170           else{
 171             nSigmaPion = TMath::Abs(pID->NumberOfSigmasITS(track,AliPID::kPion));
 172             nSigmaProton = TMath::Abs(pID->NumberOfSigmasITS(track,AliPID::kProton));
 173             nSigmaKaon = TMath::Abs(pID->NumberOfSigmasITS(track,AliPID::kKaon));
 174             nSigmaElectron = TMath::Abs(pID->NumberOfSigmasITS(track,AliPID::kElectron));
 175           }
 176           //        bool isPion = (nSigmaPion<3.0 && nSigmaProton>2.0 && nSigmaKaon>2.0);
 177           //        bool isElectron = (nSigmaElectron<2.0 && nSigmaPion>4.0 && nSigmaProton>3.0 && nSigmaKaon>3.0);
 178           //        bool isKaon = (nSigmaPion>3.0 && nSigmaProton>2.0 && nSigmaKaon<2.0);
 179           //        bool isProton = (nSigmaPion>3.0 && nSigmaProton<2.0 && nSigmaKaon>2.0);
 180           bool isPion = (nSigmaPion<3.0 && nSigmaProton>2.0 && nSigmaKaon>2.0);
 181           bool isElectron = (nSigmaElectron<2.0 && nSigmaPion>4.0 && nSigmaProton>3.0 && nSigmaKaon>3.0);
 182           bool isKaon = (nSigmaPion>3.0 && nSigmaProton>3.0 && nSigmaKaon<3.0 && track->Pt()<0.45);
 183           bool isProton = (nSigmaPion>3.0 && nSigmaProton<3.0 && nSigmaKaon>3.0 && track->Pt()<0.9);
 184
 185           bool unidentified = (!isProton && !isKaon && !isElectron && !isPion);
 186           if(cutset==1){//ITS dE/dx identification requires tighter cuts on the tracks and we don't gain much from that so we won't do it
 187             unidentified = true;
 188             isPion=false;
 189             isElectron=false;
 190             isKaon=false;
 191             isProton=false;
 192           }
 193           Float_t dEdx = track->GetTPCsignal();
 194           if(cutset==1) dEdx = track->GetITSsignal();
 195           FillHisto2D(Form("dEdxDataAll%s",cutName->Data()),track->P(),dEdx,1.0);
 196
 197           bool inPHOS = IsInPHOS(track);
 198           bool inEMCAL = IsInEMCAL(track);
 199
 200           Float_t corrBkgd=0.0;
 201           Float_t corrNotID=0.0;
 202           Float_t corrNoID = fCorrections->GetNotIDCorrectionNoPID(track->Pt());
 203           Float_t corrEff = 0.0;
 204           Float_t corrEffNoID = 0.0;
 205           if(cutset!=1){//TPC
 206             corrBkgd = fCorrections->GetBackgroundCorrectionTPC(track->Pt());
 207             corrEffNoID = fCorrections->GetTPCEfficiencyCorrectionHadron(track->Pt(),fCentBin);
 208             corrNotID = fCorrections->GetNotIDConstCorrectionTPC();
 209             corrNoID = fCorrections->GetNotIDConstCorrectionTPCNoID();
 210           }
 211           if(cutset==1){//ITS
 212             corrBkgd = fCorrections->GetBackgroundCorrectionITS(track->Pt());
 213             corrEffNoID = fCorrections->GetITSEfficiencyCorrectionHadron(track->Pt(),fCentBin);
 214             corrNotID = fCorrections->GetNotIDConstCorrectionITS();
 215             corrNoID = fCorrections->GetNotIDConstCorrectionITSNoID();
 216           }
 217           if(fDataSet==20100){
 218             FillHisto2D("fbkgdVsCentralityBin",fCentBin,corrBkgd,1.0);
 219             FillHisto2D("fnotIDVsCentralityBin",fCentBin,corrNotID,1.0);
 220             FillHisto2D("fpTcutVsCentralityBin",fCentBin,fCorrections->GetpTCutCorrectionTPC(),1.0);
 221             if(fCorrHadEtFullAcceptanceTPC>0.0) FillHisto2D("fneutralVsCentralityBin",fCentBin,1.0/fCorrHadEtFullAcceptanceTPC,1.0);
 222             if(fCorrections->GetNeutralCorrection()>0.0) FillHisto2D("ConstantCorrectionsVsCentralityBin",fCentBin,1.0/fCorrections->GetNeutralCorrection(),1.0);
 223           }
 224           Float_t et = 0.0;
 225           Float_t etNoID = Et(track->P(),track->Theta(),fgPiPlusCode,track->Charge());
 226           Float_t etpartialcorrected = 0.0;
 227           Float_t etpartialcorrectedPion = 0.0;
 228           Float_t etpartialcorrectedKaon = 0.0;
 229           Float_t etpartialcorrectedProton = 0.0;
 230           Float_t etpartialcorrectedNoID = corrNoID*corrBkgd*corrEffNoID*etNoID;
 231           FillHisto2D(Form("EtDataRaw%sNoID",cutName->Data()),track->Pt(),track->Eta(),etpartialcorrectedNoID);
 232
 233           if(isPion){
 234             FillHisto2D(Form("dEdxDataPion%s",cutName->Data()),track->P(),dEdx,1.0);
 235             et = Et(track->P(),track->Theta(),fgPiPlusCode,track->Charge());
 236             if(cutset==0){corrEff = fCorrections->GetTPCEfficiencyCorrectionPion(track->Pt(),fCentBin);}
 237             etpartialcorrected = et*corrBkgd*corrEff*corrNotID;
 238             if(corrEff>0.0&&fDataSet==20100)FillHisto2D("feffPionVsCentralityBin",fCentBin,1.0/corrEff,1.0);
 239             if(track->Charge()>0.0){
 240               FillHisto2D(Form("EtDataRaw%sPiPlus",cutName->Data()),track->Pt(),track->Eta(),et);
 241               FillHisto2D(Form("EtDataCorrected%sPiPlus",cutName->Data()),track->Pt(),track->Eta(),etpartialcorrected);
 242             }
 243             else{
 244               FillHisto2D(Form("EtDataRaw%sPiMinus",cutName->Data()),track->Pt(),track->Eta(),et);
 245               FillHisto2D(Form("EtDataCorrected%sPiMinus",cutName->Data()),track->Pt(),track->Eta(),etpartialcorrected);
 246             }
 247           }
 248           if(isKaon){
 249             FillHisto2D(Form("dEdxDataKaon%s",cutName->Data()),track->P(),dEdx,1.0);
 250             et = Et(track->P(),track->Theta(),fgKPlusCode,track->Charge());
 251             if(cutset==0){corrEff = fCorrections->GetTPCEfficiencyCorrectionKaon(track->Pt(),fCentBin);}
 252             etpartialcorrected = et*corrBkgd*corrEff*corrNotID;
 253             if(corrEff>0.0&&fDataSet==20100)FillHisto2D("feffKaonVsCentralityBin",fCentBin,1.0/corrEff,1.0);
 254
 255             if(track->Charge()>0.0){
 256               FillHisto2D(Form("EtDataRaw%sKPlus",cutName->Data()),track->Pt(),track->Eta(),et);
 257               FillHisto2D(Form("EtDataCorrected%sKPlus",cutName->Data()),track->Pt(),track->Eta(),etpartialcorrected);
 258             }
 259             else{
 260               FillHisto2D(Form("EtDataRaw%sKMinus",cutName->Data()),track->Pt(),track->Eta(),et);
 261               FillHisto2D(Form("EtDataCorrected%sKMinus",cutName->Data()),track->Pt(),track->Eta(),etpartialcorrected);
 262             }
 263           }
 264           if(isProton){
 265             FillHisto2D(Form("dEdxDataProton%s",cutName->Data()),track->P(),dEdx,1.0);
 266             et = Et(track->P(),track->Theta(),fgProtonCode,track->Charge());
 267             if(cutset==0){corrEff = fCorrections->GetTPCEfficiencyCorrectionProton(track->Pt(),fCentBin);}
 268             etpartialcorrected = et*corrBkgd*corrEff*corrNotID;
 269             if(corrEff>0.0&&fDataSet==20100)FillHisto2D("feffProtonVsCentralityBin",fCentBin,1.0/corrEff,1.0);
 270
 271             if(track->Charge()>0.0){
 272               FillHisto2D(Form("EtDataRaw%sProton",cutName->Data()),track->Pt(),track->Eta(),et);
 273               FillHisto2D(Form("EtDataCorrected%sProton",cutName->Data()),track->Pt(),track->Eta(),etpartialcorrected);
 274             }
 275             else{
 276               FillHisto2D(Form("EtDataRaw%sAntiProton",cutName->Data()),track->Pt(),track->Eta(),et);
 277               FillHisto2D(Form("EtDataCorrected%sAntiProton",cutName->Data()),track->Pt(),track->Eta(),etpartialcorrected);
 278             }
 279           }
 280           if(isElectron){
 281             FillHisto2D(Form("dEdxDataElectron%s",cutName->Data()),track->P(),dEdx,1.0);
 282           }
 283           if(unidentified){
 284             if(isPion) cerr<<"I should not be here!!  AliAnalysisHadEtReconstructed 273"<<endl;
 285             FillHisto2D(Form("dEdxDataUnidentified%s",cutName->Data()),track->P(),dEdx,1.0);
 286             et = Et(track->P(),track->Theta(),fgPiPlusCode,track->Charge());
 287             Float_t etProton = Et(track->P(),track->Theta(),fgProtonCode,track->Charge());
 288             Float_t etKaon = Et(track->P(),track->Theta(),fgKPlusCode,track->Charge());
 289             if(corrEff>0.0&&fDataSet==20100)FillHisto2D("feffHadronVsCentralityBin",fCentBin,1.0/corrEff,1.0);
 290             etpartialcorrected = et*corrBkgd*corrEffNoID*corrNotID;
 291             etpartialcorrectedPion = et*corrBkgd*corrEffNoID;
 292             etpartialcorrectedProton = etProton*corrBkgd*corrEffNoID;
 293             etpartialcorrectedKaon = etKaon*corrBkgd*corrEffNoID;
 294             FillHisto2D(Form("EtDataCorrected%sUnidentified",cutName->Data()),track->Pt(),track->Eta(),etpartialcorrected);
 295           }
 296           else{
 297             etpartialcorrectedPion = etpartialcorrected;
 298             etpartialcorrectedKaon = etpartialcorrected;
 299             etpartialcorrectedProton = etpartialcorrected;
 300           }
 301           if(!isTPC){
 302             etpartialcorrected = etpartialcorrectedNoID;//Not using PID for ITS
 303           }
 304           AddEt(et,etNoID,etpartialcorrected,etpartialcorrectedPion,etpartialcorrectedProton,etpartialcorrectedKaon,etpartialcorrectedNoID,track->Pt(),isTPC,inPHOS,inEMCAL);
 305         }
 306       }
 307     delete list;
 308   }
 309
 310   Int_t nondiff = 0;//(Int_t) AliPWG0Helper::kND;
 311   Int_t doublediff = 0;//(Int_t) AliPWG0Helper::kDD;
 312   Int_t singlediff = 0;//(Int_t) AliPWG0Helper::kSD;
 313   if(fDataSet!=20100){
 314     nondiff = (Int_t) AliPWG0Helper::kND;
 315     doublediff = (Int_t) AliPWG0Helper::kDD;
 316     singlediff = (Int_t) AliPWG0Helper::kSD;
 317   }
 318   if(eventtype == nondiff && fGoodEvent){
 319     FillHisto1D("RecoHadEtFullAcceptanceTPCND",GetCorrectedHadEtFullAcceptanceTPC(),1.0);
 320     FillHisto1D("RecoPiKPEtFullAcceptanceTPCND",GetCorrectedPiKPEtFullAcceptanceTPC(),1.0);
 321     FillHisto1D("RecoTotEtFullAcceptanceTPCND",GetCorrectedTotEtFullAcceptanceTPC(),1.0);
 322     FillHisto1D("RecoHadEtFullAcceptanceITSND",GetCorrectedHadEtFullAcceptanceITS(),1.0);
 323     FillHisto1D("RecoPiKPEtFullAcceptanceITSND",GetCorrectedPiKPEtFullAcceptanceITS(),1.0);
 324     FillHisto1D("RecoTotEtFullAcceptanceITSND",GetCorrectedTotEtFullAcceptanceITS(),1.0);
 325     FillHisto1D("RecoHadEtFullAcceptanceTPCNoPIDND",GetCorrectedHadEtFullAcceptanceTPCNoPID(),1.0);
 326     FillHisto1D("RecoPiKPEtFullAcceptanceTPCNoPIDND",GetCorrectedPiKPEtFullAcceptanceTPCNoPID(),1.0);
 327     FillHisto1D("RecoTotEtFullAcceptanceTPCNoPIDND",GetCorrectedTotEtFullAcceptanceTPCNoPID(),1.0);
 328     FillHisto1D("RecoHadEtFullAcceptanceITSNoPIDND",GetCorrectedHadEtFullAcceptanceITSNoPID(),1.0);
 329     FillHisto1D("RecoPiKPEtFullAcceptanceITSNoPIDND",GetCorrectedPiKPEtFullAcceptanceITSNoPID(),1.0);
 330     FillHisto1D("RecoTotEtFullAcceptanceITSNoPIDND",GetCorrectedTotEtFullAcceptanceITSNoPID(),1.0);
 331   }
 332   if(eventtype == doublediff){
 333     FillHisto1D("RecoHadEtFullAcceptanceTPCDD",GetCorrectedHadEtFullAcceptanceTPC(),1.0);
 334     FillHisto1D("RecoPiKPEtFullAcceptanceTPCDD",GetCorrectedPiKPEtFullAcceptanceTPC(),1.0);
 335     FillHisto1D("RecoTotEtFullAcceptanceTPCDD",GetCorrectedTotEtFullAcceptanceTPC(),1.0);
 336     FillHisto1D("RecoPiKPEtFullAcceptanceITSDD",GetCorrectedPiKPEtFullAcceptanceITS(),1.0);
 337     FillHisto1D("RecoHadEtFullAcceptanceITSDD",GetCorrectedHadEtFullAcceptanceITS(),1.0);
 338     FillHisto1D("RecoTotEtFullAcceptanceITSDD",GetCorrectedTotEtFullAcceptanceITS(),1.0);
 339     FillHisto1D("RecoHadEtFullAcceptanceTPCNoPIDDD",GetCorrectedHadEtFullAcceptanceTPCNoPID(),1.0);
 340     FillHisto1D("RecoPiKPEtFullAcceptanceTPCNoPIDDD",GetCorrectedPiKPEtFullAcceptanceTPCNoPID(),1.0);
 341     FillHisto1D("RecoTotEtFullAcceptanceTPCNoPIDDD",GetCorrectedTotEtFullAcceptanceTPCNoPID(),1.0);
 342     FillHisto1D("RecoHadEtFullAcceptanceITSNoPIDDD",GetCorrectedHadEtFullAcceptanceITSNoPID(),1.0);
 343     FillHisto1D("RecoPiKPEtFullAcceptanceITSNoPIDDD",GetCorrectedPiKPEtFullAcceptanceITSNoPID(),1.0);
 344     FillHisto1D("RecoTotEtFullAcceptanceITSNoPIDDD",GetCorrectedTotEtFullAcceptanceITSNoPID(),1.0);
 345   }
 346   if(eventtype == singlediff){
 347     FillHisto1D("RecoPiKPEtFullAcceptanceTPCSD",GetCorrectedPiKPEtFullAcceptanceTPC(),1.0);
 348     FillHisto1D("RecoHadEtFullAcceptanceTPCSD",GetCorrectedHadEtFullAcceptanceTPC(),1.0);
 349     FillHisto1D("RecoTotEtFullAcceptanceTPCSD",GetCorrectedTotEtFullAcceptanceTPC(),1.0);
 350     FillHisto1D("RecoHadEtFullAcceptanceITSSD",GetCorrectedHadEtFullAcceptanceITS(),1.0);
 351     FillHisto1D("RecoPiKPEtFullAcceptanceITSSD",GetCorrectedPiKPEtFullAcceptanceITS(),1.0);
 352     FillHisto1D("RecoTotEtFullAcceptanceITSSD",GetCorrectedTotEtFullAcceptanceITS(),1.0);
 353     FillHisto1D("RecoHadEtFullAcceptanceTPCNoPIDSD",GetCorrectedHadEtFullAcceptanceTPCNoPID(),1.0);
 354     FillHisto1D("RecoPiKPEtFullAcceptanceTPCNoPIDSD",GetCorrectedPiKPEtFullAcceptanceTPCNoPID(),1.0);
 355     FillHisto1D("RecoTotEtFullAcceptanceTPCNoPIDSD",GetCorrectedTotEtFullAcceptanceTPCNoPID(),1.0);
 356     FillHisto1D("RecoHadEtFullAcceptanceITSNoPIDSD",GetCorrectedHadEtFullAcceptanceITSNoPID(),1.0);
 357     FillHisto1D("RecoPiKPEtFullAcceptanceITSNoPIDSD",GetCorrectedPiKPEtFullAcceptanceITSNoPID(),1.0);
 358     FillHisto1D("RecoTotEtFullAcceptanceITSNoPIDSD",GetCorrectedTotEtFullAcceptanceITSNoPID(),1.0);
 359   }
 360   if(fGoodEvent){
 361     FillHisto1D("RecoPiKPEtFullAcceptanceTPC",GetCorrectedPiKPEtFullAcceptanceTPC(),1.0);
 362     FillHisto1D("RecoPiKPEtFullAcceptanceITS",GetCorrectedPiKPEtFullAcceptanceITS(),1.0);
 363     FillHisto1D("RecoPiKPEtFullAcceptanceTPCNoPID",GetCorrectedPiKPEtFullAcceptanceTPCNoPID(),1.0);
 364     FillHisto1D("RecoPiKPEtFullAcceptanceITSNoPID",GetCorrectedPiKPEtFullAcceptanceITSNoPID(),1.0);
 365     FillHisto1D("RecoHadEtFullAcceptanceTPC",GetCorrectedHadEtFullAcceptanceTPC(),1.0);
 366     FillHisto1D("RecoTotEtFullAcceptanceTPC",GetCorrectedTotEtFullAcceptanceTPC(),1.0);
 367     FillHisto1D("RecoHadEtFullAcceptanceTPCAssumingPion",GetCorrectedHadEtFullAcceptanceTPCAssumingPion(),1.0);
 368     FillHisto1D("RecoHadEtFullAcceptanceTPCAssumingProton",GetCorrectedHadEtFullAcceptanceTPCAssumingProton(),1.0);
 369     FillHisto1D("RecoHadEtFullAcceptanceTPCAssumingKaon",GetCorrectedHadEtFullAcceptanceTPCAssumingKaon(),1.0);
 370     FillHisto1D("RecoHadEtEMCALAcceptanceTPC",GetCorrectedHadEtEMCALAcceptanceTPC(),1.0);
 371     FillHisto1D("RecoTotEtEMCALAcceptanceTPC",GetCorrectedTotEtEMCALAcceptanceTPC(),1.0);
 372     FillHisto1D("RecoHadEtPHOSAcceptanceTPC",GetCorrectedHadEtPHOSAcceptanceTPC(),1.0);
 373     FillHisto1D("RecoTotEtPHOSAcceptanceTPC",GetCorrectedTotEtPHOSAcceptanceTPC(),1.0);
 374     FillHisto1D("RecoHadEtFullAcceptanceTPCNoPID",GetCorrectedHadEtFullAcceptanceTPCNoPID(),1.0);
 375     FillHisto1D("RecoTotEtFullAcceptanceTPCNoPID",GetCorrectedTotEtFullAcceptanceTPCNoPID(),1.0);
 376     FillHisto1D("RecoHadEtEMCALAcceptanceTPCNoPID",GetCorrectedHadEtEMCALAcceptanceTPCNoPID(),1.0);
 377     FillHisto1D("RecoTotEtEMCALAcceptanceTPCNoPID",GetCorrectedTotEtEMCALAcceptanceTPCNoPID(),1.0);
 378     FillHisto1D("RecoHadEtPHOSAcceptanceTPCNoPID",GetCorrectedHadEtPHOSAcceptanceTPCNoPID(),1.0);
 379     FillHisto1D("RecoTotEtPHOSAcceptanceTPCNoPID",GetCorrectedTotEtPHOSAcceptanceTPCNoPID(),1.0);
 380     FillHisto1D("RecoHadEtFullAcceptanceITS",GetCorrectedHadEtFullAcceptanceITS(),1.0);
 381     FillHisto1D("RecoHadEtFullAcceptanceITSAssumingPion",GetCorrectedHadEtFullAcceptanceITSAssumingPion(),1.0);
 382     FillHisto1D("RecoHadEtFullAcceptanceITSAssumingProton",GetCorrectedHadEtFullAcceptanceITSAssumingProton(),1.0);
 383     FillHisto1D("RecoHadEtFullAcceptanceITSAssumingKaon",GetCorrectedHadEtFullAcceptanceITSAssumingKaon(),1.0);
 384     FillHisto1D("RecoTotEtFullAcceptanceITS",GetCorrectedTotEtFullAcceptanceITS(),1.0);
 385     FillHisto1D("RecoHadEtEMCALAcceptanceITS",GetCorrectedHadEtEMCALAcceptanceITS(),1.0);
 386     FillHisto1D("RecoTotEtEMCALAcceptanceITS",GetCorrectedTotEtEMCALAcceptanceITS(),1.0);
 387     FillHisto1D("RecoHadEtPHOSAcceptanceITS",GetCorrectedHadEtPHOSAcceptanceITS(),1.0);
 388     FillHisto1D("RecoTotEtPHOSAcceptanceITS",GetCorrectedTotEtPHOSAcceptanceITS(),1.0);
 389     FillHisto1D("RecoHadEtFullAcceptanceITSNoPID",GetCorrectedHadEtFullAcceptanceITSNoPID(),1.0);
 390     FillHisto1D("RecoTotEtFullAcceptanceITSNoPID",GetCorrectedTotEtFullAcceptanceITSNoPID(),1.0);
 391     FillHisto1D("RecoHadEtEMCALAcceptanceITSNoPID",GetCorrectedHadEtEMCALAcceptanceITSNoPID(),1.0);
 392     FillHisto1D("RecoTotEtEMCALAcceptanceITSNoPID",GetCorrectedTotEtEMCALAcceptanceITSNoPID(),1.0);
 393     FillHisto1D("RecoHadEtPHOSAcceptanceITSNoPID",GetCorrectedHadEtPHOSAcceptanceITSNoPID(),1.0);
 394     FillHisto1D("RecoTotEtPHOSAcceptanceITSNoPID",GetCorrectedTotEtPHOSAcceptanceITSNoPID(),1.0);
 395
 396     FillHisto1D("RecoRawEtFullAcceptanceTPC",GetRawEtFullAcceptanceTPC(),1.0);
 397     FillHisto1D("RecoRawEtEMCALAcceptanceTPC",GetRawEtEMCALAcceptanceTPC(),1.0);
 398     FillHisto1D("RecoRawEtPHOSAcceptanceTPC",GetRawEtPHOSAcceptanceTPC(),1.0);
 399     FillHisto1D("RecoRawEtFullAcceptanceTPCNoPID",GetRawEtFullAcceptanceTPCNoPID(),1.0);
 400     FillHisto1D("RecoRawEtEMCALAcceptanceTPCNoPID",GetRawEtEMCALAcceptanceTPCNoPID(),1.0);
 401     FillHisto1D("RecoRawEtPHOSAcceptanceTPCNoPID",GetRawEtPHOSAcceptanceTPCNoPID(),1.0);
 402     FillHisto1D("RecoRawEtFullAcceptanceITS",GetRawEtFullAcceptanceITS(),1.0);
 403     FillHisto1D("RecoRawEtEMCALAcceptanceITS",GetRawEtEMCALAcceptanceITS(),1.0);
 404     FillHisto1D("RecoRawEtPHOSAcceptanceITS",GetRawEtPHOSAcceptanceITS(),1.0);
 405     FillHisto1D("RecoRawEtFullAcceptanceITSNoPID",GetRawEtFullAcceptanceITSNoPID(),1.0);
 406     FillHisto1D("RecoRawEtEMCALAcceptanceITSNoPID",GetRawEtEMCALAcceptanceITSNoPID(),1.0);
 407     FillHisto1D("RecoRawEtPHOSAcceptanceITSNoPID",GetRawEtPHOSAcceptanceITSNoPID(),1.0);
 408     if(fCentBin>-1){//if we have Pb+Pb and found a centrality bin
 409       FillHisto1D(Form("RecoHadEtFullAcceptanceTPCCB%i",fCentBin),GetCorrectedHadEtFullAcceptanceTPC(),1.0);
 410       FillHisto1D(Form("RecoTotEtFullAcceptanceTPCCB%i",fCentBin),GetCorrectedTotEtFullAcceptanceTPC(),1.0);
 411       FillHisto1D(Form("RecoHadEtFullAcceptanceITSCB%i",fCentBin),GetCorrectedHadEtFullAcceptanceITS(),1.0);
 412       FillHisto1D(Form("RecoTotEtFullAcceptanceITSCB%i",fCentBin),GetCorrectedTotEtFullAcceptanceITS(),1.0);
 413       FillHisto1D(Form("RecoRawEtFullAcceptanceTPCCB%i",fCentBin),GetRawEtFullAcceptanceTPC(),1.0);
 414       FillHisto1D(Form("RecoRawEtFullAcceptanceITSCB%i",fCentBin),GetRawEtFullAcceptanceITS(),1.0);
 415     }
 416   }
 417   delete pID;
 418   delete strTPC;
 419   delete strITS;
 420   delete strTPCITS;
 421   // cout<<"Reconstructed pi/k/p et "<<GetCorrectedPiKPEtFullAcceptanceTPC()<<endl;
 422   return 1;
 423 }
 424 void AliAnalysisHadEtReconstructed::AddEt(Float_t rawEt, Float_t rawEtNoPID, Float_t corrEt, Float_t corrEtPion, Float_t corrEtProton, Float_t corrEtKaon, Float_t corrEtNoPID, Float_t pt, Bool_t IsTPC, Bool_t InPHOS, Bool_t InEMCAL) {//Adding Et to each of the variables that tracks et event by event
 425   if(pt>=AliAnalysisHadEt::fgPtTPCCutOff && IsTPC){//TPC tracks
 426     //adding to the raw Et
 427     fRawEtFullAcceptanceTPC += rawEt;
 428     if(InPHOS)fRawEtPHOSAcceptanceTPC += rawEt;
 429     if(InEMCAL)fRawEtEMCALAcceptanceTPC += rawEt;
 430     fRawEtFullAcceptanceTPCNoPID += rawEtNoPID;
 431     if(InPHOS)fRawEtPHOSAcceptanceTPCNoPID += rawEtNoPID;
 432     if(InEMCAL)fRawEtEMCALAcceptanceTPCNoPID += rawEtNoPID;
 433     //adding to the corrected Et
 434     fCorrectedHadEtFullAcceptanceTPC += corrEt;//the pi/k/p et
 435     fCorrectedHadEtFullAcceptanceTPCAssumingPion += corrEtPion;
 436     fCorrectedHadEtFullAcceptanceTPCAssumingProton += corrEtProton;
 437     fCorrectedHadEtFullAcceptanceTPCAssumingKaon += corrEtKaon;
 438     if(InPHOS)fCorrectedHadEtPHOSAcceptanceTPC += corrEt;
 439     if(InEMCAL)fCorrectedHadEtEMCALAcceptanceTPC += corrEt;
 440     fCorrectedHadEtFullAcceptanceTPCNoPID += corrEtNoPID;
 441     if(InPHOS)fCorrectedHadEtPHOSAcceptanceTPCNoPID += corrEtNoPID;
 442     if(InEMCAL)fCorrectedHadEtEMCALAcceptanceTPCNoPID += corrEtNoPID;
 443   }
 444   //if(pt<AliAnalysisHadEt::fgPtTPCCutOff &&pt>=AliAnalysisHadEt::fgPtITSCutOff && !IsTPC){//ITS tracks
 445   //If we use standalone tracks - not pure standalone tracks - the only tracks we get are ones that were missed by the TPC+ITS tracking.  Therefore we don't need to add a momentum cut-off
 446   if(pt<AliAnalysisHadEt::fgPtTPCCutOff && !IsTPC){//ITS tracks
 447     //adding to the raw Et
 448     fRawEtFullAcceptanceITS += rawEt;
 449     if(InPHOS)fRawEtPHOSAcceptanceITS += rawEt;
 450     if(InEMCAL)fRawEtEMCALAcceptanceITS += rawEt;
 451     fRawEtFullAcceptanceITSNoPID += rawEtNoPID;
 452     if(InPHOS)fRawEtPHOSAcceptanceITSNoPID += rawEtNoPID;
 453     if(InEMCAL)fRawEtEMCALAcceptanceITSNoPID += rawEtNoPID;
 454     //adding to the corrected Et
 455     fCorrectedHadEtFullAcceptanceITS += corrEt;
 456     fCorrectedHadEtFullAcceptanceITSAssumingPion += corrEtPion;
 457     fCorrectedHadEtFullAcceptanceITSAssumingProton += corrEtProton;
 458     fCorrectedHadEtFullAcceptanceITSAssumingKaon += corrEtKaon;
 459     if(InPHOS)fCorrectedHadEtPHOSAcceptanceITS += corrEt;
 460     if(InEMCAL)fCorrectedHadEtEMCALAcceptanceITS += corrEt;
 461     fCorrectedHadEtFullAcceptanceITSNoPID += corrEtNoPID;
 462     if(InPHOS)fCorrectedHadEtPHOSAcceptanceITSNoPID += corrEtNoPID;
 463     if(InEMCAL)fCorrectedHadEtEMCALAcceptanceITSNoPID += corrEtNoPID;
 464   }
 465 }
 466
 467 Bool_t AliAnalysisHadEtReconstructed::IsInPHOS(AliESDtrack *track){//This function will need to be elaborated on later to include PHOS dead channels
 468   if(!track){
 469     cout<<"Error: Track does not exist!!"<<endl;
 470     return kFALSE;
 471   }
 472   return   TMath::Abs(track->Eta()) < fCuts->GetGeometryPhosEtaAccCut()//in eta acceptance
 473     && track->Phi()*180.0/TMath::Pi() > fCuts->GetGeometryPhosPhiAccMinCut()//greater than the minimum phi
 474     && track->Phi()*180.0/TMath::Pi() < fCuts->GetGeometryPhosPhiAccMaxCut();//less than the maximum phi
 475 }
 476 Bool_t AliAnalysisHadEtReconstructed::IsInEMCAL(AliESDtrack *track){//This function will need to be elaborated on later to include EMCAL dead channels
 477   if(!track){
 478     cout<<"Error: Track does not exist!!"<<endl;
 479     return kFALSE;
 480   }
 481   return   TMath::Abs(track->Eta()) < fCuts->GetGeometryEmcalEtaAccCut()//in eta acceptance
 482     && track->Phi()*180.0/TMath::Pi() > fCuts->GetGeometryEmcalPhiAccMinCut()//greater than the minimum phi
 483     && track->Phi()*180.0/TMath::Pi() < fCuts->GetGeometryEmcalPhiAccMaxCut();//less than the maximum phi
 484 }
 485 Bool_t AliAnalysisHadEtReconstructed::CheckGoodVertex(AliVParticle* track)
 486 { // check vertex
 487
 488   Float_t bxy = 999.;
 489   Float_t bz = 999.;
 490   if(!track){
 491     AliError("ERROR: no track");
 492     return kFALSE;
 493   }
 494   AliESDtrack *esdTrack = dynamic_cast<AliESDtrack*>(track);
 495   if(!esdTrack){
 496     AliError("ERROR: no track");
 497     return kFALSE;
 498   }
 499   esdTrack->GetImpactParametersTPC(bxy,bz);
 500
 501   bool status = (TMath::Abs(track->Xv()) < fCuts->GetReconstructedVertexXCut()) &&
 502     (TMath::Abs(track->Yv()) < fCuts->GetReconstructedVertexYCut()) &&
 503     (TMath::Abs(track->Zv()) < fCuts->GetReconstructedVertexZCut()) &&
 504     (TMath::Abs(bxy) < fCuts->GetReconstructedIPxyCut()) &&
 505     (TMath::Abs(bz) < fCuts->GetReconstructedIPzCut());
 506
 507   return status;
 508 }
 509
 510 void AliAnalysisHadEtReconstructed::Init()
 511 { // Init
 512   AliAnalysisHadEt::Init();
 513   if(fCorrections){
 514     fCorrTotEtFullAcceptanceTPC = fCorrections->GetConstantCorrections(kTRUE,fgPtTPCCutOff,"Full");
 515     fCorrTotEtFullAcceptanceITS = fCorrections->GetConstantCorrections(kTRUE,fgPtITSCutOff,"Full");
 516     fCorrHadEtFullAcceptanceTPC = fCorrections->GetConstantCorrections(kFALSE,fgPtTPCCutOff,"Full");
 517     fCorrHadEtFullAcceptanceITS = fCorrections->GetConstantCorrections(kFALSE,fgPtITSCutOff,"Full");
 518     fCorrPiKPEtFullAcceptanceTPC = fCorrections->GetConstantCorrections(kFALSE,fgPtTPCCutOff,"PiKP");
 519     fCorrPiKPEtFullAcceptanceITS = fCorrections->GetConstantCorrections(kFALSE,fgPtITSCutOff,"PiKP");
 520     fCorrTotEtEMCALAcceptanceTPC = fCorrections->GetConstantCorrections(kTRUE,fgPtTPCCutOff,"EMCAL");
 521     fCorrTotEtEMCALAcceptanceITS = fCorrections->GetConstantCorrections(kTRUE,fgPtITSCutOff,"EMCAL");
 522     fCorrHadEtEMCALAcceptanceTPC = fCorrections->GetConstantCorrections(kFALSE,fgPtTPCCutOff,"EMCAL");
 523     fCorrHadEtEMCALAcceptanceITS = fCorrections->GetConstantCorrections(kFALSE,fgPtITSCutOff,"EMCAL");
 524     fCorrTotEtPHOSAcceptanceTPC = fCorrections->GetConstantCorrections(kTRUE,fgPtTPCCutOff,"PHOS");
 525     fCorrTotEtPHOSAcceptanceITS = fCorrections->GetConstantCorrections(kTRUE,fgPtITSCutOff,"PHOS");
 526     fCorrHadEtPHOSAcceptanceTPC = fCorrections->GetConstantCorrections(kFALSE,fgPtTPCCutOff,"PHOS");
 527     fCorrHadEtPHOSAcceptanceITS = fCorrections->GetConstantCorrections(kFALSE,fgPtITSCutOff,"PHOS");
 528   }
 529   else{
 530     cout<<"Warning!  You have not set corrections.  Your code will crash.  You have to set the corrections."<<endl;
 531   }
 532 }
 533
 534 void AliAnalysisHadEtReconstructed::ResetEventValues(){//resetting event by event et's
 535   AliAnalysisHadEt::ResetEventValues();
 536   fCorrectedHadEtFullAcceptanceTPCNoPID=0.0;
 537   fCorrectedHadEtFullAcceptanceITSNoPID=0.0;
 538   fCorrectedHadEtEMCALAcceptanceTPCNoPID=0.0;
 539   fCorrectedHadEtEMCALAcceptanceITSNoPID=0.0;
 540   fCorrectedHadEtPHOSAcceptanceTPCNoPID=0.0;
 541   fCorrectedHadEtPHOSAcceptanceITSNoPID=0.0;
 542   fCorrectedHadEtFullAcceptanceTPC=0.0;
 543   fCorrectedHadEtFullAcceptanceITS=0.0;
 544   fCorrectedHadEtFullAcceptanceTPCAssumingPion=0.0;
 545   fCorrectedHadEtFullAcceptanceITSAssumingPion=0.0;
 546   fCorrectedHadEtFullAcceptanceTPCAssumingProton=0.0;
 547   fCorrectedHadEtFullAcceptanceITSAssumingProton=0.0;
 548   fCorrectedHadEtFullAcceptanceTPCAssumingKaon=0.0;
 549   fCorrectedHadEtFullAcceptanceITSAssumingKaon=0.0;
 550   fCorrectedHadEtEMCALAcceptanceTPC=0.0;
 551   fCorrectedHadEtEMCALAcceptanceITS=0.0;
 552   fCorrectedHadEtPHOSAcceptanceTPC=0.0;
 553   fCorrectedHadEtPHOSAcceptanceITS=0.0;
 554   fRawEtFullAcceptanceTPC=0.0;
 555   fRawEtFullAcceptanceITS=0.0;
 556   fRawEtEMCALAcceptanceTPC=0.0;
 557   fRawEtEMCALAcceptanceITS=0.0;
 558   fRawEtPHOSAcceptanceTPC=0.0;
 559   fRawEtPHOSAcceptanceITS=0.0;
 560   fRawEtFullAcceptanceTPCNoPID=0.0;
 561   fRawEtFullAcceptanceITSNoPID=0.0;
 562   fRawEtEMCALAcceptanceTPCNoPID=0.0;
 563   fRawEtEMCALAcceptanceITSNoPID=0.0;
 564   fRawEtPHOSAcceptanceTPCNoPID=0.0;
 565   fRawEtPHOSAcceptanceITSNoPID=0.0;
 566
 567   if(TMath::Abs(fCorrTotEtFullAcceptanceTPC)<1e-3){
 568     if (fConfigFile.Length()) {
 569       cout<<"Warning: Rereading fCorrections file..."<<endl;
 570       gROOT->LoadMacro(fConfigFile);
 571       fCorrections = (AliAnalysisHadEtCorrections *) gInterpreter->ProcessLine("ConfigHadEtAnalysis()");
 572       fCorrTotEtFullAcceptanceTPC = fCorrections->GetConstantCorrections(kTRUE,fgPtTPCCutOff,"Full");
 573       fCorrTotEtFullAcceptanceITS = fCorrections->GetConstantCorrections(kTRUE,fgPtITSCutOff,"Full");
 574       fCorrHadEtFullAcceptanceTPC = fCorrections->GetConstantCorrections(kFALSE,fgPtTPCCutOff,"Full");
 575       fCorrHadEtFullAcceptanceITS = fCorrections->GetConstantCorrections(kFALSE,fgPtITSCutOff,"Full");
 576       fCorrPiKPEtFullAcceptanceTPC = fCorrections->GetConstantCorrections(kFALSE,fgPtTPCCutOff,"PiKP");
 577       fCorrPiKPEtFullAcceptanceITS = fCorrections->GetConstantCorrections(kFALSE,fgPtITSCutOff,"PiKP");
 578       fCorrTotEtEMCALAcceptanceTPC = fCorrections->GetConstantCorrections(kTRUE,fgPtTPCCutOff,"EMCAL");
 579       fCorrTotEtEMCALAcceptanceITS = fCorrections->GetConstantCorrections(kTRUE,fgPtITSCutOff,"EMCAL");
 580       fCorrHadEtEMCALAcceptanceTPC = fCorrections->GetConstantCorrections(kFALSE,fgPtTPCCutOff,"EMCAL");
 581       fCorrHadEtEMCALAcceptanceITS = fCorrections->GetConstantCorrections(kFALSE,fgPtITSCutOff,"EMCAL");
 582       fCorrTotEtPHOSAcceptanceTPC = fCorrections->GetConstantCorrections(kTRUE,fgPtTPCCutOff,"PHOS");
 583       fCorrTotEtPHOSAcceptanceITS = fCorrections->GetConstantCorrections(kTRUE,fgPtITSCutOff,"PHOS");
 584       fCorrHadEtPHOSAcceptanceTPC = fCorrections->GetConstantCorrections(kFALSE,fgPtTPCCutOff,"PHOS");
 585       fCorrHadEtPHOSAcceptanceITS = fCorrections->GetConstantCorrections(kFALSE,fgPtITSCutOff,"PHOS");
 586     }
 587     else{cerr<<"Uh-oh!  Unable to open configuration file!"<<endl;}
 588   }
 589 }
 590 void AliAnalysisHadEtReconstructed::CreateHistograms(){//Creating histograms and adding them to the output TList
 591
 592   if(kDoTriggerChecksOnly){return;}//In this case we are just after trigger efficiencies and don't care about the ET reconstructed.
 593   //TString *strTPC = new TString("TPC");
 594   TString *strITS = new TString("ITS");
 595   TString *strTPCITS = new TString("TPCITS");
 596   for(Int_t i=0;i<2;i++){
 597     TString *cutName = NULL;
 598     Float_t maxPtdEdx = 10;
 599     Float_t mindEdx = 35;
 600     Float_t maxdEdx = 150.0;
 601     switch(i){
 602     case 0:
 603       cutName = strTPCITS;
 604       break;
 605     case 1:
 606       cutName = strITS;
 607       maxPtdEdx = 5;
 608       maxdEdx = 500.0;
 609       break;
 610       //not deleting this completely since we might want to add it back in later
 611       //     case 2:
 612       //       cutName = strTPC;
 613       //       maxPtdEdx = 5;
 614       //       maxdEdx = 500.0;
 615       //       break;
 616     default:
 617       cerr<<"Error:  cannot make histograms!"<<endl;
 618       return;
 619     }
 620
 621     CreateEtaPtHisto2D(Form("EtDataRaw%sPiPlus",cutName->Data()),"Raw reconstructed E_{T} from identified #pi^{+}");
 622     CreateEtaPtHisto2D(Form("EtDataRaw%sPiMinus",cutName->Data()),"Raw reconstructed E_{T} from identified #pi^{-}");
 623     CreateEtaPtHisto2D(Form("EtDataRaw%sKPlus",cutName->Data()),"Raw reconstructed E_{T} from identified K^{+}");
 624     //     CreateEtaPtHisto2D(Form("EtDataRaw%sEMinus",cutName->Data()),"Raw reconstructed E_{T} from identified e^{-}");
 625     //     CreateEtaPtHisto2D(Form("EtDataRaw%sEPlus",cutName->Data()),"Raw reconstructed E_{T} from identified e^{+}");
 626     CreateEtaPtHisto2D(Form("EtDataRaw%sKMinus",cutName->Data()),"Raw reconstructed E_{T} from identified K^{-}");
 627     CreateEtaPtHisto2D(Form("EtDataRaw%sProton",cutName->Data()),"Raw reconstructed E_{T} from identified p");
 628     CreateEtaPtHisto2D(Form("EtDataRaw%sAntiProton",cutName->Data()),"Raw reconstructed E_{T} from identified #bar{p}");
 629     CreateEtaPtHisto2D(Form("EtDataRaw%sUnidentified",cutName->Data()),"Raw reconstructed E_{T} from unidentified particles using real mass");
 630     CreateEtaPtHisto2D(Form("EtDataRaw%sNoID",cutName->Data()),"Raw reconstructed E_{T} from unidentified particles using real mass");
 631
 632     CreateEtaPtHisto2D(Form("EtDataCorrected%sPiPlus",cutName->Data()),"Corrected reconstructed E_{T} from identified #pi^{+}");
 633     CreateEtaPtHisto2D(Form("EtDataCorrected%sPiMinus",cutName->Data()),"Corrected reconstructed E_{T} from identified #pi^{-}");
 634     CreateEtaPtHisto2D(Form("EtDataCorrected%sKPlus",cutName->Data()),"Corrected reconstructed E_{T} from identified K^{+}");
 635     //     CreateEtaPtHisto2D(Form("EtDataCorrected%sEMinus",cutName->Data()),"Corrected reconstructed E_{T} from identified e^{-}");
 636     //     CreateEtaPtHisto2D(Form("EtDataCorrected%sEPlus",cutName->Data()),"Corrected reconstructed E_{T} from identified e^{+}");
 637     CreateEtaPtHisto2D(Form("EtDataCorrected%sKMinus",cutName->Data()),"Corrected reconstructed E_{T} from identified K^{-}");
 638     CreateEtaPtHisto2D(Form("EtDataCorrected%sProton",cutName->Data()),"Corrected reconstructed E_{T} from identified p");
 639     CreateEtaPtHisto2D(Form("EtDataCorrected%sAntiProton",cutName->Data()),"Corrected reconstructed E_{T} from identified #bar{p}");
 640     CreateEtaPtHisto2D(Form("EtDataCorrected%sUnidentified",cutName->Data()),"Corrected reconstructed E_{T} from unidentified particles using real mass");
 641     CreateEtaPtHisto2D(Form("EtDataCorrected%sNoID",cutName->Data()),"Corrected reconstructed E_{T} from unidentified particles using real mass");
 642
 643
 644     CreateEtaPtHisto2D(Form("EtNData%sPiPlus",cutName->Data()),"Number of reconstructed #pi^{+}");
 645     CreateEtaPtHisto2D(Form("EtNData%sPiMinus",cutName->Data()),"Number of reconstructed #pi^{-}");
 646     CreateEtaPtHisto2D(Form("EtNData%sKPlus",cutName->Data()),"Number of reconstructed K^{+}");
 647     CreateEtaPtHisto2D(Form("EtNData%sKMinus",cutName->Data()),"Number of reconstructed K^{-}");
 648     CreateEtaPtHisto2D(Form("EtNData%sProton",cutName->Data()),"Number of reconstructed p");
 649     CreateEtaPtHisto2D(Form("EtNData%sAntiProton",cutName->Data()),"Number of reconstructed #bar{p}");
 650     CreateEtaPtHisto2D(Form("EtNData%sUnidentified",cutName->Data()),"Number of Reconstructed unidentified particles");
 651
 652     CreateHisto2D(Form("dEdxDataAll%s",cutName->Data()),"dE/dx for all particles","momentum (GeV/c)","dE/dx",400,0.0,maxPtdEdx,200,mindEdx,maxdEdx);
 653     CreateHisto2D(Form("dEdxDataPion%s",cutName->Data()),"dE/dx for #pi^{#pm}","momentum (GeV/c)","dE/dx",400,0.0,maxPtdEdx,200,mindEdx,maxdEdx);
 654     CreateHisto2D(Form("dEdxDataKaon%s",cutName->Data()),"dE/dx for K^{#pm}","momentum (GeV/c)","dE/dx",400,0.0,maxPtdEdx,200,mindEdx,maxdEdx);
 655     CreateHisto2D(Form("dEdxDataProton%s",cutName->Data()),"dE/dx for p(#bar{p})","momentum (GeV/c)","dE/dx",400,0.0,maxPtdEdx,200,mindEdx,maxdEdx);
 656     CreateHisto2D(Form("dEdxDataElectron%s",cutName->Data()),"dE/dx for e^{#pm}","momentum (GeV/c)","dE/dx",400,0.0,maxPtdEdx,200,mindEdx,maxdEdx);
 657     CreateHisto2D(Form("dEdxDataUnidentified%s",cutName->Data()),"dE/dx for unidentified particles","momentum (GeV/c)","dE/dx",400,0.0,maxPtdEdx,200,mindEdx,maxdEdx);
 658   }
 659
 660   Float_t minEt = 0.0;
 661   Float_t maxEt = 100.0;
 662   Float_t minEtPiKP = 0.0;
 663   Float_t maxEtPiKP = 50.0;
 664   if(fDataSet==20100) maxEt=4000.0;
 665   Int_t nbinsEt = 200;
 666   char histoname[200];
 667   char histotitle[200];
 668   char xtitle[50];
 669   TString *ytitle = new TString("Number of events");
 670   TString *sTPC = new TString("TPC");
 671   TString *sITS = new TString("ITS");
 672   TString *sTPCpt = new TString("0.15");
 673   TString *sITSpt = new TString("0.10");
 674   TString *sPID = new TString("");
 675   TString *sNoPID = new TString("NoPID");
 676   TString *sNoPIDString = new TString(", No PID");
 677   TString *sHadEt = new TString("HadEt");
 678   TString *sPiKPEt = new TString("PiKPEt");
 679   TString *sRawEt = new TString("RawEt");
 680   TString *sTotEt = new TString("TotEt");
 681   TString *sTotEtString = new TString("total E_{T}");
 682   TString *sHadEtString = new TString("hadronic E_{T}");
 683   TString *sPiKPEtString = new TString("E_{T}^{#pi,K,p}");
 684   TString *sRawEtString = new TString("raw E_{T}");
 685   TString *sFull = new TString("Full");
 686   TString *sEMCAL = new TString("EMCAL");
 687   TString *sPHOS = new TString("PHOS");
 688
 689   for(int tpc = 0;tpc<2;tpc++){
 690     for(int hadet = 0;hadet<3;hadet++){
 691       for(int type = 0;type<3;type++){
 692         for(int pid = 0;pid<2;pid++){
 693           TString *detector = NULL;
 694           TString *partid = NULL;
 695           TString *et = sHadEt;
 696           TString *acceptance = NULL;
 697           TString *ptstring = NULL;
 698           TString *partidstring = NULL;
 699           TString *etstring = sHadEtString;
 700           if(tpc==1) {detector = sTPC; ptstring = sTPCpt;}
 701           else{detector = sITS; ptstring = sITSpt;}
 702           if(pid==1){partid = sPID; partidstring = sPID;}
 703           else{partid = sNoPID; partidstring = sNoPIDString;}
 704           if(hadet==1) {et = sHadEt; etstring = sHadEtString;}
 705           if(hadet==0){et = sTotEt; etstring = sTotEtString;}
 706           if(hadet==2){et = sRawEt; etstring = sRawEtString;}
 707           switch(type){
 708           case 0:
 709             acceptance = sFull;
 710             break;
 711           case 1:
 712             acceptance = sEMCAL;
 713             break;
 714           case 2:
 715             acceptance = sPHOS;
 716             break;
 717           default:
 718             acceptance = sFull;
 719           }
 720           snprintf(histoname,200,"Reco%s%sAcceptance%s%s",et->Data(),acceptance->Data(),detector->Data(),partid->Data());
 721           snprintf(histotitle,200,"Reconstructed %s with %s acceptance for p_{T}>%s GeV/c%s",etstring->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data());
 722           snprintf(xtitle,50,"Reconstructed %s",etstring->Data());
 723           CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEt,maxEt);
 724           if(type==0){//full acceptance only
 725             snprintf(xtitle,50,"Reconstructed %s",etstring->Data());
 726
 727             snprintf(histoname,200,"Reco%s%sAcceptance%s%sND",et->Data(),acceptance->Data(),detector->Data(),partid->Data());
 728             snprintf(histotitle,200,"Reconstructed non-diffractive events %s with %s acceptance for p_{T}>%s GeV/c%s",etstring->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data());
 729             CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEt,maxEt);
 730             snprintf(histoname,200,"Reco%s%sAcceptance%s%sDD",et->Data(),acceptance->Data(),detector->Data(),partid->Data());
 731             snprintf(histotitle,200,"Reconstructed doubly-diffractive events %s with %s acceptance for p_{T}>%s GeV/c%s",etstring->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data());
 732             CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEt,maxEt);
 733             snprintf(histoname,200,"Reco%s%sAcceptance%s%sSD",et->Data(),acceptance->Data(),detector->Data(),partid->Data());
 734             snprintf(histotitle,200,"Reconstructed singly-diffractive events %s with %s acceptance for p_{T}>%s GeV/c%s",etstring->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data());
 735             CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEt,maxEt);
 736
 737
 738
 739           snprintf(histoname,200,"Reco%s%sAcceptance%s%s",sPiKPEt->Data(),acceptance->Data(),detector->Data(),partid->Data());
 740           snprintf(histotitle,200,"Reconstructed %s with %s acceptance for p_{T}>%s GeV/c%s",sPiKPEtString->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data());
 741           snprintf(xtitle,50,"Reconstructed %s",sPiKPEtString->Data());
 742           CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEtPiKP,maxEtPiKP);
 743
 744             snprintf(xtitle,50,"Reconstructed %s",sPiKPEt->Data());
 745             snprintf(histoname,200,"Reco%s%sAcceptance%s%sND",sPiKPEt->Data(),acceptance->Data(),detector->Data(),partid->Data());
 746             snprintf(histotitle,200,"Reconstructed non-diffractive events %s with %s acceptance for p_{T}>%s GeV/c%s",sPiKPEtString->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data());
 747             CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEtPiKP,maxEtPiKP);
 748             snprintf(histoname,200,"Reco%s%sAcceptance%s%sDD",sPiKPEt->Data(),acceptance->Data(),detector->Data(),partid->Data());
 749             snprintf(histotitle,200,"Reconstructed doubly-diffractive events %s with %s acceptance for p_{T}>%s GeV/c%s",sPiKPEtString->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data());
 750             CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEtPiKP,maxEtPiKP);
 751             snprintf(histoname,200,"Reco%s%sAcceptance%s%sSD",sPiKPEt->Data(),acceptance->Data(),detector->Data(),partid->Data());
 752             snprintf(histotitle,200,"Reconstructed singly-diffractive events %s with %s acceptance for p_{T}>%s GeV/c%s",sPiKPEtString->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data());
 753             CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEtPiKP,maxEtPiKP);//et
 754
 755           }
 756           if(fDataSet==20100 && type ==0 &&pid==1){//If this is Pb+Pb and full acceptance with pid
 757             Int_t width = 5;
 758             if(fNCentBins<21) width = 10;
 759             for(Int_t i=0;i<fNCentBins;i++){
 760               snprintf(histoname,200,"Reco%s%sAcceptance%s%sCB%i",et->Data(),acceptance->Data(),detector->Data(),partid->Data(),i);
 761               snprintf(histotitle,200,"Reconstructed %s with %s acceptance for p_{T}>%s GeV/c%s for centrality %i-%i",etstring->Data(),acceptance->Data(),ptstring->Data(),partidstring->Data(),i*width,(i+1)*width);
 762               snprintf(xtitle,50,"Reconstructed %s",etstring->Data());
 763               CreateHisto1D(histoname,histotitle,xtitle,ytitle->Data(),nbinsEt*2,minEt,maxEt);
 764             }
 765           }
 766         }
 767       }
 768     }
 769   }
 770   CreateHisto1D("RecoHadEtFullAcceptanceTPCAssumingPion","Reconstructing E_{T}^{had} with full acceptance for p_{T}>0.15 GeV/c assuming pions","Reconstructed E_{T}^{had}","dN_{eve}/dE_{T}^{had}",nbinsEt*2,minEt,maxEt);
 771   CreateHisto1D("RecoHadEtFullAcceptanceTPCAssumingProton","Reconstructing E_{T}^{had} with full acceptance for p_{T}>0.15 GeV/c assuming protons","Reconstructed E_{T}^{had}","dN_{eve}/dE_{T}^{had}",nbinsEt*2,minEt,maxEt);
 772   CreateHisto1D("RecoHadEtFullAcceptanceTPCAssumingKaon","Reconstructing E_{T}^{had} with full acceptance for p_{T}>0.15 GeV/c assuming kaons","Reconstructed E_{T}^{had}","dN_{eve}/dE_{T}^{had}",nbinsEt*2,minEt,maxEt);
 773   CreateHisto1D("RecoHadEtFullAcceptanceITSAssumingPion","Reconstructing E_{T}^{had} with full acceptance for p_{T}>0.10 GeV/c assuming pions","Reconstructed E_{T}^{had}","dN_{eve}/dE_{T}^{had}",nbinsEt*2,minEt,maxEt);
 774   CreateHisto1D("RecoHadEtFullAcceptanceITSAssumingProton","Reconstructing E_{T}^{had} with full acceptance for p_{T}>0.10 GeV/c assuming protons","Reconstructed E_{T}^{had}","dN_{eve}/dE_{T}^{had}",nbinsEt*2,minEt,maxEt);
 775   CreateHisto1D("RecoHadEtFullAcceptanceITSAssumingKaon","Reconstructing E_{T}^{had} with full acceptance for p_{T}>0.10 GeV/c assuming kaons","Reconstructed E_{T}^{had}","dN_{eve}/dE_{T}^{had}",nbinsEt*2,minEt,maxEt);
 776
 777   //Cross checks that corrections are applied correctly
 778   if(fDataSet==20100){
 779     CreateHisto2D("fbkgdVsCentralityBin","f_{bkgd} vs centrality bin","centrality bin","f_{bkgd}",21,-1.5,19.5,200,0.7,1.05);//
 780     CreateHisto2D("feffPionVsCentralityBin","Pion efficiency vs centrality bin","centrality bin","pion efficiency",21,-1.5,19.5,200,0,1.2);//
 781     CreateHisto2D("feffHadronVsCentralityBin","Hadron efficiency vs centrality bin","centrality bin","hadron efficiency",21,-1.5,19.5,200,0,1.2);//
 782     CreateHisto2D("feffKaonVsCentralityBin","Kaon efficiency vs centrality bin","centrality bin","kaon efficiency",21,-1.5,19.5,200,0,1.2);//
 783     CreateHisto2D("feffProtonVsCentralityBin","Proton efficiency vs centrality bin","centrality bin","proton efficiency",21,-1.5,19.5,200,0,1.2);//
 784     CreateHisto2D("fnotIDVsCentralityBin","f_{notID} vs centrality bin","centrality bin","f_{notID}",21,-1.5,19.5,50,0.95,1.05);//
 785     CreateHisto2D("fpTcutVsCentralityBin","f_{pTcut} vs centrality bin","centrality bin","f_{pTcut}",21,-1.5,19.5,50,0.95,1.05);
 786     CreateHisto2D("fneutralVsCentralityBin","f_{neutral} vs centrality bin","centrality bin","f_{neutral}",21,-1.5,19.5,50,0.5,1.00);
 787     CreateHisto2D("ConstantCorrectionsVsCentralityBin","constant corrections vs centrality bin","centrality bin","constant corrections",21,-1.5,19.5,50,0.5,1.00);
 788   }
 789
 790   delete sTPC;
 791   delete sITS;
 792   delete sTPCpt;
 793   delete sITSpt;
 794   delete sPID;
 795   delete sNoPID;
 796   delete sNoPIDString;
 797   delete sPiKPEt;
 798   delete sHadEt;
 799   delete sTotEt;
 800   delete sTotEtString;
 801   delete sHadEtString;
 802   delete sPiKPEtString;
 803   delete sFull;
 804   delete sEMCAL;
 805   delete sPHOS;
 806
 807 }