if(!recParam){
AliFatal("Reconstruction parameters for EMCAL not set!");
}
-
- fCutX = recParam->GetTrkCutX();
- fCutY = recParam->GetTrkCutY();
- fCutZ = recParam->GetTrkCutZ();
- fMaxDist = recParam->GetTrkCutR();
- fCutAngle = recParam->GetTrkCutAngle();
- fCutAlphaMin = recParam->GetTrkCutAlphaMin();
- fCutAlphaMax = recParam->GetTrkCutAlphaMax();
- fCutNITS = recParam->GetTrkCutNITS();
- fCutNTPC = recParam->GetTrkCutNTPC();
+ else{
+ fCutX = recParam->GetTrkCutX();
+ fCutY = recParam->GetTrkCutY();
+ fCutZ = recParam->GetTrkCutZ();
+ fMaxDist = recParam->GetTrkCutR();
+ fCutAngle = recParam->GetTrkCutAngle();
+ fCutAlphaMin = recParam->GetTrkCutAlphaMin();
+ fCutAlphaMax = recParam->GetTrkCutAlphaMax();
+ fCutNITS = recParam->GetTrkCutNITS();
+ fCutNTPC = recParam->GetTrkCutNTPC();
+ }
}
//
for (Int_t i = 0; i < nClusters; i++) {
AliEMCALRecPoint *cluster = (AliEMCALRecPoint*)clusters->At(i);
if (!cluster) continue;
- if (cluster->GetClusterType() != AliESDCaloCluster::kEMCALClusterv1) continue;
+ if (cluster->GetClusterType() != AliVCluster::kEMCALClusterv1) continue;
AliEMCALMatchCluster *matchCluster = new AliEMCALMatchCluster(i, cluster);
fClusters->AddLast(matchCluster);
}
//
// TEMP
- Bool_t isTrue = kFALSE;
+ //Bool_t isTrue = kFALSE;
// if (tr->GetSeedLabel() == cl->Label()) {
// isTrue = kTRUE;
// }
Double_t r=0.;
cout.setf(ios::fixed);
cout.precision(5);
- if (isTrue) cout << "Init : " << rt << ' ' << x << ' ' << y << ' ' << z << endl;
+ //if (isTrue) cout << "Init : " << rt << ' ' << x << ' ' << y << ' ' << z << endl;
for (i = 0; i < fNPropSteps; i++) {
r = rt + (rc - rt) * ((Double_t)(i+1)/(Double_t)fNPropSteps);
if (!tr->PropagateTo(r, x0, rho)){
return distance;
}
tr->GetXYZ(pos);
- if (isTrue) cout << "Step : " << r << ' ' << x << ' ' << y << ' ' << z << endl;
+ // if (isTrue) cout << "Step : " << r << ' ' << x << ' ' << y << ' ' << z << endl;
}
- if (isTrue) cout << "Clstr: " << rc << ' ' << cl->X() << ' ' << cl->Y() << ' ' << cl->Z() << endl;
+ //if (isTrue) cout << "Clstr: " << rc << ' ' << cl->X() << ' ' << cl->Y() << ' ' << cl->Z() << endl;
}
else {
// when no steps are used, no correction makes sense