OK, I find the core problem is not getting the proper detector parameters
form the db. By added a Print("db") statement I find:
BrDetectorVolume: TPM1 - TPM1
Detector Volume TPM1, Type :
Size: 0.00 0.00 0.00
Position 0.00 0.00 0.00
(0.000000,0.000000,0.000000)
( 0.000000, 0.000000, 0.000000)
( 0.000000, 0.000000, 0.000000)
( 0.000000, 0.000000, 0.000000)
(0.000000,0.000000,0.000000)
( 0.000000, 0.000000, 0.000000)
( 0.000000, 0.000000, 0.000000)
( 0.000000, 0.000000, 0.000000)
Now if I can only figure out what the correct name is...
...steve
Stephen J. Sanders wrote:
> Hi,
> I've been trying unsuccessfully to get the tpm1 vertex using the
> BrTPMTrackVertexModule. This worked about a month ago, but
> now I'm not getting anything reasonable. By forcing a printout of
> some of the diagnostic information, it looks like the problem is
> with the LocalToGlobal transformation, which always returns
> zero:
>
> //sjs if(DebugLevel()>2){
> cout<< "volume name:" << fFrontVolume->GetName() << endl;
> cout<< "Vertex0: position" << FrontTrack_p->GetPos() << endl;
> cout<< "Vertex0: alpha" << FrontTrack_p->GetSlope() << endl;
> //sjs }
>
> fFrontVolume->LocalToGlobal(FrontTrack_p->GetPos(),xg,0);
> // fFrontVolume->LocalToGlobal(FrontTrack_p->GetAlpha(),ag,1);
> fFrontVolume->LocalToGlobal(FrontTrack_p->GetSlope(),ag,1);
> //sjs if(DebugLevel()>2){
> cout<< "Vertex0: xg[0]" << xg[0] << endl;
> cout<< "Vertex0: xg[1]" << xg[1] << endl;
> cout<< "Vertex0: xg[2]" << xg[2] << endl;
> cout<< "Vertex0: ag[0]" << ag[0] << endl;
> cout<< "Vertex0: ag[1]" << ag[1] << endl;
> cout<< "Vertex0: ag[2]" << ag[2] << endl;
> //sjs }
>
> results in the following output:
>
> volume name:TPM1
> Vertex0: position(-9.944769, -3.670153, 0.000000)
> Vertex0: alpha(-0.237948, -0.030677, 0.970793)
> Vertex0: xg[0]0.000000
> Vertex0: xg[1]0.000000
> Vertex0: xg[2]0.000000
> Vertex0: ag[0]0.000000
> Vertex0: ag[1]0.000000
> Vertex0: ag[2]0.000000
>
> Is anyone currently using this module for vertex determination? Any
> ideas on what might
> be the problem?
>
> ...steve
>
This archive was generated by hypermail 2b30 : Thu May 02 2002 - 22:43:07 EDT