Ignore:
Timestamp:
07/18/01 09:37:32 (23 years ago)
Author:
bigongia
Message:
Some debugging printout added. They are activated using option verbose_level 4
C.Bigongiari
:wq
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Simulation/Detector/ReflectorII/reflector.c

    r725 r870  
    3333static long writep = 0L;        /*  write ptr (on rflfile)      */
    3434
    35 
    3635extern CerHeader *cheadp;       /*  var inited in header.c      */
    3736extern RflHeader *rheadp;       /*  var inited in header.c      */
     
    4342static int ProcessEvent(CerHeader *cheadp, FILE *cerfile, FILE *rflfile);
    4443
     44FILE *chkf = NULL;     /**********************/
     45long myloop;     /**********************/
     46
    4547void main(void)
    4648{   long event = 0L;            /*  event counter       */
     
    4850    /*  Read init & geometry parms, init files and vars */
    4951    init(NULL);
     52    chkf = fopen("check", "w");     /**********************/
    5053
    5154    /*  Processing loop */
     
    6164    Log(RFLF_CLOSE_LOG);
    6265    fclose(rflfile);
     66    fclose(chkf);     /**********************/
    6367
    6468    /*  Clean memory and exit  */
     
    7882    extern void makeOmegaI(float theta, float phi);
    7983
     84
    8085    /*  Various counters: phs = absphs + refphs[0..3] + cphs  */
    8186    long phs,           /*  Number of incoming photons  */
     
    134139                    1.f - Photons[ph].u*Photons[ph].u
    135140                        - Photons[ph].v*Photons[ph].v));
     141
    136142            /*  Check absorption  */
    137 
    138143            if (absorption(wlen, Photons[ph].h, theta))
    139144                absphs++;
     
    141146            /*  Check reflection  */
    142147            else if (0 != (ref_type =
    143                            ph2cph(&Photons[ph], &CPhotons[cphs]))) 
     148                           ph2cph(&Photons[ph], &CPhotons[cphs])))
    144149                refphs[ref_type-1]++;
    145 
    146150            else    /*  Photon passed   */
    147151            {
     152            Debug("Ph %d\t%f\t%f\t%f\t%f\t%f\n",ph,
     153                   Photons[ph].x,Photons[ph].y,
     154                   Photons[ph].u,Photons[ph].v,theta);
     155            Debug("CPh %d\t%d\t%f\t%f\n\n",cphs,ph,CPhotons[cphs].x,
     156                   CPhotons[cphs].y);
     157
     158
    148159                /*  Update first/last arrival times  */
    149160                if (first > CPhotons[cphs].t) first = CPhotons[cphs].t;
     
    166177
    167178                }   /*  if overflow  */
    168             }  /*  else (=Photon passed)  */
     179            }   /*  else (=Photon passed)  */
    169180        }   /*  end of loop inside datablock    */
    170181    }   /*  end of loop on datablocks           */
     
    191202        rheadp->CPhotons     = (long) overflow * NR_OF_CPHOTONS + cphs;
    192203        fwrite(rheadp, sizeof(RflHeader), 1, rflfile);
     204        for (myloop=0; myloop<sizeof(RflHeader)/4; myloop++)
     205            fprintf(chkf, "%e ", *((float *)rheadp+myloop));
     206        fputc('\n', chkf);
    193207
    194208        /*  If there was an overflow, append data from tempfile   */
     
    252266            /*  Ok: set found at TRUE and exit loop  */
    253267            if (strncmp(cheadp->EVTH, "EVTH", 4) == 0
    254              && first_Event <= cheadp->EvtNumber
    255              &&                cheadp->EvtNumber <= last_Event
     268             && first_Event <= (long)cheadp->EvtNumber
     269             &&                (long)cheadp->EvtNumber <= last_Event
    256270             &&    low_Ecut <= cheadp->Etotal
    257271             &&                cheadp->Etotal <= high_Ecut)
Note: See TracChangeset for help on using the changeset viewer.