New histogramas and new cuts for QA
authorpgonzale <pedro.gonzalez.zamora@cern.ch>
Mon, 24 Nov 2014 17:01:14 +0000 (18:01 +0100)
committerpgonzale <pedro.gonzalez.zamora@cern.ch>
Mon, 24 Nov 2014 17:02:09 +0000 (18:02 +0100)
PWGGA/GammaConv/AliAnalysisTaskGammaConvDalitzV1.cxx
PWGGA/GammaConv/AliAnalysisTaskGammaConvDalitzV1.h
PWGGA/GammaConv/macros/AddTask_GammaConvDalitzQAV1_pPb.C

index 02edbf1..88f2c72 100644 (file)
@@ -1,7 +1,7 @@
 /**************************************************************************
  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  *                                                                        *
- * Author: Pedro Gonz??lez, Pedro Ladr??n de Guevara, Ernesto L??pez Torres, *
+ * Author: Pedro Gonzalez, Pedro Ladron de Guevara, Ernesto Lopez Torres, *
  *         Eulogio Serradilla, Ana Marin, Friederike Bock                 *
  * Version 2                                                              *
  *                                                                        *
@@ -103,13 +103,21 @@ AliAnalysisTaskGammaConvDalitzV1::AliAnalysisTaskGammaConvDalitzV1():
        hESDDalitzElectronAfterPhi(NULL),
        hESDDalitzPositronAfterPhi(NULL),
        hESDDalitzElectronAfterNClsITS(NULL),
+       hESDDalitzElectronAfterNClsITSPCut(NULL),
        hESDDalitzPositronAfterNClsITS(NULL),
+       hESDDalitzPositronAfterNClsITSPCut(NULL),
        hESDDalitzElectronAfterNFindClsTPC(NULL),
+       hESDDalitzElectronAfterNFindClsTPCPCut(NULL),
        hESDDalitzPositronAfterNFindClsTPC(NULL),
+       hESDDalitzPositronAfterNFindClsTPCPCut(NULL),
        hESDDalitzElectronAfterNClsTPC(NULL),
+       hESDDalitzElectronAfterNClsTPCPCut(NULL),
        hESDDalitzPositronAfterNClsTPC(NULL),
+       hESDDalitzPositronAfterNClsTPCPCut(NULL),
        hESDDalitzElectronAfterNCrossedRowsTPC(NULL),
+       hESDDalitzElectronAfterNCrossedRowsTPCPCut(NULL),
        hESDDalitzPositronAfterNCrossedRowsTPC(NULL),
+       hESDDalitzPositronAfterNCrossedRowsTPCPCut(NULL),
        hESDDalitzPosEleAfterDCAxy(NULL),
        hESDDalitzPosEleAfterDCAz(NULL),
        hESDDalitzElectronAfterTPCdEdxVsP(NULL),
@@ -275,13 +283,21 @@ AliAnalysisTaskGammaConvDalitzV1::AliAnalysisTaskGammaConvDalitzV1( const char*
        hESDDalitzElectronAfterPhi(NULL),
        hESDDalitzPositronAfterPhi(NULL),
        hESDDalitzElectronAfterNClsITS(NULL),
+       hESDDalitzElectronAfterNClsITSPCut(NULL),
        hESDDalitzPositronAfterNClsITS(NULL),
+       hESDDalitzPositronAfterNClsITSPCut(NULL),
        hESDDalitzElectronAfterNFindClsTPC(NULL),
+       hESDDalitzElectronAfterNFindClsTPCPCut(NULL),
        hESDDalitzPositronAfterNFindClsTPC(NULL),
+       hESDDalitzPositronAfterNFindClsTPCPCut(NULL),
        hESDDalitzElectronAfterNClsTPC(NULL),
+       hESDDalitzElectronAfterNClsTPCPCut(NULL),
        hESDDalitzPositronAfterNClsTPC(NULL),
+       hESDDalitzPositronAfterNClsTPCPCut(NULL),
        hESDDalitzElectronAfterNCrossedRowsTPC(NULL),
+       hESDDalitzElectronAfterNCrossedRowsTPCPCut(NULL),
        hESDDalitzPositronAfterNCrossedRowsTPC(NULL),
+       hESDDalitzPositronAfterNCrossedRowsTPCPCut(NULL),
        hESDDalitzPosEleAfterDCAxy(NULL),
        hESDDalitzPosEleAfterDCAz(NULL),
        hESDDalitzElectronAfterTPCdEdxVsP(NULL),
@@ -558,13 +574,21 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
                hESDDalitzElectronAfterPhi                              = new TH1F*[fnCuts];
                hESDDalitzPositronAfterPhi                              = new TH1F*[fnCuts];
                hESDDalitzElectronAfterNClsITS                  = new TH1F*[fnCuts];
+               hESDDalitzElectronAfterNClsITSPCut              = new TH1F*[fnCuts];
                hESDDalitzPositronAfterNClsITS                  = new TH1F*[fnCuts];
+               hESDDalitzPositronAfterNClsITSPCut              = new TH1F*[fnCuts];
                hESDDalitzElectronAfterNFindClsTPC              = new TH2F*[fnCuts];
+               hESDDalitzElectronAfterNFindClsTPCPCut          = new TH1F*[fnCuts];
                hESDDalitzPositronAfterNFindClsTPC              = new TH2F*[fnCuts];
+               hESDDalitzPositronAfterNFindClsTPCPCut          = new TH1F*[fnCuts];
                hESDDalitzElectronAfterNClsTPC                  = new TH2F*[fnCuts];
+               hESDDalitzElectronAfterNClsTPCPCut              = new TH1F*[fnCuts];
                hESDDalitzPositronAfterNClsTPC                  = new TH2F*[fnCuts];
-               hESDDalitzElectronAfterNCrossedRowsTPC  = new TH2F*[fnCuts];
-               hESDDalitzPositronAfterNCrossedRowsTPC  = new TH2F*[fnCuts];
+               hESDDalitzPositronAfterNClsTPCPCut              = new TH1F*[fnCuts];
+               hESDDalitzElectronAfterNCrossedRowsTPC          = new TH2F*[fnCuts];
+               hESDDalitzElectronAfterNCrossedRowsTPCPCut      = new TH1F*[fnCuts];
+               hESDDalitzPositronAfterNCrossedRowsTPC          = new TH2F*[fnCuts];
+               hESDDalitzPositronAfterNCrossedRowsTPCPCut      = new TH1F*[fnCuts];
                hESDDalitzPosEleAfterDCAxy                              = new TH2F*[fnCuts];
                hESDDalitzPosEleAfterDCAz                               = new TH2F*[fnCuts];
                hESDDalitzElectronAfterTPCdEdxVsP               = new TH2F*[fnCuts];
@@ -720,26 +744,52 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
                        hESDDalitzElectronAfterNClsITS[iCut]  = new TH1F("ESD_DalitzElectron_After_NClsITS","ESD_DalitzElectron_After_NClsITS",7,0.,7.);
                        fQAFolder[iCut]->Add(hESDDalitzElectronAfterNClsITS[iCut]);
                        
+                       hESDDalitzElectronAfterNClsITSPCut[iCut]  = new TH1F("ESD_DalitzElectron_After_NClsITS_PCut","ESD_DalitzElectron_After_NClsITS_PCut",7,0.,7.);
+                       fQAFolder[iCut]->Add(hESDDalitzElectronAfterNClsITSPCut[iCut]);
+                       
                        hESDDalitzPositronAfterNClsITS[iCut]  = new TH1F("ESD_DalitzPositron_After_NClsITS","ESD_DalitzPositron_After_NClsITS",7,0.,7.);
                        fQAFolder[iCut]->Add(hESDDalitzPositronAfterNClsITS[iCut]);
                        
+                       hESDDalitzPositronAfterNClsITSPCut[iCut]  = new TH1F("ESD_DalitzPositron_After_NClsITS_PCut","ESD_DalitzPositron_After_NClsITS_PCut",7,0.,7.);
+                       fQAFolder[iCut]->Add(hESDDalitzPositronAfterNClsITSPCut[iCut]);
+                                               
                        hESDDalitzElectronAfterNFindClsTPC[iCut]  = new TH2F("ESD_DalitzElectron_After_NFindClsTPC","ESD_DalitzElectron_After_NFindClsTPC",60,0,1.5,kPtBins,binsPtDummy);
                        fQAFolder[iCut]->Add(hESDDalitzElectronAfterNFindClsTPC[iCut]);
                        
+                       hESDDalitzElectronAfterNFindClsTPCPCut[iCut]  = new TH1F("ESD_DalitzElectron_After_NFindClsTPC_PCut","ESD_DalitzElectron_After_NFindClsTPC_PCut",60,0,1.5);
+                       fQAFolder[iCut]->Add(hESDDalitzElectronAfterNFindClsTPCPCut[iCut]);
+                       
+                       
                        hESDDalitzPositronAfterNFindClsTPC[iCut]  = new TH2F("ESD_DalitzPositron_After_NFindClsTPC","ESD_DalitzPositron_After_NFindClsTPC",60,0,1.5,kPtBins,binsPtDummy);
                        fQAFolder[iCut]->Add(hESDDalitzPositronAfterNFindClsTPC[iCut]);
                        
+                       hESDDalitzPositronAfterNFindClsTPCPCut[iCut]  = new TH1F("ESD_DalitzPositron_After_NFindClsTPC_PCut","ESD_DalitzPositron_After_NFindClsTPC_PCut",60,0,1.5);
+                       fQAFolder[iCut]->Add(hESDDalitzPositronAfterNFindClsTPCPCut[iCut]);
+                       
+                       
                        hESDDalitzElectronAfterNClsTPC[iCut]  = new TH2F("ESD_DalitzElectron_After_NClsTPC","ESD_DalitzElectron_After_NClsTPC",200,0,200,kPtBins,binsPtDummy);
                        fQAFolder[iCut]->Add(hESDDalitzElectronAfterNClsTPC[iCut]);
                        
+                       hESDDalitzElectronAfterNClsTPCPCut[iCut]  = new TH1F("ESD_DalitzElectron_After_NClsTPC_PCut","ESD_DalitzElectron_After_NClsTPC_PCut",200,0,200);
+                       fQAFolder[iCut]->Add(hESDDalitzElectronAfterNClsTPCPCut[iCut]);
+                       
                        hESDDalitzPositronAfterNClsTPC[iCut]  = new TH2F("ESD_DalitzPositron_After_NClsTPC","ESD_DalitzPositron_After_NClsTPC",200,0,200,kPtBins,binsPtDummy);
                        fQAFolder[iCut]->Add(hESDDalitzPositronAfterNClsTPC[iCut]);
+                       
+                       hESDDalitzPositronAfterNClsTPCPCut[iCut]  = new TH1F("ESD_DalitzPositron_After_NClsTPC_PCut","ESD_DalitzPositron_After_NClsTPC_PCut",200,0,200);
+                       fQAFolder[iCut]->Add(hESDDalitzPositronAfterNClsTPCPCut[iCut]);
                                        
                        hESDDalitzPositronAfterNCrossedRowsTPC[iCut] =  new TH2F("ESD_DalitzPositron_After_NCrossedRowsTPC","ESD_DalitzPositron_After_NCrossedRowsTPC",165,-0.5,164.5,kPtBins,binsPtDummy);
                        fQAFolder[iCut]->Add(hESDDalitzPositronAfterNCrossedRowsTPC[iCut]);
                        
+                       hESDDalitzPositronAfterNCrossedRowsTPCPCut[iCut] =  new TH1F("ESD_DalitzPositron_After_NCrossedRowsTPC_PCut","ESD_DalitzPositron_After_NCrossedRowsTPC_PCut",165,-0.5,164.5);
+                       fQAFolder[iCut]->Add(hESDDalitzPositronAfterNCrossedRowsTPCPCut[iCut]);
+                                               
                        hESDDalitzElectronAfterNCrossedRowsTPC[iCut] =  new TH2F("ESD_DalitzElectron_After_NCrossedRowsTPC","ESD_DalitzElectron_After_NCrossedRowsTPC",165,-0.5,164.5,kPtBins,binsPtDummy);
                        fQAFolder[iCut]->Add(hESDDalitzElectronAfterNCrossedRowsTPC[iCut]);
+                       
+                       hESDDalitzElectronAfterNCrossedRowsTPCPCut[iCut] =  new TH1F("ESD_DalitzElectron_After_NCrossedRowsTPC_PCut","ESD_DalitzElectron_After_NCrossedRowsTPC_PCut",165,-0.5,164.5);
+                       fQAFolder[iCut]->Add(hESDDalitzElectronAfterNCrossedRowsTPCPCut[iCut]);
                                        
                        hESDDalitzPosEleAfterDCAxy[iCut] = new TH2F("ESD_DalitzPosEle_After_DCAxy","ESD_DalitzPosEle_After_DCAxy",124,-0.62,0.62,kPtBins,binsPtDummy);
                        fQAFolder[iCut]->Add(hESDDalitzPosEleAfterDCAxy[iCut]);
@@ -2685,10 +2735,20 @@ void AliAnalysisTaskGammaConvDalitzV1::FillElectronQAHistos(AliAODConversionPhot
        dcaToVertexZNeg  = bNeg[1];
                
        if( electronVgamma->P() > 0.3 &&  electronVgamma->P() < 0.45 ){ 
+         
                hESDDalitzElectronAfterEtaPCut[fiCut]->Fill( electronVgamma->Eta() );
+               hESDDalitzElectronAfterNClsITSPCut[fiCut]->Fill( NumClsITSNeg );
+               hESDDalitzElectronAfterNFindClsTPCPCut[fiCut]->Fill(clsToFNeg);
+               hESDDalitzElectronAfterNClsTPCPCut[fiCut]->Fill( NumClsTPCNeg );
+               hESDDalitzElectronAfterNCrossedRowsTPCPCut[fiCut]->Fill( nCrossedRowsTPCNeg );
        }
        if( positronVgamma->P() > 0.3 &&  positronVgamma->P() < 0.45 ){ 
+         
                hESDDalitzPositronAfterEtaPCut[fiCut]->Fill( positronVgamma->Eta() );
+               hESDDalitzPositronAfterNClsITSPCut[fiCut]->Fill( NumClsITSPos );
+               hESDDalitzPositronAfterNFindClsTPCPCut[fiCut]->Fill(clsToFPos);
+               hESDDalitzPositronAfterNClsTPCPCut[fiCut]->Fill( NumClsTPCPos );
+               hESDDalitzPositronAfterNCrossedRowsTPCPCut[fiCut]->Fill( nCrossedRowsTPCPos );
        }
        
        hESDDalitzElectronAfterPt[fiCut]->Fill(  electronVgamma->Pt()  );
index 2f2a9cd..31d4b43 100644 (file)
@@ -126,13 +126,21 @@ class AliAnalysisTaskGammaConvDalitzV1: public AliAnalysisTaskSE
                TH1F                                                                    **hESDDalitzElectronAfterPhi;
                TH1F                                                                    **hESDDalitzPositronAfterPhi;
                TH1F                                                                    **hESDDalitzElectronAfterNClsITS;
+               TH1F                                                                    **hESDDalitzElectronAfterNClsITSPCut;
                TH1F                                                                    **hESDDalitzPositronAfterNClsITS;
+               TH1F                                                                    **hESDDalitzPositronAfterNClsITSPCut;
                TH2F                                                                    **hESDDalitzElectronAfterNFindClsTPC;
+               TH1F                                                                    **hESDDalitzElectronAfterNFindClsTPCPCut;
                TH2F                                                                    **hESDDalitzPositronAfterNFindClsTPC;
+               TH1F                                                                    **hESDDalitzPositronAfterNFindClsTPCPCut;
                TH2F                                                                    **hESDDalitzElectronAfterNClsTPC;
+               TH1F                                                                    **hESDDalitzElectronAfterNClsTPCPCut;
                TH2F                                                                    **hESDDalitzPositronAfterNClsTPC;
+               TH1F                                                                    **hESDDalitzPositronAfterNClsTPCPCut;
                TH2F                                                                    **hESDDalitzElectronAfterNCrossedRowsTPC;
+               TH1F                                                                    **hESDDalitzElectronAfterNCrossedRowsTPCPCut;
                TH2F                                                                    **hESDDalitzPositronAfterNCrossedRowsTPC;
+               TH1F                                                                    **hESDDalitzPositronAfterNCrossedRowsTPCPCut;
                TH2F                                                                    **hESDDalitzPosEleAfterDCAxy;
                TH2F                                                                    **hESDDalitzPosEleAfterDCAz;
                TH2F                                                                    **hESDDalitzElectronAfterTPCdEdxVsP;
index 3dc452f..d0b762f 100644 (file)
@@ -72,7 +72,7 @@ void AddTask_GammaConvDalitzQAV1_pPb(  Int_t trainConfig = 1,
    
    if(trainConfig == 9 || trainConfig == 10 ){
    cutnumberPhoton = "160084001001500000000";   //Offline  V0 finder 
-   } else if ( trainConfig == 54  ){   // Warning!!!   Those trains must run in stand-alone mode
+   } else if ( trainConfig >= 54 && trainConfig <= 60  ){   // Warning!!!   Those trains must run in stand-alone mode
    cutnumberPhoton = "000700040000500000000";  
    } else {
    cutnumberPhoton = "060084001000500000000";   //Online  V0 finder //change         
@@ -434,6 +434,30 @@ if( trainConfig == 1 ) {  // No eta shift |Y| < 0.8
   
         eventCutArray[0]="8000011"; photonCutArray[0] = "000000040000500000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny    + new psiPair Cut    0.60, 0.0 0.12   + 0   cm < Rconv  < 180 cm   //open cuts
        
+} else if( trainConfig == 55 ) {
+      
+       eventCutArray[0]="8000011"; photonCutArray[0] = "002000040000500000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Rconv > 5 cm
+       
+} else if( trainConfig == 56 ) {
+      
+       eventCutArray[0]="8000011"; photonCutArray[0] = "009000040000500000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Rconv > 7.5 cm
+       
+} else if( trainConfig == 57 ) {
+      
+       eventCutArray[0]="8000011"; photonCutArray[0] = "005000040000500000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Rconv > 10 cm
+       
+} else if( trainConfig == 58 ) {
+      
+       eventCutArray[0]="8000011"; photonCutArray[0] = "008000040000500000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Rconv > 12.5 cm
+       
+} else if( trainConfig == 59 ) {
+      
+       eventCutArray[0]="8000011"; photonCutArray[0] = "006000040000500000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Rconv > 20 cm
+       
+} else if( trainConfig == 60 ) {
+      
+       eventCutArray[0]="8000011"; photonCutArray[0] = "007000040000500000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Rconv > 35 cm
+       
 }
        
 
@@ -479,7 +503,8 @@ if( trainConfig == 1 ) {  // No eta shift |Y| < 0.8
 
          if (  ( trainConfig >= 1 && trainConfig <= 6 ) || trainConfig == 9  ||  trainConfig == 11  || trainConfig == 13 || trainConfig == 14 || trainConfig == 16 || trainConfig == 17 || trainConfig == 18 || trainConfig == 19 || trainConfig == 20 || trainConfig == 21 || trainConfig == 22 || trainConfig == 23 ||
                  trainConfig == 28 || trainConfig == 29 || trainConfig == 30 ||  trainConfig == 31  || trainConfig == 32 || trainConfig == 33 || trainConfig == 37 || trainConfig == 38 || trainConfig == 39 || trainConfig == 40 || trainConfig == 41 || trainConfig == 41 || trainConfig == 43 || trainConfig == 44 ||
-                 trainConfig == 45 || trainConfig == 46 || trainConfig == 47 ||  trainConfig == 48  || trainConfig == 49 || trainConfig == 50 || trainConfig == 51 || trainConfig == 52 || trainConfig == 54 ) {
+                 trainConfig == 45 || trainConfig == 46 || trainConfig == 47 ||  trainConfig == 48  || trainConfig == 49 || trainConfig == 50 || trainConfig == 51 || trainConfig == 52 || trainConfig == 54 || trainConfig == 55 || trainConfig == 56 || trainConfig == 57 || trainConfig == 58 || trainConfig == 59 || 
+                 trainConfig == 60 ) {
            
            if (doWeighting){
              if (generatorName.CompareTo("DPMJET")==0){