That's what I was doing actually, but was concerned about things like the accuracy of the calculated cross sections from using millicharge coupling, etc. tailored for the A'. I'll eventually need to adjust the model anyway for other states, since they have different diagrams, but the actual yields are especially important to get right at this point.


On further investigation, it looks like the equation I was wondering about is the (differential) cross section calculation where the Weizsaecker-Williams approximation is applied. I think this would be ok to use if I set epsilon=1 and change any numerical factors already calculated in it to reflect the correct mass in the Mandelstam variables etc., but I'll make sure nothing else is affected from the theory as well, since I'll basically be treating TM like a virtual photon.


Thanks

Brad



From: [log in to unmask] <[log in to unmask]> on behalf of Mathew Graham <[log in to unmask]>
Sent: Monday, February 16, 2015 1:46 PM
To: Bradley T Yale; Holly Vance
Cc: [log in to unmask]
Subject: Re: No clusters (and Madgraph question)
 

I don't know the details of the couplings routine, but the simplest thing to do is to generate prompt "A'" events and, in the conversion from lhe to stdhep, add the vertex displacements according to whatever model you like.  

The daughter energies appear above the threshold (strongly forward scattered), but I will explore this further. Apart from the clustering issue, it seems as though my signal is decaying too promptly as well (at the target), but I think I know why. This is a first pass at trying to create true muonium events, specifically the 1S state with the Bremsstrahlung mechanism similar to the A', and I'm currently trying to learn how to implement the correct model in its production.


This is a good time to ask someone who is familiar with MadGraph:

According to proc_card.dat within the mock.xml procedure, the model used to generate the A' is within /u/group/hps/production/mc/MadGraph/ap/MG_mini/Models/fixedAP (correct me if I'm wrong). In the file 'couplings.f' within that directory, under the section labeled 'UserMode couplings', I don't immediately recognize the 'fullcoupling' equation. Could somebody point me to where it is in literature so that I could learn more about what the procedure is doing with it? I would think that this section would need to be modified for a different signal, beyond setting the mass differently and 'epsilon=1' or something. Also please correct me if I'm way off...


Brad



From: Holly Vance <[log in to unmask]>
Sent: Friday, February 13, 2015 2:53 PM
To: Bradley T Yale
Cc: [log in to unmask]
Subject: Re: No clusters
 
Hi Bradley,

It seems that you are using the most updated clusterer if this is what you are using. 

I am not sure what type of signal you have generated. Have you checked the thresholds you are using in clustering? Is it possible that they are too high? There's a timing cut (which you may want to set to false, initially), and you can try setting the crystal, seed, and cluster thresholds to something very low. 

If you are still having trouble, you might try looking at the generated file in JAS3 event display to see what type of energies and positions your hits are. 

-Holly

On Fri, Feb 13, 2015 at 1:51 PM, Bradley T Yale <[log in to unmask]> wrote:

Hey everyone,


When doing recon using HPS2014OfflineRecon (both with and without pileup), I'm getting ecal hits but no clusters.

Could I get a stdhep file that should definitely produce clusters under those preset cuts? I was testing my own signal files, and need to make sure that something didn't happen on the generation end.

Also, if a different steering file for recon is now primarily used for this purpose, or an updated clusterer beyond ReconClusterer, that would be good to know as well.


Thanks

Brad


NOTE: This message was trained as non-spam. If this is wrong, please correct the training as soon as possible.
Spam
Not spam
Forget previous vote


Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1




Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1



Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1