2 #include "TVirtualMCApplication.h"
3 #include "TGeoMaterial.h"
4 #include "TGeoManager.h"
6 #include "TFlukaCerenkov.h"
9 #include "TFlukaCodes.h"
11 #include "Fdimpar.h" //(DIMPAR) fluka include
12 #include "Ftrackr.h" //(TRACKR) fluka common
13 #include "Fltclcm.h" //(LTCLCM) fluka common
14 #include "Fpaprop.h" //(PAPROP) fluka common
17 # define endraw endraw_
19 # define endraw ENDRAW
22 void endraw(Int_t& icode, Int_t& mreg, Double_t& rull, Double_t& xsco, Double_t& ysco, Double_t& zsco)
24 TFluka* fluka = (TFluka*) gMC;
25 // nothing to do if particle in dummy region
26 if (mreg == fluka->GetDummyRegion()) return;
27 Int_t verbosityLevel = fluka->GetVerbosityLevel();
28 Bool_t debug = (verbosityLevel >= 3)? kTRUE : kFALSE;
29 Int_t mlttc = (icode==kKASKADinelarecoil) ? TRACKR.lt2trk : TRACKR.lt1trk; //LTCLCM.mlatm1;
30 fluka->SetCaller(kENDRAW);
35 fluka->SetMreg(mreg, mlttc);
38 Int_t ipt = fluka->PDGFromId(TRACKR.jtrack);
40 if (debug) printf("Unknown particle %5d %5d \n", TRACKR.jtrack, icode);
45 if (TRACKR.jtrack == -1) {
46 // Handle quantum efficiency the G3 way
47 if (debug) printf("endraw: Cerenkov photon depositing energy: %d %e\n", mreg, rull);
48 TGeoMaterial* material = (gGeoManager->GetCurrentVolume())->GetMaterial();
49 TFlukaCerenkov* cerenkov = dynamic_cast<TFlukaCerenkov*> (material->GetCerenkovProperties());
51 Double_t eff = (cerenkov->GetQuantumEfficiency(rull));
52 if (gRandom->Rndm() > eff) {
58 TVirtualMCStack* cppstack = fluka->GetStack();
61 cout << "ENDRAW For icode=" << icode
62 << " track=" << TRACKR.ispusr[mkbmx2-1] << " pdg=" << fluka->PDGFromId(TRACKR.jtrack)
63 << " edep="<< edep << " mreg=" << mreg << endl;
66 // check region lattice consistency (debug Ernesto)
67 // *****************************************************
69 Int_t volId = fluka->CurrentVolID(nodeId);
70 Int_t crtlttc = gGeoManager->GetCurrentNodeId()+1;
71 if(debug && mreg != volId && !gGeoManager->IsOutside() ) {
72 cout << " endraw: track=" << TRACKR.ispusr[mkbmx2-1] << " pdg=" << fluka->PDGFromId(TRACKR.jtrack)
73 << " icode=" << icode << " gNstep=" << fluka->GetNstep() << endl
74 << " fluka mreg=" << mreg << " mlttc=" << mlttc << endl
75 << " TGeo volId=" << volId << " crtlttc=" << crtlttc << endl
76 << " common TRACKR lt1trk=" << TRACKR.lt1trk << " lt2trk=" << TRACKR.lt2trk << endl
77 << " common LTCLCM newlat=" << LTCLCM.newlat << " mlatld=" << LTCLCM.mlatld << endl
78 << " mlatm1=" << LTCLCM.mlatm1 << " mltsen=" << LTCLCM.mltsen << endl
79 << " mltsm1=" << LTCLCM.mltsm1 << " mlattc=" << LTCLCM.mlattc << endl;
80 if( mlttc == crtlttc ) cout << " *************************************************************" << endl;
82 // *****************************************************
83 if (icode != kEMFSCOstopping1 && icode != kEMFSCOstopping2) {
84 fluka->SetIcode((FlukaProcessCode_t)icode);
86 if (icode == kKASKADelarecoil && TRACKR.ispusr[mkbmx2-5]) {
87 // Elastic recoil and in stuprf npprmr > 0,
88 // the secondary being loaded is actually still the interacting particle
89 cppstack->SetCurrentTrack( TRACKR.ispusr[mkbmx2-4] );
91 cppstack->SetCurrentTrack(TRACKR.ispusr[mkbmx2-1] );
94 (TVirtualMCApplication::Instance())->Stepping();
98 // For icode 21,22 the particle has fallen below thresshold.
99 // This has to be signalled to the StepManager()
101 cppstack->SetCurrentTrack( TRACKR.ispusr[mkbmx2-1] );
102 fluka->SetRull(edep);
103 fluka->SetIcode((FlukaProcessCode_t) icode);
104 (TVirtualMCApplication::Instance())->Stepping();
105 fluka->SetTrackIsNew(kFALSE);
106 fluka->SetIcode((FlukaProcessCode_t)icode);
108 (TVirtualMCApplication::Instance())->Stepping();
111 } // end of extern "C"