Int_t clsId_ESD = track->GetEMCALcluster();
Int_t clsId = -9999;
- if(phi>1.0 && phi<3.5)clsId = TrMatch(track);
+ if(phi>1.0 && phi<3.5 && fabs(eta)<0.75)clsId = TrMatch(track);
//cout << "match ID" << clsId << " ; " << clsId_Tender << endl;
Double_t trkPos[3];
int matchID = -1;
Double_t MaxR = 0.05;
+ Double_t MaxEta = 0.025;
+ Double_t MaxPhi = 0.05;
Double_t magF = fESD->GetMagneticField();
Double_t magSign = 1.0;
Double_t delEmcphi = clsPosVec.Phi()-trkPosVec.Phi(); // track cluster matching
Double_t delEmceta = clsPosVec.Eta()-trkPosVec.Eta(); // track cluster matching
+ if(fabs(delEmcphi)>0.05)continue;
+ if(fabs(delEmceta)>0.025)continue;
+
double rmatch = sqrt(pow(delEmcphi,2)+pow(delEmceta,2));
- if(rmatch<MaxR)
+ //if(rmatch<MaxR)
+ if(fabs(delEmcphi)<MaxPhi && fabs(delEmceta)<MaxEta )
{
matchID = icl;
- MaxR = rmatch;
+ //MaxR = rmatch;
+ MaxPhi = fabs(delEmcphi);
+ MaxEta = fabs(delEmceta);
}
}
else if(configIndex==3) hfecuts->SetMinNClustersTPC(80);
else if(configIndex==4) hfecuts->SetMinNClustersTPC(85);
else if(configIndex==5) hfecuts->SetMinNClustersTPC(115);
- else if(configIndex==6) hfecuts->SetMinNClustersTPC(120); //Minimum number of clusters on TPC
+ else if(configIndex==6) hfecuts->SetMinNClustersTPC(120);
+ else if(configIndex==80) hfecuts->SetMinNClustersTPC(130);
+ else if(configIndex==81) hfecuts->SetMinNClustersTPC(140);//Minimum number of clusters on TPC
else hfecuts->SetMinNClustersTPC(100); //Minimum number of clusters on TPC
if(configIndex==7) hfecuts->SetMinNClustersTPCPID(70);
//ITS
if(configIndex==13) hfecuts->SetCutITSpixel(AliHFEextraCuts::kBoth); //Require at least one cluster on SPD
+ else if(configIndex==82) hfecuts->SetCutITSpixel(AliHFEextraCuts::kFirst);
else hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny); //Require at least one cluster on SPD
//hfecuts->SetCutITSdrift(AliHFEextraCuts::kAny); //Require at least one cluster on SDD
hfecuts->SetCheckITSLayerStatus(kFALSE);
//DCA cut included in the analysis 12 March 2014
- if(configIndex==18) hfecuts->SetMaxImpactParam(2,4);
- else if (configIndex==19) hfecuts->SetMaxImpactParam(0.5,1);
+ if(configIndex==18) hfecuts->SetMaxImpactParam(2,5);
+ else if(configIndex==19) hfecuts->SetMaxImpactParam(0.5,1);
+ else if(configIndex==83) hfecuts->SetMaxImpactParam(0.1,0.2);
else hfecuts->SetMaxImpactParam(1,2); //DCA to vertex
//Event Selection
//partner cuts
if(configIndex==29) task->SetAdditionalCuts(0.3,80);
- if(configIndex==30) task->SetAdditionalCuts(0.5,80);
- if(configIndex==31) task->SetAdditionalCuts(0.7,80);
- if(configIndex==32) task->SetAdditionalCuts(0.9,80);
+ else if(configIndex==30) task->SetAdditionalCuts(0.5,80);
+ else if(configIndex==31) task->SetAdditionalCuts(0.7,80);
+ else if(configIndex==32) task->SetAdditionalCuts(0.9,80);
- if(configIndex==33) task->SetAdditionalCuts(0,60);
- if(configIndex==34) task->SetAdditionalCuts(0,70);
- if(configIndex==35) task->SetAdditionalCuts(0,90);
- if(configIndex==36) task->SetAdditionalCuts(0,100);
+ else if(configIndex==33) task->SetAdditionalCuts(0,60);
+ else if(configIndex==34) task->SetAdditionalCuts(0,70);
+ else if(configIndex==35) task->SetAdditionalCuts(0,90);
+ else if(configIndex==36) task->SetAdditionalCuts(0,100);
- task->SetAdditionalCuts(0,80);
+ else task->SetAdditionalCuts(0,80);
//eta cuts
///_______________________________________________________________________________________________________________
/// New configurations for random cuts -- March, 05, 2014 -- Values in the macro "Random_configurations.C"
-
+ /*
if (configIndex==80){
hfecuts->SetMinNClustersTPC(86);
hfecuts->SetMinNClustersTPCPID(76);
Double_t params[0]=-1.15;
pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
}
+ */
///_______________________________________________________________________________________________________________