upated for flattering

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
index 809cc2b7357087ddbcc7f6397fb8ef9e2928f438..b23975bb35e7c5e97f083826b52d9afda67cee4e 100644 (file)
--- 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 ==" <<fTrigger<< endl;
                 centralitypass = kFALSE;
@@ -1213,9 +1226,13 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::CheckCentrality(AliAODEvent* event, Bool_t
             }//2011
         }//2011 Central
     }
+    if (centralitypass){
+        fMultCorAfterCorrCut->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 5507bb9a556b7223f0dd0a9339bc9fe78a74008b..5408bd92a73a8633f2431432488ebcc7e0d4290d 100644 (file)
--- 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 cd43cda1d70ae5971b67fee8d35f1456242eaefb..586ab62cdfc171d8fa7733e56b48fa1d71159e73 100644 (file)
--- a/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C
+++ b/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C
@@ -1,5 +1,3 @@
-\r
-\r
 ///////////////////////////////////////////////////////////////////\r
 //                                                               //\r
 // AddTaskFlowTPCEMCalQCSP macro                                 //\r
@@ -83,7 +81,7 @@ AliAnalysisTaskFlowTPCEMCalQCSP*  AddTaskFlowTPCEMCalQCSP(
     }\r
     taskHFE->SetTrigger(Trigger);\r
     \r
-    if(Trigger==0){\r
+    if(Trigger==0 || Trigger==4){\r
         TFile *fFlat=TFile::Open(histoflatname.Data());\r
         TCanvas *c=fFlat->Get("cintegral");\r
         TH1F *hfl=(TH1F*)c->FindObject("hint");\r
@@ -370,5 +368,4 @@ AliAnalysisTaskFlowTPCEMCalQCSP* ConfigHFEemcalMod(Bool_t useMC,Int_t minTPCCuls
     \r
 }\r
 \r
-//_____________________________________________________________________________\r
-\r
+//_____________________________________________________________________________
\ No newline at end of file