[Brahms-dev-l] RE: analysis macro

Hi Ramiro,

    thanks for your reply. I have a question about the triggers.
What are the mappings between the various triggers and the hexadecimal numbers
For example you use for pp at 200 GeV for year 05,

G_fTrigger&0x0200 for trigger 2
G_fTrigger&0x2000 for trigger 6
G_fTrigger&0x1000 for trigger 5 (This mibias trigger, right?)

For CuCu at 200 GeV, I was told that I need to use

G_fTrigger&0x0800

for minbias trigger. Do you have any idea what I should use for
trigger 2? Is it the same as for pp, i.e,

G_fTrigger&0x0200

Thanks in advance,

Selemon,

> -----Original Message-----
> From: Ramiro Debbe [mailto:debbe@rcf.rhic.bnl.gov]
> Sent: Mon 10/3/2005 12:16 PM
> To: Bekele, Selemon
> Subject: Re: analysis macro
>
> Hi, Selemon,
>
>
> On Oct 3, 2005, at 12:47 PM, Bekele, Selemon wrote:
>
>>
>> Hi Ramiro,
>>
>>    thanks for your reply. Just to be sure I understand your answers:
>> ------------(1)
>> ###########################
>> In the expression  " 1/nEvents * 1/pt *d2N/dptdy "
>>
>> is nEvents the same number as
>>
>>      mNEvents[vertexBin] += mDstAnalysis->R_fScaledown5[0];
>>
>> used to normalize the spectra?
>>
>> #############################
>>
> Yes mNEevents[.... is the normalization I fill it this way:
>          if ((G_fTrigger[0]&0x1000)!=0) {
>            fNEvents[vertexBin] += R_fScaledown5[0]/R_fScaledown2[0];
>          }
>
>> ------------- (2)
>> ################################
>> Suppose I want to pick up trigger 2 events and look at charged 
>> hadrons.
>> Is what I do below OK?
>>
>>  if(G_fTrigger[0]&0x0200)!=0){
>>  ....
>>  }
>>
>> mOneDimPt2->Fill(pt,  mDstAnalysis->R_fScaledown2[0]);
>> mYieldRaw2[vertexBin]->Fill(eta, pt,  mDstAnalysis->R_fScaledown2[0]);
>> ################################
>>
>> or I should do instead
>>
>> ---------------(3)
>> ################################
>>
>>  if(G_fTrigger[0]&0x0200)!=0){
>>  ....
>>  }
>>
>>  mOneDimPt2->Fill(pt);
>>  mYieldRaw2[vertexBin]->Fill(eta, pt);
>> ################################
>>
> If you have both scaledowns (5 and 2) in the Events array, then you are
> right. Again nowadays I use the array
> hYieldRaw[pid][vertexbin] and pid runs from 0 to 3
>
> Ramiro
>> I know you are very busy, thanks for taking the time to answer
>> my questions.
>>
>>
>> -----Original Message-----
>> From: Ramiro Debbe [mailto:debbe@rcf.rhic.bnl.gov]
>> Sent: Mon 10/3/2005 11:19 AM
>> To: Bekele, Selemon
>> Subject: Re: analysis macro
>>
>> Hi Selemon,
>> The first two lines pickup exclusive trigger 6 and 2 respectively
>>
>> The histograms OneDimPt6 etc are left from the beginning of the
>> project, they were used for diagnostic purposes. I also used YieldRaw2
>> and 6 to look at charged hadrons. You can ignore them, I should clean
>> up my scripts.
>> For the arrays that I really use the scaledowns are carried in the
>> array that contains the trigger 5 normalization.
>> The plan was to leave the array that carry raw counts free of any
>> correction, I never got to do that because I still correct for 
>> tracking
>> and RICH efficiency per track. That makes my statistical errors some
>> 2-3% bigger.
>> That is also something I will fix in the near future, right now I'm
>> trying to explain some problems in pp positives at 4 degrees.
>>
>> Ramiro
>>
>> On Oct 3, 2005, at 11:05 AM, Bekele, Selemon wrote:
>>
>>>
>>> Hi Ramiro,
>>>
>>>      hope you are doing fine. I have a few questions regarding
>>> some of the lines in the analysis macro you use:
>>>
>>> what do the if conditions mean?
>>>
>>> if(G_fTrigger[0]&0x2000)!=0){
>>> ....
>>> }
>>> if(G_fTrigger[0]&0x0200)!=0){
>>> ....
>>> }
>>>
>>> why do you do, e.g.,
>>>
>>> mOneDimPt6->Fill(pt,  mDstAnalysis->R_fScaledown6[0]);
>>> mYieldRaw6[vertexBin]->Fill(eta, pt,  
>>> mDstAnalysis->R_fScaledown6[0]);
>>> mOneDimPt2->Fill(pt,  mDstAnalysis->R_fScaledown2[0]);
>>> mYieldRaw2[vertexBin]->Fill(eta, pt,  
>>> mDstAnalysis->R_fScaledown2[0]);
>>>
>>> instead of
>>>
>>> mOneDimPt6->Fill(pt)
>>> mYieldRaw6[vertexBin]->Fill(eta, pt)
>>> mOneDimPt2->Fill(pt);
>>> mYieldRaw2[vertexBin]->Fill(eta, pt);
>>>
>>> Thanks for the help
>>>
>>> Selemon
>>>
>>
>>
>
>


_______________________________________________
Brahms-dev-l mailing list
Brahms-dev-l@lists.bnl.gov
http://lists.bnl.gov/mailman/listinfo/brahms-dev-l