including QA histogram (E/p vs. E) to survey material in front of the EMCAL
authormcosenti <mcosenti@cern.ch>
Wed, 24 Sep 2014 17:03:48 +0000 (14:03 -0300)
committermcosenti <mcosenti@cern.ch>
Wed, 24 Sep 2014 17:04:48 +0000 (14:04 -0300)
PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h

index 7e5ee2f..b969dc8 100644 (file)
@@ -147,7 +147,8 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton() :
   fTrackPtEta(0),     
   fTrackPtEtaCut(0),
   fMaxCellEPhi(0),
-  fDetaDphiFromTM(0)
+  fDetaDphiFromTM(0),
+  fEoverPvsE(0)
 {
   // Default constructor.
   for(Int_t i = 0; i < 12;    i++)  fGeomMatrix[i] =  0;
@@ -258,7 +259,8 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton(const char *name) :
   fTrackPtEta(0),     
   fTrackPtEtaCut(0),   
   fMaxCellEPhi(0),
-  fDetaDphiFromTM(0)
+  fDetaDphiFromTM(0),
+  fEoverPvsE(0)
 {
   // Constructor
 
@@ -456,6 +458,10 @@ void AliAnalysisTaskEMCALIsoPhoton::UserCreateOutputObjects()
   fDetaDphiFromTM->Sumw2();
   fQAList->Add(fDetaDphiFromTM);
 
+  fEoverPvsE = new TH2F("fEoverPvsE","E^{clus}/p^{track} vs E^{clus} (80<TPCsignal<100);E^{clus} [GeV];E^{clus}/p^{track} [c^{-1}]",fNBinsPt, fPtBinLowEdge,fPtBinHighEdge,100,0,2);
+  fEoverPvsE->Sumw2();
+  fQAList->Add(fEoverPvsE);
+
   PostData(1, fOutputList);
   PostData(2, fQAList);
 }
@@ -1482,6 +1488,14 @@ void AliAnalysisTaskEMCALIsoPhoton::FillQA()
       fTrackPtPhiCut->Fill(t->Pt(), t->Phi());
       fTrackPtEtaCut->Fill(t->Pt(), t->Eta());
     }
+    if(t->GetTPCsignal()<80 || t->GetTPCsignal()>100)
+      continue;
+    if(t->GetEMCALcluster()<=0 || t->GetEMCALcluster()>nclus)
+      continue;
+    AliVCluster *c = dynamic_cast<AliVCluster*>(clusters->At(t->GetEMCALcluster()));
+    if(!c)
+      continue;
+    fEoverPvsE->Fill(c->E(),c->E()/t->P());
   }
   for(int ic=0;ic<nclus;ic++){
     AliVCluster *c = dynamic_cast<AliVCluster*>(clusters->At(ic));
index 5e8074a..d69014e 100644 (file)
@@ -188,6 +188,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE {
   TH2F        *fTrackPtEtaCut;    //!selected tracks pt vs. eta in "triggered event"
   TH2F        *fMaxCellEPhi;      //!max cell energy vs. cell phi
   TH2F        *fDetaDphiFromTM;   //!dphi vs deta of track->GetEMCALcluster() clusters
+  TH2F        *fEoverPvsE;        //!E/p for tracks with 80<TPCsignal<100 vs cluster E (check material)
 
 
   AliAnalysisTaskEMCALIsoPhoton(const AliAnalysisTaskEMCALIsoPhoton&); // not implemented