GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » General » Fieldmaps for the half current solenoid
Re: Fieldmaps for the half current solenoid [message #15300 is a reply to message #15294] Thu, 22 August 2013 09:47 Go to previous messageGo to previous message
Simone Esch is currently offline  Simone Esch
Messages: 37
Registered: August 2010
Location: Forschungszentrum Jülich...
continuous participant

From: *ext.kfa-juelich.de
Hello Prometeusz, Hello Pandas,

thank you for organizing the files for the half current solenoid maps.
I have just a question about there usage and implementation in my PandaRoot version. I saw that you just replaced them, but I would like to have automaticly the right files in the classe (otherwise I will screw it up for sure). And I would like to understand which one I really need.

I was told, that the Dipol- and Transitionmaps for a beammomenta of e.g. 1.6 Gev/c are an interpolation of the maps of 1.5 GeV/c and 4.06 GeV/c. So I would expect that I have to deliver maps of two momenta for the transition field for each moemtasetting of the solenoid.
You have uploaded all these files which matched with my expectation.


But if I have a look in the classes (PndTransMap and PndDipolMap), I see that the classes just load one and not two maps, what I would expect if they do an interpolation.
(the following code is out if PndTransMap)

  fType = 4;
   TString Suffix="";
   FairRunSim *fRun= FairRunSim::Instance();
   if(fRun) fBeamMom= fRun->GetBeamMom();
   
   if(fBeamMom< 3)Suffix=".0150" ;
   else if (fBeamMom< 6.0 && fBeamMom >= 3.0)Suffix=".0406";
   else if (fBeamMom< 10.0 && fBeamMom >= 6.0 )Suffix=".0890" ;
   else if (fBeamMom< 13.0 && fBeamMom >= 10.0)Suffix=".1191";
   else if (fBeamMom> 13.0) Suffix=".1500";
   
   
   TString NewName=mapName;
   NewName=mapName+Suffix;
   SetName(NewName.Data());
   TString dir = getenv("VMCWORKDIR");
   fFileName = dir + "/input/" + NewName;
   if ( fileType[0] == 'R' ) fFileName += ".root";
   else                      fFileName += ".dat";



Whith this code I would expect to just need a transition map for 1.5 GeV/c, and not 4.06 GeV/c.

So I am not sure what to do now. Which fils do I have to change for 1.6 and 1.9 GeV/c beammomenta?

Best regards

Simone

 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: [NOTSUPPORTED]Problem in Installation of external package dec 2013
Next Topic: Target Specifications
Goto Forum:
  


Current Time: Tue Dec 07 00:58:04 CET 2021

Total time taken to generate the page: 0.02435 seconds