Small update
authorrbailhac <rbailhac@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 11 Nov 2013 14:04:08 +0000 (14:04 +0000)
committerrbailhac <rbailhac@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 11 Nov 2013 14:04:08 +0000 (14:04 +0000)
PWGHF/hfe/AliAnalysisTaskHFEQA.cxx
PWGHF/hfe/AliHFEpidTRD.cxx
PWGHF/hfe/macros/AddTaskHFEQA.C
PWGHF/hfe/macros/configs/PbPb/ConfigHFEQA.C

index 8fe52c5..0a33d7c 100644 (file)
@@ -343,6 +343,8 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
 \r
   AliMCEvent *mcEvent = MCEvent();\r
  \r
+\r
+  AliDebug(2,"MC info");\r
   // MC info\r
   Bool_t mcthere = kTRUE;\r
   if(fAODAnalysis) {\r
@@ -375,7 +377,7 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
   ////////////////////////////////////\r
   // Number of contributors\r
   ///////////////////////////////////\r
-   \r
+  AliDebug(2,"Number of contributors");\r
   Int_t ncontribVtx = 0;\r
   if(fAODAnalysis) {\r
     AliAODEvent *fAOD = dynamic_cast<AliAODEvent *>(fInputEvent);\r
@@ -399,6 +401,7 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
       ncontribVtx = priVtx->GetNContributors();\r
     }\r
   }\r
+  AliDebug(2,Form("Number of contributors %d",ncontribVtx));\r
 \r
 \r
   /////////////////////////////////\r
@@ -424,6 +427,7 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
     if((90.0< cntr) && (cntr<100.0)) binct = 10.5;\r
     if(binct > 11.0) return;\r
   }\r
+  AliDebug(2,Form("Centrality %f with %s",binct,fCentralityEstimator.Data()));\r
  \r
   //////////////////////\r
   // run number\r
@@ -453,7 +457,7 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
   //////////\r
   // PID\r
   //////////\r
\r
+  AliDebug(2,"PID response");\r
   AliPIDResponse *pidResponse = fInputHandler->GetPIDResponse();\r
   if(!pidResponse){\r
     AliDebug(2,"No PID response set");\r
@@ -467,7 +471,7 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
   //////////////////\r
   // Event cut\r
   //////////////////\r
-  //printf("Event cut\n");\r
+  AliDebug(2,"Event cut");\r
   if(!fHFECuts->CheckEventCuts("fEvRecCuts", fInputEvent)) {\r
     AliDebug(2,"Does not pass the event cut");\r
     PostData(1, fListHist);\r
@@ -478,41 +482,40 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
   // Loop over track\r
   //////////////////////////\r
   Int_t nbtracks = fInputEvent->GetNumberOfTracks();\r
-  //printf("There are %d tracks \n",nbtracks);\r
+  AliDebug(2,Form("Number of tracks %d",nbtracks));\r
   for(Int_t k = 0; k < nbtracks; k++){\r
       \r
     AliVTrack *track = (AliVTrack *) fInputEvent->GetTrack(k);\r
     if(!track) continue;\r
     Double_t pt = track->Pt();     \r
 \r
-    //printf("test 0\n");\r
+    AliDebug(2,"test 0\n");\r
     \r
     // RecKine: ITSTPC cuts  \r
     if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepRecKineITSTPC + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-    //printf("test 1\n");\r
+    AliDebug(2,"test 1\n");\r
 \r
     // RecPrim\r
     if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepRecPrim + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-    //printf("test 2\n");\r
+    AliDebug(2,"test 2\n");\r
 \r
     // HFEcuts: ITS layers cuts\r
     if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepHFEcutsITS + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-    //printf("test 3\n");\r
+    AliDebug(2,"test 3\n");\r
 \r
     // HFE cuts: TOF and mismatch flag\r
     if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepHFEcutsTOF + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-    //printf("test 4\n");\r
+    AliDebug(2,"test 4\n");\r
 \r
     // HFE cuts: TPC PID cleanup\r
     if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepHFEcutsTPC + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-    //printf("test 5\n");\r
+    AliDebug(2,"test 5\n");\r
 \r
     // HFEcuts: Nb of tracklets TRD0\r
     if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepHFEcutsTRD + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-    //printf("test 6\n");\r
-\r
+    \r
     AliDebug(2,"Survived");\r
-    //printf("Survived track cuts\n");\r
+    \r
 \r
     ////////////////////////\r
     // Apply PID\r
@@ -539,16 +542,20 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
     if(fPIDTPConly->IsSelected(&hfetrack,0x0,"recTrackCont",0x0)) {\r
       fTPConly->Fill(pt);\r
     }\r
+    AliDebug(2,"TPC only PID\n");\r
+           \r
     \r
     // Complete PID TPC TOF \r
     if(fPIDTOFTPC->IsSelected(&hfetrack,0x0,"recTrackCont",fPIDqaTOFTPC)) {\r
       fTOFTPC->Fill(pt);\r
+      AliDebug(2,"TOF TPC PID\n");\r
     }\r
     \r
     // Complete PID TPC TRD \r
     if(fPIDTPCTRD->IsSelected(&hfetrack,0x0,"recTrackCont",fPIDqaTPCTRD)) {\r
       fTPCTRD->Fill(pt);\r
     }\r
+    AliDebug(2,"TPC TRD PID\n");\r
 \r
 \r
 \r
@@ -558,6 +565,7 @@ void AliAnalysisTaskHFEQA::UserExec(Option_t */*option*/)
        fTPCEMCal->Fill(pt);\r
       }\r
     }\r
+    AliDebug(2,"TPC EMCal PID\n");\r
     \r
     \r
   }\r
index bca6a60..6278eba 100644 (file)
@@ -297,6 +297,7 @@ Int_t AliHFEpidTRD::IsSelected2D(const AliHFEpidObject *track, AliHFEpidQAmanage
 
   if(fkPIDResponse->IdentifiedAsElectronTRD(track->GetRecTrack(),fElectronEfficiency,centrality,AliTRDPIDResponse::kLQ2D)){
       AliDebug(2, Form("Electron effi: %f %i %i %f\n", fElectronEfficiency,track->GetCentrality(),centralitybin,centrality));
+      if(pidqa) pidqa->ProcessTrack(track, AliHFEpid::kTRDpid, AliHFEdetPIDqa::kAfterPID); 
       return 11;
   } else return 211;
 
index a8acd61..9c478f7 100644 (file)
@@ -1,4 +1,4 @@
-AliAnalysisTask *AddTaskHFEQA(Bool_t useMC, Bool_t isAOD, Int_t icollisionsystem = 2, Int_t icent = 2){
+AliAnalysisTask *AddTaskHFEQA(Bool_t useMC, Bool_t isAOD, Int_t icollisionsystem = 2, Int_t icent = 2,Int_t debuglevel = 4){
 
    // Name
   TString appendixx("HFEQA");
@@ -14,6 +14,7 @@ AliAnalysisTask *AddTaskHFEQA(Bool_t useMC, Bool_t isAOD, Int_t icollisionsystem
   AliAnalysisTaskHFEQA *task = ConfigHFEQA(useMC,isAOD,icollisionsystem,icent);  
 
   mgr->AddTask(task);
+  mgr->AddClassDebug("AliAnalysisTaskHFEQA",debuglevel);
 
   TString containerName = mgr->GetCommonFileName();
   containerName += ":";
index 5862e6f..909aa5e 100644 (file)
@@ -96,8 +96,8 @@ AliAnalysisTaskHFEQA* ConfigHFEQA(Bool_t useMC, Bool_t isAOD, Int_t icollisionsy
   AliHFEpidTRD *trdpid = pidTPCTRD->GetDetPID(AliHFEpid::kTRDpid);
   trdpid->SetTRD2DPID();
   trdpid->SetElectronEfficiency(0.80);   // efficiency
-  trdpid->SetNTracklets(5);      // ntracklets threshold
-  trdpid->SetCutNTracklets(5, kTRUE);
+  trdpid->SetNTracklets(6);      // ntracklets threshold
+  trdpid->SetCutNTracklets(6, kTRUE);
 
   // change E/p cuts
   AliHFEpidEMCAL *emcpid = pidTPCEMCal->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);