ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Mods/src/PhotonTreeWriter.cc
(Generate patch)

Comparing UserCode/MitPhysics/Mods/src/PhotonTreeWriter.cc (file contents):
Revision 1.78 by veverka, Fri Dec 13 02:58:46 2013 UTC vs.
Revision 1.79 by mingyang, Mon Dec 16 16:52:36 2013 UTC

# Line 131 | Line 131 | PhotonTreeWriter::PhotonTreeWriter(const
131    fElectronMVAWeights_Subdet2Pt20ToInf(""),
132  
133    fTheRhoType(RhoUtilities::DEFAULT),
134 +
135 +  fdor9rescale                   (false),
136 +  fp0b                           (0.),
137 +  fp1b                           (1.),
138 +  fp0e                           (0.),
139 +  fp1e                           (1.),
140 +
141    fProcessedEvents(0)
142  
143   {
# Line 903 | Line 910 | void PhotonTreeWriter::Process()
910                                         fPhfixph, fPhfixele, fTracks, fPV,
911                                         fPFCands, rho, fFillClusterArrays,
912                                         fPhotons, fPFSuperClusters,
913 <                                       fElectrons, fConversions, bsp,
913 >                                       fElectrons, fConversions, bsp,fIsData, fdor9rescale, fp0b, fp1b,fp0e,fp1e,
914                                         fApplyElectronVeto, realVtx);
915      fDiphotonEvent->photons[1].SetVars(phSoft, conv2, ele2, pfsc2, phgen2,
916                                         fPhfixph, fPhfixele, fTracks, fPV,
917                                         fPFCands, rho, fFillClusterArrays,
918                                         fPhotons, fPFSuperClusters,
919                                         fElectrons, fConversions,
920 <                                       bsp, fApplyElectronVeto, realVtx);
920 >                                       bsp,fIsData, fdor9rescale, fp0b, fp1b,fp0e,fp1e, fApplyElectronVeto, realVtx);
921      
922      Float_t ph1ecor    = fDiphotonEvent->photons[0].Ecor();
923      Float_t ph1ecorerr = fDiphotonEvent->photons[0].Ecorerr();
# Line 1132 | Line 1139 | void PhotonTreeWriter::Process()
1139      fSinglePhoton->SetVars(ph, conv, ele, pfsc, phgen, fPhfixph, fPhfixele,
1140                             fTracks, fPV, fPFCands, rho, fFillClusterArrays,
1141                             fPhotons, fPFSuperClusters, fElectrons, fConversions,
1142 <                           bsp, fApplyElectronVeto);
1142 >                           bsp, fIsData, fdor9rescale, fp0b, fp1b,fp0e,fp1e, fApplyElectronVeto);
1143      hCiCTupleSingle->Fill();
1144    }
1145  
# Line 1422 | Line 1429 | PhotonTreeWriterPhoton<NClus>::SetVars(c
1429                                         const SuperClusterCol* scs,
1430                                         const ElectronCol* els,
1431                                         const DecayParticleCol *convs,
1432 <                                       const BaseVertex *bs,
1432 >                                       const BaseVertex *bs,bool isdata,
1433 >                                       bool dor9rescale, double p0b, double p1b, double  p0e,  double  p1e,
1434                                         Bool_t applyElectronVeto,
1435                                         const Vertex* realVtx)
1436   {
# Line 1515 | Line 1523 | PhotonTreeWriterPhoton<NClus>::SetVars(c
1523      // -----------------------------------------------------
1524      // PF-CiC4 Debug Stuff
1525      std::vector<double> debugVals;
1526 <    PhotonTools::PassCiCPFIsoSelection(p, vtx, fPFCands, vtxCol, rho, 20., &debugVals);
1526 >    PhotonTools::PassCiCPFIsoSelection(p, vtx, fPFCands, vtxCol, rho, 20., dor9rescale,p0b,p1b,p0e,p1e,&debugVals);
1527      if( debugVals.size() == 13 ) {
1528        pfcic4_tIso1   = debugVals[0];
1529        pfcic4_tIso2   = debugVals[1];
# Line 1577 | Line 1585 | PhotonTreeWriterPhoton<NClus>::SetVars(c
1585      combiso2 = -99.;
1586    }
1587  
1588 +  if(dor9rescale &&  !isdata){
1589 +    if(s->AbsEta()<1.5){
1590 +      r9 = p0b + p1b * r9;
1591 +    }else{
1592 +      r9 = p0e + p1e * r9;
1593 +    }
1594 +  }
1595 +  
1596    // TODO: fix the bug with supercluster index
1597    scindex = PhotonTreeWriter::IndexOfNearestSuperClusterInCollection(s, scs);
1598    /// DEBUG

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines