Hi, Let me give a comment. sigmaang, sigmamult, sigmat are used to parameterise the overall sigmam^2, see e.g. my Ph.d. appendix D (and probably Djamels thesis as well). This also describes how to extract the values by fitting the data.... Cheers, Peter On Fri, 16 Sep 2005, Bekele, Selemon wrote: > Hi Hans, > > I found the following piece of code in > > fstree2datamap.cxx,v > > which I assume is the official version of the code in CVS from which the excutable > > "fstree2datamap" is made. > > //______________________________________________________________________ > void SetupTof2Pid(TString sett) > { > Float_t pi = mass[0]; > Float_t p = mass[2]; > > Float_t sigmaAng = 0; > Float_t sigmaMult = 0; > Float_t sigmaTof = 0; > > if (sett == "4B608") { > pi = 0.0073251; p = 0.854308; > sigmaAng = 0.0001; sigmaMult = 0.0223507; sigmaTof = 0.0015291; > } > else if (sett == "4A608") { > pi = 0.00868369; p = 0.858988; > sigmaAng = 0.0001; sigmaMult = 0.0319085; sigmaTof = 0.00151653; > } > else if (sett == "6B861") { > pi = 0.00694888; p = 0.849384; > sigmaAng = 0.0001; sigmaMult = 0.0191514; sigmaTof = 0.00156334; > } > else if (sett == "6A861") { > pi = 0.00694888; p = 0.849384; > sigmaAng = 0.0001; sigmaMult = 0.0191514; sigmaTof = 0.00156334; > } > else if (sett == "8B861") { > pi = 0.00467955; p = 0.849995; > sigmaAng = 0.00818456; sigmaMult = 0.000100182; sigmaTof = 0.00153155; > } > else if (sett == "8A861") { > pi = 0.00467955; p = 0.849995; > sigmaAng = 0.00818456; sigmaMult = 0.000100182; sigmaTof = 0.00153155; > } > > > fPidPi2 ->SetParameters(sigmaAng, sigmaMult, sigmaTof, pi); > fPidP2 ->SetParameters(sigmaAng, sigmaMult, sigmaTof, p); > > // diagnostic > fPiLow2 ->SetParameters(sigmaAng, sigmaMult, sigmaTof, pi, -kMaxSigm2); > fPiHigh2->SetParameters(sigmaAng, sigmaMult, sigmaTof, pi, kMaxSigm2); > > fPLow2 ->SetParameters(sigmaAng, sigmaMult, sigmaTof, p, -kMaxSigm2); > fPHigh2->SetParameters(sigmaAng, sigmaMult, sigmaTof, p, kMaxSigm2); > > } > > > ################################################ > > > My question is, how are the hardwired numbers for > > pi > p > sigmaAng > sigmaMult > sigmaTof > > are calculated? How does one get these values for other field settings? > > Thanks for the help, > > Selemon > > > -----Original Message----- > From: Hans Hjersing Dalsgaard [mailto:canute@fys.ku.dk] > Sent: Thu 9/15/2005 3:41 PM > To: Bekele, Selemon > Subject: RE: [Brahms-dev-l] fstree2datamap > > Hi Selemon, > > The problem is that when Truls revised the code, he focused on making it > compile. The TOF PID in FS is handled by three functions (and a fourth > that sets up hardcoded values for H2) , 'TofPid', 'TofPid2' and 'RichPid' > (I think). The RICH pid is not implemented and the H1 pid is not set up so > the only thing 'working' at the moment is the H2 pid, handled by > 'Tof2Pid'. There is, however, a problem with this function, namely that it > contains a cut in the delta(dy) which is quite reasonable but not working, > so by removing this cut you will get not-empty histograms. At the moment I > am trying to write some TOF pid in the FS that makes sense but progress is > slow. If you can help me I would be grateful. I have no idea whatsoever > about the RICH pid, as it is apparently not generated by 'dst2tree'... > Regards, > > Hans > > > Hi Casper, Hans > > > > thanks for your prompt reply. Your suggestion worked. Now I am able to > > make > > a microDst for the fs. Please have a look at the plot > > > > http://www4.rcf.bnl.gov/~tigist/BbVsZdcVertex.gif > > > > The above plot tells me the microDst is OK. Now when I try to make an > > "fstree2datamap", > > all the histograms are empty. I tried all I can to understand what is > > going on but have > > not succeded yet. This is what I do: > > > > ~/>fstree2datamap -a FS_accMap -c 0 -C 90 -b fstree -o fsdataH1_ -p H1 > > -S 4A2442 > > > > which is all one has to do to get a data map. Do you guys have any idea > > why > > the histograms turn out to be empty? > > > > Thanks, > > > > -----Original Message----- > > From: Casper Nygaard [mailto:caz@fys.ku.dk] > > Sent: Thu 9/15/2005 2:11 AM > > To: Bekele, Selemon > > Subject: Re: [Brahms-dev-l] fstree2datamap > > > > The fstree2datamap doesnt search for the .root part of the name. You > > can easily circumvent the problem by making a symbolic link to the > > file; e.g. ln -s fstree4B442.root fstree4B442 > > > > This should help. > > > > caz > > > > "Bekele, Selemon" <bekeleku@ku.edu> said: > > > >> > >> > >> Hi, > >> > >> I was trying to make data maps using the excutable > >> "fstree2datamap". > >> I have files the files > >> > >> accMap4B2442.root > >> fstree4B2442.root > >> > >> in the same directory from which I run "fstree2datamap" as shown > >> below: > >> > >> #################################3 > >> > >> tigist@rcas0032:~/MyBrahms/brahms_app/tml_app/analyze_app/anatools> > > fstree2datamap -b fstree -S 4B2442 -p H1 -c 0 -C 20 -a FS_accMap -o > > fsdataH1_ > >> > >> AcceptanceMaps::LoadMaps : Loading maps from FS_accMap4B2442.root > >> Defining 11 vertex bin limits: -50 -40 -30 -20 -10 0 10 20 30 40 50 > >> ...loading maps for h, pi, k, p, d, . > >> file : FS_accMap4B2442.root contains geometry for run 13147. > >> ----------------------------------------------------- > >> Selecting pi,K,p, setting 4 B 2442 > >> ----------------------------------------------------- > >> --> found [1] Tree files > >> Error in <TFile::TFile>: file fstree4B2442 does not exist > >> Error in <TFile::TFile>: file fstree4B2442 does not exist > >> Error in <TFile::TFile>: file fstree4B2442 does not exist > >> Error in <TFile::TFile>: file fstree4B2442 does not exist > >> --- Vertex binning from map: > >> Nbins : 10 > >> Range : -50 | 50 > >> Looping over 0 events... > >> pion efficiency map: bin [140] > >> kaon efficiency map: bin [140] > >> proton efficiency map: bin [140] > >> > >> ################################################## > >> > >> I expect the file "fstree4B2442.root" to be picked up by the > >> program. > > I wanted to look at > >> the source code to see what I might need to modify in my input but > >> could > > not find it easily. > >> Any ideas? > >> > >> Thanks, > >> > >> Selemon > >> > >> > >> _______________________________________________ > >> Brahms-dev-l mailing list > >> Brahms-dev-l@lists.bnl.gov > >> http://lists.bnl.gov/mailman/listinfo/brahms-dev-l > >> > > > > > > > > -- > > *************************************************************************** > > --------------------------------- > > / Casper Nygaard, stud.scient / "In the end it is not the amount > > of / Frederiksborgvej 54, 3.tv / years in your life, but the > > amount / 2400 Copenhagen NV / of life in your years > > that counts" / tlf#: (+45)27125518) / - > > Abraham Lincoln --------------------------------- > > *************************************************************************** > > > > > > > > > > _______________________________________________ > > Brahms-dev-l mailing list > > Brahms-dev-l@lists.bnl.gov > > http://lists.bnl.gov/mailman/listinfo/brahms-dev-l > > > > > > _______________________________________________ > Brahms-dev-l mailing list > Brahms-dev-l@lists.bnl.gov > http://lists.bnl.gov/mailman/listinfo/brahms-dev-l > -- :-) --------------------------- )-: Peter H L Christiansen pchristi@nbi.dk / (+41)764870425 :-D --------------------------- \-: _______________________________________________ Brahms-dev-l mailing list Brahms-dev-l@lists.bnl.gov http://lists.bnl.gov/mailman/listinfo/brahms-dev-lReceived on Fri Sep 16 12:11:38 2005
This archive was generated by hypermail 2.1.8 : Fri Sep 16 2005 - 12:11:43 EDT
