/**************************************************************************
* 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 *
* *
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),
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),
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];
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]);
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() );
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
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
+
}
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){