Added histo to study matching eff vs eta and phi
authordainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 12 Nov 2010 15:55:54 +0000 (15:55 +0000)
committerdainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 12 Nov 2010 15:55:54 +0000 (15:55 +0000)
PWG1/ITS/AliAnalysisTaskITSTrackingCheck.cxx
PWG1/ITS/AliAnalysisTaskITSTrackingCheck.h

index 78492b5..d7f83c3 100644 (file)
@@ -121,6 +121,7 @@ fHistxlocVSmodSDDok(0),
 fHistxlocSDDall(0),
 fHistzlocSDDall(0),
 fHistPhiTPCInAcc(0),
+fHistEtaTPCInAcc(0),
 fHistPtTPC(0),
 fHistPtTPCInAcc(0),
 fHistPtTPCInAccMCtwoSPD(0),
@@ -152,6 +153,20 @@ fHistPtITSMI5InAcc(0),
 fHistPtITSMI6InAcc(0),
 fHistPtITSMISPDInAcc(0),
 fHistPtITSMIoneSPDInAcc(0),
+fHistPhiITSMI2InAcc(0),
+fHistPhiITSMI3InAcc(0),
+fHistPhiITSMI4InAcc(0),
+fHistPhiITSMI5InAcc(0),
+fHistPhiITSMI6InAcc(0),
+fHistPhiITSMISPDInAcc(0),
+fHistPhiITSMIoneSPDInAcc(0),
+fHistEtaITSMI2InAcc(0),
+fHistEtaITSMI3InAcc(0),
+fHistEtaITSMI4InAcc(0),
+fHistEtaITSMI5InAcc(0),
+fHistEtaITSMI6InAcc(0),
+fHistEtaITSMISPDInAcc(0),
+fHistEtaITSMIoneSPDInAcc(0),
 fHistPtITSMI2InAccFake(0),
 fHistPtITSMI3InAccFake(0),
 fHistPtITSMI4InAccFake(0),
@@ -320,6 +335,7 @@ fHistxlocVSmodSDDok(0),
 fHistxlocSDDall(0),
 fHistzlocSDDall(0),
 fHistPhiTPCInAcc(0),
+fHistEtaTPCInAcc(0),
 fHistPtTPC(0),
 fHistPtTPCInAcc(0),
 fHistPtTPCInAccMCtwoSPD(0),
@@ -351,6 +367,20 @@ fHistPtITSMI5InAcc(0),
 fHistPtITSMI6InAcc(0),
 fHistPtITSMISPDInAcc(0),
 fHistPtITSMIoneSPDInAcc(0),
+fHistPhiITSMI2InAcc(0),
+fHistPhiITSMI3InAcc(0),
+fHistPhiITSMI4InAcc(0),
+fHistPhiITSMI5InAcc(0),
+fHistPhiITSMI6InAcc(0),
+fHistPhiITSMISPDInAcc(0),
+fHistPhiITSMIoneSPDInAcc(0),
+fHistEtaITSMI2InAcc(0),
+fHistEtaITSMI3InAcc(0),
+fHistEtaITSMI4InAcc(0),
+fHistEtaITSMI5InAcc(0),
+fHistEtaITSMI6InAcc(0),
+fHistEtaITSMISPDInAcc(0),
+fHistEtaITSMIoneSPDInAcc(0),
 fHistPtITSMI2InAccFake(0),
 fHistPtITSMI3InAccFake(0),
 fHistPtITSMI4InAccFake(0),
@@ -738,6 +768,11 @@ void AliAnalysisTaskITSTrackingCheck::UserCreateOutputObjects()
   fHistPhiTPCInAcc->SetMinimum(0);
   fOutput->Add(fHistPhiTPCInAcc);
 
+  fHistEtaTPCInAcc = new TH1F("fHistEtaTPCInAcc","Eta distribution of TPC tracks; #eta; N tracks",100, -1.5, +1.5);
+  fHistEtaTPCInAcc->Sumw2();
+  fHistEtaTPCInAcc->SetMinimum(0);
+  fOutput->Add(fHistEtaTPCInAcc);
+
   fHistPhiITSMIokbadoutinz6InAcc = new TH1F("fHistPhiITSMIokbadoutinz6InAcc","Azimuthal distribution of ITSMI tracks with 6 layers OK; #phi; N tracks",100,0,2.*3.1415);
   fHistPhiITSMIokbadoutinz6InAcc->Sumw2();
   fHistPhiITSMIokbadoutinz6InAcc->SetMinimum(0);
@@ -867,6 +902,77 @@ void AliAnalysisTaskITSTrackingCheck::UserCreateOutputObjects()
   fHistPtITSMIoneSPDInAcc->SetMinimum(0);
   fOutput->Add(fHistPtITSMIoneSPDInAcc);
 
+  fHistPhiITSMI6InAcc = new TH1F("fHistPhiITSMI6InAcc","phi distribution of ITSMI6 tracks; #phi; N tracks",100,0,2*3.1415);
+  fHistPhiITSMI6InAcc->Sumw2();
+  fHistPhiITSMI6InAcc->SetMinimum(0);
+  fOutput->Add(fHistPhiITSMI6InAcc);
+  
+  fHistPhiITSMI5InAcc = new TH1F("fHistPhiITSMI5InAcc","phi distribution of ITSMI5 tracks; #phi; N tracks",100,0,2*3.1415);
+  fHistPhiITSMI5InAcc->Sumw2();
+  fHistPhiITSMI5InAcc->SetMinimum(0);
+  fOutput->Add(fHistPhiITSMI5InAcc);
+  
+  fHistPhiITSMI4InAcc = new TH1F("fHistPhiITSMI4InAcc","phi distribution of ITSMI4 tracks; #phi; N tracks",100,0,2*3.1415);
+  fHistPhiITSMI4InAcc->Sumw2();
+  fHistPhiITSMI4InAcc->SetMinimum(0);
+  fOutput->Add(fHistPhiITSMI4InAcc);
+  
+  fHistPhiITSMI3InAcc = new TH1F("fHistPhiITSMI3InAcc","phi distribution of ITSMI3 tracks; #phi; N tracks",100,0,2*3.1415);
+  fHistPhiITSMI3InAcc->Sumw2();
+  fHistPhiITSMI3InAcc->SetMinimum(0);
+  fOutput->Add(fHistPhiITSMI3InAcc);
+  
+  fHistPhiITSMI2InAcc = new TH1F("fHistPhiITSMI2InAcc","phi distribution of ITSMI2 tracks; #phi; N tracks",100,0,2*3.1415);
+  fHistPhiITSMI2InAcc->Sumw2();
+  fHistPhiITSMI2InAcc->SetMinimum(0);
+  fOutput->Add(fHistPhiITSMI2InAcc);
+  
+  fHistPhiITSMISPDInAcc = new TH1F("fHistPhiITSMISPDInAcc","phi distribution of ITSMISPD tracks; #phi; N tracks",100,0,2*3.1415);
+  fHistPhiITSMISPDInAcc->Sumw2();
+  fHistPhiITSMISPDInAcc->SetMinimum(0);
+  fOutput->Add(fHistPhiITSMISPDInAcc);
+  
+  fHistPhiITSMIoneSPDInAcc = new TH1F("fHistPhiITSMIoneSPDInAcc","phi distribution of ITSMISPD tracks; #phi; N tracks",100,0,2*3.1415);
+  fHistPhiITSMIoneSPDInAcc->Sumw2();
+  fHistPhiITSMIoneSPDInAcc->SetMinimum(0);
+  fOutput->Add(fHistPhiITSMIoneSPDInAcc);
+
+  fHistEtaITSMI6InAcc = new TH1F("fHistEtaITSMI6InAcc","eta distribution of ITSMI6 tracks; #eta; N tracks",100,-1.5,+1.5);
+  fHistEtaITSMI6InAcc->Sumw2();
+  fHistEtaITSMI6InAcc->SetMinimum(0);
+  fOutput->Add(fHistEtaITSMI6InAcc);
+  
+  fHistEtaITSMI5InAcc = new TH1F("fHistEtaITSMI5InAcc","eta distribution of ITSMI5 tracks; #eta; N tracks",100,-1.5,+1.5);
+  fHistEtaITSMI5InAcc->Sumw2();
+  fHistEtaITSMI5InAcc->SetMinimum(0);
+  fOutput->Add(fHistEtaITSMI5InAcc);
+  
+  fHistEtaITSMI4InAcc = new TH1F("fHistEtaITSMI4InAcc","eta distribution of ITSMI4 tracks; #eta; N tracks",100,-1.5,+1.5);
+  fHistEtaITSMI4InAcc->Sumw2();
+  fHistEtaITSMI4InAcc->SetMinimum(0);
+  fOutput->Add(fHistEtaITSMI4InAcc);
+  
+  fHistEtaITSMI3InAcc = new TH1F("fHistEtaITSMI3InAcc","eta distribution of ITSMI3 tracks; #eta; N tracks",100,-1.5,+1.5);
+  fHistEtaITSMI3InAcc->Sumw2();
+  fHistEtaITSMI3InAcc->SetMinimum(0);
+  fOutput->Add(fHistEtaITSMI3InAcc);
+  
+  fHistEtaITSMI2InAcc = new TH1F("fHistEtaITSMI2InAcc","eta distribution of ITSMI2 tracks; #eta; N tracks",100,-1.5,+1.5);
+  fHistEtaITSMI2InAcc->Sumw2();
+  fHistEtaITSMI2InAcc->SetMinimum(0);
+  fOutput->Add(fHistEtaITSMI2InAcc);
+  
+  fHistEtaITSMISPDInAcc = new TH1F("fHistEtaITSMISPDInAcc","eta distribution of ITSMISPD tracks; #eta; N tracks",100,-1.5,+1.5);
+  fHistEtaITSMISPDInAcc->Sumw2();
+  fHistEtaITSMISPDInAcc->SetMinimum(0);
+  fOutput->Add(fHistEtaITSMISPDInAcc);
+  
+  fHistEtaITSMIoneSPDInAcc = new TH1F("fHistEtaITSMIoneSPDInAcc","eta distribution of ITSMISPD tracks; #eta; N tracks",100,-1.5,+1.5);
+  fHistEtaITSMIoneSPDInAcc->Sumw2();
+  fHistEtaITSMIoneSPDInAcc->SetMinimum(0);
+  fOutput->Add(fHistEtaITSMIoneSPDInAcc);
+
+
   fHistPtITSMI6InAccFake = new TH1F("fHistPtITSMI6InAccFake","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
   fHistPtITSMI6InAccFake->Sumw2();
   fHistPtITSMI6InAccFake->SetMinimum(0);
@@ -1781,7 +1887,10 @@ void AliAnalysisTaskITSTrackingCheck::UserExec(Option_t *)
          Double_t ptTPCinnerwall=TMath::Sqrt(pTPCinnerwall[0]*pTPCinnerwall[0]+pTPCinnerwall[1]*pTPCinnerwall[1]);
          fHistPtTPCInAccWithPtTPCAtInnerWall->Fill(ptTPCinnerwall);
          if(track->GetStatus()&AliESDtrack::kITSrefit && !(track->GetStatus()&AliESDtrack::kTOFout)) fHistDeltaPtTPC->Fill(trackTPC->P(),track->P()-trackTPC->P());
-         fHistPhiTPCInAcc->Fill(track->Phi());  
+         if(track->Pt()>1) {
+           fHistPhiTPCInAcc->Fill(track->Phi());  
+           fHistEtaTPCInAcc->Fill(track->Eta());
+         }  
          if(isPrimary) {
            fHistPtTPCInAccP->Fill(track->Pt());
            if(pdgTrk==321) fHistPtTPCInAccPfromStrange->Fill(track->Pt());
@@ -1827,36 +1936,64 @@ void AliAnalysisTaskITSTrackingCheck::UserExec(Option_t *)
       if(itsfindableAcc) {
        if(nclsITS==6) {
          fHistPtITSMI6InAcc->Fill(track->Pt());
+         if(track->Pt()>1) {
+           fHistPhiITSMI6InAcc->Fill(track->Phi());
+           fHistEtaITSMI6InAcc->Fill(track->Eta());
+         }
          if(isFake) fHistPtITSMI6InAccFake->Fill(track->Pt());
          if(isPrimary) {fHistPtITSMI6InAccP->Fill(track->Pt());} else {fHistPtITSMI6InAccS->Fill(track->Pt());}  
        }
        if(nclsITS==5) {
          fHistPtITSMI5InAcc->Fill(track->Pt());
+         if(track->Pt()>1) {
+           fHistPhiITSMI5InAcc->Fill(track->Phi());
+           fHistEtaITSMI5InAcc->Fill(track->Eta());
+         }
          if(isFake) fHistPtITSMI5InAccFake->Fill(track->Pt());
          if(isPrimary) {fHistPtITSMI5InAccP->Fill(track->Pt());} else {fHistPtITSMI5InAccS->Fill(track->Pt());}  
        }
        if(nclsITS==4) {
          fHistPtITSMI4InAcc->Fill(track->Pt());
+         if(track->Pt()>1) {
+           fHistPhiITSMI4InAcc->Fill(track->Phi());
+           fHistEtaITSMI4InAcc->Fill(track->Eta());
+         }
          if(isFake) fHistPtITSMI4InAccFake->Fill(track->Pt());
          if(isPrimary) {fHistPtITSMI4InAccP->Fill(track->Pt());} else {fHistPtITSMI4InAccS->Fill(track->Pt());}  
        }
        if(nclsITS==3) {
          fHistPtITSMI3InAcc->Fill(track->Pt());
+         if(track->Pt()>1) {
+           fHistPhiITSMI3InAcc->Fill(track->Phi());
+           fHistEtaITSMI3InAcc->Fill(track->Eta());
+         }
          if(isFake) fHistPtITSMI3InAccFake->Fill(track->Pt());
          if(isPrimary) {fHistPtITSMI3InAccP->Fill(track->Pt());} else {fHistPtITSMI3InAccS->Fill(track->Pt());}  
        }
        if(nclsITS==2) {
          fHistPtITSMI2InAcc->Fill(track->Pt());
+         if(track->Pt()>1) {
+           fHistPhiITSMI2InAcc->Fill(track->Phi());
+           fHistEtaITSMI2InAcc->Fill(track->Eta());
+         }
          if(isFake) fHistPtITSMI2InAccFake->Fill(track->Pt());
          if(isPrimary) {fHistPtITSMI2InAccP->Fill(track->Pt());} else {fHistPtITSMI2InAccS->Fill(track->Pt());}  
        }
        if(track->HasPointOnITSLayer(0) && track->HasPointOnITSLayer(1)) {
          fHistPtITSMISPDInAcc->Fill(track->Pt());
+         if(track->Pt()>1) {
+           fHistPhiITSMISPDInAcc->Fill(track->Phi());
+           fHistEtaITSMISPDInAcc->Fill(track->Eta());
+         }
          if(isFake) fHistPtITSMISPDInAccFake->Fill(track->Pt());
          if(isPrimary) {fHistPtITSMISPDInAccP->Fill(track->Pt());} else {fHistPtITSMISPDInAccS->Fill(track->Pt());}  
        }
        if(track->HasPointOnITSLayer(0) || track->HasPointOnITSLayer(1)) {
          fHistPtITSMIoneSPDInAcc->Fill(track->Pt());
+         if(track->Pt()>1) {
+           fHistPhiITSMIoneSPDInAcc->Fill(track->Phi());
+           fHistEtaITSMIoneSPDInAcc->Fill(track->Eta());
+         }
          Float_t chi2redITS = track->GetITSchi2()/track->GetNcls(0);
          if(track->Pt()>0.18 && track->Pt()<0.25)  fHistITSRedChi2NonFakePt02->Fill(chi2redITS);
          if(track->Pt()>0.4 && track->Pt()<0.6)  fHistITSRedChi2NonFakePt05->Fill(chi2redITS);
@@ -1874,7 +2011,7 @@ void AliAnalysisTaskITSTrackingCheck::UserExec(Option_t *)
        if(nclsokbadoutinzITS==6) fHistPtITSMIokbadoutinz6InAcc->Fill(track->Pt());
        if(nclsokbadoutinzITS==5) fHistPtITSMIokbadoutinz5InAcc->Fill(track->Pt());
        if(nclsokbadoutinzITS==4) fHistPtITSMIokbadoutinz4InAcc->Fill(track->Pt());
-       if(nclsokbadoutinzITS==6) fHistPhiITSMIokbadoutinz6InAcc->Fill(track->Phi());  
+       if(nclsokbadoutinzITS==6 && track->Pt()>1) fHistPhiITSMIokbadoutinz6InAcc->Fill(track->Phi());  
       }
     }
 
index 33e02c7..d5d1483 100644 (file)
@@ -102,6 +102,7 @@ class AliAnalysisTaskITSTrackingCheck : public AliAnalysisTaskSE
   TH1F         *fHistxlocSDDall; //! output hist
   TH1F         *fHistzlocSDDall; //! output hist
   TH1F         *fHistPhiTPCInAcc; //! output hist
+  TH1F         *fHistEtaTPCInAcc; //! output hist
   TH1F         *fHistPtTPC; //! output hist
   TH1F         *fHistPtTPCInAcc; //! output hist
   TH1F         *fHistPtTPCInAccMCtwoSPD; //! output hist
@@ -133,6 +134,20 @@ class AliAnalysisTaskITSTrackingCheck : public AliAnalysisTaskSE
   TH1F         *fHistPtITSMI6InAcc; //! output hist
   TH1F         *fHistPtITSMISPDInAcc; //! output hist
   TH1F         *fHistPtITSMIoneSPDInAcc; //! output hist
+  TH1F         *fHistPhiITSMI2InAcc; //! output hist
+  TH1F         *fHistPhiITSMI3InAcc; //! output hist
+  TH1F         *fHistPhiITSMI4InAcc; //! output hist
+  TH1F         *fHistPhiITSMI5InAcc; //! output hist
+  TH1F         *fHistPhiITSMI6InAcc; //! output hist
+  TH1F         *fHistPhiITSMISPDInAcc; //! output hist
+  TH1F         *fHistPhiITSMIoneSPDInAcc; //! output hist
+  TH1F         *fHistEtaITSMI2InAcc; //! output hist
+  TH1F         *fHistEtaITSMI3InAcc; //! output hist
+  TH1F         *fHistEtaITSMI4InAcc; //! output hist
+  TH1F         *fHistEtaITSMI5InAcc; //! output hist
+  TH1F         *fHistEtaITSMI6InAcc; //! output hist
+  TH1F         *fHistEtaITSMISPDInAcc; //! output hist
+  TH1F         *fHistEtaITSMIoneSPDInAcc; //! output hist
   TH1F         *fHistPtITSMI2InAccFake; //! output hist
   TH1F         *fHistPtITSMI3InAccFake; //! output hist
   TH1F         *fHistPtITSMI4InAccFake; //! output hist
@@ -250,7 +265,7 @@ class AliAnalysisTaskITSTrackingCheck : public AliAnalysisTaskSE
   Int_t MakeITSflag(AliESDtrack *track) const;
   Bool_t IsSelectedCentrality() const;
 
-  ClassDef(AliAnalysisTaskITSTrackingCheck,11); // ITS tracks analysis
+  ClassDef(AliAnalysisTaskITSTrackingCheck,12); // ITS tracks analysis
 };
 
 #endif