Physics Log

Time:
 11:56 am
Date:  Friday, July 22
Projects currently working on:
   
Efficiency calculations

Issues:
   
I have now rerun all the efficiency data for the comparison of a known two particle system with a known three particle system, for both the low and high field data. However, completing the calculations for all the efficiencies continues to elude me.

Problems solved:
   
All the data is completed for the low and high field data of the K+, K-, and P efficiency calculations.

Problems unresolved:
    I have thoroughly checked my code and cannot find any problems with it. The calculations for the K+ efficiency is fine and turns out as it did before, around 90% in the region with enough data for accurate measurements (0 to 40 degrees). However, I am having serious problems with the K- and P efficiencies. It seems there are far too many K- events present. Without any cuts, the events I am getting for tlkm (the K- theta distribution) for PK+K- is 29822, whereas with tlkm for PK+X the number is 408740, which is approximately 20 times as many!!
    I don't understand how I managed to get all the calculations to work out when I originally calculated the efficiencies about a month ago. I am going to have to retrace my steps and hope I can solve this problem fast. Getting all the data run and completed by next week will not be hard, but if I continue to have problems like this that are going to take an indefinite amount of time to solve, actually producing the results from the data on time may get nerve-racking.






Time:
1:54 pm
Date: Wednesday, July 20
Projects currently working on:
    
Efficiency calculations
    Working on completing the webpage
    Analysing low field deuterium data
     Complete documentation of data files and code on my userspace

Issues:
    The efficiency calculations have run into some trouble and so I have decided to run the data over. The main problem is that the contents of the .root files within the efficiency folder do not seem to contain the correct data. This is bizarre because I was able to complete my efficiency calculations before and posted the results to my webpage, but now going back I cannot reproduce these results. Hopefully once the data has been redone, the calculations should be easy to complete.
       I am hoping to be finished with the documentation of my userspace soon which should help clarify the use of each program and what each data file contains. This will be especially useful in the future if other people need access to my account and the code I have written.

Problems solved:
   
The main part of the webpage has been set up now and all that is needed from here are modifications.
     There are a few remaining areas I have to document for my userspace and once those are done I need to post the documentation on the webpage.

Problems unresolved:
   
The efficiency data is still being processed so until that has finished running I will not be able to post the results of the efficiency calculations on the website. For the moment the old efficiency calculations are still up.







Time:
2:23 pm
Date: Tuesday, June 28
Projects currently working on:
   
High-field data analysis for the reaction: gamma d -> p k+ k- X
     More efficiency results for the reaction: gamma p -> p k+ k-

Issues:
   
Although I have still not received a single e-mail from anyone, despite it being a week now, I have at least been added to the CLAS group. However, attempting to run my program to read the CLAS  data has still proven unsuccessful. Also, the only files for the low-field hydrogen data are actually bad runs. None of the good runs that I need to access are available, which means I will have to ask Dr. Mibe about this as well, and hope he will eventually respond to e-mail.

Problems solved:
  
I have FINALLY been added to the CLAS group which means I should be able to access the files I need.

Problems unresolved:
   
The correct .B00 files are unavailable for the low-field hydrogen runs. As well, currently my code for the high-field D2_full data does not work.






Time:
3:05 pm
Date: Monday, June 27
Projects currently working on:
   
High-field data analysis for the reaction: gamma d -> p k+ k- X
    More efficiency results for the reaction: gamma p -> p k+ k-

Issues:
   
I have still received no word from Dr. Tedeschi, Dr. Mibe or Wang Xu about the status of both my CLAS group account or my code. Until I receive an e-mail I will be unable to proceed.







Time:
4:35 pm
Date: Friday, June 17
Projects currently working on:
   
Efficiency for the reaction: gamma p -> p k+ k-              COMPLETED
    High-field data analysis for the reaction: gamma d -> p k+ k- X

Issues:
   
I was finally able to resolve the issues I had with attempting to divide graphs, and I am now able to, for the most part, manipulate histograms completely, including the axes and the labels.

Problems solved:
   
I finished up the efficiency plots completely, and analyzed them.  [Report on efficiency plots]

Problems unresolved:
   
I am still currently working on the high-field data analysis.




Time: 2:17 pm
Date: Wednesday, June 15
Projects currently working on:
   
Efficiency for the reaction: gamma p -> p k+ k-
    High-field data analysis for the reaction: gamma d -> p k+ k- X

Issues:
   
I have been messing around with the new dmep server and making sure that I can get my files to run correctly on the new machine.
    Also, I am still attempting to sort out the issues regarding dividing the two histograms. Although Prof. Gao has said that an actual plot of this efficiency is not required, now that I am so close to completing it I might as well go ahead and finish it. Also, it has required me to go into the Tree structure branch by branch and edit data which should be very useful in the future, since I can then directly manipulate the histograms I am editing. Kevin Kramer has also been a help in this regard, since he has written some code on a similar topic.
    The graph I have below for PXK- as compared to PK+K- still needs to be done for K- as X and for P as X. I am in the process of getting these done, but there have been some issues with Zombie files that are stopping me from manipulating the created .root files.
    Also, I have begun work on the high-field data but have not made large amounts of progress yet, due to trying to solve the above problems.

Problems solved:
   
I have succeeded in getting the make files to work on the new dmep server, and am currently attempting to run some of the programs on it. If this works out it will be EXTREMELY helpful because I will finally be able to run multiple jobs at the same time without worrying about getting my access revoked!
   
Problems unresolved:
   
I have yet to remedy the problems with the Zombie files. It seems that the .root files I have been creating have not been closed once I am done with them. I am attempting to rerun the program on both the smeagol and dmep servers in the hope that one of them will produce .root files I can use. Once this is remedied it should take me no time at all to create all the graphs.
    I need to talk to Kevin as there are a few problems with my code for creating the two histograms. Once this is fixed,  and I fix the .root problem I have it should be very simple for me to create a complete set of efficiency graphs.
    There is also the general project of getting the high-field data analyzed.






Time:
3:15 pm
Date: Monday, June 13
Projects currently working on:
    
Efficiency for the reaction:  gamma p -> p k+ k-.

Issues:
   
I have created plots of the theta distribution (after the required proton acceptance, k- acceptance, and missing mass X cuts) for both the missing mass X in the reaction gamma p -> p X k-  and the calculated missing mass of K+ in the reaction gamma p -> p k+ k-.
Unlike before, these seem to have turned out alright (on the left is for PK+K- and on the right is for PXK-)

Missing mass distributions

    The issue is that the two plots have different scaling for the x-axis and the bin widths are different. Therefore, when I attempt to divide one graph by the other in order to obtain a plot of efficiency as dependent on the solid angle theta, the plot produced is incorrect.

Efficiency

    If however, I choose to calculate the efficiency in my head (dividing one plot by the other), the efficiency seems viable. It appears that the efficiency is around 60% which at least seems somewhat valid.
    I need to be able to change my TTree ppk into a histogram and then I should be able to rescale the x-axis for both missing mass plots as well as create the same number of bins for each plot. The other advantage of changing it to a histogram is that I should be able to edit the properties of the plots and add labels.

Problems solved:
   
I have managed to solve the problem of turning the TTree pkk into a histogram using the h->Draw("var>>A,"cuts") command where A is the histogram TH1F to be created.

Problems unresolved:
   
Despite creating the histograms, and defining the bin numbers and axis limits within them, ROOT continues to automatically assign the bins and x-axis limits, and so the plot for efficiency remains as it was, and still needs to be fixed. I have e-mailed Jason about this and hope to hear back from him soon.





Time: 3:06 pm
Date: Monday, June 13
Projects currently working on:
   
At the moment I am trying to set up this Physics log in order to keep track of my progress, and also to inform Prof. Haiyan Gao of what I am up to.

Problems solved:
   
Well I guess in one form or another I have solved the problem of this Physics log. At first I thought I was going to have to use straight HTML, which would have sucked, but instead I managed to download NVU which facilitates web-design on Linux machines.

Problems unresolved:
   
Everything to do with the efficiency!!