From acd6f996c79a0fd5e896300d11e133e6e13f145c Mon Sep 17 00:00:00 2001 From: ssakai Date: Thu, 19 Sep 2013 12:25:24 +0000 Subject: [PATCH] upated for flattering --- PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx | 39 ++++++++++++------- PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h | 10 ++--- PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C | 7 +--- 3 files changed, 32 insertions(+), 24 deletions(-) diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx index 809cc2b7357..b23975bb35e 100644 --- a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx +++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx @@ -186,6 +186,8 @@ AliAnalysisTaskFlowTPCEMCalQCSP::AliAnalysisTaskFlowTPCEMCalQCSP(const char *nam ,fMultCorAfterVZTRKComp(0) ,fCentralityBeforePileup(0) ,fCentralityAfterVZTRK(0) +,fCentralityAfterCorrCut(0) +,fMultCorAfterCorrCut(0) ,EPVz(0) ,EPTPCp(0) ,EPTPCn(0) @@ -292,6 +294,8 @@ AliAnalysisTaskFlowTPCEMCalQCSP::AliAnalysisTaskFlowTPCEMCalQCSP() ,fMultCorAfterVZTRKComp(0) ,fCentralityBeforePileup(0) ,fCentralityAfterVZTRK(0) +,fCentralityAfterCorrCut(0) +,fMultCorAfterCorrCut(0) ,EPVz(0) ,EPTPCp(0) ,EPTPCn(0) @@ -379,7 +383,10 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserExec(Option_t*) if(fTrigger==3){ if(!(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kMB) ) return; } - + if(fTrigger==4){ + if(!(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & (AliVEvent::kCentral | AliVEvent::kSemiCentral))) return; + } + //---------------CENTRALITY AND EVENT SELECTION----------------------- Int_t fNOtrks = fAOD->GetNumberOfTracks(); @@ -775,11 +782,11 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::SelectPhotonicElectron(Int_t itrack,const Bool_t fFlagLS=kFALSE, fFlagULS=kFALSE; Double_t openingAngle = -999.; Double_t ptcutonmasshighpt = track->Pt(); - + ptAsso = trackAsso->Pt(); Short_t chargeAsso = trackAsso->Charge(); Short_t charge = track->Charge(); - + nsigma = fPID->GetPIDResponse() ? fPID->GetPIDResponse()->NumberOfSigmasTPC(trackAsso, AliPID::kElectron) : 1000; //80 @@ -815,8 +822,8 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::SelectPhotonicElectron(Int_t itrack,const if(fFlagULS) fInvmassULS1->Fill(mass); if(ptcutonmasshighpt >= 8.){ - if(fFlagLS) fInvmassLS1highpt->Fill(mass); - if(fFlagULS) fInvmassULS1highpt->Fill(mass); + if(fFlagLS) fInvmassLS1highpt->Fill(mass); + if(fFlagULS) fInvmassULS1highpt->Fill(mass); } @@ -928,7 +935,7 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserCreateOutputObjects() fInvmassULS1 = new TH1F("fInvmassULS1", "Inv mass of ULS (e,e); mass(GeV/c^2); counts;", 1000,0,1.0); fOutputList->Add(fInvmassULS1); - + fInvmassLS1highpt = new TH1F("fInvmassLS1highpt", "Inv mass of LS (e,e); mass(GeV/c^2) highpt; counts;", 1000,0,1.0); fOutputList->Add(fInvmassLS1highpt); @@ -949,6 +956,9 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserCreateOutputObjects() fCentralityAfterVZTRK = new TH1F("fCentralityAfterVZTRK", "fCentralityAfterVZTRK Pass", 101, -1, 100); fOutputList->Add(fCentralityAfterVZTRK); + + fCentralityAfterCorrCut = new TH1F("fCentralityAfterCorrCut", "fCentralityAfterCorrCut Pass", 101, -1, 100); + fOutputList->Add(fCentralityAfterCorrCut); fPhi = new TH1F("fPhi", "#phi distribution", 100, -.5, 7); fOutputList->Add(fPhi); @@ -980,6 +990,9 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserCreateOutputObjects() fMultCorAfterVZTRKComp = new TH2F("fMultCorAfterVZTRKComp", "TPC vs Global multiplicity (After V0-TRK); Global multiplicity; TPC multiplicity", 100, 0, 3000, 100, 0, 3000); fOutputList->Add(fMultCorAfterVZTRKComp); + fMultCorAfterCorrCut = new TH2F("fMultCorAfterCorrCut", "TPC vs Global multiplicity (After CorrCut); Global multiplicity; TPC multiplicity", 100, 0, 3000, 100, 0, 3000); + fOutputList->Add(fMultCorAfterCorrCut); + fMultvsCentr = new TH2F("fMultvsCentr", "Multiplicity vs centrality; centrality; Multiplicity", 100, 0., 100, 100, 0, 3000); fOutputList->Add(fMultvsCentr); @@ -1056,7 +1069,7 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserCreateOutputObjects() //---------------------------------------------------------------------------- - + if(fPhiminusPsi){ Int_t binsvElectH[8]={ 600, 200, 200 ,100, 100, 100, 10, 10}; //pt, E/p,TPCnSigma,M20,M02,Disp Phi-psiV0A ,Phi-PsiV0C,eta (commented) Double_t xminvElectH[8]={0, 0, -10 , 0, 0, 0, 0, 0}; @@ -1198,7 +1211,7 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::CheckCentrality(AliAODEvent* event, Bool_t }//...after centrality selectrion //============================================================================================================================ if(fMultCut){ - if(fTrigger==1){ + if(fTrigger==1 || fTrigger==4){ if(! (multTPC > (-36.73 + 1.48*multGlob) && multTPC < (62.87 + 1.78*multGlob))){ // cout <<" Trigger ==" <Fill(multGlob, multTPC); + fCentralityAfterCorrCut->Fill(fCentrality); + }//...after CORR CUT //=================================All cuts are passed==================++++================================================== //=================================Now Centrality flattening for central trigger==================++++================================================== - if(fTrigger==0){ + if(fTrigger==0 || fTrigger==4){ if(!IsEventSelectedForCentrFlattening(fCentrality)){ centralitypass = kFALSE; fCentralityNoPassForFlattening->Fill(fCentrality); @@ -1338,7 +1355,3 @@ Bool_t AliAnalysisTaskFlowTPCEMCalQCSP::IsEventSelectedForCentrFlattening(Float_ } //--------------------------------------------------------------------------- - - - - diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h index 5507bb9a556..5408bd92a73 100644 --- a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h +++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h @@ -1,5 +1,3 @@ - - #ifndef ALIANALYSISTASKFLOWTPCEMCALQCSP_H #define ALIANALYSISTASKFLOWTPCEMCALQCSP_H @@ -78,7 +76,7 @@ public: void SetHistoForCentralityFlattening(TH1F *h,Double_t minCentr,Double_t maxCentr,Double_t centrRef=0.,Int_t switchTRand=0); Bool_t IsEventSelectedForCentrFlattening(Float_t centvalue); - + @@ -164,6 +162,8 @@ private: TH2F *fMultCorAfterVZTRKComp;//!after cent comp TH1F *fCentralityBeforePileup;//!cent chneck TH1F *fCentralityAfterVZTRK;//!cent chneck2 + TH1F *fCentralityAfterCorrCut;//!cent chneck2 + TH2F *fMultCorAfterCorrCut;//!after cent comp TH1D *EPVz;//v0cep TH1D *EPTPCp;//tpcep TH1D *EPTPCn;//!tpcep @@ -173,7 +173,7 @@ private: TH1F *fCentralityNoPassForFlattening; //! QA histogram of events that do not pass centrality cut for flattening TH1F *fInvmassLS1highpt; //LS Invmass for all rec par high pt TH1F *fInvmassULS1highpt;//ULS Invmass for all rec par high pt - + AliAnalysisTaskFlowTPCEMCalQCSP(const AliAnalysisTaskFlowTPCEMCalQCSP&); // not implemented @@ -184,5 +184,3 @@ private: #endif - - diff --git a/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C b/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C index cd43cda1d70..586ab62cdfc 100644 --- a/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C +++ b/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C @@ -1,5 +1,3 @@ - - /////////////////////////////////////////////////////////////////// // // // AddTaskFlowTPCEMCalQCSP macro // @@ -83,7 +81,7 @@ AliAnalysisTaskFlowTPCEMCalQCSP* AddTaskFlowTPCEMCalQCSP( } taskHFE->SetTrigger(Trigger); - if(Trigger==0){ + if(Trigger==0 || Trigger==4){ TFile *fFlat=TFile::Open(histoflatname.Data()); TCanvas *c=fFlat->Get("cintegral"); TH1F *hfl=(TH1F*)c->FindObject("hint"); @@ -370,5 +368,4 @@ AliAnalysisTaskFlowTPCEMCalQCSP* ConfigHFEemcalMod(Bool_t useMC,Int_t minTPCCuls } -//_____________________________________________________________________________ - +//_____________________________________________________________________________ \ No newline at end of file -- 2.43.0