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-lReceived on Fri Sep 16 11:30:22 2005
This archive was generated by hypermail 2.1.8 : Fri Sep 16 2005 - 11:30:32 EDT
