Corrections for the TGeo interface. (A. Gheata)
[u/mrichter/AliRoot.git] / TFluka / endraw.cxx
CommitLineData
fa3d1cc7 1#include <Riostream.h>
8fd1d27e 2#include "TVirtualMCApplication.h"
df14e7fb 3#include "TGeoMaterial.h"
4#include "TGeoManager.h"
18e0cabb 5#include <TParticle.h>
df14e7fb 6#include "TFlukaCerenkov.h"
a7bb59a2 7
fa3d1cc7 8#include "TFluka.h"
d566901f 9#include "TFlukaCodes.h"
a7bb59a2 10
7dac99f1 11#include "Fdimpar.h" //(DIMPAR) fluka include
12#include "Ftrackr.h" //(TRACKR) fluka common
d566901f 13#include "Fltclcm.h" //(LTCLCM) fluka common
14#include "Fpaprop.h" //(PAPROP) fluka common
18e0cabb 15
fa3d1cc7 16#ifndef WIN32
17# define endraw endraw_
18#else
19# define endraw ENDRAW
20#endif
21extern "C" {
22void endraw(Int_t& icode, Int_t& mreg, Double_t& rull, Double_t& xsco, Double_t& ysco, Double_t& zsco)
23{
4b81a588 24 TFluka* fluka = (TFluka*) gMC;
d566901f 25 // nothing to do if particle in dummy region
26 if (mreg == fluka->GetDummyRegion()) return;
2bc4c610 27 Int_t verbosityLevel = fluka->GetVerbosityLevel();
2bc78182 28 Bool_t debug = (verbosityLevel >= 3)? kTRUE : kFALSE;
4aba9d66 29 Int_t mlttc = (icode==kKASKADinelarecoil) ? TRACKR.lt2trk : TRACKR.lt1trk; //LTCLCM.mlatm1;
d566901f 30 fluka->SetCaller(kENDRAW);
4b81a588 31 fluka->SetRull(rull);
32 fluka->SetXsco(xsco);
33 fluka->SetYsco(ysco);
34 fluka->SetZsco(zsco);
d566901f 35 fluka->SetMreg(mreg, mlttc);
4aba9d66 36
022d4275 37 Float_t edep = rull;
38
df14e7fb 39 if (TRACKR.jtrack == -1) {
18e0cabb 40 // Handle quantum efficiency the G3 way
2bc4c610 41 if (debug) printf("endraw: Cerenkov photon depositing energy: %d %e\n", mreg, rull);
df14e7fb 42 TGeoMaterial* material = (gGeoManager->GetCurrentVolume())->GetMaterial();
df14e7fb 43 TFlukaCerenkov* cerenkov = dynamic_cast<TFlukaCerenkov*> (material->GetCerenkovProperties());
44 if (cerenkov) {
18e0cabb 45 Double_t eff = (cerenkov->GetQuantumEfficiency(rull));
46 if (gRandom->Rndm() > eff) {
47 edep = 0.;
48 }
df14e7fb 49 }
50 }
2bc78182 51
18e0cabb 52 TVirtualMCStack* cppstack = fluka->GetStack();
18e0cabb 53
54 if (debug) {
4aba9d66 55 cout << "ENDRAW For icode=" << icode
18e0cabb 56 << " track=" << TRACKR.ispusr[mkbmx2-1] << " pdg=" << fluka->PDGFromId(TRACKR.jtrack)
4aba9d66 57 << " edep="<< edep << " mreg=" << mreg << endl;
18e0cabb 58 }
59
4aba9d66 60 // check region lattice consistency (debug Ernesto)
61 // *****************************************************
62 Int_t nodeId;
63 Int_t volId = fluka->CurrentVolID(nodeId);
64 Int_t crtlttc = gGeoManager->GetCurrentNodeId()+1;
4b1c3976 65 if(debug && mreg != volId && !gGeoManager->IsOutside() ) {
4aba9d66 66 cout << " endraw: track=" << TRACKR.ispusr[mkbmx2-1] << " pdg=" << fluka->PDGFromId(TRACKR.jtrack)
67 << " icode=" << icode << " gNstep=" << fluka->GetNstep() << endl
68 << " fluka mreg=" << mreg << " mlttc=" << mlttc << endl
69 << " TGeo volId=" << volId << " crtlttc=" << crtlttc << endl
70 << " common TRACKR lt1trk=" << TRACKR.lt1trk << " lt2trk=" << TRACKR.lt2trk << endl
71 << " common LTCLCM newlat=" << LTCLCM.newlat << " mlatld=" << LTCLCM.mlatld << endl
72 << " mlatm1=" << LTCLCM.mlatm1 << " mltsen=" << LTCLCM.mltsen << endl
73 << " mltsm1=" << LTCLCM.mltsm1 << " mlattc=" << LTCLCM.mlattc << endl;
74 if( mlttc == crtlttc ) cout << " *************************************************************" << endl;
75 }
76 // *****************************************************
77
78
d566901f 79 if (icode != kEMFSCOstopping1 && icode != kEMFSCOstopping2) {
80 fluka->SetIcode((FlukaProcessCode_t)icode);
022d4275 81 fluka->SetRull(edep);
18e0cabb 82 if (icode == kKASKADelarecoil && TRACKR.ispusr[mkbmx2-5]) {
4aba9d66 83 // Elastic recoil and in stuprf npprmr > 0,
84 // the secondary being loaded is actually still the interacting particle
85 cppstack->SetCurrentTrack( TRACKR.ispusr[mkbmx2-4] );
86 // cout << "endraw elastic recoil track=" << TRACKR.ispusr[mkbmx2-1] << " parent=" << TRACKR.ispusr[mkbmx2-4]
87 // << endl;
18e0cabb 88 }
89 else
4aba9d66 90 cppstack->SetCurrentTrack(TRACKR.ispusr[mkbmx2-1] );
022d4275 91 (TVirtualMCApplication::Instance())->Stepping();
92 } else {
2bc78182 93 //
d566901f 94 // For icode 21,22 the particle has fallen below thresshold.
2bc78182 95 // This has to be signalled to the StepManager()
96 //
18e0cabb 97 cppstack->SetCurrentTrack( TRACKR.ispusr[mkbmx2-1] );
a5b73d69 98 fluka->SetRull(edep);
d566901f 99 fluka->SetIcode((FlukaProcessCode_t) icode);
022d4275 100 (TVirtualMCApplication::Instance())->Stepping();
101 fluka->SetTrackIsNew(kFALSE);
d566901f 102 fluka->SetIcode((FlukaProcessCode_t)icode);
a5b73d69 103 fluka->SetRull(0.);
022d4275 104 (TVirtualMCApplication::Instance())->Stepping();
18e0cabb 105// cppstack->SetCurrentTrack( saveTrackId );
106
022d4275 107 }
fa3d1cc7 108} // end of endraw
109} // end of extern "C"
110