/[PAMELA software]/calo/flight/CaloBragg/src/CaloBragg.cpp
ViewVC logotype

Diff of /calo/flight/CaloBragg/src/CaloBragg.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.14 by mocchiut, Thu Oct 6 11:52:57 2011 UTC revision 1.15 by biancoa, Wed Nov 23 14:33:32 2011 UTC
# Line 565  void CaloBragg::LoadParam(){ Line 565  void CaloBragg::LoadParam(){
565    
566    elem[3][0] = 9.012182;  //Be 4    elem[3][0] = 9.012182;  //Be 4
567    elem[3][1] = 10.01353;  //10Be (Isotope) (most stable)    elem[3][1] = 10.01353;  //10Be (Isotope) (most stable)
568    elem[3][2] = -1.;    elem[3][2] = 7.01693;   //9Be no EC in space?
569    elem[3][3] = -1.;    elem[3][3] = -1.;
570    elem[3][4] = -1.;    elem[3][4] = -1.;
571    elem[3][5] = -1.;    elem[3][5] = -1.;
572    elem[3][6] = -1.;    elem[3][6] = -1.;
573    
574    elem[4][0] = 11.0093;   //B 5    elem[4][0] = 11.00930;   //B 5
575    elem[4][1] = 10.01294;  //10B  (Isotope)    elem[4][1] = 10.01294;  //10B  (Isotope)
576    elem[4][2] = -1.;    elem[4][2] = -1.;
577    elem[4][3] = -1.;    elem[4][3] = -1.;
# Line 1169  if(debug)  printf("loopze inizio max  %d Line 1169  if(debug)  printf("loopze inizio max  %d
1169    
1170    //loop isotopi    //loop isotopi
1171    while ( elem[inucl-1][isotope] > 0. ){      while ( elem[inucl-1][isotope] > 0. ){  
1172    
1173              if( fiso != -1 ){
1174                      isotope=fiso;
1175                      if(debug) printf("In Loopze - Isotope N %d",isotope);
1176              }
1177      Massa = elem[inucl-1][isotope]*MassP;      Massa = elem[inucl-1][isotope]*MassP;
1178    
1179      //loop energia      //loop energia
# Line 1196  if(debug)  printf("loopze inizio max  %d Line 1201  if(debug)  printf("loopze inizio max  %d
1201                          matrixchi2[inucl][isotope][iene2][2]=chi2[2];//piani saltati nel chi2                          matrixchi2[inucl][isotope][iene2][2]=chi2[2];//piani saltati nel chi2
1202              if( fene > 0. ) break;              if( fene > 0. ) break;
1203           } else {                 } else {      
1204                          matrixchi2[inucl][isotope][iene2][0]=1000.;//valore chi2 per questo z a questa energia                          matrixchi2[inucl][isotope][iene2][0]=numeric_limits<Float_t>::max();//valore chi2 per questo z a questa energia
1205                          matrixchi2[inucl][isotope][iene2][1]=1000.;//energia per questo chi2                          matrixchi2[inucl][isotope][iene2][1]=numeric_limits<Float_t>::max();//energia per questo chi2
1206                          matrixchi2[inucl][isotope][iene2][2]=1000.;//piani saltati nel chi2                          matrixchi2[inucl][isotope][iene2][2]=numeric_limits<Float_t>::max();//piani saltati nel chi2
1207                  break;                  break;
1208        }        }
1209    
1210      }//fine loop energia      }//fine loop energia
1211    
1212      if( fiso != -1 ){
1213              if(debug) printf("exited form isotopes loop");
1214              break;
1215      }
1216    
1217    isotope++; //incremento il contatore isotopi    isotope++; //incremento il contatore isotopi
1218    }//fine loop isotopi    }//fine loop isotopi
1219    isotope=0; //riazzero il contatore isotopi    isotope=0; //riazzero il contatore isotopi
# Line 1213  if(debug)  printf("loopze inizio max  %d Line 1223  if(debug)  printf("loopze inizio max  %d
1223    isotope=0;//non dovrebbe servire    isotope=0;//non dovrebbe servire
1224    
1225    //Emi    //Emi
1226    for (Int_t nu=zmin; nu<max; nu++){   for (Int_t nu=zmin; nu<max; nu++){
1227      while(elem[nu-1][isotope]> 0.){     if( fiso != -1 ){
1228              isotope=fiso;
1229              if(debug) printf("In Loopze EMI - Isotope N %d",isotope);
1230       }
1231         while(elem[nu-1][isotope]> 0.){
1232          for (Int_t en=0; en<nostep; en++){          for (Int_t en=0; en<nostep; en++){
1233          if((matrixchi2[nu][isotope][en][0]<bestchi2[0]) && (matrixchi2[nu][isotope][en][0] >0.)){          if((matrixchi2[nu][isotope][en][0]<bestchi2[0]) && (matrixchi2[nu][isotope][en][0] >0.)){
1234            bestchi2[0]= matrixchi2[nu][isotope][en][0];// chi2            bestchi2[0]= matrixchi2[nu][isotope][en][0];// chi2
# Line 1223  if(debug)  printf("loopze inizio max  %d Line 1237  if(debug)  printf("loopze inizio max  %d
1237            bestchi2[3]= matrixchi2[nu][isotope][en][2];// totale piani saltati            bestchi2[3]= matrixchi2[nu][isotope][en][2];// totale piani saltati
1238        bestchi2[4]= (Float_t)isotope; //isotopo        bestchi2[4]= (Float_t)isotope; //isotopo
1239          }          }
1240        }       }
1241    
1242      if( fiso != -1 ){
1243        if(debug) printf("exited form isotopes loop");
1244        break;
1245      }
1246    
1247      isotope++;      isotope++;
1248      }      }
1249     isotope=0;     isotope=0;
1250    }   }
1251    
1252  };//endloopze  };//endloopze
1253    
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.23