70 likes | 197 Views
Jupiter/Satellites update. A.Miyamoto 2006-06-12. J4OptDet. Purpose: Monitor particle distributions at arbitrary points in a detector Especially for background studies Mechanism Add J4OptDet, J4OptDetSD, J4OptDetHit classes in lcexp.
E N D
Jupiter/Satellites update A.Miyamoto 2006-06-12
J4OptDet • Purpose: • Monitor particle distributions at arbitrary points in a detector • Especially for background studies • Mechanism • Add J4OptDet, J4OptDetSD, J4OptDetHit classes in lcexp. • Optional detectors are introduced in IR/VTX/IT by a call to J4OptDet in Assemble • Number of detectors , geometries, etc are defined by a geometry data file. • Hit information = J4VHit + Detector ID • Detector ID is defined by a geometry data file. • Any call to ProcessHit create a hit disregarding energy deposit or particle ID. • Need to define a thin geometry to save a data size. • Hit created by a same mother particle should be resolved in user analysis • S4OptDetExactHit is used to save hit information. • No MCTruth module is created. Use direct access to JSFJ4 buffer.
Howto define optional detectors • Make sure not to overlap with existing volumes. • Sample definition J4IR.OptDet.N: 4 : Number of optional detector in IR region ( in this case 4 ) J4IR.OptDet.ID: 31 32 33 34 : Detector ID’s assigned to the exact hit J4IR.OptDet.0.ShapeID: 0 : ShapeID of 0th detector (0=Tubs, 1=Cons) J4IR.OptDet.0.Geometry: 0.0 16.19 398.8 398.801 : Geometry parameter for 0th detetor Tubs: rmin, rmax, zmin, zmax ( in cm) Cons: rmin1, rmax1, rmin2, rmax2, zmin, zmax J4IR.OptDet.0.VisAtt:true : Visibility flag ( true or false ) J4IR.OptDet.0.Color: 1.0 0.5 0.5 1.0 : Color of 0th detector Indicates 0th detector Detector description line repeats • To define in VTX/IR, use parameter name J4VTX.OptDet…. J4IT.OptDet…
A sample code to access hits JSFJ4 *jsfj4=(JSFJ4*)gJSF->FindModule("JSFJ4"); JSFJ4Buf *jsfj4buf=(JSFJ4Buf*)jsfj4->EventBuf(); TObjArray *opthits=jsfj4buf->GetComponent("OptDetHits"); std::cerr << "Number of opt hits is " << opthits->GetEntries() << std::endl; TIter onext(opthits); S4OptDetExactHit *ohit; while( (ohit=(S4OptDetExactHit*)onext() ) ) { ohit->DebugPrint(); std::cerr << ohit->GetDetID(); // Detector ID assigned by the input data << ohit->GetDetSerialID(); // Serial number from 0 to Max-1. // All OptDet has a serial number disregarding IR/VTX/IT, … }
FCAL/BCAL • Sensitive detectors are defined. • Energy Deposit, DE weighted positions are calculated by J4IRXXXSD • Parameters for segmentation • J4IR.YYY.NDIVR: Number of radial division • J4IR.yyy.NDIVPhi: Number of phi division ( decided by center line ) yyy = FCAL.Front, FCAL.Tail, BCAL • Typical values for ~1cm is defined in jun06 geometry file. Cell boundary is not aligned
FCAL/BCAL PreHit • J4IRFCAL(BCAL)PreHit is defined • No posthit classes • Uses of J4HistoryKeeper are turned off by J4IR.FCAL.KeepHistory:NoJ4IR.BCAL.KeepHistory:No
Jun06 data production • After fixing initial problems, • started since Jun 9(Friday) • Using Release 1.29 version • Produced data available • Gamma: 1, 2, 5, 10, 20, 50, 100, 200GeV 1K each • K0L: 1,2,5,10,20,50,100,200 GeV ( 5 and 10 GeV sample, 10K ) • Pi0: in progress: using 1x1cm segmentation • UDS91: trying to make 10K sample by 2x2, 1x1 seg. • CAIN sample O(10) larger than previously. • No FCAL/BCAL prehit. • Generate both +Z/-Z side