Added histogram for ITS-PID cluster selection

diff --git a/PWG1/ITS/AliAnalysisTaskITSTrackingCheck.cxx b/PWG1/ITS/AliAnalysisTaskITSTrackingCheck.cxx
index 874fd5062a732481b88ef431bdb6ba85c18747f4..8a133a8e10689a08ce0a8fbb5596abae338b085d 100644 (file)
--- a/PWG1/ITS/AliAnalysisTaskITSTrackingCheck.cxx
+++ b/PWG1/ITS/AliAnalysisTaskITSTrackingCheck.cxx
@@ -150,6 +150,7 @@ fHistPtITSMI5InAcc(0),
 fHistPtITSMI6InAcc(0),
 fHistPtITSMISPDInAcc(0),
 fHistPtITSMIoneSPDInAcc(0),
+fHistPtITSMIoneSPDthreeSDDSSDInAcc(0),
 fHistPtITSTPCsel(0),
 fHistPtITSTPCselP(0),
 fHistPtITSTPCselS(0),
@@ -332,6 +333,7 @@ fHistPtITSMI5InAcc(0),
 fHistPtITSMI6InAcc(0),
 fHistPtITSMISPDInAcc(0),
 fHistPtITSMIoneSPDInAcc(0),
+fHistPtITSMIoneSPDthreeSDDSSDInAcc(0),
 fHistPtITSTPCsel(0),
 fHistPtITSTPCselP(0),
 fHistPtITSTPCselS(0),
@@ -832,6 +834,11 @@ void AliAnalysisTaskITSTrackingCheck::UserCreateOutputObjects()
   fHistPtITSMIoneSPDInAcc->SetMinimum(0);
   fOutput->Add(fHistPtITSMIoneSPDInAcc);
 
+  fHistPtITSMIoneSPDthreeSDDSSDInAcc = new TH1F("fHistPtITSMIoneSPDthreeSDDSSDInAcc","pt distribution of ITSMI tracks (>0 in SPD, >2 in SDD+SSD); p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
+  fHistPtITSMIoneSPDthreeSDDSSDInAcc->Sumw2();
+  fHistPtITSMIoneSPDthreeSDDSSDInAcc->SetMinimum(0);
+  fOutput->Add(fHistPtITSMIoneSPDthreeSDDSSDInAcc);
+
   fHistPtITSTPCsel = new TH1F("fHistPtITSTPCsel","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
   fHistPtITSTPCsel->Sumw2();
   fHistPtITSTPCsel->SetMinimum(0);
@@ -1489,11 +1496,13 @@ void AliAnalysisTaskITSTrackingCheck::UserExec(Option_t *)
    
 
     Int_t nclsITS = track->GetNcls(0);
+    Int_t nclsSDDSSD = 0;
     Int_t nclsokbadoutinzITS = 0;
     Bool_t outInZ=kFALSE;
     Bool_t skipTrack=kFALSE;
 
     for(Int_t layer=0; layer<6; layer++) {
+      if(layer>=2 && track->HasPointOnITSLayer(layer)) nclsSDDSSD++;
       if(layer==0 && !track->HasPointOnITSLayer(1)) continue;
       if(layer==1 && !track->HasPointOnITSLayer(0)) continue;
       track->GetITSModuleIndexInfo(layer,idet,status,xloc,zloc);
@@ -1751,6 +1760,8 @@ void AliAnalysisTaskITSTrackingCheck::UserExec(Option_t *)
        }
        if(track->HasPointOnITSLayer(0) || track->HasPointOnITSLayer(1)) {
          fHistPtITSMIoneSPDInAcc->Fill(track->Pt());
+         if(nclsSDDSSD>=3) fHistPtITSMIoneSPDthreeSDDSSDInAcc->Fill(track->Pt());
+
          if(isPrimary) {fHistPtITSMIoneSPDInAccP->Fill(track->Pt());} else {fHistPtITSMIoneSPDInAccS->Fill(track->Pt());}  
        }
        if(nclsokbadoutinzITS==6) fHistPtITSMIokbadoutinz6InAcc->Fill(track->Pt());

diff --git a/PWG1/ITS/AliAnalysisTaskITSTrackingCheck.h b/PWG1/ITS/AliAnalysisTaskITSTrackingCheck.h
index 1a34ac4f6e7bf6f26c94df10b1d9a669f21dd4a0..65fbbe3a6f11c0d9c9c033b5ed96e2388f5a38ea 100644 (file)
--- a/PWG1/ITS/AliAnalysisTaskITSTrackingCheck.h
+++ b/PWG1/ITS/AliAnalysisTaskITSTrackingCheck.h
@@ -130,6 +130,7 @@ class AliAnalysisTaskITSTrackingCheck : public AliAnalysisTaskSE
   TH1F         *fHistPtITSMI6InAcc; //! output hist
   TH1F         *fHistPtITSMISPDInAcc; //! output hist
   TH1F         *fHistPtITSMIoneSPDInAcc; //! output hist
+  TH1F         *fHistPtITSMIoneSPDthreeSDDSSDInAcc; //! output hist
   TH1F         *fHistPtITSTPCsel; //! output hist
   TH1F         *fHistPtITSTPCselP; //! output hist
   TH1F         *fHistPtITSTPCselS; //! output hist