| bae04e0a |
1 | class LMEECutLib { |
| 2 | |
| 3 | public: |
| 4 | static enum LMMECutSet { |
| 5 | |
| 6 | kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight, |
| 7 | kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight, |
| 8 | kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight, |
| 9 | kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight, |
| 10 | kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4, |
| 11 | kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4, |
| 12 | kPbPb2011_pidITSTPC_trkSPDfirst_3, // (cutset w/o pairing) |
| 13 | kPbPb2011_pidTPC_trkSPDfirst_3, // (cutset w/o pairing) |
| 14 | kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_2_loose, // (cutset w/o pairing) |
| 15 | kPbPb2011_pidITSTPCTOFif_trkSPDfirst_2_loose, // (cutset w/o pairing) |
| 16 | kPbPb2011_pidTPCTOF_trkSPDorSDD_2_loose, // (cutset w/o pairing) |
| 17 | kPbPb2011_pidTPCTOF_trkSPDfirst_2_loose, // (cutset w/o pairing) |
| 18 | kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1, |
| 19 | kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1, // Cutset for Technical Preliminaries for QM2014 (no prefilter used!) |
| 20 | kPbPb2011_pidTPCTOF_trkSPDorSDD_1, |
| 21 | kPbPb2011_pidTPCTOF_trkSPDfirst_1, |
| 22 | kPbPb2011PID_ITSTPCTOFif2, // (NO FULL CUTSET) |
| 23 | kPbPb2011PID_TPCTOF3, // (NO FULL CUTSET) |
| 24 | kPbPb2011TRK_SDDfirstSPDnone, // (NO FULL CUTSET) complimentary tracks, strictly without SPD, to be combined with others! |
| 25 | kPbPb2011TRK_SPDfirst, // (NO FULL CUTSET) main track selection, with SPD first |
| 26 | kPbPb2011TRK_SDDfirstSPDnone4cls, // (NO FULL CUTSET) complimentary tracks, strictly without SPD, to be combined with others! |
| 27 | kPbPb2011TRK_SPDfirst5cls, // (NO FULL CUTSET) main track selection, with SPD first |
| 28 | kPbPb2011_TPCITS_TOFif1, |
| 29 | kPbPb2011_TPCTOF_Semi2, // changed PairCutsAna from PhiV to OpeningAngle. prefilter cuts renewed (if applicable) |
| 30 | // following cutsets are not complete anymore! |
| 31 | kPbPb2011_TPCTOF_Semi1, // old prefilter cuts (leg & pair), some are confusing |
| 32 | kPbPb2011NoPID, // pairing disabled in config |
| 33 | kPbPb2011TPCandTOF, // this was the final one activated by Christoph! |
| 34 | kPbPb2011TPCandTOFHPT, |
| 35 | kPbPb2011TPC, //TOF required, more relaxed cut on TPC |
| 36 | kPbPb2011TPCandTOFwide, //TOF required, more relaxed cut on TPC |
| 37 | kPbPb2011TPCorTOF, |
| 38 | kpp2010TPCandTOF, |
| 39 | kpp2010TPCorTOF, |
| 40 | kCUTSETMAX |
| 41 | }; |
| 42 | |
| 43 | static enum LMMECentSel { |
| 44 | kPbPb2011Central, |
| 45 | kPbPb2011MidCentral, |
| 46 | kPbPb2011SemiCentral, |
| 47 | kPbPb2011Peripheral, |
| 48 | kCENTSELMAX |
| 49 | }; |
| 50 | |
| 51 | //char* LMEECutNames[kCUTSETMAX] = { "PbPb2011TPCandTOF","PbPb2011TPCorTOF"}; |
| 52 | |
| 53 | |
| 54 | LMEECutLib() {} |
| 55 | |
| 56 | AliDielectronEventCuts* GetEventCuts(Int_t cutSet); |
| 57 | AliAnalysisCuts* GetCentralityCuts(Int_t centSel); |
| 58 | AliDielectronTrackRotator* GetTrackRotator(Int_t cutSet); |
| 59 | AliDielectronMixingHandler* GetMixingHandler(Int_t cutSet); |
| 60 | |
| 61 | AliAnalysisCuts* GetPairCutsAna(Int_t cutSet, Bool_t tooglePC=kFALSE); |
| 62 | AliAnalysisCuts* GetPairCutsPre(Int_t cutSet); |
| 63 | |
| 64 | AliAnalysisCuts* GetPIDCutsAna(Int_t cutSet); |
| 65 | AliAnalysisCuts* GetPIDCutsPre(Int_t cutSet); |
| 66 | |
| 67 | AliAnalysisCuts* GetTrackCutsAna(Int_t cutSet); |
| 68 | AliAnalysisCuts* GetTrackCutsPre(Int_t cutSet); |
| 69 | AliAnalysisCuts* GetESDTrackCutsAna(Int_t cutSet); |
| 70 | |
| 71 | |
| 72 | // Note: event cuts are identical for all analysis 'cutDefinition's that run together! |
| 73 | // the selection is hardcoded in the AddTask, currently to 'kPbPb2011_TPCTOF_Semi1' |
| 74 | AliDielectronEventCuts* GetEventCuts(Int_t cutSet) { |
| 75 | AliDielectronEventCuts* eventCuts = 0x0; |
| 76 | switch (cutSet) { |
| 77 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 78 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 79 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 80 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 81 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 82 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 83 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 84 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 85 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_2_loose: |
| 86 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_2_loose: |
| 87 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_2_loose: |
| 88 | case kPbPb2011_pidTPCTOF_trkSPDfirst_2_loose: |
| 89 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 90 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 91 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 92 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 93 | case kPbPb2011_TPCITS_TOFif1: |
| 94 | case kPbPb2011_TPCTOF_Semi2: |
| 95 | case kPbPb2011_TPCTOF_Semi1: |
| 96 | case kPbPb2011NoPID: |
| 97 | case kPbPb2011TPCandTOF : |
| 98 | case kPbPb2011TPCandTOFHPT: |
| 99 | case kPbPb2011TPC : |
| 100 | case kPbPb2011TPCandTOFwide : |
| 101 | case kPbPb2011TPCorTOF : |
| 102 | case kpp2010TPCandTOF : |
| 103 | case kpp2010TPCorTOF : |
| 104 | //Basic Event Cuts for pp and Pb-Pb, additional cuts may be in the AddTask |
| 105 | eventCuts=new AliDielectronEventCuts("eventCuts","Vertex Track && |vtxZ|<10 && ncontrib>0"); |
| 106 | eventCuts->SetVertexType(AliDielectronEventCuts::kVtxSPD); // AOD |
| 107 | //eventCuts->SetVertexType(AliDielectronEventCuts::kVtxTPC); // AOD |
| 108 | // eventCuts->SetCentralityRange(0.0,80.0); |
| 109 | eventCuts->SetRequireVertex(); |
| 110 | eventCuts->SetMinVtxContributors(1); |
| 111 | eventCuts->SetVertexZ(-10.,10.); |
| 112 | break; |
| 113 | default: cout << "No Event Cut defined" << endl; |
| 114 | } |
| 115 | return eventCuts; |
| 116 | } |
| 117 | |
| 118 | |
| 119 | //Selection of relatively 'flat' centralities |
| 120 | AliAnalysisCuts* GetCentralityCuts(Int_t centSel) { |
| 121 | AliDielectronVarCuts* centCuts = 0x0; |
| 122 | switch (centSel) { |
| 123 | case kPbPb2011Central: |
| 124 | centCuts = new AliDielectronVarCuts("centCuts","CentralityPbPb2011Central"); |
| 125 | centCuts->AddCut(AliDielectronVarManager::kCentrality,0.,10.); |
| 126 | break; |
| 127 | case kPbPb2011MidCentral: |
| 128 | centCuts = new AliDielectronVarCuts("centCuts","CentralityPbPb2011MidCentral"); |
| 129 | centCuts->AddCut(AliDielectronVarManager::kCentrality,10.,20.); |
| 130 | break; |
| 131 | case kPbPb2011SemiCentral: |
| 132 | centCuts = new AliDielectronVarCuts("centCuts","CentralityPbPb2011SemiCentral"); |
| 133 | centCuts->AddCut(AliDielectronVarManager::kCentrality,20.,50.); |
| 134 | break; |
| 135 | case kPbPb2011Peripheral: |
| 136 | centCuts = new AliDielectronVarCuts("centCuts","CentralityPbPb2011Peripheral"); |
| 137 | centCuts->AddCut(AliDielectronVarManager::kCentrality,50.,90.); |
| 138 | break; |
| 139 | default: cout << "No Centrality selected" << endl; |
| 140 | } |
| 141 | return centCuts; |
| 142 | } |
| 143 | |
| 144 | |
| 145 | //Basic track rotator settings from J/Psi, more investigation needed |
| 146 | AliDielectronTrackRotator* GetTrackRotator(Int_t cutSet) { |
| 147 | AliDielectronTrackRotator* trackRotator = 0x0; |
| 148 | switch (cutSet) { |
| 149 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 150 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 151 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 152 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 153 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 154 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 155 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 156 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 157 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 158 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 159 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 160 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 161 | case kPbPb2011_TPCITS_TOFif1: |
| 162 | case kPbPb2011_TPCTOF_Semi2: |
| 163 | case kPbPb2011_TPCTOF_Semi1: |
| 164 | case kPbPb2011NoPID: |
| 165 | case kPbPb2011TPCandTOF : |
| 166 | case kPbPb2011TPCandTOFHPT: |
| 167 | case kPbPb2011TPC : |
| 168 | case kPbPb2011TPCandTOFwide : |
| 169 | case kPbPb2011TPCorTOF : |
| 170 | case kpp2010TPCandTOF : |
| 171 | case kpp2010TPCorTOF : |
| 172 | trackRotator = new AliDielectronTrackRotator(); |
| 173 | trackRotator->SetIterations(20); |
| 174 | trackRotator->SetConeAnglePhi(TMath::Pi()/180*165); |
| 175 | trackRotator->SetStartAnglePhi(TMath::Pi()); |
| 176 | break; |
| 177 | default: cout << "No Rotator defined" << endl; |
| 178 | } |
| 179 | return trackRotator; |
| 180 | } |
| 181 | |
| 182 | |
| 183 | AliDielectronMixingHandler* GetMixingHandler(Int_t cutSet) { |
| 184 | AliDielectronMixingHandler* mixingHandler = 0x0; |
| 185 | switch (cutSet) { |
| 186 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 187 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 188 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 189 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 190 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 191 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 192 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 193 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 194 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 195 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 196 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 197 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 198 | case kPbPb2011_TPCITS_TOFif1: |
| 199 | case kPbPb2011_TPCTOF_Semi2: |
| 200 | case kPbPb2011_TPCTOF_Semi1: |
| 201 | mixingHandler = new AliDielectronMixingHandler; |
| 202 | mixingHandler->AddVariable(AliDielectronVarManager::kZvPrim,"-10,-5,0,5,10"); |
| 203 | mixingHandler->AddVariable(AliDielectronVarManager::kCentrality,"0,5,10,20,30,50,80"); |
| 204 | // now using TPC event plane, uncorrected. (also, the old phi range was wrong, now same effective binning.) |
| 205 | mixingHandler->AddVariable(AliDielectronVarManager::kTPCrpH2uc, 6, TMath::Pi()/-2., TMath::Pi()/2.); |
| 206 | mixingHandler->SetDepth(15); |
| 207 | mixingHandler->SetMixType(AliDielectronMixingHandler::kAll); |
| 208 | break; |
| 209 | //[...] |
| 210 | default: cout << "No Mixer defined" << endl; |
| 211 | } |
| 212 | return mixingHandler; |
| 213 | } |
| 214 | |
| 215 | |
| 216 | |
| 217 | //Pair Cuts for Analysis step - take care of logic - inverted compared to other PairCuts!! |
| 218 | // cuts = SELECTION!!! |
| 219 | AliAnalysisCuts* GetPairCutsAna(Int_t cutSet, Bool_t togglePC /*=kFALSE*/) { |
| 220 | cout << " >>>>>>>>>>>>>>>>>>>>>> GetPairCutsAna() >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 221 | AliAnalysisCuts* pairCuts=0x0; |
| 222 | switch (cutSet) { |
| 223 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 224 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 225 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 226 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 227 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 228 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 229 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 230 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 231 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 232 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 233 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 234 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 235 | cout << "No Pair Cuts used - ok " << endl; // since 18.02.2014 |
| 236 | break; |
| 237 | |
| 238 | case kPbPb2011_TPCITS_TOFif1: |
| 239 | case kPbPb2011_TPCTOF_Semi2: |
| 240 | // AliDielectronVarCuts* pairCutsPhivGood =new AliDielectronVarCuts("pairCutsPhivGood","pairCutsPhivGood"); |
| 241 | // pairCutsPhivGood->AddCut(AliDielectronVarManager::kPhivPair, 0.0, 2.0); |
| 242 | AliDielectronVarCuts* pairCutsOpAngGood =new AliDielectronVarCuts("pairCutsOpAngGood","pairCutsOpAngGood"); |
| 243 | pairCutsOpAngGood->AddCut(AliDielectronVarManager::kOpeningAngle, 0.05, 999.); // in upgrade: 0.05 |
| 244 | AliDielectronVarCuts* pairCutsInvM =new AliDielectronVarCuts("pairCutsInvM","pairCutsInvM"); |
| 245 | pairCutsInvM->AddCut(AliDielectronVarManager::kM, 0.0, 0.02); // in upgrade: 0.01 |
| 246 | AliDielectronVarCuts* pairCutsInvMgood =new AliDielectronVarCuts("pairCutsInvMgood","pairCutsInvMgood"); |
| 247 | pairCutsInvMgood->AddCut(AliDielectronVarManager::kM, 0.02, 99999.); |
| 248 | |
| 249 | AliDielectronCutGroup* pairCutsCG =new AliDielectronCutGroup("pairCutsCG","pairCutsCG",AliDielectronCutGroup::kCompAND); |
| 250 | pairCutsCG->AddCut(pairCutsInvM); |
| 251 | pairCutsCG->AddCut(pairCutsOpAngGood); |
| 252 | // pairCutsCG->AddCut(pairCutsPhivGood); |
| 253 | |
| 254 | AliDielectronCutGroup* pairCutsCG2 =new AliDielectronCutGroup("pairCutsCG2","pairCutsCG2",AliDielectronCutGroup::kCompOR); |
| 255 | pairCutsCG2->AddCut(pairCutsInvMgood); |
| 256 | pairCutsCG2->AddCut(pairCutsCG); |
| 257 | pairCuts = pairCutsCG2; |
| 258 | break; |
| 259 | |
| 260 | case kPbPb2011_TPCTOF_Semi1: |
| 261 | //[...] // PhiV and InvMass |
| 262 | default: cout << "No Pair Cuts defined " << endl; |
| 263 | } |
| 264 | return pairCuts; |
| 265 | } |
| 266 | |
| 267 | |
| 268 | //Pair Cuts for PREFILTER step |
| 269 | // cuts = REJECTION!!! |
| 270 | AliAnalysisCuts* GetPairCutsPre(Int_t cutSet) { |
| 271 | cout << " >>>>>>>>>>>>>>>>>>>>>> GetPairCutsPre() >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 272 | AliAnalysisCuts* pairCuts=0x0; |
| 273 | switch (cutSet) { |
| 274 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 275 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 276 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 277 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 278 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 279 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 280 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 281 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 282 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 283 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 284 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 285 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 286 | case kPbPb2011_TPCITS_TOFif1: |
| 287 | case kPbPb2011_TPCTOF_Semi2: |
| 288 | AliDielectronVarCuts* pairCutsInvM =new AliDielectronVarCuts("pairCutsInvM","pairCutsInvM"); |
| 289 | pairCutsInvM->AddCut(AliDielectronVarManager::kM, 0.0, 0.02); // in upgrade: 0.01 |
| 290 | AliDielectronVarCuts* pairCutsOpAng =new AliDielectronVarCuts("pairCutsOpAng","pairCutsOpAng"); |
| 291 | pairCutsOpAng->AddCut(AliDielectronVarManager::kOpeningAngle, 0.0, 0.05); // in upgrade: 0.05 |
| 292 | |
| 293 | AliDielectronCutGroup* pairCutsCG =new AliDielectronCutGroup("pairCutsCG","pairCutsCG",AliDielectronCutGroup::kCompAND); |
| 294 | pairCutsCG->AddCut(pairCutsInvM); |
| 295 | pairCutsCG->AddCut(pairCutsOpAng); |
| 296 | //pairCutsCG->AddCut(pairCutsPhiv); |
| 297 | pairCuts = pairCutsCG; |
| 298 | break; |
| 299 | |
| 300 | case kPbPb2011_TPCTOF_Semi1: |
| 301 | //[...] // PhiV and InvMass |
| 302 | default: cout << "No Prefilter Pair Cuts defined " << endl; |
| 303 | } |
| 304 | return pairCuts; |
| 305 | } |
| 306 | |
| 307 | |
| 308 | |
| 309 | AliAnalysisCuts* GetPIDCutsAna(Int_t cutSet) { |
| 310 | cout << " >>>>>>>>>>>>>>>>>>>>>> GetPIDCutsAna() >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 311 | AliAnalysisCuts* pidCuts=0x0; |
| 312 | |
| 313 | //----------------------------------------------- |
| 314 | // Define different PID Cuts, that are used later |
| 315 | //----------------------------------------------- |
| 316 | // PID cuts depend on TPC_inner_p, if not specified |
| 317 | // PID cut ranges correspond to global momentum P |
| 318 | // check it again!!! |
| 319 | //----------------------------------------------- |
| 320 | |
| 321 | // |
| 322 | // |
| 323 | //TPC: electron inclusion asymmetric |
| 324 | // pion exclusion 3sigma |
| 325 | //TOF: electron inclusion 3sigma in region where p,K cross electrons in TPC |
| 326 | AliDielectronPID *pidTPCTOF_Semi1 = new AliDielectronPID("pidTPCTOF_Semi1","pidTPCTOF_Semi1"); |
| 327 | pidTPCTOF_Semi1->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -1.5, 3., 0. ,100., kFALSE); |
| 328 | pidTPCTOF_Semi1->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3. , 3., 0. ,100., kTRUE); |
| 329 | pidTPCTOF_Semi1->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -3. , 3., 0. ,1.7 , kFALSE); |
| 330 | // |
| 331 | // |
| 332 | // LOOSE PID TPC+TOF |
| 333 | AliDielectronPID *pidTPCTOF_Semi_LOOSE = new AliDielectronPID("pidTPCTOF_Semi_LOOSE","pidTPCTOF_Semi_LOOSE"); |
| 334 | pidTPCTOF_Semi_LOOSE->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-12. ,20. , 0. ,100., kFALSE); |
| 335 | pidTPCTOF_Semi_LOOSE->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -3. , 3. , 0. ,1.7 , kFALSE); |
| 336 | // |
| 337 | // |
| 338 | // PID TPC only |
| 339 | AliDielectronPID *pidTPC_3 = new AliDielectronPID("pidTPC_3","pidTPC_3"); |
| 340 | pidTPC_3->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -1.5, 3., 0. ,100., kFALSE); |
| 341 | pidTPC_3->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3. , 3., 0. ,100., kTRUE); |
| 342 | |
| 343 | |
| 344 | //TPC: electron inclusion asymmetric |
| 345 | // pion exclusion 3sigma |
| 346 | //ITS: electron inclusion asymmetric in region where p,K cross electrons in TPC |
| 347 | //TOF: electron inclusion 3sigma in similar region - BUT ONLY IF AVAILABLE |
| 348 | AliDielectronPID *pidTPCITS_TOFif1 = new AliDielectronPID("pidTPCITS_TOFif1","pidTPCITS_TOFif1"); |
| 349 | pidTPCITS_TOFif1->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -1.5, 3. , 0. ,100., kFALSE); |
| 350 | pidTPCITS_TOFif1->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3. , 3. , 0. ,100., kTRUE); |
| 351 | pidTPCITS_TOFif1->AddCut(AliDielectronPID::kITS,AliPID::kElectron, -4. , 1. , 0. ,1.5 , kFALSE); |
| 352 | pidTPCITS_TOFif1->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -3. , 3. , 0. ,1.7 , kFALSE, AliDielectronPID::kIfAvailable); |
| 353 | // |
| 354 | // |
| 355 | //TPC: electron inclusion asymmetric |
| 356 | // pion exclusion 3sigma |
| 357 | //ITS: electron inclusion asymmetric OVER FULL MOMENTUM RANGE |
| 358 | //TOF: electron inclusion 3sigma - BUT ONLY IF AVAILABLE |
| 359 | AliDielectronPID *pidTPCITS_TOFif2 = new AliDielectronPID("pidTPCITS_TOFif2","pidTPCITS_TOFif2"); |
| 360 | pidTPCITS_TOFif2->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -1.5, 3. , 0. ,100., kFALSE); |
| 361 | pidTPCITS_TOFif2->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3. , 3. , 0. ,100., kTRUE); |
| 362 | pidTPCITS_TOFif2->AddCut(AliDielectronPID::kITS,AliPID::kElectron, -4. , 1. , 0. ,100., kFALSE); |
| 363 | pidTPCITS_TOFif2->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -3. , 3. , 0. ,100., kFALSE, AliDielectronPID::kIfAvailable); |
| 364 | // |
| 365 | // |
| 366 | // LOOSE PID ITS+TPC+TOFif |
| 367 | AliDielectronPID *pidTPCITS_TOFif_LOOSE = new AliDielectronPID("pidTPCITS_TOFif_LOOSE","pidTPCITS_TOFif_LOOSE"); |
| 368 | pidTPCITS_TOFif_LOOSE->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-12. ,20. , 0. ,100., kFALSE); |
| 369 | pidTPCITS_TOFif_LOOSE->AddCut(AliDielectronPID::kITS,AliPID::kElectron,-10. ,20. , 0. ,100., kFALSE); |
| 370 | pidTPCITS_TOFif_LOOSE->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -3. , 3. , 0. ,100., kFALSE, AliDielectronPID::kIfAvailable); |
| 371 | // |
| 372 | // |
| 373 | // PID ITS+TPC |
| 374 | AliDielectronPID *pidTPCITS_3 = new AliDielectronPID("pidTPCITS_3","pidTPCITS_3"); |
| 375 | pidTPCITS_3->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -1.5, 3. , 0. ,100., kFALSE); |
| 376 | pidTPCITS_3->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3. , 3. , 0. ,100., kTRUE); |
| 377 | pidTPCITS_3->AddCut(AliDielectronPID::kITS,AliPID::kElectron, -4. , 1. , 0. ,100., kFALSE); |
| 378 | // |
| 379 | // |
| 380 | // tighter PID ITS+TPC+TOFif |
| 381 | // ITS only up to momentum where proton contamination is seen in TPC signal |
| 382 | AliDielectronPID *pidTPCITS_TOFif56 = new AliDielectronPID("pidTPCITS_TOFif56","pidTPCITS_TOFif56"); |
| 383 | pidTPCITS_TOFif56->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -1.5, 2.5, 0. ,100., kFALSE); |
| 384 | pidTPCITS_TOFif56->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3. , 3. , 0. ,100., kTRUE); |
| 385 | pidTPCITS_TOFif56->AddCut(AliDielectronPID::kITS,AliPID::kElectron, -4. , 0.5, 0. , 2., kFALSE); |
| 386 | pidTPCITS_TOFif56->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -2. , 2. , 0. ,100., kFALSE, AliDielectronPID::kIfAvailable); |
| 387 | |
| 388 | |
| 389 | //Apply ITS cuts (see Hongyan's talks): |
| 390 | //3 sigma inclusion of electrons in TPC |
| 391 | //3 sigma exclusion of pions in TPC |
| 392 | //3 sigma inclusion of electrons in ITS,TOF for p<1.5GeV, where p,K contamination |
| 393 | //TOF only IF available! |
| 394 | AliDielectronPID *pidTPCandITSTOF = new AliDielectronPID("pidTPCandITSTOF","pidTPCandITSTOF");//"TPC-TOF-ITS" |
| 395 | pidTPCandITSTOF->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -3.,3., 0.,100., kFALSE); |
| 396 | pidTPCandITSTOF->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3.,3., 0.,100., kTRUE); |
| 397 | pidTPCandITSTOF->AddCut(AliDielectronPID::kITS,AliPID::kElectron, -6.,3., 0.,1.5, kFALSE ); |
| 398 | pidTPCandITSTOF->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -3.,3., 0.,1.5, kFALSE,AliDielectronPID::kIfAvailable ); |
| 399 | |
| 400 | // eta range: |
| 401 | AliDielectronVarCuts *etaRange090 = new AliDielectronVarCuts("etaRange090","etaRange090"); |
| 402 | etaRange090->AddCut(AliDielectronVarManager::kEta, -0.90, 0.90); |
| 403 | AliDielectronVarCuts *etaRange084 = new AliDielectronVarCuts("etaRange084","etaRange084"); |
| 404 | etaRange084->AddCut(AliDielectronVarManager::kEta, -0.84, 0.84); |
| 405 | AliDielectronVarCuts *etaRange076 = new AliDielectronVarCuts("etaRange076","etaRange076"); |
| 406 | etaRange076->AddCut(AliDielectronVarManager::kEta, -0.76, 0.76); |
| 407 | // pt range: |
| 408 | AliDielectronVarCuts *ptRange400to3500 = new AliDielectronVarCuts("ptRange400to3500","ptRange400to3500"); |
| 409 | ptRange400to3500->AddCut(AliDielectronVarManager::kPt, .4, 3.5); |
| 410 | |
| 411 | |
| 412 | //----------------------------------------------- |
| 413 | // Now see what Config actually loads and assemble final cuts |
| 414 | //----------------------------------------------- |
| 415 | switch (cutSet) { |
| 416 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: // tighter "ITSTPCTOFif" PID |
| 417 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 418 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 419 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 420 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 421 | cgPIDCutsAna->AddCut(etaRange076); |
| 422 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 423 | cgPIDCutsAna->AddCut(pidTPCITS_TOFif56); |
| 424 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); |
| 425 | pidCuts = cgPIDCutsAna; |
| 426 | break; |
| 427 | |
| 428 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 429 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 430 | cgPIDCutsAna->AddCut(etaRange076); |
| 431 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 432 | cgPIDCutsAna->AddCut(pidTPCITS_3); |
| 433 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); |
| 434 | pidCuts = cgPIDCutsAna; |
| 435 | break; |
| 436 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 437 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 438 | cgPIDCutsAna->AddCut(etaRange090); |
| 439 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 440 | cgPIDCutsAna->AddCut(pidTPC_3); |
| 441 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); |
| 442 | pidCuts = cgPIDCutsAna; |
| 443 | break; |
| 444 | |
| 445 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_2_loose: // loose "ITSTPCTOFif" PID - for 2D contamination study |
| 446 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_2_loose: |
| 447 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 448 | cgPIDCutsAna->AddCut(etaRange076); |
| 449 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 450 | cgPIDCutsAna->AddCut(pidTPCITS_TOFif_LOOSE); |
| 451 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); |
| 452 | pidCuts = cgPIDCutsAna; |
| 453 | break; |
| 454 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_2_loose: |
| 455 | case kPbPb2011_pidTPCTOF_trkSPDfirst_2_loose: |
| 456 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 457 | cgPIDCutsAna->AddCut(etaRange090); |
| 458 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 459 | cgPIDCutsAna->AddCut(pidTPCTOF_Semi_LOOSE); |
| 460 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); |
| 461 | pidCuts = cgPIDCutsAna; |
| 462 | break; |
| 463 | |
| 464 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: // regular "ITSTPCTOFif" PID |
| 465 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 466 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 467 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 468 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 469 | cgPIDCutsAna->AddCut(etaRange076); |
| 470 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 471 | cgPIDCutsAna->AddCut(pidTPCITS_TOFif2); |
| 472 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); |
| 473 | pidCuts = cgPIDCutsAna; |
| 474 | break; |
| 475 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 476 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 477 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 478 | cgPIDCutsAna->AddCut(etaRange090); |
| 479 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 480 | cgPIDCutsAna->AddCut(pidTPCTOF_Semi1); |
| 481 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); |
| 482 | pidCuts = cgPIDCutsAna; |
| 483 | break; |
| 484 | |
| 485 | case kPbPb2011_TPCITS_TOFif1: |
| 486 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 487 | cgPIDCutsAna->AddCut(etaRange084); // was 0.84 -> not ideal |
| 488 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 489 | cgPIDCutsAna->AddCut(pidTPCITS_TOFif1); |
| 490 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); |
| 491 | pidCuts = cgPIDCutsAna; |
| 492 | break; |
| 493 | case kPbPb2011_TPCTOF_Semi2: |
| 494 | case kPbPb2011_TPCTOF_Semi1: |
| 495 | AliDielectronCutGroup* cgPIDCutsAna = new AliDielectronCutGroup("cgPIDCutsAna","cgPIDCutsAna",AliDielectronCutGroup::kCompAND); |
| 496 | cgPIDCutsAna->AddCut(etaRange084); |
| 497 | cgPIDCutsAna->AddCut(ptRange400to3500); |
| 498 | cgPIDCutsAna->AddCut(pidTPCTOF_Semi1); |
| 499 | cgPIDCutsAna->AddCut(GetTrackCutsAna(cutSet)); // for 'kPbPb2011_TPCTOF_Semi1', this was called in the Config |
| 500 | pidCuts = cgPIDCutsAna; |
| 501 | break; |
| 502 | //[...] |
| 503 | case kPbPb2011TPCandTOFHPT: |
| 504 | //test Hongyan's cut |
| 505 | AliDielectronCutGroup* cgSecondTrackFilterPIDTPC1 = new AliDielectronCutGroup("cgPIDTPC1","cgPIDTPC1",AliDielectronCutGroup::kCompAND); |
| 506 | cgSecondTrackFilterPIDTPC1->AddCut(etaRange084); |
| 507 | cgSecondTrackFilterPIDTPC1->AddCut(ptRange400to3500); |
| 508 | cgSecondTrackFilterPIDTPC1->AddCut(pidTPCandITSTOF); |
| 509 | pidCuts = cgSecondTrackFilterPIDTPC1; |
| 510 | break; |
| 511 | //[...] |
| 512 | default: cout << "No Analysis PID Cut defined " << endl; |
| 513 | } |
| 514 | return pidCuts; |
| 515 | } |
| 516 | |
| 517 | |
| 518 | //Make/Tighten track Cuts that are *NOT* already |
| 519 | //done in the AOD production |
| 520 | //**IMPORTANT**: For AODs, select FilterBit |
| 521 | //the method is ignored for ESDs |
| 522 | |
| 523 | AliAnalysisCuts* GetTrackCutsAna(Int_t cutSet) { |
| 524 | cout << " >>>>>>>>>>>>>>>>>>>>>> GetTrackCutsAna() >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 525 | AliDielectronCutGroup* trackCuts=0x0; |
| 526 | switch (cutSet) { |
| 527 | |
| 528 | |
| 529 | //---------- |
| 530 | // these MAIN settings have to combine different track selections: |
| 531 | //---------- |
| 532 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 533 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_2_loose: |
| 534 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_2_loose: |
| 535 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 536 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 537 | // combine typical and new trackcuts with "kCompOR" condition: |
| 538 | cgTrackCutsAnaSPDorSDD = new AliDielectronCutGroup("cgTrackCutsAnaSPDorSDD","cgTrackCutsAnaSPDorSDD",AliDielectronCutGroup::kCompOR); |
| 539 | cgTrackCutsAnaSPDorSDD->AddCut(GetTrackCutsAna(kPbPb2011TRK_SPDfirst)); // typical trackcuts with requirement of SPD |
| 540 | cgTrackCutsAnaSPDorSDD->AddCut(GetTrackCutsAna(kPbPb2011TRK_SDDfirstSPDnone)); // new additional trackcuts with SDD instead of SPD |
| 541 | trackCuts = cgTrackCutsAnaSPDorSDD; |
| 542 | break; |
| 543 | |
| 544 | //---------- |
| 545 | // these MAIN settings just load the main track selection directly below: |
| 546 | //---------- |
| 547 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 548 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 549 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 550 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_2_loose: |
| 551 | case kPbPb2011_pidTPCTOF_trkSPDfirst_2_loose: |
| 552 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 553 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 554 | //---------- |
| 555 | case kPbPb2011TRK_SPDfirst: // main track selection, now closer to what Hongyan does... |
| 556 | AliDielectronVarCuts* trackCutsAOD =new AliDielectronVarCuts("trackCutsAOD","trackCutsAOD"); |
| 557 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0); |
| 558 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0); |
| 559 | trackCutsAOD->AddCut(AliDielectronVarManager::kNclsITS, 4.0, 100.0); // means at least 2 with PID |
| 560 | trackCutsAOD->AddCut(AliDielectronVarManager::kTPCchi2Cl, 0.0, 4.0); |
| 561 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCr, 100.0, 160.0); |
| 562 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCfCross, 0.8, 1.1); // lower limit 0.8 in most filterbits! // 1.1 since 26.02.2014 |
| 563 | AliDielectronTrackCuts *trackCutsDiel = new AliDielectronTrackCuts("trackCutsDiel","trackCutsDiel"); |
| 564 | trackCutsDiel->SetAODFilterBit(1<<4); // (=16) filterbit 4! //GetStandardITSTPCTrackCuts2011(kFALSE); loose DCA, 2D cut |
| 565 | trackCutsDiel->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst); |
| 566 | |
| 567 | cgTrackCutsAnaSPDfirst = new AliDielectronCutGroup("cgTrackCutsAnaSPDfirst","cgTrackCutsAnaSPDfirst",AliDielectronCutGroup::kCompAND); |
| 568 | cgTrackCutsAnaSPDfirst->AddCut(trackCutsDiel); |
| 569 | cgTrackCutsAnaSPDfirst->AddCut(trackCutsAOD); |
| 570 | trackCuts = cgTrackCutsAnaSPDfirst; |
| 571 | break; |
| 572 | |
| 573 | case kPbPb2011TRK_SDDfirstSPDnone: // complimentary tracks, strictly without SPD, to be combined with others! |
| 574 | AliDielectronVarCuts* trackCutsAOD =new AliDielectronVarCuts("trackCutsAOD","trackCutsAOD"); |
| 575 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0); |
| 576 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0); |
| 577 | trackCutsAOD->AddCut(AliDielectronVarManager::kNclsITS, 3.0, 100.0); // means at least 3 with PID |
| 578 | trackCutsAOD->AddCut(AliDielectronVarManager::kTPCchi2Cl, 0.0, 4.0); |
| 579 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCr, 100.0, 160.0); |
| 580 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCfCross, 0.8, 1.1); // lower limit 0.8 in most filterbits! // 1.1 since 26.02.2014 |
| 581 | AliDielectronTrackCuts *trackCutsDiel = new AliDielectronTrackCuts("trackCutsDiel","trackCutsDiel"); |
| 582 | trackCutsDiel->SetAODFilterBit(1<<6); // GetStandardITSTPCTrackCuts2011(kTRUE), SPD none, SDD first |
| 583 | |
| 584 | cgTrackCutsAnaSDDnoSPD = new AliDielectronCutGroup("cgTrackCutsAnaSDDnoSPD","cgTrackCutsAnaSDDnoSPD",AliDielectronCutGroup::kCompAND); |
| 585 | cgTrackCutsAnaSDDnoSPD->AddCut(trackCutsDiel); |
| 586 | cgTrackCutsAnaSDDnoSPD->AddCut(trackCutsAOD); |
| 587 | trackCuts = cgTrackCutsAnaSDDnoSPD; |
| 588 | break; |
| 589 | |
| 590 | //---------- |
| 591 | // MAIN settings - combined trackset - variation 1 |
| 592 | //---------- |
| 593 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 594 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 595 | // combine typical and new trackcuts with "kCompOR" condition: |
| 596 | cgTrackCutsAnaSPDorSDD = new AliDielectronCutGroup("cgTrackCutsAnaSPDorSDD","cgTrackCutsAnaSPDorSDD",AliDielectronCutGroup::kCompOR); |
| 597 | cgTrackCutsAnaSPDorSDD->AddCut(GetTrackCutsAna(kPbPb2011TRK_SPDfirst5cls)); // typical trackcuts with requirement of SPD |
| 598 | cgTrackCutsAnaSPDorSDD->AddCut(GetTrackCutsAna(kPbPb2011TRK_SDDfirstSPDnone4cls)); // new additional trackcuts with SDD instead of SPD |
| 599 | trackCuts = cgTrackCutsAnaSPDorSDD; |
| 600 | break; |
| 601 | |
| 602 | //---------- |
| 603 | // MAIN settings - single trackset - variation 1 |
| 604 | //---------- |
| 605 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 606 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 607 | //---------- |
| 608 | case kPbPb2011TRK_SPDfirst5cls: // main track selection, 5+ ITS clusters |
| 609 | AliDielectronVarCuts* trackCutsAOD =new AliDielectronVarCuts("trackCutsAOD","trackCutsAOD"); |
| 610 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0); |
| 611 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0); |
| 612 | trackCutsAOD->AddCut(AliDielectronVarManager::kNclsITS, 5.0, 100.0); // means at least 3 with PID |
| 613 | trackCutsAOD->AddCut(AliDielectronVarManager::kTPCchi2Cl, 0.0, 4.0); |
| 614 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCr, 100.0, 160.0); |
| 615 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCfCross, 0.8, 1.1); |
| 616 | AliDielectronTrackCuts *trackCutsDiel = new AliDielectronTrackCuts("trackCutsDiel","trackCutsDiel"); |
| 617 | trackCutsDiel->SetAODFilterBit(1<<4); // (=16) filterbit 4! //GetStandardITSTPCTrackCuts2011(kFALSE); loose DCA, 2D cut |
| 618 | trackCutsDiel->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst); |
| 619 | |
| 620 | cgTrackCutsAnaSPDfirst = new AliDielectronCutGroup("cgTrackCutsAnaSPDfirst","cgTrackCutsAnaSPDfirst",AliDielectronCutGroup::kCompAND); |
| 621 | cgTrackCutsAnaSPDfirst->AddCut(trackCutsDiel); |
| 622 | cgTrackCutsAnaSPDfirst->AddCut(trackCutsAOD); |
| 623 | trackCuts = cgTrackCutsAnaSPDfirst; |
| 624 | break; |
| 625 | |
| 626 | case kPbPb2011TRK_SDDfirstSPDnone4cls: // complimentary tracks, 4+ ITS clusters, strictly without SPD, to be combined with others! |
| 627 | AliDielectronVarCuts* trackCutsAOD =new AliDielectronVarCuts("trackCutsAOD","trackCutsAOD"); |
| 628 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0); |
| 629 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0); |
| 630 | trackCutsAOD->AddCut(AliDielectronVarManager::kNclsITS, 4.0, 100.0); // means at least 4 with PID |
| 631 | trackCutsAOD->AddCut(AliDielectronVarManager::kTPCchi2Cl, 0.0, 4.0); |
| 632 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCr, 100.0, 160.0); |
| 633 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCfCross, 0.8, 1.1); |
| 634 | AliDielectronTrackCuts *trackCutsDiel = new AliDielectronTrackCuts("trackCutsDiel","trackCutsDiel"); |
| 635 | trackCutsDiel->SetAODFilterBit(1<<6); // GetStandardITSTPCTrackCuts2011(kTRUE), SPD none, SDD first |
| 636 | |
| 637 | cgTrackCutsAnaSDDnoSPD = new AliDielectronCutGroup("cgTrackCutsAnaSDDnoSPD","cgTrackCutsAnaSDDnoSPD",AliDielectronCutGroup::kCompAND); |
| 638 | cgTrackCutsAnaSDDnoSPD->AddCut(trackCutsDiel); |
| 639 | cgTrackCutsAnaSDDnoSPD->AddCut(trackCutsAOD); |
| 640 | trackCuts = cgTrackCutsAnaSDDnoSPD; |
| 641 | break; |
| 642 | |
| 643 | // ========== |
| 644 | case kPbPb2011_TPCITS_TOFif1: |
| 645 | case kPbPb2011_TPCTOF_Semi2: // no pt and eta ranges in the trackcuts anymore! |
| 646 | AliDielectronVarCuts* trackCutsAOD =new AliDielectronVarCuts("trackCutsAOD","trackCutsAOD"); |
| 647 | // trackCutsAOD->AddCut(AliDielectronVarManager::kEta, -0.84, 0.84); // eta commented out later. (05.02.2014) |
| 648 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0); |
| 649 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0); |
| 650 | trackCutsAOD->AddCut(AliDielectronVarManager::kNclsITS, 3.0, 100.0); |
| 651 | trackCutsAOD->AddCut(AliDielectronVarManager::kTPCchi2Cl, 0.0, 3.5); |
| 652 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCr, 110.0, 160.0); |
| 653 | trackCutsAOD->AddCut(AliDielectronVarManager::kNFclsTPCfCross, 0.8, 1.0); |
| 654 | AliDielectronTrackCuts *trackCutsDiel = new AliDielectronTrackCuts("trackCutsDiel","trackCutsDiel"); |
| 655 | trackCutsDiel->SetAODFilterBit(16); //does nothing for ESDs // 16=2^4 -> filter bit 4! |
| 656 | trackCutsDiel->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst); //function in AliDielectronTrackCuts |
| 657 | |
| 658 | cgTrackCutsAna = new AliDielectronCutGroup("cgTrackCutsAna","cgTrackCutsAna",AliDielectronCutGroup::kCompAND); |
| 659 | cgTrackCutsAna->AddCut(trackCutsDiel); |
| 660 | cgTrackCutsAna->AddCut(trackCutsAOD); |
| 661 | trackCuts = cgTrackCutsAna; |
| 662 | break; |
| 663 | |
| 664 | case kPbPb2011_TPCTOF_Semi1: |
| 665 | //[...] |
| 666 | default: cout << "No Analysis Track Cut defined " << endl; |
| 667 | } |
| 668 | return trackCuts; |
| 669 | } |
| 670 | |
| 671 | |
| 672 | |
| 673 | //Relaxed PID cuts for additional rejectin step, do not use blindly |
| 674 | AliAnalysisCuts* GetPIDCutsPre(Int_t cutSet) { |
| 675 | cout << " >>>>>>>>>>>>>>>>>>>>>> GetPIDCutsPre() >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 676 | AliAnalysisCuts* pidCuts=0x0; |
| 677 | switch (cutSet) { |
| 678 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 679 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 680 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_2_loose: |
| 681 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_2_loose: |
| 682 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_2_loose: |
| 683 | case kPbPb2011_pidTPCTOF_trkSPDfirst_2_loose: |
| 684 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 685 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 686 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 687 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 688 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 689 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 690 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 691 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 692 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 693 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 694 | case kPbPb2011_TPCITS_TOFif1: |
| 695 | case kPbPb2011_TPCTOF_Semi2: |
| 696 | |
| 697 | // eta range: |
| 698 | AliDielectronVarCuts *etaRangePre1 = new AliDielectronVarCuts("etaRangePre1","etaRangePre1"); |
| 699 | etaRangePre1->AddCut(AliDielectronVarManager::kEta,-0.9,0.9); |
| 700 | // pt range: |
| 701 | AliDielectronVarCuts *ptRangePre1 = new AliDielectronVarCuts("ptRangePre1","ptRangePre1"); |
| 702 | ptRangePre1->AddCut(AliDielectronVarManager::kPt, .2, 3.5); // 0.2 is realistic. turnon at ~180MeV |
| 703 | //AliDielectronVarCuts *ptRangePre2 = new AliDielectronVarCuts("ptRangePre2","ptRangePre2"); |
| 704 | //ptRangePre2->AddCut(AliDielectronVarManager::kPt, .4, 3.5); |
| 705 | //AliDielectronVarCuts *ptRangePre3 = new AliDielectronVarCuts("ptRangePre3","ptRangePre3"); |
| 706 | //ptRangePre3->AddCut(AliDielectronVarManager::kPt, 0.05, 1.5); |
| 707 | |
| 708 | AliDielectronCutGroup* cgITSTPCTOFpre = new AliDielectronCutGroup("cgITSTPCTOFpre","cgITSTPCTOFpre",AliDielectronCutGroup::kCompAND); |
| 709 | AliDielectronPID *pidITSTPCTOFpre = new AliDielectronPID("pidITSTPCTOFpre","pidITSTPCTOFpre"); |
| 710 | pidITSTPCTOFpre->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -3. , 3., 0. ,100., kFALSE); |
| 711 | pidITSTPCTOFpre->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3. , 3., 0. ,100., kTRUE); |
| 712 | // ITS will be used: |
| 713 | // pidITSTPCTOFpre->AddCut(AliDielectronPID::kITS,AliPID::kElectron, -3. , 3., 0. ,1.7 , kFALSE); |
| 714 | // TOF will be used if available, and with pt instead of p: |
| 715 | // pidITSTPCTOFpre->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -3. , 3., 0.4,100., kFALSE, |
| 716 | // AliDielectronPID::kIfAvailable, AliDielectronVarManager::kPt); |
| 717 | cgITSTPCTOFpre->AddCut(pidITSTPCTOFpre); |
| 718 | cgITSTPCTOFpre->AddCut(etaRangePre1); |
| 719 | cgITSTPCTOFpre->AddCut(ptRangePre1); |
| 720 | cgITSTPCTOFpre->AddCut(GetTrackCutsAna(cutSet)); |
| 721 | |
| 722 | // AliDielectronCutGroup* cgTPCpre = new AliDielectronCutGroup("cgTPCpre","cgTPCpre",AliDielectronCutGroup::kCompAND); |
| 723 | // AliDielectronPID *pidTPCpre = new AliDielectronPID("pidTPCpre","pidTPCpre"); |
| 724 | // pidTPCpre->AddCut(AliDielectronPID::kTPC,AliPID::kElectron, -3. , 3., 0. ,100., kFALSE); |
| 725 | // pidTPCpre->AddCut(AliDielectronPID::kTPC,AliPID::kPion, -3. , 3., 0. ,100., kTRUE); |
| 726 | // // TOF will be used if available, and with pt instead of p: |
| 727 | // pidTPCpre->AddCut(AliDielectronPID::kTOF,AliPID::kElectron, -3. , 3., 0.4,5. , kFALSE, |
| 728 | // AliDielectronPID::kIfAvailable, AliDielectronVarManager::kPt); |
| 729 | // cgTPCpre->AddCut(pidTPCpre); |
| 730 | // cgITSTPCTOFpre->AddCut(etaRangePre1); |
| 731 | // cgTPCpre->AddCut(ptRangePre2); |
| 732 | // cgTPCpre->AddCut(GetTrackCutsAna(cutSet)); |
| 733 | |
| 734 | // AliDielectronCutGroup* cgITSSA = new AliDielectronCutGroup("cgITSSA","cgITSSA",AliDielectronCutGroup::kCompAND); |
| 735 | // AliDielectronPID *pidITSSA = new AliDielectronPID("pidITSSA","pidITSSA"); |
| 736 | // pidITSSA->AddCut(AliDielectronPID::kITS,AliPID::kElectron, -3., 3.); |
| 737 | // // this means that very many pions will be used for rejection! |
| 738 | // cgITSSA->AddCut(pidITSSA); |
| 739 | // cgITSTPCTOFpre->AddCut(etaRangePre1); |
| 740 | // cgITSSA->AddCut(ptRangePre3); |
| 741 | // cgITSSA->AddCut(GetTrackCutsPre(cutSet)); |
| 742 | |
| 743 | AliDielectronCutGroup* cgInitialTrackFilter = new AliDielectronCutGroup("cgInitialTrackFilter","cgInitialTrackFilter",AliDielectronCutGroup::kCompOR); |
| 744 | cgInitialTrackFilter->AddCut(GetPIDCutsAna(cutSet)); // in case the prefilter cuts do not include all needed global tracks. |
| 745 | cgInitialTrackFilter->AddCut(cgITSTPCTOFpre); |
| 746 | //cgInitialTrackFilter->AddCut(cgTPCpre); |
| 747 | //cgInitialTrackFilter->AddCut(cgITSSA); |
| 748 | pidCuts = cgInitialTrackFilter; // kCompOR works!!! <- checked with 'SetNoPairing()' and commented out 'GetPIDCutsAna(selectedPID)' |
| 749 | //cout << " ========== pidCuts prefilter: ========== " << endl; |
| 750 | //pidCuts->Print(); |
| 751 | break; |
| 752 | |
| 753 | case kPbPb2011_TPCTOF_Semi1: |
| 754 | //[...] |
| 755 | default: cout << "No Prefilter PID Cut defined " << endl; |
| 756 | } |
| 757 | return pidCuts; |
| 758 | } |
| 759 | |
| 760 | |
| 761 | //Possibly different cut sets for Prefilter step |
| 762 | //Not used at the moment |
| 763 | AliAnalysisCuts* GetTrackCutsPre(Int_t cutSet) { |
| 764 | cout << " >>>>>>>>>>>>>>>>>>>>>> GetTrackCutsPre() >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 765 | AliDielectronCutGroup* trackCuts=0x0; |
| 766 | switch (cutSet) { |
| 767 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 768 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 769 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_2_loose: |
| 770 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_2_loose: |
| 771 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_2_loose: |
| 772 | case kPbPb2011_pidTPCTOF_trkSPDfirst_2_loose: |
| 773 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 774 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 775 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 776 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 777 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 778 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 779 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 780 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 781 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 782 | AliDielectronVarCuts* trackCutsAOD =new AliDielectronVarCuts("trackCutsAOD","trackCutsAOD"); |
| 783 | trackCutsAOD->AddCut(AliDielectronVarManager::kEta,-0.84,0.84); |
| 784 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0); |
| 785 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0); |
| 786 | trackCutsAOD->AddCut(AliDielectronVarManager::kNclsITS, 3.0, 100.0); |
| 787 | AliDielectronTrackCuts *trackCutsDiel = new AliDielectronTrackCuts("trackCutsDiel","trackCutsDiel"); |
| 788 | trackCutsDiel->SetAODFilterBit(1); //does nothing for ESDs, ITSSA(???) // maybe use FilterBit(2) instead! |
| 789 | |
| 790 | cgTrackCutsPre = new AliDielectronCutGroup("cgTrackCutsPre","cgTrackCutsPre",AliDielectronCutGroup::kCompAND); |
| 791 | cgTrackCutsPre->AddCut(trackCutsDiel); |
| 792 | cgTrackCutsPre->AddCut(trackCutsAOD); |
| 793 | trackCuts = cgTrackCutsPre; |
| 794 | break; |
| 795 | |
| 796 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 797 | AliDielectronVarCuts* trackCutsAOD =new AliDielectronVarCuts("trackCutsAOD","trackCutsAOD"); |
| 798 | trackCutsAOD->AddCut(AliDielectronVarManager::kEta,-0.9,0.9); |
| 799 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0); |
| 800 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0); |
| 801 | trackCutsAOD->AddCut(AliDielectronVarManager::kNclsITS, 4.0, 100.0); |
| 802 | AliDielectronTrackCuts *trackCutsDiel = new AliDielectronTrackCuts("trackCutsDiel","trackCutsDiel"); |
| 803 | trackCutsDiel->SetAODFilterBit(1); //does nothing for ESDs, ITSSA(???) // maybe use FilterBit(2) instead! |
| 804 | trackCutsDiel->SetRequireITSRefit(kTRUE); //function in AliDielectronTrackCuts |
| 805 | trackCutsDiel->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst); //function in AliDielectronTrackCuts |
| 806 | |
| 807 | cgTrackCutsPre = new AliDielectronCutGroup("cgTrackCutsPre","cgTrackCutsPre",AliDielectronCutGroup::kCompAND); |
| 808 | cgTrackCutsPre->AddCut(trackCutsDiel); |
| 809 | cgTrackCutsPre->AddCut(trackCutsAOD); |
| 810 | trackCuts = cgTrackCutsPre; |
| 811 | break; |
| 812 | |
| 813 | case kPbPb2011_TPCITS_TOFif1: |
| 814 | case kPbPb2011_TPCTOF_Semi2: // no pt ranges in the trackcuts anymore! |
| 815 | AliDielectronVarCuts* trackCutsAOD =new AliDielectronVarCuts("trackCutsAOD","trackCutsAOD"); |
| 816 | trackCutsAOD->AddCut(AliDielectronVarManager::kEta,-0.84,0.84); |
| 817 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0); |
| 818 | trackCutsAOD->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0); |
| 819 | trackCutsAOD->AddCut(AliDielectronVarManager::kNclsITS, 3.0, 100.0); |
| 820 | AliDielectronTrackCuts *trackCutsDiel = new AliDielectronTrackCuts("trackCutsDiel","trackCutsDiel"); |
| 821 | trackCutsDiel->SetAODFilterBit(1); //does nothing for ESDs, ITSSA(???) // maybe use FilterBit(2) instead! |
| 822 | trackCutsDiel->SetRequireITSRefit(kTRUE); //function in AliDielectronTrackCuts |
| 823 | trackCutsDiel->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst); //function in AliDielectronTrackCuts |
| 824 | |
| 825 | cgTrackCutsPre = new AliDielectronCutGroup("cgTrackCutsPre","cgTrackCutsPre",AliDielectronCutGroup::kCompAND); |
| 826 | cgTrackCutsPre->AddCut(trackCutsDiel); |
| 827 | cgTrackCutsPre->AddCut(trackCutsAOD); |
| 828 | trackCuts = cgTrackCutsPre; |
| 829 | break; |
| 830 | |
| 831 | case kPbPb2011_TPCTOF_Semi1: |
| 832 | //[...] |
| 833 | default: cout << "No Prefilter Track Cut defined " << endl; |
| 834 | } |
| 835 | return trackCuts; |
| 836 | } |
| 837 | |
| 838 | |
| 839 | |
| 840 | //******************************************************************************* |
| 841 | //******************************************************************************* |
| 842 | //** ESD TRACK CUTS TUNED FOR AGREEMENT BETWEEN AODS AND ESDS ****************** |
| 843 | //** NOT NECESSARILY 100% OPTIMIZED FOR DIEL-ANALYSIS ****************** |
| 844 | //******************************************************************************* |
| 845 | //******************************************************************************* |
| 846 | |
| 847 | //WHEN RUNNING ON ESDs: LOAD Default Cuts for AODs |
| 848 | AliAnalysisCuts* GetESDTrackCutsAna(Int_t cutSet) { |
| 849 | //cout << " >>>>>>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 850 | cout << " >>>>>>>>>>>>>>>>>>>>>> GetESDTrackCutsAna() >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 851 | //cout << " >>>>>>>>>>>>>>>>>>>>>> Setting ESD Track Cuts >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 852 | //cout << " >>>>>>>>>>>>>>>>>>>>>> ( do we run on ESD?! ) >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 853 | //cout << " >>>>>>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>>>>> " << endl; |
| 854 | AliESDtrackCuts* esdTrackCutsH = 0x0; |
| 855 | switch (cutSet) { |
| 856 | case kPbPb2011_pidITSTPC_trkSPDfirst_3: |
| 857 | case kPbPb2011_pidTPC_trkSPDfirst_3: |
| 858 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_2_loose: |
| 859 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_2_loose: |
| 860 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_2_loose: |
| 861 | case kPbPb2011_pidTPCTOF_trkSPDfirst_2_loose: |
| 862 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPD5orSDD4cls_6_tight: |
| 863 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst5cls_6_tight: |
| 864 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDorSDD_5_tight: |
| 865 | case kPbPb2011_pidITS2gevTPCTOFif_trkSPDfirst_5_tight: |
| 866 | case kPbPb2011_pidITSTPCTOFif_trkSPD5orSDD4cls_4: |
| 867 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst5cls_4: |
| 868 | case kPbPb2011_pidITSTPCTOFif_trkSPDorSDD_1: |
| 869 | case kPbPb2011_pidITSTPCTOFif_trkSPDfirst_1: |
| 870 | case kPbPb2011_pidTPCTOF_trkSPDorSDD_1: |
| 871 | case kPbPb2011_pidTPCTOF_trkSPDfirst_1: |
| 872 | case kPbPb2011_TPCITS_TOFif1: |
| 873 | case kPbPb2011_TPCTOF_Semi2: |
| 874 | case kPbPb2011_TPCTOF_Semi1: |
| 875 | case kPbPb2011NoPID: |
| 876 | case kPbPb2011TPCandTOF : |
| 877 | case kPbPb2011TPCandTOFHPT: |
| 878 | case kPbPb2011TPC : |
| 879 | case kPbPb2011TPCandTOFwide : |
| 880 | case kPbPb2011TPCorTOF : |
| 881 | case kpp2010TPCandTOF : |
| 882 | case kpp2010TPCorTOF : |
| 883 | // standard cuts with very loose DCA: Bit4 (Int: 16), AOD095&115 |
| 884 | |
| 885 | esdTrackCutsH = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kFALSE); |
| 886 | esdTrackCutsH->SetMaxDCAToVertexXY(2.4); |
| 887 | esdTrackCutsH->SetMaxDCAToVertexZ(3.2); |
| 888 | esdTrackCutsH->SetDCAToVertex2D(kTRUE); |
| 889 | |
| 890 | //The cuts below should be the onyl ones that are missing |
| 891 | //explicitely in the TrackCutsAna method |
| 892 | //To be sure, StandardITSTPCTrackCuts is loaded however |
| 893 | /* |
| 894 | esdTrackCutsH = new AliESDtrackCuts(); |
| 895 | esdTrackCutsH->SetAcceptKinkDaughters(kFALSE); |
| 896 | //Not done so far via dielectron cuts: |
| 897 | */ |
| 898 | /* |
| 899 | esdTrackCuts->SetDCAToVertex2D(kFALSE); |
| 900 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); |
| 901 | esdTrackCuts->SetMaxChi2PerClusterITS(36); |
| 902 | */ |
| 903 | |
| 904 | break; |
| 905 | default: cout << "No ESD Track Cut defined " << endl; |
| 906 | } |
| 907 | return esdTrackCutsH; |
| 908 | } |
| 909 | |
| 910 | |
| 911 | }; |
git://git.uio.no / u / mrichter / AliRoot.git / blame