Coding violations corrected.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 7 Jul 2010 14:36:26 +0000 (14:36 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 7 Jul 2010 14:36:26 +0000 (14:36 +0000)
PWG1/AliAnalysisTaskQASym.cxx
PWG1/AliAnalysisTaskQASym.h

index 904e33f..832061e 100644 (file)
@@ -104,8 +104,6 @@ ClassImp(AliAnalysisTaskQASym)
     ,fVertexX(0)
     ,fVertexY(0)
     ,fVertexZ(0)
-    ,test(0)
-  
     ,fRecDcaPosPhi(0)
     ,fRecDcaNegPhi(0)
     ,fRecPtPosPhi(0)
@@ -155,13 +153,6 @@ ClassImp(AliAnalysisTaskQASym)
     ,fEtavPt(0)  
     ,fCompareTPCparam(0)
     ,fITSlayer(0)
-    
-    ,sdca(0)
-    ,xy(0)
-    ,z(0)
-    ,xvertexcor(0)
-    ,yvertexcor(0)  
     ,fCuts(0)
 
 {
@@ -211,8 +202,7 @@ void AliAnalysisTaskQASym::UserCreateOutputObjects()
   Double_t pt = 20.;
 
   fHists = new TList();
-  //  test   = new TNtuple("test","test",  
-  //                     "pt:phi:theta:x:y:z:charge");
+
   fHistRECpt   = new TH1F("fHistRECpt", 
                          " p_{T}",
                          200, 0., pt);
@@ -1116,7 +1106,8 @@ void AliAnalysisTaskQASym::UserExec(Option_t *)
     sigmapt= sqrt(sigmapt);
     sigmapt= sigmapt *(tpcP->Pt()*tpcP->Pt()); 
 
-    if(sigmapt == 0.)continue;
+    if(TMath::Abs(sigmapt) < 1.e-10) continue;
+
     fsigmaPt->Fill(TMath::Log10(sigmapt));
  
 
@@ -1155,10 +1146,12 @@ void AliAnalysisTaskQASym::UserExec(Option_t *)
    
    
     //------------------- 
+    Float_t xvertexcor = 0.;
+    Float_t yvertexcor = 0.;
 
     xvertexcor = tpcP->Xv() - vertex->GetX(); // coordinate corrected for vertex position
     yvertexcor = tpcP->Yv() - vertex->GetY(); // "
-    sdca = (tpcP->Py()*xvertexcor - tpcP->Px()*yvertexcor)/tpcP->Pt();
+    Double_t sdca = (tpcP->Py()*xvertexcor - tpcP->Px()*yvertexcor)/tpcP->Pt();
 
 
     fqPtRec[cas]->Fill(tpcP->Charge()/tpcP->Pt());
@@ -1181,6 +1174,8 @@ void AliAnalysisTaskQASym::UserExec(Option_t *)
     fDca->Fill(sdca);
     fRecQPtPhi->Fill(tpcP->Charge()/tpcP->Pt(), phiIn);
 
+    Float_t xy = 0.;
+    Float_t  z = 0.;
 
     tpcP->GetImpactParameters(xy,z);
     fDiffDcaD->Fill(sdca+xy);
index 92c4f7e..9757ed8 100644 (file)
@@ -1,12 +1,11 @@
-#ifndef AliAnalysisTaskQASym_cxx
-#define AliAnalysisTaskQASym_cxx
+#ifndef ALIANALYSISTASKQASYM_H
+#define ALIANALYSISTASKQASYM_H
  
 
 class TH1F;
 class TH2F;
 class TH3F;
 class TList;
-class TNtuple;
 
 class AliESDEvent;
 class AliESDtrack;
@@ -15,7 +14,6 @@ class AliESDtrackCuts;
 
 #include "AliAnalysisTaskSE.h"
 #include "TFile.h"
-#include "TNtuple.h"
 
 class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
  public:
@@ -37,9 +35,9 @@ class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
   
  private:
 
-  Int_t       fTrackType;
-  Int_t       fStandAlone; // needed for ITS tracks
-  Bool_t      fFieldOn;
+  Int_t       fTrackType;       // track type
+  Int_t       fStandAlone;      // needed for ITS tracks
+  Bool_t      fFieldOn;         // field status
 
   TList       *fHists;          // List of histos
 
@@ -82,7 +80,7 @@ class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
 
   // 2D: all measures as function of z of first trackpoint 
   TH2F        *fRecPtPosVz;     //pt-zfirst of pos tracks
-  TH2F        *fRecPtNegVz;      //pt-zfirst of neg tracks
+  TH2F        *fRecPtNegVz;     //pt-zfirst of neg tracks
   TH2F        *fRecEtaPosVz;    //eta-zfirst of pos tracks
   TH2F        *fRecEtaNegVz;    //eta-zfirst of neg tracks
   TH2F        *fRecPhiPosVz;    //phi-zfirst of pos tracks
@@ -110,7 +108,6 @@ class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
   TH1F * fVertexX;             // x of vertex
   TH1F * fVertexY;             // y of vertex
   TH1F * fVertexZ;             // z of vertex
-  TNtuple * test;
 
   //new
   TH2F        *fRecDcaPosPhi;     //dca-phi for pos.
@@ -158,13 +155,11 @@ class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
   //  TH3F        *fRecDcaPhiPtPosEtaNeg; //dca-pt-phi for pos tracks and neg eta
   //  TH3F        *fRecDcaPhiPtNegEtaNeg; //dca-pt-phi for neg tracks and neg eta
 
-  TH2F        *fEtavPt;
-  TH2F        *fCompareTPCparam;
-  TH1F        *fITSlayer;
+  TH2F        *fEtavPt;                 // eta vs pt 
+  TH2F        *fCompareTPCparam;        // TPC param
+  TH1F        *fITSlayer;               // ITS layer
 
-  Double_t  sdca;
-  Float_t xy, z, xvertexcor, yvertexcor;
-  AliESDtrackCuts* fCuts;                      // List of cuts
+  AliESDtrackCuts* fCuts;               // List of cuts
 
   // sectors of TPC 
   TH1F        *fRecPtTpcSector[18];     //pt for TPC sectors