photrans::SimpleTransport Class Reference

#include <SimpleTransport.h>

Inheritance diagram for photrans::SimpleTransport:

List of all members.

Public Member Functions
	SimpleTransport (const char *version)
void	Update (const cfg::Config &c)
virtual	~SimpleTransport ()
jobc::Result	Reco (edm::EventHandle &evt)
Private Attributes
float	fMeVToNPhot
float	fCollectionEff
	conversion factor MeV->NPhoton
float	fAttenLength
	normalization factor
geo::Geometry *	fGeo
	the attenuation length in cm
int	msglevel

---

Constructor & Destructor Documentation

SimpleTransport::SimpleTransport (  const char *  version  )

Definition at line 31 of file SimpleTransport.cxx. References jobc::Module::SetCfgVersion(). : PhotonTransporter(), jobc::Module("SimpleTransport"), fMeVToNPhot(0), fAttenLength(0), fCollectionEff(0), fGeo(0) { this->SetCfgVersion(version); }

SimpleTransport::~SimpleTransport (   )  [virtual]

Definition at line 43 of file SimpleTransport.cxx. { if (fGeo) delete fGeo; }

---

Member Function Documentation

jobc::Result SimpleTransport::Reco (  edm::EventHandle &  evt  )  [virtual]

Reimplemented from jobc::Module. Definition at line 61 of file SimpleTransport.cxx. References assert_jobc, edm::EventHandle::DAQ(), edm::EventHandle::DetSim(), fAttenLength, fCollectionEff, fGeo, fMeVToNPhot, edm::EventHandle::Header(), geo::Geometry::Instance(), msglevel, photrans::PhotonSignal::NPhoton(), edm::EventHandle::Raw(), photrans::PhotonTransporter::SeedOnEvent(), photrans::PhotonSignal::SetPlaneCell(), and photrans::PhotonSignal::SetTrackId(). { if (!fGeo){ std::vector<const rawdata::DAQHeader*> header; short int det = rawdata::kFar; try{ evt.DAQ().Get("./",header); } catch(edm::Exception e){ std::cerr << "Error retrieving daq header, while looking "<< "in CalHit::Raw(), using default fardet geometry" << std::endl; } if(header.size() > 0) det = header[0]->DetId(); fGeo = &geo::Geometry::Instance(evt.Header().Run(), det); } //Setup the random numbers SeedOnEvent(evt); vector<const sim::FLSHitList*> hitlist(0); vector<photrans::PhotonSignal> photlist(0); assert_jobc(evt.DetSim().Get("./",hitlist), "No FLSHitLists found in DetSim() folder!"); int nhits = 0; int planeId = -1; int cellId = -1; float fGeVToMeV = 1000.; //geo::CellUniqueId lastuid; for (unsigned int i=0; i<hitlist.size(); ++i) { // geo::CellUniqueId uid; float edep = 0.; float time = 0.; float dist=0.; unsigned short trkid = 0; for (unsigned int j=0; j< hitlist[i]->fHits.size(); ++j) { ++nhits; // uid = hitlist[i]->fHits[j].fId; //if(!(lastuid==uid)) if(planeId!=hitlist[i]->fHits[j].fPlaneId || cellId!=hitlist[i]->fHits[j].fCellId) try{ //const geo::CellGeo cell = fGeo->IdToCell(uid, &planeId, &cellId); planeId=hitlist[i]->fHits[j].fPlaneId; cellId=hitlist[i]->fHits[j].fCellId; //cout <<"plane = "<<planeId<<endl; //cout << "("<<planeId<<","<<cellId<<") (" << hitlist[i]->fHits[j].fPlaneId<<","<<hitlist[i]->fHits[j].fCellId<<")"<<endl; // if (&cell == 0) continue; // cout << "cell = " << cellId << endl; // cellHalfLength = cell.HalfL(); // cout << "cell Half length = " << cellHalfLength << endl; // lastuid=uid; }catch(...){printf("bad photon geometry\n");continue;} edep = hitlist[i]->fHits[j].fEdep; time = hitlist[i]->fHits[j].fT; // double pos1[3] = {0.,0.,hitlist[i]->fHits[j].fEntryZ()}; // double pos2[3]; // cell.WorldToLocal(pos1,pos2); dist = fGeo-> Plane((unsigned int)planeId). Cell((unsigned int)cellId). DistToReadOut(0.5*(hitlist[i]->fHits[j].fEntryZ + hitlist[i]->fHits[j].fExitZ)); trkid = hitlist[i]->fHits[j].fTrackId; photrans::PhotonSignal photon(time,0.,int(edep*fMeVToNPhot*fGeVToMeV*fCollectionEff*exp(-dist/fAttenLength))); photon.SetPlaneCell(planeId, cellId); photon.SetTrackId(trkid); if(photon.NPhoton() > 0.) photlist.push_back(photon); } } if(msglevel>=INFO)cout << "Read " << nhits << " FLS hits" << endl; if(msglevel>=INFO)cout << "Created " << photlist.size() << " photon signals" << endl; for (unsigned int i=0; i<photlist.size(); ++i) evt.Raw().Put(photlist[i],""); photlist.clear(); return jobc::kPassed; }

void SimpleTransport::Update (  const cfg::Config &  c  )  [virtual]

Implements cfg::Observer. Definition at line 51 of file SimpleTransport.cxx. References fAttenLength, fCollectionEff, fMeVToNPhot, and msglevel. { c("MeVToNPhot"). Get(fMeVToNPhot); c("CollectionEff").Get(fCollectionEff); c("AttenLength"). Get(fAttenLength); c("MessageLevel"). Get(msglevel); }

---

Member Data Documentation

float photrans::SimpleTransport::fAttenLength [private]

normalization factor Definition at line 42 of file SimpleTransport.h. Referenced by Reco(), and Update().

float photrans::SimpleTransport::fCollectionEff [private]

conversion factor MeV->NPhoton Definition at line 41 of file SimpleTransport.h. Referenced by Reco(), and Update().

geo::Geometry* photrans::SimpleTransport::fGeo [private]

the attenuation length in cm Definition at line 43 of file SimpleTransport.h. Referenced by Reco().

float photrans::SimpleTransport::fMeVToNPhot [private]

Definition at line 40 of file SimpleTransport.h. Referenced by Reco(), and Update().

int photrans::SimpleTransport::msglevel [private]

Definition at line 45 of file SimpleTransport.h. Referenced by Reco(), and Update().

---

The documentation for this class was generated from the following files:

• SimpleTransport.h
• SimpleTransport.cxx

---