X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ANALYSIS%2FAliAnalysisTaskPIDqa.cxx;h=5e1d623f43aa1796d05ade72fc762c2e5e2490cb;hb=9120829d30a51779aa89933dd52a0700b5b0f7aa;hp=e1eff32802d3071d504bf8dbdcd475cf2d42f451;hpb=f25dc3ef023b507057ca7f5878751cb1c61e47b1;p=u%2Fmrichter%2FAliRoot.git

diff --git a/ANALYSIS/AliAnalysisTaskPIDqa.cxx b/ANALYSIS/AliAnalysisTaskPIDqa.cxx
index e1eff32802d..5e1d623f43a 100644
--- a/ANALYSIS/AliAnalysisTaskPIDqa.cxx
+++ b/ANALYSIS/AliAnalysisTaskPIDqa.cxx
@@ -381,13 +381,13 @@ void AliAnalysisTaskPIDqa::FillITSqa()
     }
     
     
-    for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+    for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
       TH2 *h=(TH2*)theList->At(ispecie);
       if (!h) continue;
       Double_t nSigma=fPIDResponse->NumberOfSigmasITS(track, (AliPID::EParticleType)ispecie);
       h->Fill(mom,nSigma);
     }
-    TH2 *h=(TH2*)theList->At(AliPID::kSPECIES);
+    TH2 *h=(TH2*)theList->At(AliPID::kSPECIESC);
     if (h) {
       Double_t sig=track->GetITSsignal();
       h->Fill(mom,sig);
@@ -428,15 +428,35 @@ void AliAnalysisTaskPIDqa::FillTPCqa()
     if ( nCrossedRowsTPC<70 || ratioCrossedRowsOverFindableClustersTPC<.8 ) continue;
     
     Double_t mom=track->GetTPCmomentum();
-    
-    for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+    // the default scenario
+    for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
       TH2 *h=(TH2*)fListQAtpc->At(ispecie);
       if (!h) continue;
       Double_t nSigma=fPIDResponse->NumberOfSigmasTPC(track, (AliPID::EParticleType)ispecie);
       h->Fill(mom,nSigma);
     }
+    // the "hybrid" scenario
+    for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
+      TH2 *h=(TH2*)fListQAtpc->At(ispecie+AliPID::kSPECIESC);
+      if (!h) continue;
+      Double_t nSigma=fPIDResponse->NumberOfSigmasTPC(track, (AliPID::EParticleType)ispecie, AliTPCPIDResponse::kdEdxHybrid);
+      h->Fill(mom,nSigma);
+    }
+    
+    // the "OROC" scenario
+    for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
+      TH2 *h=(TH2*)fListQAtpc->At(ispecie+2*AliPID::kSPECIESC);
+      if (!h) continue;
+      Double_t nSigma=fPIDResponse->NumberOfSigmasTPC(track, (AliPID::EParticleType)ispecie, AliTPCPIDResponse::kdEdxOROC);
+      //TSpline3* spline = fPIDResponse->GetTPCResponse().GetCurrentResponseFunction();
+      //std::cout<<ispecie<<" "<<nSigma<<" phi:"<<track->Phi()<<". "<<std::endl;
+      //if (spline) {cout<<spline->GetName()<<endl;}
+      //else {cout<<"NULL spline"<<endl;}
+      h->Fill(mom,nSigma);
+    }
     
-    TH2 *h=(TH2*)fListQAtpc->At(AliPID::kSPECIES);
+    TH2 *h=(TH2*)fListQAtpc->At(3*AliPID::kSPECIESC);
+
     if (h) {
       Double_t sig=track->GetTPCsignal();
       h->Fill(mom,sig);
@@ -529,7 +549,7 @@ void AliAnalysisTaskPIDqa::FillTOFqa()
 
     Double_t mom=track->P();
 
-    for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+    for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
       TH2 *h=(TH2*)fListQAtof->At(ispecie);
       if (!h) continue;
       Double_t nSigma=fPIDResponse->NumberOfSigmasTOF(track, (AliPID::EParticleType)ispecie);
@@ -875,7 +895,7 @@ void AliAnalysisTaskPIDqa::FillTPCTOFqa()
     Double_t mom=track->P();
     Double_t momTPC=track->GetTPCmomentum();
 
-    for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+    for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
       //TOF nSigma
       Double_t nSigmaTOF=fPIDResponse->NumberOfSigmasTOF(track, (AliPID::EParticleType)ispecie);
       Double_t nSigmaTPC=fPIDResponse->NumberOfSigmasTPC(track, (AliPID::EParticleType)ispecie);
@@ -885,11 +905,11 @@ void AliAnalysisTaskPIDqa::FillTPCTOFqa()
       if (h && TMath::Abs(nSigmaTOF)<3.) h->Fill(momTPC,nSigmaTPC);
 
       //TOF after TPC cut
-      h=(TH2*)fListQAtpctof->At(ispecie+AliPID::kSPECIES);
+      h=(TH2*)fListQAtpctof->At(ispecie+AliPID::kSPECIESC);
       if (h && TMath::Abs(nSigmaTPC)<3.) h->Fill(mom,nSigmaTOF);
 
       //EMCAL after TOF and TPC cut
-      h=(TH2*)fListQAtpctof->At(ispecie+2*AliPID::kSPECIES);
+      h=(TH2*)fListQAtpctof->At(ispecie+2*AliPID::kSPECIESC);
       if (h && TMath::Abs(nSigmaTOF)<3. && TMath::Abs(nSigmaTPC)<3. ){
 
 	Int_t nMatchClus = track->GetEMCALcluster();
@@ -929,7 +949,7 @@ void AliAnalysisTaskPIDqa::SetupITSqa()
   TVectorD *vX=MakeLogBinning(200,.1,30);
   
   //ITS+TPC tracks
-  for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
     TH2F *hNsigmaP = new TH2F(Form("hNsigmaP_ITS_%s",AliPID::ParticleName(ispecie)),
                               Form("ITS n#sigma %s vs. p;p [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
                               vX->GetNrows()-1,vX->GetMatrixArray(),
@@ -943,7 +963,7 @@ void AliAnalysisTaskPIDqa::SetupITSqa()
   fListQAits->Add(hSig);
 
   //ITS Standalone tracks
-  for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
     TH2F *hNsigmaPSA = new TH2F(Form("hNsigmaP_ITSSA_%s",AliPID::ParticleName(ispecie)),
 				Form("ITS n#sigma %s vs. p;p [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
 				vX->GetNrows()-1,vX->GetMatrixArray(),
@@ -957,7 +977,7 @@ void AliAnalysisTaskPIDqa::SetupITSqa()
   fListQAitsSA->Add(hSigSA);
   
   //ITS Pure Standalone tracks
-  for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
     TH2F *hNsigmaPPureSA = new TH2F(Form("hNsigmaP_ITSPureSA_%s",AliPID::ParticleName(ispecie)),
 				    Form("ITS n#sigma %s vs. p;p [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
 				    vX->GetNrows()-1,vX->GetMatrixArray(),
@@ -982,20 +1002,39 @@ void AliAnalysisTaskPIDqa::SetupTPCqa()
   
   TVectorD *vX=MakeLogBinning(200,.1,30);
   
-  for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
     TH2F *hNsigmaP = new TH2F(Form("hNsigmaP_TPC_%s",AliPID::ParticleName(ispecie)),
                               Form("TPC n#sigma %s vs. p;p [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
                               vX->GetNrows()-1,vX->GetMatrixArray(),
                               200,-10,10);
     fListQAtpc->Add(hNsigmaP);
   }
+
+  // the "hybrid" scenario
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
+    TH2F *hNsigmaP = new TH2F(Form("hNsigmaP_TPC_%s_Hybrid",AliPID::ParticleName(ispecie)),
+                              Form("TPC n#sigma %s vs. p (Hybrid gain scenario);p [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
+                              vX->GetNrows()-1,vX->GetMatrixArray(),
+                              200,-10,10);
+    fListQAtpc->Add(hNsigmaP);
+  }
+   
+  // the "OROC high" scenario
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
+    TH2F *hNsigmaP = new TH2F(Form("hNsigmaP_TPC_%s_OROChigh",AliPID::ParticleName(ispecie)),
+                              Form("TPC n#sigma %s vs. p (OROChigh gain scenario);p [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
+                              vX->GetNrows()-1,vX->GetMatrixArray(),
+                              200,-10,10);
+    fListQAtpc->Add(hNsigmaP);
+  }
+  
   
   
   TH2F *hSig = new TH2F("hSigP_TPC",
                         "TPC signal vs. p;p [GeV]; TPC signal [arb. units]",
                         vX->GetNrows()-1,vX->GetMatrixArray(),
                         300,0,300);
-  fListQAtpc->Add(hSig);
+  fListQAtpc->Add(hSig); //3*AliPID::kSPECIESC
 
   delete vX;  
 }
@@ -1028,7 +1067,7 @@ void AliAnalysisTaskPIDqa::SetupTOFqa()
   
   TVectorD *vX=MakeLogBinning(200,.1,30);
 
-  for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
     TH2F *hNsigmaP = new TH2F(Form("hNsigmaP_TOF_%s",AliPID::ParticleName(ispecie)),
                               Form("TOF n#sigma %s vs. p;p [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
                               vX->GetNrows()-1,vX->GetMatrixArray(),
@@ -1158,7 +1197,7 @@ void AliAnalysisTaskPIDqa::SetupTPCTOFqa()
   TVectorD *vX=MakeLogBinning(200,.1,30);
 
   //TPC signals after TOF cut
-  for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
     TH2F *hNsigmaP = new TH2F(Form("hNsigmaP_TPC_TOF_%s",AliPID::ParticleName(ispecie)),
                               Form("TPC n#sigma %s vs. p (after TOF 3#sigma cut);p_{TPC} [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
                               vX->GetNrows()-1,vX->GetMatrixArray(),
@@ -1167,7 +1206,7 @@ void AliAnalysisTaskPIDqa::SetupTPCTOFqa()
   }
 
   //TOF signals after TPC cut
-  for (Int_t ispecie=0; ispecie<AliPID::kSPECIES; ++ispecie){
+  for (Int_t ispecie=0; ispecie<AliPID::kSPECIESC; ++ispecie){
     TH2F *hNsigmaP = new TH2F(Form("hNsigmaP_TOF_TPC_%s",AliPID::ParticleName(ispecie)),
                               Form("TOF n#sigma %s vs. p (after TPC n#sigma cut);p [GeV]; n#sigma",AliPID::ParticleName(ispecie)),
                               vX->GetNrows()-1,vX->GetMatrixArray(),