/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/analysissubroutines.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/F77/analysissubroutines.f

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

revision 1.46 by mocchiut, Fri Jan 17 10:14:38 2014 UTC revision 1.47 by pam-ts, Wed Jun 4 07:57:04 2014 UTC
# Line 1046  c$$$      PFAy = 'COG4'!PFA Line 1046  c$$$      PFAy = 'COG4'!PFA
1046           ipy=npl(VIEW(icy))           ipy=npl(VIEW(icy))
1047                    
1048           if( (nplanes-ipx+1).ne.ip )then           if( (nplanes-ipx+1).ne.ip )then
1049              print*,'xyzpam: ***WARNING*** cluster ',icx              print*,'xyzpam: ***WARNING*** cluster icx=',icx
1050       $           ,' does not belong to plane: ',ip       $           ,' belongs to plane ',(nplanes-ipx+1)            
1051         $           ,' and not ',ip
1052              icx = -1*icx              icx = -1*icx
1053              return              return
1054           endif           endif
1055           if( (nplanes-ipy+1).ne.ip )then           if( (nplanes-ipy+1).ne.ip )then
1056              print*,'xyzpam: ***WARNING*** cluster ',icy              print*,'xyzpam: ***WARNING*** cluster icy=',icy
1057       $           ,' does not belong to plane: ',ip       $           ,' belongs to plane ',(nplanes-ipy+1)            
1058              icy = -1*icy       $           ,' and not ',ip
1059                 icy = -1*icy
1060              return              return
1061           endif           endif
1062    
# Line 1070  c$$$      PFAy = 'COG4'!PFA Line 1072  c$$$      PFAy = 'COG4'!PFA
1072           zm(ip) = zPAM           zm(ip) = zPAM
1073           xm_A(ip) = 0.D0           xm_A(ip) = 0.D0
1074           ym_A(ip) = 0.D0           ym_A(ip) = 0.D0
1075             zm_A(ip) = 0.D0
1076           xm_B(ip) = 0.D0           xm_B(ip) = 0.D0
1077           ym_B(ip) = 0.D0           ym_B(ip) = 0.D0
1078             zm_B(ip) = 0.D0
1079    
1080  c         zv(ip) = zPAM  c         zv(ip) = zPAM
1081    
# Line 1079  c         zv(ip) = zPAM Line 1083  c         zv(ip) = zPAM
1083    
1084           ipy=npl(VIEW(icy))           ipy=npl(VIEW(icy))
1085           if( (nplanes-ipy+1).ne.ip )then           if( (nplanes-ipy+1).ne.ip )then
1086              print*,'xyzpam: ***WARNING*** cluster ',icy              print*,'xyzpam: ***WARNING*** cluster icy=',icy
1087       $           ,' does not belong to plane: ',ip       $           ,' belongs to plane ',(nplanes-ipy+1)            
1088         $           ,' and not ',ip
1089              icy = -1*icy              icy = -1*icy
1090              return              return
1091           endif           endif
# Line 1100  c$$$         zm(ip) = (zPAM_A+zPAM_B)/2. Line 1105  c$$$         zm(ip) = (zPAM_A+zPAM_B)/2.
1105           zm(ip) = zPAM           zm(ip) = zPAM
1106           xm_A(ip) = xPAM_A           xm_A(ip) = xPAM_A
1107           ym_A(ip) = yPAM_A           ym_A(ip) = yPAM_A
1108             zm_A(ip) = zPAM_A
1109           xm_B(ip) = xPAM_B           xm_B(ip) = xPAM_B
1110           ym_B(ip) = yPAM_B           ym_B(ip) = yPAM_B
1111             zm_B(ip) = zPAM_B
1112    
1113  c         zv(ip) = (zPAM_A+zPAM_B)/2.  c         zv(ip) = (zPAM_A+zPAM_B)/2.
1114                    
# Line 1110  c         zv(ip) = (zPAM_A+zPAM_B)/2. Line 1117  c         zv(ip) = (zPAM_A+zPAM_B)/2.
1117           ipx=npl(VIEW(icx))           ipx=npl(VIEW(icx))
1118    
1119           if( (nplanes-ipx+1).ne.ip )then           if( (nplanes-ipx+1).ne.ip )then
1120              print*,'xyzpam: ***WARNING*** cluster ',icx              print*,'xyzpam: ***WARNING*** cluster icx=',icx
1121       $           ,' does not belong to plane: ',ip       $           ,' belongs to plane ',(nplanes-ipx+1)            
1122         $           ,' and not ',ip
1123              icx = -1*icx              icx = -1*icx
1124              return              return
1125           endif           endif
# Line 1131  c$$$         zm(ip) = (zPAM_A+zPAM_B)/2. Line 1139  c$$$         zm(ip) = (zPAM_A+zPAM_B)/2.
1139           zm(ip) = zPAM           zm(ip) = zPAM
1140           xm_A(ip) = xPAM_A           xm_A(ip) = xPAM_A
1141           ym_A(ip) = yPAM_A           ym_A(ip) = yPAM_A
1142             zm_A(ip) = zPAM_A
1143           xm_B(ip) = xPAM_B           xm_B(ip) = xPAM_B
1144           ym_B(ip) = yPAM_B           ym_B(ip) = yPAM_B
1145             zm_B(ip) = zPAM_B
1146                    
1147  c         zv(ip) = (zPAM_A+zPAM_B)/2.  c         zv(ip) = (zPAM_A+zPAM_B)/2.
1148    
# Line 1153  c         zv(ip) = (zPAM_A+zPAM_B)/2. Line 1163  c         zv(ip) = (zPAM_A+zPAM_B)/2.
1163           zm(ip) = z_mech_sensor(nplanes-ip+1,il,is)*1000./1.d4           zm(ip) = z_mech_sensor(nplanes-ip+1,il,is)*1000./1.d4
1164           xm_A(ip) = 0.           xm_A(ip) = 0.
1165           ym_A(ip) = 0.           ym_A(ip) = 0.
1166             zm_A(ip) = 0.
1167           xm_B(ip) = 0.           xm_B(ip) = 0.
1168           ym_B(ip) = 0.           ym_B(ip) = 0.
1169             zm_B(ip) = 0.
1170    
1171  c         zv(ip) = z_mech_sensor(nplanes-ip+1,il,is)*1000./1.d4  c         zv(ip) = z_mech_sensor(nplanes-ip+1,il,is)*1000./1.d4
1172    
Legend:
Removed from v.1.46  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.23