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

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

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

revision 1.10 by pam-fi, Mon May 14 11:03:06 2007 UTC revision 1.13 by pam-fi, Thu May 24 13:29:09 2007 UTC
# Line 164  c      print*,pfastrips Line 164  c      print*,pfastrips
164        end        end
165    
166  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***
167        real function ris_eta(ic,angle)  c      real function riseta(ic,angle)
168          real function riseta(iview,angle)
169  *--------------------------------------------------------------  *--------------------------------------------------------------
170  *     this function returns the average spatial resolution  *     this function returns the average spatial resolution
171  *     (in cm) for the ETA algorithm (function pfaeta(ic,angle))  *     (in cm) for the ETA algorithm (function pfaeta(ic,angle))
172  *     it calls:  *     it calls:
173  *     - risx_eta2(angle)  *     - risxeta2(angle)
174  *     - risy_eta2(angle)  *     - risyeta2(angle)
175  *     - risx_eta3(angle)  *     - risxeta3(angle)
176  *     - risx_eta4(angle)  *     - risxeta4(angle)
177  *     according to the angle  *     according to the angle
178  *--------------------------------------------------------------  *--------------------------------------------------------------
179        include 'commontracker.f'        include 'commontracker.f'
180        include 'level1.f'        include 'level1.f'
181        include 'calib.f'        include 'calib.f'
182    
183        ris_eta = 0        riseta = 0
184    
185        if(mod(int(VIEW(ic)),2).eq.1)then !Y-view  c      if(mod(int(VIEW(ic)),2).eq.1)then !Y-view
186          if(mod(iview,2).eq.1)then !Y-view
187                
188    
189           if( abs(angle).ge.e2fay.and.abs(angle).le.e2tay )then           if( abs(angle).ge.e2fay.and.abs(angle).le.e2tay )then
190              ris_eta = risy_eta2(angle)              riseta = risyeta2(angle)
191           elseif( abs(angle).ge.e3fay.and.abs(angle).le.e3tay )then           elseif( abs(angle).ge.e3fay.and.abs(angle).le.e3tay )then
192              ris_eta = risy_cog(angle) !ATTENZIONE!!              riseta = risy_cog(angle) !ATTENZIONE!!
193           elseif( abs(angle).ge.e4fay.and.abs(angle).le.e4tay )then           elseif( abs(angle).ge.e4fay.and.abs(angle).le.e4tay )then
194              ris_eta = risy_cog(angle) !ATTENZIONE!!              riseta = risy_cog(angle) !ATTENZIONE!!
195           else           else
196              ris_eta = risy_cog(angle)              riseta = risy_cog(angle)
197           endif                       endif            
198    
199        else                      !X-view        else                      !X-view
200    
201           if( abs(angle).ge.e2fax.and.abs(angle).le.e2tax )then           if( abs(angle).ge.e2fax.and.abs(angle).le.e2tax )then
202              ris_eta = risx_eta2(angle)              riseta = risxeta2(angle)
203           elseif( abs(angle).ge.e3fax.and.abs(angle).le.e3tax )then           elseif( abs(angle).ge.e3fax.and.abs(angle).le.e3tax )then
204              ris_eta = risx_eta3(angle)              riseta = risxeta3(angle)
205           elseif( abs(angle).ge.e4fax.and.abs(angle).le.e4tax )then           elseif( abs(angle).ge.e4fax.and.abs(angle).le.e4tax )then
206              ris_eta = risx_eta4(angle)              riseta = risxeta4(angle)
207           else           else
208              ris_eta = risx_cog(angle)              riseta = risx_cog(angle)
209           endif                       endif            
210                            
211        endif        endif
212    
213          print*,'---- ',riseta,iview,angle
214    
215   100  return   100  return
216        end        end
217    
# Line 781  c         print*,'## ',sl2,sl1,sc,sr1,sr Line 785  c         print*,'## ',sl2,sl1,sc,sr1,sr
785                 if((sl2+sl1+sc+sr1).ne.0)                 if((sl2+sl1+sc+sr1).ne.0)
786       $              COG = (sr1-sl1-2*sl2)/(sl2+sl1+sc+sr1)       $              COG = (sr1-sl1-2*sl2)/(sl2+sl1+sc+sr1)
787              elseif(sl2.le.sr2)then              elseif(sl2.le.sr2)then
788                  if((sl2+sl1+sc+sr1).ne.0)                  if((sr2+sl1+sc+sr1).ne.0)
789       $              COG = (2*sr2+sr1-sl1)/(sl2+sl1+sc+sr1)       $              COG = (2*sr2+sr1-sl1)/(sr2+sl1+sc+sr1)
790              endif              endif
791           else           else
792              print*,'function COG(NCOG,IC) ==> WARNING!! NCOG=',NCOG              print*,'function COG(NCOG,IC) ==> WARNING!! NCOG=',NCOG
# Line 814  c         print*,'-------' Line 818  c         print*,'-------'
818                 COG = COG + ipos*CLSIGNAL(i)                 COG = COG + ipos*CLSIGNAL(i)
819                 SGN = SGN + CLSIGNAL(i)                 SGN = SGN + CLSIGNAL(i)
820                 mu = mu + 1                 mu = mu + 1
821                 print*,ipos,CLSIGNAL(i)  c               print*,ipos,CLSIGNAL(i)
822              else              else
823                 goto 10                 goto 10
824              endif              endif
# Line 827  c         print*,'-------' Line 831  c         print*,'-------'
831                 COG = COG + ipos*CLSIGNAL(i)                 COG = COG + ipos*CLSIGNAL(i)
832                 SGN = SGN + CLSIGNAL(i)                 SGN = SGN + CLSIGNAL(i)
833                 mu = mu + 1                 mu = mu + 1
834                 print*,ipos,CLSIGNAL(i)  c               print*,ipos,CLSIGNAL(i)
835              else              else
836                 goto 20                 goto 20
837              endif              endif
838           enddo           enddo
839   20      continue   20      continue
840           if(SGN.le.0)then           if(SGN.le.0)then
841  c            print*,'cog(0,ic) --> ic, dedx ',ic,SGN              print*,'cog(0,ic) --> ic, dedx ',ic,SGN
842              print*,(CLSIGNAL(i)/CLSIGMA(i),i=istart,istop)              print*,(CLSIGNAL(i)/CLSIGMA(i),i=istart,istop)
843              print*,(CLSIGNAL(i),i=istart,istop)              print*,(CLSIGNAL(i),i=istart,istop)
844  c            print*,'cog(0,ic) --> NOT EVALUATED '  c            print*,'cog(0,ic) --> NOT EVALUATED '
# Line 1189  c         if(BAD(VIEW(ic),nvk(MAXS(ic)), Line 1193  c         if(BAD(VIEW(ic),nvk(MAXS(ic)),
1193    
1194  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***
1195    
1196        FUNCTION risx_eta2(x)        FUNCTION risxeta2(x)
1197    
1198        DOUBLE PRECISION V( 1)        DOUBLE PRECISION V( 1)
1199        INTEGER NPAR, NDIM, IMQFUN, I, J        INTEGER NPAR, NDIM, IMQFUN, I, J
# Line 1276  c         if(BAD(VIEW(ic),nvk(MAXS(ic)), Line 1280  c         if(BAD(VIEW(ic),nvk(MAXS(ic)),
1280     20 CONTINUE     20 CONTINUE
1281        IF (IMQFUN .EQ. 2) HQUADF = VCONST * EXP (HQUADF)        IF (IMQFUN .EQ. 2) HQUADF = VCONST * EXP (HQUADF)
1282                
1283        risx_eta2=HQUADF* 1e-4        risxeta2=HQUADF* 1e-4
1284    
1285        END        END
1286    
1287  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***
1288        FUNCTION risx_eta3(x)        FUNCTION risxeta3(x)
1289        DOUBLE PRECISION V( 1)        DOUBLE PRECISION V( 1)
1290        INTEGER NPAR, NDIM, IMQFUN, I, J        INTEGER NPAR, NDIM, IMQFUN, I, J
1291        DOUBLE PRECISION HQDJ, VV, VCONST        DOUBLE PRECISION HQDJ, VV, VCONST
# Line 1367  c         if(BAD(VIEW(ic),nvk(MAXS(ic)), Line 1371  c         if(BAD(VIEW(ic),nvk(MAXS(ic)),
1371     20 CONTINUE     20 CONTINUE
1372        IF (IMQFUN .EQ. 2) HQUADF = VCONST * EXP (HQUADF)        IF (IMQFUN .EQ. 2) HQUADF = VCONST * EXP (HQUADF)
1373                
1374        risx_eta3 = HQUADF* 1e-4        risxeta3 = HQUADF* 1e-4
1375    
1376        END        END
1377  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***
1378        FUNCTION risx_eta4(x)        FUNCTION risxeta4(x)
1379        DOUBLE PRECISION V( 1)        DOUBLE PRECISION V( 1)
1380        INTEGER NPAR, NDIM, IMQFUN, I, J        INTEGER NPAR, NDIM, IMQFUN, I, J
1381        DOUBLE PRECISION HQDJ, VV, VCONST        DOUBLE PRECISION HQDJ, VV, VCONST
# Line 1457  c         if(BAD(VIEW(ic),nvk(MAXS(ic)), Line 1461  c         if(BAD(VIEW(ic),nvk(MAXS(ic)),
1461     20 CONTINUE     20 CONTINUE
1462        IF (IMQFUN .EQ. 2) HQUADF = VCONST * EXP (HQUADF)        IF (IMQFUN .EQ. 2) HQUADF = VCONST * EXP (HQUADF)
1463                
1464        risx_eta4=HQUADF* 1e-4        risxeta4=HQUADF* 1e-4
1465    
1466        END        END
1467  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***
1468        FUNCTION risy_eta2(x)        FUNCTION risyeta2(x)
1469        DOUBLE PRECISION V( 1)        DOUBLE PRECISION V( 1)
1470        INTEGER NPAR, NDIM, IMQFUN, I, J        INTEGER NPAR, NDIM, IMQFUN, I, J
1471        DOUBLE PRECISION HQDJ, VV, VCONST        DOUBLE PRECISION HQDJ, VV, VCONST
# Line 1529  c         if(BAD(VIEW(ic),nvk(MAXS(ic)), Line 1533  c         if(BAD(VIEW(ic),nvk(MAXS(ic)),
1533     20 CONTINUE     20 CONTINUE
1534        IF (IMQFUN .EQ. 2) HQUADF = VCONST * EXP (HQUADF)        IF (IMQFUN .EQ. 2) HQUADF = VCONST * EXP (HQUADF)
1535    
1536        risy_eta2=HQUADF* 1e-4        risyeta2=HQUADF* 1e-4
1537    
1538        END        END
1539  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.23