Linear gradient inside radiators. Height in the radiator used to find temperature

[u/mrichter/AliRoot.git] / HMPID / AliHMPIDTracker.cxx

diff --git a/HMPID/AliHMPIDTracker.cxx b/HMPID/AliHMPIDTracker.cxx
index 528ea538ff27475375c49875977940e6ba9b7eb7..837164c91e310cfffe080795949c5e3851274f9b 100644 (file)
--- a/HMPID/AliHMPIDTracker.cxx
+++ b/HMPID/AliHMPIDTracker.cxx
@@ -256,13 +256,13 @@ Int_t AliHMPIDTracker::Recon(AliESDEvent *pEsd,TObjArray *pClus,TObjArray *pNmea
 
     FillResiduals(hmpTrk,bestHmpCluster,kFALSE);
  
 
     FillResiduals(hmpTrk,bestHmpCluster,kFALSE);
  

+    Int_t iRad     = pParam->Radiator(yRa);                                                      //evaluate the radiator involved

     
     //evaluate nMean
     if(tsRight){
      if(pNmean->GetEntries()==21) {                                                              //for backward compatibility
        nmean=((TF1*)pNmean->At(3*ipCh))->Eval(ts);                                               //C6F14 Nmean for this chamber
      } else {
     
     //evaluate nMean
     if(tsRight){
      if(pNmean->GetEntries()==21) {                                                              //for backward compatibility
        nmean=((TF1*)pNmean->At(3*ipCh))->Eval(ts);                                               //C6F14 Nmean for this chamber
      } else {

-       Int_t iRad     = pParam->Radiator(yRa);                                                   //evaluate the radiator involved

        if(iRad < 0) {
        nmean = -1;
        } else {
        if(iRad < 0) {
        nmean = -1;
        } else {


@@ -376,12 +376,14 @@ Int_t AliHMPIDTracker::ReconHiddenTrk(AliESDEvent *pEsd,TObjArray *pClus,TObjArr
 
     Double_t yRa = yMip;                                                                        //just an approx...
     Double_t nmean;
 
     Double_t yRa = yMip;                                                                        //just an approx...
     Double_t nmean;

+
+    Int_t iRad     = pParam->Radiator(yRa);                                                   //evaluate the radiator involved
+    

     //evaluate nMean
     if(tsRight){
     if(pNmean->GetEntries()==21) {                                                              //for backward compatibility
       nmean=((TF1*)pNmean->At(3*chMip))->Eval(ts);                                              //C6F14 Nmean for this chamber
     } else {
     //evaluate nMean
     if(tsRight){
     if(pNmean->GetEntries()==21) {                                                              //for backward compatibility
       nmean=((TF1*)pNmean->At(3*chMip))->Eval(ts);                                              //C6F14 Nmean for this chamber
     } else {

-      Int_t iRad     = pParam->Radiator(yRa);                                                   //evaluate the radiator involved

       if(iRad < 0) {
        nmean = -1;
       } else {
       if(iRad < 0) {
        nmean = -1;
       } else {